Next Article in Journal
The Impact of the Nitric Oxide (NO)/Soluble Guanylyl Cyclase (sGC) Signaling Cascade on Kidney Health and Disease: A Preclinical Perspective
Next Article in Special Issue
Mathematical Models in the Description of Pregnane X Receptor (PXR)-Regulated Cytochrome P450 Enzyme Induction
Previous Article in Journal
ARID1B as a Potential Therapeutic Target for ARID1A-Mutant Ovarian Clear Cell Carcinoma
Previous Article in Special Issue
The Role of Genetic Polymorphisms in Chronic Pain Patients
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 19
Issue 6
10.3390/ijms19061711
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus

by
Mathula Lancelot Ngwenya
1,
Wanping Chen
2,
Albert Kotze Basson
1,
Jabulani Siyabonga Shandu
1,
Jae-Hyuk Yu
3,
David R. Nelson
4,* and
Khajamohiddin Syed
1,*
1
Department of Biochemistry and Microbiology, Faculty of Science and Agriculture, University of Zululand, KwaDlangezwa 3886, South Africa
2
College of Food Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
3
Department of Bacteriology, University of Wisconsin-Madison, 3155 MSB, 1550 Linden Drive, Madison, WI 53706, USA
4
Department of Microbiology, Immunology and Biochemistry, University of Tennessee Health Science Center, Memphis, TN 38163, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2018, 19(6), 1711; https://doi.org/10.3390/ijms19061711
Submission received: 14 March 2018 / Revised: 2 May 2018 / Accepted: 14 May 2018 / Published: 9 June 2018
(This article belongs to the Special Issue Cytochromes P450: Drug Metabolism and Bioactivation)
Browse Figures
Versions Notes

Abstract

:
While the Zygomycete fungus Conidiobolus coronatus primarily infects insects, it can be pathogenic to mammals as well, including humans. High variability in the treatment of this fungal infection with currently available drugs, including azole drugs is a very common phenomenon. Azoles bind to the cytochrome P450 monooxygenases (P450s/CYP) including CYP51, a sterol 14-α-demethylase, inhibiting the synthesis of cell membrane ergosterol and thus leading to the elimination of infecting fungi. Despite P450’s role as a drug target, to date, no information on C. coronatus P450s has been reported. Genome-wide data mining has revealed the presence of 142 P450s grouped into 12 families and 21 subfamilies in C. coronatus. Except for CYP51, the remaining 11 P450 families are new (CYP5854-CYP5864). Despite having a large number of P450s among entomopathogenic fungi, C. coronatus has the lowest number of P450 families, which suggests blooming P450s. Further analysis has revealed that 79% of the same family P450s is tandemly positioned, suggesting that P450 tandem duplication led to the blooming of P450s. The results of this study; i.e., unravelling the C. coronatus P450 content, will certainly help in designing experiments to understand P450s’ role in C. coronatus physiology, including a highly variable response to azole drugs with respect to P450s.

1. Introduction

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins that are ubiquitously distributed in organisms belonging to different biological kingdoms [1]. P450s are well known for their stereo- and regio-specific enzymatic reactions [2] and thus play a key role in organisms’ physiology, both in primary and secondary metabolism [3]. Because of their important catalytic activities, P450s’ role in drug metabolism and drug discovery, in the generation of commercial products including antibiotics, in bioremediation, and in the production of biofuels has been explored [3].
Analysis of fungal species genomes has revealed the presence of a large number of P450s, with few exceptions. It is now well known that the number of P450s and diversity of P450s in a fungal species is largely affected by the fungal species’ lifestyle or their ecological niche [4,5,6,7]. This indicates that P450s play a key role in fungal species’ physiology in terms of their primary and secondary metabolism. Fungal P450s’ role in fungal species’ physiology and their biotechnological application has been thoroughly reviewed [8,9,10,11]. Among all fungal P450 applications, the use of P450s as drug targets against fungal pathogens is well explored [12,13] because of their binding affinity to azole drugs, which are widely used for treating fungal infections [14]. Most of the azole drugs bind and inhibit CYP51, a sterol 14α-demethylase, thus stopping the synthesis of cell membrane ergosterol, which is an indispensable component of fungal membranes [12]. Ergosterol depletion affects not only the structure of the membrane, but also several of its functions, including its hormone-like role that stimulates the growth and proliferation of fungal cells that eventually lead to the death of fungi [14,15,16]. Fungal pathogens were found to develop drug resistance [17] to anti-fungal drugs such as azole drugs and other drugs belonging to the classes of allylamine, echinocandin, nucleoside analog, and polyene [18]. Understanding the fungal pathogens’ physiology at the molecular level will help to unravel novel drug targets and develop novel drugs.
Among fungal pathogens, Conidiobolus coronatus, an entomopathogen (pathogen of insects) developed the ability to infect animals (horses, sheep, and dogs), including humans [19]. C. coronatus lives in soil and decaying matter, especially dead leaves [20], in addition to the feces of amphibians, sheep, insects, and horses [19]. C. coronatus has a worldwide distribution, especially in the tropical rain forests of Africa. It has been found in the United Kingdom [21], on the eastern coast of the United States [22], in India [23,24], the western region of Africa, and recently in South Africa [25].
C. coronatus causes conidiobolomycosis, an entomophthoramycosis in humans, and is considered a rare and neglected disease [19]. The definitive route of infection has not been established, but it is believed that C. coronatus enters the human body via inhalation of spores or from minor trauma [26]. Conidiobolomycosis is a disease of the nasal submucosa and paranasal sinuses, which slowly spreads to the nasal skin, glabella, cheek, upper lip, and pharynx. Rarely, contiguous lymph nodes may be involved [24,27,28,29,30]. Infection of the nasal and paranasal mucocutaneous tissue is known as rhinoentomophthoramycosis. Recently, C. coronatus infection in the vagina has been reported [22], indicating that this fungus is capable of infecting other sites of the human body as well.
C. coronatus is capable of producing mycotoxins that are known to be very toxic to insects [31,32]. This parasitic fungus will attack normal or stressed insects, gaining access by penetrating the host’s inner skin, and the host will eventually die because of tissue demolition, nutrient depreciation, and the production of toxins [33]. This entomopathogenic fungus has the potential of becoming a biological insect control, owing to its ability to infect and kill a varied range of insects [34].
In order to understand C. coronatus’s physiology at the molecular level, genome sequencing of this fungus has been carried out [35]. However, the genome sequencing study was focused on analyzing the evolution of cell wall digestion enzymes and no information on C. coronatus P450s has been presented [35]. Interestingly, conidiobolomycosis treatment employing amphotericin B, azoles and/or potassium iodide was found to be highly variable [19]. In a case study, a patient infected with C. coronatus showed no response when treated with amphotericin B; azoles (ketoconazole, fluconazole and itraconazole) and potassium iodide individually, and a combination of azole drug therapies were needed to treat the patient [36]. Sometimes treatment includes both the surgical removal of infected tissue and the use of different drugs [37]. High variability in susceptibility to azole drugs by C. coronatus was widely observed [19], including C. coronatus species that are resistant to azole drugs [38].
Since P450s were found to play a key role in organisms’ physiology, including serving as azole drug targets [12,14,15,16,17], it is necessary to unravel the P450 content of this fungus to be able to design experiments to understand P450s role in C. coronatus physiology and the variability in susceptibility of this fungus to azole drugs with respect to P450s.

2. Results and Discussion

2.1. C. coronatus Has the Highest Number of P450s among Entomopathogenic Fungi

Genome data mining, identification, and annotation of P450s in C. coronatus revealed the presence of 142 P450s in its genome (Figure 1 and Table 1). C. coronatus was found to have the highest number of P450s compared to other entomopathogenic and animal (including human) pathogenic fungi (Table 2).

2.2. Large Number of New P450s Found in C. coronatus

C. coronatus P450s can be grouped into 12 P450 families and 21 P450 subfamilies (Table 1). Except for CYP51, the most conserved family in fungi [43], the remaining P450 families are new. It is quite interesting that C. coronatus, apart from CYP51, does not share a P450 family with other entomopathogenic fungi listed in Table 2. Interestingly, C. coronatus lacks CYP61, the P450 involved in fungal membrane ergosterol synthesis [44], compared with pathogenic fungi listed in Table 2. Among 142 P450s, one P450 (protein ID: 71373) is not assigned to a P450 family because of a short amino acid sequence and low sequence identity to other P450s. In addition, one P450 (protein ID: 2292) belonging to the CYP5854 family has not been assigned to a subfamily because of a short amino acid sequence and low identity to other subfamilies. P450s that are not full-length (<300 amino acids) were indicated with word “fragments” after their P450 subfamily (Figure 1 and Tables S1 and S2). The new P450 families identified in C. coronatus include CYP5854-CYP5864 (Table 1). Despite having a large number of P450s, C. coronatus has the lowest number of P450 families compared to entomopathogenic and animal (including human) pathogenic fungi (Table 2), suggesting P450 blooms (a few families with many genes) in C. coronatus. Phylogenetic analysis revealed that CYP5859 P450s show a close relationship with CYP5855 P450s and seem to have evolved from CYP5855 P450s (Figure 1).

2.3. P450 Signature Motifs EXXR and CXG Are Conserved in C. coronatus P450s

Among different motifs that are characteristic of P450s, two motifs, namely EXXR in the K-helix and FXXGXRXCXG (also known as CXG) in the heme-binding domain, are found to be conserved in P450s, with few exceptions [45,46,47,48]. These two motifs are largely explored in the identification of P450s across biological kingdoms [48]. Analysis of EXXR and CXG motifs in C. coronatus P450s revealed that all P450s have both signature motifs (Table S3). Only a few P450s that are named as fragments (as mentioned in Section 2.2) did not have one of these motifs (Table S3). Thus, these P450s named as fragments represent pseudo-P450s. Analysis revealed that C. coronatus P450s have highly conserved amino acids, such as glutamic acid (E) and cysteine (C), in the motifs EXXR and CXG, with the exception of CYP5857A4, which has glycine instead of glutamic acid in the EXXR motif (Table S3). Non-conservation of glutamic acid in EXXR motifs is not new and has been reported for other P450s as well [48]. Analysis of amino acid patterns for the P450 families, CYP5854-CYP5856, revealed family-specific amino acid patterns at the EXXR and CXG motifs (Figure 2). CYP5854 and CYP5855 have the same E-S-M-R amino acid pattern, whereas CYP5856 has an E-T/V-M-R amino acid pattern in the EXXR motif (Figure 2). Comparative analysis with P450s across biological kingdoms [48,49] revealed that the CYP5854 and CYP5855 families’ EXXR amino acid pattern matched that of the CYP94 family, whereas the CYP5856 family EXXR amino acid patterns partially matched those of the CYP53, CYP92, and CYP176 P450 families, indicating a possible evolutionary relation between the families. Contrary to the EXXR motifs, the amino acid patterns at the CXG motifs of the CYP5854-CYP5856 families did not match any P450s [48,49]. This strongly supports the hypothesis that each P450 family has its unique signature amino acid patterns at the EXXR and CXG motifs [48].

2.4. P450 Blooming in C. coronatus

Comparative analysis of P450s revealed that among 11 P450 families, five P450 families contributed 88% of P450s to the total number of P450s in C. coronatus (Figure 2 and Table 1). This indicates that the five P450 families, namely CYP5854, CYP5855, CYP5856, CYP5858, and CYP5857, are highly bloomed in C. coronatus. The highest number of member P450s were found in CYP5854 (48 P450s) followed by CYP5855 (36 P450s), CYP5856 (20 P450s), CYP5858 (11 P450s), and CYP5857 (10 P450s) (Figure 3). Blooming of certain P450 families in an organism indicates that these family members play a key role in an organism’s physiology and thus they are expanded [4,50]. In some Basidiomycetes, P450 families were found to be expanded in order to help those organisms to colonize on wood and wood-derived materials [4]. At present the reason or need for blooming of these five P450 families in C. coronatus is not clear. It is noteworthy that none of the P450 families is expanded in other entomopathogenic fungi listed in Table 2.

2.5. Extensive Tandem Duplications Led to P450 Blooming in C. coronatus

P450 blooming is only possible when certain P450s belong to the same P450 family duplicated in large numbers in an organism [50]. The duplicated P450s (paralogs) can be identified by observing their location in the genome, intron-exon organization, and percentage identity at protein level, as described elsewhere [4,6,13]. Tandem arrangement of P450s (P450s located one behind the other) belonging to the same family is a good indication of P450 duplications. Analysis of the C. coronatus genome revealed that 79% of P450s (112 P450s) belonging to the same P450 family are tandemly located (Figure 4 and Table S5).
As shown in Figure 4, the five P450 families bloomed in C. coronatus were found to have the highest number of tandemly duplicated P450s. This indicates that P450 tandem duplications led to the blooming of these five P450 families in C. coronatus. The CYP5855 family was found to have the highest number of tandemly duplicated P450s (33 P450s), followed by the CYP5844 family with 32 P450s, the CYP5856 family with 17 P450s, the CYP5858 family with 11 P450s, and the CYP5857 family with 10 P450s. The number of tandemly duplicated P450s in the respective families was shown in Figure 4. In addition to the five P450 families bloomed in C. coronatus, three P450 families were found to have P450s that are tandemly duplicated, namely CYP5861 (2 P450s), CYP5859 (5 P450s) and CYP5860 (2 P450s) (Figure 4 and Table S5). It is interesting to note that P450s that are tandemly located on different scaffolds (portion of the genome sequence reconstructed from end-sequenced whole-genome shotgun clones) belong to the same subfamily (Table S5). This strongly indicates that these P450s indeed evolved by tandem duplications.

2.6. C. coronatus P450s Have Same Gene Structure

Another characteristic of tandemly duplicated genes is that they have the same gene structure; i.e., the same number and size of introns and exons. Fungal P450s with the same gene structure have previously been reported and the authors suggested that this is a typical characteristic of gene duplications [4,6,13]. In this study, gene structure analysis was carried out for P450 families such as CYP5854, CYP5857, and CYP5860 (Figure 5 and Figure 6), since most of the members in these families are full-length P450s.
As shown in Figure 5, the size of the exons in CYP5854 P450s is conserved. The number of introns and exons in some P450s is also highly conserved. P450s that are short contain the same size exons as full-length P450s, indicating their origin from full-length P450s. This was clearly observed for some CYP5854 P450s such as CYP5854-fragment2, which originated from CYP5854A2, CYP5854A2-fragment4, which originated from CYP5854-fragment3 and CYP5854E1, which originated from CYP5854E2 (Figure 5).
The phenomenon of the same gene structure was also observed in CYP5857 and CYP5860 P450s (Figure 6). The origin of P450s from other P450s in the CYP5857 family was observed as well, where CYP5857A5 P450s seem to have originated from CYP5857A1 (Figure 6). The two members of the CYP5860 family have the same gene structure (Figure 6).

2.7. C. coronatus Has the Lowest P450 Diversity among Entomopathogenic Fungi

Another characteristic of tandem duplication of P450s is that organisms have the lowest P450 diversity percentage [7,49,51]. The P450 diversity percentage is a good indication of P450 family blooming. The lowest P450 diversity means that certain P450 families are bloomed in an organism. Analysis of the P450 diversity percentage revealed that C. coronatus has only 8% P450 diversity, despite having a large number of P450s in its genome (Figure 7).
Comparative analysis of the P450 diversity percentage with other entomopathogenic fungi revealed that C. coronatus has the lowest diversity (Figure 7). This further indicates that extensive P450 duplication led to blooming of certain P450 families and thus the lowest P450 diversity percentage.

2.8. Functional Prediction of C. coronatus P450s

Except for CYP51 P450, all P450s in C. coronatus are new, implying that functional analysis based on characterized homologs is not possible. However, as indicated earlier, CYP51 of C. coronatus can be predicted to be involved in ergosterol biosynthesis based on the characterized homolog P450 [43]. It is interesting to see the role of bloomed P450s in C. coronatus, as these families are unique to this fungus compared to other entomopathogenic fungi. It is interesting to note that another P450, CYP61, involved in fungal membrane ergosterol synthesis [44] is not present in C. coronatus.

3. Materials and Methods

3.1. Genome Data Mining for P450s

The genome sequence of C. coronatus NRRL28638 v1.0 has been published [35] and is available at the JGI MycoCosm database (Available online: https://genome.jgi.doe.gov/Conco1/Conco1.home.html), which was data-mined for P450s. The data mining of P450s was carried out as described elsewhere [4,6,52], with slight modifications. Briefly, the C. coronatus genome was mined for P450s using InterPro code “IPR001128”. The hit protein sequences were downloaded and subjected to the NCBI Batch Web CD-Search Tool [53]. Proteins that grouped under the P450 superfamily were selected for further analysis.

3.2. Annotation of P450s

The above selected P450 proteins were subjected to BLAST analysis against all named fungal P450s at the Cytochrome P450 Homepage [54] to identify a homolog P450. The name of the homolog P450 and percentage of identity to the homolog P450 were recorded (Table S1). Based on the percentage identity to homologous P450s, a family and subfamily were assigned to C. coronatus P450s. This was done by following the rules set by the International P450 Nomenclature Committee [55,56,57]; i.e., >40% identity for a family and >55% identity for a subfamily. C. coronatus P450s that had <40% and <55% identity to the named fungal P450s were assigned to new families and new subfamilies, respectively (Table S1). Percentage identity among C. coronatus P450s was also used to assign P450 families and subfamilies. Information on homolog P450s that are used in the naming of C. coronatus P450s was listed in Table S1 along with the C. coronatus P450s identified and annotated in this study. C. coronatus P450 protein sequences along with their amino acid length are presented in Table S2.

3.3. Phylogenetic Analysis of P450s

The phylogenetic tree of C. coronatus P450s was constructed as described previously [51], with slight modifications. Briefly, the protein sequences were aligned by MUSCLE embedded in MEGA 7 [58]. Then, the best-fit substitution model for alignment was determined by the IQ-TREE web server (Available online: http://iqtree.cibiv.univie.ac.at/) [59]. Finally, the tree was constructed in MEGA7 by the maximum likelihood method, along with the best-fit substitution model and 100 bootstrap replications [60]. iTOL was used to view and highlight the tree [61].

3.4. Analysis of EXXR and CXG Motifs

An analysis of EXXR and CXG motifs in C. coronatus P450s was carried out as described elsewhere [48]. Briefly, a multiple sequence alignment of P450s was carried out using Clustal Omega [62] and then a manual search was performed for the presence of EXXR and CXG motifs (Table S3). The WebLogo for EXXR and CXG motifs for the P450 families, CYP5854-CYP5856, was deduced as described elsewhere [48] following the WebLogo program [63]. WebLogos were deduced using the default parameter with the amino acid sequence type. The criteria for the selection of P450 families for EXXR and CXG amino acid pattern analysis are based on the number of P450s in a P450 family [49]. Three P450 families, CYP5854-CYP5856, have more than 15 P450s for analysis (Table S4) compared to other new P450 families, thus they qualify for EXXR and CXG motif analysis.

3.5. Identification of Tandemly Duplicated P450s

Tandem P450 gene duplications were assessed as described elsewhere [4,6]. Briefly, the physical location of P450 genes, such as scaffold number, start and end point of the P450 gene on the DNA strand, was recorded by scanning the C. coronatus genome using the P450 protein ID (Table S5). P450s that are tandemly arranged and belong to the same family are regarded as duplicated P450s. The percentage P450 duplications are calculated as the percentage of P450s duplicated in the total number of P450s.

3.6. Gene Structure Analysis

Gene structure analysis of P450s was carried out as described elsewhere [4,6,7,13]. Briefly, the C. coronatus genome was mined using P450 protein ID. P450 gene structure graphics downloaded from the C. coronatus genome database were aligned in such a way that exons with the same size would align together. The length of exons was noted as an indication of possible gene duplication, if P450s showed conservation in the size of exons. P450s whose gene structures contained the same sizes of exons or introns were presented in the figures.

3.7. P450 Diversity Percentage Analysis

P450 diversity percentage analysis was carried out as described elsewhere [7,49,51]. Briefly, the P450 diversity percentage in C. coronatus was measured as a percentage contribution of the number of P450 families in the total number of P450s.

3.8. Comparative Analysis of P450s

Entomopathogenic and animal (including human) pathogenic fungi P450s were retrieved from published articles [7,39,40,41,42] and used for comparative analysis with C. coronatus P450s.

4. Conclusions

This study is the first of its kind on an in silico analysis of P450s in C. coronatus. Except for CYP51 P450, a conserved P450 in fungi and a primary target of azole drugs, all P450s were found to be new in C. coronatus. Unprecedented blooming of novel P450s in C. coronatus compared to other fungal species was observed. It would be interesting to observe the role of bloomed P450s in C. coronatus physiology, considering this fungus has a broad host range and is capable of infecting humans and other animals, in addition to insects.

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/1422-0067/19/6/1711/s1.

Author Contributions

K.S. conceived and designed the experiments; M.L.N., W.C., A.K.B., J.S.S., J-H.Y., D.R.N., and K.S. performed the experiments, analysed the data and contributed analysis tools. M.L.N., W.C., A.K.B., J.S.S., J-H.Y., D.R.N., and K.S. were involved in writing the manuscript. All authors reviewed and approved the manuscript.

Acknowledgments

Mathula Lancelot Ngwenya, Albert Kotze Basson, Jabulani Siyabonga Shandu and Khajamohiddin Syed express sincere gratitude to the University of Zululand Research Committee for funding. Khajamohiddin Syed is grateful to the National Research Foundation, South Africa for a research grant (Grant No: 114159). The authors also want to extend their thanks to Barbara Bradley, Pretoria, South Africa for English language editing.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; the collection, analysis, or interpretation of data, the writing of the manuscript or the decision to publish the results.

References

  1. Nelson, D.R. Cytochrome P450 diversity in tree of life. Biochim. Biophys. Acta 2018, 1866, 141–154. [Google Scholar] [CrossRef] [PubMed]
  2. Bernhardt, R. Cytochromes P450 as versatile biocatalysts. J. Biotechnol. 2006, 124, 128–145. [Google Scholar] [CrossRef] [PubMed]
  3. Montellano, P.R.O. Cytochrome P450: Structure, Mechanism, and Biochemistry, 4th ed.; Springer: Basel, Switzerland, 2015; pp. 1–912. ISBN 978-3-319-12108-6. [Google Scholar]
  4. Syed, K.; Shale, K.; Pagadala, N.S.; Tuszynski, J. Systematic identification and evolutionary analysis of catalytically versatile cytochrome P450 monooxygenase families enriched in model basidiomycete fungi. PLoS ONE 2014, 9, e86683. [Google Scholar] [CrossRef] [PubMed]
  5. Kgosiemang, I.K.R.; Mashele, S.S.; Syed, K. Comparative genomics and evolutionary analysis of cytochrome P450 monooxygenases in fungal subphylum Saccharomycotina. J. Pure Appl. Microbiol. 2014, 8, 291–302. [Google Scholar]
  6. Qhanya, L.B.; Matowane, G.; Chen, W.; Sun, Y.; Letsimo, E.M.; Parvez, M.; Yu, J.H.; Mashele, S.S.; Syed, K. Genome-wide annotation and comparative analysis of cytochrome P450 monooxygenases in Basidiomycete biotrophic plant pathogens. PLoS ONE 2015, 10, e0142100. [Google Scholar] [CrossRef] [PubMed]
  7. Matowane, R.G.; Wieteska, L.; Bamal, H.D.; Kgosiemang, I.K.R.; Van Wyk, M.; Manume, N.A.; Abdalla, S.M.H.; Mashele, S.S.; Gront, D.; Syed, K. In silico analysis of cytochrome P450 monooxygenases in chronic granulomatous infectious fungus Sporothrix schenckii: Special focus on CYP51. Biochim. Biophys. Acta 2018, 1866, 166–177. [Google Scholar] [CrossRef] [PubMed]
  8. Sakaki, T. Practical application of cytochrome P450. Biol. Pharm. Bull. 2012, 35, 844–849. [Google Scholar] [CrossRef] [PubMed]
  9. Girhard, M.; Bakkes, P.J.; Mahmoud, O.; Urlacher, V.B. P450 Biotechnology. In Cytochrome P450: Structure, Mechanism, and Biochemistry, 4th ed.; Montellano, P.R.O., Ed.; Springer: Basel, Switzerland, 2015; pp. 451–520. ISBN 978-3-319-12108-6. [Google Scholar]
  10. Durairaj, P.; Hur, J.S.; Yun, H. Versatile biocatalysis of fungal cytochrome P450 monooxygenases. Microb. Cell Fact. 2016, 15, 125. [Google Scholar] [CrossRef] [PubMed]
  11. Shin, J.; Kim, J.E.; Lee, Y.W.; Son, H. Fungal cytochrome P450s and the P450 complement (CYPome) of Fusarium graminearum. Toxins 2018, 10, 112. [Google Scholar] [CrossRef] [PubMed]
  12. Kelly, S.L.; Kelly, D.E. Microbial cytochromes P450: Biodiversity and biotechnology. Where do cytochromes P450 come from, what do they do and what can they do for us? Phil. Trans. R. Soc. B 2013, 368. [Google Scholar] [CrossRef] [PubMed]
  13. Jawallapersand, P.; Mashele, S.S.; Kovačič, L.; Stojan, J.; Komel, R.; Pakala, S.B.; Kraševec, N.; Syed, K. Cytochrome P450 monooxygenase CYP53 family in fungi: Comparative structural and evolutionary analysis and its role as a common alternative anti-fungal drug target. PLoS ONE 2014, 9, e107209. [Google Scholar] [CrossRef] [PubMed]
  14. Sheehan, D.J.; Hitchcock, C.A.; Sibley, C.M. Current and emerging azole antifungal agents. Clin. Microbiol. Rev. 1999, 12, 40–79. [Google Scholar] [PubMed]
  15. Shapiro, R.S.; Robbins, N.; Cowen, L.E. Regulatory circuitry governing fungal development, drug resistance, and disease. Microbiol. Mol. Biol. Rev. 2011, 75, 213–267. [Google Scholar] [CrossRef] [PubMed]
  16. Hargrove, T.Y.; Wawrzak, Z.; Lamb, D.C.; Guengerich, F.P.; Lepesheva, G.I. Structure-fuctional characterization of cytochrome P450 sterol 14α-demethylase (CYP51B) from Aspergillus fumigatus and molecular basis for the development of antifungal drugs. J. Biol. Chem. 2015, 290, 23916–23934. [Google Scholar] [CrossRef] [PubMed]
  17. Sanglard, D. Emerging threats in antifungal-resistant fungal pathogens. Front. Med. 2016, 3, 11. [Google Scholar] [CrossRef] [PubMed]
  18. Ghannoum, M.A.; Rice, L.B. Antifungal agents: Mode of action, mechanisms of resistance, and correlation of these mechanisms with bacterial resistance. Clin. Microbiol. Rev. 1999, 12, 501–517. [Google Scholar] [PubMed]
  19. Shaikh, N.; Hussain, K.A.; Petraitiene, R.; Schuetz, A.N.; Walsh, T.J. Entomophthoramycosis: A neglected tropical mycosis. Clin. Microbiol. Infect. 2016, 22, 688–694. [Google Scholar] [CrossRef] [PubMed]
  20. Page, R.M.; Humber, R.A. Phototropism in Conidiobolus coronatus. Mycologia 1973, 65, 335–354. [Google Scholar] [CrossRef] [PubMed]
  21. Smith, M.F.; Callaghan, A.A. Quantitative survey of Conidiobolus and Basidiobolus in soils and litter. Trans. Br. Mycol. Soc. 1987, 89, 179–185. [Google Scholar] [CrossRef]
  22. Subramanian, C.; Sobel, J.D. A case of Conidiobolus coronatus in the vagina. Med. Mycol. 2011, 49, 427–429. [Google Scholar] [CrossRef] [PubMed]
  23. Chowdhary, A.; Randhawa, H.S.; Khan, Z.U.; Ahmad, S.; Khanna, G.; Gupta, R.; Chakravarti, A.; Roy, P. Rhinoentomophthoromycosis due to Conidiobolus coronatus. A case report and an overview of the disease in India. Med. Mycol. 2010, 48, 870–879. [Google Scholar] [CrossRef] [PubMed]
  24. Dutta, S.; Sarkar, S.; Linka, U.; Dora, S. Conidiobolomycosis: A case report of rare fungal infection from the eastern India. Indian Dermatol. Online J. 2015, 6, 393–395. [Google Scholar] [CrossRef] [PubMed]
  25. Hoogendijk, C.F.; Pretorius, E.; Marx, J.; van Heerden, W.E.P.; Imhof, A.; Schneemann, M. Detection of villous conidia of Conidiobolus coronatus in a blood sample by scanning electron microscopy investigation. Ultrastruct. Pathol. 2006, 30, 53–58. [Google Scholar] [CrossRef] [PubMed]
  26. Gugnani, H.C. A review of zygomycosis due to Basidiobolus ranarum. Eur. J. Epidemiol. 1999, 15, 923–929. [Google Scholar] [CrossRef] [PubMed]
  27. Rippon, J.W. Medical Mycology: The Pathogenic Fungi and the Pathogenic Actinomycetes, 3rd ed.; W.B. Saunders Co.: Philadelphia, PA, USA, 1988; pp. 506–531. ISBN 0-7216-2444-8. [Google Scholar]
  28. Kwon-Chung, K.J.; Bennett, J.E. Medical Mycology; Lea and Febiger: Philadelphia, PA, USA, 1992; p. 866. ISBN 0-8121-1463-9. [Google Scholar]
  29. Ribes, J.A.; Vanover-Sams, C.L.; Baker, D.J. Zygomycetes in human disease. Clin. Microbiol. Rev. 2000, 13, 236–301. [Google Scholar] [CrossRef] [PubMed]
  30. Ellis, D.H. Systemic zygomycosis. In Topley and Wilson’ Microbiology and Microbial Infections. Medical Mycology, 10th ed.; Merz, W.G., Hay, R.J., Eds.; Edward Arnold: London, UK, 2005; pp. 659–686. [Google Scholar]
  31. Paszkiewicz, M.; Tyma, M.; Ligęza-Żuber, M.; Włóka, E.; Boguś, M.I.; Stepnowski, P. Mycotoxin production by entomopathogenic fungus Conidiobolus coronatus. Int. J. Environ. Agric. Res. 2017, 3, 33–40. [Google Scholar]
  32. Paszkiewicz, M.; Tyma, M.; Ligęza-Żuber, M.; Włóka, E.; Boguś, M.I.; Stepnowski, P. Trichothecenes production by entomopathogenic fungus Conidiobolus Coronatus. Adv. Toxicol. Toxic Eff. 2016, 1, 7–14. [Google Scholar] [CrossRef]
  33. Bogus, M.I.; Maria, S. Histopathology of Conidiobolus coronatus (entomophthorales) infection in Galleria mellonella (lepidoptera) larvae. Acta Parasitol. 2000, 45, 48–54. [Google Scholar]
  34. Malinowski, H. Entomopathogenic fungi as insecticides in protecting the forest. Prog. Plant Prot. 2009, 49, 865–873. [Google Scholar]
  35. Chang, Y.; Wang, S.; Sekimoto, S.; Aerts, A.L.; Choi, C.; Clum, A.; LaButti, K.M.; Lindquist, E.A.; Yee Ngan, C.; Ohm, R.A.; et al. Phylogenomic analyses indicate that early fungi evolved digesting cell walls of algal ancestors of land plants. Genome Biol. Evol. 2015, 7, 1590–1601. [Google Scholar] [CrossRef] [PubMed]
  36. Valle, A.C.; Wanke, B.; Lazéra, M.S.; Monteiro, P.C.; Viegas, M.L. Entomophthoramycosis by Conidiobolus coronatus. Report of a case successfully treated with the combination of itraconazole and fluconazole. Rev. Inst. Med. Trop. Sao Paulo 2001, 43, 233–236. [Google Scholar] [CrossRef] [PubMed]
  37. Choon, S.E.; Kang, J.; Neafie, R.C.; Ragsdale, B.; Klassen-Fischer, M.; Carlson, J.A. Conidiobolomycosis in a young Malaysian woman showing chronic localized fibrosing leukocytoclastic vasculitis: A case report and meta-analysis focusing on clinicopathologic and therapeutic correlations with outcome. Am. J. Dermatopathol. 2012, 34, 511–522. [Google Scholar] [CrossRef] [PubMed]
  38. Guarro, J.; Aguilar, C.; Pujol, I. In-vitro antifungal susceptibilities of Basidiobolus and Conidiobolus spp. strains. J. Antimicrob. Chemother. 1999, 44, 557–560. [Google Scholar] [CrossRef] [PubMed]
  39. Xiao, G.; Ying, S.H.; Zheng, P.; Wang, Z.L.; Zhang, S.; Xie, X.Q.; Shang, Y.; Leger, R.J.S.; Zhao, G.P.; Wang, C.; et al. Genomic perspectives on the evolution of fungal entomopathogenicity in Beauveria bassiana. Sci. Rep. 2012, 2, 483. [Google Scholar] [CrossRef] [PubMed]
  40. Zheng, P.; Xia, Y.; Xiao, G.; Xiong, C.; Hu, X.; Zhang, S.; Zheng, H.; Huang, Y.; Zhou, Y.; Wang, S.; et al. Genome sequence of the insect pathogenic fungus Cordyceps militaris, a valued traditional Chinese medicine. Genome Biol. 2011, 12, R116. [Google Scholar] [CrossRef] [PubMed]
  41. Gao, Q.; Jin, K.; Ying, S.H.; Zhang, Y.; Xiao, G.; Shang, Y.; Duan, Z.; Hu, X.; Xie, X.Q.; Zhou, G.; et al. Genome sequencing and comparative transcriptomics of the model entomopathogenic fungi Metarhizium anisopliae and M. acridum. PLoS Genet. 2011, 7, e1001264. [Google Scholar] [CrossRef] [PubMed]
  42. Cuomo, C.A.; Rodriguez-Del Valle, N.; Perez-Sanchez, L.; Abouelleil, A.; Goldberg, J.; Young, S.; Zeng, Q.; Birren, B.W. Genome sequence of the pathogenic fungus Sporothrix schenckii (ATCC 58251). Genome Announc. 2014, 2. [Google Scholar] [CrossRef] [PubMed]
  43. Lepesheva, G.I.; Waterman, M.R. Sterol 14α-demethylase cytochrome P450 (CYP51), a P450 in all biological kingdoms. Biochim. Biophys. Acta 2007, 1770, 467–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kelly, S.L.; Lamb, D.C.; Baldwin, B.C.; Corran, A.J.; Kelly, D.E. Characterization of Saccharomyces cerevisiae CYP61, Sterol Δ22-desaturase, and inhibition by azole antifungal agents. J. Biol. Chem. 1997, 272, 9986–9988. [Google Scholar] [CrossRef] [PubMed]
  45. Gotoh, O. Substrate recognition sites in cytochrome P450 family 2 (CYP2) proteins inferred from comparative analyses of amino acid and coding nucleotide sequences. J. Biol. Chem. 1992, 267, 83–90. [Google Scholar] [PubMed]
  46. Sirim, D.; Widmann, M.; Wagner, F.; Pleiss, J. Prediction and analysis of the modular structure of cytochrome P450 monooxygenases. BMC Struct. Biol. 2010, 10, 34. [Google Scholar] [CrossRef] [PubMed]
  47. Gricman, Ł.; Vogel, C.; Pleiss, J. Conservation analysis of class-specific positions in cytochrome P450 monooxygenases: Functional and structural relevance. Proteins 2014, 82, 491–504. [Google Scholar] [CrossRef] [PubMed]
  48. Syed, K.; Mashele, S.S. Comparative analysis of P450 signature motifs EXXR and CXG in the large and diverse kingdom of fungi: Identification of evolutionarily conserved amino acid patterns characteristic of P450 family. PLoS ONE 2014, 9, e95616. [Google Scholar] [CrossRef] [PubMed]
  49. Sello, M.M.; Jafta, N.; Nelson, D.R.; Chen, W.; Yu, J.H.; Parvez, M.; Kgosiemang, I.K.R.; Monyaki, R.; Raselemane, S.C.; Qhanya, L.B.; et al. Diversity and evolution of cytochrome P450 monooxygenases in Oomycetes. Sci. Rep. 2015, 5, 11572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Feyereisen, R. Arthropod CYPomes illustrate the tempo and mode in P450 evolution. Biochim. Biophys. Acta 2011, 1814, 19–28. [Google Scholar] [CrossRef] [PubMed]
  51. Parvez, M.; Qhanya, L.B.; Mthakathi, N.T.; Kgosiemang, I.K.R.; Bamal, H.D.; Pagadala, N.S.; Xie, T.; Yang, H.; Chen, H.; Theron, C.W.; et al. Molecular evolutionary dynamics of cytochrome P450 monooxygenases across kingdoms: Special focus on mycobacterial P450s. Sci. Rep. 2016, 6, 33099. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Syed, K.; Shale, K.; Nazir, K.N.H.; Krasevec, N.; Mashele, S.S.; Pagadala, N.S. Genome-wide identification, annotation and characterization of novel thermostable cytochrome P450 monooxygenases from the thermophilic biomass-degrading fungi Thielavia terrestris and Myceliophthora thermophila. Genes Genom. 2014, 36, 321–333. [Google Scholar] [CrossRef]
  53. Marchler-Bauer, A.; Bo, Y.; Han, L.; He, J.; Lanczycki, C.J.; Lu, S.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; et al. CDD/SPARCLE: Functional classification of proteins via subfamily domain architectures. Nucleic Acids Res. 2016, 45, D200–D203. [Google Scholar] [CrossRef] [PubMed]
  54. Nelson, D.R. The cytochrome p450 homepage. Hum. Genom. 2009, 4, 59–65. [Google Scholar]
  55. Nelson, D.R.; Kamataki, T.; Waxman, D.J.; Guengerich, F.P.; Estabrook, R.W.; Feyereisen, R.; Gonzalez, F.J.; Coon, M.J.; Gunsalus, I.C.; Gotoh, O.; et al. The P450 superfamily: Update on new sequences, gene mapping, accession numbers, early trivial names of enzymes, and nomenclature. DNA Cell Biol. 1993, 12, 1–51. [Google Scholar] [CrossRef] [PubMed]
  56. Nelson, D.R. Cytochrome P450 nomenclature. Methods Mol. Biol. 1998, 107, 15–24. [Google Scholar] [PubMed]
  57. Nelson, D.R. Cytochrome P450 Nomenclature, 2004. Methods Mol. Biol. 2006, 320, 1–10. [Google Scholar] [CrossRef] [PubMed]
  58. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed]
  59. Trifinopoulos, J.; Nguyen, L.T.; von Haeseler, A.; Minh, B.Q. W-IQ-TREE: A fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 2016, 44, W232–W235. [Google Scholar] [CrossRef] [PubMed]
  60. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed]
  61. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef] [PubMed]
  62. Sievers, F.; Higgins, D.G. Clustal Omega, accurate alignment of very large numbers of sequences. Methods Mol. Biol. 2014, 1079, 105–116. [Google Scholar] [CrossRef] [PubMed]
  63. Crooks, G.E.; Hon, G.; Chandonia, J.M.; Brenner, S.E. WebLogo: A sequence logo generator. Genome Res. 2004, 14, 1188–1190. [Google Scholar] [CrossRef] [PubMed]
Ijms 19 01711 g001 550
Figure 1. Phylogenetic analysis of C. coronatus P450s. P450 families bloomed in C. coronatus are highlighted in different colors. Bootstrap replications used to evaluate each node of the phylogenetic tree are presented in Figure S1.
Figure 1. Phylogenetic analysis of C. coronatus P450s. P450 families bloomed in C. coronatus are highlighted in different colors. Bootstrap replications used to evaluate each node of the phylogenetic tree are presented in Figure S1.
Ijms 19 01711 g001
Ijms 19 01711 g002 550
Figure 2. Analysis of amino acid patterns at the EXXR and CXG motifs in P450 families CYP5854-CYP5856. P450 protein sequences used to deduce EXXR and CXG signature sequences were presented in Table S4.
Figure 2. Analysis of amino acid patterns at the EXXR and CXG motifs in P450 families CYP5854-CYP5856. P450 protein sequences used to deduce EXXR and CXG signature sequences were presented in Table S4.
Ijms 19 01711 g002
Ijms 19 01711 g003 550
Figure 3. Comparative analysis of P450 family members in C. coronatus. The P450 family name and percentage in the total number of P450s are shown in the figure. The total number of P450s in a P450 family is presented in Table 1.
Figure 3. Comparative analysis of P450 family members in C. coronatus. The P450 family name and percentage in the total number of P450s are shown in the figure. The total number of P450s in a P450 family is presented in Table 1.
Ijms 19 01711 g003
Ijms 19 01711 g004 550
Figure 4. Analysis of tandem duplications of P450s in C. coronatus. The numbers next to bars indicate the number of P450s located on the scaffold.
Figure 4. Analysis of tandem duplications of P450s in C. coronatus. The numbers next to bars indicate the number of P450s located on the scaffold.
Ijms 19 01711 g004
Ijms 19 01711 g005 550
Figure 5. Gene-structure analysis of P450s in the CYP5854 family. Each P450 gene structure is presented with the P450 name, exons (red color bars), and introns (gap between bars). The size of the exons (number at the top of the red bars) and introns (number in the gaps) is also shown in the figure. Some P450s’ evolution was deduced from other P450s and their origin is indicated with dotted lines. Blue regions indicate an untranslated DNA region. The location of the respective P450s is indicated with their scaffold number.
Figure 5. Gene-structure analysis of P450s in the CYP5854 family. Each P450 gene structure is presented with the P450 name, exons (red color bars), and introns (gap between bars). The size of the exons (number at the top of the red bars) and introns (number in the gaps) is also shown in the figure. Some P450s’ evolution was deduced from other P450s and their origin is indicated with dotted lines. Blue regions indicate an untranslated DNA region. The location of the respective P450s is indicated with their scaffold number.
Ijms 19 01711 g005
Ijms 19 01711 g006 550
Figure 6. Gene structure analysis of P450s in CYP5857 and CYP5860 family. Each P450 gene structure is presented with the P450 name, exons (red color bars) and introns (gap between bars). The size of exons (number at the top of the red bars) and introns (number in the gaps) is shown in the figure. Some P450s’ evolution was deduced from other P450s and their origin is indicated with dotted lines. The location of the respective P450s is indicated with their scaffold number.
Figure 6. Gene structure analysis of P450s in CYP5857 and CYP5860 family. Each P450 gene structure is presented with the P450 name, exons (red color bars) and introns (gap between bars). The size of exons (number at the top of the red bars) and introns (number in the gaps) is shown in the figure. Some P450s’ evolution was deduced from other P450s and their origin is indicated with dotted lines. The location of the respective P450s is indicated with their scaffold number.
Ijms 19 01711 g006
Ijms 19 01711 g007 550
Figure 7. Comparative analysis of P450 diversity percentage among entomopathogenic fungi. The corresponding values are presented in Table S6.
Figure 7. Comparative analysis of P450 diversity percentage among entomopathogenic fungi. The corresponding values are presented in Table S6.
Ijms 19 01711 g007
Table
Table 1. Annotation (family and subfamily assignment) of C. coronatus P450s.
Table 1. Annotation (family and subfamily assignment) of C. coronatus P450s.
P450 FamilyP450 SubfamilyNo. of P450s in a Family
Subfamily NameNo. of P450s
CYP51F11
CYP5854A848
B4
C1
D1
E34
CYP5855 *A536
B2
C12
D15
E1
1
CYP5856A520
B15
CYP5857A1010
CYP5858A1111
CYP5859A55
CYP5860A33
CYP5861A33
CYP5862A22
CYP5863A11
CYP5864A11
CYP-fragment #-11
12 Families21 subfamilies142142
* Belonging to the same family but not assigned to a subfamily because of short amino acid sequence; - not applicable; # not assigned to family because of short amino acid sequences (284 amino acids) and only showed 24% identity to named P450s.
Table
Table 2. Comparative P450 analysis in entomopathogenic and animal (including human) pathogenic fungi.
Table 2. Comparative P450 analysis in entomopathogenic and animal (including human) pathogenic fungi.
FungusHostP450 CountNo. of P450 FamiliesReference
Entomopathogenic fungus
Beauveria bassianaArthropods (termites, thrips, whiteflies, aphids, and beetles)8349[39]
Cordyceps militarisButterflies and caterpillars5737[40]
Metarhizium acridum≥200 insects10067[41]
Metarhizium anisopliae (formerly Metarhizium robertsii)Locusts12387[41]
Animal including human pathogen
Sporothrix schenckiiHumans and other animals (cats, dogs, rodents, squirrels, horses, and birds)4032[7,42]
Animals (including humans) and entomopathogenic fungus
Conidiobolus coronatusInsects, humans, and other animals (horses, sheep and dogs)14212This work

Share and Cite

MDPI and ACS Style

Ngwenya, M.L.; Chen, W.; Basson, A.K.; Shandu, J.S.; Yu, J.-H.; Nelson, D.R.; Syed, K. Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus. Int. J. Mol. Sci. 2018, 19, 1711. https://doi.org/10.3390/ijms19061711

AMA Style

Ngwenya ML, Chen W, Basson AK, Shandu JS, Yu J-H, Nelson DR, Syed K. Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus. International Journal of Molecular Sciences. 2018; 19(6):1711. https://doi.org/10.3390/ijms19061711

Chicago/Turabian Style

Ngwenya, Mathula Lancelot, Wanping Chen, Albert Kotze Basson, Jabulani Siyabonga Shandu, Jae-Hyuk Yu, David R. Nelson, and Khajamohiddin Syed. 2018. "Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus" International Journal of Molecular Sciences 19, no. 6: 1711. https://doi.org/10.3390/ijms19061711

APA Style

Ngwenya, M. L., Chen, W., Basson, A. K., Shandu, J. S., Yu, J. -H., Nelson, D. R., & Syed, K. (2018). Blooming of Unusual Cytochrome P450s by Tandem Duplication in the Pathogenic Fungus Conidiobolus coronatus. International Journal of Molecular Sciences, 19(6), 1711. https://doi.org/10.3390/ijms19061711

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop