Enhancing Bioinformatics and Genomics Courses: Building Capacity and Skills via Lab Meeting Activities: Fostering a Culture of Critical Capacities to Read, Write, Communicate and Engage in Rigorous Scientific Exchanges.

Reading, writing, publishing, and publicly presenting scientific works are vital for a young researcher's profile building and career development. Generally, the traditional educational curricula do not offer training possibilities to learn and practice how to prepare, write, and present scientific works. These are rather a part of lab meeting activities in research groups. The lack of such training is more critical in some developing countries because this adds to the rare opportunities to discuss and become involved in the exchanges on state of the art scientific literature. Here the authors relate their experience in introducing a weekly 1-day lab meeting in the framework of two previously organized 3-month courses on "Bioinformatics and Genome Analyses". The main activities which are developed during these lab meetings include scientific literature follow up as well as preparing and presenting oral and written scientific reviews. These activities prove to be useful for a student's self-confidence building, for enhancing their active participation during the lectures and practical sessions, as well as for the positive impact on running the whole course program. Incorporation of such lab meeting activities in the course program significantly improves the capacity building of the participants, their analytical and critical reading of scientific literature, as well as communication skills. In this work it is shown how to proceed with the different steps involved in the implementation of lab meeting activities, and to recommend their regular institution in similar courses.

Designing and running an advanced Bioinformatics and genome analyses course in Tunisia.

Genome data, with underlying new knowledge, are accumulating at exponential rate thanks to ever-improving sequencing technologies and the parallel development of dedicated efficient Bioinformatics methods and tools. Advanced Education in Bioinformatics and Genome Analyses is to a large extent not accessible to students in developing countries where endeavors to set up Bioinformatics courses concern most often only basic levels. Here, we report a pioneering pilot experience concerning the design and implementation, from scratch, of a three-months advanced and extensive course in Bioinformatics and Genome Analyses in the Institut Pasteur de Tunis. Most significantly the outcome of the course was upgrading the participants' skills in Bioinformatics and Genome Analyses to recognized international standards. Here we detail the different steps involved in the implementation of this course as well as the topics covered in the program. The description of this pilot experience might be helpful for the implementation of other similar educational projects, notably in developing countries, aiming to go beyond basics and providing young researchers with high-level skills.

Genome Data Exploration Using Correspondence Analysis.

Recent developments of sequencing technologies that allow the production of massive amounts of genomic and genotyping data have highlighted the need for synthetic data representation and pattern recognition methods that can mine and help discovering biologically meaningful knowledge included in such large data sets. Correspondence analysis (CA) is an exploratory descriptive method designed to analyze two-way data tables, including some measure of association between rows and columns. It constructs linear combinations of variables, known as factors. CA has been used for decades to study high-dimensional data, and remarkable inferences from large data tables were obtained by reducing the dimensionality to a few orthogonal factors that correspond to the largest amount of variability in the data. Herein, I review CA and highlight its use by considering examples in handling high-dimensional data that can be constructed from genomic and genetic studies. Examples in amino acid compositions of large sets of species (viruses, phages, yeast, and fungi) as well as an example related to pairwise shared orthologs in a set of yeast and fungal species, as obtained from their proteome comparisons, are considered. For the first time, results show striking segregations between yeasts and fungi as well as between viruses and phages. Distributions obtained from shared orthologs show clusters of yeast and fungal species corresponding to their phylogenetic relationships. A direct comparison with the principal component analysis method is discussed using a recently published example of genotyping data related to newly discovered traces of an ancient hominid that was compared to modern human populations in the search for ancestral similarities. CA offers more detailed results highlighting links between modern humans and the ancient hominid and their characterizations. Compared to the popular principal component analysis method, CA allows easier and more effective interpretation of results, particularly by the ability of relating individual patterns with their corresponding characteristic variables.

Inferring Orthologs: Open Questions and Perspectives.

With the increasing number of sequenced genomes and their comparisons, the detection of orthologs is crucial for reliable functional annotation and evolutionary analyses of genes and species. Yet, the dynamic remodeling of genome content through gain, loss, transfer of genes, and segmental and whole-genome duplication hinders reliable orthology detection. Moreover, the lack of direct functional evidence and the questionable quality of some available genome sequences and annotations present additional difficulties to assess orthology. This article reviews the existing computational methods and their potential accuracy in the high-throughput era of genome sequencing and anticipates open questions in terms of methodology, reliability, and computation. Appropriate taxon sampling together with combination of methods based on similarity, phylogeny, synteny, and evolutionary knowledge that may help detecting speciation events appears to be the most accurate strategy. This review also raises perspectives on the potential determination of orthology throughout the whole species phylogeny.

Complete DNA sequence of Kuraishia capsulata illustrates novel genomic features among budding yeasts (Saccharomycotina).

The numerous yeast genome sequences presently available provide a rich source of information for functional as well as evolutionary genomics but unequally cover the large phylogenetic diversity of extant yeasts. We present here the complete sequence of the nuclear genome of the haploid-type strain of Kuraishia capsulata (CBS1993(T)), a nitrate-assimilating Saccharomycetales of uncertain taxonomy, isolated from tunnels of insect larvae underneath coniferous barks and characterized by its copious production of extracellular polysaccharides. The sequence is composed of seven scaffolds, one per chromosome, totaling 11.4 Mb and containing 6,029 protein-coding genes, ~13.5% of which being interrupted by introns. This GC-rich yeast genome (45.7%) appears phylogenetically related with the few other nitrate-assimilating yeasts sequenced so far, Ogataea polymorpha, O. parapolymorpha, and Dekkera bruxellensis, with which it shares a very reduced number of tRNA genes, a novel tRNA sparing strategy, and a common nitrate assimilation cluster, three specific features to this group of yeasts. Centromeres were recognized in GC-poor troughs of each scaffold. The strain bears MAT alpha genes at a single MAT locus and presents a significant degree of conservation with Saccharomyces cerevisiae genes, suggesting that it can perform sexual cycles in nature, although genes involved in meiosis were not all recognized. The complete absence of conservation of synteny between K. capsulata and any other yeast genome described so far, including the three other nitrate-assimilating species, validates the interest of this species for long-range evolutionary genomic studies among Saccharomycotina yeasts.

Detection and characterization of megasatellites in orthologous and nonorthologous genes of 21 fungal genomes.

Megasatellites are large DNA tandem repeats, originally described in Candida glabrata, in protein-coding genes. Most of the genes in which megasatellites are found are of unknown function. In this work, we extended the search for megasatellites to 20 additional completely sequenced fungal genomes and extracted 216 megasatellites in 203 out of 142,121 genes, corresponding to the most exhaustive description of such genetic elements available today. We show that half of the megasatellites detected encode threonine-rich peptides predicted to be intrinsically disordered, suggesting that they may interact with several partners or serve as flexible linkers. Megasatellite motifs were clustered into several families. Their distribution in fungal genes shows that different motifs are found in orthologous genes and similar motifs are found in unrelated genes, suggesting that megasatellite formation or spreading does not necessarily track the evolution of their host genes. Altogether, these results suggest that megasatellites are created and lost during evolution of fungal genomes, probably sharing similar functions, although their primary sequences are not necessarily conserved.

Investigation of secreted protein transcripts as early biomarkers for type 1 diabetes in the mouse model.

Type 1 diabetes (T1D) represents a serious health burden in the world, complicated by the fact that disease onset can be preceded by a long time period without evident clinical signs. It would be then of critical importance to detect the disease in its early stages. In this direction, we seek here to identify early preinflammatory markers for autoimmune diabetes, mining our previously reported transcriptome data relevant to distinct early sub-phenotypes in the NOD mouse, associated with early insulin autoantibodies (E-IAA). More specifically we focus on secreted or transmembrane protein transcripts, identifying in this category 71 differentially expressed transcripts which are regulated at the early preinflammatory stages of T1D in the pancreatic lymph nodes (PLN). Following the expression patterns of these 71 transcripts, correspondence analysis (a multivariate analysis method) reveals a clear-cut segregation of the individual samples according to the early subphenotype used. Thus the 71 transcripts coding for secreted proteins constitute a candidate-set of predictive biomarkers for the development of autoimmune damage of the ß cells of the pancreas. The majority of these genes have human orthologs and accordingly they represent potential candidate biomarkers for the human disease. In addition, for predictive purposes, the analysis reveals the possibility to reduce significantly the size of the candidate-set in practice, with various genes displaying identical expression profiles.

SuperPartitions: detection and classification of orthologs.

The proper detection of orthologs is crucial for evolutionary studies of genes and species. Despite large efforts to solve this problem the methodological situation appears unsettled to a large extent and the "quest for orthologs" is still an ongoing task in large-scale genome comparisons. Here, we introduce a simple operational framework for the detection of orthologs and their classification. The operational framework relies on well-established principles, optimizing their implementation for the considered purposes, and chaining components in coherent procedures: 1) We take advantage of the efficiency and simplicity of the Reciprocal Best Hit (RBH) detections, remedying (by design) the drawback concerning the limitations in terms of 1:1 detections. The procedure is based on the partitioning of Reciprocal Best Hits, with the further merging of partitions including members of the same paralogous classes ("SuperPartition of Orthologs" (SPOs)). 2) We then resort to the conservation profiles of the obtained clusters, allowing simple detection of SPOs containing duplicated members. Based on accepted evolutionary principles, such members can be further tagged as in-paralogs (co-orthologs) or out-paralogs. The method is illustrated and validated by extensive genomic analyses. The performances of the overall approach are characterized in global terms for three sets of species (Chlamydiae, Mycobacteria, Aspergilli), showing that at least 75% of the sets of orthologs contain at most one protein from a given species. The sets including more than one protein from a given species are shown to contain in-paralogs in proportions varying from 28% to 58%. The characterizations also show that the large majority of SPOs are associated with ancestral motifs, and accordingly not prone to chaining effects that might be triggered by multi-domain proteins. Further the SPO formulation is compared to other similarity based ortholog detection methods. Beyond core common results, significant differences are observed between various methods, which can be accounted for to a large extent on conceptual grounds, relative to the different merging schemes involved. Such comparisons highlight a major advantage of the SPO approach concerning the proper clustering of associated paralogs, which appear to be often dispatched spuriously into distinct orthologous classes. Finally the perspectives for future applications and elaborations of SPO-based compositional analyses are discussed.

Use of fine-needle aspiration for diagnosis of Mycobacterium ulcerans infection.

Noninvasive methods for the bacteriological diagnosis of early-stage Mycobacterium ulcerans infection are not available. It was recently shown that fine-needle aspiration (FNA) could be used for diagnosing M. ulcerans infection in ulcerative lesions. We report that FNA is an appropriate sampling method for diagnosing M. ulcerans infection in nonulcerative lesions.

Promiscuous DNA in the nuclear genomes of hemiascomycetous yeasts.

Abstract Transfer of fragments of mtDNA to the nuclear genome is a general phenomenon that gives rise to NUMTs (NUclear sequences of MiTochondrial origin). We present here the first comparative analysis of the NUMT content of entirely sequenced species belonging to a monophyletic group, the hemiascomycetous yeasts (Candida glabrata, Kluyveromyces lactis, Kluyveromyces thermotolerans, Debaryomyces hansenii and Yarrowia lipolytica, along with the updated NUMT content of Saccharomyces cerevisiae). This study revealed a huge diversity in NUMT number and organization across the six species. Debaryomyces hansenii harbors the highest number of NUMTs (145), half of which are distributed in numerous large mosaics of up to eight NUMTs arising from multiple noncontiguous mtDNA fragments inserted at the same chromosomal locus. Most NUMTs, in all species, are found within intergenic regions including seven NUMTs in pseudogenes. However, five NUMTs overlap a gene, suggesting a positive impact of NUMTs on protein evolution. Contrary to the other species, K. lactis and K. thermotolerans harbor only a few diverged NUMTs, suggesting that mitochondrial transfer to the nuclear genome has decreased or ceased in these phylogenetic branches. The dynamics of NUMT acquisition and loss are illustrated here by their species-specific distribution.

Protection against Mycobacterium ulcerans lesion development by exposure to aquatic insect saliva.

BACKGROUND: Buruli ulcer is a severe human skin disease caused by Mycobacterium ulcerans. This disease is primarily diagnosed in West Africa with increasing incidence. Antimycobacterial drug therapy is relatively effective during the preulcerative stage of the disease, but surgical excision of lesions with skin grafting is often the ultimate treatment. The mode of transmission of this Mycobacterium species remains a matter of debate, and relevant interventions to prevent this disease lack (i) the proper understanding of the M. ulcerans life history traits in its natural aquatic ecosystem and (ii) immune signatures that could be correlates of protection. We previously set up a laboratory ecosystem with predatory aquatic insects of the family Naucoridae and laboratory mice and showed that (i) M. ulcerans-carrying aquatic insects can transmit the mycobacterium through bites and (ii) that their salivary glands are the only tissues hosting replicative M. ulcerans. Further investigation in natural settings revealed that 5%-10% of these aquatic insects captured in endemic areas have M. ulcerans-loaded salivary glands. In search of novel epidemiological features we noticed that individuals working close to aquatic environments inhabited by insect predators were less prone to developing Buruli ulcers than their relatives. Thus we set out to investigate whether those individuals might display any immune signatures of exposure to M. ulcerans-free insect predator bites, and whether those could correlate with protection. METHODS AND FINDINGS: We took a two-pronged approach in this study, first investigating whether the insect bites are protective in a mouse model, and subsequently looking for possibly protective immune signatures in humans. We found that, in contrast to control BALB/c mice, BALB/c mice exposed to Naucoris aquatic insect bites or sensitized to Naucoris salivary gland homogenates (SGHs) displayed no lesion at the site of inoculation of M. ulcerans coated with Naucoris SGH components. Then using human serum samples collected in a Buruli ulcer-endemic area (in the Republic of Benin, West Africa), we assayed sera collected from either ulcer-free individuals or patients with Buruli ulcers for the titre of IgGs that bind to insect predator SGH, focusing on those molecules otherwise shown to be retained by M. ulcerans colonies. IgG titres were lower in the Buruli ulcer patient group than in the ulcer-free group. CONCLUSIONS: These data will help structure future investigations in Buruli ulcer-endemic areas, providing a rationale for research into human immune signatures of exposure to predatory aquatic insects, with special attention to those insect saliva molecules that bind to M. ulcerans.

Reductive evolution and niche adaptation inferred from the genome of Mycobacterium ulcerans, the causative agent of Buruli ulcer.

Mycobacterium ulcerans is found in aquatic ecosystems and causes Buruli ulcer in humans, a neglected but devastating necrotic disease of subcutaneous tissue that is rampant throughout West and Central Africa. Here, we report the complete 5.8-Mb genome sequence of M. ulcerans and show that it comprises two circular replicons, a chromosome of 5632 kb and a virulence plasmid of 174 kb. The plasmid is required for production of the polyketide toxin mycolactone, which provokes necrosis. Comparisons with the recently completed 6.6-Mb genome of Mycobacterium marinum revealed >98% nucleotide sequence identity and genome-wide synteny. However, as well as the plasmid, M. ulcerans has accumulated 213 copies of the insertion sequence IS2404, 91 copies of IS2606, 771 pseudogenes, two bacteriophages, and multiple DNA deletions and rearrangements. These data indicate that M. ulcerans has recently evolved via lateral gene transfer and reductive evolution from the generalist, more rapid-growing environmental species M. marinum to become a niche-adapted specialist. Predictions based on genome inspection for the production of modified mycobacterial virulence factors, such as the highly abundant phthiodiolone lipids, were confirmed by structural analyses. Similarly, 11 protein-coding sequences identified as M. ulcerans-specific by comparative genomics were verified as such by PCR screening a diverse collection of 33 strains of M. ulcerans and M. marinum. This work offers significant insight into the biology and evolution of mycobacterial pathogens and is an important component of international efforts to counter Buruli ulcer.

Evolution of proteomes: fundamental signatures and global trends in amino acid compositions.

BACKGROUND: The evolutionary characterization of species and lifestyles at global levels is nowadays a subject of considerable interest, particularly with the availability of many complete genomes. Are there specific properties associated with lifestyles and phylogenies? What are the underlying evolutionary trends? One of the simplest analyses to address such questions concerns characterization of proteomes at the amino acids composition level. RESULTS: In this work, amino acid compositions of a large set of 208 proteomes, with significant number of representatives from the three phylogenetic domains and different lifestyles are analyzed, resorting to an appropriate multidimensional method: Correspondence analysis. The analysis reveals striking discrimination between eukaryotes, prokaryotic mesophiles and hyperthemophiles-themophiles, following amino acid usage. In sharp contrast, no similar discrimination is observed for psychrophiles. The observed distributional properties are compared with various inferred chronologies for the recruitment of amino acids into the genetic code. Such comparisons reveal correlations between the observed segregations of species following amino acid usage, and the separation of amino acids following early or late recruitment. CONCLUSION: A simple description of proteomes according to amino acid compositions reveals striking signatures, with sharp segregations or on the contrary non-discriminations following phylogenies and lifestyles. The distribution of species, following amino acid usage, exhibits a discrimination between [high GC]-[high optimal growth temperatures] and [low GC]-[moderate temperatures] characteristics. This discrimination appears to coincide closely with the separation of amino acids following their inferred early or late recruitment into the genetic code. Taken together the various results provide a consistent picture for the evolution of proteomes, in terms of amino acid usage.

Genome trees from conservation profiles.

The concept of the genome tree depends on the potential evolutionary significance in the clustering of species according to similarities in the gene content of their genomes. In this respect, genome trees have often been identified with species trees. With the rapid expansion of genome sequence data it becomes of increasing importance to develop accurate methods for grasping global trends for the phylogenetic signals that mutually link the various genomes. We therefore derive here the methodological concept of genome trees based on protein conservation profiles in multiple species. The basic idea in this derivation is that the multi-component "presence-absence" protein conservation profiles permit tracking of common evolutionary histories of genes across multiple genomes. We show that a significant reduction in informational redundancy is achieved by considering only the subset of distinct conservation profiles. Beyond these basic ideas, we point out various pitfalls and limitations associated with the data handling, paving the way for further improvements. As an illustration for the methods, we analyze a genome tree based on the above principles, along with a series of other trees derived from the same data and based on pair-wise comparisons (ancestral duplication-conservation and shared orthologs). In all trees we observe a sharp discrimination between the three primary domains of life: Bacteria, Archaea, and Eukarya. The new genome tree, based on conservation profiles, displays a significant correspondence with classically recognized taxonomical groupings, along with a series of departures from such conventional clusterings.

A human-curated annotation of the Candida albicans genome.

Recent sequencing and assembly of the genome for the fungal pathogen Candida albicans used simple automated procedures for the identification of putative genes. We have reviewed the entire assembly, both by hand and with additional bioinformatic resources, to accurately map and describe 6,354 genes and to identify 246 genes whose original database entries contained sequencing errors (or possibly mutations) that affect their reading frame. Comparison with other fungal genomes permitted the identification of numerous fungus-specific genes that might be targeted for antifungal therapy. We also observed that, compared to other fungi, the protein-coding sequences in the C. albicans genome are especially rich in short sequence repeats. Finally, our improved annotation permitted a detailed analysis of several multigene families, and comparative genomic studies showed that C. albicans has a far greater catabolic range, encoding respiratory Complex 1, several novel oxidoreductases and ketone body degrading enzymes, malonyl-CoA and enoyl-CoA carriers, several novel amino acid degrading enzymes, a variety of secreted catabolic lipases and proteases, and numerous transporters to assimilate the resulting nutrients. The results of these efforts will ensure that the Candida research community has uniform and comprehensive genomic information for medical research as well as for future diagnostic and therapeutic applications.

Aspergillus fumigatus: saprophyte or pathogen?

Large-scale genome comparisons have shown that no gene sets are shared exclusively by both Aspergillus fumigatus and any other human pathogen sequenced to date, such as Candida or Cryptococcus species. By contrast, and in agreement with the environmental occurrence of this fungus in decaying vegetation, the enzymatic machinery required by a fungus to colonize plant substrates has been found in the A. fumigatus genome. In addition, the proteome of this fungus contains numerous efflux pumps, including >100 major facilitators that help the fungus to resist either natural aggressive molecules present in the environment or antifungal drugs in humans. Environment sensing, counteracting reactive oxidants, and retrieving essential nutriments from the environment are general metabolic traits that are associated with the growth of the saprotrophic mold A. fumigatus in an unfriendly environment such as its human host.

Continued colonization of the human genome by mitochondrial DNA.

Integration of mitochondrial DNA fragments into nuclear chromosomes (giving rise to nuclear DNA sequences of mitochondrial origin, or NUMTs) is an ongoing process that shapes nuclear genomes. In yeast this process depends on double-strand-break repair. Since NUMTs lack amplification and specific integration mechanisms, they represent the prototype of exogenous insertions in the nucleus. From sequence analysis of the genome of Homo sapiens, followed by sampling humans from different ethnic backgrounds, and chimpanzees, we have identified 27 NUMTs that are specific to humans and must have colonized human chromosomes in the last 4-6 million years. Thus, we measured the fixation rate of NUMTs in the human genome. Six such NUMTs show insertion polymorphism and provide a useful set of DNA markers for human population genetics. We also found that during recent human evolution, Chromosomes 18 and Y have been more susceptible to colonization by NUMTs. Surprisingly, 23 out of 27 human-specific NUMTs are inserted in known or predicted genes, mainly in introns. Some individuals carry a NUMT insertion in a tumor-suppressor gene and in a putative angiogenesis inhibitor. Therefore in humans, but not in yeast, NUMT integrations preferentially target coding or regulatory sequences. This is indeed the case for novel insertions associated with human diseases and those driven by environmental insults. We thus propose a mutagenic phenomenon that may be responsible for a variety of genetic diseases in humans and suggest that genetic or environmental factors that increase the frequency of chromosome breaks provide the impetus for the continued colonization of the human genome by mitochondrial DNA.

Novel transporters from hemiascomycete yeasts.

We screened nearly one thousand random sequenced targets obtained by partial sequencing of 13 hemiascomycete genomes identified by higher amino acid sequence similarity to a non-Saccharomyces cerevisiae protein than to a S. Cerevisiae protein. Among those sequences we have identified 36 novel phylogenetic clusters of putative transporters which, according to the Transport Commission system (TC-DB, 2002; http:// tcdb.ucsd.edu/tcdb), do not belong to acknowledged S. Cerevisiae protein families [De Hertogh et al.: Funct. Integr. Genomics 2002;2:154-170; http://cbi.labri.u-bordeaux.fr/Genolevures]. These novel hemiascomycete transporters comprise 3 channels, 23 secondary transporters, 5 primary transporters and 5 membrane proteins of unknown function.

A novel design of whole-genome microarray probes for Saccharomyces cerevisiae which minimizes cross-hybridization.

BACKGROUND: Numerous DNA microarray hybridization experiments have been performed in yeast over the last years using either synthetic oligonucleotides or PCR-amplified coding sequences as probes. The design and quality of the microarray probes are of critical importance for hybridization experiments as well as subsequent analysis of the data. RESULTS: We present here a novel design of Saccharomyces cerevisiae microarrays based on a refined annotation of the genome and with the aim of reducing cross-hybridization between related sequences. An effort was made to design probes of similar lengths, preferably located in the 3'-end of reading frames. The sequence of each gene was compared against the entire yeast genome and optimal sub-segments giving no predicted cross-hybridization were selected. A total of 5660 novel probes (more than 97% of the yeast genes) were designed. For the remaining 143 genes, cross-hybridization was unavoidable. Using a set of 18 deletant strains, we have experimentally validated our cross-hybridization procedure. Sensitivity, reproducibility and dynamic range of these new microarrays have been measured. Based on this experience, we have written a novel program to design long oligonucleotides for microarray hybridizations of complete genome sequences. CONCLUSIONS: A validated procedure to predict cross-hybridization in microarray probe design was defined in this work. Subsequently, a novel Saccharomyces cerevisiae microarray (which minimizes cross-hybridization) was designed and constructed. Arrays are available at Eurogentec S. A. Finally, we propose a novel design program, OliD, which allows automatic oligonucleotide design for microarrays. The OliD program is available from authors.

Expressed sequence tag analysis of the human pathogen Paracoccidioides brasiliensis yeast phase: identification of putative homologues of Candida albicans virulence and pathogenicity genes.

Paracoccidioides brasiliensis, a thermodimorphic fungus, is the causative agent of the prevalent systemic mycosis in Latin America, paracoccidioidomycosis. We present here a survey of expressed genes in the yeast pathogenic phase of P. brasiliensis. We obtained 13,490 expressed sequence tags from both 5' and 3' ends. Clustering analysis yielded the partial sequences of 4,692 expressed genes that were functionally classified by similarity to known genes. We have identified several Candida albicans virulence and pathogenicity homologues in P. brasiliensis. Furthermore, we have analyzed the expression of some of these genes during the dimorphic yeast-mycelium-yeast transition by real-time quantitative reverse transcription-PCR. Clustering analysis of the mycelium-yeast transition revealed three groups: (i) RBT, hydrophobin, and isocitrate lyase; (ii) malate dehydrogenase, contigs Pb1067 and Pb1145, GPI, and alternative oxidase; and (iii) ubiquitin, delta-9-desaturase, HSP70, HSP82, and HSP104. The first two groups displayed high mRNA expression in the mycelial phase, whereas the third group showed higher mRNA expression in the yeast phase. Our results suggest the possible conservation of pathogenicity and virulence mechanisms among fungi, expand considerably gene identification in P. brasiliensis, and provide a broader basis for further progress in understanding its biological peculiarities.