The mitochondrial genome of the thermal dimorphic fungus Penicillium marneffei is more closely related to those of molds than yeasts

An Overlooked and Underrated Endemic Mycosis-Talaromycosis and the Pathogenic Fungus Talaromyces marneffei

Talaromycosis is an invasive mycosis endemic in tropical and subtropical Asia and is caused by the pathogenic fungus Talaromyces marneffei. Approximately 17,300 cases of T. marneffei infection are diagnosed annually, and the reported mortality rate is extremely high (~1/3). Despite the devastating impact of talaromycosis on immunocompromised individuals, particularly HIV-positive persons, and the increase in reported occurrences in HIV-uninfected persons, diagnostic and therapeutic approaches for talaromycosis have received far too little attention worldwide. In 2021, scientists living in countries where talaromycosis is endemic raised a global demand for it to be recognized as a neglected tropical disease. Therefore, T. marneffei and the infectious disease induced by this fungus must be treated with concern. T. marneffei is a thermally dimorphic saprophytic fungus with a complicated mycological growth process that may produce various cell types in its life cycle, including conidia, hyphae, and yeast, all of which are associated with its pathogenicity. However, understanding of the pathogenic mechanism of T. marneffei has been limited until recently. To achieve a holistic view of T. marneffei and talaromycosis, the current knowledge about talaromycosis and research breakthroughs regarding T. marneffei growth biology are discussed in this review, along with the interaction of the fungus with environmental stimuli and the host immune response to fungal infection. Importantly, the future research directions required for understanding this serious infection and its causative pathogenic fungus are also emphasized to identify solutions that will alleviate the suffering of susceptible individuals worldwide.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

Mp1p homologues as virulence factors in Aspergillus fumigatus

Recently, we showed that Mp1p is an important virulence factor of Talaromyces marneffei, a dimorphic fungus phylogenetically closely related to Aspergillus fumigatus. In this study, we investigated the virulence properties of the four Mp1p homologues (Afmp1p, Afmp2p, Afmp3p, and Afmp4p) in A. fumigatus using a mouse model. All mice died 7 days after challenge with wild-type A. fumigatus QC5096, AFMP1 knockdown mutant, AFMP2 knockdown mutant and AFMP3 knockdown mutant and 28 days after challenge with AFMP4 knockdown mutant (P<.0001). Only 11% of mice died 30 days after challenge with AFMP1-4 knockdown mutant (P<.0001). For mice challenge with AFMP1-4 knockdown mutant, lower abundance of fungal elements was observed in brains, kidneys, and spleens compared to mice challenge with QC5096 at day 4 post-infection. Fungal counts in brains of mice challenge with QC5096 or AFMP4 knockdown mutant were significantly higher than those challenge with AFMP1-4 knockdown mutant (P<.01 and P<.05). Fungal counts in kidneys of mice challenge with QC5096 or AFMP4 knockdown mutant were significantly higher than those challenge with AFMP1-4 knockdown mutant (P<.001 and P<.001) and those of mice challenge with QC5096 were significantly higher than those challenge with AFMP4 knockdown mutant (P<.05). There is no difference among the survival rates of wild-type A. fumigatus, AFMP4 knockdown mutant and AFMP1-4 knockdown mutant, suggesting that Mp1p homologues in A. fumigatus do not mediate its virulence via improving its survival in macrophage as in the case in T. marneffei. Afmp1p, Afmp2p, Afmp3p, and Afmp4p in combination are important virulence factors of A. fumigatus.

Novel Partitivirus Enhances Virulence of and Causes Aberrant Gene Expression in Talaromyces marneffei

Talaromyces marneffei is the most important thermal dimorphic fungus causing systemic mycosis in Southeast Asia. We report the discovery of a novel partitivirus, Talaromyces marneffeipartitivirus-1 (TmPV1). TmPV1 was detected in 7 (12.7%) of 55 clinical T. marneffei isolates. Complete genome sequencing of the seven TmPV1 isolates revealed two double-stranded RNA (dsRNA) segments encoding RNA-dependent RNA polymerase (RdRp) and capsid protein, respectively. Phylogenetic analysis showed that TmPV1 occupied a distinct clade among the members of the genus Gammapartitivirus Transmission electron microscopy confirmed the presence of isometric, nonenveloped viral particles of 30 to 45 nm in diameter, compatible with partitiviruses, in TmPV1-infected T. marneffei Quantitative reverse transcription-PCR (qRT-PCR) demonstrated higher viral load of TmPV1 in the yeast phase than in the mycelial phase of T. marneffei Two virus-free isolates, PM1 and PM41, were successfully infected by purified TmPV1 using protoplast transfection. Mice challenged with TmPV1-infected T. marneffei isolates showed significantly shortened survival time (P < 0.0001) and higher fungal burden in organs than mice challenged with isogenic TmPV1-free isolates. Transcriptomic analysis showed that TmPV1 causes aberrant expression of various genes in T. marneffei, with upregulation of potential virulence factors and suppression of RNA interference (RNAi)-related genes. This is the first report of a mycovirus in a thermally dimorphic fungus. Further studies are required to ascertain the mechanism whereby TmPV1 enhances the virulence of T. marneffei in mice and the potential role of RNAi-related genes in antiviral defense in T. marneffeiIMPORTANCETalaromyces marneffei (formerly Penicillium marneffei) is the most important thermal dimorphic fungus in Southeast Asia, causing highly fatal systemic penicilliosis in HIV-infected and immunocompromised patients. We discovered a novel mycovirus, TmPV1, in seven clinical isolates of T. marneffei TmPV1 belongs to the genus Gammapartitivirus of the family Partitiviridae We showed that TmPV1 enhanced the virulence of T. marneffei in mice, with shortened survival time and higher fungal burden in the organs of mice challenged with TmPV1-infected T. marneffei isolates than in those of mice challenged with virus-free isogenic isolates. Transcriptomics analysis showed that TmPV1 altered the expression of genes involved in various cellular processes in T. marneffei, with upregulation of potential virulence factors and suppression of RNAi machinery which may be involved in antiviral defense. This is the first report of a mycovirus in a thermal dimorphic fungus. The present results offer insights into mycovirus-fungus interactions and pathogenesis of thermal dimorphic fungi.

Taxonomy and evolution of Aspergillus, Penicillium and Talaromyces in the omics era - Past, present and future

Aspergillus, Penicillium and Talaromyces are diverse, phenotypically polythetic genera encompassing species important to the environment, economy, biotechnology and medicine, causing significant social impacts. Taxonomic studies on these fungi are essential since they could provide invaluable information on their evolutionary relationships and define criteria for species recognition. With the advancement of various biological, biochemical and computational technologies, different approaches have been adopted for the taxonomy of Aspergillus, Penicillium and Talaromyces; for example, from traditional morphotyping, phenotyping to chemotyping (e.g. lipotyping, proteotypingand metabolotyping) and then mitogenotyping and/or phylotyping. Since different taxonomic approaches focus on different sets of characters of the organisms, various classification and identification schemes would result. In view of this, the consolidated species concept, which takes into account different types of characters, is recently accepted for taxonomic purposes and, together with the lately implemented 'One Fungus - One Name' policy, is expected to bring a more stable taxonomy for Aspergillus, Penicillium and Talaromyces, which could facilitate their evolutionary studies. The most significant taxonomic change for the three genera was the transfer of Penicillium subgenus Biverticillium to Talaromyces (e.g. the medically important thermally dimorphic 'P. marneffei' endemic in Southeast Asia is now named T. marneffei), leaving both Penicillium and Talaromyces as monophyletic genera. Several distantly related Aspergillus-like fungi were also segregated from Aspergillus, making this genus, containing members of both sexual and asexual morphs, monophyletic as well. In the current omics era, application of various state-of-the-art omics technologies is likely to provide comprehensive information on the evolution of Aspergillus, Penicillium and Talaromyces and a stable taxonomy will hopefully be achieved.

Comparative mitogenomic and phylogenetic characterization on the complete mitogenomes of Talaromyces (Penicillium) marneffei

We report the complete mitochondrial DNA (mtDNA) sequences of four Talaromyces marneffei strains and performed comparative genomic and phylogenetic analyses. The gene orders of the four mtDNAs are identical to the previously published mtDNA of strain PM1 (nad4l, nad5, nad2, atp9, cob, nad1, nad4, atp8, atp6, nad6, cox3, rps, cox1, nad3, cox2). Phylogenetic analysis showed that the four mtDNAs were clustered with that of PM1 with high bootstrap support. Compared to mtDNA of PM1, the only non-synonymous mutation was located in nad2 (T505M) of strain PM26. Synonymous single nucleotide polymorphisms were observed at eight positions in the four mtDNAs.

Expression and characterization of a Talaromyces marneffei active phospholipase B expressed in a Pichia pastoris expression system

Phospholipase B is a virulence factor for several clinically important pathogenic fungi, including Candida albicans, Cryptococcus neoformans and Aspergillus fumigatus, but its role in the thermally dimorphic fungus Talaromyces marneffei remains unclear. Here, we provide the first report of the expression of a novel phospholipase gene, designated TmPlb1, from T. marneffei in the eukaryotic expression system of Pichia pastoris GS115. Sensitive real-time quantitative reverse-transcription PCR (qRT-PCR) demonstrated that the expression of TmPlb1 increased 1.85-fold in the yeast phase compared with the mycelial phase. TmPlb1 contains an open reading frame (ORF) of 732 bp that encodes a protein of 243 amino acids. The conserved serine, aspartate and histidine catalytic triad and the G-X-S-X-G domain of TmPLB1 provide the structural basis for its molecular activity. The ORF of TmPlb1 was successfully cloned into a pPIC9K vector containing an α-mating factor secretion signal that allowed the secretory expression of TmPLB1 in P. pastoris. The heterologous protein expression began 12 h after methanol induction and peaked at 96 h. Through analysis with SDS-polyacrylamide gel electrophoresis (SDS-PAGE), western blotting and mass spectrometry, we confirmed that TmPLB1 was successfully expressed. Through Ni-affinity chromatography, TmPLB1 was highly purified, and its concentration reached 240.4 mg/L of culture medium. With specific substrates, the phospholipase A1 and phospholipase A2 activities of TmPLB1 were calculated to be 5.96 and 1.59 U/mg, respectively. The high purity and activity of the TmPLB1 obtained here lay a solid foundation for further investigation.

Mp1p Is a Virulence Factor in Talaromyces (Penicillium) marneffei

Background

Talaromyces marneffei is an opportunistic dimorphic fungus prevalent in Southeast Asia. We previously demonstrated that Mp1p is an immunogenic surface and secretory mannoprotein of T. marneffei. Since Mp1p is a surface protein that can generate protective immunity, we hypothesized that Mp1p and/or its homologs are virulence factors.

Methodology/Principal Findings

We examined the pathogenic roles of Mp1p and its homologs in a mouse model. All mice died 21 and 30 days after challenge with wild-type T. marneffei PM1 and MP1 complemented mutant respectively. None of the mice died 60 days after challenge with MP1 knockout mutant (P<0.0001). Seventy percent of mice died 60 days after challenge with MP1 knockdown mutant (P<0.0001). All mice died after challenge with MPLP1 to MPLP13 knockdown mutants, suggesting that only Mp1p plays a significant role in virulence. The mean fungal loads of PM1 and MP1 complemented mutant in the liver, lung, kidney and spleen were significantly higher than those of the MP1 knockout mutant. Similarly, the mean load of PM1 in the liver, lung and spleen were significantly higher than that of the MP1 knockdown mutant. Histopathological studies showed an abundance of yeast in the kidney, spleen, liver and lung with more marked hepatic and splenic necrosis in mice challenged with PM1 compared to MP1 knockout and MP1 knockdown mutants. Likewise, a higher abundance of yeast was observed in the liver and spleen of mice challenged with MP1 complemented mutant compared to MP1 knockout mutant. PM1 and MP1 complemented mutant survived significantly better than MP1 knockout mutant in macrophages at 48 hours (P<0.01) post-infection. The mean fungal counts of Pichia pastoris GS115-MP1 in the liver (P<0.001) and spleen (P<0.05) of mice were significantly higher than those of GS115 at 24 hours post-challenge.

Conclusions/Significance

Mp1p is a key virulence factor of T. marneffei. Mp1p mediates virulence by improving the survival of T. marneffei in macrophages.

Talaromyces (Penicillium) marneffei is an opportunistic thermal dimorphic fungus most prevalent in Southeast Asia. Our team has previously shown that Mp1p, a protein encoded by the MP1 gene, is an immunogenic surface and secretory protein of T. marneffei. In this study, we showed that mice challenged with T. marneffei with the MP1 gene died but those challenged with T. marneffei without the MP1 gene did not die. There was also significantly higher fungal load and more necrosis in organs of mice challenged with T. marneffei with the MP1 gene than T. marneffei without the MP1 gene. Furthermore, T. marneffei with the MP1 gene survived better in macrophages than T. marneffei without the MP1 gene and Pichia pastoris with the MP1 gene survived in mice better than P. pastoris without the MP1 gene. Our data support that Mp1p is a key virulence factor of T. marneffei and Mp1p mediates virulence by improving the survival of T. marneffei in macrophages.

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry for rapid identification of mold and yeast cultures of Penicillium marneffei

BACKGROUND

Penicillium marneffei is the most important thermal dimorphic fungus causing systemic mycosis in HIV-infected and other immunocompromised patients in Southeast Asia. However, laboratory diagnosis of penicilliosis, which relies on microscopic morphology and mycelial-to-yeast conversion, is time-consuming and expertise-dependent, thus delaying diagnosis and treatment. Although matrix -assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS) is useful for identification of various medically important fungi, its performance for identification of P. marneffei is less clear.

RESULTS

We evaluated the performance of the Bruker MALDI-TOF MS system for identification of mold and yeast cultures of 59 clinical strains and the type strain of P. marneffei using the direct transfer method, with results compared to four phylogenetically closely related species, P. brevi-compactum, P. chrysogenum, Talaromyces aurantiacus and T. stipitatus. Using the Bruker original database combined with BDAL v4.0.0.1 and Filamentous Fungi Library 1.0, MALDI-TOF MS failed to identify the 60 P. marneffei strains grown in mold and yeast phase (identified as P. funiculosum and P. purpurogenum with scores <1.7 respectively). However, when the combined database was expanded with inclusion of spectra from 21 P. marneffei strains in mold and/or yeast phase, all the remaining 39 P. marneffei strains grown in mold or phase were correctly identified to the species level with score >2.0. The MS spectra of P. marneffei exhibited significant difference to those of P. brevi-compactum, P. chrysogenum, T. aurantiacus and T. stipitatus. However, MALDI-TOF MS failed to identify these four fungi to the species level using the combined database with or without spectra from P. marneffei.

CONCLUSIONS

MALDI-TOF MS is useful for rapid identification of both yeast and mold cultures of P. marneffei and differentiation from related species. However, accurate identification to the species level requires database expansion using P. marneffei strains.

Talaromyces (Penicillium) marneffei infection in non-HIV-infected patients

Talaromyces (Penicillium) marneffei is an important pathogenic thermally dimorphic fungus causing systemic mycosis in Southeast Asia. The clinical significance of T. marneffei became evident when the human immunodeficiency virus (HIV)/acquired immunodeficiency syndrome epidemic arrived in Southeast Asia in 1988. Subsequently, a decline in the incidence of T. marneffei infection among HIV-infected patients was seen in regions with access to highly active antiretroviral therapy and other control measures for HIV. Since the 1990s, an increasing number of T. marneffei infections have been reported among non-HIV-infected patients with impaired cell-mediated immunity. Their comorbidities included primary adult-onset immunodeficiency due to anti-interferon-gamma autoantibodies and secondary immunosuppressive conditions including other autoimmune diseases, solid organ and hematopoietic stem cell transplantations, T-lymphocyte-depleting immunsuppressive drugs and novel anti-cancer targeted therapies such as anti-CD20 monoclonal antibodies and kinase inhibitors. Moreover, improved immunological diagnostics identified more primary immunodeficiency syndromes associated with T. marneffei infection in children. The higher case-fatality rate of T. marneffei infection in non-HIV-infected than HIV-infected patients might be related to delayed diagnosis due to the lack of clinical suspicion. Correction of the underlying immune defects and early use of antifungals are important treatment strategies. Clinicians should be familiar with the changing epidemiology and clinical management of T. marneffei infection among non-HIV-infected patients.

Polyketides, toxins and pigments in Penicillium marneffei

Penicillium marneffei (synonym: Talaromyces marneffei) is the most important pathogenic thermally dimorphic fungus in China and Southeastern Asia. The HIV/AIDS pandemic, particularly in China and other Southeast Asian countries, has led to the emergence of P. marneffei infection as an important AIDS-defining condition. Recently, we published the genome sequence of P. marneffei. In the P. marneffei genome, 23 polyketide synthase genes and two polyketide synthase-non-ribosomal peptide synthase hybrid genes were identified. This number is much higher than those of Coccidioides immitis and Histoplasma capsulatum, important pathogenic thermally dimorphic fungi in the Western world. Phylogenetically, these polyketide synthase genes were distributed evenly with their counterparts found in Aspergillus species and other fungi, suggesting that polyketide synthases in P. marneffei did not diverge from lineage-specific gene duplication through a recent expansion. Gene knockdown experiments and ultra-high performance liquid chromatography-photodiode array detector/electrospray ionization-quadruple time of flight-mass spectrometry analysis confirmed that at least four of the polyketide synthase genes were involved in the biosynthesis of various pigments in P. marneffei, including melanin, mitorubrinic acid, mitorubrinol, monascorubrin, rubropunctatin, citrinin and ankaflavin, some of which were mycotoxins and virulence factors of the fungus.

High variability of mitochondrial gene order among fungi

From their origin as an early alpha proteobacterial endosymbiont to their current state as cellular organelles, large-scale genomic reorganization has taken place in the mitochondria of all main eukaryotic lineages. So far, most studies have focused on plant and animal mitochondrial (mt) genomes (mtDNA), but fungi provide new opportunities to study highly differentiated mtDNAs. Here, we analyzed 38 complete fungal mt genomes to investigate the evolution of mtDNA gene order among fungi. In particular, we looked for evidence of nonhomologous intrachromosomal recombination and investigated the dynamics of gene rearrangements. We investigated the effect that introns, intronic open reading frames (ORFs), and repeats may have on gene order. Additionally, we asked whether the distribution of transfer RNAs (tRNAs) evolves independently to that of mt protein-coding genes. We found that fungal mt genomes display remarkable variation between and within the major fungal phyla in terms of gene order, genome size, composition of intergenic regions, and presence of repeats, introns, and associated ORFs. Our results support previous evidence for the presence of mt recombination in all fungal phyla, a process conspicuously lacking in most Metazoa. Overall, the patterns of rearrangements may be explained by the combined influences of recombination (i.e., most likely nonhomologous and intrachromosomal), accumulated repeats, especially at intergenic regions, and to a lesser extent, mobile element dynamics.

The biosynthetic pathway for a thousand-year-old natural food colorant and citrinin in Penicillium marneffei

Monascorubrin and its derivatives are polyketides used as natural colorants for a wide range of food for more than one thousand years. Since the biosynthetic pathway for this ancient chemical compound is unknown and genome sequence unavailable for any Monascus species, monascorubrin production has relied on extraction from fungal cultures of Monascus species. In vitro synthesis and genetic manipulation are not possible. Here we report the polyketide gene cluster and pathway for monascorubrin biosynthesis in Penicillium marneffei, a diffusible red pigment-producing, thermal dimorphic fungus, taking advantage of available genome sequence and faster growth rate than Monascus species. We also documented that the red pigment of P. marneffei is a mixture of more than 16 chemical compounds, which are amino acid conjugates of monascorubrin and rubropunctatin, and showed that this polyketide gene cluster and pathway are also responsible for biosynthesis of ankaflavin and citrinin, a mycotoxin with nephrotoxic activity in mammals. The present study on elucidation of the biosynthetic pathway of monascorubrin is a proof-of-the-concept study that serves as a cornerstone for future studies on monascorubrin biosynthesis pathway dissection in Monascus species.

Proteome profiling of the dimorphic fungus Penicillium marneffei extracellular proteins and identification of glyceraldehyde-3-phosphate dehydrogenase as an important adhesion factor for conidial attachment

Despite being the most important thermal dimorphic fungus causing systemic mycosis in Southeast Asia, the pathogenic mechanisms of Penicillium marneffei remain largely unknown. By comparing the extracellular proteomes of P. marneffei in mycelial and yeast phases, we identified 12 differentially expressed proteins among which glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and heat shock protein 60 (HSP60) were found to be upregulated in mycelial and yeast phases respectively. Based on previous findings in other pathogens, we hypothesized that these two extracellular proteins may be involved in adherence during P. marneffei-host interaction. Using inhibition assays with recombinant GAPDH (rGAPDH) proteins and anti-rGAPDH sera, we demonstrated that adhesion of P. marneffei conidia to fibronectin and laminin was inhibited by rGAPDH or rabbit anti-rGAPDH serum in a dose-dependent manner. Similarly, a dose-dependent inhibition of conidial adherence to A549 pneumocytes by rGAPDH or rabbit anti-rGAPDH serum was observed, suggesting that P. marneffei GAPDH can mediate binding of conidia to human extracellular matrix proteins and pneumocytes. However, HSP60 did not exhibit similar inhibition on conidia adherence, and neither GAPDH norHSP60 exhibited inhibition on adherence to J774 or THP-1 macrophage cell lines. This report demonstrates GAPDH as an adherence factor in P. marneffei by mediating conidia adherence to host bronchoalveolar epithelium during the early establishment phase of infection.

Characterization and phylogenetic analysis of the mitochondrial genome of Glarea lozoyensis indicates high diversity within the order Helotiales

BACKGROUND

Glarea lozoyensis is a filamentous fungus used for the industrial production of non-ribosomal peptide pneumocandin B0. In the scope of a whole genome sequencing the complete mitochondrial genome of the fungus has been assembled and annotated. It is the first one of the large polyphyletic Helotiaceae family. A phylogenetic analysis was performed based on conserved proteins of the oxidative phosphorylation system in mitochondrial genomes.

RESULTS

The total size of the mitochondrial genome is 45,038 bp. It contains the expected 14 genes coding for proteins related to oxidative phosphorylation,two rRNA genes, six hypothetical proteins, three intronic genes of which two are homing endonucleases and a ribosomal protein rps3. Additionally there is a set of 33 tRNA genes. All genes are located on the same strand. Phylogenetic analyses based on concatenated mitochondrial protein sequences confirmed that G. lozoyensis belongs to the order of Helotiales and that it is most closely related to Phialocephala subalpina. However, a comparison with the three other mitochondrial genomes known from Helotialean species revealed remarkable differences in size, gene content and sequence. Moreover, it was found that the gene order found in P. subalpina and Sclerotinia sclerotiorum is not conserved in G. lozoyensis.

CONCLUSION

The arrangement of genes and other differences found between the mitochondrial genome of G. lozoyensis and those of other Helotiales indicates a broad genetic diversity within this large order. Further mitochondrial genomes are required in order to determine whether there is a continuous transition between the different forms of mitochondrial genomes or G. lozoyensis belongs to a distinct subgroup within Helotiales.

Identification of microRNA-like RNAs in mycelial and yeast phases of the thermal dimorphic fungus Penicillium marneffei

Background

Penicillium marneffei is the most important thermal dimorphic fungus causing systemic mycosis in China and Southeast Asia. While miRNAs are increasingly recognized for their roles in post-transcriptional regulation of gene expression in animals and plants, miRNAs in fungi were less well studied and their potential roles in fungal dimorphism were largely unknown. Based on P. marneffei genome sequence, we hypothesize that miRNA-like RNAs (milRNAs) may be expressed in the dimorphic fungus.

Methodology/Principal Findings

We attempted to identify milRNAs in P. marneffei in both mycelial and yeast phase using high-throughput sequencing technology. Small RNAs were more abundantly expressed in mycelial than yeast phase. Sequence analysis revealed 24 potential milRNA candidates, including 17 candidates in mycelial and seven in yeast phase. Two genes, dcl-1 and dcl-2, encoding putative Dicer-like proteins and the gene, qde-2, encoding Argonaute-like protein, were identified in P. marneffei. Phylogenetic analysis showed that dcl-2 of P. marneffei was more closely related to the homologues in other thermal dimorphic pathogenic fungi than to Penicillium chrysogenum and Aspergillus spp., suggesting the co-evolution of dcl-2 among the thermal dimorphic fungi. Moreover, dcl-2 demonstrated higher mRNA expression levels in mycelial than yeast phase by 7 folds (P<0.001). Northern blot analysis confirmed the expression of two milRNAs, PM-milR-M1 and PM-milR-M2, only in mycelial phase. Using dcl-1^KO, dcl-2^KO, dcl^DKO and qde-2^KO deletion mutants, we showed that the biogenesis of both milRNAs were dependent on dcl-2 but not dcl-1 or qde-2. The mRNA expression levels of three predicted targets of PM-milR-M1 were upregulated in knockdown strain PM-milR-M1^KD, supporting regulatory function of milRNAs.

Conclusions/Significance

Our findings provided the first evidence for differential expression of milRNAs in different growth phases of thermal dimorphic fungi and shed light on the evolution of fungal proteins involved in milRNA biogenesis and possible role of post-transcriptional control in governing thermal dimorphism.

Penicillium marneffei is the most important thermal dimorphic pathogenic fungus in Southeast Asia. Despite findings on diverse genes and mechanisms involved in dimorphic switching, the key to signally pathways governing the switch is still unknown. Since miRNAs are important regulatory molecules in eukaryotes, we attempt to define if miRNAs are expressed in different growth phases of P. marneffei. Using high-throughput sequencing, we identified 24 potential milRNA candidates in P. marneffei, which were more abundantly expressed in mycelial than yeast phase. Two genes, dcl-1 and dcl-2, encoding Dicer-like proteins and the gene, qde-2, encoding Argonaute-like protein, were also identified. Phylogenetic analysis showed that dcl-2 of P. marneffei was more closely related to the homologues in other thermal dimorphic pathogenic fungi than to Penicillium chrysogenum and Aspergillus spp.. dcl-2 demonstrated higher mRNA levels in mycelial than yeast phase. Northern blot analysis confirmed expression of two milRNAs, PM-milR-M1 and PM-milR-M2, only in mycelial phase, whose expression was dependent on dcl-2 but not dcl-1 or qde-2. The mRNA levels of three predicted targets of PM-milR-M1 were upregulated in knockdown strain PM-milR-M1^KD, supporting its regulatory function. This study represents the first discovery of milRNAs in thermal dimorphic fungi, with differential expression in different growth phases.

Romance of the three domains: how cladistics transformed the classification of cellular organisms

Cladistics is a biological philosophy that uses genealogical relationship among species and an inferred sequence of divergence as the basis of classification. This review critically surveys the chronological development of biological classification from Aristotle through our postgenomic era with a central focus on cladistics. In 1957, Julian Huxley coined cladogenesis to denote splitting from subspeciation. In 1960, the English translation of Willi Hennig’s 1950 work, Systematic Phylogenetics, was published, which received strong opposition from pheneticists, such as numerical taxonomists Peter Sneath and Robert Sokal, and evolutionary taxonomist, Ernst Mayr, and sparked acrimonious debates in 1960–1980. In 1977–1990, Carl Woese pioneered in using small subunit rRNA gene sequences to delimitate the three domains of cellular life and established major prokaryotic phyla. Cladistics has since dominated taxonomy. Despite being compatible with modern microbiological observations, i.e. organisms with unusual phenotypes, restricted expression of characteristics and occasionally being uncultivable, increasing recognition of pervasiveness and abundance of horizontal gene transfer has challenged relevance and validity of cladistics. The mosaic nature of eukaryotic and prokaryotic genomes was also gradually discovered. In the mid-2000s, high-throughput and whole-genome sequencing became routine and complex geneologies of organisms have led to the proposal of a reticulated web of life. While genomics only indirectly leads to understanding of functional adaptations to ecological niches, computational modeling of entire organisms is underway and the gap between genomics and phenetics may soon be bridged. Controversies are not expected to settle as taxonomic classifications shall remain subjective to serve the human scientist, not the classified.

Sequencing of mitochondrial genomes of nine Aspergillus and Penicillium species identifies mobile introns and accessory genes as main sources of genome size variability

BACKGROUND

The genera Aspergillus and Penicillium include some of the most beneficial as well as the most harmful fungal species such as the penicillin-producer Penicillium chrysogenum and the human pathogen Aspergillus fumigatus, respectively. Their mitochondrial genomic sequences may hold vital clues into the mechanisms of their evolution, population genetics, and biology, yet only a handful of these genomes have been fully sequenced and annotated.

RESULTS

Here we report the complete sequence and annotation of the mitochondrial genomes of six Aspergillus and three Penicillium species: A. fumigatus, A. clavatus, A. oryzae, A. flavus, Neosartorya fischeri (A. fischerianus), A. terreus, P. chrysogenum, P. marneffei, and Talaromyces stipitatus (P. stipitatum). The accompanying comparative analysis of these and related publicly available mitochondrial genomes reveals wide variation in size (25-36 Kb) among these closely related fungi. The sources of genome expansion include group I introns and accessory genes encoding putative homing endonucleases, DNA and RNA polymerases (presumed to be of plasmid origin) and hypothetical proteins. The two smallest sequenced genomes (A. terreus and P. chrysogenum) do not contain introns in protein-coding genes, whereas the largest genome (T. stipitatus), contains a total of eleven introns. All of the sequenced genomes have a group I intron in the large ribosomal subunit RNA gene, suggesting that this intron is fixed in these species. Subsequent analysis of several A. fumigatus strains showed low intraspecies variation. This study also includes a phylogenetic analysis based on 14 concatenated core mitochondrial proteins. The phylogenetic tree has a different topology from published multilocus trees, highlighting the challenges still facing the Aspergillus systematics.

CONCLUSIONS

The study expands the genomic resources available to fungal biologists by providing mitochondrial genomes with consistent annotations for future genetic, evolutionary and population studies. Despite the conservation of the core genes, the mitochondrial genomes of Aspergillus and Penicillium species examined here exhibit significant amount of interspecies variation. Most of this variation can be attributed to accessory genes and mobile introns, presumably acquired by horizontal gene transfer of mitochondrial plasmids and intron homing.

First discovery of two polyketide synthase genes for mitorubrinic acid and mitorubrinol yellow pigment biosynthesis and implications in virulence of Penicillium marneffei

Background

The genome of P. marneffei, the most important thermal dimorphic fungus causing respiratory, skin and systemic mycosis in China and Southeast Asia, possesses 23 polyketide synthase (PKS) genes and 2 polyketide synthase nonribosomal peptide synthase hybrid (PKS-NRPS) genes, which is of high diversity compared to other thermal dimorphic pathogenic fungi. We hypothesized that the yellow pigment in the mold form of P. marneffei could also be synthesized by one or more PKS genes.

Methodology/Principal Findings

All 23 PKS and 2 PKS-NRPS genes of P. marneffei were systematically knocked down. A loss of the yellow pigment was observed in the mold form of the pks11 knockdown, pks12 knockdown and pks11pks12 double knockdown mutants. Sequence analysis showed that PKS11 and PKS12 are fungal non-reducing PKSs. Ultra high performance liquid chromatography-photodiode array detector/electrospray ionization-quadruple time of flight-mass spectrometry (MS) and MS/MS analysis of the culture filtrates of wild type P. marneffei and the pks11 knockdown, pks12 knockdown and pks11pks12 double knockdown mutants showed that the yellow pigment is composed of mitorubrinic acid and mitorubrinol. The survival of mice challenged with the pks11 knockdown, pks12 knockdown and pks11pks12 double knockdown mutants was significantly better than those challenged with wild type P. marneffei (P<0.05). There was also statistically significant decrease in survival of pks11 knockdown, pks12 knockdown and pks11pks12 double knockdown mutants compared to wild type P. marneffei in both J774 and THP1 macrophages (P<0.05).

Conclusions/Significance

The yellow pigment of the mold form of P. marneffei is composed of mitorubrinol and mitorubrinic acid. This represents the first discovery of PKS genes responsible for mitorubrinol and mitorubrinic acid biosynthesis. pks12 and pks11 are probably responsible for sequential use in the biosynthesis of mitorubrinol and mitorubrinic acid. Mitorubrinol and mitorubrinic acid are virulence factors of P. marneffei by improving its intracellular survival in macrophages.

Penicillium marneffei is the most important thermal dimorphic fungus causing respiratory, skin and systemic mycosis in China and Southeast Asia. Its genome possesses a large number of polyketide synthase (PKS) genes, which should be responsible for synthesis of secondary metabolites such as pigments, antibiotics and mycotoxins. Using state-of-the-art gene knockdown and ultra high performance liquid chromatography-photodiode array detector/electrospray ionization-quadruple time of flight-mass spectrometry technologies, we discovered that the yellow pigment of P. marneffei was composed of mitorubrinol and mitorubrinic acid and was synthesized by two PKS genes, named pks12 and pks11. This represents the first discovery of PKS genes responsible for mitorubrinol and mitorubrinic acid biosynthesis, in which pks12 and pks11 are probably responsible for sequential use in the biosynthesis of mitorubrinol and mitorubrinic acid. Using a mouse model and human and mouse macrophage cell line models for P. marneffei infection, we also discovered that mitorubrinol and mitorubrinic acid are virulence factors of P. marneffei by improving its intracellular survival in macrophages.

Phylogenetic analysis of the complete mitochondrial genome of Madurella mycetomatis confirms its taxonomic position within the order Sordariales

BACKGROUND

Madurella mycetomatis is the most common cause of human eumycetoma. The genus Madurella has been characterized by overall sterility on mycological media. Due to this sterility and the absence of other reliable morphological and ultrastructural characters, the taxonomic classification of Madurella has long been a challenge. Mitochondria are of monophyletic origin and mitochondrial genomes have been proven to be useful in phylogenetic analyses.

RESULTS

The first complete mitochondrial DNA genome of a mycetoma-causative agent was sequenced using 454 sequencing. The mitochondrial genome of M. mycetomatis is a circular DNA molecule with a size of 45,590 bp, encoding for the small and the large subunit rRNAs, 27 tRNAs, 11 genes encoding subunits of respiratory chain complexes, 2 ATP synthase subunits, 5 hypothetical proteins, 6 intronic proteins including the ribosomal protein rps3. In phylogenetic analyses using amino acid sequences of the proteins involved in respiratory chain complexes and the 2 ATP synthases it appeared that M. mycetomatis clustered together with members of the order Sordariales and that it was most closely related to Chaetomium thermophilum. Analyses of the gene order showed that within the order Sordariales a similar gene order is found. Furthermore also the tRNA order seemed mostly conserved.

CONCLUSION

Phylogenetic analyses of fungal mitochondrial genomes confirmed that M. mycetomatis belongs to the order of Sordariales and that it was most closely related to Chaetomium thermophilum, with which it also shared a comparable gene and tRNA order.

Comprehensive transcription analysis of human pathogenic fungus Penicillium marneffei in mycelial and yeast cells

Penicillium marneffei is a dimorphic fungus that responds to changes in temperature. We performed a comparative analysis of gene expression in mycelial- and yeast-phase P. marneffei cells using high-throughput DNA microarrays. A total of 1,884 differentially expressed genes with annotations in the gene ontology (GO) database were identified between P. marneffei mycelial and yeast cells. These differentially expressed genes mainly belong to 18 categories in the organism's ontology, including reproduction, immunity, metabolism, signaling, etc. Bioinformatics suggests that these differentially expressed genes may help explain the resistance to adverse environments and the virulence of P. marneffei. Nine genes from the results of the DNA microarray experiment were selected for further analysis with real-time quantitative PCR to validate the differential expression. Our data provide a global description of the transcriptional response accompanying adaptation to high temperature in yeast.

Draft genome sequence of Penicillium marneffei strain PM1

Penicillium marneffei is the most important thermal dimorphic, pathogenic fungus endemic in China and Southeast Asia and is particularly important in HIV-positive patients. We report the 28,887,485-bp draft genome sequence of P. marneffei, which contains its complete mitochondrial genome, sexual cycle genes, a high diversity of Mp1p homologues, and polyketide synthase genes.

Complete mitochondrial genome of compactin-producing fungus Penicillium solitum and comparative analysis of Trichocomaceae mitochondrial genomes

We determined the complete mitochondrial genome sequence of the compactin-producing fungus Penicillium solitum strain 20-01. The 28 601-base pair circular-mapping DNA molecule encodes a characteristic set of mitochondrial proteins and RNA genes and is intron-free. All 46 protein- and RNA-encoding genes are located on one strand and apparently transcribed in one direction. Comparative analysis of this mtDNA and previously sequenced but unannotated mitochondrial genomes of several medically and industrially important species of the Aspergillus/Penicillium group revealed their extensive similarity in terms of size, gene content and sequence, which is also reflected in the almost perfect conservation of mitochondrial gene order in Penicillium and Aspergillus. Phylogenetic analysis based on concatenated mitochondrial protein sequences confirmed the monophyletic origin of Eurotiomycetes.

Approaches to Fungal Genome Annotation

Fungal genome annotation is the starting point for analysis of genome content. This generally involves the application of diverse methods to identify features on a genome assembly such as protein-coding and non-coding genes, repeats and transposable elements, and pseudogenes. Here we describe tools and methods leveraged for eukaryotic genome annotation with a focus on the annotation of fungal nuclear and mitochondrial genomes. We highlight the application of the latest technologies and tools to improve the quality of predicted gene sets. The Broad Institute eukaryotic genome annotation pipeline is described as one example of how such methods and tools are integrated into a sequencing center's production genome annotation environment.

Cordyceps militaris (Hypocreales: Cordycipitaceae): transcriptional analysis and molecular characterization of cox1 and group I intron with putative LAGLIDADG endonuclease

The full-length cytochrome c oxidase subunit I gene (cox1) containing a group I intron was isolated from an important medical fungus Cordyceps militaris (Cordycipitaceae). The open reading frame (ORF) of 1,593 nucleotides encoded a predicted protein COX1 of 530 amino acids. The group I intron encoded a putative homing endonuclease (HE) with two LAGLIDADG motifs. RT-PCR and Northern analysis showed a mature transcript of spliced cox1. Both 5'exon-intron and intron-3'exon junctions were also found by RT-PCR, suggesting the possible presence of unspliced cox1 RNA in total RNA. Sequence comparison by BLASTn showed that the coding region of cox1 (CRcox1) of C. militaris had significant similarities to those of related species (such as Cordyceps bassiana and C. brongniartii), while the intron had no significant homologous sequences of Cordycipitaceae fungi in NCBI database. The phylogenetic tree based on the CRcox1 confirmed the present taxonomic status of related species, but the cox1 introns were phylogenetically distinct. Compared to C. bassiana and C. brongniartii, the cox1 intron of C. militaris exhibited specific splicing site and different intronic ORF. The analysis of the folding RNA structures of the known cox1 introns from Cordyceps species showed different base pairs and conserved regions (P1-P10) in their structures. The present results provide useful information on the studies of cox1 intron splicing and Cordyceps evolution.

Classical and alternative components of the mitochondrial respiratory chain in pathogenic fungi as potential therapeutic targets

The frequency of opportunistic fungal infection has increased drastically, mainly in patients who are immunocompromised due to organ transplant, leukemia or HIV infection. In spite of this, only a few classes of drugs with a limited array of targets, are available for antifungal therapy. Therefore, more specific and less toxic drugs with new molecular targets is desirable for the treatment of fungal infections. In this context, searching for differences between mitochondrial mammalian hosts and fungi in the classical and alternative components of the mitochondrial respiratory chain may provide new potential therapeutic targets for this purpose.

Characterization of the cytochrome b (cyt b) gene from Monilinia species causing brown rot of stone and pome fruit and its significance in the development of QoI resistance

BACKGROUND

Quinone outside inhibitor (QoI) resistance as a consequence of point mutations in the cytochrome b (cyt b) gene has been reported in numerous plant pathogenic fungi. To examine the potential for QoI resistance development in those Monilinia species causing brown rot of stone and pome fruits [Monilinia fructicola (G Winter) Honey, M. laxa (Aderhold & Ruhland) Honey and M. fructigena (Aderhold & Ruhland) Honey], an examination was made of the sequence and exon/intron structure of their cyt b genes for the presence of any point mutations and/or introns commonly associated with resistance to QoIs in fungal plant pathogens.

RESULTS

None of the point mutations typically linked to QoI resistance was present in any of the Monilinia isolates examined. Furthermore, the cyt b genes from M. fructicola and M. laxa, but not M. fructigena, possessed a group-I-like intron directly after codon 143. Based on the results obtained, a simple PCR assay using a single primer pair was developed, allowing discrimination between the three Monilinia species without the need for culturing.

CONCLUSIONS

Results suggest that resistance to QoI fungicides based on the G143A mutation is not likely to occur in M. fructicola or M. laxa. Conversely, M. fructigena may be at higher risk for developing QoI resistance owing to the absence of a G143-associated intron.

Recent dermatophyte divergence revealed by comparative and phylogenetic analysis of mitochondrial genomes

Background

Dermatophytes are fungi that cause superficial infections of the skin, hair, and nails. They are the most common agents of fungal infections worldwide. Dermatophytic fungi constitute three genera, Trichophyton, Epidermophyton, and Microsporum, and the evolutionary relationships between these genera are epidemiologically important. Mitochondria are considered to be of monophyletic origin and mitochondrial sequences offer many advantages for phylogenetic studies. However, only one complete dermatophyte mitochondrial genome (E. floccosum) has previously been determined.

Results

The complete mitochondrial DNA sequences of five dermatophyte species, T. rubrum (26,985 bp), T. mentagrophytes (24,297 bp), T. ajelloi (28,530 bp), M. canis (23,943 bp) and M. nanum (24,105 bp) were determined. These were compared to the E. floccosum sequence. Mitochondrial genomes of all 6 species were found to harbor the same set of genes arranged identical order indicating that these dermatophytes are closely related. Genome size differences were largely due to variable lengths of non-coding intergenic regions and the presence/absence of introns. Phylogenetic analyses based on complete mitochondrial genomes reveals that the divergence of the dermatophyte clade was later than of other groups of pathogenic fungi.

Conclusion

This is the first systematic comparative genomic study on dermatophytes, a highly conserved and recently-diverged lineage of ascomycota fungi. The data reported here provide a basis for further exploration of interrelationships between dermatophytes and will contribute to the study of mitochondrial evolution in higher fungi.

The mitochondrial view of Blastocladiella emersonii

The mitochondrial genome of the chytrid Blastocladiella emersonii was sequenced and annotated, revealing the complete set of oxidative phosphorylation genes and tRNAs/rRNAs necessary for the translation process. Phylogenetic reconstructions reinforce the proposal of the new phylum Blastocladiomycota. Evidences of gene duplication due to inserted elements suggest shared susceptibility to gene invasion/exchange between chytrids and zygomycetes. The gene content of B. emersonii is very similar to Allomyces macrogynus but the content of intronic and changeable elements is much lower suggesting a stronger resistance to this kind of exchange. In addition, a total of 401 potential nuclear transcripts encoding mitochondrial proteins were obtained after B. emersonii EST database scanning using Saccharomyces cerevisiae, Homo sapiens and Arabidopsis thaliana data as probes and TargetP tool to find N-terminal mitochondrial signal in translated sequences.

Complete sequence and analysis of the mitochondrial genome of Hemiselmis andersenii CCMP644 (Cryptophyceae)

BACKGROUND

Cryptophytes are an enigmatic group of unicellular eukaryotes with plastids derived by secondary (i.e., eukaryote-eukaryote) endosymbiosis. Cryptophytes are unusual in that they possess four genomes-a host cell-derived nuclear and mitochondrial genome and an endosymbiont-derived plastid and 'nucleomorph' genome. The evolutionary origins of the host and endosymbiont components of cryptophyte algae are at present poorly understood. Thus far, a single complete mitochondrial genome sequence has been determined for the cryptophyte Rhodomonas salina. Here, the second complete mitochondrial genome of the cryptophyte alga Hemiselmis andersenii CCMP644 is presented.

RESULTS

The H. andersenii mtDNA is 60,553 bp in size and encodes 30 structural RNAs and 36 protein-coding genes, all located on the same strand. A prominent feature of the genome is the presence of a approximately 20 Kbp long intergenic region comprised of numerous tandem and dispersed repeat units of between 22-336 bp. Adjacent to these repeats are 27 copies of palindromic sequences predicted to form stable DNA stem-loop structures. One such stem-loop is located near a GC-rich and GC-poor region and may have a regulatory function in replication or transcription. The H. andersenii mtDNA shares a number of features in common with the genome of the cryptophyte Rhodomonas salina, including general architecture, gene content, and the presence of a large repeat region. However, the H. andersenii mtDNA is devoid of inverted repeats and introns, which are present in R. salina. Comparative analyses of the suite of tRNAs encoded in the two genomes reveal that the H. andersenii mtDNA has lost or converted its original trnK(uuu) gene and possesses a trnS-derived 'trnK(uuu)', which appears unable to produce a functional tRNA. Mitochondrial protein coding gene phylogenies strongly support a variety of previously established eukaryotic groups, but fail to resolve the relationships among higher-order eukaryotic lineages.

CONCLUSION

Comparison of the H. andersenii and R. salina mitochondrial genomes reveals a number of cryptophyte-specific genomic features, most notably the presence of a large repeat-rich intergenic region. However, unlike R. salina, the H. andersenii mtDNA does not possess introns and lacks a Lys-tRNA, which is presumably imported from the cytosol.

Agrobacterium tumefaciens-mediated transformation as a tool for insertional mutagenesis in the fungus Penicillium marneffei

Penicillium marneffei is an opportunistic fungal pathogen of humans, causing respiratory, skin, and systemic mycosis in south-east Asia. Here we describe the transformation of P. marneffei with Agrobacterium tumefaciens, and the optimization of the transformation procedure. Transformations in different combinations between A. tumefaciens stains (LBA4404 and EHA105) and binary vectors (pCB309A, pBI129A, and pCaMBIA1312A) showed that EHA105/pBI129A were the most efficient partners. Southern blot analysis suggested that 87.5% of transformants obtained with this protocol displayed single hybridization bands, indicating a single insert of T-DNA in each of the transformants. Unique hybridization patterns, along with thermal asymmetric interlaced PCR (TAIL-PCR) analysis of T-DNA insertion sites, suggested that A. tumefaciens-mediated transformation may be a powerful tool for insertional mutagenesis in P. marneffei. Several mutants with altered phenotypes were obtained during the construction of the mutant library, indicating the usefulness of the approach for functional genetic analysis in this important fungal pathogen.

Dothideomycete plant interactions illuminated by genome sequencing and EST analysis of the wheat pathogen Stagonospora nodorum

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

The mitochondrial genome from the thermal dimorphic fungus Paracoccidioides brasiliensis

We present here the sequence of the mitochondrial DNA of the pathogenic thermodimorphic fungus Paracoccidioides brasiliensis, agent of an endemic disease in most South American countries. The sequenced genome has 71 334 bp and is organized as a circular molecule with two gaps of unknown size flanking the middle exon of the nad5 gene. We located genes coding for the three subunits of the ATP synthase (atp6, atp8 and atp9), the apocytochrome b (cob), three subunits of the cytochrome c oxidase enzyme complex (cox1, cox2 and cox3), seven subunits of the reduced nicotinamide adenine dinucleotide ubiquinone oxidoreductase (nad1, nad2, nad3, nad4, nad5, nad6 and nad4L) and the large (rnl) and small (rns) subunits of ribosomal RNA. Two maturases and a ribosomal protein (rms5) are located inside introns. Twenty-five tRNAs were identified with acceptors for all 20 amino acids. Seven polypurine/polypyrimidine tracts (140-240 bp) have been found in this genome. All genes are in the same orientation over the genome, while their order is closest to the mitochondrial genomes from Penicillium marneffei and Aspergillus nidulans.

MP1 homologue-based multilocus sequence system for typing the pathogenic fungus Penicillium marneffei: a novel approach using lineage-specific genes

A highly reproducible and discriminative typing system is essential for better understanding of the epidemiology of Penicillium marneffei, the most important thermal dimorphic fungus causing respiratory, skin, and systemic mycosis in Southeast Asia. The sequences of 11 housekeeping genes were identical among 10 strains of P. marneffei, but those of MP1 and its 13 homologues, a novel superfamily of mannoproteins in the subdivision Pezizomycotina of Ascomycetes, mostly species of Penicillium and Aspergillus, showed significant variations. Therefore, a multilocus sequence typing (MLST) system for P. marneffei was constructed using MP1 (549 bp) and the four of its homologues (MPLP4 [337 bp], MPLP7 [347 bp], MPLP10 [546 bp], and MPLP13 [422 bp]) that showed the greatest variations. Among the 2,201 bp of the five loci, 183 polymorphic sites were observed in 44 strains of P. marneffei. The median number of alleles at each locus was five (range, 5 [MPLP4, MPLP7, and MPLP13] to 15 [MPLP10]). Four of the five genes had nonsynonymous substitution/synonymous substitution (d(n)/d(s)) ratios of >1. A total of 35 different sequence types (STs) were assigned to the 44 P. marneffei isolates, with 28 of the 35 STs identified only once. The discriminatory power was 0.9884. MP1 and its homologues were better than housekeeping genes for MLST in P. marneffei. Due to their more rapid evolutionary rates, lineage-specific genes may be better candidates than housekeeping genes for sequence-based typing, especially in microbes that evolve slowly or have evolved recently.

Novel group I introns encoding a putative homing endonuclease in the mitochondrial cox1 gene of Scleractinian corals

Analyses of mitochondrial sequences revealed the existence of a group I intron in the cytochrome oxidase subunit 1 (cox1) gene in 13 of 41 genera (20 out of 73 species) of corals conventionally assigned to the suborder Faviina. With one exception, phylogenies of the coral cox1 gene and its intron were concordant, suggesting at most two insertions and many subsequent losses. The coral introns were inferred to encode a putative homing endonuclease with a LAGLI-DADG motif as reported for the cox1 group I intron in the sea anemone Metridium senile. However, the coral and sea anemone cox1 group I introns differed in several aspects, such as the intron insertion site and sequence length. The coral cox1 introns most closely resemble the mitochondrial cox1 group I introns of a sponge species, which also has the same insertion site. The coral introns are also more similar to the introns of several fungal species than to that of the sea anemone (although the insertion site differs in the fungi). This suggests either a horizontal transfer between a sponge and a coral or independent transfers from a similar fungal donor (perhaps one with an identical insertion site that has not yet been discovered). The common occurrence of this intron in corals strengthens the evidence for an elevated abundance of group I introns in the mitochondria of anthozoans.

Prospects for fungus identification using CO1 DNA barcodes, with Penicillium as a test case

DNA barcoding systems employ a short, standardized gene region to identify species. A 648-bp segment of mitochondrial cytochrome c oxidase 1 (CO1) is the core barcode region for animals, but its utility has not been tested in fungi. This study began with an examination of patterns of sequence divergences in this gene region for 38 fungal taxa with full CO1 sequences. Because these results suggested that CO1 could be effective in species recognition, we designed primers for a 545-bp fragment of CO1 and generated sequences for multiple strains from 58 species of Penicillium subgenus Penicillium and 12 allied species. Despite the frequent literature reports of introns in fungal mitochondrial genomes, we detected introns in only 2 of 370 Penicillium strains. Representatives from 38 of 58 species formed cohesive assemblages with distinct CO1 sequences, and all cases of sequence sharing involved known species complexes. CO1 sequence divergences averaged 0.06% within species, less than for internal transcribed spacer nrDNA or beta-tubulin sequences (BenA). CO1 divergences between species averaged 5.6%, comparable to internal transcribed spacer, but less than values for BenA (14.4%). Although the latter gene delivered higher taxonomic resolution, the amplification and alignment of CO1 was simpler. The development of a barcoding system for fungi that shares a common gene target with other kingdoms would be a significant advance.

Nucleotide sequence-based analysis for determining the molecular epidemiology of Penicillium marneffei

The dimorphic fungus, Penicillium marneffei, is an emerging opportunistic pathogen endemic in Southeast Asia, especially for those with impaired cellular immunity such as human immunodeficiency virus-infected persons. A discriminatory and reproducible method based on the analysis of nucleotide sequences would facilitate epidemiologic investigations of this fungus. Twenty-four clinical or environmental isolates of P. marneffei obtained from China, Thailand, and Vietnam were analyzed by nucleotide sequence analysis. A total of 3,803 bp, consisting of eight nuclear gene fragments (transcription factor [AbaA], catalase [CpeA]], homodomain transcription factor [StlA], isocitrate lyase [Icl1], polyaromatic amino acid biosynthesis [PAA], NADH-dependent glutamate synthase [NGS], lovastatin nonaketide synthase [LNS], a cell wall mannoprotein [MP1], and a gene fragment of the cytochrome oxidase subunit 1 gene [COX1] of the P. marneffei mitochondrial genome) were amplified by PCR and then sequenced. No polymorphic sites within the Cox1 gene fragment were observed. Likewise, no nucleotide sequence polymorphisms were observed for three gene fragments: StlA, AbaA, and NGS. Seven single-nucleotide polymorphisms were observed for three gene fragments, Icl1, CpeA, and PAA, providing only a low degree of discriminatory power (D = 0.747). In contrast, the gene fragment for an antigenic cell wall glycoprotein, MP1, a useful immunologic marker for infection, was observed to be highly polymorphic with 12 different MP1 types (D = 0.887). Single-nucleotide polymorphisms were observed at 21 different locations in the MP1 gene fragment. Indels of 3, 21, 24, and 42 bp were observed and were in frame for protein translation. The relatively high degree of MP1 polymorphisms suggests the sequence is rapidly evolving in order to evade host immune responses. After all polymorphic gene sequences were combined, a high degree of genetic variation was observed (D = 0.949) for a total of 16 different haploid sequence types with 11 genotypes represented by single isolates. Phylogenetic analysis detected clusters composed of isolates obtained only from China or Thailand, as well as clusters with a combination of isolates from these two countries, indicating some mixing or common descent. Identical sequences were observed for isolates passed in vitro for 8 weeks, suggesting good reproducibility. The low degree of nucleotide diversity in housekeeping and regulatory genes suggests the recent emergence and spread as a species or an evolutionary bottleneck. In summary, multilocus sequence typing demonstrated a high degree of discriminatory power and reproducibility and may provide a robust and reliable adjunct method for genotyping isolates of P. marneffei and facilitating interlaboratory comparisons.

Putative cross-kingdom horizontal gene transfer in sponge (Porifera) mitochondria

BACKGROUND

The mitochondrial genome of Metazoa is usually a compact molecule without introns. Exceptions to this rule have been reported only in corals and sea anemones (Cnidaria), in which group I introns have been discovered in the cox1 and nad5 genes. Here we show several lines of evidence demonstrating that introns can also be found in the mitochondria of sponges (Porifera).

RESULTS

A 2,349 bp fragment of the mitochondrial cox1 gene was sequenced from the sponge Tetilla sp. (Spirophorida). This fragment suggests the presence of a 1143 bp intron. Similar to all the cnidarian mitochondrial introns, the putative intron has group I intron characteristics. The intron is present in the cox1 gene and encodes a putative homing endonuclease. In order to establish the distribution of this intron in sponges, the cox1 gene was sequenced from several representatives of the demosponge diversity. The intron was found only in the sponge order Spirophorida. A phylogenetic analysis of the COI protein sequence and of the intron open reading frame suggests that the intron may have been transmitted horizontally from a fungus donor.

CONCLUSION

Little is known about sponge-associated fungi, although in the last few years the latter have been frequently isolated from sponges. We suggest that the horizontal gene transfer of a mitochondrial intron was facilitated by a symbiotic relationship between fungus and sponge. Ecological relationships are known to have implications at the genomic level. Here, an ecological relationship between sponge and fungus is suggested based on the genomic analysis.

A novel approach for screening immunogenic proteins inPenicillium marneffeiusing the ΔAFMP1ΔAFMP2deletion mutant ofAspergillus fumigatus

Using serum from guinea-pigs immunized with a DeltaAFMP1DeltaAFMP2 deletion mutant of Aspergillus fumigatus to screen a cDNA library of A. fumigatus, we cloned a novel immunogenic 57-kDa protein in A. fumigatus. We also cloned its 55-kDa homologue in Penicillium marneffei, which was possibly related to amino acid biosynthesis and metabolism, with homologues present only in the subphylum Pezizomycotina of Ascomycota. The recombinant 55-kDa protein of P. marneffei reacted strongly with guinea-pig serum immunized with P. marneffei and with the sera of patients with P. marneffei infection. A similar approach could be applied to immunogenic protein screening in other microorganisms for serological diagnosis, epidemiological studies and the study of vaccines.

Genomic and experimental evidence for a potential sexual cycle in the pathogenic thermal dimorphic fungus Penicillium marneffei

All meiotic genes (except HOP1) and genes encoding putative pheromone processing enzymes, pheromone receptors and pheromone response pathways proteins in Aspergillus fumigatus and Aspergillus nidulans and a putative MAT-1 alpha box mating-type gene were present in the Penicillium marneffei genome. A putative MAT-2 high-mobility group mating-type gene was amplified from a MAT-1 alpha box mating-type gene-negative P. marneffei strain. Among 37 P. marneffei patient strains, MAT-1 alpha box and MAT-2 high-mobility group mating-type genes were present in 23 and 14 isolates, respectively. We speculate that P. marneffei can potentially be a heterothallic fungus that does not switch mating type.

Penicillium marneffei infection and recent advances in the epidemiology and molecular biology aspects

Penicillium marneffei infection is an important emerging public health problem, especially among patients infected with human immunodeficiency virus in the areas of endemicity in southeast Asia, India, and China. Within these regions, P. marneffei infection is regarded as an AIDS-defining illness, and the severity of the disease depends on the immunological status of the infected individual. Early diagnosis by serologic and molecular assay-based methods have been developed and are proving to be important in diagnosing infection. The occurrence of natural reservoirs and the molecular epidemiology of P. marneffei have been studied; however, the natural history and mode of transmission of the organism remain unclear. Soil exposure, especially during the rainy season, has been suggested to be a critical risk factor. Using a highly discriminatory molecular technique, multilocus microsatellite typing, to characterize this fungus, several isolates from bamboo rats and humans were shown to share identical multilocus genotypes. These data suggest either that transmission of P. marneffei may occur from rodents to humans or that rodents and humans are coinfected from common environmental sources. These putative natural cycles of P. marneffei infection need further investigation. Studies on the fungal genetics of P. marneffei have been focused on the characterization of genetic determinants that may play important roles in asexual development, mycelial-to-yeast phase transition, and the expression of antigenic determinants. Molecular studies have identified several genes involved in germination, hyphal development, conidiogenesis, and yeast cell polarity. A number of functionally important genes, such as the malate synthase- and catalase-peroxidase protein-encoding genes, have been identified as being upregulated in the yeast phase. Future investigations pertaining to the roles of these genes in host-fungus interactions may provide the key knowledge to understanding the pathogenicity of P. marneffei.

The complete DNA sequence of the mitochondrial genome of the dermatophyte fungus Epidermophyton floccosum

We report here the complete nucleotide sequence of the 30.9-kb mitochondrial genome of the dermatophyte fungus Epidermophyton floccosum. All genes are encoded on the same DNA strand and include seven subunits of the reduced nicotinamide adenine dinucleotide ubiquinone oxireductase (nad1, nad2, nad3, nad4, nad4L, nad5, and nad6), three subunits of cytochrome oxidase (cox1, cox2, and cox3), apocytochrome b (cob), three subunits of ATP synthase (atp6, atp8, and atp9), the small and large ribosomal RNAs (rns and rnl), and 25 tRNAs. A ribosomal protein gene (rps5) is present as an intronic ORF in the large ribosomal subunit. The genes coding for cob and cox1 carry one intron and nad5 carries two introns with ORFs. The mtDNA of E. floccosum has the same gene order as Trichophyton rubrum mtDNA, with the exception of some tRNA genes. Maximum likelihood phylogenetic analysis confirms T. rubrum as a close relative of E. floccosum. This is the first complete mitochondrial sequence of a species of the order Onygenales. This sequence is available under GenBank accession number AY916130.

AFMP2 encodes a novel immunogenic protein of the antigenic mannoprotein superfamily in Aspergillus fumigatus

We cloned the Aspergillus fumigatus mannoprotein 2 (AFMP2) gene, which encodes a novel immunogenic protein (Afmp2p) of the antigenic mannoprotein superfamily, in A. fumigatus. Sequence analysis revealed that Afmp2p has 510 amino acid residues, with a predicted molecular mass of 51.5 kDa. Afmp2p has a putative N-terminal signal peptide, a putative C-terminal glycosylphosphatidylinositol membrane attachment signal sequence, and an upstream GAA cleavage site commonly used for cytoplasmic membrane attachment and implicated in fungal cell wall assembly. Upstream of the GAA cleavage site, Afmp2p contains a 302-amino-acid serine- and threonine-rich region as a site for potential O-glycosylation. Within this serine- and threonine-rich region, 13 repeats of ETSTPCE(T)(n) were observed. Western blot analysis of Afmp2p in A. fumigatus fungal cell lysate and culture supernatant and immunogold staining and electron microscopy showed that Afmp2p is predominantly secreted into the culture supernatant, whereas only minimal amounts can be detected in the cell lysate and cell wall. Finally, it was observed that patients with aspergilloma and invasive aspergillosis due to A. fumigatus develop a specific antibody response against recombinant Afmp2p. The abundance of Afmp2p in secreted form, its minimal cross-reactivity with Afmp1p, and the presence of an antibody response against Afmp2p in patients with A. fumigatus infections suggest that Afmp2p is a good candidate for complementing Afmp1p in serodiagnosis of A. fumigatus infections.