Gene replacement and ectopic integration in the zygomycete Mucor circinelloides

Genetic Manipulation in Mucorales and New Developments to Study Mucormycosis

The study of the Mucoralean fungi physiology is a neglected field that the lack of effective genetic tools has hampered in the past. However, the emerging fungal infection caused by these fungi, known as mucormycosis, has prompted many researchers to study the pathogenic potential of Mucorales. The main reasons for this current attraction to study mucormycosis are its high lethality, the lack of effective antifungal drugs, and its recent increased incidence. The most contemporary example of the emergence character of mucormycosis is the epidemics declared in several Asian countries as a direct consequence of the COVID-19 pandemic. Fortunately, this pressure to understand mucormycosis and develop new treatment strategies has encouraged the blossoming of new genetic techniques and methodologies. This review describes the history of genetic manipulation in Mucorales, highlighting the development of methods and how they allowed the main genetic studies in these fungi. Moreover, we have emphasized the recent development of new genetic models to study mucormycosis, a landmark in the field that will configure future research related to this disease.

Heterotrimeric G-alpha subunits Gpa11 and Gpa12 define a transduction pathway that control spore size and virulence in Mucor circinelloides

Mucor circinelloides is one of the causal agents of mucormycosis, an emerging and high mortality rate fungal infection produced by asexual spores (sporangiospores) of fungi that belong to the order Mucorales. M. circinelloides has served as a model genetic system to understand the virulence mechanism of this infection. Although the G-protein signaling cascade plays crucial roles in virulence in many pathogenic fungi, its roles in Mucorales are yet to be elucidated. Previous study found that sporangiospore size and calcineurin are related to the virulence in Mucor, in which larger spores are more virulent in an animal mucormycosis model and loss of a calcineurin A catalytic subunit CnaA results in larger spore production and virulent phenotype. The M. circinelloides genome is known to harbor twelve gpa (gpa1 to gpa12) encoding G-protein alpha subunits and the transcripts of the gpa11 and gpa12 comprise nearly 72% of all twelve gpa genes transcript in spores. In this study we demonstrated that loss of function of Gpa11 and Gpa12 led to larger spore size associated with reduced activation of the calcineurin pathway. Interestingly, we found lower levels of the cnaA mRNAs in sporangiospores from the Δgpa12 and double Δgpa11/Δgpa12 mutant strains compared to wild-type and the ΔcnaA mutant had significantly lower gpa11 and gpa12 mRNA levels compared to wild-type. However, in contrast to the high virulence showed by the large spores of ΔcnaA, the spores from Δgpa11/Δgpa12 were avirulent and produced lower tissue invasion and cellular damage, suggesting that the gpa11 and gpa12 define a signal pathway with two branches. One of the branches controls spore size through regulation of calcineurin pathway, whereas virulences is controlled by an independent pathway. This virulence-related regulatory pathway could control the expression of genes involved in cellular responses important for virulence, since sporangiospores of Δgpa11/Δgpa12 were less resistant to oxidative stress and phagocytosis by macrophages than the ΔcnaA and wild-type strains. The characterization of this pathway could contribute to decipher the signals and mechanism used by Mucorales to produce mucormycosis.

Functional analysis of the Phycomyces carRA gene encoding the enzymes phytoene synthase and lycopene cyclase

Phycomyces carRA gene encodes a protein with two domains. Domain R is characterized by red carR mutants that accumulate lycopene. Domain A is characterized by white carA mutants that do not accumulate significant amounts of carotenoids. The carRA-encoded protein was identified as the lycopene cyclase and phytoene synthase enzyme by sequence homology with other proteins. However, no direct data showing the function of this protein have been reported so far. Different Mucor circinelloides mutants altered at the phytoene synthase, the lycopene cyclase or both activities were transformed with the Phycomyces carRA gene. Fully transcribed carRA mRNA molecules were detected by Northern assays in the transformants and the correct processing of the carRA messenger was verified by RT-PCR. These results showed that Phycomyces carRA gene was correctly expressed in Mucor. Carotenoids analysis in these transformants showed the presence of ß-carotene, absent in the untransformed strains, providing functional evidence that the Phycomyces carRA gene complements the M. circinelloides mutations. Co-transformation of the carRA cDNA in E. coli with different combinations of the carotenoid structural genes from Erwinia uredovora was also performed. Newly formed carotenoids were accumulated showing that the Phycomyces CarRA protein does contain lycopene cyclase and phytoene synthase activities. The heterologous expression of the carRA gene and the functional complementation of the mentioned activities are not very efficient in E. coli. However, the simultaneous presence of both carRA and carB gene products from Phycomyces increases the efficiency of these enzymes, presumably due to an interaction mechanism.

4-dihydrotrisporin-dehydrogenase, an enzyme of the sex hormone pathway of Mucor mucedo: purification, cloning of the corresponding gene, and developmental expression

The NADP-dependent 4-dihydrotrisporin-dehydrogenase is a (-) mating-type-specific enzyme in the pathway from beta-carotene to trisporic acid. This substance and its isomers and derivatives represent the general system of sexual communication in zygomycetes. The (-) mating type of Mucor mucedo was stimulated by trisporic acid and the enzyme was purified by ion exchange and affinity chromatography. Several peptides of the 26-kDa protein, digested with trypsin, were sequenced by mass spectrometry. Oligonucleotides based on protein sequence data were used for PCR amplification of genomic DNA. The primary PCR fragment was sequenced and the complete gene, TSP2, was isolated. A labeled TSP2 hybridization probe detects a single-copy gene in the genome of M. mucedo. Northern blot analysis with RNAs from different growth stages reveals that the expression of the gene depends on the developmental stage of the mycelium in both mating types of M. mucedo. At the enzyme level, activity is found exclusively in the (-) mating type. However, renaturation of proteins in sodium dodecyl sulfate-containing gels revealed the TSP2 gene product in both mating types. Analyzing the protein sequence places the enzyme in the short chain dehydrogenase superfamily. Thus, it has an evolutionary origin distinct from that of the previously isolated 4-dihydromethyltrisporate dehydrogenase, which belongs to the aldo/keto reductase superfamily. Apart from the TSP2 genes in the three sequenced zygomycetous genomes (Phycomyces blakesleeanus, Rhizopus oryzae, and Mucor circinelloides), the closest relative is the Myxococcus xanthus CsgA gene product, which is also a short chain dehydrogenase, involved in C signaling and fruiting body formation.

Heterologous expression of astaxanthin biosynthesis genes in Mucor circinelloides

Most Mucor species accumulate beta-carotene as the main carotenoid. The crtW and crtZ astaxanthin biosynthesis genes from Agrobacterium aurantiacum were placed under the control of Mucor circinelloides expression signals. Expression vectors containing the bacterial genes were constructed, and PEG-mediated transformations were performed on a selected M. circinelloides strain. Transformants that exhibited altered carotene production were isolated and analyzed. Southern analysis showed that all plasmids behave as autoreplicative elements. Northern analysis detected the actual heterologous transcription products, whereas thin layer chromatography and high-performance liquid chromatography studies revealed the presence of new carotenoid compounds and intermediates among the transformants.

Blakeslea trispora genes for carotene biosynthesis

We cloned the carB and carRA genes involved in beta-carotene biosynthesis from overproducing and wild-type strains of Blakeslea trispora. The carB gene has a length of 1,955 bp, including two introns of 141 and 68 bp, and encodes a protein of 66.4 kDa with phytoene dehydrogenase activity. The carRA gene contains 1,894 bp, with a single intron of 70 bp, and encodes a protein of 69.6 kDa with separate domains for lycopene cyclase and phytoene synthase. The estimated transcript sizes for carB and carRA were 1.8 and 1.9 kb, respectively. CarB from the beta-carotene-overproducing strain B. trispora F-744 had an S528R mutation and a TAG instead of a TAA stop codon. The overproducing strain also had a P143S mutation in CarRA. Both B. trispora genes could complement mutations in orthologous genes in Mucor circinelloides and could be used to construct transformed strains of M. circinelloides that produced higher levels of beta-carotene than did the nontransformed parent. The results show that these genes are conserved across the zygomycetes and that the B. trispora carB and carRA genes are functional and potentially useable to increase carotenoid production.

Genetic damage following introduction of DNA in Phycomyces

Introduction of plasmids in Phycomyces blakesleeanus caused extensive changes in the exogenous DNA and in the resident genome. Plasmids with a bacterial gene for geneticin resistance under a Phycomyces promoter were either injected into immature sporangia or incubated with spheroplasts. An improved method produced about one viable spheroplast per cell. Colonies resistant to geneticin were rare and only about 0.1% of their spores grew in the presence of geneticin. The transformation frequency was very low, < or =1 transformed colony per million spheroplasts or per microg DNA. Few nuclei in the transformants contained exogenous DNA, as shown by a selective procedure that sampled single nuclei from heterokaryons. The exogenous DNA was not integrated into the genome and no stable transformants were obtained. The plasmids were replicated in the recipient cells, but their DNA sequences were modified by deletions and rearrangements and the transformed phenotype was eventually lost. The spores developed in injected sporangia were often inviable; a genetic test showed that spore death was caused by impaired nuclear proliferation and induction of lethal mutations. About one-fourth of the viable spores from injected sporangia formed abnormal colonies with obvious changes in shape, texture, or color. The abnormalities that could be investigated were due to dominant mutations. The results indicate that incoming DNA is not only attacked, but signals a situation of stress that leads to increased mutation and nuclear and cellular death.

A multicopy vector system for genetic studies in Mucor circinelloides and other zygomycetes

Transformation of Mucor circinelloides is routinely achieved by using a plasmid containing the wild-type leuA gene to complement the leucine requirement of an auxotrophic host strain. As is the case for other zygomycetes, the transforming DNA is usually not integrated into the genome of M. circinelloides, but is maintained as an autonomously replicating plasmid. However, even under selective conditions, the plasmid is segregationally unstable, resulting in a rather low number of cells carrying the plasmid. We report here on a new transformation vector based on a dominant selection marker conferring resistance to geneticin, which allows for plasmid maintenance in high copy numbers. The vector was also used to transform Mucor rouxii and Rhizomucor pusillus, and should therefore be a valuable tool for gene expression studies in zygomycetes. The functionality and regulatory properties of the promoter of the M. circinelloides gpd1 gene (which codes for glyceraldehyde-3P-dehydrogenase) were demonstrated in R. pusillus using geneticin selection. In this work, we have also determined the molecular basis of the Leu(-) phenotype of the M. circinelloides host strain R7B. The leucine requirement is due to a single point mutation in the leuA gene that results in the replacement of a glutamic acid by a lysine residue.

Development of a system for integrative and stable transformation of the zygomycete Rhizopus oryzae by Agrobacterium-mediated DNA transfer

Two transformation systems, based on the use of CaCl(2)/PEG and Agrobacterium tumefaciens, respectively, were developed for the zygomycete Rhizopus oryzae. Irrespective of the selection marker used, a pyr4 marker derived from R. niveus or a dominant amdS(+) marker from Aspergillus nidulans, and irrespective of the configuration of the transforming DNA (linear or circular), the transformants obtained with the CaCl(2)/PEG transformation method were found to carry multiple copies of tandemly linked vector molecules, which failed to integrate into the genomic DNA. Furthermore, these transformants displayed low mitotic stability. In contrast, transformants obtained by Agrobacterium-mediated transformation were mitotically stable, even under non-selective conditions. Detailed analysis of these transformants revealed that the transforming DNA had integrated into the genome of R. oryzae at a single locus in independently obtained transformants. In addition, truncation of the transforming DNA was observed, resulting in the integration of the R. niveus pyr4 marker gene, but not the second gene located on the transferred DNA. Modification of the transforming DNA, resulting in partial resistance to restriction enzyme digestion, was observed in transformants obtained with the CaCl(2)/PEG transformation method, suggesting that a specific genome defence mechanism may exist in R. oryzae. It is likely that the unique mechanism used by A. tumefaciens to deliver its transferred DNA to its hosts facilitates bypass of the host defence mechanisms, thus allowing the DNA to integrate into the chromosomal genome.

Transformation of Mucor miehei results in plasmid deletion and phenotypic instability

Mucor miehei transformants resistant to geneticin have been obtained by treatment of protoplasts with different plasmids and by Agrobacterium-mediated DNA transfer. All transformants exhibited a very unstable phenotype whose maintenance required continuous selective pressure. Molecular analysis of transformants showed that the plasmid DNA was extensively modified and maintained in very low amount. Our results indicate that M. miehei reluctance to transformation is due to different causes, including the coenocytic nature of its mycelium and the existence of specific mechanisms for the detection and elimination of foreign DNA.

Cloning of glyceraldehyde-3-phosphate dehydrogenase-encoding genes in Mucor circinelloides (Syn. racemosus) and use of the gpd1 promoter for recombinant protein production

Three genes (gpd1, gpd2, and gpd3) encoding glyceraldehyde-3-phosphate dehydrogenase were isolated from the dimorphic zygomycete Mucor circinelloides by PCR using degenerated primers. Transcription of gpd1 could be detected during vegetative growth under both aerobic and anaerobic conditions, whereas neither gpd2 nor gpd3 transcription was detected, indicating that gpd1 is the major transcribed gpd gene. The transcription of gpd1 was regulated by carbon source. The gpd1 promoter was successfully used for recombinant expression of genes of both homologous (crgA encoding a regulator of carotene biosynthesis) and heterologous (gox1 from Aspergillus niger encoding glucose oxidase; GOX) nature. Growth of a gox1 transformant strain resulted in the secretion of enzymatically active GOX. The potential advantages of using a dimorphic fungus for heterologous protein production are discussed.

Isolation and use of a homologous histone H4 promoter and a ribosomal DNA region in a transformation vector for the oil-producing fungus Mortierella alpina

Mortierella alpina was transformed successfully to hygromycin B resistance by using a homologous histone H4 promoter to drive gene expression and a homologous ribosomal DNA region to promote chromosomal integration. This is the first description of transformation in this commercially important oleaginous organism. Two pairs of histone H3 and H4 genes were isolated from this fungus. Each pair consisted of one histone H3 gene and one histone H4 gene, transcribed divergently from an intergenic promoter region. The pairs of encoded histone H3 or H4 proteins were identical in amino acid sequence. At the DNA level, each histone H3 or H4 open reading frame showed 97 to 99% identity to its counterpart but the noncoding regions had little sequence identity. Unlike the histone genes from other filamentous fungi, all four M. alpina genes lacked introns. During normal vegetative growth, transcripts from the two histone H4 genes were produced at approximately the same level, indicating that either histone H4 promoter could be used in transformation vectors. The generation of stable, hygromycin B-resistant transformants required the incorporation of a homologous ribosomal DNA region into the transformation vector to promote chromosomal integration.

Isolation, characterization and transformation, by autonomous replication, of Mucor circinelloides OMPdecase-deficient mutants

Pyrimidine auxotrophs of Mucor circinelloides were isolated after mutagenesis with nitrosoguanidine and selected for resistance to 5-fluoroorotate. These mutants were genetically and biochemically characterized and found to be deficient either in orotidine-5'-monophosphate decarboxylase (OMPdecase) activity or in orotate phosphoribosyl transferase (OPRTase) activity. Different circular DNA molecules containing the homologous pyrG gene were used to transform a representative OMPdecase-deficient strain to uracil prototrophy. Southern analysis, as well as mitotic stability analysis of the transformants, showed that the transforming DNA is always maintained extrachromosomally. The smallest fragment tested that retained both the capacity to complement the pyrG4 mutation and the ability to be maintained extrachromosomally when cloned in a suitable vector is a 1.85 kb M. circinelloides genomic DNA fragment. This fragment consists of the pyrG coding region flanked by 606 nucleotides at the 5' and 330 nucleotides at the 3' ends, respectively. Sequence analysis reveals that it does not share any element in common with another M. circinelloides genomic DNA fragment which also promotes autonomous replication in this organism, except those related to transcription. Furthermore, it differs from elements which have been shown to be involved in autonomous replication in other fungal systems. An equivalent plasmid harbouring the heterologous Phycomyces blakesleeanus pyrG gene yielded lower transformation rates, but the transforming DNA was also maintained extrachromosomally. Our results suggest that autonomous replication in M. circinelloides may be driven by elements normally present in nuclear coding genes.

Electrophoretic karyotype of Mucor circinelloides

Contour-clamped homogeneous electric field (CHEF) gel electrophoresis was used to separate chromosomal size DNA molecules of two Mucor circinelloides strains. Electrophoretic karyotypes revealed the presence of eight distinct bands for the M. circinelloides f. lusitanicus strain, and four, presumably multiple, bands for the M. circinelloides f. gryseo-cyanus strain. The approximate sizes of the resolved chromosomal DNA bands ranged from 2.3 to 8.1 Mb, giving estimated genome sizes of 38.7 and 32.6 Mb, respectively. Hybridisation techniques were used to assign the leuA gene to a chromosome.