Pycallingcards: an integrated environment for visualizing, analyzing, and interpreting Calling Cards data.

MOTIVATION: Unraveling the transcriptional programs that control how cells divide, differentiate, and respond to their environments requires a precise understanding of transcription factors' (TFs) DNA-binding activities. Calling cards (CC) technology uses transposons to capture transient TF binding events at one instant in time and then read them out at a later time. This methodology can also be used to simultaneously measure TF binding and mRNA expression from single-cell CC and to record and integrate TF binding events across time in any cell type of interest without the need for purification. Despite these advantages, there has been a lack of dedicated bioinformatics tools for the detailed analysis of CC data. RESULTS: We introduce Pycallingcards, a comprehensive Python module specifically designed for the analysis of single-cell and bulk CC data across multiple species. Pycallingcards introduces two innovative peak callers, CCcaller and MACCs, enhancing the accuracy and speed of pinpointing TF binding sites from CC data. Pycallingcards offers a fully integrated environment for data visualization, motif finding, and comparative analysis with RNA-seq and ChIP-seq datasets. To illustrate its practical application, we have reanalyzed previously published mouse cortex and glioblastoma datasets. This analysis revealed novel cell-type-specific binding sites and potential sex-linked TF regulators, furthering our understanding of TF binding and gene expression relationships. Thus, Pycallingcards, with its user-friendly design and seamless interface with the Python data science ecosystem, stands as a critical tool for advancing the analysis of TF functions via CC data. AVAILABILITY AND IMPLEMENTATION: Pycallingcards can be accessed on the GitHub repository: https://github.com/The-Mitra-Lab/pycallingcards.

Zinc cluster transcription factors frequently activate target genes using a non-canonical half-site binding mode.

Gene expression changes are orchestrated by transcription factors (TFs), which bind to DNA to regulate gene expression. It remains surprisingly difficult to predict basic features of the transcriptional process, including in vivo TF occupancy. Existing thermodynamic models of TF function are often not concordant with experimental measurements, suggesting undiscovered biology. Here, we analyzed one of the most well-studied TFs, the yeast zinc cluster Gal4, constructed a Shea-Ackers thermodynamic model to describe its binding, and compared the results of this model to experimentally measured Gal4p binding in vivo. We found that at many promoters, the model predicted no Gal4p binding, yet substantial binding was observed. These outlier promoters lacked canonical binding motifs, and subsequent investigation revealed Gal4p binds unexpectedly to DNA sequences with high densities of its half site (CGG). We confirmed this novel mode of binding through multiple experimental and computational paradigms; we also found most other zinc cluster TFs we tested frequently utilize this binding mode, at 27% of their targets on average. Together, these results demonstrate a novel mode of binding where zinc clusters, the largest class of TFs in yeast, bind DNA sequences with high densities of half sites.

A Suite of New Strain Construction Vectors for Gene Expression Knockdown in Budding Yeast.

Numerous tools for gene expression knockdown have been developed and characterized in the model organism Saccharomyces cerevisiae and extended to facilitate studies in multicellular models. To comparatively evaluate the efficacy of these approaches, we systematically applied seven such published constitutive and inducible knockdown strategies to a panel of essential genes encoding nuclear-localized proteins. In this effort, we created the CEAS (C-SWAT for Essential Allele Strains) collection, a suite of tagging vectors for improved utility and ease of strain construction. Of particular note, we adapted an improved auxin inducible degron (AID) protein degradation strategy previously available only in mammalian tissue culture for one-step strain construction in budding yeast by leveraging both the C-SWAT system and CRISPR/Cas9 editing. Taken together, this work presents a toolbox for endogenous gene expression knockdown and allows us to make recommendations on the efficacy and applicability of these tools for the perturbation of essential genes.

Quantitative analysis of transcription factor binding and expression using calling cards reporter arrays.

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 KBT for sites that are positioned 'out of phase'. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a 'TF collective', an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

Homotypic cooperativity and collective binding are determinants of bHLH specificity and function.

Eukaryotic cells express transcription factor (TF) paralogues that bind to nearly identical DNA sequences in vitro but bind at different genomic loci and perform different functions in vivo. Predicting how 2 paralogous TFs bind in vivo using DNA sequence alone is an important open problem. Here, we analyzed 2 yeast bHLH TFs, Cbf1p and Tye7p, which have highly similar binding preferences in vitro, yet bind at almost completely nonoverlapping target loci in vivo. We dissected the determinants of specificity for these 2 proteins by making a number of chimeric TFs in which we swapped different domains of Cbf1p and Tye7p and determined the effects on in vivo binding and cellular function. From these experiments, we learned that the Cbf1p dimer achieves its specificity by binding cooperatively with other Cbf1p dimers bound nearby. In contrast, we found that Tye7p achieves its specificity by binding cooperatively with 3 other DNA-binding proteins, Gcr1p, Gcr2p, and Rap1p. Remarkably, most promoters (63%) that are bound by Tye7p do not contain a consensus Tye7p binding site. Using this information, we were able to build simple models to accurately discriminate bound and unbound genomic loci for both Cbf1p and Tye7p. We then successfully reprogrammed the human bHLH NPAS2 to bind Cbf1p in vivo targets and a Tye7p target intergenic region to be bound by Cbf1p. These results demonstrate that the genome-wide binding targets of paralogous TFs can be discriminated using sequence information, and provide lessons about TF specificity that can be applied across the phylogenetic tree.

Inositol polyphosphates regulate and predict yeast pseudohyphal growth phenotypes.

Pseudohyphal growth is a nutrient-regulated program in which budding yeast form multicellular filaments of elongated and connected cells. Filamentous growth is required for virulence in pathogenic fungi and provides an informative model of stress-responsive signaling. The genetics and regulatory networks modulating pseudohyphal growth have been studied extensively, but little is known regarding the changes in metabolites that enable pseudohyphal filament formation. Inositol signaling molecules are an important class of metabolite messengers encompassing highly phosphorylated and diffusible inositol polyphosphates (InsPs). We report here that the InsP biosynthesis pathway is required for wild-type pseudohyphal growth. Under nitrogen-limiting conditions that can induce filamentation, InsPs exhibit characteristic profiles, distinguishing the InsP7 pyrophosphate isoforms 1PP-InsP5 and 5PP-InsP5. Deletion and overexpression analyses of InsP kinases identify elevated levels of 5PP-InsP5 relative to 1PP-InsP5 in mutants exhibiting hyper-filamentous growth. Overexpression of KCS1, which promotes formation of inositol pyrophosphates, is sufficient to drive pseudohyphal filamentation on medium with normal nitrogen levels. We find that the kinases Snf1p (AMPK), Kss1p, and Fus3p (MAPKs), required for wild-type pseudohyphal growth, are also required for wild-type InsP levels. Deletion analyses of the corresponding kinase genes indicate elevated InsP3 levels and an absence of exaggerated 5PP-InsP5 peaks in trace profiles from snf1Δ/Δ and kss1Δ/Δ mutants exhibiting decreased pseudohyphal filamentation. Elevated 5PP-InsP5:1PP-InsP5 ratios are present in the hyperfilamentous fus3 deletion mutant. Collectively, the data identify the presence of elevated 5PP-InsP5 levels relative to other inositol pyrophosphates as an in vivo marker of hyper-filamentous growth, while providing initial evidence for the regulation of InsP signaling by pseudohyphal growth kinases.

Formation of an Intramolecular Periplasmic Disulfide Bond in TcpP Protects TcpP and TcpH from Degradation in Vibrio cholerae.

UNLABELLED: TcpP and ToxR coordinately regulate transcription of toxT, the master regulator of numerous virulence factors in Vibrio cholerae. TcpP and ToxR are membrane-localized transcription factors, each with a periplasmic domain containing two cysteines. In ToxR, these cysteines form an intramolecular disulfide bond and a cysteine-to-serine substitution affects activity. We determined that the two periplasmic cysteines of TcpP also form an intramolecular disulfide bond. Disruption of this intramolecular disulfide bond by mutation of either cysteine resulted in formation of intermolecular disulfide bonds. Furthermore, disruption of the intramolecular disulfide bond in TcpP decreased the stability of TcpP. While the decreased stability of TcpP-C207S resulted in a nearly complete loss of toxT activation and cholera toxin (CT) production, the second cysteine mutant, TcpP-C218S, was partially resistant to proteolytic degradation and maintained ∼50% toxT activation capacity. TcpP-C218S was also TcpH independent, since deletion of tcpH did not affect the stability of TcpP-C218S, whereas wild-type TcpP was degraded in the absence of TcpH. Finally, TcpH was also unstable when intramolecular disulfides could not be formed in TcpP, suggesting that the single periplasmic cysteine in TcpH may assist with disulfide bond formation in TcpP by interacting with the periplasmic cysteines of TcpP. Consistent with this finding, a TcpH-C114S mutant was unable to stabilize TcpP and was itself unstable. Our findings demonstrate a periplasmic disulfide bond in TcpP is critical for TcpP stability and virulence gene expression. IMPORTANCE: The Vibrio cholerae transcription factor TcpP, in conjunction with ToxR, regulates transcription of toxT, the master regulator of numerous virulence factors in Vibrio cholerae. TcpP is a membrane-localized transcription factor with a periplasmic domain containing two cysteines. We determined that the two periplasmic cysteines of TcpP form an intramolecular disulfide bond and disruption of the intramolecular disulfide bond in TcpP decreased the stability of TcpP and reduced virulence gene expression. Normally TcpH, another membrane-localized periplasmic protein, protects TcpP from degradation. However, we found that TcpH was also unstable when intramolecular disulfides could not be formed in TcpP, indicating that the periplasmic cysteines of TcpP are required for functional interaction with TcpH and that this interaction is required for both TcpP and TcpH stability.

Large-Scale Analysis of Kinase Signaling in Yeast Pseudohyphal Development Identifies Regulation of Ribonucleoprotein Granules.

Yeast pseudohyphal filamentation is a stress-responsive growth transition relevant to processes required for virulence in pathogenic fungi. Pseudohyphal growth is controlled through a regulatory network encompassing conserved MAPK (Ste20p, Ste11p, Ste7p, Kss1p, and Fus3p), protein kinase A (Tpk2p), Elm1p, and Snf1p kinase pathways; however, the scope of these pathways is not fully understood. Here, we implemented quantitative phosphoproteomics to identify each of these signaling networks, generating a kinase-dead mutant in filamentous S. cerevisiae and surveying for differential phosphorylation. By this approach, we identified 439 phosphoproteins dependent upon pseudohyphal growth kinases. We report novel phosphorylation sites in 543 peptides, including phosphorylated residues in Ras2p and Flo8p required for wild-type filamentous growth. Phosphoproteins in these kinase signaling networks were enriched for ribonucleoprotein (RNP) granule components, and we observe co-localization of Kss1p, Fus3p, Ste20p, and Tpk2p with the RNP component Igo1p. These kinases localize in puncta with GFP-visualized mRNA, and KSS1 is required for wild-type levels of mRNA localization in RNPs. Kss1p pathway activity is reduced in lsm1Δ/Δ and pat1Δ/Δ strains, and these genes encoding P-body proteins are epistatic to STE7. The P-body protein Dhh1p is also required for hyphal development in Candida albicans. Collectively, this study presents a wealth of data identifying the yeast phosphoproteome in pseudohyphal growth and regulatory interrelationships between pseudohyphal growth kinases and RNPs.

The yeast Sks1p kinase signaling network regulates pseudohyphal growth and glucose response.

The yeast Saccharomyces cerevisiae undergoes a dramatic growth transition from its unicellular form to a filamentous state, marked by the formation of pseudohyphal filaments of elongated and connected cells. Yeast pseudohyphal growth is regulated by signaling pathways responsive to reductions in the availability of nitrogen and glucose, but the molecular link between pseudohyphal filamentation and glucose signaling is not fully understood. Here, we identify the glucose-responsive Sks1p kinase as a signaling protein required for pseudohyphal growth induced by nitrogen limitation and coupled nitrogen/glucose limitation. To identify the Sks1p signaling network, we applied mass spectrometry-based quantitative phosphoproteomics, profiling over 900 phosphosites for phosphorylation changes dependent upon Sks1p kinase activity. From this analysis, we report a set of novel phosphorylation sites and highlight Sks1p-dependent phosphorylation in Bud6p, Itr1p, Lrg1p, Npr3p, and Pda1p. In particular, we analyzed the Y309 and S313 phosphosites in the pyruvate dehydrogenase subunit Pda1p; these residues are required for pseudohyphal growth, and Y309A mutants exhibit phenotypes indicative of impaired aerobic respiration and decreased mitochondrial number. Epistasis studies place SKS1 downstream of the G-protein coupled receptor GPR1 and the G-protein RAS2 but upstream of or at the level of cAMP-dependent PKA. The pseudohyphal growth and glucose signaling transcription factors Flo8p, Mss11p, and Rgt1p are required to achieve wild-type SKS1 transcript levels. SKS1 is conserved, and deletion of the SKS1 ortholog SHA3 in the pathogenic fungus Candida albicans results in abnormal colony morphology. Collectively, these results identify Sks1p as an important regulator of filamentation and glucose signaling, with additional relevance towards understanding stress-responsive signaling in C. albicans.

Genetic networks inducing invasive growth in Saccharomyces cerevisiae identified through systematic genome-wide overexpression.

The budding yeast Saccharomyces cerevisiae can respond to nutritional and environmental stress by implementing a morphogenetic program wherein cells elongate and interconnect, forming pseudohyphal filaments. This growth transition has been studied extensively as a model signaling system with similarity to processes of hyphal development that are linked with virulence in related fungal pathogens. Classic studies have identified core pseudohyphal growth signaling modules in yeast; however, the scope of regulatory networks that control yeast filamentation is broad and incompletely defined. Here, we address the genetic basis of yeast pseudohyphal growth by implementing a systematic analysis of 4909 genes for overexpression phenotypes in a filamentous strain of S. cerevisiae. Our results identify 551 genes conferring exaggerated invasive growth upon overexpression under normal vegetative growth conditions. This cohort includes 79 genes lacking previous phenotypic characterization. Pathway enrichment analysis of the gene set identifies networks mediating mitogen-activated protein kinase (MAPK) signaling and cell cycle progression. In particular, overexpression screening suggests that nuclear export of the osmoresponsive MAPK Hog1p may enhance pseudohyphal growth. The function of nuclear Hog1p is unclear from previous studies, but our analysis using a nuclear-depleted form of Hog1p is consistent with a role for nuclear Hog1p in repressing pseudohyphal growth. Through epistasis and deletion studies, we also identified genetic relationships with the G2 cyclin Clb2p and phenotypes in filamentation induced by S-phase arrest. In sum, this work presents a unique and informative resource toward understanding the breadth of genes and pathways that collectively constitute the molecular basis of filamentation.

A profile of differentially abundant proteins at the yeast cell periphery during pseudohyphal growth.

Yeast filamentous growth is a stress response to conditions of nitrogen deprivation, wherein yeast colonies form pseudohyphal filaments of elongated and connected cells. As proteins mediating adhesion and transport are required for this growth transition, we expect that the protein complement at the yeast cell periphery plays a critical and tightly regulated role in pseudohyphal filamentation. To identify proteins differentially abundant at the yeast cell periphery during pseudohyphal growth, we generated quantitative proteomic profiles of plasma membrane protein preparations under conditions of vegetative growth and filamentation. By isobaric tags for relative and absolute quantification chemistry and two-dimensional liquid chromatography-tandem mass spectrometry, we profiled 2463 peptides and 356 proteins, identifying 11 differentially abundant proteins that localize to the yeast cell periphery. This protein set includes Ylr414cp, herein renamed Pun1p, a previously uncharacterized protein localized to the plasma membrane compartment of Can1. Pun1p abundance is doubled under conditions of nitrogen stress, and deletion of PUN1 abolishes filamentous growth in haploids and diploids; pun1Delta mutants are noninvasive, lack surface-spread filamentation, grow slowly, and exhibit impaired cell adhesion. Conversely, overexpression of PUN1 results in exaggerated cell elongation under conditions of nitrogen stress. PUN1 contributes to yeast nitrogen signaling, as pun1Delta mutants misregulate amino acid biosynthetic genes during nitrogen stress. By chromatin immunoprecipitation and reverse transcription-PCR, we find that the filamentous growth factor Mss11p directly binds the PUN1 promoter and regulates its transcription. In total, this study provides the first profile of differential protein abundance during pseudohyphal growth, identifying a previously uncharacterized membrane compartment of Can1 protein required for wild-type nitrogen signaling and filamentous growth.