Mapping yeast transcriptional networks

"Hit-and-Run" leaves its mark: catalyst transcription factors and chromatin modification

Understanding how transcription factor (TF) binding is related to gene regulation is a moving target. We recently uncovered genome‐wide evidence for a “Hit‐and‐Run” model of transcription. In this model, a master TF “hits” a target promoter to initiate a rapid response to a signal. As the “hit” is transient, the model invokes recruitment of partner TFs to sustain transcription over time. Following the “run”, the master TF “hits” other targets to propagate the response genome‐wide. As such, a TF may act as a “catalyst” to mount a broad and acute response in cells that first sense the signal, while the recruited TF partners promote long‐term adaptive behavior in the whole organism. This “Hit‐and‐Run” model likely has broad relevance, as TF perturbation studies across eukaryotes show small overlaps between TF‐regulated and TF‐bound genes, implicating transient TF‐target binding. Here, we explore this “Hit‐and‐Run” model to suggest molecular mechanisms and its biological relevance.

Genome-Wide Identification of the Target Genes of AP2-O, a Plasmodium AP2-Family Transcription Factor

Stage-specific transcription is a fundamental biological process in the life cycle of the Plasmodium parasite. Proteins containing the AP2 DNA-binding domain are responsible for stage-specific transcriptional regulation and belong to the only known family of transcription factors in Plasmodium parasites. Comprehensive identification of their target genes will advance our understanding of the molecular basis of stage-specific transcriptional regulation and stage-specific parasite development. AP2-O is an AP2 family transcription factor that is expressed in the mosquito midgut-invading stage, called the ookinete, and is essential for normal morphogenesis of this stage. In this study, we identified the genome-wide target genes of AP2-O by chromatin immunoprecipitation-sequencing and elucidate how this AP2 family transcription factor contributes to the formation of this motile stage. The analysis revealed that AP2-O binds specifically to the upstream genomic regions of more than 500 genes, suggesting that approximately 10% of the parasite genome is directly regulated by AP2-O. These genes are involved in distinct biological processes such as morphogenesis, locomotion, midgut penetration, protection against mosquito immunity and preparation for subsequent oocyst development. This direct and global regulation by AP2-O provides a model for gene regulation in Plasmodium parasites and may explain how these parasites manage to control their complex life cycle using a small number of sequence-specific AP2 transcription factors.

Although malarial parasites have a complex life cycle, they harbor only 30 transcription factors in their genome. The majority of these transcription factors belong to a single family referred to as the AP2 family. Our previous study suggested that stage-specific AP2 family transcription factors have critical roles in maintaining the Plasmodium parasite life cycle. However, it remains fairly elusive as to how these transcription factors regulate each stage. AP2-O is an AP2 family transcription factor that is expressed during the mosquito midgut-invading stage, the ookinete, and is essential for normal development of this stage. In the present study, we identified the entire set of AP2-O target genes to elucidate how this AP2 family transcription factor contributes to the formation of this stage. Our results showed that AP2-O directly regulates 10% of the parasite genome and is involved in the whole process of mosquito midgut-invasion by ookinetes. The global and comprehensive regulation by the AP2 family transcription factor that we revealed provides a model for transcriptional regulation of this parasite and may explain how malarial parasites regulate their complex life cycle using a small number of sequence-specific transcription factors.

Role of intragenic binding of cAMP responsive protein (CRP) in regulation of the succinate dehydrogenase genes Rv0249c-Rv0247c in TB complex mycobacteria

Bacterial pathogens adapt to changing environments within their hosts, and the signaling molecule adenosine 3', 5'-cyclic monophosphate (cAMP) facilitates this process. In this study, we characterized in vivo DNA binding and gene regulation by the cAMP-responsive protein CRP in M. bovis BCG as a model for tuberculosis (TB)-complex bacteria. Chromatin immunoprecipitation followed by deep-sequencing (ChIP-seq) showed that CRP associates with ∼900 DNA binding regions, most of which occur within genes. The most highly enriched binding region was upstream of a putative copper transporter gene (ctpB), and crp-deleted bacteria showed increased sensitivity to copper toxicity. Detailed mutational analysis of four CRP binding sites upstream of the virulence-associated Rv0249c-Rv0247c succinate dehydrogenase genes demonstrated that CRP directly regulates Rv0249c-Rv0247c expression from two promoters, one of which requires sequences intragenic to Rv0250c for maximum expression. The high percentage of intragenic CRP binding sites and our demonstration that these intragenic DNA sequences significantly contribute to biologically relevant gene expression greatly expand the genome space that must be considered for gene regulatory analyses in mycobacteria. These findings also have practical implications for an important bacterial pathogen in which identification of mutations that affect expression of drug target-related genes is widely used for rapid drug resistance screening.

Controllability analysis of transcriptional regulatory networks reveals circular control patterns among transcription factors

Transcriptional regulation is the most committed type of regulation in living cells where transcription factors (TFs) control the expression of their target genes and TF expression is controlled by other TFs forming complex transcriptional regulatory networks that can be highly interconnected. Here we analyze the topology and organization of nine transcriptional regulatory networks for E. coli, yeast, mouse and human, and we evaluate how the structure of these networks influences two of their key properties, namely controllability and stability. We calculate the controllability for each network as a measure of the organization and interconnectivity of the network. We find that the number of driver nodes nD needed to control the whole network is 64% of the TFs in the E. coli transcriptional regulatory network in contrast to only 17% for the yeast network, 4% for the mouse network and 8% for the human network. The high controllability (low number of drivers needed to control the system) in yeast, mouse and human is due to the presence of internal loops in their regulatory networks where the TFs regulate each other in a circular fashion. We refer to these internal loops as circular control motifs (CCM). The E. coli transcriptional regulatory network, which does not have any CCMs, shows a hierarchical structure of the transcriptional regulatory network in contrast to the eukaryal networks. The presence of CCMs also has influence on the stability of these networks, as the presence of cycles can be associated with potential unstable steady-states where even small changes in binding affinities can cause dramatic rearrangements of the state of the network.

ChIP-seq and in vivo transcriptome analyses of the Aspergillus fumigatus SREBP SrbA reveals a new regulator of the fungal hypoxia response and virulence

The Aspergillus fumigatus sterol regulatory element binding protein (SREBP) SrbA belongs to the basic Helix-Loop-Helix (bHLH) family of transcription factors and is crucial for antifungal drug resistance and virulence. The latter phenotype is especially striking, as loss of SrbA results in complete loss of virulence in murine models of invasive pulmonary aspergillosis (IPA). How fungal SREBPs mediate fungal virulence is unknown, though it has been suggested that lack of growth in hypoxic conditions accounts for the attenuated virulence. To further understand the role of SrbA in fungal infection site pathobiology, chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) was used to identify genes under direct SrbA transcriptional regulation in hypoxia. These results confirmed the direct regulation of ergosterol biosynthesis and iron uptake by SrbA in hypoxia and revealed new roles for SrbA in nitrate assimilation and heme biosynthesis. Moreover, functional characterization of an SrbA target gene with sequence similarity to SrbA identified a new transcriptional regulator of the fungal hypoxia response and virulence, SrbB. SrbB co-regulates genes involved in heme biosynthesis and demethylation of C4-sterols with SrbA in hypoxic conditions. However, SrbB also has regulatory functions independent of SrbA including regulation of carbohydrate metabolism. Loss of SrbB markedly attenuates A. fumigatus virulence, and loss of both SREBPs further reduces in vivo fungal growth. These data suggest that both A. fumigatus SREBPs are critical for hypoxia adaptation and virulence and reveal new insights into SREBPs' complex role in infection site adaptation and fungal virulence.

Despite improvements in diagnostics and antifungal drug treatments, mortality rates from invasive mold infections remain high. Defining the fungal adaptation and growth mechanisms at the infection site microenvironment is one research focus that is expected to improve treatment of established invasive fungal infections. The Aspergillus fumigatus transcription factor SrbA is a major regulator of the fungal response to hypoxia found at sites of invasive fungal growth in vivo. In this study, new insights into how SrbA mediates hypoxia adaptation and virulence were revealed through identification of direct transcriptional targets of SrbA under hypoxic conditions. A major novel finding from these studies is the identification of a critical role in fungal hypoxia adaptation and virulence of an SrbA target gene, srbB, which is also in the SREBP family. SrbB plays a major role in regulation of heme biosynthesis and carbohydrate metabolism early in the response to hypoxia. The discovery of SrbA-dependent regulation of srbB gene expression, and the target genes they regulate opens new avenues to understand how SREBPs and their target genes mediate adaptation to the in vivo infection site microenvironment and responses to current antifungal therapies.

Hit-and-run transcriptional control by bZIP1 mediates rapid nutrient signaling in Arabidopsis

The dynamic nature of gene regulatory networks allows cells to rapidly respond to environmental change. However, the underlying temporal connections are missed, even in kinetic studies, as transcription factor (TF) binding within at least one time point is required to identify primary targets. The TF-regulated but unbound genes are dismissed as secondary targets. Instead, we report that these genes comprise transient TF-target interactions most relevant to rapid signal transduction. We temporally perturbed a master TF (Basic Leucine Zipper 1, bZIP1) and the nitrogen (N) signal it transduces and integrated TF regulation and binding data from the same cell samples. Our enabling approach could identify primary TF targets based solely on gene regulation, in the absence of TF binding. We uncovered three classes of primary TF targets: (i) poised (TF-bound but not TF-regulated), (ii) stable (TF-bound and TF-regulated), and (iii) transient (TF-regulated but not TF-bound), the largest class. Unexpectedly, the transient bZIP1 targets are uniquely relevant to rapid N signaling in planta, enriched in dynamic N-responsive genes, and regulated by TF and N signal interactions. These transient targets include early N responders nitrate transporter 2.1 and NIN-like protein 3, bound by bZIP1 at 1-5 min, but not at later time points following TF perturbation. Moreover, promoters of these transient targets are uniquely enriched with cis-regulatory motifs coinherited with bZIP1 binding sites, suggesting a recruitment role for bZIP1. This transient mode of TF action supports a classic, but forgotten, "hit-and-run" transcription model, which enables a "catalyst TF" to activate a large set of targets within minutes of signal perturbation.

Discovering functions of unannotated genes from a transcriptome survey of wild fungal isolates

Most fungal genomes are poorly annotated, and many fungal traits of industrial and biomedical relevance are not well suited to classical genetic screens. Assigning genes to phenotypes on a genomic scale thus remains an urgent need in the field. We developed an approach to infer gene function from expression profiles of wild fungal isolates, and we applied our strategy to the filamentous fungus Neurospora crassa. Using transcriptome measurements in 70 strains from two well-defined clades of this microbe, we first identified 2,247 cases in which the expression of an unannotated gene rose and fell across N. crassa strains in parallel with the expression of well-characterized genes. We then used image analysis of hyphal morphologies, quantitative growth assays, and expression profiling to test the functions of four genes predicted from our population analyses. The results revealed two factors that influenced regulation of metabolism of nonpreferred carbon and nitrogen sources, a gene that governed hyphal architecture, and a gene that mediated amino acid starvation resistance. These findings validate the power of our population-transcriptomic approach for inference of novel gene function, and we suggest that this strategy will be of broad utility for genome-scale annotation in many fungal systems. IMPORTANCE Some fungal species cause deadly infections in humans or crop plants, and other fungi are workhorses of industrial chemistry, including the production of biofuels. Advances in medical and industrial mycology require an understanding of the genes that control fungal traits. We developed a method to infer functions of uncharacterized genes by observing correlated expression of their mRNAs with those of known genes across wild fungal isolates. We applied this strategy to a filamentous fungus and predicted functions for thousands of unknown genes. In four cases, we experimentally validated the predictions from our method, discovering novel genes involved in the metabolism of nutrient sources relevant for biofuel production, as well as colony morphology and starvation resistance. Our strategy is straightforward, inexpensive, and applicable for predicting gene function in many fungal species.

Saccharomyces genome database provides new regulation data

The Saccharomyces Genome Database (SGD; http://www.yeastgenome.org) is the community resource for genomic, gene and protein information about the budding yeast Saccharomyces cerevisiae, containing a variety of functional information about each yeast gene and gene product. We have recently added regulatory information to SGD and present it on a new tabbed section of the Locus Summary entitled 'Regulation'. We are compiling transcriptional regulator-target gene relationships, which are curated from the literature at SGD or imported, with permission, from the YEASTRACT database. For nearly every S. cerevisiae gene, the Regulation page displays a table of annotations showing the regulators of that gene, and a graphical visualization of its regulatory network. For genes whose products act as transcription factors, the Regulation page also shows a table of their target genes, accompanied by a Gene Ontology enrichment analysis of the biological processes in which those genes participate. We additionally synthesize information from the literature for each transcription factor in a free-text Regulation Summary, and provide other information relevant to its regulatory function, such as DNA binding site motifs and protein domains. All of the regulation data are available for querying, analysis and download via YeastMine, the InterMine-based data warehouse system in use at SGD.

Influence of RpoS, cAMP-receptor protein, and ppGpp on expression of the opgGH operon and osmoregulated periplasmic glucan content of Salmonella enterica serovar Typhimurium

A transcriptional fusion (opgG1::MudJ) to the opgGH operon of Salmonella enterica serovar Typhimurium (S. Typhimurium) LT2, isolated by resistance to mecillinam, was used to study the influence of global regulators RpoS, ppGpp, and cAMP/cAMP-receptor protein (CRP) on expression of the opgGH operon and osmoregulated periplasmic glucan (OPG) content. Neither high growth medium osmolarity nor absence of ppGpp or CRP had important effects on opgG1::MudJ expression in exponential cultures. However, under the same conditions, OPG content was strongly decreased by high osmolarity or cAMP/CRP defectiveness, and reduced to a half by lack of ppGpp. In stationary cultures, high osmolarity as well as CRP loss caused significant descents in opgG1::MudJ expression that were compensated by inactivation of RpoS sigma factor. No effect of RpoS inactivation on OPG content was observed. It is concluded that opgGH expression in S. Typhimurium is only slightly affected by high osmolarity, but is inversely modulated by RpoS level. On the other hand, osmolarity and the cAMP/CRP global regulatory system appear to control OPG content, either directly or indirectly, mainly at the post-transcriptional level.

Protective ability of subcellular extracts from Salmonella Enteritidis and from a rough isogenic mutant against salmonellosis in mice

We evaluated the efficacy of surface components enriched hot saline extracts (HE) from parental and two isogenic rough mutant strains of Salmonella Enteritidis as subcellular vaccine candidates. By a randomized mutagenesis approach from a clinical isolate of S. Enteritidis there were selected two rough mutants defective in LPS synthesis (R1 and R2 mutants). The mutations mapped to the wcaI gene and gmd gene, respectively, of the O-antigen gene cluster involved in O-antigen synthesis. BALB/c mice received intraperitoneally one single dose of 30 microg of HE from parental and mutant strains, and the protection against a lethal infection with S. Enteritidis was determined. In contrast to the wild type extract, immunization with rough extracts did not induce any distress symptoms in the mice. HE extract from wild type and R1 strains induced the highest immunogenic response with respect IFN-gamma eliciting splenic cells, in contrast with HE-R2. These results correlated with the obtained levels of protection. Thus, at day 63 post-infection, HE from parental strain rendered an 80% level of protection; HE-R1 conferred a 60% level of protection, whereas HE-R2 did not protect the mice. Any of the antigenic extracts elicited systemic IgG1 and IgG2a responses, although these antibodies did not, however, correlate with protection. These results put forward the importance of cellular immune response mediated by IFN-gamma in protection against salmonellosis. The significantly different protective capacity between HE extracts from both rough mutants suggest that other factors independent of the O-chain, like outer membrane proteins and fimbrial antigens, may be involved in protection. In summary, the HE is a good candidate acellular extract for evaluation of its protective ability against salmonellosis following vaccination in poultry.