Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Fungal Biol., 30 July 2025

Sec. Fungal Genomics and Evolution

Volume 6 - 2025 | https://doi.org/10.3389/ffunb.2025.1633301

Orthrus: a Pumilio-family gene involved in fruiting body and dark stipe development in Coprinopsis cinerea

Benedek SzathmriBenedek Szathmári1Balzs Blint&#x;Balázs Bálint1†Botond HegedüsBotond Hegedüs1Mt VirghMáté Virágh1Zhihao Hou,Zhihao Hou1,2Xiao-Bin LiuXiao-Bin Liu1Hongli WuHongli Wu1Csenge FldiCsenge Földi1Julien Gagneur,Julien Gagneur3,4Johann PromeuschelJohann Promeuschel3rpd Csernetics&#x;Árpád Csernetics1‡Lszl G. Nagy,&#x;*László G. Nagy1,5‡*
  • 1Synthetic and Systems Biology Unit, Institute of Biochemistry, HUN-REN Biological Research Centre Szeged, Szeged, Hungary
  • 2Doctoral School of Biology, Faculty of Science and Informatics, University of Szeged, Szeged, Hungary
  • 3School of Computation, Information and Technology, Technical University of Munich, Garching, Germany
  • 4Munich Center for Machine Learning, Munich, Germany
  • 5Korea University, Seoul, Republic of Korea

Fruiting bodies of mushroom-forming fungi (Agaricomycetes) are complex multicellular structures whose formation is regulated by a developmental program that dynamically responds to environmental changes, such as light intensity. However, the genetic architecture and regulation of this developmental program are poorly known. Here, we characterize a novel Pumilio family gene, ort2, which influences fruiting body development, particularly the formation of dark stipes, a light-dependent alternative developmental trajectory. Phylogenetic analysis of this RNA-binding protein family in fungi revealed a distinct subfamily structure, with high conservation of each subfamily within Agaricomycetes. Reverse genetics experiments in the model species Coprinopsis cinerea revealed that ort2 disruptants produced fruiting bodies, but were deficient in dark stipe formation, whereas the overexpression mutants produced significantly more dark stipes. The gene was named after Orthrus, the two-headed dog of classical mythology, based on rare but reproducible branching fruiting body phenotypes observed upon overexpression. Our findings reveal fruiting-related functions for ort2, a novel conserved RNA-binding protein, and may serve as a novel entry point for understanding the molecular basis of dark stipe development.

Introduction

Fungal fruiting bodies are sexual reproductive structures that facilitate spore dispersal (Hibbett, 2007). They are complex multicellular (Knoll, 2011) – with a three-dimensional structure, in which not all the cells are in direct contact with the environment. They exhibit all key complex multicellular traits, such as intercellular communication, and a development that is orchestrated by a genetically encoded program (Kües, 2000; Nagy et al., 2018, 2023). They represent the most morphologically complex structures that fungi produce and have been a subject of intense research interest (Virágh et al., 2022).

Fruiting body development starts with the appearance of a dense hyphal aggregate, the primary hyphal knot, which appears in darkness, and turns into a secondary hyphal knot in the presence of light (Kües and Navarro-González, 2015). Hyphal knots represent the first three-dimensional developmental stage that fulfills all criteria of complex multicellularity (Knoll, 2011). Tissue differentiation primarily occurs during the primordium stages (P1–P5) and requires alternating light conditions (Kües, 2000; Kües and Navarro-González, 2015). While morphological aspects of the developmental process are well-described, their genetic underpinnings are poorly known, even though the development of fruiting bodies has been a central question in mycology.

Recent genomic and transcriptomic studies have significantly advanced our understanding of fruiting body morphogenesis in Agaricomycetes and complemented data from classical genetic studies (Kües, 2000; Ohm et al., 2011). Fruiting body development appears to utilize a deeply conserved yet flexible genetic toolkit (Krizsán et al., 2019; Nagy et al., 2023). Several genes, including those that encode transcription factors, cell wall remodeling enzymes, and signaling components, have been predicted to contribute to fruiting body development; however, the regulatory architecture remains incompletely resolved, since experimental characterization of most genes is lacking (Földi et al., 2024). Recent studies uncovered several key regulators of fruiting body development in model species such as Coprinopsis cinerea (Sakamoto et al., 2024; Masuda et al., 2016; Muraguchi et al., 2008; Muraguchi and Kamada, 1998) and Schizophyllum commune (Pelkmans et al., 2017; Beijen et al., 2025; Vonk et al., 2024), or commercially produced taxa, e.g., Lentinula edodes (Miyazaki et al., 2004) and Agaricus bisporus (Pelkmans et al., 2016). Most of these are transcription factors that bind DNA in a sequence-specific manner and regulate the expression of genes. However, the regulation of development likely also involves post-transcriptional and post-translational layers. The genes involved in these are currently unknown in mushroom-forming fungi.

In typical agaricomycete fruiting bodies, development proceeds from primordial stages to sporulation through stipe elongation and cap expansion. A characteristic, environment-dependent alteration of this developmental trajectory results in the formation of dark stipes or etiolated fruiting bodies (Tsusué, 1969). They are best studied in C. cinerea, which forms dark stipes in darkness; however, the phenomenon exists in several other species, often referred to by different names. Dark stipes are elongate structures that bear underdeveloped pileus and stipe tissues at their tip. They elongate via growth of the organ called nodulus (basal plectenchyma, basal shaft) (Kües, 2000; Chaisaena, 2009). Developmental stages up to P4 primordium can undergo etiolation when transferred to constant darkness (Subba, 2022). At the tip, the development of caps and stipes is arrested in the absence of light (Kües, 2000; Terashima et al., 2005). Upon light exposure, tissues at the tip can resume development and grow into a mature fruiting body (Tsusué, 1969); however, the original primordium is frequently aborted, and secondary primordia emerge from the dark stipe (Kües, 2000). In addition to darkness, elevated CO2 levels can also promote the formation of etiolated fruiting bodies (Subba, 2022).

High CO2 concentration and darkness trigger the formation of elongate fruiting bodies in other species, too. For example, Flammulina velutipes, F. filiformis (Yan et al., 2019), and Pleurotus eryngii (Zadražil, 1975) are cultivated under elevated CO2 conditions, to promote elongation and repress cap development. For these popular species, the elongated forms with smaller caps (called pinhead mushrooms in Flammulina species) are generally regarded to be of higher quality (Yan et al., 2019).

In earlier terminology, C. cinerea dark stipes are called pseudorhizae (Kües, 2000), since in nature, they appear as root-like structures within the substrate. The primary function of pseudorhizae is to push cap and stipe initials towards the light from below the surface of the substrate (e.g., within dung) (Buller, 1958; Clémençon, 2012; Kües, 2000). In a broader perspective, the developmental program of dark stipes may represent an adaptation to escape enclosed, dark places, such as the interior of the substrate or cavities. The term pseudorhiza is widely used in mycology to describe similar, substrate-embedded, elongated forms in several species (Buller, 1958). Some of them (e.g., Hymenopellis radicata, Mycena galericulata) are known to be functionally analogous (Buller, 1958) to C. cinerea pseudorhizae, but their developmental homology has not been proven.

Our current understanding of the molecular biology underlying dark stipe development is incomplete, limited to photoreception. The dst1 gene is a WC1 photoreceptor, while dst2 encodes a FAD/FMN-containing dehydrogenase (Subba, 2022). The inactivation of both dst genes results in a ‘blind’ phenotype, which refers to the formation of dark stipes even under standard (i.e., alternating light) fruiting conditions (Chaisaena, 2009).

Pumilio (Pum/PUF) constitutes a eukaryotic family of RNA-binding proteins with extensively studied members across animals, plants, and fungi alike. Pumilio proteins contribute to a broad spectrum of biological processes, such as embryonic patterning, stem cell maintenance, along with ribosome and mitochondrial biogenesis, among others (Nishanth and Simon, 2020). They play a critical role in post-transcriptional regulation, typically by destabilizing mRNA targets bound in a sequence-specific manner. mRNA destabilization is mediated via mechanisms such as inhibition of translation initiation, deadenylation, or microRNA activity (Nishanth and Simon, 2020). PUF RNA regulatory networks evolve dynamically; comparative analyses of predicted and experimentally validated targets across fungal species demonstrated recurrent rewiring (Wilinski et al., 2017).

In the fungal kingdom, the Pumilio proteins of Saccharomyces cerevisiae are the most extensively characterized. Yeast PUF proteins bind distinct subsets of transcripts, with defined recognition motifs (Gerber et al., 2004). Puf3p, one of the best-studied members of the family, localizes to the outer mitochondrial membrane and plays a role in mitochondrial biogenesis and motility (García-Rodríguez et al., 2007). The deletion of PUF3, as well as hyperphosphorylation of the encoded protein (Bhondeley and Liu, 2020) increases, whereas its overexpression decreases mitochondrial function, indicating that Puf3p is a negative regulator. Additionally, Puf3p was also reported to contribute to adaptation to metabolic changes via switching the balance in translational flux between mitochondrial and cytosolic ribosome biogenesis (Wang et al., 2019). The ortholog in Schizosaccharomyces pombe, Puf3, was shown to function via the Ccr4-Not complex (Webster et al., 2019). An analysis comprising 42 fungal (four basidiomycete) species indicates the mitochondrion-related function is an apomorphy of the Saccharomycotina clade (Jiang et al., 2010). The study of Wilinski et al. (2017) revealed that the Neurospora crassa Puf3p ortholog binds a similar motif to that of S. cerevisiae Puf3, but with a cytosine in a central position in half of the target sequences, instead of an adenine. In Cryptococcus neoformans, the Puf3 ortholog Pum1 serves as a pleiotropic regulator of sexual development, including filamentation and basidium maturation (Liu et al., 2018). In this species, the maintenance of the vegetative phase requires the autorepression of PUM1 (Kaur and Panepinto, 2016), with a longer, yeast phase-specific transcript bearing the recognition motif of its own product. Pumilio proteins have been reported to have dynamic expression patterns during fruiting body formation of complex multicellular Agaricomycetes (Krizsán et al., 2019; Merényi et al., 2020; Nagy et al., 2023); however, whether these proteins contribute to complex development has not been experimentally studied until now.

In this study, we analyze Pumilio protein-encoding genes in C. cinerea and select one member, referred to as ort2, that is orthologous to S. cerevisiae PUF3 and C. neoformans PUM1 for functional investigation. We show that Pumilio proteins form distinct and conserved subfamilies in the Agaricomycetes, and that ort2 influences fruiting body development in this species. We detected the largest impact of ort2 on dark stipe development and attempted to decipher the mechanistic bases by which it influences development.

Materials and methods

Strain, buffers and culture conditions

We used the self-compatible, para-aminobenzoic acid auxotrophic AmutBmut1 pab1-1 #326 strain of C. cinerea for genetic modification and as a control in all experiments. The fungus was maintained on solid YMG (yeast-malt-glucose) medium. All media and buffers used for culturing and transformation (YMG, minimal medium, regeneration medium, top agar, PEG/CaCl2 solution, MM and MMC buffers) were prepared as described by Dörnte and Kües (2012). Oidiation was induced by culturing at 37°C under constant white light. After transformation, we selected for the complementation of PABA auxotrophy using synthetic minimal medium.

Phylogenetic analyses

Homologs of ort2 were retrieved by blastp (e < 10-5) and an InterPro search (search terms: IPR001313, IPR033133) from a proteome dataset comprising 13 agaricomycete species, two additional basidiomycetes (Cryptococcus neoformans, Ustilago maydis) (Nagy et al., 2023) complemented by two ascomycete yeasts, Saccharomyces cerevisiae and Schizosaccharomyces pombe. The proteome of Phanerochaete chrysosporium is included in two versions. For our model species, C. cinerea version 2 was used. We refer to C. cinerea proteins based on the annotation of Hegedüs et al. (2025) (https://mushroomdb.brc.hu/); the IDs we supply can be found in the JGI annotation as gene (model) names (https://mycocosm.jgi.doe.gov/Copcin2/Copcin2.home.html). Homologs fitting the HMM of IPR040000 were excluded, as these are orthologs of the nucleolar protein Nop9 of S. cerevisiae (C. cinerea protein: 438228). Retained homologs were aligned by PRANK v.170427 (Löytynoja, 2014) by using the -F parameter and a maximum likelihood guide tree inferred in preliminary analyses by RAxML version 8.2.12 (Stamatakis, 2014), from a multiple sequence alignment generated by Mafft-linsi v7.453 (Katoh and Standley, 2013). The PRANK alignment was used for tree inference and bootstrapping by RAxML version 8.2.12, under the WAG model with gamma-distributed rate heterogeneity. Branch support was assessed by 500 non-parametric rapid bootstrap replicates. Trees were visualized in FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree).

Design of molecular tools

We intended to perform gene deletion through Cas9 cleavage followed by homologous repair, and overexpression via ectopic integration. We designed two crRNAs and a plasmid (p.Δort2) as a repair template, and another plasmid (p.ort2OE) was constructed to generate overexpressing strains (Supplementary Figure S2). The plasmids, primers and crRNAs were designed using SnapGene version 4.3.11 (https://www.snapgene.com/). The sequences of primers and crRNAs are listed in Supplementary Table S1. Both plasmids possess a backbone derived from pUC19 and carry the positive selection marker gene pab1-1 with its promoter and terminator regions (CcPAB1) amplified from the pMA412 vector (Stanley et al., 2014). p.Δort2 carries the marker gene between homologous arms (1 kb sequences upstream and downstream of the target sequence). p.ort2OE carries the ort2 gene driven by the promoter of a heat shock protein gene, Cc.hsp3 (ID: 423239), given that the promoter of its ortholog has been shown to be useful for overexpression in Schizophyllum commune (Ohm et al., 2011).

Plasmid assembly and preparation

Fragments were amplified by Phusion Plus DNA polymerase (Thermo Fisher Scientific) using overlapping Gibson primers. Plasmids were assembled using the Gibson Assembly Cloning Kit (New England Biolabs, USA), following the manufacturer’s instructions. The ligation mixture was introduced into competent NEB5α cells following the protocol of Rabe and Cepko (2020), and the plasmids were purified from the cultures of colonies for which the presence and proper composition of the plasmids had been confirmed, using the NucleoBond Xtra Midi kit (Macherey-Nagel, Germany).

Protoplasting oidia

Oidia were harvested from 5-day-old cultures grown under constant light, by washing with distilled water, filtering through a 40 µm cell strainer, and pelleting by centrifugation (1000 g, 10 min). Pellets were washed with 5 ml MM buffer, pelleted again, and resuspended in 1 ml sterile-filtered MM buffer containing ~20% w/v VinoTaste Pro (Novozymes, Denmark). Protoplasting was carried out by shaking oidia at 70–90 rpm in an incubator at 36°C for 3–4 hours. Once the majority of oidia turned into spherical protoplasts, the process was terminated by adding 5 ml MMC buffer, and the protoplasts were collected (1000 g, 15 min).

PEG-mediated transformation of C. cinerea protoplasts

The pelleted protoplasts were resuspended in 100 µl MMC buffer. We added 10 µg plasmid DNA; in the case of the knockout transformation, the two RNP complexes [assembled as described by Pareek et al. (2022)]; and for both transformations, 30 µl PEG-CaCl2 solution. The cells were incubated on ice for 30 min, then 0.5 ml PEG-CaCl2 solution was added, and the mixture was incubated at room temperature for 10–20 minutes. After adding 20 ml MMC buffer, the solution was gently mixed, and the cells were pelleted (1000 g, 15 minutes). The cells were resuspended in 1 ml MMC buffer, mixed with top agar, and poured onto regeneration medium. Colonies were collected onto selective minimal medium plates. Transformants were checked by PCR, amplicon sequencing, and RT-qPCR. To obtain DNA template for the PCR reactions, vegetative mycelium was lysed thermally as described by Dörnte and Kües (2012). Amplicon sequencing was conducted commercially (Microsynth, Switzerland) after purification using the NucleoSpin Gel and PCR Clean-up kit (Macherey-Nagel, Germany). To purify mixed-genotype cultures, we prepared single-spore cultures from oidia.

RNA purification and RT-qPCR

Cultures grown on cellophane were homogenized in a mortar with liquid nitrogen. RNA extraction was performed using the Quick-RNA MiniPrep kit (Zymo Research, USA) following the manufacturer’s instructions. We screened the overexpressing mutant candidates by RT-qPCR. The colonies of the transformants and the wild-type controls were exposed to 42°C heat shock for 2 hrs before sample collection. Reverse transcription was carried out with the RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific), and qPCR reactions were performed using SsoAdvanced Universal SYBR Green Supermix (BioRad), following the manufacturers’ instructions. RNA extracted from the wild-type strain was used as a control, and the β-tubulin gene (protein ID: 393528) was chosen as the reference gene. Expression changes were calculated using the formula FC = 2−ΔΔCt (Livak and Schmittgen, 2001). qPCR was conducted in six technical replicates: two reactions per experiment, repeated three times starting from the same RNA extracts.

Phenotyping of fruiting bodies

Since light- and dark-grown colonies of C. cinerea produce fruiting bodies of different morphology, we investigated fruiting under both conditions. Strains were inoculated onto ‘half-sugar’ YMG (containing 2 g/L glucose) (Muraguchi et al., 2015), incubated in darkness at 28°C for six days, then half of the plates were transferred to alternating light conditions (12 hrs light, 12 hrs dark; Panasonic FL40SS-ENW/37 light bulb), while the rest were kept in darkness. Both groups were further incubated at 28°C for two weeks. After several smaller-scale phenotyping experiments, fifty plates of the deletion mutants and one hundred of the selected overexpression strains were phenotyped and compared to a wild-type control group. For the overexpression mutants, after the additional two-week incubation, we quantified fruiting, along with making morphological observations. We counted etiolated fruiting bodies (dark stipes) longer than 1 cm per Petri dish formed by dark-grown colonies, and structures that had at least started to elongate (P4 and older) in the case of plates exposed to alternating light. Phenotyping experiments were also conducted in jars containing 150 ml ‘half-sugar’ YMG with 6 g/L agar, with the same protocol.

Assessment of relative mtDNA abundance

Following the method described by Venegas et al. (2011) and Galeota-Sprung et al. (2022), we performed qPCR on a nuclear (β-tubulin gene, JGI protein ID: 393528) and a mitochondrially encoded gene (NADH dehydrogenase subunit 1, protein ID: 384263), using DNA templates extracted from an ort2OE, four Δort2 mutants, and the wild-type strain, by using the Blood & Cell Culture DNA Maxi Kit (Qiagen, Germany). We used five extracts from all the strains; and for each sample, two technical replicate qPCR reactions were performed. qPCR was carried out in the same way as described above.

Functional studies

Research Group Julien Gagneur (Munich Centre for Machine Learning) provided us with a list of putative Ort2 recognition motifs and mRNA targets, predicted by their species-aware DNA language model (Karollus et al., 2024). We utilized the hits belonging to the most significant motif. The list contains 772 Okayama 7 genes; 655 remain when converted to AmutBmut identifiers. A Gene Ontology (GO) term enrichment analysis was performed on the AmutBmut gene set. A FIMO (Grant et al., 2011) search was also conducted for UTRs of C. cinerea (Hegedüs et al., 2025), using the S. cerevisiae Puf3p recognition motif PWM as described by Gerber et al. (2004). We aimed to select the top 100 most significant hits; however, due to the identical p-values among several entries, a total of 116 were included. We also generated three mitochondrial gene sets. Okayama 7 genes with a mitochondrial GO annotation were downloaded from QuickGO (https://www.ebi.ac.uk/QuickGO) and converted to AmutBmut identifiers. Mitochondrial gene sets were also obtained based on predicted subcellular localization. We employed TargetP 2.0 (Almagro Armenteros et al., 2019) and DeepLoc 2.0 (Thumuluri et al., 2022) to identify proteins potentially localizing to mitochondria. Hits with 0.7 or higher probability scores were considered significant. An overlap analysis was conducted on the five gene sets. An Euler diagram was drawn using the eulerr R package. For predicting the subcellular localization of Ort2, we used DeepLoc 2.0 (Thumuluri et al., 2022) as implemented at https://services.healthtech.dtu.dk/ and WolfPSort (Horton et al., 2007) at https://wolfpsort.hgc.jp. Both software were run with default parameters and the setting ‘slow’ for DeepLoc 2.0 and ‘fungi’ for WolfPsort.

Results

Pumilio proteins form five conserved subfamilies

Through BLASTp and InterPro searches, we collected the Pumilio protein sequences from the proteomes of 13 agaricomycete species, two additional basidiomycete species (Cryptococcus neoformans, Ustilago maydis) (Nagy et al., 2023), and two ascomycete yeasts, Saccharomyces cerevisiae and Schizosaccharomyces pombe. We identified 4 to 7 Pumilio proteins in the examined species of Basidiomycota, which is similar to copy numbers reported for Ascomycota (Matsushima et al., 2007). In Coprinopsis cinerea, six proteins that possess Pumilio-family RNA-binding repeat domains (IPR001313) were found. C. cinerea 438228 was excluded from the analyses, since its sequence is highly divergent, contains another domain (IPR040000), and, as an ortholog of the nucleolar Nop9 of Saccharomyces cerevisiae, it might function differently. Nop9 orthologs also formed a distant subfamily in previous studies (Najdrová et al., 2020).

The Pumilio family formed 5 clearly discernible and strongly supported (maximum likelihood bootstrap > 70%) subfamilies, each of which contained S. cerevisiae and Sch. pombe sequences, indicating conservation of major Pumilio subfamilies retained across Asco- and Basidiomycota (Figure 1; Supplementary Figure S1). Each subfamily harbors a single protein of each agaricomycete species, with few exceptions (e.g., two Pleurotus ostreatus paralogs in clade puf4-puf5). Subfamily topologies mostly reflect the species tree (Naranjo‐Ortiz and Gabaldón, 2019): agaricomycete proteins group together, while more basal branches represent the two ascomycetes and the two basidiomycete species outside Agaricomycetes. This indicates highly conserved evolutionary trajectories and low overall divergence within subfamilies. Subfamilies differ from each other considerably, as indicated by the long branches separating them. C. cinerea 354109 (CopciAB_ort2) grouped with Puf3 of Sch. pombe, Puf3p of S. cerevisiae, and Pum1 of Cryptococcus neoformans (Figure 1), whose function has been investigated (Webster et al., 2019; Liu et al., 2018), especially that of S. cerevisiae Puf3p (García-Rodríguez et al., 2007).

Figure 1
Phylogenetic tree highlighting different clades in various colors: green for “mcp2-mpf1” clade, purple for “puf3” clade, blue for “puf1/puf2” clade, orange for “puf6” clade, and gray for “puf4/puf5” clade. Surrounding graphs show data variability across conditions VM, HK, P1, P2, YFB, and MFB with species names like CopciAB, Schpo, Sacce, and Cryne annotated in matching colors.

Figure 1. Unrooted maximum likelihood gene tree of fungal Pumilio proteins and expression dynamics of Coprinopsis cinerea Pumilio genes. The phylogenetic analysis comprises protein sequences from 17 fungal species. Pumilio proteins form five conserved subfamilies. Subfamilies are named after Sch. pombe proteins. Clade colors match those of the expression diagrams of the corresponding C. cinerea genes. In each clade, only terminal nodes belonging to C. cinerea, Cr. neoformans, S. cerevisiae, and Sch. pombe are labeled by JGI identifiers or common names. For more details, see Supplementary Figure S1. For the expression plots, we used Krizsán et al.’s RNA-Seq dataset (2019). The developmental stages are as follows: VM, vegetative mycelium; HK, hyphal knot; P1, stage 1 primordium; P2, stage 2 primordium; YFB, young fruiting body; MFB, mature fruiting body. In earlier stages (VM to P2), bulk expression values are shown (indicated by filled circles), in more developed stages (YFB and MFB), tissue-specific expression is plotted. In the YFB stage, squares, triangles, and diamonds indicate expression in stipe, cap, and gills, respectively. In the MFB stage, squares similarly indicate expression in the stipe tissue, while empty triangles mark combined cap and gill expression. Expression values are expressed as CPM on all plots.

We inspected the expression dynamics of the five C. cinerea genes encoding Pumilio proteins in Krizsán et al.’s RNA-Seq dataset (2019) (Figure 1). The expression of C. cinerea 354109 (CopciAB_ort2) peaks in the hyphal knot (HK), the earliest stage of fruiting body formation (Kües and Navarro-González, 2015), suggesting it has functions in the early developmental events. We note that Xie et al.’s data (2020) indicate a maximum expression in vegetative mycelium; however, both datasets agree that its expression decreases as the developmental program progresses. For the phylogenetic and expression-related reasons discussed above, we selected C. cinerea 354109, later named ort2, for experimental characterization.

Disruption of ort2 affects fruiting body development

Using the RNP-mediated CRISPR-Cas9 genome editing technique (Pareek et al., 2022), we generated four Δort2 disruptant strains (see the plasmid in Supplementary Figure S2a). The lack of ort2 expression was demonstrated by RT-qPCR (Supplementary Table S2). Though we used two crRNAs, instead of complete deletion, we could achieve only gene disruption. The disruptant genomes were not screened for off-target integrations, given that the phenotypic similarity of the independent disruptant strains implies the observed characters likely result from the intended disruption of ort2 rather than off-target effects.

Phenotyping of the four Δort2 disruptants revealed no differences in vegetative growth patterns and oidiation relative to the parent strain (not shown), whereas fruiting body and dark stipe development seem to be impaired to a degree (Figure 2), the latter process considerably. Grown under alternating light conditions in jars, stage 1 primordia appeared at least 3 days later on disruptant cultures than on the wild-type ones, and no Δort2 fruiting bodies reached maturity (i.e., failed to produce basidiospores and deliquesce). We counted structures whose stipes had already started to elongate, including P4 and P5 stage primordia (Kües and Navarro-González, 2015), as well as young and mature fruiting bodies, together on four wild-type cultures and 16 Δort2 jar cultures (Figure 2a). On average, on wild-type cultures 10.5; while on Δort2 cultures 1.875 such structures formed, which is a significant difference (p = 0.0009, two-sample t-test); although we acknowledge the limited power of the analysis due to the low number of data points and the presence of outliers. We conducted a larger-scale experiment with 20 wild-type and 20 mutant jar cultures. The wild-type cultures formed 16.05 ± 9.56 fruiting bodies (P4 to MFB), while none of the disruptants produced any fruiting bodies or initials, indicating the decreased readiness of the mutants for fruiting compared to the parent strain. No consistent difference was observed in the tissue structure of primordia; however, many of them exhibited diverse abnormalities (Figure 2b). When cultured under constant darkness on plates, wild-type cultures produced dark stipes, whereas the Δort2 cultures did not. Across four wild-type jar cultures, a total of 34 dark stipes were counted, approximately 20 of which exceeded 1 cm in length/height. In contrast, on 16 Δort2 jar cultures, only three dark stipe-like structures were observed altogether (all measuring less than 1 cm), along with some rudimentary structures that failed to elongate. This indicates that Δort2 strains are capable of producing dark stipes, though to a considerably lower extent.

Figure 2
Violin plot on the left compares the number of fruiting bodies between WT and Δort2, showing a significant difference (p = 0.0009). Images on the right depict eight different samples of fruiting bodies, each with a scale bar indicating size, ranging from 1 millimeter to 2 millimeters.

Figure 2. Phenotypic characterization of Δort2 disruptants. (a) Violin plot: total number of fruiting bodies in stages from P4 to MFB on light-grown jar cultures of the wild-type and Δort2 strains. (b) Cross-sections of wild-type (left panel) and Δort2 (right panel) primordia.

Overexpression of ort2

Utilizing an ectopically integrating plasmid (Supplementary Figure S2b), additional copies of ort2 were introduced into the genome of the AmutBmut strain, under the control of the hsp3 heat shock promoter. RT-qPCR was performed to measure ort2 expression in mycelia grown at 28°C (Supplementary Table S3). We selected the strain with the highest expression change (ort2OE(10); FC = 4.14 ± 2.79), along with another overexpression strain (ort2OE(7); FC = 1.28 ± 0.65) for preliminary phenotyping. After confirming the phenotypes were consistent across the independent mutants in preliminary experiments, ort2OE(10) was characterized further in larger-scale experiments.

According to our observations, the overexpression of ort2 did not affect vegetative growth characteristics and oidiation, however, we observed an effect on fruiting. Under alternating light, the ort2OE strains are able to complete the fruiting program; from hyphal knots to mature fruiting bodies, all developmental stages were observed, as well as sclerotia and dark stipes. For plates grown under alternating light conditions, we counted structures whose stipes had already started to elongate (P4 to MFB), together. On 50-50 Petri dish cultures, we found that the number of fruiting bodies reaching stage 4 primordium or older was similar between the wild type and overexpression strains (p = 0.2354, two-sample t-test) (Figure 3a). Mature, deliquescent fruiting bodies were not observed in the wild-type strain during the phenotyping, whereas, on 50 ort2OE Petri-dish cultures, altogether 8, and on 4 cultures grown in jars, 6 fruiting bodies were able to complete development by the end of the experiment. In the case of dark-grown colonies, we counted etiolated fruiting bodies taller than 1 cm, as a proxy for the number of dark stipes. The overexpression mutant produced significantly more dark stipes than the wild-type strain (p = 7.23 × 10-9, two-sample t-test) (Figures 3b, c). The average number of these structures per Petri dish was approximately 7.5 times higher in the case of the ort2OE mutant than in the wild-type strain. This, combined with the disruptant phenotypes, suggests that ort2 influences the development of dark stipes.

Figure 3
Panel a shows a violin plot comparing the number of fruiting bodies between WT and ort2OE groups with a non-significant p-value of 0.2354. Panel b displays another violin plot, comparing the number of dark stipes over one centimeter between the same groups, showing a significant p-value of 7.23 x 10^-9. Panel c presents two agar plates with fungal growth, labeled WT and ort2OE. Panel d includes close-up images of dark stipe characteristics for WT and ort2OE, with visible morphological differences.

Figure 3. Phenotypic characterization of the ort2OE mutants. (a) Violin plots illustrating the total number of fruiting bodies (P4 to MFB) on light-grown Petri dish cultures of the wild-type and ort2OE strains. (b) Violin plots of the number of dark stipes exceeding 1 cm in length, formed by the wild-type and ort2OE strains. (c) Comparison of WT and ort2OE dark-grown cultures. (d) An overexpression-related abnormality; the reason we chose ort2 (two-headed Orthrus) as a name for C. cinerea 354109. The branched fruiting bodies were more pronounced in dark stipes (pictured). For other overexpression-related abnormalities, see Supplementary Figure S3.

In addition, the ort2 overexpression strains exhibit diverse developmental abnormalities, which occur at low frequencies (and therefore cannot be subjected to statistical tests) (Supplementary Figure S3). For dark-grown colonies, these phenotypes include the intertwining of etiolated fruiting bodies, forming amorphous structures (Supplementary Figure S3d), and the formation of robust prostrate dark stipes. Some of these reached lengths of up to 20 cm, nearly encircling the dish (Supplementary Figures S3c-d), a phenomenon that is uncommon in the wild-type strain. Branched fruiting body morphologies, occurring on both dark-grown colonies and those exposed to alternating light, constitute another low-frequency developmental abnormality observed upon ort2 overexpression (Figure 3d; Supplementary Figure S3a). We documented fruiting bodies with branched stipes, as well as ones bearing protrusions on their caps (or tips). The branched fruiting bodies and snake-like prostrate dark stipes led us to name C. cinerea 354109 ort2 after Orthrus, the two-headed, snake-tailed dog from classical mythology.

Mechanistic insights into the function of ort2 are limited

We attempted obtaining mechanistic insights into the role of Ort2 in C. cinerea, by (1) functional enrichment analyses of putative targets inferred by a large language model (Karollus et al., 2024) and (2) by specifically addressing whether it is involved in mitochondrial biogenesis, a function that was postulated to be a Saccharomycotina-specific innovation (Jiang et al., 2010).

Using a language model (Karollus et al., 2024), reconstruction probability arrays were generated for the C. cinerea Okayama 7 genome (Stajich et al., 2010) (300 bps downstream of STOP codons, since Pumilio proteins primarily bind the 3’-UTRs of transcripts (Nishanth and Simon, 2020). Three motifs were identified that significantly matched the S. cerevisiae Puf3p motif (using TF-MoDISco and TOMTOM). The most significant motif (Figure 4a) was found in 772 C. cinerea Okayama 7 transcripts, 655 of which were found in the AmutBmut strainThe information content of the reconstruction probabilities of motif sites was significantly higher than that of sites containing the shuffled motif (p = 4.02×10-166, Mann-Whitney U test), suggesting the DNA language model recognizes the specific nucleotide arrangement defining the motif, not just general sequence composition. This supports the motif’s robustness and increased confidence in future use. It is worth noting that the central position in the predicted C. cinerea motif – where it differs from the S. cerevisiae Puf3p motif (A) – can similarly be occupied by a cytosine in the experimentally determined binding site of the Neurospora crassa ortholog (Wilinski et al., 2017) (Figure 4a), supporting our prediction. We analyzed GO enrichment within the 655 C. cinerea AmutBmut target genes. This analysis did not reveal a clear signal, with only five significantly enriched terms (adjusted p-value < 0.05, Supplementary Table S4). The cause of the lack of signal in this dataset may be manifold, though we note a clear functional enrichment was not found for N. crassa either, perhaps meaning that these proteins bind functionally diverse transcripts. Nevertheless, our GO results did not comprise terms related to the mitochondrion, in agreement with previous results on filamentous fungi (Jiang et al., 2010).

Figure 4
Panel A shows sequence logos of DNA motifs with nucleotide heights indicating their relative frequency. Panel B is a bar chart displaying relative mitochondrial DNA abundance for different genetic variants, with ort2OE showing the highest level. Panel C is a Venn diagram illustrating the overlap of different datasets, labeled mt_GO, mt_DeepLoc, mt_targetP, targets_LM, and targets_FIMO, with various intersections quantified.

Figure 4. In silico prediction of Ort2 target genes. (a) S. cerevisiae Puf3p recognition motif (based on (Gerber et al., 2004)) (left) and Ort2 recognition motif predicted by the language model (right). Nucleotide frequencies of the aligned sites were transformed to information content. (b) qPCR measurement of mtDNA levels, relative to a nuclear gene, in the Δort2 and ort2OE strains, compared to a wild-type control. The relative mtDNA abundance of the wild type corresponds to 1 on the plot. Due to the necessity of averaging during the calculation of FC values, standard deviations of FC values cannot be assessed. See raw data in Supplementary Table S5. (c) Euler diagram of five C. cinerea protein sets: those with mitochondrial GO annotations (mt_GO), proteins with mitochondrial localization predicted by DeepLoc 2.0 (mt_DeepLoc), proteins with mitochondrial localization predicted by TergetP 2.0 (mt_TargetP), Ort2 targets predicted by Karollus et al.’s language model (targets_LM), and Ort2 targets predicted using FIMO (targets_FIMO). See gene sets in Supplementary Table S6. Euler plots are model-based, and the model might not include low-cardinality intersections. For the precise overlap analysis, see Supplementary Figure S4.

To verify this finding, we first compared mtDNA abundance in the Δort2, ort2OE, and wild-type strains, using quantitative PCR targeting a nuclear (β-tubulin) and a mitochondrially encoded gene (NADH dehydrogenase subunit 1, Figure 4b; Supplementary Table S5). If Ort2 regulated mitochondrial biogenesis, one would expect overexpression and deletion to have opposite effects on mtDNA abundance. However, no such differences were observed. The measured mtDNA levels in the disruptants were not significantly different from those in the wild type. Although the selected overexpression strain displayed a more pronounced difference, this change occurred in the opposite direction than one would expect based on the role of S. cerevisiae Puf3p as a negative regulator of mitochondrial biogenesis (García-Rodríguez et al., 2007).

As an alternative approach, we compared candidate Ort2 targets inferred by in silico methods to mitochondrial proteins of C. cinerea (1) Gene Ontology (GO) annotations, (2) subcellular protein localization predicted by DeepLoc 2.0, and (3) by TargetP 2.0. The three mitochondrial gene sets considerably overlap (Figure 4c), validating these predictions. Mitochondrial protein lists are provided in Supplementary Table S6. However, when these gene sets were compared with 655 putative Ort2 targets predicted by the language model, and with 116 C. cinerea mRNA targets predicted with FIMO (using the S. cerevisiae Puf3p recognition motif as query), we obtained small overlaps (Figure 4c; Supplementary Table S6). Only one comparison, between the set of genes with mitochondrial GO annotation and the list of Ort2 targets predicted by the language model, revealed a significant overlap: 14.77% of the LM hits are in the GO set (p < 10-5, Fisher’s, exact test).

Prediction of the subcellular localization of ort2 using in silico tools suggested that ort2 is cytoplasmic and nuclear. DeepLoc 2.0 yielded scores of 0.7543 and 0.5299 for cytoplasmic and nuclear localization, respectively, with much lower (>0.09) scores for other localizations. WolfPsort ranked a nuclear localization first (score: 13), followed by cytoplasmic (6), mitochondrial (4), and peroxisomal (3). We hypothesize that ort2 may shuttle between the cytoplasm and the nucleus, possibly depending on posttranslational modifications. This is consistent with the behavior of other Pumilio family proteins (Gerber et al., 2004). Overall, we interpret these results as an absence of evidence for a mitochondrial function of Ort2 is in line with the data of Jiang et al. (2010) on filamentous Ascomycota. It remains to be tested whether the predicted mRNA targets are biologically significant and how these contribute to fruiting and dark stipe development remains unclear.

Discussion

In this study, we examined the Pumilio family of RNA-binding proteins in mushroom-forming fungi and addressed the role of a conserved member of the family, ort2, selected based on its expression patterns during mushroom formation. Our experiments show that ort2 influences both fruiting in the light and dark stipe development, with a stronger effect on the latter. Since no effect on vegetative morphology was observed, the role of ort2 is likely to be specific to sexual development.

Phylogenetic analyses indicated that the Pumilio family displays a remarkable internal structure, with clearly discernible clades, which may be interpreted as subfamilies. Each of these clades showed a high level of sequence conservation, and orthology relationships were conserved with a single copy from most species. This is in contrast with the much more dynamic evolution of the family at the scale of all eukaryotes (Najdrová et al., 2020), suggesting that within fungi, there is a higher level of Pumilio conservation. On the other hand, our data indicate that the conservation of function might not follow the conservation of sequence and orthology. Despite being a close ortholog of S. cerevisiae PUF3, our results do not support ort2 having a similar, mitochondrion-related function, which may be due to fast evolutionary turnover of function, despite a high level of sequence conservation, as has been described for RNA-regulatory networks (Wilinski et al., 2017).

The deletion of ort2 resulted in deficient dark stipe development, while its overexpression significantly increased their number. At this point, we lack a mechanistic explanation for the described qualitative and quantitative phenotypes. Dark stipe formation is a complex developmental program that remains largely unexplored. Hitherto, molecular insights into dark stipes have been limited exclusively to photoreception (Subba, 2022; Chaisaena, 2009). Given that our ort2 mutants are not ‘blind’, and dark stipe formation is more profoundly affected than fruiting body formation in light, we suggest that it acts independently of light-regulated pathways and that further investigation of ort2 might illuminate other aspects of dark stipe formation. Since dark stipes arise as a result of the growth of the basal plectenchyma (nodulus) (Terashima et al., 2005), it is tempting to hypothesize that ort2 contributes to nodulus functioning; nonetheless, further investigation is needed. We acknowledge the lack of a supported mechanistic explanation behind the described phenotypes and propose conducting transcriptomic analyses of ort2 mutants to gain more knowledge about the regulatory function of Ort2.

Our findings might have various practical applications. The overexpression of ort2 resulted in the higher abundance of dark stipes. There are popular edible species whose etiolated fruiting bodies are consumed – one such species, F. velutipes is among those with the greatest economic importance (Rugolo et al., 2016). Given our results, overexpression of the F. velutipes ort2 ortholog may be worth exploring from a commercial perspective. If the function is conserved, these findings may be applicable in mushroom cultivation. Furthermore, in our experiments fruiting bodies of the ort2 overexpression strain more frequently reached maturity than the parent strain. The unpredictability and varying success of fructification of different fungal species, even in a controlled environment (Prayook et al., 2013), pose a major challenge for the developmental biology of fruiting bodies. In our experience, even fruiting bodies of the C. cinerea AmutBmut strain frequently arrest their development in a laboratory environment. We suggest that for future fruiting body research, especially for studies focusing on late stages or the entire developmental process, the ort2OE strain might be a more suitable model system than the AmutBmut parent strain.

Ort2 is orthologous with Pum1 of Cr. neoformans, and, though RNA regulatory networks rewire quickly during evolution (Wilinski et al., 2017), a functional comparison between the two genes is of interest. Pum1 was postulated to be a pleiotropic regulator that positively influences basidium maturation and is involved in maintaining filamentous growth in Cr. neoformans (Liu et al., 2018; Kaur and Panepinto, 2016). It is notable that, at the highest level, both ort2 and PUM1 influence sexual development. The vastly different sexual development of these two species - Cr. neoformans produces simple basidial lawns (Kwon-Chung, 1976), whereas C. cinerea basidia are enclosed in complex fruiting bodies - currently precludes a more detailed comparison, however, it is possible that these genes regulate similar pathways or mechanisms. Further studies and comparisons of target gene sets in C. cinerea with those of Cr. neoformans are required to address these questions.

Ort2 is the first RNA-binding protein with a documented role in fruiting body development. Given that Pumilio proteins have been shown to possess developmental significance in animals and plants (Nishanth and Simon, 2020), further characterization of RNA regulatory networks in fungi may illuminate novel aspects of complex multicellularity in fungi. We anticipate that the characterization of ort2, reported here, will be useful in charting this new realm of fungal developmental regulation. Experimental identification of Ort2 targets might provide mechanistic explanations and an interesting future research direction.

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 author.

Author contributions

BS: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Software, Validation, Visualization, Writing – original draft, Writing – review & editing. BB: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing – review & editing. BH: Data curation, Formal Analysis, Investigation, Methodology, Software, Visualization, Writing – review & editing. MV: Investigation, Methodology, Writing – review & editing. ZH: Investigation, Methodology, Software, Writing – review & editing. X-BL: Investigation, Methodology, Writing – review & editing. HW: Investigation, Methodology, Software, Writing – review & editing. CF: Investigation, Methodology, Writing – review & editing. JG: Methodology, Software, Writing – review & editing. JP: Formal Analysis, Methodology, Software, Visualization, Writing – review & editing. ÁC: Conceptualization, Data curation, Investigation, Methodology, Project administration, Software, Supervision, Validation, Writing – original draft, Writing – review & editing. LN: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. They acknowledge support by the National Research Development and Innovation Office (Grant No. OTKA 142188), the European Research Council (grant no. 101086900 to L.G.N.), and the Hungarian Research Network (KSZF-65/2023, to LN.). This research work was conducted with the support of the National Academy of Scientist Education Program of the National Biomedical Foundation under the sponsorship of the Hungarian Ministry of Culture and Innovation. Supported by the professional funding of the University Researcher Scholarship Program of the Hungarian Ministry of Culture and Innovation, financed by the National Research, Development and Innovation Fund.

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.

The reviewer MR declared a shared affiliation with the author BB to the handling editor at the time of review.

Generative AI statement

The author(s) declare that Generative AI was used in the creation of this manuscript. AI was used solely to detect grammatical errors and to enhance the clarity of the text.

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/ffunb.2025.1633301/full#supplementary-material

Supplementary Figure S1 | Rooted gene tree of fungal Pumilio proteins, with bootstrap values. The phylogenetic analysis comprises protein sequences from 17 fungal species. All terminal nodes are labeled. The five clades are colored differently.

Supplementary Figure S2 | Plasmids used for genetic modification of C. cinerea. a. We used plasmid p.Δort2 to generate disruptants. Abbreviations: CcPAB1 cassette - selection marker amplified from the pMA412 vector (Stanley et al., 2014), UHA and DHA - upstream and downstream homologous arms, and pUC19 backbone. b. We used plasmid p.ort2OE to generate overexpression mutants. Parts of the construct: CcPAB1 cassette, hsp3 (CopciAB_423239) promoter, ort2 without its own promoter, and pUC19 backbone.

Supplementary Figure S3 | Low-frequency developmental abnormalities of the ort2 overexpression mutants. a. Branched fruiting bodies. Top row: fruiting bodies formed under alternating light/dark conditions; bottom row: etiolated fruiting bodies. Two photographs on the left represent the wild-type strain, the other four are ort2OE fruiting bodies. b-c. Prostrate etiolated fruiting bodies formed by the ort2OE strains. Blue arrows follow a dark stipe through the edge of the dish. d. ort2OE dark stipes usually intertwine and form amorphous structures. Note that at the base of the structure in the picture, individual dark stipes could not be distinguished.

Supplementary Figure S4 | UpSet plot for the overlap analysis of five C. cinerea AmutBmut gene sets: genes with mitochondrial GO annotations (mt_GO), proteins with mitochondrial localization predicted by DeepLoc 2.0 (mt_DeepLoc), proteins with mitochondrial localization predicted by TergetP 2.0 (mt_TargetP), Ort2 targets predicted by Karollus et al.’s language model (targets_LM), and Ort2 targets predicted using FIMO (targets_FIMO).

Supplementary Table S1 | Sequences of crRNAs used for the disruption through Cas9 cleavage and primers used for plasmid construction, transformant screening, and qPCR.

Supplementary Table S2 | Raw qPCR data for the examination of ort2 expression in Δort2 disruptants.

Supplementary Table S3 | Raw qPCR data and fold-changes in ort2 expression in ort2OE strains relative to the parent strain.

Supplementary Table S4 | Gene Ontology terms significantly (p < 0.05) enriched among putative Ort2 targets belonging to the most significant recognition motif predicted by the species-aware language model.

Supplementary Table S5 | Raw data and calculation of relative mtDNA abundance in ort2 mutants.

Supplementary Table S6 | The five gene sets used for functional studies: genes with mitochondrial GO annotations (mt_GO), proteins with mitochondrial localization predicted by DeepLoc 2.0 (mt_DeepLoc), proteins with mitochondrial localization predicted by TergetP 2.0 (mt_TargetP), Ort2 targets predicted by Karollus et al.’s language model (targets_LM), and Ort2 targets predicted using FIMO (targets_FIMO). targets_LM and mt_GO are provided with both Okayama and AmutBmut IDs, too. In both cases, the datasets containing Okayama identifiers are the original, and they were converted to AmutBmut, to align with the other three gene sets, resulting in the loss of some genes. For the comparison between targets_LM and mt_GO, Fisher’s exact test was performed using the lists containing Okayama IDs; for comparisons across all five datasets, the C. cinerea AmutBmut lists were used.

References

Almagro A., Juan J., Salvatore M., Emanuelsson O., Winther O., Heijne G. v., et al. (2019). Detecting sequence signals in targeting peptides using deep learning. Life Sci. Alliance 2, e2019004295. doi: 10.26508/lsa.201900429, PMID: 31570514

PubMed Abstract | Crossref Full Text | Google Scholar

Beijen E. P. W., Maanen M. H.v., Bergh E. S.v. d., Brouns R., Marian I. M., Vries T. J.de, et al. (2025). Transcription factors fst8, ftr3 and gat1 are regulators of the immune system of the mushroom schizophyllum commune. Fungal Genet. Biol. 179, 103987. doi: 10.1016/j.fgb.2025.103987, PMID: 40311743

PubMed Abstract | Crossref Full Text | Google Scholar

Bhondeley M. and Liu Z. (2020). Mitochondrial biogenesis is positively regulated by casein kinase I hrr25 through phosphorylation of puf3 in. Saccharomyces Cerevisiae. Genet. 215, 463–825. doi: 10.1534/genetics.120.303191, PMID: 32317286

PubMed Abstract | Crossref Full Text | Google Scholar

Buller A. H. R. (1958). Researches on fungi. Vol (VI. New York: Hafner Publishing Co).

Google Scholar

Chaisaena W. (2009). Light effects on fruiting body development of wildtype in comparison to light-insensitive mutant strains of the basidiomycete coprinopsis cinerea, grazing of mites (Tyrophagus putrescentiae) on the strains and production of volatile organic compounds during fruiting body development Doctoral thesis. (Göttingen, Germany: Georg-August-University Göttingen). doi: 10.53846/goediss-2341

Crossref Full Text | Google Scholar

Clémençon H. (2012). Cytology and Plectology of the Hymenomycetes. 2nd revised edition. J. Cramer. 343–346.

Google Scholar

Dörnte B. and Kües U. (2012). Reliability in transformation of the basidiomycete coprinopsis cinerea. Curr. Trends Biotechnol. Pharm. 6, 340–555.

Google Scholar

Földi C., Merényi Z., Balázs B., Csernetics Á., Miklovics N., Wu H., et al. (2024). Snowball: A novel gene family required for developmental patterning of fruiting bodies of mushroom-forming fungi (Agaricomycetes). MSystems 9. doi: 10.1128/msystems.01208-23, PMID: 38334416

PubMed Abstract | Crossref Full Text | Google Scholar

Galeota-Sprung B., Fernandez A., and Sniegowski P. (2022). Changes to the mtDNA copy number during yeast culture growth. R. Soc. Open Sci. 9. doi: 10.1098/rsos.211842, PMID: 35814911

PubMed Abstract | Crossref Full Text | Google Scholar

García-Rodríguez L. J., Gay A. C., and Pon L. A. (2007). Puf3p, a pumilio family RNA binding protein, localizes to mitochondria and regulates mitochondrial biogenesis and motility in budding yeast. J. Cell Biol. 176, 197–2075. doi: 10.1083/jcb.200606054, PMID: 17210948

PubMed Abstract | Crossref Full Text | Google Scholar

Gerber AndréP., Herschlag D., and Brown P. O. (2004). Extensive association of functionally and cytotopically related MRNAs with puf family RNA-binding proteins in yeast. PloS Biol. 2, e795. doi: 10.1371/journal.pbio.0020079, PMID: 15024427

PubMed Abstract | Crossref Full Text | Google Scholar

Grant C. E., Bailey T. L., and Noble W. S. (2011). FIMO: scanning for occurrences of a given motif. Bioinformatics 27, 1017–1185. doi: 10.1093/bioinformatics/btr064, PMID: 21330290

PubMed Abstract | Crossref Full Text | Google Scholar

Hegedüs B., Sahu N., Bálint B., Haridas S., Bense V., Merényi Z., et al. (2025). Morphogenesis, starvation, and light responses in a mushroom-forming fungus revealed by long-read sequencing and extensive expression profiling. Cell Genomics. 5, 100853. doi: 10.1016/j.xgen.2025.100853, PMID: 40262612

PubMed Abstract | Crossref Full Text | Google Scholar

Hibbett D. S. (2007). After the gold rush, or before the flood? Evolutionary morphology of mushroom-forming fungi (Agaricomycetes) in the early 21st century. Mycological Res. 111, 1001–1018. doi: 10.1016/j.mycres.2007.01.012, PMID: 17964768

PubMed Abstract | Crossref Full Text | Google Scholar

Horton P., Park K.-J., Obayashi T., Fujita N., Harada H., Adams-Collier C. J., et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, W585–W587. doi: 10.1093/nar/gkm259, PMID: 17517783

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang H., Guan W., and Gu Z. (2010). Tinkering evolution of post-transcriptional RNA regulons: puf3p in fungi as an example. PloS Genet. 6, e10010305. doi: 10.1371/journal.pgen.1001030, PMID: 20661438

PubMed Abstract | Crossref Full Text | Google Scholar

Karollus A., Hingerl J., Gankin D., Grosshauser M., Klemon K., and Gagneur J. (2024). Species-aware DNA language models capture regulatory elements and their evolution. Genome Biol. 25, 835. doi: 10.1186/s13059-024-03221-x, PMID: 38566111

PubMed Abstract | Crossref Full Text | Google Scholar

Katoh K. and Standley D. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780. doi: 10.1093/molbev/mst010, PMID: 23329690

PubMed Abstract | Crossref Full Text | Google Scholar

Kaur J. N. and Panepinto J. C. (2016). Morphotype-specific effector functions of cryptococcus neoformans PUM1. Sci. Rep. 6, 236385. doi: 10.1038/srep23638, PMID: 27008977

PubMed Abstract | Crossref Full Text | Google Scholar

Knoll A. H. (2011). The multiple origins of complex multicellularity. Annu. Rev. Earth Planetary Sci. 39, 217–239. doi: 10.1146/annurev.earth.031208.100209

Crossref Full Text | Google Scholar

Krizsán K., Almási Éva, Merényi Z., Sahu N., Virágh Máté, Kószó Tamás, et al. (2019). Transcriptomic atlas of mushroom development reveals conserved genes behind complex multicellularity in fungi. Proc. Natl. Acad. Sci. United States America 116, 7409–7418. doi: 10.1073/pnas.1817822116, PMID: 30902897

PubMed Abstract | Crossref Full Text | Google Scholar

Kües, U (2000). Life history and developmental processes in the basidiomycete coprinus cinereus. Microbiol. Mol. Biol. Rev. 64, 316–353. doi: 10.1128/MMBR.64.2.316-353.2000, PMID: 10839819

PubMed Abstract | Crossref Full Text | Google Scholar

Kües U. and Navarro-González M. (2015). How do agaricomycetes shape their fruiting bodies? 1. How do agaricomycetes shape their fruiting bodies? 1. Morphological aspects of development. Fungal Biol. Rev. 29 (2), 63–97. doi: 10.1016/j.fbr.2015.05.001

Crossref Full Text | Google Scholar

Kwon-Chung K. J. (1976). Morphogenesis of filobasidiella neoformans, the sexual state of cryptococcus neoformans. Mycologia 68, 821–833. doi: 10.1080/00275514.1976.12019959, PMID: 790172

PubMed Abstract | Crossref Full Text | Google Scholar

Liu L., He G.-J., Chen L., Zheng J., Chen Y., Shen L., et al. (2018). Genetic basis for coordination of meiosis and sexual structure maturation in cryptococcus neoformans. ELife 7. doi: 10.7554/eLife.38683, PMID: 30281018

PubMed Abstract | Crossref Full Text | Google Scholar

Livak K. J. and Schmittgen T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods 25, 402–485. doi: 10.1006/meth.2001.1262, PMID: 11846609

PubMed Abstract | Crossref Full Text | Google Scholar

Löytynoja A. (2014). Phylogeny-aware alignment with PRANK. Methods Mol. Biol. 155–170. doi: 10.1007/978-1-62703-646-7_10, PMID: 24170401

PubMed Abstract | Crossref Full Text | Google Scholar

Masuda R., Iguchi N., Tukuta K., Nagoshi T., Kemuriyama K., and Muraguchi H. (2016). The coprinopsis cinereaTup1 homologue cag1 is required for gill formation during fruiting body morphogenesis. Biol. Open 5, 1844–1525. doi: 10.1242/bio.021246, PMID: 27815245

PubMed Abstract | Crossref Full Text | Google Scholar

Matsushima Y., Sugiura R., and Kuno T. (2007). Studies on fungal pumilio gene family through mining multiple genome-scale data sets. Kobe J. Med. Sci. 53, 163–175. doi: 10.24546/81000115, PMID: 17932456

PubMed Abstract | Crossref Full Text | Google Scholar

Merényi Z., Prasanna A. N., Wang Z., Kovács Károly, Hegedüs B., Bálint Balázs, et al. (2020). Unmatched level of molecular convergence among deeply divergent complex multicellular fungi. Mol. Biol. Evol. 37, 2228–2405. doi: 10.1093/molbev/msaa077, PMID: 32191325

PubMed Abstract | Crossref Full Text | Google Scholar

Miyazaki Y., Sakuragi Y., Yamazaki T., and Shishido K. (2004). Target genes of the developmental regulator PRIB of the mushroom. Lentinula Edodes. Bioscience Biotechnology Biochem. 68, 1898–19055. doi: 10.1271/bbb.68.1898, PMID: 15388965

PubMed Abstract | Crossref Full Text | Google Scholar

Muraguchi H., Fujita T., Kishibe Y., Konno K., Ueda N., Nakahori K., et al. (2008). The exp1 gene essential for pileus expansion and autolysis of the inky cap mushroom coprinopsis cinerea (Coprinus cinereus) encodes an HMG protein. Fungal Genet. Biol. 45, 890–965. doi: 10.1016/j.fgb.2007.11.004, PMID: 18164224

PubMed Abstract | Crossref Full Text | Google Scholar

Muraguchi H. and Kamada T. (1998). The ich1 gene of the mushroom coprinus cinereus is essential for pileus formation in fruiting. Development 125, 3133–3415. doi: 10.1242/dev.125.16.3133, PMID: 9671586

PubMed Abstract | Crossref Full Text | Google Scholar

Muraguchi H., Umezawa K., Niikura M., Yoshida M., Kozaki T., Ishii K., et al. (2015). Strand-specific RNA-seq analyses of fruiting body development in coprinopsis cinerea. PloS One 10, e0141586. doi: 10.1371/journal.pone.0141586, PMID: 26510163

PubMed Abstract | Crossref Full Text | Google Scholar

Nagy László G., Kovács Gábor M, and Krizsán K. (2018). Complex multicellularity in fungi: evolutionary convergence, single origin, or both? Biol. Rev. 93, 1778–1945. doi: 10.1111/brv.12418, PMID: 29675836

PubMed Abstract | Crossref Full Text | Google Scholar

Nagy L. G., Vonk P. J., Künzler M., Földi C., Virágh M., Ohm R. A., et al. (2023). Lessons on fruiting body morphogenesis from genomes and transcriptomes of agaricomycetes. Stud. Mycology. 104, 1–85. doi: 10.3114/sim.2022.104.01, PMID: 37351542

PubMed Abstract | Crossref Full Text | Google Scholar

Najdrová Vladimíra, Stairs C. W., Vinopalová M., Voleman Luboš, and Doležal P. (2020). The evolution of the puf superfamily of proteins across the tree of eukaryotes. BMC Biol. 18, 775. doi: 10.1186/s12915-020-00814-3, PMID: 32605621

PubMed Abstract | Crossref Full Text | Google Scholar

Naranjo-Ortiz M. A. and Gabaldón T. (2019). Fungal evolution: diversity, taxonomy and phylogeny of the fungi. Biol. Rev. 94, 2101–2375. doi: 10.1111/brv.12550, PMID: 31659870

PubMed Abstract | Crossref Full Text | Google Scholar

Nishanth M. J. and Simon B. (2020). Functions, mechanisms and regulation of pumilio/Puf family RNA binding proteins: A comprehensive review. Mol. Biol. Rep. 47, 785–8075. doi: 10.1007/s11033-019-05142-6, PMID: 31643042

PubMed Abstract | Crossref Full Text | Google Scholar

Ohm R. A., de Jong J. F., Bekker C. de, Wösten H. A.B., and Lugones L. G. (2011). Transcription factor genes of schizophyllum commune involved in regulation of mushroom formation. Mol. Microbiol. 81, 1433–1455. doi: 10.1111/j.1365-2958.2011.07776.x, PMID: 21815946

PubMed Abstract | Crossref Full Text | Google Scholar

Pareek M., Hegedüs B., Hou Z., Csernetics Árpád, Wu H., Virágh Máté, et al. (2022). Preassembled cas9 ribonucleoprotein-mediated gene deletion identifies the carbon catabolite repressor and its target genes in coprinopsis cinerea. Appl. Environ. Microbiol. 88. doi: 10.1128/aem.00940-22, PMID: 36374019

PubMed Abstract | Crossref Full Text | Google Scholar

Pelkmans J. F., Patil M. B., Gehrmann T., Reinders M. J.T., Wösten H. A.B., and Lugones L. G. (2017). Transcription factors of schizophyllum commune involved in mushroom formation and modulation of vegetative growth. Sci. Rep. 7, 3105. doi: 10.1038/s41598-017-00483-3, PMID: 28331193

PubMed Abstract | Crossref Full Text | Google Scholar

Pelkmans J. F., Vos A. M., Scholtmeijer K., Hendrix Ed, Baars J. J.P., Gehrmann T., et al. (2016). The transcriptional regulator C2h2 accelerates mushroom formation in agaricus bisporus. Appl. Microbiol. Biotechnol. 100, 7151–7595. doi: 10.1007/s00253-016-7574-9, PMID: 27207144

PubMed Abstract | Crossref Full Text | Google Scholar

Prayook S., Ananyaphon P., and Preecha P. (2013). Initiation of fruiting body development of the medicinal mushroom phellinus linteus from Cambodia. Afr. J. Microbiol. Res. 7, 2885–2925. doi: 10.5897/AJMR12.1464

Crossref Full Text | Google Scholar

Rabe B. A. and Cepko C. (2020). A simple enhancement for gibson isothermal assembly. bioRxiv. doi: 10.1101/2020.06.14.150979

Crossref Full Text | Google Scholar

Rugolo M., Levin L., and Lechner B. E. (2016). Flammulina velutipes: an option for ‘Alperujo’ Use. Rev. Iberoamericana Micología 33, 242–475. doi: 10.1016/j.riam.2015.12.001, PMID: 27468641

PubMed Abstract | Crossref Full Text | Google Scholar

Sakamoto Y., Sato S., Yoshida H., Takahashi M., Osakabe K., and Muraguchi H. (2024). The exp2 gene, which encodes a protein with two zinc finger domains, regulates cap expansion and autolysis in coprinopsis cinerea. Microbiological Res. 283, 127695. doi: 10.1016/j.micres.2024.127695, PMID: 38554651

PubMed Abstract | Crossref Full Text | Google Scholar

Stajich J. E., Wilke S. K., Ahrén D., Au C. H., Birren B. W., Borodovsky M., et al. (2010). Insights into evolution of multicellular fungi from the assembled chromosomes of the mushroom coprinopsis cinerea (Coprinus cinereus). Proc. Natl. Acad. Sci. 107, 11889–11894. doi: 10.1073/pnas.1003391107, PMID: 20547848

PubMed Abstract | Crossref Full Text | Google Scholar

Stamatakis A. (2014). RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033, PMID: 24451623

PubMed Abstract | Crossref Full Text | Google Scholar

Stanley C. E., Stöckli M., Swaay D. v., Sabotič J., Kallio P. T., Künzler M., et al. (2014). Probing bacterial–fungal interactions at the single cell level. Integr. Biol. 6, 935–455. doi: 10.1039/C4IB00154K, PMID: 25144657

PubMed Abstract | Crossref Full Text | Google Scholar

Subba S. (2022). Fruiting body development of coprinopsis cinerea: cytology and regulatory factors. Georg-August-University Göttingen. doi: 10.53846/goediss-9177

Crossref Full Text | Google Scholar

Terashima K., Yuki K., Muraguchi H., Akiyama M., and Kamada T. (2005). The dst1 gene involved in mushroom photomorphogenesis of coprinus cinereus encodes a putative photoreceptor for blue light. Genetics 171, 101–185. doi: 10.1534/genetics.104.040048, PMID: 15956671

PubMed Abstract | Crossref Full Text | Google Scholar

Thumuluri V., Armenteros JoséJ. A., Johansen A. R., Nielsen H., and Winther O. (2022). DeepLoc 2.0: multi-label subcellular localization prediction using protein language models. Nucleic Acids Res. 50, W228–W234. doi: 10.1093/nar/gkac278, PMID: 35489069

PubMed Abstract | Crossref Full Text | Google Scholar

Tsusué Y. M. (1969). Experimental control of fruit-body formation in coprinus macrorhizus. Development Growth Differentiation 11, 164–178. doi: 10.1111/j.1440-169X.1969.00164.x, PMID: 5370065

PubMed Abstract | Crossref Full Text | Google Scholar

Venegas V., Wang J., Dimmock D., and Wong L.-J. (2011). Real-time quantitative PCR analysis of mitochondrial DNA content. Curr. Protoc. Hum. Genet. 68. doi: 10.1002/0471142905.hg1907s68, PMID: 21234878

PubMed Abstract | Crossref Full Text | Google Scholar

Virágh Máté, Merényi Z., Csernetics Árpád, Földi C., Sahu N., Liu X.-B., et al. (2022). Evolutionary morphogenesis of sexual fruiting bodies in basidiomycota: toward a new evo-devo synthesis. Microbiol. Mol. Biol. Rev. 86. doi: 10.1128/MMBR.00019-21, PMID: 34817241

PubMed Abstract | Crossref Full Text | Google Scholar

Vonk P. J., van der Poel M. J. P., Niemeijer ZoéE., and Ohm R. A. (2024). The regulatory network of the white collar complex during early mushroom development in schizophyllum commune. Microbiological Res. 284, 127736. doi: 10.1016/j.micres.2024.127736, PMID: 38663231

PubMed Abstract | Crossref Full Text | Google Scholar

Wang Z., Sun X., Wee J., Guo X., and Gu Z. (2019). Novel Insights into Global Translational Regulation through Pumilio Family RNA-Binding Protein Puf3p Revealed by Ribosomal Profiling. Curr. Genet. 65, 201–125. doi: 10.1007/s00294-018-0862-4, PMID: 29951697

PubMed Abstract | Crossref Full Text | Google Scholar

Webster M. W., Stowell J. A. W., and Passmore L. A. (2019). RNA-binding proteins distinguish between similar sequence motifs to promote targeted deadenylation by ccr4-not. ELife 8. doi: 10.7554/eLife.40670, PMID: 30601114

PubMed Abstract | Crossref Full Text | Google Scholar

Wilinski D., Buter N., Klocko A. D., Lapointe C. P., Selker E. U., Gasch A. P., et al. (2017). Recurrent rewiring and emergence of RNA regulatory networks. Proc. Natl. Acad. Sci. 114. doi: 10.1073/pnas.1617777114, PMID: 28320951

PubMed Abstract | Crossref Full Text | Google Scholar

Xie Y., Chang J., and Kwan H. S. (2020). Carbon metabolism and transcriptome in developmental paths differentiation of a homokaryotic coprinopsis cinerea strain. Fungal Genet. Biol. 143, 103432. doi: 10.1016/j.fgb.2020.103432, PMID: 32681999

PubMed Abstract | Crossref Full Text | Google Scholar

Yan J.-J., Tong Z.-J., Liu Y.-Y., Li Y.-N., Zhao C., Mukhtar I., et al. (2019). Comparative transcriptomics of flammulina filiformis suggests a high CO2 concentration inhibits early pileus expansion by decreasing cell division control pathways. Int. J. Mol. Sci. 20, 59235. doi: 10.3390/ijms20235923, PMID: 31775357

PubMed Abstract | Crossref Full Text | Google Scholar

Zadražil F. (1975). Influence of CO2 concentration on the mycelium growth of three pleurotus species. Eur. J. Appl. Microbiol. 1, 327–335. doi: 10.1007/BF01382692

Crossref Full Text | Google Scholar

Keywords: fruiting body formation, dark stipe, RNA-binding proteins, Pumilio, language model

Citation: Szathmári B, Bálint B, Hegedüs B, Virágh M, Hou Z, Liu X-B, Wu H, Földi C, Gagneur J, Promeuschel J, Csernetics Á and Nagy LG (2025) Orthrus: a Pumilio-family gene involved in fruiting body and dark stipe development in Coprinopsis cinerea. Front. Fungal Biol. 6:1633301. doi: 10.3389/ffunb.2025.1633301

Received: 22 May 2025; Accepted: 11 July 2025;
Published: 30 July 2025.

Edited by:

Isheng Jason Tsai, Academia Sinica, Taiwan

Reviewed by:

György Sipos, University of Sopron, Hungary
Hajime Muraguchi, Akita Prefectural University, Japan
Marjatta Raudaskoski, University of Turku, Finland

Copyright © 2025 Szathmári, Bálint, Hegedüs, Virágh, Hou, Liu, Wu, Földi, Gagneur, Promeuschel, Csernetics and Nagy. 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: László G. Nagy, bG5hZ3lAZnVuZ2Vub21lbGFiLmNvbQ==

Present address: Balázs Bálint, Medical Bioinformatics Centre, University of Turku, Turku, Finland

These authors have contributed equally

Disclaimer: 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.