1. Introduction

In eukaryotes, gene transcription is regulated by dynamic changes in chromatin. The posttranslational modifications of core histones, including acetylation, methylation, and ubiquitination, represent major mechanisms by which cells alter the chromatin structural properties and regulate gene transcription [1, 2]. Among them, the monoubiquitination of a lysine (K) residue on the C-terminal of histone H2B (H2Bub1) is a conserved modification that occurs on H2B K120 residue in Homo sapiens and K123 residue in Saccharomyces cerevisiae [3, 4]. H2Bub1 is enriched at regions of active transcription but plays roles in both gene activation and repression [57]. In addition, H2Bub1 is required for di- and trimethylation of H3 K4 and H3 K79, subsequently modulating chromatin accessibility [812].

The ubiquitin conjugase (E2) Rad6 and the ubiquitin ligase (E3) Bre1 are responsible for H2Bub1 in S. cerevisiae [1315]. In addition, H2Bub1 is regulated by additional factors in yeast and other eukaryotes, among which the conserved polymerase-associated factor 1 (Paf1) complex (Paf1C) is the prominent one [4, 1619] Paf1C is a multi-functional protein complex, which impacts RNA synthesis at multiple stages [2028]. Paf1C consists of the subunits Paf1, Ctr9, Cdc73, Rtf1, and Leo1, and the five subunits are stably associated within the complex in S. cerevisiae [22, 2931]. In contrast, Rtf1 appears not to be stably associated with Paf1C in human cells, despite the Paf1C is structurally and functionally conserved [3234]. Interestingly, the histone modification domain (HMD) within Rtf1 is both necessary and sufficient for stimulating H2Bub1 in yeast [4, 18]. Expression of the Rtf1 HMD alone restores H2Bub1 levels in S. cerevisiae mutants deleted for the RTF1 gene or all five Paf1C subunits-encoding genes [3, 4, 18]. These studies show that Rtf1 is the only Paf1C subunit that is strictly required for deposition of H2Bub1 in vivo. However, little is known about its role in human fungal pathogens.

Cryptococcus neoformans, the top-ranked fungus in the WHO fungal pathogen priority list, is a globally distributed opportunistic fungal pathogen that can cause life-threatening cryptococcosis [35, 36]. The mortality rate of cryptococcosis is alarmingly high, especially in patients with HIV infection, in whom it ranges from 41% to 61% [37, 38]. C. neoformans can be classified into two serotypes: the serotype A C. neoformans and the serotype D C. deneoformans. Both C. neoformans and C. deneoformans undergo yeast-to-hypha transition under inducing conditions, which has been shown to be associated with fungal virulence [39, 40]. Thus, deciphering the regulatory mechanisms on fungal morphogenesis and pathogenesis in C. neoformans is critical for comprehensive understanding of the nature of pathogen and combating against cryptococcal infection.

In our previous study, we characterized the subunits of complex-associated with Set1 (COMPASS) and found that COMPASS-mediated H3K4 methylation (H3K4me) affects yeast- to-hypha transition and virulence in both C. neoformans and C. deneoformans [41]. We also preliminarily showed that H2Bub1 is required for COMPASS-mediated H3K4me by deletion of RAD6 and RTF1 in C. neoformans and C. deneoformans [41]. However, the Here, we set out to characterize the roles of Rtf1 in facilitating global H2Bub1 and to gain comprehensive insights into the epigenetic regulation on fungal morphogenesis and pathogenesis in the human fungal pathogen C. neoformans.

2. Materials and methods

2.1. Strains, culture conditions, and microscopy examination

Strains used in this study are listed in the Supplemental Table S1. C. deneoformans and C. neoformans strains were maintained on YPD medium unless specified otherwise. Transformants obtained from transient CRISPR-Cas9 coupled with electroporation (TRACE) were selected on YPD medium with 100 μg/mL of nourseothricin, 100 μg/mL of G418, or 200 μg/mL of hygromycin.

Strains for phenotypic assays were grown overnight in liquid YPD medium at 30°C with shaking. The cells were washed with sterile water, adjusted to an optical density at 600 nm (OD600) of 3.0, and serially diluted. For filamentation tests, aliquots (3 μL) of cell suspensions (OD600 = 3.0) were spotted onto V8 plates and cultured at room temperature in the dark. For morphological examinations, all strains were examined under a stereoscope. For spotting assays, aliquots (3 μL) of serial dilutions starting from OD600 = 3.0 were spotted onto agar medium with supplements and cultured under the noted conditions.

2.2. Gene manipulation

Cryptococcal genes were deleted following the TRACE protocol [42, 43]. In brief, a deletion construct with approximately 1 kb of homologous arms flanking a target gene and the dominant marker was cloned through fusion PCR. This construct was mixed with PCR products of CAS9 and a relevant guide RNA (gRNA), and the mixture was introduced into recipient strains by electroporation as described previously [43]. Resulting yeast colonies were screened by two rounds of diagnostic PCR. The first round of PCR was performed to detect the integration of the construct into the corresponding locus of the target gene. The second round of PCR was performed to confirm knockout of the target fragment. All primers used to make gene deletion mutants are listed in the Supplemental Table S1.

For gene complementation, the ORFs plus approximately 1.0 kb of their upstream regions were amplified by PCR and cloned into vectors through T5 exonuclease-dependent assembly as previously described [44]. For gene overexpression with inducible or constitutively active promoters, the constructs were obtained by amplifying the entire ORF by PCR and cloning the PCR products into vectors at the downstream of CTR4, TEF1, or GPD1 promoter. All plasmids were confirmed by restriction enzyme digestion and sequencing. The confirmed constructs, together with PCR products of CAS9 and gRNA targeting the Safe Haven locus [45, 46], were introduced into recipient Cryptococcus strains. The transformants were passaged once per day for five days and cultured on selection plates to obtain stable transformants. Then, two rounds of diagnostic PCR were performed to confirm the integration and orientation of constructs into the Safe Haven locus. All primers and plasmids used for gene complementation and overexpression are listed in the Supplemental Table S1.

2.3. Protein extraction and western blotting

Proteins were extracted from Cryptococcus cells according to a previously described method [47]. Aliquots of proteins were separated on 4%-to-12% gradient SDS-PAGE gels and then transferred to a polyvinylidene difluoride membrane for Western blot analyses. Antibodies used in this study are listed in the Supplemental Table S1. For co-immunoprecipitation assays coupled with mass spectrometry (CoIP/MS), whole cell extracts of experimental strains were incubated with FLAG-trap (Sigma) according to the manufacturer’ s instructions. Proteins in the eluted samples were loaded in SDS-PAGE gel, digested, and analysed by the proteome facility centre of Institute of Microbiology, Chinese Academy of Sciences.

2.4. RNA extraction and qPCR assays

Cryptococcus strains were cultured in liquid YPD with shaking at 220 rpm at 30 °C overnight, or on solid V8 medium at room temperature in the dark for 24 h. The cultures were collected, flash frozen in liquid nitrogen, and lyophilized for 24 h. Total RNA was isolated with the PureLink RNA Mini Kit (Invitrogen), and first strand cDNA was synthesized using the GoScript Reverse Transcription System (Promega) following the manufacturer’s instructions. The Power SYBR Green system (Invitrogen) was used for RT-PCR. All the primers used here are listed in the Supplemental Table S1. Relative transcript levels were determined using the ΔΔCt method as described previously. Three biological replicates were included for all tests. Statistical significance was determined using a Student’s t-test. Differences for which p < 0.05 were considered statistically significant.

2.5. RNA-seq and data analysis

For RNA-seq analyses, strains were cultured in YPD liquid medium at 30 °C overnight. The cells were washed with ddH2O and spotted on V8 medium to stimulate unisexual reproduction. The level and integrity of RNA in each sample were evaluated using a Qubit RNA Assay Kit on a Qubit 2.0 Fluorometer (Life Technologies, CA, USA) and RNA Nano 6000 Assay Kit with the Bioanalyzer 2100 system (Agilent Technologies, CA, USA), respectively. RNA purity was assessed using a Nano Photometer spectrophotometer (IMPLEN, CA, USA). The transcriptome libraries were generated using the VAHTS mRNA-seq v2 Library Prep Kit (Vazyme Biotech Co., Ltd, Nanjing, China) according to the manufacturer’s instructions.

The transcriptome libraries were sequenced by Annoroad Gene Technology Co., Ltd (Beijing, China) on an Illumina platform. For RNA-seq analysis, the quality of sequenced clean data was analyzed using FastQC software. Subsequently, sequences from approximately 2 GB of clean data for each sample were mapped to the genome sequence of C. deneoformans XL280α using STAT. Gene expression levels were measured in transcripts per million (TPM) by Stringtie to determine unigenes. All unigenes were subsequently aligned against the well- annotated genome of JEC21, which served as the parent strain to generate XL280α through a cross with B3501α. The differential expression of genes (DEGs) was assessed using DEseq2 of the R package and defined based on the fold change criterion (log2|fold-change | > 1.0, adjusted p value < 0.05).

2.6. Virulence trait assays

Strains for examining virulence factors were grown overnight in liquid YPD at 30 °C with shaking. The overnight cultures were washed with sterile water, adjusted to OD600 = 3.0, and serially diluted. For thermal tolerance, melanin production, and capsule formation assay, aliquots (3 μL) of serially diluted cell suspensions were spotted onto YPD plate, L-dopamine media, and low-iron media [48], respectively. Thermal tolerance was test at 30, 37, and 39 °C; melanin production was tested at 30 °C in the dark; capsule formation was tested at 37 °C with 5% CO2. All assays were repeated at least three times.

2.7. Murine models of cryptococcosis

Intranasal infection model

Female Balb/C mice of 8–10 weeks old were purchased from the Laboratory Animal Center of Zhengzhou University, China. Cryptococcal strains were inoculated in 3 mL of liquid YPD medium with the initial OD600=0.2 (approximately 106 cell/mL) and incubated for 15 hr at 30 °C with shaking. Prior to intranasal infection, cells were washed with sterile saline three times and adjusted to the final concentration of 2×106 cell/mL. Once the mice were sedated with ketamine and xylazine via intraperitoneal injection, 50 μL of the cell suspension (1×105 cells per mouse) were inoculated intranasally as previously described [41, 4951]. Mice were monitored daily for disease progression. Animals were euthanized at 10 DPI, and lungs were dissected for fungal burden quantification.

Intravenous infection model

Prior to intravenous infections, cryptococcal cells were washed with sterile saline three times and adjusted to the final concentration of 2×106 cell/mL. Mice were sedated with Isoflurane. 50 μL of the cell suspension (1×105 cells per mouse) were injected intravenously as previously described [41, 4951]. After DPI 7, animals were euthanized, and the brain, lungs, kidneys, and spleens were dissected.

For fungal burden quantifications, dissected organs were homogenized in 2 mL of cold sterile PBS. Tissue suspensions were serially diluted in PBS and plated onto YNB agar medium and incubated at 30°C for 2 days before counting the CFUs.

2.8. DAPI staining

DAPI (4’,6-diamidino-2-phenylindole) staining assays were performed as previously described [52]. Briefly, yeast cells or hyphae were collected and fixed with 3.7% formaldehyde and permeabilized in 1% Triton X-100. The cells were then washed three times with PBS and incubated in 2 μg/mL DAPI before being dropped onto a glass slide for fluorescent microscopic observation.

2.9. Data availability

All RNA-seq data are going to be available at the NCBI (SUB14425795).

3. Results

3.1. Rtf1-mediated global H2Bub1 regulates cryptococcal yeast-to-hypha transition

PAF1C subunit Rtf1 functions at the interface between Paf1C and Rad6/Bre1, and is required for deposition of H2Bub1 in all the eukaryotic species examined [53]. We showed that Rtf1 is also required for H2Bub1 and subsequent COMPASS-mediated H3K4me in the C. deneoformans reference strain XL280α background (Figure 1A and B) [41]. Interestingly, loss of H2Bub1 through deleting RTF1 blocked unisexual yeast-to-hypha transition in C. deneoformans (Supplementary Figure S1A) [41]. To establish the role of Rtf1 in regulating cryptococcal filamentation during bisexual mating, we obtained RTF1 deletion mutant in the C. deneoformans reference strain XL280a background through spore dissection from cross between rtf1Δα and XL280a, and conducted bisexual cross assay under mating-inducing condition on V8 media. The mating hyphae during unilateral mating between rtf1Δα and XL280awere produced at a slightly reduced level compared to mating between reference partners XL280α and XL280a, while filamentation was significantly reduced during bilateral mating between rtf1Δα and rtf1Δa(Figure 1C).

Rtf1 regulates cryptococcal bisexual mating by facilitating H2Bub1 via the HMD domain.

(A) Immunoblot analysis of H2Bub1 in C. deneoformans wild-type XL280, rtf1Δ, and rtf1Δ/RTF1 strains. (B) Immunoblot analysis of H3K4me (including H3K4me1, H3K4me2, and H3K4me3) in C. deneoformans wild-type XL280, rtf1Δ, and rtf1Δ/RTF1 strains. (C) Colony morphology of cells during bisexual mating between indicated strains. The same volume of cells with opposite mating type at OD600 = 3 were mixed, and 3 μL of mixtures were spotted and cultured on V8 media for 2 days at room temperature in dark. (D) Schematic diagram of pheromone-dependent signaling pathway and transcript levels of genes involved in pheromone signaling. The mating cells were prepared and cultured on V8 media following the same protocol as the colony morphology assay. After 24 h, cells were collected for total RNA extraction and qPCR. (E) Domain structure of Rtf1 homologs in indicated eukaryotes. Cn, C. deneoformans, Ca, Candida albicans, Af, Aspergillus fumigatus, Fg, Fusarium graminearum, Nc, Neurospora crassa, Sp, Schizosaccharomyces pombe, Sc, S. cerevisiae, Dm, Drosophila melanogaster, Hs, Homo sapiens. (F and G) Neighbor- joining tree of Rtf1 homologs and their corresponding HMD in indicated eukaryotes. (H) Schematic diagram of overexpressing constructs of Rtf1, HMD domain, and Plus3 domain. The constitutive promoter of TEF1 gene was used to drive gene expression. (I) Immunoblot analysis of H2Bub1 in strains expressing the indicated proteins. (J) Unisexual hyphal formation of indicated strains on V8 media.

During bisexual mating in C. deneoformans, mating pheromone (MF) is produced in cells and secreted through the transporter Ste6 [54]. Secreted pheromone induces mating response by binding to the compatible receptor on the cell surface of opposite mating type (Mfa to Ste3α or Mfα to Ste3a) [55, 56]. In addition, Mat2, which is a direct downstream transcription factor of the Cpk1 MAPK pathway, regulates the transcription of genes encoding the above-mentioned pheromone sensing proteins [57] (Figure 1D). Given the bisexual mating hyphae reduction caused by RTF1 deletion, we further investigated the effects of RTF1 deletion on genes involved in bisexual mating at transcript level via qPCR. In comparison to the mating- suppressing condition (YPD media), the transcript levels of MAT2, MFα2, STE6, and STE3α were all highly induced under mating-inducing condition (V8 media). However, these inductions were significantly impaired by deletion of RTF1 (Figure 1D). These results strongly indicated that Rtf1 facilitates H2Bub1 and regulates the expression of genes involved in fungal morphogenesis in C. deneoformans.

3.2. Ectopic expression of HMD restores global H2Bub1 levels and cryptococcal yeast-to-hypha transition

As the key subunit of Paf1C in mediating histone H2Bub1, Rtf1 is conserved across eukaryotes and consists of two conserved domains, a histone modification domain HMD and a domain that contains three highly conserved positively charged residues (Plus3) (Figure 1E). It is worth noting that Rtf1 protein and Plus3 domain in C. neoformans is evolutionally close to higher eukaryotes, such as H. sapiens and Drosophila melanogaster (Figure 1F and Supplementary Figure S1B and C), while the HMD domain is distant from higher eukaryotes (Figure 1G). To further dissect the roles of Rtf1 HMD and Plus3 in facilitating histone H2Bub1 in C. deneoformans, we constructed the truncated versions of Rtf1 that encode HMD (53-145) or Plus3 (227-333) with a nuclear localization sequence (NLS) added in their N terminus, respectively, driven by the constitutive promoter and tagged with FLAG (Figure 1H). Interestingly, overexpression of HMD domain itself significantly promoted H2Bub1 to an even higher level in the rtf1Δ strain, compared to that in WT strain and the strain overexpressing the full length of RTF1 (Figure 1I), while overexpression of the Plus3 failed to restore H2Bub1 (Figure 1I). These results demonstrated that HMD alone is sufficient to facilitate the global H2Bub1 level in C. deneoformans.

Our previous studies have demonstrated that H2Bub1 is positively related to the filamentation in C. neoformans [41]. Consistently, overexpressing either the full length of RTF1 or the HMD domain, but not Plus3, promoted the filamentation in rtf1Δ strain (Figure 1J). To gain an overview of effects on gene expression patterns by the overexpression of HMD domain, we conducted transcriptome profiling by RNA-seq under filamentation-inducing condition (on V8 media). The results showed that the expression levels of 668 genes were significantly changed due to the RTF1 deletion compared to the WT on V8 media (|log2FC| > 1, adjusted P value < 0.05), with 308 genes significantly upregulated and 360 genes downregulated (Figure 2A, Supplementary Data S1). It is worth noting that the downregulated genes are mainly enriched in GO categories related to sexual reproduction, pheromone-dependent signaling, and filamentous growth (Supplementary Figure S2). Strikingly, overexpression of HMD domain alone in rtf1Δ strain successfully restored the expression of these genes to similar levels as those in wild-type XL280 strain, while overexpression of Plus3 domain failed to do so (Figure 2A). In particular, the expression levels of marker genes of filamentous growth (ZNF2 and CFL1) [57, 58]and genes involved in sexual reproduction and pheromone-dependent signaling (MFα, STE3α, and STE6) as shown in Figure 1D were restored by overexpressing HMD domain alone in the rtf1Δ background (Figure 2B, Supplementary Figure S3). These findings from transcriptome analyses were further confirmed by qPCR (Figure 2C). In addition, the downregulated genes in rtf1Δ/Plus3 cells were significantly enriched in common GO categories as the downregulated genes in rtf1Δ cells (Supplementary Figure S2), relative to the wild-type XL280 strain. Together, these results strongly suggested that HMD domain is sufficient to facilitate global H2Bub1 to promote expression of genes associated with filamentation.

The expression of HMD domain alone rescues the downregulation of genes involved in pheromone signaling and filamentous growth due to the deletion of RTF1.

(A) Volcano plots of differentially expressed genes in rtf1Δ, rtf1Δ/HMD, and rtf1Δ/Plus3 relative to wild-type XL280, respectively, during unisexual mating on V8 media. Genes involved in pheromone signaling and filamentous growth were indicated. (B) Reads coverage of indicated gene loci in XL280, rtf1Δ, rtf1Δ/HMD, and rtf1Δ/Plus3 strains. Reads coverage at TEF1 locus served as control. (C) qPCR quantification of transcript levels of indicated genes in XL280, rtf1Δ, rtf1Δ/HMD, and rtf1Δ/Plus3 strains cultured on V8. The transcript level of indicated genes were relative to its transcript level in XL280 cells cultured in YPD.

3.3. HMD is sufficient to facilitate global H2Bub1 and the consequent yeast-to- hypha transition

The full length Rtf1 or HMD domain should properly translocate into the nucleus to facilitate histone H2Bub1. To confirm the function of Rtf1 and HMD domain in facilitating H2Bub1, we artificially intervened their sub-cellular localizations and evaluated the effects of non-nuclear (cell membrane) and nuclear localizations on H2Bub1, H3K43me, and filamentation. To achieve cell membrane localization, we fused the full length Rtf1 and HMD with a cell membrane RGS2-mNeonGreen tag [59, 60] (Figure 3A), and introduced the constructs into the rtf1Δ strain, respectively. As indicated by the mNeonGreen fluorescence, HMD domain and the originally nuclear-localized full length Rtf1 translocated to cell membrane after fusing with the RGS2-mNeongreen tag (Figure 3B). Both nuclear-localized Rtf1 and HMD domain restored the levels of H2Bub1, H3K4me, and filamentation (Figure 3C and D). In contrast, the non-nuclear-localized full length Rtf1 or HMD domain failed to restore the levels of H2Bub1, H3K4me, or filamentation in the rtf1Δ strain (Figure 3C and D). These results further supported the role of Rtf1 HMD domain in facilitating H2Bub1.

The HMD domain alone is sufficient to facilitate global H2Bub1 and restore hyphal formation in rtf1Δ strain.

(A) Schematic diagram of constructs expressing mNG-labelled Rtf1 or HMD domain with cell membrane tag (RGS2) or NLS. (B) Fluorescence analysis of sub-cellular localizations of Rtf1 and HMD with RGS2 and NLS. (C) Immunoblot analysis of H2Bub1 and H3K4me in strains expressing the indicated proteins. (D) Colony morphology of indicated strains during unisexual mating on V8. (E) ClustalW multiple amino acid sequence alignment of the HMD domain in the indicated eukaryotes. The E95 and F118 residues in C.de neoformans were indicated with asterisk and dot, respectively. Cn, C. deneoformans, Ca, Candida albicans, Af, Aspergillus fumigatus, Fg, Fusarium graminearum, Nc, Neurospora crassa, Sp, Schizosaccharomyces pombe, Sc, S. cerevisiae, Dm, Drosophila melanogaster, Hs, Homo sapiens. (F and G) The 3D structure of C. deneoformans HMD domain predicted by SWISS-MODEL with the 3D structure of S. cerevisiae HMD domain (5emx) as the template. The conserved E95 residue was indicated in red. (H) Immunoblot analysis of H2Bub1 and H3K4me in strains expressing the indicated proteins. (I) Colony morphology of indicated strains during unisexual mating on V8.

Rtf1 HMD domain is conserved from various eukaryotic species, and the residue of glutamine at position 95 (E95, Figure 3E, F and G) has been shown to be critical for the function of Rtf1 [28]. It is noteworthy that the residue of phenylalanine at position 118 (F118) in C. neoformans is as conserved as the residue in S. cerevisiae (Figure 3E), which is critical for H2Bub1 in yeast, although it is not conserved in other eukaryotic species [28]. To investigate their roles in cryptococcal Rtf1 HMD domain, we constructed site-mutated alleles of full length Rtf1 (Rtf1E95A and Rtf1F118A) and HMD domain (HMDE95A and HMDF118A) and introduced them into the rtf1Δ strain, respectively. Both Rtf1E95A and HMDE95A failed to restored H2Bub1 and H3K4me levels in the rtf1Δ strain, while Rtf1F118A and HMDF118A partially restored H2Bub1 and H3K4me levels (Figure 3H). In consistent with the histone modification outputs, the mutants expressing Rtf1E95A or HMDE95A showed non-filamentous phenotypes similar as the staring rtf1Δ strain, while mutants expressing Rtf1F118A or HMDF118A produced more filaments than the starting rtf1Δ strain (Figure 3I). Together, these results demonstrated that Rtf1 HMD domain itself is sufficient to facilitate H2Bub1 and consequent cryptococcal filamentation with E95 as a critical conserved residue.

3.4. Roles of the global H2Bub1 level in cryptococcal virulence factor production

To investigate the role of HMD-mediated H2Bub1 in cryptococcal virulence, we constructed RTF1 deletion strain and mutants overexpressing the full length Rtf1, HMD domain, Plus3 domain, or mutated alleles of Rtf1E95A and HMD E95A, respectively, in the clinically isolated serotype A C. neoformans H99 strain background. Consistent with what we observed in the serotype D C. deneoformans, deletion of RTF1 abolished H2Bub1 and H3K4me, and overexpressing the full length of Rtf1 and HMD domain alone, but not the Plus3 domain, Rtf1E95A or HMD E95A, successfully restored H2Bub1 and H3K4me (Figure 4A). Next, we investigated whether Rtf1 HMD domain is involved in the production of major virulence factors in vitro and pathogenicity in murine models of cryptococcosis. As shown in Figure 4B, the rtf1Δ mutant appeared to produce wild-type levels of capsule, which was much less than the pas3Δ control strain under capsule-inducing media, indicating that capsule production is not affected by RTF1 deletion in C. neoformans.

Rtf1 and HMD domain regulate the production of virulence factors in C. neoformans.

(A) Immunoblot analysis of H2Bub1 and H3K4me in strains expressing the indicated proteins in C. neoformans. (B) Capsule production in the indicated strains on capsule-inducing media. The capsule-overproducing strain pas3Δ was used as control [52]. (C) Thermal tolerance of indicated strains. (D) melanin production of indicated strains on L-DOPA media.

In contrast to the minor role in capsule production, the global H2Bub1 level caused by RTF1 deletion played critical roles in thermal tolerance and melanin production in C. neoformans. The rtf1Δ mutant showed severe growth defect at 39°C, and overexpression of the full length Rtf1, but not the Plus3 domain, Rtf1E95A or HMD E95A, partially restored the thermal tolerance of the rtf1Δ mutant (Figure 4C). Interestingly, overexpression of HMD domain alone restored the growth defect of rtf1Δ mutant at 39°C to a level that was worse than the expression of full length of Rtf1 (Figure 4C). Furthermore, the rtf1Δ mutant was incapable to produce melanin, and only the full length of Rtf1 restored its melanin production, while the HMD domain alone failed to do so (Figure 4D). These results strongly indicate that the full length of Rtf1, but not only the levels of global H2Bub1, is required to regulate thermal tolerance and melanin production in C. neoformans.

3.5. Roles of the HMD-mediated H2Bub1 in regulating cryptococcal pathogenicity

To further investigate the role of HMD-mediated H2Bub1 in the pathogenicity of C. neoformans, we tested the fungal burdens and survival rates of wild-type, rtf1Δ, and complemented strains in intranasal and intravenous murine models of cryptococcosis (Figure 5A). Our results showed that both intranasally and intravenously infected lungs by the rtf1Δ mutant had significantly reduced fungal burden compared to lungs infected by wild-type, full length RTF1-complemented, or HMD-complemented strains (Figure 5B and C). The lungs infected by Plus3-, Rtf1E95A- or HMD E95A-complemented strain had comparable fungal burden relative to the rtf1Δ-infected lungs (Figure 5B and C). Similar trends in effects on fungal burden were observed in other intravenously infected organs, including brain, kidney, and spleen (Figure 5D, E and F). In consistent with the fungal burden analysis, the pathogenicity of rtf1Δ mutant in the intravenous model of cryptococcosis were significantly attenuated compared to the wild-type strain and strains complemented with the full length of Rtf1 or the HMD domain, while the Plus3, Rtf1E95A, or HMD E95A failed to complement the attenuated virulence of the rtf1Δ mutant (Figure 5G). Together, our findings suggest that HMD-mediated H2Bub1 is essential for the successful survival and proliferation of C. neoformans during infection.

Rtf1 and HMD domain regulate cryptococcal pathogenicity in murine models of cryptococcosis.

(A) Schematic diagram of the intranasal and intravenous infection models of cryptococcosis. The inoculum and time for detecting fungal burden in these two infection models were indicated, respectively. (B) The fungal burden of lungs infected by indicated strains through intranasal infection. (C-F) The fungal burden of lungs, brain, kidney, and spleen infected by indicated strains through intravenous infection. (G) The survival curve of animals infected by indicated strains through intravenous infection. The inoculum was the same as the fungal burden assay for intravenous infection. (H) Schematic diagram of the working model depicting the role of Rtf1 and its HMD domain in regulating fungal morphogenesis in C. deneoformans.

4. Discussion

In this study, we investigated the role of Rtf1 in promoting H2Bub1 and consequently regulating cryptococcus yeast-to-hypha transition and virulence. Here, we demonstrated that the global H2Bub1 plays pleiotropic roles in the sexual reproduction, morphogenesis, melanin production, thermal tolerance, and pathogenicity of C. deneoformans and C. neoformans. Interestingly, the Rtf1 HMD domain alone is sufficient to facilitate global H2Bub1 and subsequent H3K4me. The HMD domain could fully restore the deficiencies on filamentation in vitro and pathogenicity in a murine model of cryptococcosis. Our results fit a model in which Rtf1 facilitates the global H2Bub1 and subsequent H3K4me levels, in order to promote expression of genes involved in morphogenesis and pathogenicity in C. neoformans.

Paf1C was first identified as the RNA polymerase II transcriptional regulator functioning in transcription elongation, and is also required for Rad6/Bre1-mediated H2Bub1 and subsequent H3K4me. Whether and how these two roles interplay with each other remain unclear. In yeast, Paf1C contains five highly conserved core subunits Paf1, Leo1, Ctr9, Cdc73, and Rtf1, which is stably associated with the other subunits within Paf1C. However, the human core Paf1C was shown to interact with RNA polymerase II, in the absence of human Rtf1 homolog, indicating the dispensable role of human Rtf1 in the function of Paf1C [61]. The Paf1C is conserved and consists of five subunits in C. neoformans [41]. To investigate the association of Rtf1 with Paf1C in C. neoformans, we conducted co-immunoprecipitation coupled with mass spectrometry (CoIP/MS) assays. None of the other four Paf1C subunits could be detected with either FLAG-tagged full length of Rtf1, HMD, or Plus3 as bait (Supplementary Data S2), strongly indicating that Rtf1 is not stably associated with the Paf1C in C. neoformans.

Rtf1 is critical for H2Bub1 levels, and its deletion abolishes global H2Bub1 in both yeast and humans. It is reasoned that Rtf1 may play dual roles in regulating elongation of gene transcription and deposition of H2Bub1. To gain a comprehensive insight into the function of Rtf1 in eukaryotes, we investigated the role of Rtf1 in human fungal pathogens C. deneoformans and C. neoformans that belong to Basidiomycota. Besides its conserved functions in facilitating global H2Bub1, we also found it is required for fungal morphogenesis and pathogenicity. We showed that HMD domain alone is sufficient to restore cryptococcal filamentation and virulence (Figure 1, 3 and Figure 5), concomitant with the restoration of global H2Bub1 levels (Figure 1, 3 and 4). Given that HMD domain alone lacks regions of full length Rtf1 required for its interactions with other Paf1C subunits and transcribed regions of genes [4, 62], our results on the HMD domain support for a model in which the function of Rtf1 in regulating H2Bub1 is uncoupled from interaction with other subunits of Paf1C, and it is required for cryptococcal development and virulence (Figure 5H). Biochemical and biophysical studies on the association of Rtf1 with Paf1C subunits and Rad6/Bre1 would provide further insights into its mode of action in regulating establishment and deposition of H2Bub1.

Rtf1 contains two conserved domains Plus3 and HMD (Figure 1E). The Plus3 domain has been shown to interact with single-stranded DNA, indicating a role for Rtf1 in the elongation bubble during transcription elongation [62]. In addition, Plus3 may also function in facilitating proper positioning of H2Bub1, especially in regions of actively transcribed genes [4]. Here, we showed that the Plus3 domain alone has no effects on global H2Bub1, while the HMD domain alone could facilitate H2Bub1 (Figure 1A). Moreover, the rtf1Δ mutant showed reduced thermal tolerance with growth defects at 37°C and 39°C, compared to the wild-type strain (Figure 4C). The full length of Rtf1 or HMD domain alone fully restored the growth defect of rtf1Δ mutant at 37°C; However, both of them only partially restored the growth defect at 39°C (Figure 4C). In addition, only the full length of Rtf1 restored melanin production in the rtf1Δ mutant (Figure 4D), although the HMD domain alone fully restored the global H2Bub1 levels in C. neoformans (Figure 4A). There are two possibilities that may lead to these observations: (1) H2Bub1 was not properly depositioned with expression of only HMD domain, although the global level of H2Bub1 seems normal; (2) production of the virulence factors may require functions of Rtf1 and/or Paf1C in transcription elongation, which is absent in HMD- complemented strain [4, 18]. These results on the HMD domain in regulating virulence factors provide insights into the function of full-length Rtf1 and interactions with other subunits of Paf1C. A detailed comparison of H2Bub1 occupancies across chromosomes between cells expressing the full length of Rtf1 and HMD alone would be of great interest. In addition, further studies are required to uncover the roles of Paf1C in facilitating proper deposition of H2Bub1 to regulate fungal morphogenesis and pathogenicity.

Acknowledgements

We thank the Zhao lab for their continued interest and ideas. We thank the Big Data Center and Bioinformatics Center at Department of Veterinary Medicine, Henan Agricultural University for providing high performance computing service. This work was supported by National Natural Science Foundation of China (no. 32373093 and 30900880 to Zhao Y) and Henan Agricultural University (no. 30500946 to Zhao Y).

Disclosure statement

No potential conflict of interest was reported by the author(s).

(A) RTF1 gene regulates fungal filamentation in C. deneoformans. (B) The neighbor-joining tree of Plus3 domain in the indicated eukaryotes. (C) ClustalW multiple amino acid sequence alignment of Plus3 domain in the indicated eukaryotes. Cn, C. deneoformans, Ca, Candida albicans, Af, Aspergillus fumigatus, Fg, Fusarium graminearum, Nc, Neurospora crassa, Sp, Schizosaccharomyces pombe, Sc, S. cerevisiae, Dm, Drosophila melanogaster, Hs, Homo sapiens.

Overexpression of Plus3 domain alone failed to rescue the downregulation of genes enriched in sexual reproduction, pheromone signaling, and filamentous growth due to RTF1 deletion.

(A) Venn diagram of significantly enriched (p-value < 0.05) GO terms of downregulated genes in rtf1Δ, rtf1Δ/ΗΜD, and rtf1Δ/Plus3 relative to the wild-type XL280 strain on V8. (B) The 19 GO terms of downregulated genes shared by rtf1Δ and rtf1Δ/Plus3 strains.

HMD domain alone successfully restores the expression of MAT2 and MFα1 in rtf1Δ strain.

Strains, plasmids, primers and antibodies used in this study.

Supplementary Data S1. List of differentially expressed genes in rtf1Δ, rtf1Δ/ΗΜD, and rtf1Δ/Plus3 relative to the wild-type XL280 strain on V8.

Supplementary Data S2. Potential interacting proteins with Rtf1, HMD, and Plus3 identified by CoIP/MS.