Evolution of yeast noncoding RNAs reveals an alternative mechanism for widespread intron loss

Inter-species association mapping links splice site evolution to METTL16 and SNRNP27K

Eukaryotic genes are interrupted by introns that are removed from transcribed RNAs by splicing. Patterns of splicing complexity differ between species, but it is unclear how these differences arise. We used inter-species association mapping with Saccharomycotina species to correlate splicing signal phenotypes with the presence or absence of splicing factors. Here, we show that variation in 5' splice site sequence preferences correlate with the presence of the U6 snRNA N6-methyladenosine methyltransferase METTL16 and the splicing factor SNRNP27K. The greatest variation in 5' splice site sequence occurred at the +4 position and involved a preference switch between adenosine and uridine. Loss of METTL16 and SNRNP27K orthologs, or a single SNRNP27K methionine residue, was associated with a preference for +4 U. These findings are consistent with splicing analyses of mutants defective in either METTL16 or SNRNP27K orthologs and models derived from spliceosome structures, demonstrating that inter-species association mapping is a powerful orthogonal approach to molecular studies. We identified variation between species in the occurrence of two major classes of 5' splice sites, defined by distinct interaction potentials with U5 and U6 snRNAs, that correlates with intron number. We conclude that variation in concerted processes of 5' splice site selection by U6 snRNA is associated with evolutionary changes in splicing signal phenotypes.

Interactions between Macrophages and Mast Cells in the Female Reproductive System

Mast cells (MCs) and macrophages (Mϕs) are innate immune cells that differentiate from early common myeloid precursors and reside in all body tissues. MCs have a unique capacity to neutralize/degrade toxic proteins, and they are hypothesized as being able to adopt two alternative polarization profiles, similar to Mϕs, with distinct or even opposite roles. Mϕs are very plastic phagocytic cells that are devoted to the elimination of senescent/anomalous endogenous entities (to maintain tissue homeostasis), and to the recognition and elimination of exogenous threats. They can adopt several functional phenotypes in response to microenvironmental cues, whose extreme profiles are the inflammatory/killing phenotype (M1) and the anti-inflammatory/healing phenotype (M2). The concomitant and abundant presence of these two cell types and the partial overlap of their defensive and homeostatic functions leads to the hypothesis that their crosstalk is necessary for the optimal coordination of their functions, both under physiological and pathological conditions. This review will examine the relationship between MCs and Mϕs in some situations of homeostatic regulation (menstrual cycle, embryo implantation), and in some inflammatory conditions in the same organs (endometriosis, preeclampsia), in order to appreciate the importance of their cross-regulation.

Evolution of intron splicing towards optimized gene expression is based on various Cis- and Trans-molecular mechanisms

Splicing expands, reshapes, and regulates the transcriptome of eukaryotic organisms. Despite its importance, key questions remain unanswered, including the following: Can splicing evolve when organisms adapt to new challenges? How does evolution optimize inefficiency of introns’ splicing and of the splicing machinery? To explore these questions, we evolved yeast cells that were engineered to contain an inefficiently spliced intron inside a gene whose protein product was under selection for an increased expression level. We identified a combination of mutations in Cis (within the gene of interest) and in Trans (in mRNA-maturation machinery). Surprisingly, the mutations in Cis resided outside of known intronic functional sites and improved the intron’s splicing efficiency potentially by easing tight mRNA structures. One of these mutations hampered a protein’s domain that was not under selection, demonstrating the evolutionary flexibility of multi-domain proteins as one domain functionality was improved at the expense of the other domain. The Trans adaptations resided in two proteins, Npl3 and Gbp2, that bind pre-mRNAs and are central to their maturation. Interestingly, these mutations either increased or decreased the affinity of these proteins to mRNA, presumably allowing faster spliceosome recruitment or increased time before degradation of the pre-mRNAs, respectively. Altogether, our work reveals various mechanistic pathways toward optimizations of intron splicing to ultimately adapt gene expression patterns to novel demands.

An experimental evolution study involving an inefficiently spliced intron reveals that the splicing machinery, introns, and RNA quality control factors evolve in Cis and in Trans when cells optimize their transcriptome to new challenges.

Exon junction complex components Y14 and Mago still play a role in budding yeast

Since their divergence from Pezizomycotina, the mRNA metabolism of budding yeasts have undergone regressive evolution. With the dramatic loss of introns, a number of quality control mechanisms have been simplified or lost during evolution, such as the exon junction complex (EJC). We report the identification of the core EJC components, Mago, Y14, and eIF4A3, in at least seven Saccharomycotina species, including Yarrowia lipolytica. Peripheral factors that join EJC, either to mediate its assembly (Ibp160 or Cwc22), or trigger downstream processes, are present in the same species, forming an evolutionary package. Co-immunoprecipitation studies in Y. lipolytica showed that Mago and Y14 have retained the capacity to form heterodimers, which successively bind to the peripheral factors Upf3, Aly/REF, and Pym. Phenotypes and RNA-Seq analysis of EJC mutants showed evidence of Y14 and Mago involvement in mRNA metabolism. Differences in unspliced mRNA levels suggest that Y14 binding either interferes with pre-mRNA splicing or retains mRNA in the nucleus before their export and translation. These findings indicate that yeast could be a relevant model for understanding EJC function.

Sense-antisense gene overlap is probably a cause for retaining the few introns in Giardia genome and the implications

Background

It is widely accepted that the last eukaryotic common ancestor and early eukaryotes were intron-rich and intron loss dominated subsequent evolution, thus the presence of only very few introns in some modern eukaryotes must be the consequence of massive loss. But it is striking that few eukaryotes were found to have completely lost introns. Despite extensive research, the causes of massive intron losses remain elusive. Actually the reverse question -- how the few introns can be retained under the evolutionary selection pressure of intron loss -- is equally significant but was rarely studied, except that it was conjectured that the essential functions of some introns prevent their loss. The situation that extremely few (eight) spliceosome-mediated cis-spliced introns present in the relatively simple genome of Giardia lamblia provides an excellent opportunity to explore this question.

Results

Our investigation found three types of distribution patterns of the few introns in the intron-containing genes: ancient intron in ancient gene, later-evolved intron in ancient gene, and later-evolved intron in later-evolved gene, which can reflect to some extent the dynamic evolution of introns in Giardia. Without finding any special features or functional importance of these introns responsible for their retention, we noticed and experimentally verified that some intron-containing genes form sense-antisense gene pairs with transcribable genes on their complementary strands, and that the introns just reside in the overlapping regions.

Conclusions

In Giardia’s evolution, despite constant evolutionary selection pressure of intron loss, intron gain can still occur in both ancient and later-evolved genes, but only a few introns are retained; at least the evolutionary retention of some of the introns might not be due to the functional constraint of the introns themselves but the causes outside of introns, such as the constraints imposed by other genomic functional elements overlapping with the introns. These findings can not only provide some clues to find new genomic functional elements -- in the areas overlapping with introns, but suggest that “functional constraint” of introns may not be necessarily directly associated with intron loss and gain, and that the real functions are probably still outside of our current knowledge.

Reviewers

This article was reviewed by Mikhail Gelfand, Michael Gray, and Igor Rogozin.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0226-5) contains supplementary material, which is available to authorized users.

The fungal snoRNAome

Small nucleolar RNAs (snoRNAs) are essential players in the rRNA biogenesis due to their involvement in the nucleolytic processing of the precursor and the subsequent guidance of nucleoside modifications. Within the kingdom Fungi, merely a few species-specific surveys have explored their snoRNA repertoire. However, the wide range of the snoRNA landscape spanning all major fungal lineages has not been mapped so far, mainly because of missing tools for automatized snoRNA detection and functional analysis. For the first time, we report here a comprehensive inventory of fungal snoRNAs together with a functional analysis and an in-depth investigation of their evolutionary history including innovations, deletions, and target switches. This large-scale analysis, incorporating more than 120 snoRNA families with more than 7700 individual snoRNA sequences, catalogs and clarifies the landscape of fungal snoRNA families, assigns functions to previously orphan snoRNAs, and increases the number of sequences by 450%. We also show that the snoRNAome is subject to ongoing rearrangements and adaptations, e.g., through lineage-specific targets and redundant guiding functions.

Protein coding genes as hosts for noncoding RNA expression

With the emergence of high-throughput sequence characterization methods and the subsequent improvements in gene annotations, it is becoming increasingly clear that a large proportion of eukaryotic protein-coding genes (as many as 50% in human) serve as host genes for non-coding RNA genes. Amongst the most extensively characterized embedded non-coding RNA genes, small nucleolar RNAs and microRNAs represent abundant families. Encoded individually or clustered, in sense or antisense orientation with respect to their host and independently expressed or dependent on host expression, the genomic characteristics of embedded genes determine their biogenesis and the extent of their relationship with their host gene. Not only can host genes and the embedded genes they harbour be co-regulated and mutually modulate each other, many are functionally coupled playing a role in the same cellular pathways. And while host-non-coding RNA relationships can be highly conserved, mechanisms have been identified, and in particular an association with transposable elements, allowing the appearance of copies of non-coding genes nested in host genes, or the migration of embedded genes from one host gene to another. The study of embedded non-coding genes and their relationship with their host genes increases the complexity of cellular networks and provides important new regulatory links that are essential to properly understand cell function.

Phylogenetic distribution of plant snoRNA families

Background

Small nucleolar RNAs (snoRNAs) are one of the most ancient families amongst non-protein-coding RNAs. They are ubiquitous in Archaea and Eukarya but absent in bacteria. Their main function is to target chemical modifications of ribosomal RNAs. They fall into two classes, box C/D snoRNAs and box H/ACA snoRNAs, which are clearly distinguished by conserved sequence motifs and the type of chemical modification that they govern. Similarly to microRNAs, snoRNAs appear in distinct families of homologs that affect homologous targets. In animals, snoRNAs and their evolution have been studied in much detail. In plants, however, their evolution has attracted comparably little attention.

Results

In order to chart the phylogenetic distribution of individual snoRNA families in plants, we applied a sophisticated approach for identifying homologs of known plant snoRNAs across the plant kingdom. In response to the relatively fast evolution of snoRNAs, information on conserved sequence boxes, target sequences, and secondary structure is combined to identify additional snoRNAs. We identified 296 families of snoRNAs in 24 species and traced their evolution throughout the plant kingdom. Many of the plant snoRNA families comprise paralogs. We also found that targets are well-conserved for most snoRNA families.

Conclusions

The sequence conservation of snoRNAs is sufficient to establish homologies between phyla. The degree of this conservation tapers off, however, between land plants and algae. Plant snoRNAs are frequently organized in highly conserved spatial clusters. As a resource for further investigations we provide carefully curated and annotated alignments for each snoRNA family under investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3301-2) contains supplementary material, which is available to authorized users.

Identification of Non-Coding RNAs in the Candida parapsilosis Species Group

The Candida CTG clade is a monophyletic group of fungal species that translates CTG as serine, and includes the pathogens Candida albicans and Candida parapsilosis. Research has typically focused on identifying protein-coding genes in these species. Here, we use bioinformatic and experimental approaches to annotate known classes of non-coding RNAs in three CTG-clade species, Candida parapsilosis, Candida orthopsilosis and Lodderomyces elongisporus. We also update the annotation of ncRNAs in the C. albicans genome. The majority of ncRNAs identified were snoRNAs. Approximately 50% of snoRNAs (including most of the C/D box class) are encoded in introns. Most are within mono- and polycistronic transcripts with no protein coding potential. Five polycistronic clusters of snoRNAs are highly conserved in fungi. In polycistronic regions, splicing occurs via the classical pathway, as well as by nested and recursive splicing. We identified spliceosomal small nuclear RNAs, the telomerase RNA component, signal recognition particle, RNase P RNA component and the related RNase MRP RNA component in all three genomes. Stem loop IV of the U2 spliceosomal RNA and the associated binding proteins were lost from the ancestor of C. parapsilosis and C. orthopsilosis, following the divergence from L. elongisporus. The RNA component of the MRP is longer in C. parapsilosis, C. orthopsilosis and L. elongisporus than in S. cerevisiae, but is substantially shorter than in C. albicans.

Novel Intronic RNA Structures Contribute to Maintenance of Phenotype in Saccharomyces cerevisiae

The Saccharomyces cerevisiae genome has undergone extensive intron loss during its evolutionary history. It has been suggested that the few remaining introns (in only 5% of protein-coding genes) are retained because of their impact on function under stress conditions. Here, we explore the possibility that novel noncoding RNA structures (ncRNAs) are embedded within intronic sequences and are contributing to phenotype and intron retention in yeast. We employed de novo RNA structure prediction tools to screen intronic sequences in S. cerevisiae and 36 other fungi. We identified and validated 19 new intronic RNAs via RNA sequencing (RNA-seq) and RT-PCR. Contrary to the common belief that excised introns are rapidly degraded, we found that, in six cases, the excised introns were maintained intact in the cells. In another two cases we showed that the ncRNAs were further processed from their introns. RNA-seq analysis confirmed that introns in ribosomal protein genes are more highly expressed when they contain predicted RNA structures. We deleted the novel intronic RNA structure within the GLC7 intron and showed that this region, rather than the intron itself, is responsible for the cell’s ability to respond to salt stress. We also showed a direct association between the in cis presence of the intronic RNA and GLC7 expression. Overall, these data support the notion that some introns may have been maintained in the genome because they harbor functional RNA structures.

FgPrp4 Kinase Is Important for Spliceosome B-Complex Activation and Splicing Efficiency in Fusarium graminearum

PRP4 encodes the only kinase among the spliceosome components. Although it is an essential gene in the fission yeast and other eukaryotic organisms, the Fgprp4 mutant was viable in the wheat scab fungus Fusarium graminearum. Deletion of FgPRP4 did not block intron splicing but affected intron splicing efficiency in over 60% of the F. graminearum genes. The Fgprp4 mutant had severe growth defects and produced spontaneous suppressors that were recovered in growth rate. Suppressor mutations were identified in the PRP6, PRP31, BRR2, and PRP8 orthologs in nine suppressor strains by sequencing analysis with candidate tri-snRNP component genes. The Q86K mutation in FgMSL1 was identified by whole genome sequencing in suppressor mutant S3. Whereas two of the suppressor mutations in FgBrr2 and FgPrp8 were similar to those characterized in their orthologs in yeasts, suppressor mutations in Prp6 and Prp31 orthologs or FgMSL1 have not been reported. Interestingly, four and two suppressor mutations identified in FgPrp6 and FgPrp31, respectively, all are near the conserved Prp4-phosphorylation sites, suggesting that these mutations may have similar effects with phosphorylation by Prp4 kinase. In FgPrp31, the non-sense mutation at R464 resulted in the truncation of the C-terminal 130 aa region that contains all the conserved Prp4-phosphorylation sites. Deletion analysis showed that the N-terminal 310-aa rich in SR residues plays a critical role in the localization and functions of FgPrp4. We also conducted phosphoproteomics analysis with FgPrp4 and identified S289 as the phosphorylation site that is essential for its functions. These results indicated that FgPrp4 is critical for splicing efficiency but not essential for intron splicing, and FgPrp4 may regulate pre-mRNA splicing by phosphorylation of other components of the tri-snRNP although itself may be activated by phosphorylation at S289.

Evaluation of the mechanisms of intron loss and gain in the social amoebae Dictyostelium

Background

Spliceosomal introns are a common feature of eukaryotic genomes. To approach a comprehensive understanding of intron evolution on Earth, studies should look beyond repeatedly studied groups such as animals, plants, and fungi. The slime mold Dictyostelium belongs to a supergroup of eukaryotes not covered in previous studies.

Results

We found 441 precise intron losses in Dictyostelium discoideum and 202 precise intron losses in Dictyostelium purpureum. Consistent with these observations, Dictyostelium discoideum was found to have significantly more copies of reverse transcriptase genes than Dictyostelium purpureum. We also found that the lost introns are significantly further from the 5′ end of genes than the conserved introns. Adjacent introns were prone to be lost simultaneously in Dictyostelium discoideum. In both Dictyostelium species, the exonic sequences flanking lost introns were found to have a significantly higher GC content than those flanking conserved introns. Together, these observations support a reverse-transcription model of intron loss in which intron losses were caused by gene conversion between genomic DNA and cDNA reverse transcribed from mature mRNA. We also identified two imprecise intron losses in Dictyostelium discoideum that may have resulted from genomic deletions. Ninety-eight putative intron gains were also observed. Consistent with previous studies of other lineages, the source sequences were found in only a small number of cases, with only two instances of intron gain identified in Dictyostelium discoideum.

Conclusions

Although they diverged very early from animals and fungi, Dictyostelium species have similar mechanisms of intron loss.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0567-y) contains supplementary material, which is available to authorized users.

Intron evolution in Saccharomycetaceae

Introns in protein-coding genes are very rare in hemiascomycetous yeast genomes. It has been suggested that these species have experienced extensive intron loss during their evolution from the postulated intron-rich fungal ancestor. However, no intron-devoid yeast species have been identified and some of the introns remaining within the genomes of intron-poor species, such as Saccharomyces cerevisiae, appear to be beneficial during growth under stress conditions. In order to reveal the pattern of intron retention within intron-poor yeast species and better understand the mechanisms of intron evolution, we generated a comprehensive set of 250 orthologous introns in the 20 species that comprise the Saccharomycetaceae, by analyzing RNA deep-sequencing data and alignments of intron-containing genes. Analysis of these intron sets shows that intron loss is at least two orders of magnitude more frequent than intron gain. Fine mapping of intron positions shows that intron sliding is rare, and that introns are almost always removed without changing the primary sequence of the encoded protein. The latter finding is consistent with the prevailing view that homologous recombination between reverse-transcribed mature mRNAs and the corresponding genomic locus is the primary mechanism of intron loss. However, we also find evidence that loss of a small number of introns is mediated by micro-homology, and that the number of intron losses is diminished in yeast species that have lost the microhomology end joining and nonhomologous end joining machinery.

Imprecise intron losses are less frequent than precise intron losses but are not rare in plants

Abstract

In this study, we identified 19 intron losses, including 11 precise intron losses (PILs), six imprecise intron losses (IILs), one de-exonization, and one exon deletion in tomato and potato, and 17 IILs in Arabidopsis thaliana. Comparative analysis of related genomes confirmed that all of the IILs have been fixed during evolution. Consistent with previous studies, our results indicate that PILs are a major type of intron loss. However, at least in plants, IILs are unlikely to be as rare as previously reported.

Reviewers

This article was reviewed by Jun Yu and Zhang Zhang. For complete reviews, see the Reviewers’ Reports section.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0056-7) contains supplementary material, which is available to authorized users.

smFRET studies of the 'encounter' complexes and subsequent intermediate states that regulate the selectivity of ligand binding

The selectivity with which a biomolecule can bind its cognate ligand when confronted by the vast array of structurally similar, competing ligands that are present in the cell underlies the fidelity of some of the most fundamental processes in biology. Because they collectively comprise one of only a few methods that can sensitively detect the 'encounter' complexes and subsequent intermediate states that regulate the selectivity of ligand binding, single-molecule fluorescence, and particularly single-molecule fluorescence resonance energy transfer (smFRET), approaches have revolutionized studies of ligand-binding reactions. Here, we describe a widely used smFRET strategy that enables investigations of a large variety of ligand-binding reactions, and discuss two such reactions, aminoacyl-tRNA selection during translation elongation and splice site selection during spliceosome assembly, that highlight both the successes and challenges of smFRET studies of ligand-binding reactions. We conclude by reviewing a number of emerging experimental and computational approaches that are expanding the capabilities of smFRET approaches for studies of ligand-binding reactions and that promise to reveal the mechanisms that control the selectivity of ligand binding with unprecedented resolution.

Reverse transcriptase and intron number evolution

BACKGROUND

Introns are universal in eukaryotic genomes and play important roles in transcriptional regulation, mRNA export to the cytoplasm, nonsense-mediated decay as both a regulatory and a splicing quality control mechanism, R-loop avoidance, alternative splicing, chromatin structure, and evolution by exon-shuffling.

METHODS

Sixteen complete fungal genomes were used 13 of which were sequenced and annotated by JGI. Ustilago maydis, Cryptococcus neoformans, and Coprinus cinereus (also named Coprinopsis cinerea) were from the Broad Institute. Gene models from JGI-annotated genomes were taken from the GeneCatalog track that contained the best representative gene models. Varying fractions of the GeneCatalog were manually curated by external users. For clarity, we used the JGI unique database identifier.

RESULTS

The last common ancestor of eukaryotes (LECA) has an estimated 6.4 coding exons per gene (EPG) and evolved into the diverse eukaryotic life forms, which is recapitulated by the development of a stem cell. We found a parallel between the simulated reverse transcriptase (RT)-mediated intron loss and the comparative analysis of 16 fungal genomes that spanned a wide range of intron density. Although footprints of RT (RTF) were dynamic, relative intron location (RIL) to the 5'-end of mRNA faithfully traced RT-mediated intron loss and revealed 7.7 EPG for LECA. The mode of exon length distribution was conserved in simulated intron loss, which was exemplified by the shared mode of 75 nt between fungal and Chlamydomonas genomes. The dominant ancient exon length was corroborated by the average exon length of the most intron-rich genes in fungal genomes and consistent with ancient protein modules being ~25 aa. Combined with the conservation of a protein length of 400 aa, the earliest ancestor of eukaryotes could have 16 EPG. During earlier evolution, Ascomycota's ancestor had significantly more 3'-biased RT-mediated intron loss that was followed by dramatic RTF loss. There was a down trend of EPG from more conserved to less conserved genes. Moreover, species-specific genes have higher exon-densities, shorter exons, and longer introns when compared to genes conserved at the phylum level. However, intron length in species-specific genes became shorter than that of genes conserved in all species after genomes experiencing drastic intron loss. The estimated EPG from the most frequent exon length is more than double that from the RIL method.

CONCLUSIONS

This implies significant intron loss during the very early period of eukaryotic evolution. De novo gene-birth contributes to shorter exons, longer introns, and higher exon-density in species-specific genes relative to conserved genes.

Conservation and divergence of transcriptional coregulations between box C/D snoRNA and ribosomal protein genes in Ascomycota

Coordinated assembly of the ribosome is essential for proper translational activity in eukaryotic cells. It is therefore critical to coordinate the expression of components of ribosomal programs with the cell's nutritional status. However, coordinating expression of these components is poorly understood. Here, by combining experimental and computational approaches, we systematically identified box C/D snoRNAs in four fission yeasts and found that the expression of box C/D snoRNA and ribosomal protein (RP) genes were orchestrated by a common Homol-D box, thereby ensuring a constant balance of these two genetic components. Interestingly, such transcriptional coregulations could be observed in most Ascomycota species and were mediated by different cis-regulatory elements. Via the reservation of cis elements, changes in spatial configuration, the substitution of cis elements, and gain or loss of cis elements, the regulatory networks of box C/D snoRNAs evolved to correspond with those of the RP genes, maintaining transcriptional coregulation between box C/D snoRNAs and RP genes. Our results indicate that coregulation via common cis elements is an important mechanism to coordinate expression of the RP and snoRNA genes, which ensures a constant balance of these two components.

Species and condition specific adaptation of the transcriptional landscapes in Candida albicans and Candida dubliniensis

BACKGROUND

Although Candida albicans and Candida dubliniensis are most closely related, both species behave significantly different with respect to morphogenesis and virulence. In order to gain further insight into the divergent routes for morphogenetic adaptation in both species, we investigated qualitative along with quantitative differences in the transcriptomes of both organisms by cDNA deep sequencing.

RESULTS

Following genome-associated assembly of sequence reads we were able to generate experimentally verified databases containing 6016 and 5972 genes for C. albicans and C. dubliniensis, respectively. About 95% of the transcriptionally active regions (TARs) contain open reading frames while the remaining TARs most likely represent non-coding RNAs. Comparison of our annotations with publically available gene models for C. albicans and C. dubliniensis confirmed approximately 95% of already predicted genes, but also revealed so far unknown novel TARs in both species. Qualitative cross-species analysis of these databases revealed in addition to 5802 orthologs also 399 and 49 species-specific protein coding genes for C. albicans and C. dubliniensis, respectively. Furthermore, quantitative transcriptional profiling using RNA-Seq revealed significant differences in the expression of orthologs across both species. We defined a core subset of 84 hyphal-specific genes required for both species, as well as a set of 42 genes that seem to be specifically induced during hyphal morphogenesis in C. albicans.

CONCLUSIONS

Species-specific adaptation in C. albicans and C. dubliniensis is governed by individual genetic repertoires but also by altered regulation of conserved orthologs on the transcriptional level.

Role of Ess1 in growth, morphogenetic switching, and RNA polymerase II transcription in Candida albicans

Candida albicans is a fungal pathogen that causes potentially fatal infections among immune-compromised individuals. The emergence of drug resistant C. albicans strains makes it important to identify new antifungal drug targets. Among potential targets are enzymes known as peptidyl-prolyl cis/trans isomerases (PPIases) that catalyze isomerization of peptide bonds preceding proline. We are investigating a PPIase called Ess1, which is conserved in all major human pathogenic fungi. Previously, we reported that C. albicans Ess1 is essential for growth and morphogenetic switching. In the present study, we re-evaluated these findings using more rigorous genetic analyses, including the use of additional CaESS1 mutant alleles, distinct marker genes, and the engineering of suitably-matched isogenic control strains. The results confirm that CaEss1 is essential for growth in C. albicans, but show that reduction of CaESS1 gene dosage by half (δ/+) does not interfere with morphogenetic switching. However, further reduction of CaEss1 levels using a conditional allele does reduce morphogenetic switching. We also examine the role of the linker α-helix that distinguishes C. albicans Ess1 from the human Pin1 enzyme, and present results of a genome-wide transcriptome analysis. The latter analysis indicates that CaEss1 has a conserved role in regulation of RNA polymerase II function, and is required for efficient termination of small nucleolar RNAs and repression of cryptic transcription in C. albicans.

Splicing-related features of introns serve to propel evolution

The role of spliceosomal intronic structures played in evolution has only begun to be elucidated. Comparative genomic analyses of fungal snoRNA sequences, which are often contained within introns and/or exons, revealed that about one-third of snoRNA-associated introns in three major snoRNA gene clusters manifested polymorphisms, likely resulting from intron loss and gain events during fungi evolution. Genomic deletions can clearly be observed as one mechanism underlying intron and exon loss, as well as generation of complex introns where several introns lie in juxtaposition without intercalating exons. Strikingly, by tracking conserved snoRNAs in introns, we found that some introns had moved from one position to another by excision from donor sites and insertion into target sties elsewhere in the genome without needing transposon structures. This study revealed the origin of many newly gained introns. Moreover, our analyses suggested that intron-containing sequences were more prone to sustainable structural changes than DNA sequences without introns due to intron's ability to jump within the genome via unknown mechanisms. We propose that splicing-related structural features of introns serve as an additional motor to propel evolution.

Candida albicans: A Model Organism for Studying Fungal Pathogens

Candida albicans is an opportunistic human fungal pathogen that causes candidiasis. As healthcare has been improved worldwide, the number of immunocompromised patients has been increased to a greater extent and they are highly susceptible to various pathogenic microbes and C. albicans has been prominent among the fungal pathogens. The complete genome sequence of this pathogen is now available and has been extremely useful for the identification of repertoire of genes present in this pathogen. The major challenge is now to assign the functions to these genes of which 13% are specific to C. albicans. Due to its close relationship with yeast Saccharomyces cerevisiae, an edge over other fungal pathogens because most of the technologies can be directly transferred to C. albicans from S. cerevisiae and it is amenable to mutation, gene disruption, and transformation. The last two decades have witnessed enormous amount of research activities on this pathogen that leads to the understanding of host-parasite interaction, infections, and disease propagation. Clearly, C. albicans has emerged as a model organism for studying fungal pathogens along with other two fungi Aspergillus fumigatus and Cryptococcus neoformans. Understanding its complete life style of C. albicans will undoubtedly be useful for developing potential antifungal drugs and tackling Candida infections. This will also shed light on the functioning of other fungal pathogens.

Staying on message: ensuring fidelity in pre-mRNA splicing

The faithful expression of genes requires that cellular machinery select substrates with high specificity at each step in gene expression. High specificity is particularly important at the stage of nuclear pre-mRNA splicing, during which the spliceosome selects splice sites and excises intervening introns. With low specificity, the usage of alternative sites would yield insertions, deletions and frame shifts in mRNA. Recently, biochemical, genetic and genome-wide approaches have significantly advanced our understanding of splicing fidelity. In particular, we have learned that DExD/H-box ATPases play a general role in rejecting and discarding suboptimal substrates and that these factors serve as a paradigm for proofreading NTPases in other systems. Recent advances have also defined fundamental questions for future investigations.

Sequence and analysis of the genome of the pathogenic yeast Candida orthopsilosis

Candida orthopsilosis is closely related to the fungal pathogen Candida parapsilosis. However, whereas C. parapsilosis is a major cause of disease in immunosuppressed individuals and in premature neonates, C. orthopsilosis is more rarely associated with infection. We sequenced the C. orthopsilosis genome to facilitate the identification of genes associated with virulence. Here, we report the de novo assembly and annotation of the genome of a Type 2 isolate of C. orthopsilosis. The sequence was obtained by combining data from next generation sequencing (454 Life Sciences and Illumina) with paired-end Sanger reads from a fosmid library. The final assembly contains 12.6 Mb on 8 chromosomes. The genome was annotated using an automated pipeline based on comparative analysis of genomes of Candida species, together with manual identification of introns. We identified 5700 protein-coding genes in C. orthopsilosis, of which 5570 have an ortholog in C. parapsilosis. The time of divergence between C. orthopsilosis and C. parapsilosis is estimated to be twice as great as that between Candida albicans and Candida dubliniensis. There has been an expansion of the Hyr/Iff family of cell wall genes and the JEN family of monocarboxylic transporters in C. parapsilosis relative to C. orthopsilosis. We identified one gene from a Maltose/Galactoside O-acetyltransferase family that originated by horizontal gene transfer from a bacterium to the common ancestor of C. orthopsilosis and C. parapsilosis. We report that TFB3, a component of the general transcription factor TFIIH, undergoes alternative splicing by intron retention in multiple Candida species. We also show that an intein in the vacuolar ATPase gene VMA1 is present in C. orthopsilosis but not C. parapsilosis, and has a patchy distribution in Candida species. Our results suggest that the difference in virulence between C. parapsilosis and C. orthopsilosis may be associated with expansion of gene families.

A recently evolved transcriptional network controls biofilm development in Candida albicans

A biofilm is an organized, resilient group of microbes in which individual cells acquire properties, such as drug resistance, that are distinct from those observed in suspension cultures. Here, we describe and analyze the transcriptional network controlling biofilm formation in the pathogenic yeast Candida albicans, whose biofilms are a major source of medical device-associated infections. We have combined genetic screens, genome-wide approaches, and two in vivo animal models to describe a master circuit controlling biofilm formation, composed of six transcription regulators that form a tightly woven network with ∼1,000 target genes. Evolutionary analysis indicates that the biofilm network has rapidly evolved: genes in the biofilm circuit are significantly weighted toward genes that arose relatively recently with ancient genes being underrepresented. This circuit provides a framework for understanding many aspects of biofilm formation by C. albicans in a mammalian host. It also provides insights into how complex cell behaviors can arise from the evolution of transcription circuits.

Using RNA-seq to determine the transcriptional landscape and the hypoxic response of the pathogenic yeast Candida parapsilosis

BACKGROUND

Candida parapsilosis is one of the most common causes of Candida infection worldwide. However, the genome sequence annotation was made without experimental validation and little is known about the transcriptional landscape. The transcriptional response of C. parapsilosis to hypoxic (low oxygen) conditions, such as those encountered in the host, is also relatively unexplored.

RESULTS

We used next generation sequencing (RNA-seq) to determine the transcriptional profile of C. parapsilosis growing in several conditions including different media, temperatures and oxygen concentrations. We identified 395 novel protein-coding sequences that had not previously been annotated. We removed > 300 unsupported gene models, and corrected approximately 900. We mapped the 5' and 3' UTR for thousands of genes. We also identified 422 introns, including two introns in the 3' UTR of one gene. This is the first report of 3' UTR introns in the Saccharomycotina. Comparing the introns in coding sequences with other species shows that small numbers have been gained and lost throughout evolution. Our analysis also identified a number of novel transcriptional active regions (nTARs). We used both RNA-seq and microarray analysis to determine the transcriptional profile of cells grown in normoxic and hypoxic conditions in rich media, and we showed that there was a high correlation between the approaches. We also generated a knockout of the UPC2 transcriptional regulator, and we found that similar to C. albicans, Upc2 is required for conferring resistance to azole drugs, and for regulation of expression of the ergosterol pathway in hypoxia.

CONCLUSION

We provide the first detailed annotation of the C. parapsilosis genome, based on gene predictions and transcriptional analysis. We identified a number of novel ORFs and other transcribed regions, and detected transcripts from approximately 90% of the annotated protein coding genes. We found that the transcription factor Upc2 role has a conserved role as a major regulator of the hypoxic response in C. parapsilosis and C. albicans.

Deep sequencing of small RNAs in plants: applied bioinformatics

Small RNAs, including microRNA and short-interfering RNAs, play important roles in plants. In recent years, developments in sequencing technology have enabled the large-scale discovery of sRNAs in various cells, tissues and developmental stages and in response to various stresses. This review describes the bioinformatics challenges to analysing these large datasets of short-RNA sequences and some of the solutions to those challenges.

The Candida genome database incorporates multiple Candida species: multispecies search and analysis tools with curated gene and protein information for Candida albicans and Candida glabrata

The Candida Genome Database (CGD, http://www.candidagenome.org/) is an internet-based resource that provides centralized access to genomic sequence data and manually curated functional information about genes and proteins of the fungal pathogen Candida albicans and other Candida species. As the scope of Candida research, and the number of sequenced strains and related species, has grown in recent years, the need for expanded genomic resources has also grown. To answer this need, CGD has expanded beyond storing data solely for C. albicans, now integrating data from multiple species. Herein we describe the incorporation of this multispecies information, which includes curated gene information and the reference sequence for C. glabrata, as well as orthology relationships that interconnect Locus Summary pages, allowing easy navigation between genes of C. albicans and C. glabrata. These orthology relationships are also used to predict GO annotations of their products. We have also added protein information pages that display domains, structural information and physicochemical properties; bibliographic pages highlighting important topic areas in Candida biology; and a laboratory strain lineage page that describes the lineage of commonly used laboratory strains. All of these data are freely available at http://www.candidagenome.org/. We welcome feedback from the research community at candida-curator@lists.stanford.edu.

The mechanism of the second step of pre-mRNA splicing

The molecular mechanisms of the second step of pre-mRNA splicing in yeast and higher eukaryotes are reviewed. The important elements in the pre-mRNA, the participating proteins, and the proposed secondary structures and roles of the snRNAs are described. The sequence of events in the second step is presented, focusing on the actions of the proteins in setting up and facilitating the second reaction. Mechanisms for avoiding errors in splicing are discussed.

The intronome of budding yeasts

Whatever their abundance in genomes, spliceosomal introns are the signature of eukaryotic genes. The sequence of Saccharomyces cerevisiae, achieved fifteen years ago, revealed that this yeast has very few introns, but conserved intron boundaries typical for an intron definition mechanism. With the improvement and the development of new sequencing technologies, yeast genomes have been extensively sequenced during the last decade. We took advantage of this plethora of data to compile and assess the intron content of the protein-coding genes of 13 genomes representative of the evolution of hemiascomycetous yeasts. We first observed that intron paucity is a general rule and that the fastest evolving genomes tend to lose their introns more rapidly (e.g. S. cerevisiae versus Yarrowia lipolytica). Noticeable differences were also confirmed for 5' splice sites and branch point sites (BP) as well as for the relative position of the BP. These changes seemed to be correlated with the lineage specific evolution of splicing factors.

DNA methylation regulates phenotype-dependent transcriptional activity in Candida albicans

DNA methylation is a common epigenetic signaling mechanism associated with silencing of repeated DNA and transcriptional regulation in eukaryotes. Here we report that DNA methylation in the human fungal pathogen Candida albicans is primarily localized within structural genes and modulates transcriptional activity. Major repeat sequences and multigene families are largely free of DNA methylation. Among the genes subject to DNA methylation are those associated with dimorphic transition between yeast and hyphal forms, switching between white and opaque cells, and iron metabolism. Transcriptionally repressed methylated loci showed increased frequency of C-to-T transitions during asexual growth, an evolutionarily stable pattern of repression associated mutation that could bring about genetic alterations under changing environmental or host conditions. Dynamic differential DNA methylation of structural genes may be one factor contributing to morphological plasticity that is cued by nutrition and host interaction.

Why eukaryotic cells use introns to enhance gene expression: splicing reduces transcription-associated mutagenesis by inhibiting topoisomerase I cutting activity

Background

The costs and benefits of spliceosomal introns in eukaryotes have not been established. One recognized effect of intron splicing is its known enhancement of gene expression. However, the mechanism regulating such splicing-mediated expression enhancement has not been defined. Previous studies have shown that intron splicing is a time-consuming process, indicating that splicing may not reduce the time required for transcription and processing of spliced pre-mRNA molecules; rather, it might facilitate the later rounds of transcription. Because the densities of active RNA polymerase II on most genes are less than one molecule per gene, direct interactions between the splicing apparatus and transcriptional complexes (from the later rounds of transcription) are infrequent, and thus unlikely to account for splicing-mediated gene expression enhancement.

Presentation of the hypothesis

The serine/arginine-rich protein SF2/ASF can inhibit the DNA topoisomerase I activity that removes negative supercoiling of DNA generated by transcription. Consequently, splicing could make genes more receptive to RNA polymerase II during the later rounds of transcription, and thus affect the frequency of gene transcription. Compared with the transcriptional enhancement mediated by strong promoters, intron-containing genes experience a lower frequency of cut-and-paste processes. The cleavage and religation activity of DNA strands by DNA topoisomerase I was recently shown to account for transcription-associated mutagenesis. Therefore, intron-mediated enhancement of gene expression could reduce transcription-associated genome instability.

Testing the hypothesis

Experimentally test whether transcription-associated mutagenesis is lower in intron-containing genes than in intronless genes. Use bioinformatic analysis to check whether exons flanking lost introns have higher frequencies of short deletions.

Implications of the hypothesis

The mechanism of intron-mediated enhancement proposed here may also explain the positive correlation observed between intron size and gene expression levels in unicellular organisms, and the greater number of intron containing genes in higher organisms.

Reviewers

This article was reviewed by Dr Arcady Mushegian, Dr Igor B Rogozin (nominated by Dr I King Jordan) and Dr Alexey S Kondrashov. For the full reviews, please go to the Reviewer's Reports section.

The transcriptomes of two heritable cell types illuminate the circuit governing their differentiation

The differentiation of cells into distinct cell types, each of which is heritable for many generations, underlies many biological phenomena. White and opaque cells of the fungal pathogen Candida albicans are two such heritable cell types, each thought to be adapted to unique niches within their human host. To systematically investigate their differences, we performed strand-specific, massively-parallel sequencing of RNA from C. albicans white and opaque cells. With these data we first annotated the C. albicans transcriptome, finding hundreds of novel differentially-expressed transcripts. Using the new annotation, we compared differences in transcript abundance between the two cell types with the genomic regions bound by a master regulator of the white-opaque switch (Wor1). We found that the revised transcriptional landscape considerably alters our understanding of the circuit governing differentiation. In particular, we can now resolve the poor concordance between binding of a master regulator and the differential expression of adjacent genes, a discrepancy observed in several other studies of cell differentiation. More than one third of the Wor1-bound differentially-expressed transcripts were previously unannotated, which explains the formerly puzzling presence of Wor1 at these positions along the genome. Many of these newly identified Wor1-regulated genes are non-coding and transcribed antisense to coding transcripts. We also find that 5' and 3' UTRs of mRNAs in the circuit are unusually long and that 5' UTRs often differ in length between cell-types, suggesting UTRs encode important regulatory information and that use of alternative promoters is widespread. Further analysis revealed that the revised Wor1 circuit bears several striking similarities to the Oct4 circuit that specifies the pluripotency of mammalian embryonic stem cells. Additional characteristics shared with the Oct4 circuit suggest a set of general hallmarks characteristic of heritable differentiation states in eukaryotes.