See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/268807190 Comparative Genomics and Evolutionary Analysis of Cytochrome P450 Monooxygenases in Fungal Subphylum Saccharomycotina Article  in  Journal of Pure and Applied Microbiology · November 2014 CITATIONS 6 READS 298 3 authors: Some of the authors of this publication are also working on these related projects: Molecular modelling View project Ipeleng Kgosiemang Central University of Technology 13 PUBLICATIONS   28 CITATIONS    SEE PROFILE Sam Mashele Central University of Technology 26 PUBLICATIONS   97 CITATIONS    SEE PROFILE Khajamohiddin Syed University of Zululand 51 PUBLICATIONS   1,403 CITATIONS    SEE PROFILE All content following this page was uploaded by Khajamohiddin Syed on 27 November 2014. The user has requested enhancement of the downloaded file. JOURNAL OF PURE & APPLIED MICROBIOLOGY, Nov. 2014. Vol. 8(Spl. Edn. 2), p. 291-302 * To whom all correspondence should be addressed. Tel.: +27 51 507 3425;Fax: +27 51 507 3134; Email: khajamohiddinsyed@gmail.com Comparative Genomics and Evolutionary Analysis of Cytochrome P450 Monooxygenases in Fungal Subphylum Saccharomycotina Ipeleng Kopano Rosinah Kgosiemang, Samson Sitheni Mashele and Khajamohiddin Syed* Department of Health Sciences, Faculty of Health and Environmental Sciences, Central University of Technology, Bloemfontein 9300, Free State, South Africa. (Received: 09 September 2014; accepted: 29 October 2014) Cytochrome P450 monooxygenases (P450s) are heme-thiolate enzymes and play an important role in the primary and secondary metabolism of living organisms. Genome sequencing analysis of fungal organisms revealed the presence of numerous P450s in their genomes, with few exceptions. P450s in the fungal subphylum Saccharomycotina, which contains biotechnologically important and opportunistic human pathogen yeasts, have been underexplored because there are few P450s in their genomes. In the present study we performed comparative analysis of P450s in 25 yeast species. A hundred and seventy-two P450s were found in 25 yeast species and these are grouped into 13 P450 families and 27 subfamilies. P450s ranged from a minimum of three (Saccharomyces species) to a maximum of 21 (Candida species) in the yeast genomes. Among the P450 families, the CYP52 family showed the highest number of member P450s (71) followed by CYP51 (27), CYP61 (25), CYP56 (20) and CYP501 (11). Pichia pastoris and Dekkera bruxellensis showed a novel P450 family, CYP5489, in their genome. Based on the functional properties of characterized P450s, we conclude that P450s in Saccharomycotina species possibly play a role in organisms’ physiology either in the synthesis of cellular components or in the utilization of simpler organic molecules. The ecological niches of yeast species are highly enriched with simpler organic nutrients and it is well known that yeast species utilize simpler organic nutrients as carbon source efficiently. This might have played a role in compacting yeast genomes and possibly losing a considerable number of P450s during evolution. Key words: Adaptation, Comparative genomics, Cytochrome P450 monooxygenases, Ecological niches, Evolution, Fungi, Saccharomycotina, Yeast. The fungal kingdom is the largest biological kingdom and consists of diverse organisms that are adapted to diverse ecological niches. This kingdom is classified into four phyla, namely Ascomycota, Basidiomycota, Zygomycota and Chytridiomycota1. The phylum Ascomycota is further classified into three subphyla, namely Pezizomycotina, Saccharomycotina and Taphrinomycotina2. Our interest lies in the Saccharomycotina subphylum, which comprises of yeast organisms that are well known for their potential biotechnological value, as well as opportunistic human pathogens. Table 1 shows a list of Saccharomycotina species, provides general information and indicates the importance of these species. Genome sequencing analysis of fungal species revealed the presence of numerous cytochrome P450 monooxygenases (P450s) in their genomes, with few exceptions3. P450s are heme- thiolate enzymes and their role in organisms’ primary and secondary metabolism and their potential use in biotechnology, bioremediation, pharmacology and biofuel generation has been J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 292 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA documented4. Fungal P450s occupy a special place because of their diverse catalytic activities compared to the P450s from other biological kingdoms5. Furthermore, the presence of large numbers of P450s and diverse families in fungi makes fungal P450s ideal for the study of P450s’ evolution. Since the genome sequence of fungi became available6, study on fungal P450s, particularly on functional and evolutionary aspects, has gained momentum. Study on fungal P450s revealed that basidiomycetes enriched P450s in their genome owing to duplication to help the organism to adapt to diverse ecological niches, such as wood degradation3,7. Study also revealed that ascomycetes, particularly species belonging to the subphylum Pezizomycotina, contain diverse P450 families in their genome8. Furthermore, fungal P450s were used as model P450s to identify P450 family-specific signature sequences9. Authors have suggested that these signature sequences evolved during the evolution of P450 families from a common P450 ancestor. Recent studies also revealed the presence of a large number of thermostable P450s in fungi10. P450 evolutionary studies or genome-wide comparative P450 studies on fungal P450s were based on P450s from species belonging to the phyla Basidiomycotoa and Ascomycota. In the phylum Ascomycota only P450s from the subphylum Pezizomycotina have been explored for both evolutionary and functional characterization3,8. P450s from Saccharomycotina species are under explored and to date no study on comparative P450s in Saccharomycotina species has been carried out. Studies have been limited to describing the number of P450s in Saccharomycotina species and detailing one or two P450 families. Furthermore, genome sequencing of new Saccharomycotina species necessitated the performance of comparative P450 profiling in their genomes. Considering the importance of P450s in general and the use of Saccharomycotina species in biotechnology and their pathogenicity towards animals, particularly humans, it is important to understand the role of P450s in yeast species. In the present study we performed genome-wide comparative P450s analysis in 25 yeast species and P450s in three yeast species, Pichia pastoris, Dekkera bruxellensis and Pichia anomala, were identified and annotated as per International Cytochrome P450 Nomenclature Committee11. MATERIALS AND METHODS Saccharomycotina species A total of 25 Saccharomycotina species whose genomes have been sequenced and are publicly available were used in the study. A list of the yeast species used in this study is shown in Table 1. Genome-data mining of Saccharomycotina species for P450s The 22 Saccharomycotina species P450s were downloaded from the publicly available Cytochrome P450 Homepage12. The downloaded P450s were cross checked at the web pages: http:/ /www.yeastgenome.org/13 and http:// www.candidagenome.org/14. P450s in the remaining three yeast species,i.e. P. pastoris, D. bruxellens and P. anomala, were annotated following the standard procedure described in our recent studies7,10. Briefly, the proteome of yeast species was downloaded from their genome data base (Dekkera bruxellensis CBS 2499 v2.0: http:// genome.jgi.doe.gov/ Dekbr2/Dekbr2.home.html; Pichia Pastoris: http://bioinformatics. psb.ugent.be/orcae/;Pichia anomala[currently named Wickerhamomyces anomalus]NRRL Y-366- 8 v1.0: http://genome.jgi.doe. gov/Wican1/ Wican1.home.html) and the whole proteome was subjected to functional annotation using NCBI Batch Web CD-search tool15. The proteins grouped under the P450 superfamily were selected and analyzed for the presence of the P450 family signature motifs, namely EXXR and CXG. The proteins that showed both motifs were considered authentic P450s and used in this study. Annotation and classification of P450s The proteins identified as P450 were subjected to BLAST analysis against all named fungal species at the Cytochrome P450 Homepage12. For each P450, the closest homolog was identified and based on the homology percentage, family and subfamily names were assigned. For assigning the family and subfamily names, the standard rule set by the International P450 Nomenclature Committee was followed, i.e. P450s within a family share more than 40% amino acid homology and members of subfamilies share more than 55% amino acid homology11. Furthermore, P450s that showed less than 40% homology with known P450s were assigned to a J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 293KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA new family with the help of the P450 nomenclature of Dr David R Nelson, University of Tennessee Health Science Center, Memphis, Tennessee, USA. Phylogenetic analysis of P450s Phylogenetic analysis of P450s was carried out in the same way as described in our recent publications7,10. Briefly, evolutionary analysis was carried out using the minimum evolution method16. The phylogenetic analysis was carried out using Molecular Evolutionary Genetics Analysis (MEGA 5.05) software17. Functional analysis of P450s Considering the large number of P450s used in this study and the availability of functional data, we performed a literature survey on the functional analysis of Saccharomycotina species P450s and their role in organisms’ physiology. RESULTS AND DISCUSSION P450ome of P. pastoris, D. bruxellens and P. anomala The yeast species P. pastoris and D. bruxellens showed four P450s and P. anomala showed six P450s in their genome (Table 2). Among the P450 families found in these species CYP51, CYP61 and CYP501 are common in three species. Compared to two other yeast species, P450s belonging to CYP56, CYP5205 and CYP5217 families are only present in P. anomala. Surprisingly, P. pastoris and D. bruxellens showed a novel P450 family, i.e. CYP5489, in their genome. However, a difference was found in the CYP5489 subfamily type in yeast species. P. pastoris contained P450 belonging to subfamily A and D. bruxellens contained P450 belonging to subfamily B. It is noteworthy that this P450 family is not found in other yeast species (Table 2). Comparative analysis of P450ome in yeast species A hundred and seventy-two P450s were found in 25 yeast species (Table 2). The P450s ranged from a minimum of three to a maximum of 21 in the yeast genomes; 172 P450s were grouped under 13 P450 families and 27 subfamilies (Figure 1 and Table 2). Among the P450 families, CYP52 family showed the highest number of member P450s (71) followed by CYP51 (27), CYP61 (25), CYP56 (20) and CYP501 (11) (Figure 2). The remaining P450 families showed less than 6% of member P450s (Figure 2). Fig. 1. Phylogenetic analysis of the P450s in subphyla Saccharomycotina. In total, 172 P450 sequences were included in the tree. A minimal evolution tree was constructed using the close-neighbor-interchange algorithm in MEGA (version 5.05). For ease of visual identity, P450s belonging to different families were presented in unique colors. J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 294 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA Table 1.List of yeast species selected for the study and general information and importance of the species. Species name General information and importance Reference Saccharomyces The most useful yeast, which has been instrumental in winemaking, 24,25 cerevisiae baking and brewing since ancient times. It is well known as baker’s yeast and it was the first eukaryote ever to have its genome completely sequenced. It is the microorganism behind the most common type of fermentation. It is one of the most intensively studied eukaryotic model organisms in molecular and cell biology, much like Escherichia coli as the model bacterium. It is the organism employed most in industry in terms of production volumes. Saccharomyces Non-domesticated yeast living on the bark of deciduous trees. 26 paradoxus Closest relative of S. cerevisiae hence it is an attractive model organism for the study of the population genetics and genomics of wild yeast. Saccharomyces Belongs to Saccharomyces sensustricto complex. Shows high 26 mikatae conservation of synteny with S. cerevisiae. Saccharomyces Model organism for telomere research since many of the features of 27 castellii the telomere biology of this yeast are closer to those of humans than other yeast. Saccharomyces Used in winemaking, cider fermentation and making distilled 26, 27 bayanus beverages. It is also used extensively for comparative genomic studies, including expression patterns and nucleosome profile. Saccharomyces Belongs to Saccharomyces sensustricto complex. Used in making 27, 28 kudriavzevii beer and wine. It can grow at low temperatures. Saccharomyces More efficient user of glucose than model yeast S. cerevisiae. It is 27 kluyveri used in industrial applications, such as the production of proteins, as its biomass yield is greater than that of S. cerevisiae. Lives in diverse environments, for example this yeast is a plant pathogen and has been isolated from drosophila species and tissues of an HIV- infected patient. Candida It is one of the major opportunistic pathogenic fungi causing 29, 30, 31 albicans systemic infection (candidiasis) in immuno-compromised and immuno-competent hosts. It is a diploid fungus that grows both as yeast and filamentous cells and a casual agent of opportunistic oral and genital infections (superficial infections) in humans. Candida The most prevalent opportunistic pathogenic yeast species of the 19, 31, 32 tropicalis Candida-non-albicans group. Human pathogen that is well known as medical yeast pathogen. It causes Candida-non-albican candidiasis. It is a glucose and maltose fermenting yeast. The industrial strain of C. tropicalis is engineered to produce long-chain- dicarboxylic acids. Candida The model organism of a group of so-called “flavinogenic yeasts” 32, 33 guilliermondii capable of riboflavin oversynthesis. It often causes (anamorph of onychomycosis and is rarely involved in invasive fungal Pichia guilliermondii) infections. Candida It was thought to be a primarily non-pathogenic organism of the 34 glabrata normal flora of healthy individuals. However, with the increasing population of immuno-compromised individuals, trends have shown C. glabrata to be a highly opportunistic pathogen of the urogenital tract and of the bloodstream. Candida A fungal species of the yeast family that has become a major cause 31 parapsilosis of sepsis and of wound and tissue infections in immuno- compromised patients. Among Candida species this yeast is a particular problem in neonates, transplant recipients and patients J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 295KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA receiving parenteral nutrition. Candida It is the most closely related species to C. albicans and causes 35 dubliniensis superficial infections. Compared to C. albicans, this yeast has low virulence and a longer survival rate in hosts. Candida Responsible for less than 5% of invasive infections among Candida 31, 36 lusitaniae species. Patients who undergo bone marrow transplantation and high-dose cytoreductive chemotherapy have been identified as being at risk of infections caused by this organism. Kluyveromyces Best alternative yeast for genetics and physiology. Its name comes 37 lactis from the ability to assimilate lactose and convert it to lactic acid. K. lactis is one of the main organisms grown in industry in fermenters to produce chymosin on a large scale. Chymosin is used for cheese production. This yeast has GRAS (generally regarded as safe) status. Kluyveromyces Whole genome duplicated yeast. Distant lineage to S. cerevisiae and 38 polysporus considered as best yeast to study comparative whole genome duplication event. Kluyveromyces This yeast is a protoploid or pre whole genome duplication budding 39 waltii (recently yeast descended from a lineage and did not undergo the whole named genome duplication event. Genome sequencing of this yeast as revealed how the genome of S. cerevisiae evolved via the whole Lachancea waltii) genome duplication event. Pichia stipites Capable of both aerobic and oxygen-limited fermentation and has the 40 highest known natural ability of any microbe to ferment xylose directly, converting it to ethanol. Lodderomyces Causes blood stream infection in humans. It is the second or third 31 elongisporus most common yeast species isolated from patients with bloodstream infections in Europe, Canada, Asia and Latin America. Ashbya gossypii Contains the smallest genome of a free-living eukaryote and because 41 of this is considered as model organism to study filamentous growth. It infects plants such as cotton and citrus. It is a natural overproducer of riboflavin (vitamin B2), which protects its spores against ultraviolet light. This makes it an interesting organism for industries, where genetically modified strains are still used to produce this vitamin. Debaryomyces Metabolically versatile, osmotolerant and oleaginous. It is used for 31 hansenii surface ripening of cheese and meat products. It one of the most halotolerant species capable of growing in media containing NaCl as high as 4 M. Yarrowia Most extensively studied “non-conventional” yeast (strictly aerobic) 4, 43 lipolytica and currently used as model organism for study of protein secretion, peroxisome biogenesis, dimorphism, degradation of hydrophobic substrates. It produces important metabolites and has developed intense secretory activity. Capable of growing on organic compounds and producing bio surfactants. This yeast has been used in industry (as a biocatalyst) and in academic studies. This yeast has GRAS status. Dekkera Although this yeast is distantly related to S. cerevisiae, is commonly 44 bruxellensis found in the same habitat. It competes with S. cerevisiae in fermentation, can produce high amounts of ethanol and grow without oxygen. It plays a key role in flavor development in lambic beer. However, in the wine industry, this yeast is considered spoilage yeast owing to wine spoilage by the generation of a high amount of phenolic compounds (4-ethylguaiacol and 4-ethylphenol). Table 1 Continue J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 296 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA Pichia This methylotrophic yeast is the most commonly used yeast species 45 pastoris in the production of recombinant proteins and is widely used to produce proteins for basic research and medical applications. It is a model organism for the study of peroxisomal proliferation and methanol assimilation. Wickerhamomyces It exhibits a multitude of biotechnologically important characteristics 46 anomalus in flavor enhancement, food processing, biopreservation, (formerly dairy fermentation and waste water treatment. Pichia anomala and Hansenula anomala) Table 1 Continue Among the 13 P450 families, CYP51 and CYP61 were conserved across the 25 yeast species. A single copy of CYP61 members were found in the yeast species, whereas CYP51 members were found in more than a single copy in two yeast species , Candida guillermondii and Kluyveromyces polysporus (Table 2). The presence of more than one copy of CYP51 members is not rare and species belonging to other biological kingdoms, such as plants, also showed more than once copy of CYP51 members in their genomes18. CYP52 family members are enriched in species belonging to the genus Candida and in Yarrowia lipolytica. Among candida species C. tropicalis showed the maximum number of CYP52 members (17 members) in its genome and Y. lipolytica showed 12 members in its genome (Table 2). Phylogenetic analysis of CYP52 members (Figure 1) resulted in the grouping of Y. lipolytica CYP52 members together, suggesting possible duplication of member P450s in its genome after speciation. Among yeast species only P. stipites, Lodderomyces elongisporus and Debaryomyces hansenii showed CYP52 members in their genome. It is interesting to note that P. anomala did not contain CYP52 members compared to P. stipites. The presence of CYP52 members in selected yeast species suggested that these P450s play a role in Fig. 2. Comparative analysis of P450s in subphyla Saccharomycotina. A total of 172 P450s representing all P450 families found in Saccharomycotina were used for analysis. The P450 family, number of member P450s in a family and percentage compared to overall P450 count (172 P450s) are also shown in the pie chart. J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 297KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA T ab le 2 . G en om e- w id e co m pa ra ti ve f am il y an d su b- fa m il y le ve l a na ly si s of c yt oc hr om e P 45 0 m on oo xy ge na se s (P 45 0s ) in s ub ph yl um S ac ch ar om yc ot in a. C Y P f am il y SF SC SP S M SK SC a SK l SB C A C T C G C G u C P C D C L K L K W K P PS L E A G D H Y L D B P P PA T P C P Y 51 F 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 2 1 1 1 1 1 1 1 1 27 C Y P 52 A 4 14 4 7 4 3 5 5 5 71 B 1 C 1 1 2 1 1 D 1 F 11 S 1 C Y P 56 A 1 1 1 1 1 1 1 1 1 1 1 1 20 D 1 E 1 1 1 1 1 1 1 C Y P 61 A 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 25 C Y P 50 1 A 1 1 1 1 11 B 1 C 1 D 1 1 1 E 1 F 1 C Y P 50 4 A 1 2 G 1 C Y P 54 8 P 1 1 C Y P 52 05 N S 1 1 C Y P 52 17 A 1 1 1 1 1 1 1 1 1 9 C Y P 52 23 A 1 1 C Y P 52 51 A 1 1 C Y P 52 52 A 1 1 C Y P 54 89 A 1 2 B 1 T P 3 3 3 3 3 3 3 10 21 3 10 14 10 8 5 3 4 10 10 3 9 17 4 4 6 17 2 A bb re vi at io ns : S C , Sa cc ha ro m yc es c er ev is ia e; S P, S ac ch ar om yc es p ar ad ox us ; S M , Sa cc ha ro m yc es m ik at ae ; S K , Sa cc ha ro m yc es k ud ri av ze vi i; S C a, S ac ch ar om yc es ca st el li i; S ki , Sa cc ha ro m yc es k lu yv er i; S B , Sa cc ha ro m yc es b ay an us ; C A , C an di da a lb ic an s; C T , C an di da t ro pi ca li s; C G , C an di da g la br at a; C G u, C an di da g ui ll er m on di i; C P, C an di da p ar ps il os is ; C D , C an di da d ub li ne ns is ; C L , C an di da l us it an ia e; K L , K ly ve ro m yc es l ac ti s; K W , K lu yv er om yc es w al ti i; K P, K lu yv er om yc es p ol ys po ru s; P S , P ic hi a st ip it es ; L E , L od de ro m yc es e lo ng is po ru s; P A , P ic hi a an om al a; A G , A sh b ya g o ss yp ii ; D H , D eb a ry o m yc es h a n se n ii ; Y L , Ya rr ow ia l ip ol yt ic a ; D B , D ek ke ra br ux el le ns is ; P P, P ic hi a pa st or is ; T P, t ot al P 45 0 co un t; N S , ne w s ub fa m il y. J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 298 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA Table 3. General function and substrate specificity of fungal P450s. P450s that are functionally characterized and present in Saccharomycotina species are shown in the table. P450 General function Substrate Product(s) Reference CYP51 Cell membrane Lanosterol 4,4-dimethylcholesta-8,14,24- 47 sterol trienol biosynthesis 24-methylene-24,25- 4,4-dimethylfecosterol dihydrolanosterol CYP52 Alkanes and fatty Arachidonic acid 20-Hydroxy-5,8,11,14- 48 acid icosatetraenoic acid hydroxylation Dodecane 1-Dodecanol 2-Dodecanol Hexadecane 1-Hexadecanol 2-Hexadecanol Octadecane 1-Octadecanol 2-Octadecanol Decane 1-Decanol 2-Decanol Tetradecane 1-Tetradecano 2-Tetradecanol Lauric acid 12-Hydroxydodecanoic acid 11- Hydroxydodecanoic acid Myristic acid 14-Hydroxytetradecanoic acid 13-Hydroxytetradecanoic acid Palmitoleic acid 16-Hydroxy-9-hexadecenoic acid 15-Hydroxy-9-hexadecenoic acid Palmitic acid 16-Hydroxyhexadecanoic acid 15-Hydroxyhexadecanoic acid alpha-Linoleic acid 18-hydroxy-alpha-linoleic acid 17-hydroxy-alpha-linoleic acid Linoleic acid 18-hydroxy-linoleic acid 17-hydroxy-linoleic acid Oleic acid 18-Hydroxy-9-octadecenoic acid 17-Hydroxy-9-octadecenoic acid Stearic acid 18-Hydroxyoctadecanoic acid 17-Hydroxyoctadecanoic acid CYP56 Synthesis of component N-formyl tyrosine N,N’-bisformyl dityrosine 49 of outer spore wall layer CYP61 Cell membrane sterol Ergosta 5,7,24(28)-trienol Ergosta 5,7,22,24(28)-tetraenol 50 biosynthesis CYP504 Phenylacetate Phenylacetate 2-Hydroxyphenylacetate 51 degradation 3-Hydroxyphenylacetate Homogentisate 52 3,4-Dihydroxyphenylacetate 2,4,5-Trihydroxyphenylacetate the adaptation of yeast species to particular ecological niches, e.g. utilization of specific nutrients, such as alkanes, as carbon source19,20. The distribution of other P450 families and the number of member P450s across the selected yeast species can be obtained in Table 2. One interesting observation is that the CYP548 P450 family is present only in Y. lipolytica. The number of P450s in the yeast genomes was lower compared to the fungal species from other fungal phyla, Basidiomycota, Zygomycota and Chytridiomycota, and even compared to the subphylum Pezizomycotina12, with few exceptions, suggesting the significant loss of P450s in saccharomycotina species. The phenomenon of the presence of a low number of J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 299KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA P450s in Saccharomycotina species is mentioned in earlier studies3,21,22. It has been shown that a large number of P450s, especially certain P450 families that are highly enriched in Basidiomycota, play a key role in the adaptation of Basidiomycetes to different ecological niches, especially in wood degradation, by performing extraordinary catalytic activities5,7. This suggests that P450 numbers in an organism can be directly linked to the physiology/ecological niches of the organism. Significant loss of P450s and the role of enrichment of certain members of P450 families in few Saccharomycotina species physiology are discussed in the next section. Role of P450s in Saccharomycetes physiology Based on the literature available on the functional characterization of P450s (Table 3) and the characteristic life style of yeast species (Table 1), we conclude that Saccharomycotina species retained P450s that are critical in their physiology. For example, CYP51 and CYP61 P450s play a key role in the synthesis of membrane components and CYP56 is involved in the synthesis of the outer spore wall layer (Table 3). CYP52 family members that are enriched in Candida species and Y. lipolytica are involved in organisms’ primary (oxidation of fatty acids) and secondary (oxidation of alkanes) metabolism. CYP52 P450s of C. tropicalis play a key role in the production of biotechnologically valuable long-chain dicarboxylic acids23. The presence of the highest number of CYP52 members in both yeast species gives the species an advantage to utilize alkanes as a carbon source (adaptation to different ecological niches). A deletion mutant of Y. lipolytica, where 12 CYP52 genes were deleted, was unable to utilize n-alkanes (10-16 carbon length) as a carbon source20. This strongly supports the argument that P450s play a key role in the adaptation of yeast species to diverse ecological niches. CYP504 oxidizes simple aromatic compounds such as phenylacetate and helps organisms to grow on this organic compound. Considering the functional properties of characterized P450s, we conclude that the orphan P450s in Saccharomycotina possibly play a role in organisms’ physiology, either in synthesis of cellular components or utilization of simpler organic molecules. CONCLUSIONS The ecological niches of yeast species are highly enriched with simpler organic nutrients, unlike other fungal species such as Basidiomycetes that degrade wood, a complex polymer, to gain access to nutrients. Because of the complex metabolic process involved in degradation of wood components, Basidiomycetes are enriched with catalytically diverse P450s5,7. This clearly suggests that ecological niches of organisms play a critical role in shaping their genome content. Further evidence of this phenomenon can be obtained from the fact that Cryptococcus neoformans and Tremella mesentrica contain eight and 10 P450s in their genomes compared to other basidiomycetes7. C. neoformans is a well-known human pathogen and T. mesenterica is a parasite fungus. It is shown that these non-wood-degrading Basidiomycetes lost P450s owing to their different ecological niches (life style pattern) compared to the wood- degrading basidiomycetes7. It has been suggested that yeast-forming fungi (Saccharomycotina and Taphynomycotina) have small P450omes, while mycorrhizal relationships and complex nutrient degradation seem to enhance P450 diversification21. Considering the type of ecological niches adapted by yeast species, it is clear that these species do not need a large number of P450s. The presence of abundant and simpler organic nutrients and adaptation of yeast species for efficient utilization of these simpler nutrients further compacted their genomes and they lost a considerable number of P450s in their genome during evolution. ACKNOWLEDGMENTS KS and SSM thank the Central University of Technology (CUT), Bloemfontein, Free State for a grant from the University Research and Innovation Fund. SSM also thanks the National Research Foundation (NRF), South Africa for the research grant. The authors want to thank Dr David R Nelson, University of Tennessee, USA for naming the P450 family and subfamilies. The authors also want to thank Ms Barbara Bradley, Pretoria, South Africa for English language editing. J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 300 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA REFERENCES 1. Alexopoulos, C.J., Mims, C.W., Blackwell M. Introductory mycology, 4thedn. Wiley, New York, USA. 2. Nishida, H., Sugiyama, J. Archiascomycetes: Detection of major new lineage within the Ascomycota. Mycoscience 1994; 35: 361-366. 3. Moktali, V., Park, J., Fedorova-Abrams, N.D., Choi, J., Lee, Y.H., Kang, S. Systemic and searchable classification of cytochrome P450 protein encoded by fungal and oomycetes genome. BMC Genomics 2012; 13: 525. 4. Syed, K., Yadav, J.S. P450 monooxygenases (P450ome) of the model white rot fungus Phanerochaete chrysosporium. Crit. Rev. Microbiol., 2012; 38: 339–363. 5. Syed, K., Porollo, A., Lam, Y.W., Grimmett, P.E., Yadav, J.S. CYP63A2, a catalytically versatile fungal P450 monooxygenase capable of oxidizing higher-molecular-weight polycyclic aromatic hydrocarbons, alkylphenols, and alkanes. Appl. Environ. Microbiol., 2013; 79: 2692–702. 6. Grigoriev, I.V., Cullen, D., Goodwin, S.B., Hibbett, D., Jeffries, T.W., Kubicek, C.P., Kuske, C., Magnuson, J.K., Martin, F., Spatafora, J.W., Tsang, A., Baker, S.E. Fueling the future with fungal genomics. Mycology 2011; 2: 192–209. 7. 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(1), e86683. 8. Deng, J., Carbone, I., Dean, R.A. The evolutionary history of cytochrome P450 genes in four filamentous ascomycetes. BMC Evol. Biol., 2007; 7: 30. 9. 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(4),e95616. 10. Syed, K., Shale, K., Nazir, K.H.M.N.Z., 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. 11. Nelson, D.R. Cytochrome P450 nomenclature, 2004. Methods Mol. Biol., 2006; 320: 1-10. 12. Nelson, D.R. The cytochrome P450 homepage. Hum. Genomics 2009; 4: 59–65. 13. Christie, K.R., Weng, S., Balakrishnan, R., Costanzo, M.C., Dolinski, K., Dwight, S.S., Engel, S.R., Feierbach, B., Fisk, D.G., Hirschman, J.E., Hong, E.L., Issel-Tarver, L., Nash, R., Sethuraman, A., Starr, B., Theesfeld, C.L., Andrada, R., Binkley, G., Dong, Q., Lane, C., Schroeder, M., Botstein, D., Cherry, J.M. Saccharomyces Genome Database (SGD) provides tools to identify and analyze sequences from Saccharomyces cerevisiae and related sequences from other organisms. Nucleic Acids Res., 2004; 1: D311-4. 14. Arnaud, M.B., Costanzo, M.C., Skrzypek, M.S., Binkley, G., Lane, C., Miyasato, S.R., Sherlock, G. The Candida Genome Database (CGD), a community resource for Candida albicans gene and protein information. Nucleic Acids Res., 2005; 1: D358-63. 15. Marchler-Bauer, A., Lu, S., Anderson, J.B., Chitsaz, F., Derbyshire, M.K., DeWeese-Scott, C., Fong, J.H., Geer, L.Y., Geer, R.C., Gonzales, N.R., Gwadz, M., Hurwitz, D.I., Jackson, J.D., Ke, Z., Lanczycki, C.J., Lu, F., Marchler, G.H., Mullokandov, M., Omelchenko, M.V., Robertson, C.L., Song, J.S., Thanki, N., Yamashita, R.A., Zhang, D., Zhang, N., Zheng, C., Bryant, S.H. CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res., 2011; 39: D225–9. 16. Rzhetsky, A., Nei, M. A simple method for estimating and testing minimum evolution trees. Mol. Biol. Evol., 1992; 9: 945–967. 17. Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol., 2011; 28: 2731–2739. 18. Lepesheva, G.I., Waterman, M.R. Sterol 14±- demethylase cytochrome P450 (CYP51) a P450 in all biological kingdoms. Biochim. Biophys. Acta., 2007; 177: 467–477. 19. Picataggio, S., Rohrer, T., Deanda, K., Lanning, D., Reynolds, R., Mielenz, J., Elrich, L.D. Metabolic engineering of Candida tropicalis for the production of long-chain dicarboxylic acids. Biotechnology (N Y) 1992; 10: 894-898. 20. Takai, H., Iwama, R., Kobayashi, S., Horiuchi, H., Fukuda, R., Ohta, A. Construction and characterization of a Yarrowia lipolytica mutant lacking genes encoding cytochromes P450 subfamily 52. Fungal Genet. Biol., 2012; 49: 58-64. J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 301KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA 21. Crešnar, B., Petri, S. Cytochrome P450 enzymes in the fungal kingdom. Biochim. Biophys. Acta 2011; 1814: 29–35. 22. Chen, W., Lee, M.K., Jefcoate, C., Kim, S.C., Chen, F., Yu, J.H. Fungal cytochrome P450 monooxygenases: their distribution, structure, functions, family expansion, and evolutionary origin. Genome Biol. Evol., 2014: 6: 1620-34. 23. Eschenfeldt, W.H., Zhang, Y., Samaha, H., Stols, L., Eirich, L.D., Wilson, C.R., Donnelly, M.I. Transformation of fatty acids catalyzed by cytochrome P450 monooxygenase enzymes of Candida tropicalis. Appl. Environ. Microbiol., 2003; 69: 5992-9. 24. Goffeau, A., Barrell, B.G., Bussey, H., Davis, R.W., Dujon, B., Feldmann, H., Galibert, F., Hoheisel, J.D, Jacq, C., Johnston, M., Louis, E.J., Mewes, H.W., Murakami, Y., Philippsen, P., Tettelin, H., Oliver, S.G. Life with 6000 genes. Science 1996; 274: 563-7. 25. Engel, S.R., Dietrich, F.S., Fisk, D.G., Binkley, G., Balakrishnan, R., Costanzo, M.C., Dwight, S.S., Hitz, B.C., Karra, K., Nash, R.S., Weng, S., Wong, E.D., Lloyd, P., Skrzypek, M.S., Miyasato, S.R., Simison, M., Cherry, J.M. The reference genome sequence of Saccharomyces cerevisiae: then and now. G3 (Bethesda) 2014; 4: 389-98. 26. Kellis, M., Patterson, N., Endrizzi, M., Birren, B., Lander, E.S. Sequencing and comparison of yeast species to identify genes and regulatory elements. Nature 2003; 423: 241-54. 27. Cliften, P., Sudarsanam, P., Desikan, A., Fulton, L., Fulton, B., Majors, J., Waterson, R., Cohen, B.A., Johnston, M. Finding functional features in Saccharomyces genomes by phylogenetic footprinting. Science 2003; 301: 71-76. 28. Belloch, C., Pérez-Torrado, R., González, S.S., Pérez-Ortín, J.E., García-Martínez, J., Querol, A., Barrio, E. Chimeric genomes of natural hybrids of Saccharomyces cerevisiae and Saccharomyces kudriavzevii. Appl. Environ. Microbiol., 2009; 75: 2534–2544. 29. Jones, T., Federspiel, N.A., Chibana, H., Dungan, J., Kalman, S., Magee, B.B., Newport, G., Thorstenson, Y.R., Agabian, N., Magee, P.T., Davis, R.W., Scherer, S. The diploid genome sequence of Candida albicans. Proc. Natl. Acad. Sci. USA 2004; 101: 7329-34. 30. Braun, B.R., van Het Hoog, M., d’Enfert, C., Martchenko, M., Dungan, J., Kuo, A., Inglis, D.O., Uhl, M.A., Hogues, H., Berriman, M., Lorenz, M., Levitin, A., Oberholzer, U., Bachewich, C., Harcus, D., Marcil, A., Dignard, D., Iouk, T., Zito, R., Frangeul, L., Tekaia, F., Rutherford, K., Wang, E., Munro, C.A., Bates, S., Gow, N.A., Hoyer, L.L., Köhler, G., Morschhäuser, J., Newport, G., Znaidi, S., Raymond, M., Turcotte, B., Sherlock, G., Costanzo, M., Ihmels, J., Berman, J., Sanglard, D., Agabian, N., Mitchell, A.P., Johnson, A.D., Whiteway, M., Nantel, A.. A human-curated annotation of the Candida albicans genome. PLoS Genet., 2005; 1: 36-57. 31. Butler, G., Rasmussen, M.D., Lin, M.F., Santos, M.A., Sakthikumar, S., Munro, C.A. Evolution of pathogenicity and sexual reproduction in eight Candida genomes. Nature 2009; 459: 657-62. 32. Li, B., Guo, T., Chen, Y., Xie, J., Niu, H., Liu, D., Cheng, J., Chen, X., Wu, J., Zhuang, W., Zhu, C., Ying, H. 2014. Genome sequence of Candida tropicalis no. 121, Used for RNA Production. Genome Announc., 2014; 2, pii: e00316-14. 33. Sibirny, A.A., Boretsky, Y.R.: Pichia guilliermondii. In: Yeast biotechnology: Diversity and application (Satyanarayana T., Kunze, G, ed), Springer, Netherland, 2009; pp 113-134. 34. Gabaldón, T., Martin, T., Marcet-Houben, M., Durrens, P., Bolotin-Fukuhara, M., Lespinet, O., Arnaise, S., Boisnard, S., Aguileta, G., Atanasova, R., Bouchier, C., Couloux, A., Creno, S., Almeida Cruz, J., Devillers, H., Enache- Angoulvant, A., Guitard, J., Jaouen, L., Ma, L., Marck, C., Neuvéglise, C., Pelletier, E., Pinard, A., Poulain, J., Recoquillay, J., Westhof, E., Wincker, P., Dujon, B., Hennequin, C., Fairhead, C. Comparative genomics of emerging pathogens in the Candida glabrata clade. BMC Genomics 2013; 14: 623. 35. Jackson, A.P., Gamble, J.A., Yeomans, T., Moran, G.P., Saunders, D., Harris, D., Aslett, M., Barrell, J.F., Butler, G., Citiulo, F., Coleman, D.C., de Groot, P.W., Goodwin, T.J., Quail, M.A., McQuillan, J., Munro, C.A., Pain, A., Poulter, R.T., Rajandream, M.A., Renauld, H., Spiering, M.J., Tivey, A., Gow, N.A., Barrell, B., Sullivan, D.J., Berriman, M. Comparative genomics of the fungal pathogens Candida dubliniensis and Candida albicans. Genome Res., 2009; 19: 2231-44. 36. Wingard, J.R. Importance of Candida species other than C. albicans as pathogens in oncology patients. Clin. Infect. Dis., 1995; 20: 115-125. 37. Dias, O., Gombert, A.K., Ferreira, E.C., Rocha, I. Genome-wide metabolic (re-) annotation of Kluyveromyces lactis. BMC Genomics 2012; 13: 517. 38. Scannell, D.R., Frank, A.C., Conant, G.C., Byrne, K.P., Woolfit, M., Wolfe, K.H. Independent sorting-out of thousands of duplicated gene pairs in two yeast species J PURE APPL MICROBIO, 8(SPL. EDN.), NOVEMBER 2014. 302 KGOSIEMANG et al.: P450OME OF SACCHAROMYCOTINA descended from a whole-genome duplication. Proc. Natl. Acad. Sci. USA 2007; 104: 8397- 8402. 39. Kellis, M., Birren, B.W., Lander, E.S. Proof and evolutionary analysis of ancient genome duplication in the yeast Saccharomyces cerevisiae. Nature 2004; 428: 617-24. 40. Jeffries, T.W., Grigoriev, I.V., Grimwood, J., Laplaza, J.M., Aerts, A., Salamov, A. Genome sequence of the lignocellulose-bioconverting and xylose-fermenting yeast Pichiastipitis. Nat. Biotechnol., 2007; 25, 319-26. 41. Dietrich, F.S., Voegeli, S., Brachat, S., Lerch, A., Gates, K., Steiner, S. The Ashbyagossypii genome as a tool for mapping the ancient Saccharomyces cerevisiae genome. Science 2004; 304: 304-7. 42. Dujon, B., Sherman, D., Fischer, G., Durrens, P., Casaregola, S., Lafontaine, I. Genome evolution in yeasts. Nature 2004; 430, 35-44. 43. Fickers, P., Benetti, P.H., Waché, Y., Marty, A., Mauersberger, S, Smit, M.S., Nicaud, J.M. Hydrophobic substrate utilisation by the yeast Yarrowia lipolytica, and its potential applications. FEMS Yeast Res., 2005; 5: 527- 43. 44. Piškur, J., Ling, Z., Marcet-Houben, M., Ishchuk, O.P., Aerts, A., LaButti, K., Copeland, A., Lindquist, E., Barry, K., Compagno, C., Bisson, L., Grigoriev, I.V., Gabaldón, T., Phister, T. The genome of wine yeast Dekkera bruxellensis provides a tool to explore its food- related properties. Int. J. Food Microbiol., 2012; 157: 202-9. 45. De Schutter, K., Lin, Y.C., Tiels, P., Van Hecke, A., Glinka, S., Weber-Lehmann, J., Rouzé, P., Van de Peer, Y., Callewaert, N. Genome sequence of the recombinant protein production host Pichia pastoris. Nat. Biotech., 2009; 27: 561- 566. 46. Schneider, J., Rupp, O., Trost, E., Jaenicke, S., Passoth, V., Goesmann, A., Tauch, A., Brinkrolf, K. Genome sequence of Wickerhamomyces anomalus DSM 6766 reveals genetic basis of biotechnologically important antimicrobial activities. FEMS Yeast Res., 2012; 12: 382-386. 47. Lepesheva, G.I., Waterman, M.R. CYP51—the omnipotent P450. Mol. Cell Endocrinol., 2004; 215: 165-70. 48. Huang, F.C., Peter, A., Schwab, W. Expression and characterization of CYP52 genes involved in the biosynthesis of sophorolipid and alkane metabolism from Starmerella bombicola. Appl. Environ. Microbiol., 2014; 80: 766-76. 49. Melo, N.R., Moran, G.P., Warrilow, A.G.S., Dudley, E., Smith, S.N., Sullivan, D.J., Lamb, D.C., Kelly, D.E., Coleman, D.C., Kelly,S.L. CYP56 (Dit2p) in Candida albicans: Characterization and investigation of its role in growth and antifungal drug susceptibility. Agents Chemother., 2008; 52: 3718-3724. 50. Kelly, S.L., Lamb, D.C., Corran, A.J., Baldwin, B.C., Parks, L.W., Kelly, D.E. Purification and reconstitution of activity of Saccharomyces cerevisiae P450 61, a sterol delta 22-desaturase. FEBS Lett., 1995; 377: 217-20. 51. Mingot, J.M., Penalva, M.A. Fernandez-Canon, J.M. Disruption of phacA, an Aspergillus nidulans gene encoding a novel cytochrome P450 monooxygenase catalyzingphenylacetate 2- hydroxylation, results in penicillin over production. J. Biol. Chem. 1999; 274: 14545– 14550. 52. Ferrer-Sevillano, F., Fernandez-Canon, J.M. Novel phacB-encoded cytochrome P450 monooxygenase from Aspergillus nidulans with 3-hydroxyphenylacetate 6-hydroxylase and 3, 4-dihydroxyphenylacetate 6-hydroxylase activities. Eukaryotic Cell 2007; 6: 514–520. View publication stats