Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examine interspecies differences in RNA-protein interactions using the conserved neuronal RNA-binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstitute the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

Species-aware DNA language models capture regulatory elements and their evolution

BACKGROUND

The rise of large-scale multi-species genome sequencing projects promises to shed new light on how genomes encode gene regulatory instructions. To this end, new algorithms are needed that can leverage conservation to capture regulatory elements while accounting for their evolution.

RESULTS

Here, we introduce species-aware DNA language models, which we trained on more than 800 species spanning over 500 million years of evolution. Investigating their ability to predict masked nucleotides from context, we show that DNA language models distinguish transcription factor and RNA-binding protein motifs from background non-coding sequence. Owing to their flexibility, DNA language models capture conserved regulatory elements over much further evolutionary distances than sequence alignment would allow. Remarkably, DNA language models reconstruct motif instances bound in vivo better than unbound ones and account for the evolution of motif sequences and their positional constraints, showing that these models capture functional high-order sequence and evolutionary context. We further show that species-aware training yields improved sequence representations for endogenous and MPRA-based gene expression prediction, as well as motif discovery.

CONCLUSIONS

Collectively, these results demonstrate that species-aware DNA language models are a powerful, flexible, and scalable tool to integrate information from large compendia of highly diverged genomes.

Improving the performance of supervised deep learning for regulatory genomics using phylogenetic augmentation

MOTIVATION

Supervised deep learning is used to model the complex relationship between genomic sequence and regulatory function. Understanding how these models make predictions can provide biological insight into regulatory functions. Given the complexity of the sequence to regulatory function mapping (the cis-regulatory code), it has been suggested that the genome contains insufficient sequence variation to train models with suitable complexity. Data augmentation is a widely used approach to increase the data variation available for model training, however current data augmentation methods for genomic sequence data are limited.

RESULTS

Inspired by the success of comparative genomics, we show that augmenting genomic sequences with evolutionarily related sequences from other species, which we term phylogenetic augmentation, improves the performance of deep learning models trained on regulatory genomic sequences to predict high-throughput functional assay measurements. Additionally, we show that phylogenetic augmentation can rescue model performance when the training set is down-sampled and permits deep learning on a real-world small dataset, demonstrating that this approach improves data efficiency. Overall, this data augmentation method represents a solution for improving model performance that is applicable to many supervised deep-learning problems in genomics.

AVAILABILITY AND IMPLEMENTATION

The open-source GitHub repository agduncan94/phylogenetic_augmentation_paper includes the code for rerunning the analyses here and recreating the figures.

Understanding species-specific and conserved RNA-protein interactions in vivo and in vitro

While evolution is often considered from a DNA- and protein-centric view, RNA-based regulation can also impact gene expression and protein sequences. Here we examined interspecies differences in RNA-protein interactions using the conserved neuronal RNA binding protein, Unkempt (UNK) as model. We find that roughly half of mRNAs bound in human are also bound in mouse. Unexpectedly, even when transcript-level binding was conserved across species differential motif usage was prevalent. To understand the biochemical basis of UNK-RNA interactions, we reconstituted the human and mouse UNK-RNA interactomes using a high-throughput biochemical assay. We uncover detailed features driving binding, show that in vivo patterns are captured in vitro, find that highly conserved sites are the strongest bound, and associate binding strength with downstream regulation. Furthermore, subtle sequence differences surrounding motifs are key determinants of species-specific binding. We highlight the complex features driving protein-RNA interactions and how these evolve to confer species-specific regulation.

Computational prediction and experimental validation identify functionally conserved lncRNAs from zebrafish to human

Functional studies of long noncoding RNAs (lncRNAs) have been hindered by the lack of methods to assess their evolution. Here we present lncRNA Homology Explorer (lncHOME), a computational pipeline that identifies a unique class of long noncoding RNAs (lncRNAs) with conserved genomic locations and patterns of RNA-binding protein (RBP) binding sites (coPARSE-lncRNAs). Remarkably, several hundred human coPARSE-lncRNAs can be evolutionarily traced to zebrafish. Using CRISPR-Cas12a knockout and rescue assays, we found that knocking out many human coPARSE-lncRNAs led to cell proliferation defects, which were subsequently rescued by predicted zebrafish homologs. Knocking down coPARSE-lncRNAs in zebrafish embryos caused severe developmental delays that were rescued by human homologs. Furthermore, we verified that human, mouse and zebrafish coPARSE-lncRNA homologs tend to bind similar RBPs with their conserved functions relying on specific RBP-binding sites. Overall, our study demonstrates a comprehensive approach for studying the functional conservation of lncRNAs and implicates numerous lncRNAs in regulating vertebrate physiology.

Fusarium oxysporum f. sp. niveum Pumilio 1 Regulates Virulence on Watermelon through Interacting with the ARP2/3 Complex and Binding to an A-Rich Motif in the 3' UTR of Diverse Transcripts

Fusarium oxysporum f. sp. niveum (Fon), a soilborne phytopathogenic fungus, causes watermelon Fusarium wilt, resulting in serious yield losses worldwide. However, the underlying molecular mechanism of Fon virulence is largely unknown. The present study investigated the biological functions of six FonPUFs, encoding RNA binding Pumilio proteins, and especially explored the molecular mechanism of FonPUF1 in Fon virulence. A series of phenotypic analyses indicated that FonPUFs have distinct but diverse functions in vegetative growth, asexual reproduction, macroconidia morphology, spore germination, cell wall, or abiotic stress response of Fon. Notably, the deletion of FonPUF1 attenuates Fon virulence by impairing the invasive growth and colonization ability inside the watermelon plants. FonPUF1 possesses RNA binding activity, and its biochemical activity and virulence function depend on the RNA recognition motif or Pumilio domains. FonPUF1 associates with the actin-related protein 2/3 (ARP2/3) complex by interacting with FonARC18, which is also required for Fon virulence and plays an important role in regulating mitochondrial functions, such as ATP generation and reactive oxygen species production. Transcriptomic profiling of ΔFonPUF1 identified a set of putative FonPUF1-dependent virulence-related genes in Fon, possessing a novel A-rich binding motif in the 3' untranslated region (UTR), indicating that FonPUF1 participates in additional mechanisms critical for Fon virulence. These findings highlight the functions and molecular mechanism of FonPUFs in Fon virulence. IMPORTANCE Fusarium oxysporum is a devastating plant-pathogenic fungus that causes vascular wilt disease in many economically important crops, including watermelon, worldwide. F. oxysporum f. sp. nievum (Fon) causes serious yield loss in watermelon production. However, the molecular mechanism of Fusarium wilt development by Fon remains largely unknown. Here, we demonstrate that six putative Pumilio proteins-encoding genes (FonPUFs) differentially operate diverse basic biological processes, including stress response, and that FonPUF1 is required for Fon virulence. Notably, FonPUF1 possesses RNA binding activity and associates with the actin-related protein 2/3 complex to control mitochondrial functions. Furthermore, FonPUF1 coordinates the expression of a set of putative virulence-related genes in Fon by binding to a novel A-rich motif present in the 3' UTR of a diverse set of target mRNAs. Our study disentangles the previously unexplored molecular mechanism involved in regulating Fon virulence, providing a possibility for the development of novel strategies for disease management.

Nutritional compensation of the circadian clock is a conserved process influenced by gene expression regulation and mRNA stability

Compensation is a defining principle of a true circadian clock, where its approximately 24-hour period length is relatively unchanged across environmental conditions. Known compensation effectors directly regulate core clock factors to buffer the oscillator's period length from variables in the environment. Temperature Compensation mechanisms have been experimentally addressed across circadian model systems, but much less is known about the related process of Nutritional Compensation, where circadian period length is maintained across physiologically relevant nutrient levels. Using the filamentous fungus Neurospora crassa, we performed a genetic screen under glucose and amino acid starvation conditions to identify new regulators of Nutritional Compensation. Our screen uncovered 16 novel mutants, and together with 4 mutants characterized in prior work, a model emerges where Nutritional Compensation of the fungal clock is achieved at the levels of transcription, chromatin regulation, and mRNA stability. However, eukaryotic circadian Nutritional Compensation is completely unstudied outside of Neurospora. To test for conservation in cultured human cells, we selected top hits from our fungal genetic screen, performed siRNA knockdown experiments of the mammalian orthologs, and characterized the cell lines with respect to compensation. We find that the wild-type mammalian clock is also compensated across a large range of external glucose concentrations, as observed in Neurospora, and that knocking down the mammalian orthologs of the Neurospora compensation-associated genes CPSF6 or SETD2 in human cells also results in nutrient-dependent period length changes. We conclude that, like Temperature Compensation, Nutritional Compensation is a conserved circadian process in fungal and mammalian clocks and that it may share common molecular determinants.

Regulatory basis for reproductive flexibility in a meningitis-causing fungal pathogen

Pathogenic fungi of the genus Cryptococcus can undergo two sexual cycles, involving either bisexual diploidization (after fusion of haploid cells of different mating type) or unisexual diploidization (by autodiploidization of a single cell). Here, we construct a gene-deletion library for 111 transcription factor genes in Cryptococcus deneoformans, and explore the roles of these regulatory networks in the two reproductive modes. We show that transcription factors crucial for bisexual syngamy induce the expression of known mating determinants as well as other conserved genes of unknown function. Deletion of one of these genes, which we term FMP1, leads to defects in bisexual reproduction in C. deneoformans, its sister species Cryptococcus neoformans, and the ascomycete Neurospora crassa. Furthermore, we show that a recently evolved regulatory cascade mediates pre-meiotic unisexual autodiploidization, supporting that this reproductive process is a recent evolutionary innovation. Our findings indicate that genetic circuits with different evolutionary ages govern hallmark events distinguishing unisexual and bisexual reproduction in Cryptococcus.

Post-transcriptional control of fungal cell wall synthesis

Pathogenic fungi hide from their hosts by camouflage, obscuring immunogenic cell wall components such as beta-glucan with innocuous coverings such as mannoproteins and alpha-glucan that are less readily recognised by the host. Attempts to understand how such processes are regulated have met with varying success. Typically studies focus on understanding the transcriptional response of fungi to either their reservoir environment or the host. However, such approaches do not fully address this research question, due to the layers of post-transcriptional and post-translational regulation that occur within a cell. Although in animals the impact of post-transcriptional and post-translational regulation has been well characterised, our knowledge of these processes in the fungal kingdom is more limited. Mutations in RNA-binding proteins, like Ssd1 and Candida albicans Slr1, affect cell wall composition and fungal virulence indicating that post-transcriptional regulation plays a key role in these processes. Here, we review the current state of knowledge of fungal post-transcriptional regulation, and link this to potential mechanisms of immune evasion by drawing on studies from model yeast and plant pathogenic fungi. We highlight several RNA-binding proteins that regulate cell wall synthesis and could be involved in local translation of cell wall components. Expanding our knowledge on post-transcriptional regulation in human fungal pathogens is essential to fully comprehend fungal virulence strategies and for the design of novel antifungal therapies.

A comprehensive thermodynamic model for RNA binding by the Saccharomyces cerevisiae Pumilio protein PUF4

Genomic methods have been valuable for identifying RNA-binding proteins (RBPs) and the genes, pathways, and processes they regulate. Nevertheless, standard motif descriptions cannot be used to predict all RNA targets or test quantitative models for cellular interactions and regulation. We present a complete thermodynamic model for RNA binding to the S. cerevisiae Pumilio protein PUF4 derived from direct binding data for 6180 RNAs measured using the RNA on a massively parallel array (RNA-MaP) platform. The PUF4 model is highly similar to that of the related RBPs, human PUM2 and PUM1, with one marked exception: a single favorable site of base flipping for PUF4, such that PUF4 preferentially binds to a non-contiguous series of residues. These results are foundational for developing and testing cellular models of RNA-RBP interactions and function, for engineering RBPs, for understanding the biophysical nature of RBP binding and the evolutionary landscape of RNAs and RBPs.

PUMILIO competes with AUF1 to control DICER1 RNA levels and miRNA processing

DICER1 syndrome is a cancer pre-disposition disorder caused by mutations that disrupt the function of DICER1 in miRNA processing. Studying the molecular, cellular and oncogenic effects of these mutations can reveal novel mechanisms that control cell homeostasis and tumor biology. Here, we conduct the first analysis of pathogenic DICER1 syndrome allele from the DICER1 3'UTR. We find that the DICER1 syndrome allele, rs1252940486, abolishes interaction with the PUMILIO RNA binding protein with the DICER1 3'UTR, resulting in the degradation of the DICER1 mRNA by AUF1. This single mutational event leads to diminished DICER1 mRNA and protein levels, and widespread reprogramming of miRNA networks. The in-depth characterization of the rs1252940486 DICER1 allele, reveals important post-transcriptional regulatory events that control DICER1 levels.

Expanding the binding specificity for RNA recognition by a PUF domain

The ability to design a protein to bind specifically to a target RNA enables numerous applications, with the modular architecture of the PUF domain lending itself to new RNA-binding specificities. For each repeat of the Pumilio-1 PUF domain, we generate a library that contains the 8,000 possible combinations of amino acid substitutions at residues critical for RNA contact. We carry out yeast three-hybrid selections with each library against the RNA recognition sequence for Pumilio-1, with any possible base present at the position recognized by the randomized repeat. We use sequencing to score the binding of each variant, identifying many variants with highly repeat-specific interactions. From these data, we generate an RNA binding code specific to each repeat and base. We use this code to design PUF domains against 16 RNAs, and find that some of these domains recognize RNAs with two, three or four changes from the wild type sequence.

Investigation of RNA metabolism through large-scale genetic interaction profiling in yeast

Gene deletion and gene expression alteration can lead to growth defects that are amplified or reduced when a second mutation is present in the same cells. We performed 154 genetic interaction mapping (GIM) screens with query mutants related with RNA metabolism and estimated the growth rates of about 700 000 double mutant Saccharomyces cerevisiae strains. The tested targets included the gene deletion collection and 900 strains in which essential genes were affected by mRNA destabilization (DAmP). To analyze the results, we developed RECAP, a strategy that validates genetic interaction profiles by comparison with gene co-citation frequency, and identified links between 1471 genes and 117 biological processes. In addition to these large-scale results, we validated both enhancement and suppression of slow growth measured for specific RNA-related pathways. Thus, negative genetic interactions identified a role for the OCA inositol polyphosphate hydrolase complex in mRNA translation initiation. By analysis of suppressors, we found that Puf4, a Pumilio family RNA binding protein, inhibits ribosomal protein Rpl9 function, by acting on a conserved UGUAcauUA motif located downstream the stop codon of the RPL9B mRNA. Altogether, the results and their analysis should represent a useful resource for discovery of gene function in yeast.

The Pumilio RNA-binding protein APUM24 regulates seed maturation by fine-tuning the BPM-WRI1 module in Arabidopsis

Pumilio RNA-binding proteins participate in messenger RNA (mRNA) degradation and translational repression, but their roles in plant development are largely unclear. Here, we show that Arabidopsis PUMILIO PROTEIN24 (APUM24), an atypical Pumilio-homology domain-containing protein, plays an important part in regulating seed maturation, a major stage of plant development. APUM24 is strongly expressed in maturing seeds. Reducing APUM24 expression resulted in abnormal seed maturation, wrinkled seeds, and lower seed oil contents, and APUM24 knockdown resulted in lower levels of WRINKLED 1 (WRI1), a key transcription factor controlling seed oil accumulation, and lower expression of WRI1 target genes. APUM24 reduces the mRNA stability of BTB/POZMATH (BPM) family genes, thus decreasing BPM protein levels. BPM is responsible for the 26S proteasome-mediated degradation of WRI1 and has important functions in plant growth and development. The 3' untranslated regions of BPM family genes contain putative Pumilio response elements (PREs), which are bound by APUM24. Reduced BPM or increased WRI1 expression rescued the deficient seed maturation of apum24-2 knockdown mutants, and APUM24 overexpression resulted in increased seed size and weight. Therefore, APUM24 is crucial to seed maturation through its action as a positive regulator fine-tuning the BPM-WRI1 module, making APUM24 a promising target for breeding strategies to increase crop yields.

Phosphorylation of mRNA-Binding Proteins Puf1 and Puf2 by TORC2-Activated Protein Kinase Ypk1 Alleviates Their Repressive Effects

Members of the Puf family of RNA-binding proteins typically associate via their Pumilio homology domain with specific short motifs in the 3’-UTR of an mRNA and thereby influence the stability, localization and/or efficiency of translation of the bound transcript. In our prior unbiased proteome-wide screen for targets of the TORC2-stimulated protein kinase Ypk1, we identified the paralogs Puf1/Jsn1 and Puf2 as high-confidence substrates. Earlier work by others had demonstrated that Puf1 and Puf2 exhibit a marked preference for interaction with mRNAs encoding plasma membrane-associated proteins, consistent with our previous studies documenting that a primary physiological role of TORC2-Ypk1 signaling is maintenance of plasma membrane homeostasis. Here, we show, first, that both Puf1 and Puf2 are authentic Ypk1 substrates both in vitro and in vivo. Fluorescently tagged Puf1 localizes constitutively in cortical puncta closely apposed to the plasma membrane, whereas Puf2 does so in the absence of its Ypk1 phosphorylation, but is dispersed in the cytosol when phosphorylated. We further demonstrate that Ypk1-mediated phosphorylation of Puf1 and Puf2 upregulates production of the protein products of the transcripts to which they bind, with a concomitant increase in the level of the cognate mRNAs. Thus, Ypk1 phosphorylation relieves Puf1- and Puf2-mediated post-transcriptional repression mainly by counteracting their negative effect on transcript stability. Using a heterologous protein-RNA tethering and fluorescent protein reporter assay, the consequence of Ypk1 phosphorylation in vivo was recapitulated for full-length Puf1 and even for N-terminal fragments (residues 1-340 and 143-295) corresponding to the region upstream of its dimerization domain (an RNA-recognition motif fold) encompassing its two Ypk1 phosphorylation sites (both also conserved in Puf2). This latter result suggests that alleviation of Puf1-imposed transcript destabilization does not obligatorily require dissociation of Ypk1-phosphorylated Puf1 from a transcript. Our findings add new insight about how the TORC2-Ypk1 signaling axis regulates the content of plasma membrane-associated proteins to promote maintenance of the integrity of the cell envelope.

The GUL-1 Protein Binds Multiple RNAs Involved in Cell Wall Remodeling and Affects the MAK-1 Pathway in Neurospora crassa

The Neurospora crassa GUL-1 is part of the COT-1 pathway, which plays key roles in regulating polar hyphal growth and cell wall remodeling. We show that GUL-1 is a bona fide RNA-binding protein (RBP) that can associate with 828 "core" mRNA species. When cell wall integrity (CWI) is challenged, expression of over 25% of genomic RNA species are modulated (2,628 mRNAs, including the GUL-1 mRNA). GUL-1 binds mRNAs of genes related to translation, cell wall remodeling, circadian clock, endoplasmic reticulum (ER), as well as CWI and MAPK pathway components. GUL-1 interacts with over 100 different proteins, including stress-granule and P-body proteins, ER components and components of the MAPK, COT-1, and STRIPAK complexes. Several additional RBPs were also shown to physically interact with GUL-1. Under stress conditions, GUL-1 can localize to the ER and affect the CWI pathway-evident via altered phosphorylation levels of MAK-1, interaction with mak-1 transcript, and involvement in the expression level of the transcription factor adv-1. We conclude that GUL-1 functions in multiple cellular processes, including the regulation of cell wall remodeling, via a mechanism associated with the MAK-1 pathway and stress-response.

Translational remodeling by RNA-binding proteins and noncoding RNAs

Responsible for generating the proteome that controls phenotype, translation is the ultimate convergence point for myriad upstream signals that influence gene expression. System-wide adaptive translational reprogramming has recently emerged as a pillar of cellular adaptation. As classic regulators of mRNA stability and translation efficiency, foundational studies established the concept of collaboration and competition between RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) on individual mRNAs. Fresh conceptual innovations now highlight stress-activated, evolutionarily conserved RBP networks and ncRNAs that increase the translation efficiency of populations of transcripts encoding proteins that participate in a common cellular process. The discovery of post-transcriptional functions for long noncoding RNAs (lncRNAs) was particularly intriguing given their cell-type-specificity and historical definition as nuclear-functioning epigenetic regulators. The convergence of RBPs, lncRNAs, and microRNAs on functionally related mRNAs to enable adaptive protein synthesis is a newer biological paradigm that highlights their role as "translatome (protein output) remodelers" and reinvigorates the paradigm of "RNA operons." Together, these concepts modernize our understanding of cellular stress adaptation and strategies for therapeutic development. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications Translation > Translation Regulation Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Rpb4 and Puf3 imprint and post-transcriptionally control the stability of a common set of mRNAs in yeast

Gene expression involving RNA polymerase II is regulated by the concerted interplay between mRNA synthesis and degradation, crosstalk in which mRNA decay machinery and transcription machinery respectively impact transcription and mRNA stability. Rpb4, and likely dimer Rpb4/7, seem the central components of the RNA pol II governing these processes. In this work we unravel the molecular mechanisms participated by Rpb4 that mediate the posttranscriptional events regulating mRNA imprinting and stability. By RIP-Seq, we analysed genome-wide the association of Rpb4 with mRNAs and demonstrated that it targeted a large population of more than 1400 transcripts. A group of these mRNAs was also the target of the RNA binding protein, Puf3. We demonstrated that Rpb4 and Puf3 physically, genetically, and functionally interact and also affect mRNA stability, and likely the imprinting, of a common group of mRNAs. Furthermore, the Rpb4 and Puf3 association with mRNAs depends on one another. We also demonstrated, for the first time, that Puf3 associates with chromatin in an Rpb4-dependent manner. Our data also suggest that Rpb4 could be a key element of the RNA pol II that coordinates mRNA synthesis, imprinting and stability in cooperation with RBPs.

MS2-TRIBE Evaluates Both Protein-RNA Interactions and Nuclear Organization of Transcription by RNA Editing

Both UV-cross-linking and immunoprecipitation (CLIP) and RNA editing (TRIBE) can identify the targets of RNA-binding proteins. To evaluate false-positives of CLIP and TRIBE, endogenous β-actin mRNA was tagged with MS2 stem loops, making it the only bona fide target mRNA for the MS2 capsid protein (MCP). CLIP and TRIBE detected β-actin, albeit with false-positives. False-positive CLIP signals were attributed to nonspecific antibody interactions. In contrast, putative false-positive TRIBE targets were genes spatially proximal to the β-actin gene. MCP-ADAR edited nearby nascent transcripts consistent with interchromosomal contacts observed in Hi-C. The identification of nascent contacts implies RNA regulatory proteins (e.g., splicing factors) associated with multiple nascent transcripts, forming domains of post-transcriptional activity. Repeating these results with an integrated inducible MS2 reporter indicated that MS2-TRIBE can be applied to a broad array of cells and transcripts to study spatial organization and nuclear RNA regulation.

The evolution of the Puf superfamily of proteins across the tree of eukaryotes

BACKGROUND

Eukaryotic gene expression is controlled by a number of RNA-binding proteins (RBP), such as the proteins from the Puf (Pumilio and FBF) superfamily (PufSF). These proteins bind to RNA via multiple Puf repeat domains, each of which specifically recognizes a single RNA base. Recently, three diversified PufSF proteins have been described in model organisms, each of which is responsible for the maturation of ribosomal RNA or the translational regulation of mRNAs; however, less is known about the role of these proteins across eukaryotic diversity.

RESULTS

Here, we investigated the distribution and function of PufSF RBPs in the tree of eukaryotes. We determined that the following PufSF proteins are universally conserved across eukaryotes and can be broadly classified into three groups: (i) Nop9 orthologues, which participate in the nucleolar processing of immature 18S rRNA; (ii) 'classical' Pufs, which control the translation of mRNA; and (iii) PUM3 orthologues, which are involved in the maturation of 7S rRNA. In nearly all eukaryotes, the rRNA maturation proteins, Nop9 and PUM3, are retained as a single copy, while mRNA effectors ('classical' Pufs) underwent multiple lineage-specific expansions. We propose that the variation in number of 'classical' Pufs relates to the size of the transcriptome and thus the potential mRNA targets. We further distinguished full set of PufSF proteins in divergent metamonad Giardia intestinalis and initiated their cellular and biochemical characterization.

CONCLUSIONS

Our data suggest that the last eukaryotic common ancestor (LECA) already contained all three types of PufSF proteins and that 'classical' Pufs then underwent lineage-specific expansions.

Mitochondrial Biogenesis Is Positively Regulated by Casein Kinase I Hrr25 Through Phosphorylation of Puf3 in Saccharomyces cerevisiae

Mitochondrial biogenesis requires coordinated expression of genes encoding mitochondrial proteins, which in Saccharomyces cerevisiae is achieved in part via post-transcriptional control by the Pumilio RNA-binding domain protein Puf3 Puf3 binds to the 3'-UTR of many messenger RNAs (mRNAs) that encode mitochondrial proteins, regulating their turnover, translation, and/or mitochondrial targeting. Puf3 hyperphosphorylation correlates with increased mitochondrial biogenesis; however, the kinase responsible for Puf3 phosphorylation is unclear. Here, we show that the casein kinase I protein Hrr25 negatively regulates Puf3 by mediating its phosphorylation. An hrr25 mutation results in reduced phosphorylation of Puf3 in vivo and a puf3 deletion mutation reverses growth defects of hrr25 mutant cells grown on medium with a nonfermentable carbon source. We show that Hrr25 directly phosphorylates Puf3, and that the interaction between Puf3 and Hrr25 is mediated through the N-terminal domain of Puf3 and the kinase domain of Hrr25 We further found that an hrr25 mutation reduces GFP expression from GFP reporter constructs carrying the 3'-UTR of Puf3 targets. Downregulation of GFP expression due to an hrr25 mutation can be reversed either by puf3Δ or by mutations to the Puf3-binding sites in the 3'-UTR of the GFP reporter constructs. Together, our data indicate that Hrr25 is a positive regulator of mitochondrial biogenesis by phosphorylating Puf3 and inhibiting its function in downregulating target mRNAs encoding mitochondrial proteins.

Localization elements and zip codes in the intracellular transport and localization of messenger RNAs in Saccharomyces cerevisiae

Intracellular trafficking and localization of mRNAs provide a mechanism of regulation of expression of genes with excellent spatial control. mRNA localization followed by localized translation appears to be a mechanism of targeted protein sorting to a specific cell-compartment, which is linked to the establishment of cell polarity, cell asymmetry, embryonic axis determination, and neuronal plasticity in metazoans. However, the complexity of the mechanism and the components of mRNA localization in higher organisms prompted the use of the unicellular organism Saccharomyces cerevisiae as a simplified model organism to study this vital process. Current knowledge indicates that a variety of mRNAs are asymmetrically and selectively localized to the tip of the bud of the daughter cells, to the vicinity of endoplasmic reticulum, mitochondria, and nucleus in this organism, which are connected to diverse cellular processes. Interestingly, specific cis-acting RNA localization elements (LEs) or RNA zip codes play a crucial role in the localization and trafficking of these localized mRNAs by providing critical binding sites for the specific RNA-binding proteins (RBPs). In this review, we present a comprehensive account of mRNA localization in S. cerevisiae, various types of localization elements influencing the mRNA localization, and the RBPs, which bind to these LEs to implement a number of vital physiological processes. Finally, we emphasize the significance of this process by highlighting their connection to several neuropathological disorders and cancers. This article is categorized under: RNA Export and Localization > RNA Localization.

Human protein-RNA interaction network is highly stable across mammals

BACKGROUND

RNA-binding proteins (RBPs) are crucial in modulating RNA metabolism in eukaryotes thereby controlling an extensive network of RBP-RNA interactions. Although previous studies on the conservation of RBP targets have been carried out in lower eukaryotes such as yeast, relatively little is known about the extent of conservation of the binding sites of RBPs across mammalian species.

RESULTS

In this study, we employ CLIP-seq datasets for 60 human RBPs and demonstrate that most binding sites for a third of these RBPs are conserved in at least 50% of the studied vertebrate species. Across the studied RBPs, binding sites were found to exhibit a median conservation of 58%, ~ 20% higher than random genomic locations, suggesting a significantly higher preservation of RBP-RNA interaction networks across vertebrates. RBP binding sites were highly conserved across primates with weak conservation profiles in birds and fishes. We also note that phylogenetic relationship between members of an RBP family does not explain the extent of conservation of their binding sites across species. Multivariate analysis to uncover features contributing to differences in the extents of conservation of binding sites across RBPs revealed RBP expression level and number of post-transcriptional targets to be the most prominent factors. Examination of the location of binding sites at the gene level confirmed that binding sites occurring on the 3' region of a gene are highly conserved across species with 90% of the RBPs exhibiting a significantly higher conservation of binding sites in 3' regions of a gene than those occurring in the 5'. Gene set enrichment analysis on the extent of conservation of binding sites to identify significantly associated human phenotypes revealed an enrichment for multiple developmental abnormalities.

CONCLUSIONS

Our results suggest that binding sites of human RBPs are highly conserved across primates with weak conservation profiles in lower vertebrates and evolutionary relationship between members of an RBP family does not explain the extent of conservation of their binding sites. Expression level and number of targets of an RBP are important factors contributing to the differences in the extent of conservation of binding sites. RBP binding sites on 3' ends of a gene are the most conserved across species. Phenotypic analysis on the extent of conservation of binding sites revealed the importance of lineage-specific developmental events in post-transcriptional regulatory network evolution.

Distinct RNA-binding modules in a single PUF protein cooperate to determine RNA specificity

PUF proteins, named for Drosophila Pumilio (PUM) and Caenorhabditis elegans fem-3-binding factor (FBF), recognize specific sequences in the mRNAs they bind and control. RNA binding by classical PUF proteins is mediated by a characteristic PUM homology domain (PUM-HD). The Puf1 and Puf2 proteins possess a distinct architecture and comprise a highly conserved subfamily among fungal species. Puf1/Puf2 proteins contain two types of RNA-binding domain: a divergent PUM-HD and an RNA recognition motif (RRM). They recognize RNAs containing UAAU motifs, often in clusters. Here, we report a crystal structure of the PUM-HD of a fungal Puf1 in complex with a dual UAAU motif RNA. Each of the two UAAU tetranucleotides are bound by a Puf1 PUM-HD forming a 2:1 protein-to-RNA complex. We also determined crystal structures of the Puf1 RRM domain that identified a dimerization interface. The PUM-HD and RRM domains act in concert to determine RNA-binding specificity: the PUM-HD dictates binding to UAAU, and dimerization of the RRM domain favors binding to dual UAAU motifs rather than a single UAAU. Cooperative action of the RRM and PUM-HD identifies a new mechanism by which multiple RNA-binding modules in a single protein collaborate to create a unique RNA-binding specificity.

A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

Functional and molecular characterization of the conserved Arabidopsis PUMILIO protein, APUM9

Here we demonstrate that the APUM9 RNA-binding protein and its co-factors play a role in mRNA destabilization and how this activity might regulate early plant development. APUM9 is a conserved PUF RNA-binding protein (RBP) under complex transcriptional control mediated by a transposable element (TE) that restricts its expression in Arabidopsis. Currently, little is known about the functional and mechanistic details of the plant PUF regulatory system and the biological relevance of the TE-mediated repression of APUM9 in plant development and stress responses. By combining a range of transient assays, we show here, that APUM9 binding to target transcripts can trigger their rapid decay via its conserved C-terminal RNA-binding domain. APUM9 directly interacts with DCP2, the catalytic subunit of the decapping complex and DCP2 overexpression induces rapid decay of APUM9 targeted mRNAs. We show that APUM9 negatively regulates the expression of ABA signaling genes during seed imbibition, and thereby might contribute to the switch from dormant stage to seed germination. By contrast, strong TE-mediated repression of APUM9 is important for normal plant growth in the later developmental stages. Finally, APUM9 overexpression plants show slightly enhanced heat tolerance suggesting that TE-mediated control of APUM9, might have a role not only in embryonic development, but also in plant adaptation to heat stress conditions.

Evolutionary Transition of GAL Regulatory Circuit from Generalist to Specialist Function in Ascomycetes

The Gal4 transcription factor (TF) controls gene expression by binding the DNA sequence motif CGG(N₁₁)CCG. Well studied versions regulate metabolism of glucose in Candida albicans and galactose in Saccharomyces cerevisiae. Gal4 is also found within Aspergillus species and shows a wide range of potential binding targets. Members of the CTG clade that reassigned CUG codons from leucine to serine lack the Gal80 binding domain of Gal4, and they use the TF to regulate only glycolytic genes. In this clade, the galactose catabolic pathway (also known as the Leloir pathway) genes are regulated by Rtg1/Rtg3. In the WGD species, the complete Gal4/Gal80 module is limited to regulation of the Leloir pathway, while glycolysis is controlled by Gcr1/Gcr2. This shows a switch of Gal4 from a generalist to a specialist within the ascomycetes, and the split of glucose and galactose metabolism into distinct regulatory circuits.

Toward Identifying Subnetworks from FBF Binding Landscapes in Caenorhabditis Spermatogenic or Oogenic Germlines

Metazoan PUF (Pumilio and FBF) RNA-binding proteins regulate various biological processes, but a common theme across phylogeny is stem cell regulation. In Caenorhabditis elegans, FBF (fem-3 Binding Factor) maintains germline stem cells regardless of which gamete is made, but FBF also functions in the process of spermatogenesis. We have begun to "disentangle" these biological roles by asking which FBF targets are gamete-independent, as expected for stem cells, and which are gamete-specific. Specifically, we compared FBF iCLIP binding profiles in adults making sperm to those making oocytes. Normally, XX adults make oocytes. To generate XX adults making sperm, we used a fem-3(gf) mutant requiring growth at 25°; for comparison, wild-type oogenic hermaphrodites were also raised at 25°. Our FBF iCLIP data revealed FBF binding sites in 1522 RNAs from oogenic adults and 1704 RNAs from spermatogenic adults. More than half of these FBF targets were independent of germline gender. We next clustered RNAs by FBF-RNA complex frequencies and found four distinct blocks. Block I RNAs were enriched in spermatogenic germlines, and included validated target fog-3, while Block II and III RNAs were common to both genders, and Block IV RNAs were enriched in oogenic germlines. Block II (510 RNAs) included almost all validated FBF targets and was enriched for cell cycle regulators. Block III (21 RNAs) was enriched for RNA-binding proteins, including previously validated FBF targets gld-1 and htp-1 We suggest that Block I RNAs belong to the FBF network for spermatogenesis, and that Blocks II and III are associated with stem cell functions.

RNA-binding proteins distinguish between similar sequence motifs to promote targeted deadenylation by Ccr4-Not

The Ccr4-Not complex removes mRNA poly(A) tails to regulate eukaryotic mRNA stability and translation. RNA-binding proteins contribute to specificity by interacting with both Ccr4-Not and target mRNAs, but this is not fully understood. Here, we reconstitute accelerated and selective deadenylation of RNAs containing AU-rich elements (AREs) and Pumilio-response elements (PREs). We find that the fission yeast homologues of Tristetraprolin/TTP and Pumilio/Puf (Zfs1 and Puf3) interact with Ccr4-Not via multiple regions within low-complexity sequences, suggestive of a multipartite interface that extends beyond previously defined interactions. Using a two-color assay to simultaneously monitor poly(A) tail removal from different RNAs, we demonstrate that Puf3 can distinguish between RNAs of very similar sequence. Analysis of binding kinetics reveals that this is primarily due to differences in dissociation rate constants. Consequently, motif quality is a major determinant of mRNA stability for Puf3 targets in vivo and can be used for the prediction of mRNA targets.

When a cell needs to make a particular protein, it first copies the instructions from the matching gene into a molecule known as a messenger RNA (or an mRNA for short). The more mRNA copies it makes, the more protein it can produce. A simple way to control protein production is to raise or lower the number of these mRNA messages, and living cells have lots of ways to make this happen. One method involves codes built into the mRNAs themselves. The mRNAs can carry short sequences of genetic letters that can trigger their own destruction. Known as “destabilising motifs”, these sequences attract the attention of a group of proteins called Ccr4-Not. Together these proteins shorten the end of the mRNAs, preparing the molecules for degradation. But how does Ccr4-Not choose which mRNAs to target?

Different mRNAs carry different destabilising motifs. This means that when groups of mRNAs all carry the same motif, the cell can destroy them all together. This allows the cell to switch networks of related genes off together without affecting the mRNAs it still needs. What is puzzling is that the destabilising motifs that control different groups of mRNAs can be very similar, and scientists do not yet know how Ccr4-Not can tell the difference, or what triggers it to start breaking down groups of mRNAs. To find out, Webster et al. recreated the system in the laboratory using purified molecules.

The test-tube system confirmed previous suggestions that a protein called Puf3 forms a bridge between Ccr4-Not and mRNAs. It acts as a tether, recognising a destabilising motif and linking it to Ccr4-Not. Labelling different mRNAs with two colours of fluorescent dye showed how Puf3 helps the cell to choose which to destroy. Puf3 allows Ccr4-Not to select specific mRNAs from a mixture of molecules. Puf3 could distinguish between mRNAs that differed in a single letter of genetic code. When it matched with the wrong mRNA, it disconnected much faster than when it matched with the right one, preventing Ccr4-Not from linking up.

The ability to destroy specific mRNA messages is critical for cell survival. It happens when cells divide, during immune responses such as inflammation, and in early development. Understanding the targets of tethers like Puf3 could help scientists to predict which genes will switch off and when. This could reveal genes that work together, helping to unravel their roles inside cells.

Genetic basis for coordination of meiosis and sexual structure maturation in Cryptococcus neoformans

In the human fungal pathogen Cryptococcus neoformans, sex can benefit its pathogenicity through production of meiospores, which are believed to offer both physical and meiosis-created lineage advantages for its infections. Cryptococcus sporulation occurs following two parallel events, meiosis and differentiation of the basidium, the characteristic sexual structure of the basidiomycetes. However, the circuit integrating these events to ensure subsequent sporulation is unclear. Here, we show the spatiotemporal coordination of meiosis and basidial maturation by visualizing event-specific molecules in developing basidia defined by a quantitative approach. Monitoring of gene induction timing together with genetic analysis reveals co-regulation of the coordinated events by a shared regulatory program. Two RRM family regulators, Csa1 and Csa2, are crucial components that bridge meiosis and basidial maturation, further determining sporulation. We propose that the regulatory coordination of meiosis and basidial development serves as a determinant underlying the production of infectious meiospores in C. neoformans.

RNA-mediated gene regulation is less evolvable than transcriptional regulation

Cells regulate the activity of genes in a variety of ways. For example, they regulate transcription through DNA binding proteins called transcription factors, and they regulate mRNA stability and processing through RNA binding proteins. Based on current knowledge, transcriptional regulation is more widespread and is involved in many more evolutionary adaptations than posttranscriptional regulation. The reason could be that transcriptional regulation is studied more intensely. We suggest instead that transcriptional regulation harbors an intrinsic evolutionary advantage: when mutations change transcriptional regulation, they are more likely to bring forth novel patterns of such regulation. That is, transcriptional regulation is more evolvable. Our analysis suggests a reason why a specific kind of gene regulation is especially abundant in the living world.

Much of gene regulation is carried out by proteins that bind DNA or RNA molecules at specific sequences. One class of such proteins is transcription factors, which bind short DNA sequences to regulate transcription. Another class is RNA binding proteins, which bind short RNA sequences to regulate RNA maturation, transport, and stability. Here, we study the robustness and evolvability of these regulatory mechanisms. To this end, we use experimental binding data from 172 human and fruit fly transcription factors and RNA binding proteins as well as human polymorphism data to study the evolution of binding sites in vivo. We find little difference between the robustness of regulatory protein–RNA interactions and transcription factor–DNA interactions to DNA mutations. In contrast, we find that RNA-mediated regulation is less evolvable than transcriptional regulation, because mutations are less likely to create interactions of an RNA molecule with a new RNA binding protein than they are to create interactions of a gene regulatory region with a new transcription factor. Our observations are consistent with the high level of conservation observed for interactions between RNA binding proteins and their target molecules as well as the evolutionary plasticity of regulatory regions bound by transcription factors. They may help explain why transcriptional regulation is implicated in many more evolutionary adaptations and innovations than RNA-mediated gene regulation.

Multi-omics Reveal Specific Targets of the RNA-Binding Protein Puf3p and Its Orchestration of Mitochondrial Biogenesis

Coenzyme Q (CoQ) is a redox-active lipid required for mitochondrial oxidative phosphorylation (OxPhos). How CoQ biosynthesis is coordinated with the biogenesis of OxPhos protein complexes is unclear. Here, we show that the Saccharomyces cerevisiae RNA-binding protein (RBP) Puf3p regulates CoQ biosynthesis. To establish the mechanism for this regulation, we employed a multi-omic strategy to identify mRNAs that not only bind Puf3p but also are regulated by Puf3p in vivo. The CoQ biosynthesis enzyme Coq5p is a critical Puf3p target: Puf3p regulates the abundance of Coq5p and prevents its detrimental hyperaccumulation, thereby enabling efficient CoQ production. More broadly, Puf3p represses a specific set of proteins involved in mitochondrial protein import, translation, and OxPhos complex assembly (pathways essential to prime mitochondrial biogenesis). Our data reveal a mechanism for post-transcriptionally coordinating CoQ production with OxPhos biogenesis, and they demonstrate the power of multi-omics for defining genuine targets of RBPs.

High-resolution mapping of cis-regulatory variation in budding yeast

Genetic variants affecting gene-expression levels are a major source of phenotypic variation. The approximate locations of these variants can be mapped as expression quantitative trait loci (eQTLs); however, a major limitation of eQTLs is their low resolution, which precludes investigation of the causal variants and their molecular mechanisms. Here we report RNA-seq and full genome sequences for 85 diverse isolates of the yeast Saccharomyces cerevisiae-including wild, domesticated, and human clinical strains-which allowed us to perform eQTL mapping with 50-fold higher resolution than previously possible. In addition to variants in promoters, we uncovered an important role for variants in 3'UTRs, especially those affecting binding of the PUF family of RNA-binding proteins. The eQTLs are predominantly under negative selection, particularly those affecting essential genes and conserved genes. However, applying the sign test for lineage-specific selection revealed the polygenic up-regulation of dozens of biofilm suppressor genes in strains isolated from human patients, consistent with the key role of biofilms in fungal pathogenicity. In addition, a single variant in the promoter of a biofilm suppressor, NIT3, showed the strongest genome-wide association with clinical origin. Altogether, our results demonstrate the power of high-resolution eQTL mapping in understanding the molecular mechanisms of regulatory variation, as well as the natural selection acting on this variation that drives adaptation to environments, ranging from laboratories to vineyards to the human body.

RNA regulons in cancer and inflammation

Gene expression is the fundamental driving force that coordinates normal cellular processes and adapts to dysfunctional conditions such as oncogenic development and progression. While transcription is the basal process of gene expression, RNA transcripts are both the templates that encode proteins as well as perform functions that directly regulate diverse cellular processes. All levels of gene expression require coordination to optimize available resources, but how global gene expression drives cancers or responds to disrupting oncogenic mutations is not understood. Post-transcriptional coordination is controlled by RNA regulons that are governed by RNA-binding proteins (RBPs) and noncoding RNAs (ncRNAs) that bind and regulate multiple overlapping groups of functionally related RNAs. RNA regulons have been demonstrated to affect many biological functions and diseases, and many examples are known to regulate protein production in cancer and immune cells. In this review, we discuss RNA regulons demonstrated to coordinate global post-transcriptional mechanisms in carcinogenesis and inflammation.

The rewiring of transcription circuits in evolution

The binding of transcription regulators to cis-regulatory sequences is a key step through which all cells regulate expression of their genes. Due to gains and losses of cis-regulatory sequences and changes in the transcription regulators themselves, the binding connections between regulators and their target genes rapidly change over evolutionary time and constitute a major source of biological novelty. This review covers recent work, carried out in a wide range of species, that addresses the overall extent of these evolutionary changes, their consequences, and some of the molecular mechanisms that lie behind them.

Architecture and dynamics of overlapped RNA regulatory networks

A single protein can bind and regulate many mRNAs. Multiple proteins with similar specificities often bind and control overlapping sets of mRNAs. Yet little is known about the architecture or dynamics of overlapped networks. We focused on three proteins with similar structures and related RNA-binding specificities-Puf3p, Puf4p, and Puf5p of S. cerevisiae Using RNA Tagging, we identified a "super-network" comprised of four subnetworks: Puf3p, Puf4p, and Puf5p subnetworks, and one controlled by both Puf4p and Puf5p. The architecture of individual subnetworks, and thus the super-network, is determined by competition among particular PUF proteins to bind mRNAs, their affinities for binding elements, and the abundances of the proteins. The super-network responds dramatically: The remaining network can either expand or contract. These strikingly opposite outcomes are determined by an interplay between the relative abundance of the RNAs and proteins, and their affinities for one another. The diverse interplay between overlapping RNA-protein networks provides versatile opportunities for regulation and evolution.

Role of Mitochondrial Retrograde Pathway in Regulating Ethanol-Inducible Filamentous Growth in Yeast

In yeast, ethanol is produced as a by-product of fermentation through glycolysis. Ethanol also stimulates a developmental foraging response called filamentous growth and is thought to act as a quorum-sensing molecule. Ethanol-inducible filamentous growth was examined in a small collection of wine/European strains, which validated ethanol as an inducer of filamentous growth. Wine strains also showed variability in their filamentation responses, which illustrates the striking phenotypic differences that can occur among individuals. Ethanol-inducible filamentous growth in Σ1278b strains was independent of several of the major filamentation regulatory pathways [including fMAPK, RAS-cAMP, Snf1, Rpd3(L), and Rim101] but required the mitochondrial retrograde (RTG) pathway, an inter-organellar signaling pathway that controls the nuclear response to defects in mitochondrial function. The RTG pathway regulated ethanol-dependent filamentous growth by maintaining flux through the TCA cycle. The ethanol-dependent invasive growth response required the polarisome and transcriptional induction of the cell adhesion molecule Flo11p. Our results validate established stimuli that trigger filamentous growth and show how stimuli can trigger highly specific responses among individuals. Our results also connect an inter-organellar pathway to a quorum sensing response in fungi.

Recurrent rewiring and emergence of RNA regulatory networks

Alterations in regulatory networks contribute to evolutionary change. Transcriptional networks are reconfigured by changes in the binding specificity of transcription factors and their cognate sites. The evolution of RNA-protein regulatory networks is far less understood. The PUF (Pumilio and FBF) family of RNA regulatory proteins controls the translation, stability, and movements of hundreds of mRNAs in a single species. We probe the evolution of PUF-RNA networks by direct identification of the mRNAs bound to PUF proteins in budding and filamentous fungi and by computational analyses of orthologous RNAs from 62 fungal species. Our findings reveal that PUF proteins gain and lose mRNAs with related and emergent biological functions during evolution. We demonstrate at least two independent rewiring events for PUF3 orthologs, independent but convergent evolution of PUF4/5 binding specificity and the rewiring of the PUF4/5 regulons in different fungal lineages. These findings demonstrate plasticity in RNA regulatory networks and suggest ways in which their rewiring occurs.

Calcineurin Targets Involved in Stress Survival and Fungal Virulence

Calcineurin governs stress survival, sexual differentiation, and virulence of the human fungal pathogen Cryptococcus neoformans. Calcineurin is activated by increased Ca²⁺ levels caused by stress, and transduces signals by dephosphorylating protein substrates. Herein, we identified and characterized calcineurin substrates in C. neoformans by employing phosphoproteomic TiO₂ enrichment and quantitative mass spectrometry. The identified targets include the transactivator Crz1 as well as novel substrates whose functions are linked to P-bodies/stress granules (PBs/SGs) and mRNA translation and decay, such as Pbp1 and Puf4. We show that Crz1 is a bona fide calcineurin substrate, and Crz1 localization and transcriptional activity are controlled by calcineurin. We previously demonstrated that thermal and other stresses trigger calcineurin localization to PBs/SGs. Several calcineurin targets localized to PBs/SGs, including Puf4 and Pbp1, contribute to stress resistance and virulence individually or in conjunction with Crz1. Moreover, Pbp1 is also required for sexual development. Genetic epistasis analysis revealed that Crz1 and the novel targets Lhp1, Puf4, and Pbp1 function in a branched calcineurin pathway that orchestrates stress survival and virulence. These findings support a model whereby calcineurin controls stress and virulence, at the transcriptional level via Crz1, and post-transcriptionally by localizing to PBs/SGs and acting on targets involved in mRNA metabolism. The calcineurin targets identified in this study share little overlap with known calcineurin substrates, with the exception of Crz1. In particular, the mRNA binding proteins and PBs/SGs residents comprise a cohort of novel calcineurin targets that have not been previously linked to calcineurin in mammals or in Saccharomyces cerevisiae. This study suggests either extensive evolutionary rewiring of the calcineurin pathway, or alternatively that these novel calcineurin targets have yet to be characterized as calcineurin targets in other organisms. These findings further highlight C. neoformans as an outstanding model to define calcineurin-responsive virulence networks as targets for antifungal therapy.

Calcineurin is a Ca²⁺/calmodulin-dependent protein phosphatase essential for stress survival, sexual development, and virulence of the human fungal pathogen Cryptococcus neoformans and other major pathogenic fungi of global human health relevance. However, no calcineurin substrates are known in pathogenic fungi. Employing state-of-the-art phosphoproteomic approaches we identified calcineurin substrates, including calcineurin itself and the conserved Crz1 transcriptional activator known to function in calcium signaling and stress survival. Remarkably, our study also identified novel calcineurin targets involved in RNA processing, stability, and translation, which colocalize together with calcineurin in stress granules/P-bodies upon thermal stress. These findings support a model whereby calcineurin functions in a branched pathway, via Crz1 and several of the identified novel targets, that governs transcriptional and posttranscriptional circuits to drive stress survival, sexual development, and fungal virulence. Our study underscores C. neoformans as an experimental model to define basic paradigms of calcineurin signaling in global thermostress responsive virulence networks that can be targeted for fungal therapy.

The PUF binding landscape in metazoan germ cells

PUF (Pumilio/FBF) proteins are RNA-binding proteins and conserved stem cell regulators. The Caenorhabditis elegans PUF proteins FBF-1 and FBF-2 (collectively FBF) regulate mRNAs in germ cells. Without FBF, adult germlines lose all stem cells. A major gap in our understanding of PUF proteins, including FBF, is a global view of their binding sites in their native context (i.e., their "binding landscape"). To understand the interactions underlying FBF function, we used iCLIP (individual-nucleotide resolution UV crosslinking and immunoprecipitation) to determine binding landscapes of C. elegans FBF-1 and FBF-2 in the germline tissue of intact animals. Multiple iCLIP peak-calling methods were compared to maximize identification of both established FBF binding sites and positive control target mRNAs in our iCLIP data. We discovered that FBF-1 and FBF-2 bind to RNAs through canonical as well as alternate motifs. We also analyzed crosslinking-induced mutations to map binding sites precisely and to identify key nucleotides that may be critical for FBF-RNA interactions. FBF-1 and FBF-2 can bind sites in the 5'UTR, coding region, or 3'UTR, but have a strong bias for the 3' end of transcripts. FBF-1 and FBF-2 have strongly overlapping target profiles, including mRNAs and noncoding RNAs. From a statistically robust list of 1404 common FBF targets, 847 were previously unknown, 154 were related to cell cycle regulation, three were lincRNAs, and 335 were shared with the human PUF protein PUM2.

The Role of microRNAs in the Repeated Parallel Diversification of Lineages of Midas Cichlid Fish from Nicaragua

Cichlid fishes are an ideal model system for studying biological diversification because they provide textbook examples of rapid speciation. To date, there has been little focus on the role of gene regulation during cichlid speciation. However, in recent years, gene regulation has been recognized as a powerful force linking diversification in gene function to speciation. Here, we investigated the potential role of miRNA regulation in the diversification of six cichlid species of the Midas cichlid lineage (Amphilophus spp.) inhabiting the Nicaraguan crater lakes. Using several genomic resources, we inferred 236 Midas miRNA genes that were used to predict the miRNA target sites on 8,232 Midas 3′-UTRs. Using population genomic calculations of SNP diversity, we found the miRNA genes to be more conserved than protein coding genes. In contrast to what has been observed in other cichlid fish, but similar to what has been typically found in other groups, we observed genomic signatures of purifying selection on the miRNA targets by comparing these sites with the less conserved nontarget portion of the 3′-UTRs. However, in one species pair that has putatively speciated sympatrically in crater Lake Apoyo, we recovered a different pattern of relaxed purifying selection and high genetic divergence at miRNA targets. Our results suggest that sequence evolution at miRNA binding sites could be a critical genomic mechanism contributing to the rapid phenotypic evolution of Midas cichlids.

Post-transcriptional gene regulation in the biology and virulence of Candida albicans

In the human fungal pathogen Candida albicans, remodelling of gene expression drives host adaptation and virulence. Recent studies revealed that in addition to transcription, post‐transcriptional mRNA control plays important roles in virulence‐related pathways. Hyphal morphogenesis, biofilm formation, stress responses, antifungal drug susceptibility and virulence in animal models require post‐transcriptional regulators. This includes RNA binding proteins that control mRNA localization, decay and translation, as well as the cytoplasmic mRNA decay pathway. Comprehensive understanding of how modulation of gene expression networks drives C. albicans virulence will necessitate integration of our knowledge on transcriptional and post‐transcriptional mRNA control.