Uncontrolled transposition following RNAi loss causes hypermutation and antifungal drug resistance in clinical isolates of Cryptococcus neoformans

Cryptococcus neoformans infections cause approximately 15% of AIDS-related deaths owing to a combination of limited antifungal therapies and drug resistance. A collection of clinical and environmental C. neoformans isolates were assayed for increased mutation rates via fluctuation analysis, and we identified two hypermutator C. neoformans clinical isolates with increased mutation rates when exposed to the combination of rapamycin and FK506. Sequencing of drug target genes found that Cnl1 transposon insertions conferred the majority of resistance to rapamycin and FK506 and could also independently cause resistance to 5-fluoroorotic acid and the clinically relevant antifungal 5-flucytosine. Whole-genome sequencing revealed both hypermutator genomes harbour a nonsense mutation in the RNA-interference component ZNF3 and hundreds of Cnl1 elements organized into massive subtelomeric arrays on each of the fourteen chromosomes. Quantitative trait locus mapping in 28 progeny derived from a cross between a hypermutator and wild-type identified a locus associated with hypermutation that included znf3. CRISPR editing of the znf3 nonsense mutation abolished hypermutation and restored small-interfering-RNA production. We conclude that hypermutation and drug resistance in these clinical isolates result from RNA-interference loss and accumulation of Cnl1 elements.

The master regulator MAT1-1-1 of fungal mating binds to its targets via a conserved motif in the human pathogen Aspergillus fumigatus

Mating-type transcription factors are master regulators of sexually related signal transduction pathways in fungi; however, their recognition of specific DNA sequences from target genes is widely undetermined. Here, we identified and characterized the DNA-binding sequence of the MAT1-1-1 alpha-box domain transcription factor from the human pathogen Aspergillus fumigatus. In order to explore MAT1-1-1 DNA-binding targets, we used the previously reported MAT1-1-1 binding motif from Penicillium chrysogenum, in a bioinformatics approach. We identified 18 A. fumigatus genes carrying the MAT1.1 sequence in their upstream region, among them genes for the α-pheromone precursor (PpgA), G-protein-coupled pheromone receptor (PreA), and for TomA, an unidentified protein. To validate our prediction further, quantification of transcript levels showed a decrease in expression of ppgA, tomA, and others in a MAT1-1 deletion strain. For a functional analysis of the binding sites, truncated variants of the A. fumigatus MAT1-1-1 gene were introduced into Escherichia coli for heterologous expression. The yield of recombinant protein was further optimized for the AfMAT1-1-178-235 variant that harbors an extended alpha-box domain. AfMAT1-1-178-235 bound to a subset of the most strongly upregulated genes: ppgA, preA, and tomA. The DNA-binding specificity was confirmed by testing mutated binding sequences, as well as performing competition experiments with specific and non-specific sequences. Finally, equilibrium dissociation constants of 1.83 ± 0.1 and 1.45 ± 0.26 µM were determined for AfMAT1-1-178-235 and fusion protein GST-AfMAT1-1-178-235. Collectively, these findings provide further insights into AfMAT1-1-1-mediated gene expression and imply that alpha-box domain regulators from other members of Eurotiales control fungal development in a conserved manner.

Molecular analysis of mating type loci from the mycophenolic acid producer Penicillium brevicompactum: Phylogeny and MAT protein characterization suggest a cryptic sexual life cycle

The mycophenolic acid producing ascomycete Penicillium brevicompactum is considered to be an anamorphic (asexual) species, for which a sexual cycle was never observed. However, since recent reports of otherwise asexually propagating filamentous fungi have demonstrated a sexual cycle controlled by mating type loci, we carried out a molecular analysis of mating type loci from P. brevicompactum. Using data from extensive DNA sequencing analysis, we determined the mating type loci from 22 strains derived from various type culture collections. We found 8 strains carrying a MAT1-1 locus encoding a 362 amino acid alpha domain transcription factor. The other 14 possessed a MAT1-2 locus encoding a 298 amino acid HMG domain transcription factor. cDNA analysis confirmed that both mating type loci are transcriptionally expressed. The karyotype of six selected strains, determined using contour-clamped homogeneous electric field (CHEF) electrophoresis, demonstrated distinct differences in size and numbers of chromosomes between the strains investigated. Interestingly, our phylogenetic survey of 72 strains from 11 different Penicillium species revealed that MAT genes serve as excellent molecular markers to determine phylogenetic relationships among species closely related to P. brevicompactum. Based on our sequencing results, we constructed transformation vectors for site-specific deletion of mating type loci from two selected strains of opposite mating type. Complementation strains were constructed containing both the mating type locus deletion cassette and a MAT-egfp fusion gene. These strains were used for comparative phenotypic analyses between strains containing or lacking the mating type gene. Whereas all MAT1-2 strains were indistinguishable, the MAT1-1 and MAT1-1-1 deletion strains differed distinctly. The MAT1-1-1 deletion strain produced more conidiospores on solid media, but smaller pellets in liquid media. This is probably the consequence of fewer conidial germ tubes than with the wild type mating type strain. Finally, we showed that the MAT-EGPF fusion protein is localized to the nuclei and detectable in protein samples by Western analysis. Together, our results suggest that the asexually propagating fungus P. brevicompactum might be a heterothallic species with a cryptic sexual life cycle.

Mating-Type-Specific Ribosomal Proteins Control Aspects of Sexual Reproduction in Cryptococcus neoformans

The MAT locus of Cryptococcus neoformans has a bipolar organization characterized by an unusually large structure, spanning over 100 kb. MAT genes have been characterized by functional genetics as being involved in sexual reproduction and virulence. However, classical gene replacement failed to achieve mutants for five MAT genes (RPL22, RPO41, MYO2, PRT1, and RPL39), indicating that they are likely essential. In the present study, targeted gene replacement was performed in a diploid strain for both the α and a alleles of the ribosomal genes RPL22 and RPL39 Mendelian analysis of the progeny confirmed that both RPL22 and RPL39 are essential for viability. Ectopic integration of the RPL22 allele of opposite MAT identity in the heterozygous RPL22a/rpl22αΔ or RPL22α/rpl22aΔ mutant strains failed to complement their essential phenotype. Evidence suggests that this is due to differential expression of the RPL22 genes, and an RNAi-dependent mechanism that contributes to control RPL22a expression. Furthermore, via CRISPR/Cas9 technology, the RPL22 alleles were exchanged in haploid MATα and MATa strains of C. neoformans These RPL22 exchange strains displayed morphological and genetic defects during bilateral mating. These results contribute to elucidating functions of C. neoformans essential mating type genes that may constitute a type of imprinting system to promote inheritance of nuclei of both mating types.

RNA Editing During Sexual Development Occurs in Distantly Related Filamentous Ascomycetes

RNA editing is a post-transcriptional process that modifies RNA molecules leading to transcript sequences that differ from their template DNA. A-to-I editing was found to be widely distributed in nuclear transcripts of metazoa, but was detected in fungi only recently in a study of the filamentous ascomycete Fusarium graminearum that revealed extensive A-to-I editing of mRNAs in sexual structures (fruiting bodies). Here, we searched for putative RNA editing events in RNA-seq data from Sordaria macrospora and Pyronema confluens, two distantly related filamentous ascomycetes, and in data from the Taphrinomycete Schizosaccharomyces pombe. Like F. graminearum, S. macrospora is a member of the Sordariomycetes, whereas P. confluens belongs to the early-diverging group of Pezizomycetes. We found extensive A-to-I editing in RNA-seq data from sexual mycelium from both filamentous ascomycetes, but not in vegetative structures. A-to-I editing was not detected in different stages of meiosis of S. pombe. A comparison of A-to-I editing in S. macrospora with F. graminearum and P. confluens, respectively, revealed little conservation of individual editing sites. An analysis of RNA-seq data from two sterile developmental mutants of S. macrospora showed that A-to-I editing is strongly reduced in these strains. Sequencing of cDNA fragments containing more than one editing site from P. confluens showed that at the beginning of sexual development, transcripts were incompletely edited or unedited, whereas in later stages transcripts were more extensively edited. Taken together, these data suggest that A-to-I RNA editing is an evolutionary conserved feature during fruiting body development in filamentous ascomycetes.

Transcriptome analysis of the two unrelated fungal β-lactam producers Acremonium chrysogenum and Penicillium chrysogenum: Velvet-regulated genes are major targets during conventional strain improvement programs

BACKGROUND

Cephalosporins and penicillins are the most frequently used β-lactam antibiotics for the treatment of human infections worldwide. The main industrial producers of these antibiotics are Acremonium chrysogenum and Penicillium chrysogenum, two taxonomically unrelated fungi. Both were subjects of long-term strain development programs to reach economically relevant antibiotic titers. It is so far unknown, whether equivalent changes in gene expression lead to elevated antibiotic titers in production strains.

RESULTS

Using the sequence of PcbC, a key enzyme of β-lactam antibiotic biosynthesis, from eighteen different pro- and eukaryotic microorganisms, we have constructed a phylogenetic tree to demonstrate the distant relationship of both fungal producers. To address the question whether both fungi have undergone similar genetic adaptions, we have performed a comparative gene expression analysis of wild-type and production strains. We found that strain improvement is associated with the remodeling of the transcriptional landscape in both fungi. In P. chrysogenum, 748 genes showed differential expression, while 1572 genes from A. chrysogenum are differentially expressed in the industrial strain. Common in both fungi is the upregulation of genes belonging to primary and secondary metabolism, notably those involved in precursor supply for β-lactam production. Other genes not essential for β-lactam production are downregulated with a preference for those responsible for transport processes or biosynthesis of other secondary metabolites. Transcriptional regulation was shown to be an important parameter during strain improvement in different organisms. We therefore investigated deletion strains of the major transcriptional regulator velvet from both production strains. We identified 567 P. chrysogenum and 412 A. chrysogenum Velvet target genes. In both deletion strains, approximately 50% of all secondary metabolite cluster genes are differentially regulated, including β-lactam biosynthesis genes. Most importantly, 35-57% of Velvet target genes are among those that showed differential expression in both improved industrial strains.

CONCLUSIONS

The major finding of our comparative transcriptome analysis is that strain improvement programs in two unrelated fungal β-lactam antibiotic producers alter the expression of target genes of Velvet, a global regulator of secondary metabolism. From these results, we conclude that regulatory alterations are crucial contributing factors for improved β-lactam antibiotic titers during strain improvement in both fungi.

Sexual recombination as a tool for engineering industrial Penicillium chrysogenum strains

The recent discovery and functional characterization of opposite mating-type loci in the industrial penicillin producer Penicillium chrysogenum demonstrated their regulatory role in sexual as well as asexual development. Subsequent experiments further showed that a sexual life cycle can be induced in P. chrysogenum that was for long believed to reproduce exclusively by asexual propagation. Finally, crossing of wild type and production strains resulted in the generation of recombinant ascospore isolates. We predict from these recent findings that recombinant progeny for industrial applications can be obtained by sexual crossings and discuss experimental difficulties that occur when parental strains with karyotype heterogeneity are used for mating.

Dicer-Dependent Biogenesis of Small RNAs and Evidence for MicroRNA-Like RNAs in the Penicillin Producing Fungus Penicillium chrysogenum

MicroRNAs (miRNAs) are non-coding small RNAs (sRNAs) that regulate gene expression in a wide range of eukaryotes. In this study, we analyzed regulatory sRNAs in Penicillium chrysogenum, the industrial producer of the β-lactam antibiotic penicillin. To identify sRNAs and microRNA-like RNAs (milRNAs) on a global approach, two sRNA sequencing libraries were constructed. One library was created with pooled total RNA, obtained from twelve differently grown cultures (RNA Mix), and the other with total RNA from a single submerged cultivation (∆ku70FRT2). Illumina sequencing of both RNA libraries produced 84,322,825 mapped reads. To distinguish between Dicer-dependent and independent sRNA formation, we further constructed two single dicer gene mutants (∆dcl2 and ∆dcl1) and a dicer double mutant (∆dcl2∆dcl1) and analyzed an sRNA library from the Dicer-deficient double-mutant. We identified 661 Dicer-dependent loci and in silico prediction revealed 34 milRNAs. Northern blot hybridization of two milRNAs provided evidence for mature milRNAs that are processed either in a complete or partial Dicer-dependent manner from an RNA precursor. Identified milRNAs share typical characteristics of previously discovered fungal milRNAs, like a strong preference for a 5' uracil and the typical length distribution. The detection of potential milRNA target sites in the genome suggests that milRNAs might play a role in posttranscriptional gene regulation. Our data will further increase our knowledge of sRNA dependent gene regulation processes, which is an important prerequisite to develop more effective strategies for improving industrial fermentations with P. chrysogenum.

A MAT1-2 wild-type strain from Penicillium chrysogenum: functional mating-type locus characterization, genome sequencing and mating with an industrial penicillin-producing strain

In heterothallic ascomycetes, mating is controlled by two nonallelic idiomorphs that determine the 'sex' of the corresponding strains. We recently discovered mating-type loci and a sexual life cycle in the penicillin-producing fungus, Penicillium chrysogenum. All industrial penicillin production strains worldwide are derived from a MAT1-1 isolate. No MAT1-2 strain has been investigated in detail until now. Here, we provide the first functional analysis of a MAT1-2 locus from a wild-type strain. Similar to MAT1-1, the MAT1-2 locus has functions beyond sexual development. Unlike MAT1-1, the MAT1-2 locus affects germination and surface properties of conidiospores and controls light-dependent asexual sporulation. Mating of the MAT1-2 wild type with a MAT1-1 high penicillin producer generated sexual spores. We determined the genomic sequences of parental and progeny strains using next-generation sequencing and found evidence for genome-wide recombination. SNP calling showed that derived industrial strains had an uneven distribution of point mutations compared with the wild type. We found evidence for meiotic recombination in all chromosomes. Our results point to a strategy combining the use of mating-type genes, genetics, and next-generation sequencing to optimize conventional strain improvement methods.