Murat C. Kalem
John C. Panepinto- Department of Microbiology and Immunology, Witebsky Center for Microbial Pathogenesis and Immunology, Jacobs School of Medicine and Biomedical Sciences, State University of New York (SUNY), University at Buffalo, Buffalo, NY, United States
Long non-coding RNAs (lncRNAs) are highly expressed and can modulate multiple cellular processes including transcription, splicing, translation, and many diverse signaling events. LncRNAs can act as sponges for miRNAs, RNA and DNA binding proteins, functioning as competitive endogenous RNAs. The contribution of lncRNAs to microbial pathogenesis is largely neglected in eukaryotic pathogens despite the abundance of RNA sequencing datasets encompassing conditions of stress, gene deletions and conditions that mimic the host environment. The human fungal pathogen Cryptococcus neoformans encodes 6975 (84%) protein-coding and 1359 (16%) non-protein-coding RNAs, of which 1182 (14.2%) are lncRNAs defined by a threshold of greater than 200 nucleotides in length. Here, we discuss the current state of knowledge in C. neoformans lncRNA biology. Utilizing existing RNA seq datasets, we examine trends in lncRNA expression and discuss potential implications for pathogenesis.
Introduction
Long non-coding RNAs (lncRNAs) were first thought to be transcriptional noise or junk RNAs. Although they are named for what they do not do – encode for proteins – we now know that lncRNAs have an impressive range of biological functions (Ponting et al., 2009; Rinn and Chang, 2012; Scaria and Pasha, 2013; Akhade et al., 2017; Zhang et al., 2019). The central dogma of molecular biology states that genetic information is transferred from DNA to RNA to protein, creating an incomplete narrative that an RNA molecule only has value if it encodes for a protein. This dogma, accepted as the backbone of molecular biology, constrains our vision of the vast functional diversity that exists among RNA molecules. Over the last two decades, this dogma has been challenged by non-coding RNAs and major advances in the field of RNA biology.
Non-coding RNAs include a diverse cohort of RNA molecules that can be categorized under two broad groups: housekeeping ncRNAs and regulatory ncRNAs. Housekeeping ncRNAs include ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), small nuclear RNAs (snRNAs), small nucleolar RNAs (snoRNAs), and tRNA fragments (tRFs). Regulatory ncRNAs include micro RNAs (miRNAs), small interfering RNAs (siRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) (Rinn and Chang, 2012; Akhade et al., 2017; Zhang et al., 2019). LncRNAs are classified as non-coding RNA molecules that are larger than 200 bases. Classification of the lncRNAs solely based on length is arbitrary. To broaden this definition, it is proposed that non-coding RNAs that function as primary or spliced RNAs outside of the known small ncRNA classes are also considered lncRNAs (Amaral et al., 2011). Consequently, there are some lncRNAs that do not meet the arbitrary length threshold such as BC1 and snaR (Parrott and Mathews, 2007; Kim et al., 2016).
Steps in RNA metabolism from mRNA transcription to maturation might be differentially modulated for lncRNAs. As products of RNA polymerase II transcription, lncRNAs may or may not be processed through canonical spliceosome-dependent splicing, 3’ end processing and polyadenylation. For example, human lncRNAs MALAT1 and NEAT1_2 have 3’ triple helix structures instead of poly-A tails (Brown et al., 2012; Wilusz et al., 2012). This suggests that alternatively processed lncRNAs may have unique post-transcriptional regulatory mechanisms, possibly through RNA-binding proteins that recognize and bind to these unique structural elements. Interestingly, MALAT1 3’ triple helix structures are shown to engage with ribosomes suggesting a potential role in translation regulation (Bhandari et al., 2014). This type of poly-A-independent mechanism is not uncommon, and in fact, histone mRNAs contain 3’ stem-loops instead of poly-A tails and are regulated by the binding of stem-loop binding proteins (Wang et al., 1996).
LncRNAs localize to both the cytoplasm and the nucleus, therefore they are able to regulate molecular processes across multiple cellular locations (Cabili et al., 2015; Dunagin et al., 2015; Miao et al., 2019). They can regulate gene expression at the level of transcription, can impact post-transcriptional processes, and modulate translation. LncRNAs can act as a guide for chromatin-modifying enzymes, and as decoys for RNA or DNA binding proteins or other signaling molecules (Wang and Chang, 2011; Chacko and Lin, 2013; Akhade et al., 2017). They can act as endogenous competitors, or as sponges for miRNAs, and they themselves can be precursors that are processed to form small RNAs (Wang and Chang, 2011; Sweta et al., 2019). Many lncRNAs lack a canonical open reading frame, yet there is recent evidence that some lncRNAs produce micropeptides with functional roles (Hartford and Lal, 2020). For example, micropeptide Myoregulin modulates muscle development in mice and another micropeptide Sda controls sporulation in Bacillus subtilis (Burkholder et al., 2001; Anderson et al., 2016).
A single lncRNA can exert robust and widespread molecular control in cellular processes; one great example is the lncRNA Xist which can silence an entire chromosome by recruiting necessary effector proteins (Sunwoo et al., 2015; Jachowicz et al., 2021). Studies in mammalian cells and model systems implicate lncRNAs as important players in cellular and molecular processes and demonstrate their impact on disease pathologies including cancers, neurodegenerative and autoimmune diseases (Gao et al., 2018; Wang et al., 2018; Piccolo et al., 2019; Zhou et al., 2021). Despite their importance, the contribution of lncRNAs to eukaryotic gene regulation in fungi, especially in fungal pathogens, is a neglected area of study.
Here, we aim to briefly highlight some findings on the importance of lncRNAs in pathogenic fungi and then focus on lncRNAs in the human fungal pathogen Cryptococcus neoformans. In this perspective, we curate the lncRNAs of C. neoformans and utilize the published RNA sequencing data sets to investigate lncRNA expression under conditions that mimic host microenvironments. Given their large numbers and dynamic expression across stress conditions, we are hopeful that further investigation into lncRNA function will identify lncRNA molecules involved in fungal-specific virulence factor production, stress adaptation, and pathogenesis.
Non-Coding RNAs in Fungal Pathogens
Although studies of lncRNAs in fungi are not as abundant as in other kingdoms of life, there is evidence that lncRNAs have regulatory roles in many fungal species (Chacko and Lin, 2013; Till et al., 2018; Dhingra, 2020; Li et al., 2021). A majority of the knowledge on lncRNAs in fungi comes from the model yeasts S. cerevisiae and S. pombe (Yamashita et al., 2016). In Saccharomyces cerevisiae, there are both cis-acting lncRNAs, which regulate neighboring genes at a single genetic locus, and trans-acting lncRNAs, which can exert regulatory function outside of the genetic locus that encodes it. For example, CDC28-antisense lncRNA, pHO-lncRNA, PWR1, ICR1, RME2, and RME3 are identified to be cis-acting since ectopic expression of these lncRNAs presents a loss-of-function due to the absence of lncRNA-mediated cis-regulation (Hongay et al., 2006; Bumgarner et al., 2009; Nadal-Ribelles et al., 2014; Yu et al., 2016; Till et al., 2018). Cis-acting RME2 and RME3 have unique functions in regulating meiosis (Hongay et al., 2006). Although trans-acting lncRNAs are less common in S. cerevisiae than cis-acting lncRNAs, SUT457 is an example of an S. cerevisiae trans-acting lncRNA that interacts with 12 target genes that are important for telomere organization (Kyriakou et al., 2016).
Recent studies also identified lncRNAs in the plant pathogens Fusarium graminearum and Magnaporthe oryzae, the insect pathogen Metarhizium robertsii, the model organism Neurospora crassa, and human fungal pathogens Candida auris, Aspergillus fumigatus, and Cryptococcus neoformans. In the plant pathogenic fungus Fusarium graminearum, lncRNAs exhibit dynamic expression patterns during sexual fruiting body formation, suggesting a regulatory role for lncRNAs in fungal development (Kim et al., 2018). In the insect pathogen Metarhizium robertsii, lncRNAs are regulated during heat shock, 1085 out of 1655 lncRNAs were shown to be differentially expressed, both up- and down-regulated, following heat stress (Wang et al., 2019). One specific lncRNA, mrLncRNA1, plays a role in conidial thermotolerance, and in the absence of mrLncRNA1, M. robertsii exhibits a defect in germination (Wang et al., 2019). Bioinformatic analysis of lncRNAs in Candida pathogens (C. albicans, C. tropicalis, C. parapsilosis, C. auris, and C. glabrata) demonstrated the co-regulation of lncRNA expression with certain protein-coding genes while identifying many differentially regulated lncRNAs during the infection of human epithelial cells (Hovhannisyan and Gabaldón, 2021). Another important fungal lncRNA was identified during a piggyBac transposon mutagenesis screen in C. auris that resulted in a mutant that exhibited constitutive filamentous growth (Gao et al., 2021). This mutant had a transposon insertion in a lncRNA gene named DNA damage inducible lncRNA – DINOR (Gao et al., 2021). Deletion of DINOR resulted in DNA damage exhibited as fragmented DNA and caused hyperphosphorylation of kinase Rad53 while leading to attenuated virulence in mice. DINOR expression was upregulated in the presence of DNA damaging alkylating agents, antifungals such as amphotericin B and fluconazole, and peroxide-induced oxidative stress (Gao et al., 2021). The discovery of DINOR, along with other studies highlighting the importance of lncRNAs in fungal pathogens, provides a compelling proof-of-concept that lncRNAs can have robust roles in fungal virulence.
In Cryptococcus neoformans, there is only one lncRNA with an identified function. LncRNA RZE1 was identified through a forward genetic screen in C. neoformans to identify upstream regulators of Znf2 (Chacko et al., 2015). ZNF2 deletion locks C. neoformans in the yeast form causing increased virulence. Conversely, ZNF2 overexpression results in hyperfilamentous phenotype and decreased virulence. LncRNA RZE1 exhibits mainly nuclear localization and exerts its function on the ZNF2 mRNA. RZE1 deletion results in decreased ZNF2 transcript levels and altered subcellular localization of the ZNF2 mRNA (Chacko et al., 2015). It is likely that more than one lncRNA in C. neoformans exerts biological function, and the remaining unstudied lncRNAs may represent unique fungal biology, and importantly, novel virulence regulatory mechanisms. In the rest of this perspective, we present the characterization of C. neoformans lncRNAs and explore their expression profiles in published RNA sequencing studies. We believe that studying lncRNA biology will provide novel insights into Cryptococcal pathogenesis and regulation of virulence.
Characterization of Long Non-Coding RNAs in Cryptococcus neoformans
To begin with, we analyzed the Cryptococcus neoformans genome for annotated lncRNAs using the annotations on FungiDB (Basenko et al., 2018). C. neoformans encodes for 6975 (84%) protein-coding and 1359 (16%) non-protein-coding RNAs, of which 1182 (14.2%) are lncRNAs with a 200-nucleotide length threshold. We scrutinized the various features of coding and non-coding transcripts. First, we looked at the distribution of the lncRNA encoding genes across the genome and found that they are equally distributed across 14 chromosomes with no apparent enrichment (Figure 1A). The average transcript size for coding and long non-coding RNAs revealed that lncRNAs are on average shorter than mRNAs (Figure 1B). Despite not being translated, lncRNAs can still undergo splicing and other RNA processing steps. We analyzed the number of exons in coding and non-coding transcripts. This analysis revealed that more than ~60% percent of non-coding RNAs either contain a single exon (no intron) or very few exons in comparison to mRNAs (Figure 1C). RNA sequencing evidence suggests that lncRNA introns are as efficiently spliced as most mRNAs. LncRNAs do not have an accumulation of intronic sequences. LncRNAs can be classified based on their relative location to protein-coding genes (Figure 1D). These classes are 1) sense – lncRNA overlaps another gene on the same strand, 2) antisense – lncRNA overlaps another gene on the opposite strand, 3) bidirectional – transcription of lncRNA and another gene on the opposite strand are initiated in close proximity, 4) intronic – full lncRNA is derived from an intron, 5) intergenic – lncRNA exists between multiple genes. Investigations of these different class types using GffCompare and annotations on FungiDB resulted in that majority of the lncRNAs in C. neoformans are antisense lncRNAs that have exonic overlap with the genes on the opposite strand (Pertea and Pertea, 2020). This was followed by the intergenic lncRNAs being the second most common class type in C. neoformans lncRNAs (Figure 1E and Supplementary File 1). This suggests a profound possibility that lncRNA expression in C. neoformans can possess regulation in cis and impact the expression of neighboring and overlapping genes.
Figure 1 Characterization of annotated long non-coding RNAs in Cryptococcus neoformans. C. neoformans (strain H99) non-coding RNA sequences were obtained from FungiDB. We investigated the genomic locations of lncRNAs across chromosomes (A), transcript length (B), and the number of exons (C). LncRNAs are classified based on their juxtaposition relative to protein-coding genes. We determined the lncRNA class types using the GffCompare program (D, E). We investigated the occupancy of lncRNAs by ribosomes by plotting ribosome-protected fragments (RPF) vs. RNA abundance using a 10 reads minimum threshold. Regression analysis, shown as black lines, and corresponding equations with R2 values are also presented (F). Translation efficiency was calculated by using the RPF/RNA formula and plotted for both coding and non-coding RNAs (G). Ribosome profiling data is obtained from GSE133125. Bed files supplied by the authors were utilized to obtain the coding and non-coding transcript counts.
Following this initial characterization of C. neoformans lncRNAs, we utilized the previously published ribosome profiling and RNA sequencing data sets to ask if any of the identified lncRNAs are ribosome-associated (Wallace et al., 2020). Of note, our analysis also includes some shorter lncRNAs as lncRNA length classification is arbitrary. Analysis of the ribosome protected fragments (RPFs) and RNA abundance revealed that while protein-coding transcripts have a more linear correlation between RPFs and RNA abundance (: 0.7378), long non-coding transcripts do not have as robust of a relationship between RPFs and RNA abundance (: 0.4491) (Figure 1F and Supplementary File 1). Following this, we calculated the translation efficiency (TE) by dividing the RPF read count by RNA read count. This analysis revealed that although long non-coding RNAs have lower average TE values, there are some ncRNAs with TE values as high as that of protein-coding transcripts (Figure 1G and Supplementary File 1). One can speculate multiple scenarios that could explain this. First, these lncRNAs might be misannotated and they might encode unstable peptides that are yet to be identified. Second, ribosome profiling experiments include ribosome purification and nuclease digest steps to selectively isolate the RNA fragments protected by ribosomes. Despite all efforts to only isolate ribosome-bound RNAs, there might be some RNAs that are protected from the nuclease digest by a large protein complex such as RNA-binding protein complexes. And third, these ncRNAs with high TE values might actually be bound by ribosomes and act as endogenous competitive RNAs to occupy and eliminate ribosomes from translating other transcripts. These hypotheses provide exciting possibilities for the identification of fungal-specific gene regulatory mechanisms as most lncRNAs in Cryptococcus do not seem to be conserved across kingdoms based on BLAST analysis.
Harnessing the Power of Next-Generation Sequencing to Investigate lncRNA Expression in Cryptococcus neoformans
RNA sequencing became one of the most popular and commonly utilized methods over the past two decades. This widely-appreciated technological advance not only allowed scientists to appreciate the transcriptome-wide gene expression patterns but also created this enormous pool of data that can be mined to address many neglected hypotheses - one of which is that fungal lncRNAs have dynamic expression patterns and are important in modulating virulence and virulence-associated traits in the fungal pathogen C. neoformans. We investigated the expression of lncRNAs in C. neoformans in response to infection-relevant stressors using published data sets. Some of the major stressors and unique microenvironments that environmental fungi encounter during infection include elevated temperature, growth within macrophages in the presence of oxidative stress insult, and growth within the central nervous system such as cerebrospinal fluid (CSF). For that reason, we utilized RNA sequencing data obtained from C. neoformans cells phagocytosed by activated J774A murine macrophages following 2 hours co-infection, C. neoformans isolated from corticosteroid treated New Zealand white rabbit CSF following 24 hours post intracisternal inoculation, and C. neoformans cells grown at 30°C and shifted to 37°C for one hour (Bloom et al., 2019; Yu et al., 2020).
Analysis of the expression of both coding and long non-coding RNAs revealed massive remodeling of both kinds of transcripts. The macrophage condition had 572 upregulated and 108 downregulated lncRNAs compared to YPD (Figure 2A and Supplementary File 1). The CSF condition had 423 upregulated and 128 downregulated lncRNAs compared to YPD (Figure 2B and Supplementary File 1). Under elevated temperature challenge, we identified 47 upregulated and 38 downregulated lncRNAs with a 0.05 adjusted p-value and 2-fold change threshold (Figure 2C and Supplementary File 1). These results suggest that lncRNA expression, along with mRNA expression, is dynamic under infection-relevant stress conditions and likely crucial in fine-tuning gene regulatory or signaling pathways necessary to adapt to constantly and rapidly changing environments.
Figure 2 Long non-coding RNA expression under conditions that mimic the host environment and in gene deletion mutants. We utilized the public datasets to investigate the expression of lncRNAs under conditions that are relevant to infection. These conditions include in vitro J774A macrophage infection (A), CSF (B), and 37˚C compared to 30˚C in YPD (C). The number of upregulated and downregulated genes were noted on each graph for both coding (orange) and non-coding (purple) genes. We then show the impact of RRP6 (D) and AGO1 (E) gene deletions on lncRNA expression as well as changes in the expression of small RNAs that map to full-length lncRNAs. We show the expression of antisense lncRNAs and antisense protein-coding transcripts in the rrp6Δ (F). RNA sequencing data were obtained from GSE136879, GSE121183 and GSE128009. GSE136879 data files were trimmed and mapped to the H99 genome using STAR alignment and gene expression analysis was performed using RSEM and DeSeq2. GSE121183 and GSE128009 data were analyzed for ncRNA expression using the supplementary tables provided by the authors of these studies. (G) Graphical summary/model. LncRNAs can have unique functions in distinct cellular and extracellular compartments that impact pathogenesis and immune modulation. The graphical summary is created using biorender.com.
In addition to RNA sequencing data sets from host-relevant conditions, there is an abundance of transcriptome-wide analyses from many gene deletion mutants of C. neoformans. Investigations of lncRNA expression in these mutants can yield the determination of the key factors that control lncRNA expression. Expression of lncRNAs in gene deletion mutants along with the phylogenetic conservation of deleted gene products could hint us at the evolutionarily conserved or fungi-specific mechanisms of lncRNA expression modulation.
We utilized published rrp6Δ and ago1Δ RNA sequencing data sets and determined the expression of lncRNAs (Burke et al., 2019). Rrp6 is a 3’ to 5’ exoribonuclease and a critical component of the nuclear RNA exosome which is the complex responsible for catalyzing degradation, processing, surveillance, and export of coding and non-coding RNAs (Davis and Ares, 2006; Zhang et al., 2014; Fox and Mosley, 2016). The RNA exosome complex regulates the maturation of some housekeeping ncRNAs such as rRNAs, snoRNAs, and snRNAs (Fox and Mosley, 2016). It has been shown that mutations in Rrp6 yield aberrant polyadenylation of snoRNA transcripts in yeast and these misregulated transcripts are then rapidly degraded in an Rrp6-dependent fashion (Grzechnik and Kufel, 2008). Substrates of Rrp6 include many non-coding RNAs (Chacko et al., 2015). In humans, Rrp6 deletion results in increased transcription upstream of promoters, called promoter upstream transcripts or PROMPTs. PROMPTs are stabilized in exosome depleted cells and some even have functions (Preker et al., 2008). Rrp6-mediated RNA decay has shown to be important for eliminating many non-coding RNAs that arise from pervasive transcription and transcriptional noise (Davis and Ares, 2006). Another study highlights that an Rrp6-like protein in Arabidopsis maintains DNA methylation through retention of lncRNAs, which provides an example in which Rrp6 can positively regulate lncRNA expression (Zhang et al., 2014). Analysis of lncRNA expression in the C. neoformans rrp6Δ mutant revealed that 538 lncRNAs were upregulated and 230 were downregulated in the mutant compared to wildtype (Figure 2D). This suggests that while the majority of lncRNAs are upregulated in the absence of Rrp6, the nuclear exosome might modulate lncRNA expression both negatively and positively.
Another dataset we investigated was from the Argonaute 1 (Ago1) deletion mutant (Burke et al., 2019). Ago1 binds to short RNAs and is a critical component of the RNAi pathway and small RNA mediated gene silencing (Williams and Rubin, 2002). LncRNAs can act as precursors for miRNAs and small RNAs (Arikit et al., 2013). Since the majority of the lncRNAs are antisense to other coding genes, we reasoned that small RNAs generated from lncRNAs could regulate the expression of coding genes expressed on the opposite strand through small RNA mediated gene silencing. We analyzed the expression of ncRNAs and small RNAs that map to full-length ncRNAs in the ago1Δ mutant. Interestingly, this mutant did not impact the expression of lncRNAs, however, the mutant had 14 upregulated and 78 downregulated small RNAs that map to lncRNAs (Figure 2E). This suggests that Ago1 might bind to small RNAs that are generated from lncRNAs, and in the absence of Ago1 these small RNAs are likely degraded or dysregulated. Investigation of small RNAs originating from lncRNAs in the rrp6Δ showed that 59 were upregulated and 16 were downregulated (Figure 2E).
Since many C. neoformans lncRNAs are antisense to protein-coding genes, we investigated the expression of lncRNAs and their corresponding antisense protein-coding genes in the rrp6Δ (Figure 2F). Our analysis revealed that upregulated antisense lncRNAs can both positively and negatively influence the expression of the antisense protein-coding gene. This suggests that Rrp6-mediated lncRNA expression regulation may alter the antisense gene expression in a context-dependent manner. These observations utilizing published RNA sequencing datasets will allow us to generate testable hypotheses to explore lncRNA biology in C. neoformans.
Discussion
Because evolution is an opportunistic process, we can speculate that a majority of lncRNAs likely have functions otherwise they would not exist. Even though lncRNAs are not highly conserved in terms of linear sequence homology, conserved biological functions, and even structures, are likely under evolutionary selection (Johnsson et al., 2014; Ulitsky, 2016; Ramírez-Colmenero and Fernandez-Valverde, 2020). LncRNAs DINOR and RZE1 are the best examples highlighting the functional roles of lncRNAs in fungi yet they only represent the tip of the iceberg (Chacko et al., 2015; Gao et al., 2021). Both RZE1 and DINOR are only conserved within their respective species complexes based on sequence analyses, yet this does not exclude the possibility of functional conservation without retaining sequence homology.
Advances in sequencing technologies revealed the diversity of RNA molecules and allowed us to investigate the expression of coding and non-coding transcripts under various conditions. However, post-transcriptional processing of lncRNAs is unclear, and RNA-sequencing experiments including a poly-A purification method may exclude a subset of lncRNAs that are either non-polyadenylated or are conditionally polyadenylated. With this in mind, all of the data sets utilized in our analyses are biased for polyadenylated RNAs. Optimization and widespread utilization of efficient RNase H based rRNA depletion methods utilizing antisense oligonucleotides will improve our understanding of the diverse elements within fungal transcriptomes including non-coding RNAs. Recent developments in the field revealed that RNase H based rRNA removal is efficient in C. neoformans and analysis of lncRNAs using LncPipe following an RNase H based rRNA depletion method revealed 11 novel lncRNAs (Telzrow et al., 2021). Our brief analysis of lncRNA expression utilizing published sequencing data sets is not able to determine the difference between changes in expression versus polyadenylation or deadenylation. Future work investigating the non-coding RNA processing networks under host-relevant stress conditions could reveal unique fungal biology, including pathways for lncRNA maturation and processing.
C. neoformans secretes extracellular vesicles containing components of the polysaccharide capsule - one of its virulence factors - that are important for macrophage modulation (Rodrigues et al., 2008; Oliveira et al., 2010). Additionally, C. neoformans EVs are packaged with protein and RNA cargo (Da Silva et al., 2015; Rizzo et al., 2021). The RNA cargo includes a wide variety of RNA molecules including extracellular RNAs (exRNAs), miRNAs, siRNAs and lncRNAs (Liu et al., 2020). We can speculate that the secreted RNA cargo might have an immunomodulatory role during infection or amoeba predation (Figure 2G). For example, a deployed small silencing RNA from entomopathogenic fungus Beauveria bassiana binds to host cell Ago1 and hijacks the host RNA interference machinery (Cui et al., 2019). Further studies are needed to explore the immunomodulatory mechanisms of secreted non-coding RNAs of C. neoformans.
This brief characterization and expression analysis highlight the dynamic expression patterns of lncRNAs in C. neoformans. We hypothesize that changes in lncRNA expression may result in altered nuclear and cytoplasmic post-transcriptional events that impact virulence factors (Figure 2G). We also suggest that lncRNAs play functional roles in regulating C. neoformans pathogenesis (Figure 2G). Future studies focusing on individual lncRNAs will delineate functions for each lncRNA in modulating cellular processes important for virulence and pathogenesis. Even though the majority of the C. neoformans lncRNAs are antisense to a protein-coding gene which makes it impossible to generate specific deletions using homologous recombination, the latest advances in gene manipulations utilizing CRISPR Cas9 and Cas13 will make it possible to more efficiently knockout and knockdown lncRNAs, respectively. These tools could be utilized to create lncRNA knockout and knockdown libraries in order to screen for lncRNA involvement in virulence and pathogenesis as well as to identify lncRNA mutant strains that are sensitive or resistant to current antifungals. We are confident that the studies of lncRNAs in fungi will continue to be fruitful and lead to discoveries of unique fungal biology.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Author Contributions
MCK conceptualized the piece, performed the analyses and drafted the manuscript. JCP edited the manuscript. Both authors contributed to the article and approved the submitted version.
Funding
National Institute of Health (NIH) R01AI131977 to JCP. American Heart Association Fellowship Award ID 827110 to MCK.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2022.858317/full#supplementary-material
References
Akhade, V. S., Pal, D., Kanduri, C. (2017). Long Non Coding RNA: Genome Organization and Mechanism of Action. Long Non Coding RNA Biol. 1008, 47–75. doi: 10.1007/978-981-10-5203-3_2
Amaral, P. P., Clark, M. B., Gascoigne, D. K., Dinger, M. E., Mattick, J. S. (2011). Lncrnadb: A Reference Database for Long Noncoding RNAs. Nucleic Acids Res. 39, 146–151. doi: 10.1093/nar/gkq1138
Anderson, D. M., Makarewich, C. A., Anderson, K. M., Shelton, J. M., Bezprozvannaya, S., Bassel-duby, R., et al. (2016). Widespread Control of Calcium Signaling by a Family of SERCA-Inhibiting Micropeptides. Sci. Signal. 9 (457), ra119. doi: 10.1126/scisignal.aaj1460
Arikit, S., Zhai, J., Meyers, B. C. (2013). Biogenesis and Function of Rice Small RNAs From Non-Coding RNA Precursors. Curr. Opin. Plant Biol. 16, 170–179. doi: 10.1016/j.pbi.2013.01.006
Basenko, E., Pulman, J., Shanmugasundram, A., Harb, O., Crouch, K., Starns, D., et al. (2018). FungiDB: An Integrated Bioinformatic Resource for Fungi and Oomycetes. J. Fungi 4, 39. doi: 10.3390/jof4010039
Bhandari, D., Raisch, T., Weichenrieder, O., Jonas, S., Izaurralde, E. (2014). Structural Basis for the Nanos-Mediated Recruitment of the CCR4 – NOT Complex and Translational Repression. Genes Dev. 28 (8), 888–901. doi: 10.1101/gad.237289.113
Bloom, A. L. M., Jin, R. M., Leipheimer, J., Bard, J. E., Yergeau, D., Wohlfert, E. A., et al. (2019). Thermotolerance in the Pathogen Cryptococcus Neoformans is Linked to Antigen Masking via mRNA Decay-Dependent Reprogramming. Nat. Commun. 10, 1–13. doi: 10.1038/s41467-019-12907-x
Brown, J. A., Valenstein, M. L., Yario, T. A., Tycowski, K. T., Steitz, J. A. (2012). Formation of Triple-Helical Structures by the 3 ′ -End Sequences of MALAT1 and MEN β Noncoding RNAs 2012. Proc. Natl. Acad. Sci. 109 (47), 19202–19207. doi: 10.1073/pnas.1217338109
Bumgarner, S. L., Dowell, R. D., Grisafi, P., Gifford, D. K., Fink, G. R. (2009). Toggle Involving Cis-Interfering Noncoding RNAs Controls Variegated Gene Expression in Yeast. Proc. Natl. Acad. Sci. 106 (43), 1–6. doi: 10.1073/pnas.0909641106
Burke, J. E., Longhurst, A. D., Natarajan, P., Rao, B., Liu, J., Sales-Lee, J., et al. (2019). A Non-Dicer RNase III and Four Other Novel Factors Required for RNAi-Mediated Transposon Suppression in the Human Pathogenic Yeast Cryptococcus Neoformans. G3: Genes Genomes Genet. 9 (7), 2235–2244. doi: 10.1534/g3.119.400330
Burkholder, W. F., Kurtser, I., Grossman, A. D. (2001). Replication Initiation Proteins Regulate a Developmental Checkpoint in Bacillus Subtilis. Cell. 104 (2), 269–279. doi: 10.1016/S0092-8674(01)00211-2
Cabili, M. N., Dunagin, M. C., Mcclanahan, P. D., Biaesch, A., Padovan-merhar, O., Regev, A., et al. (2015). Localization and Abundance Analysis of Human lncRNAs at Single-Cell and Single-Molecule Resolution. Genome Biol. 16 (1), 1–16. doi: 10.1186/s13059-015-0586-4
Chacko, N., Lin, X. (2013). Non-Coding RNAs in the Development and Pathogenesis of Eukaryotic Microbes. Appl. Microbiol. Biotechnol. 97, 7989–7997. doi: 10.1007/s00253-013-5160-y
Chacko, N., Zhao, Y., Yang, E., Wang, L., Cai, J. J., Lin, X. (2015). The lncRNA RZE1 Controls Cryptococcal Morphological Transition. PloS Genet. 11, 1–27. doi: 10.1371/journal.pgen.1005692
Cui, C., Wang, Y., Liu, J., Zhao, J., Sun, P., Wang, S. (2019). A Fungal Pathogen Deploys a Small Silencing RNA That Attenuates Mosquito Immunity and Facilitates Infection. Nat. Commun. 10, 1–10. doi: 10.1038/s41467-019-12323-1
Da Silva, R. P., Puccia, R., Rodrigues, M. L., Oliveira, D. L., Joffe, L. S., César, G. V., et al. (2015). Extracellular Vesicle-Mediated Export of Fungal RNA. Sci. Rep. 5 (7763), 1–12. doi: 10.1038/srep07763
Davis, C. A., Ares, M. (2006). Accumulation of Unstable Promoter-Associated Transcripts Upon Loss of the Nuclear Exosome Subunit Rrp6p in Saccharomyces Cerevisiae. Proc. Natl. Acad. Sci. U.S.A. 103, 3262–3267. doi: 10.1073/pnas.0507783103
Dhingra, S. (2020). Role of Non-Coding RNAs in Fungal Pathogenesis and Antifungal Drug Responses. Curr. Clin. Microbiol. Rep. 7, 133–141. doi: 10.1007/s40588-020-00151-7
Dunagin, M., Cabili, M. N., Rinn, J., Raj, A. (2015). Nuclear Bodies and Noncoding RNAs. Methods Mol. Biol. 1262, 3–19. doi: 10.1007/978-1-4939-2253-6_1
Fox, M. J., Mosley, A. L. (2016). Rrp6: Integrated Roles in Nuclear RNA Metabolism and Transcription Termination. Wiley Interdiscip Rev RNA. 7, 91–104. doi: 10.1002/wrna.1317
Gao, J., Chow, E. W. L., Wang, H., Xu, X., Cai, C., Song, Y., et al. (2021). LncRNA DINOR Is a Virulence Factor and Global Regulator of Stress Responses in Candida Auris. Nat. Microbiol. 6 (7), 842–851. doi: 10.1038/s41564-021-00915-x
Gao, Y., Li, S., Zhang, Z., Yu, X., Zheng, J. (2018). The Role of Long Non-Coding RNAs in the Pathogenesis of RA, SLE, and SS. Front. Med. 5 (193), 1–14. doi: 10.3389/fmed.2018.00193
Grzechnik, P., Kufel, J. (2008). Polyadenylation Linked to Transcription Termination Directs the Processing of snoRNA Precursors in Yeast. Mol. Cell 32, 247–258. doi: 10.1016/j.molcel.2008.10.003
Hartford, C. C. R., Lal, A. (2020). When Long Noncoding Becomes Protein Coding. Mol. Cell Biol. 40 (6), e00528–19. doi: 10.1128/MCB.00528-19
Hongay, C. F., Grisafi, P. L., Galitski, T., Fink, G. R. (2006). Antisense Transcription Controls Cell Fate in Saccharomyces Cerevisiae. Cell 127 (4), 735–745. doi: 10.1016/j.cell.2006.09.038
Hovhannisyan, H., Gabaldón, T. (2021). The Long non-Coding RNA Landscape of Candida Yeast Pathogens. Nat. Commun. 12, 1–13. doi: 10.1038/s41467-021-27635-4
Jachowicz, J. W., Strehle, M., Banerjee, A. K., Thai, J., Blanco, M. R., Guttman, M. (2021). Xist Spatially Amplifies SHARP Recruitment to Balance Chromosome-Wide Silencing and Specificity to the X Chromosome. bioRxiv. doi: 10.1101/2021.10.27.466149
Johnsson, P., Lipovich, L., Grandér, D., Morris, K. V. (2014). Biochimica Et Biophysica Acta Evolutionary Conservation of Long non-Coding RNAs; Sequence, Structure, Function. BBA - Gen. Subj. 1840 (3), 1063–1071. doi: 10.1016/j.bbagen.2013.10.035
Kim, M. Y., Hwang, D. W., Li, F., Choi, Y., Byun, J. W., Kim, D., et al. (2016). Detection of Intra-Brain Cytoplasmic 1 (BC1) Long Noncoding RNA Using Graphene Oxide- Fluorescence Beacon Detector. Scientific Rep. 6 (1), 1–7. doi: 10.1038/srep22552
Kim, W., Miguel-Rojas, C., Wang, J., Townsend, J. P., Trail, F. (2018). Developmental Dynamics of Long Noncoding RNA Expression During Sexual Fruiting Body Formation in Fusarium Graminearum. MBio. 9, 1–17. doi: 10.1128/mBio.01292-18
Kyriakou, D., Stavrou, E., Demosthenous, P., Angelidou, G., Luis, B. S., Boone, C., et al. (2016). Functional Characterisation of Long Intergenic non-Coding RNAs Through Genetic Interaction Profiling in Saccharomyces Cerevisiae. BMC Biol. 14 (1), 1–16. doi: 10.1186/s12915-016-0325-7
Li, J., Liu, X., Yin, Z., Hu, Z., Zhang, K. Q. (2021). An Overview on Identification and Regulatory Mechanisms of Long Non-Coding RNAs in Fungi. Front. Microbiol. 12, 1–15. doi: 10.3389/fmicb.2021.638617
Liu, M., Zhang, Z., Ding, C., Wang, T., Kelly, B., Wang, P. (2020). Transcriptomic Analysis of Extracellular RNA Governed by the Endocytic Adaptor Protein Cin1 of Cryptococcus Deneoformans. Front. Cell Infect. Microbiol. 10, 1–13. doi: 10.3389/fcimb.2020.00256
Miao, H., Wang, L., Zhan, H., Dai, J., Chang, Y., Wu, F., et al. (2019). A Long Noncoding RNA Distributed in Both Nucleus and Cytoplasm Operates in the PYCARD-Regulated Apoptosis by Coordinating the Epigenetic and Translational Regulation. PloS Genet. 15, 1–24. doi: 10.1371/journal.pgen.1008144
Nadal-Ribelles, M., Solé, C., Xu, Z., Steinmetz, L. M., de Nadal, E., Posas, F. (2014). Control of Cdc28 CDK1 by a Stress-Induced lncRNA. Mol. Cell. 53 (4), 549–561. doi: 10.1016/j.molcel.2014.01.006
Oliveira, D. L., Freire-de-Lima, C. G., Nosanchuk, J. D., Casadevall, A., Rodrigues, M. L., Nimrichter, L. (2010). Extracellular Vesicles From Cryptococcus Neoformans Modulate Macrophage Functions. Infect. Immun. 78, 1601–1609. doi: 10.1128/IAI.01171-09
Parrott, A. M., Mathews, M. B. (2007). Novel Rapidly Evolving Hominid RNAs Bind Nuclear Factor 90 and Display Tissue-Restricted Distribution. Nucleic Acids Res. 35 (18), 6249–6258. doi: 10.1093/nar/gkm668
Pertea, G., Pertea, M. (2020). GFF Utilities: GffRead and GffCompare. F1000Research 9, 1–19. doi: 10.12688/f1000research.23297.2
Piccolo, L. L., Mochizuki, H., Nagai, Y. (2019). The lncRNA Hsrω Regulates Arginine Dimethylation of Human FUS to Cause its Proteasomal Degradation in Drosophila. J. Cell Sci. 132 (20), 1–14. doi: 10.1242/jcs.236836
Ponting, C. P., Oliver, P. L., Reik, W. (2009). Evolution and Functions of Long Noncoding RNAs. Cell 136, 629–641. doi: 10.1016/j.cell.2009.02.006
Preker, P., Nielsen, J., Kammler, S., Lykke-Andersen, S., Christensen, M. S., Mapendano, C. K., et al. (2008). RNA Exosome Depletion Reveals Transcription Upstream of Active Human Promoters. Science 322, 1851–1854. doi: 10.1126/science.1164096
Ramírez-Colmenero, A., Fernandez-Valverde, S. L. (2020). Evolution of Genome-Organizing Long Non-Coding RNAs in Metazoans. Front. Genet. 11, 1512. doi: 10.3389/fgene.2020.589697
Rinn, J. L., Chang, H. Y. (2012). Genome Regulation by Long Noncoding RNAs. Annu. Rev. Biochem. 81, 145–166. doi: 10.1146/annurev-biochem-051410-092902
Rizzo, J., Taheraly, A., Janbon, G., Mycologie, D., Pasteur, I., Paris, F-. (2021). Structure, Composition and Biological Properties of Fungal Extracellular Vesicles. microLife. 2, 1–13. doi: 10.1093/femsml/uqab009
Rodrigues, M. L., Nakayasu, E. S., Oliveira, D. L., Nimrichter, L., Nosanchuk, J. D., Almeida, I. C., et al. (2008). Extracellular Vesicles Produced by Cryptococcus Neoformans Contain Protein Components Associated With Virulence. Eukaryot. Cell. 7 (1), 58–67. doi: 10.1128/EC.00370-07
Scaria, V., Pasha, A. (2013). Long Non-Coding RNAs in Infection Biology. Front. Genet. 3, 1–6. doi: 10.3389/fgene.2012.00308
Sunwoo, H., Wu, J. Y., Lee, J. T. (2015). The Xist RNA-PRC2 Complex at 20-Nm Resolution Reveals a Low Xist Stoichiometry and Suggests a Hit-and-Run Mechanism in Mouse Cells. Proc. Natl. Acad. Sci. 112 (31), E4216–E4225. doi: 10.1073/pnas.1503690112
Sweta, S., Dudnakova, T., Sudheer, S., Baker, A. H., Bhushan, R. (2019). Importance of Long Non-Coding RNAs in the Development and Disease of Skeletal Muscle and Cardiovascular Lineages. Front. Cell Dev. Biol. 7, 228. doi: 10.3389/fcell.2019.00228
Telzrow, C. L., Zwack, P. J., Righi, S. E., Dietrich, F. S., Chan, C., Owzar, K., et al. (2021). Comparative Analysis of RNA Enrichment Methods for Preparation of Cryptococcus Neoformans RNA Sequencing Libraries. G3 11 (11). doi: 10.1101/2021.03.01.433483
Till, P., Mach, R. L., Mach-Aigner, A. R. (2018). A Current View on Long Noncoding RNAs in Yeast and Filamentous Fungi. Appl. Microbiol. Biotechnol. 102, 7319–7331. doi: 10.1007/s00253-018-9187-y
Ulitsky, I. (2016). Evolution to the Rescue: Using Comparative Genomics to Understand Long non-Coding RNAs. Nat. Rev. Genet. 17 (10), 601–614. doi: 10.1038/nrg.2016.85
Wallace, E. W. J., Maufrais, C., Sales-lee, J., Tuck, L. R., Natarajan, P., De Oliveira, L., et al. (2020). Quantitative Global Studies Reveal Differential Translational Control by Start Codon Context Across the Fungal Kingdom. Nucleic Acids Res. 48 (5), 2312–2331. doi: 10.1093/nar/gkaa060
Wang, K. C., Chang, H. Y. (2011). Molecular Mechanisms of Long Noncoding RNAs. Mol. Cell. 43, 904–914. doi: 10.1016/j.molcel.2011.08.018
Wang, Z., Jiang, Y., Wu, H., Xie, X., Huang, B. (2019). Genome-Wide Identification and Functional Prediction of Long Non-Coding RNAs Involved in the Heat Stress Response in Metarhizium Robertsii. Front. Microbiol. 10, 1–11. doi: 10.3389/fmicb.2019.02336
Wang, Z. F., Whitfield, M. L., Ingledue, C., Dominski, Z., Marzluff, W. F. (1996). The Protein That Binds the 3 ‘ End of Histone mRNA: A Novel RNA-Binding Protein Required for Histone pre-mRNA Processing. Genes Dev. 10 (23), 3028–3040. doi: 10.1101/gad.10.23.3028
Wang, Z., Yang, B., Zhang, M., Guo, W., Wu, Z., Wang, Y., et al. (2018). lncRNA Epigenetic Landscape Analysis Identifies EPIC1 as an Oncogenic lncRNA That Interacts With MYC and Promotes Cell-Cycle Progression in Cancer. Cancer Cell. 33, 706–720.e9. doi: 10.1016/j.ccell.2018.03.006
Williams, R. W., Rubin, G. M. (2002). ARGONAUTE1 is Required for Efficient RNA Interference in Drosophila Embryos. Proc. Natl. Acad. Sci. 99 (10), 6889–6894. doi: 10.1073/pnas.072190799
Wilusz, J. E., Jnbaptiste, C. K., Lu, L. Y., Kuhn, C., Joshua-tor, L., Sharp, P. A. (2012). A Triple Helix Stabilizes the 3 9 Ends of Long Noncoding RNAs That Lack Poly (A) Tails. Genes Dev. 26 (21), 2392–2407. doi: 10.1101/gad.204438.112
Yamashita, A., Shichino, Y., Yamamoto, M. (2016). The Long Non-Coding RNA World in Yeasts. Biochim. Biophys. Acta - Gene Regul. Mech. 1859, 147–154. doi: 10.1016/j.bbagrm.2015.08.003
Yu, C.-H., Chen, Y., Desjardins, C. A., Tenor, J. L., Toffaletti, D. L., Giamberardino, C., et al. (2020). Landscape of Gene Expression Variation of Natural Isolates of Cryptococcus Neoformans in Response to Biologically Relevant Stresses. Microb. Genom. 6 (1). doi: 10.1099/mgen.0.000319
Yu, Y., Yarrington, R. M., Chuong, E. B., Elde, N. C., Stillman, D. J. (2016). Disruption of Promoter Memory by Synthesis of a Long Noncoding RNA 113. Proc. Natl. Acad. Sci. 113 (34), 9575–9580. doi: 10.1073/pnas.1601793113
Zhang, H., Tang, K., Qian, W., Duan, C. G., Wang, B., Zhang, H., et al. (2014). An Rrp6-Like Protein Positively Regulates Noncoding RNA Levels and DNA Methylation in Arabidopsis. Mol. Cell. 54, 418–430. doi: 10.1016/j.molcel.2014.03.019
Zhang, P., Wu, W., Chen, Q., Chen, M. (2019). Non-Coding RNAs and Their Integrated Networks. J. Integr. Bioinform. 16, 1–12. doi: 10.1515/jib-2019-0027
Keywords: lncRNA, cryptococcus, fungi, pathogenesis, genome
Citation: Kalem MC and Panepinto JC (2022) Long Non-Coding RNAs in Cryptococcus neoformans: Insights Into Fungal Pathogenesis. Front. Cell. Infect. Microbiol. 12:858317. doi: 10.3389/fcimb.2022.858317
Received: 19 January 2022; Accepted: 18 February 2022;
Published: 17 March 2022.
Edited by:
Robert T. Wheeler, University of Maine, United StatesReviewed by:
Andrew Alspaugh, Duke University, United StatesElizabeth R. Ballou, University of Birmingham, United Kingdom
Xiaorong Lin, University of Georgia, United States
Copyright © 2022 Kalem and Panepinto. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Murat C. Kalem, muratcankalem@gmail.com; John C. Panepinto, jcp25@buffalo.edu