High Prevalence of bla(CTX-M) Extended Spectrum Beta-Lactamase Genes in Klebsiella pneumoniae Isolates from a Tertiary Care Hospital: First report of bla(SHV-12), bla(SHV-31), bla(SHV-38), and bla(CTX-M-15) in Brazil

The aim of this study was to investigate the presence and prevalence of bla(TEM), bla(SHV), and bla(CTX-M) and bla(GES)-like genes, responsible for extended spectrum beta-lactamases (ESBLs) production in clinical isolates of Klebsiella pneumoniae collected from a Brazilian tertiary care hospital. Sixty-five ESBL producing K. pneumoniae isolates, collected between 2005 and 2007, were screened by polymerase chain reaction (PCR). Identification of bla genes was achieved by sequencing. Genotyping of ESBL producing K. pneumoniae was performed by the enterobacterial repetitive intergenic consensus-PCR with cluster analysis by the Dice coefficient. The presence of genes encoding ESBLs was confirmed in 59/65 (90.8%) isolates, comprising 20 bla(CTX-M-2), 14 bla(CTX-M-59), 12 bla(CTX-M-15), 9 bla(SHV-12), 1 bla(SHV-2), 1 bla(SHV-2a), 1 bla(SHV-5), and 1 bla(SHV-31) genes. The ESBL genes bla(SHV-12), bla(SHV-31), and bla(CTX-M-15), and the chromosome-encoded SHV-type beta-lactamase capable of hydrolyzing imipenem were detected in Brazil for the first time. The analysis of the enterobacterial repetitive intergenic consensus-PCR band patterns revealed a high rate of multiclonal bla(CTX-M) carrying K. pneumoniae isolates (70.8%), suggesting that dissemination of encoding plasmids is likely to be the major cause of the high prevalence of these genes among the K. pneumoniae isolates considered in this study.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale

Background: Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus) microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection. Results: Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93 Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells. Conclusions: Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.

Identification and Phylogenetic Analysis of Heme Synthesis Genes in Trypanosomatids and Their Bacterial Endosymbionts

It has been known for decades that some insect-infecting trypanosomatids can survive in culture without heme supplementation while others cannot, and that this capability is associated with the presence of a betaproteobacterial endosymbiont in the flagellate's cytoplasm. However, the specific mechanisms involved in this process remained obscure. In this work, we sequence and phylogenetically analyze the heme pathway genes from the symbionts and from their hosts, as well as from a number of heme synthesis-deficient Kinetoplastida. Our results show that the enzymes responsible for synthesis of heme are encoded on the symbiont genomes and produced in close cooperation with the flagellate host. Our evidence suggests that this synergistic relationship is the end result of a history of extensive gene loss and multiple lateral gene transfer events in different branches of the phylogeny of the Trypanosomatidae.

Dissecting the Relation between a Nuclear Receptor and GATA: Binding Affinity Studies of Thyroid Hormone Receptor and GATA2 on TSH beta Promoter

Background: Much is known about how genes regulated by nuclear receptors (NRs) are switched on in the presence of a ligand. However, the molecular mechanism for gene down-regulation by liganded NRs remains a conundrum. The interaction between two zinc-finger transcription factors, Nuclear Receptor and GATA, was described almost a decade ago as a strategy adopted by the cell to up-or down-regulate gene expression. More recently, cell-based assays have shown that the Zn-finger region of GATA2 (GATA2-Zf) has an important role in down-regulation of the thyrotropin gene (TSH beta) by liganded thyroid hormone receptor (TR). Methodology/Principal Findings: In an effort to better understand the mechanism that drives TSH beta down-regulation by a liganded TR and GATA2, we have carried out equilibrium binding assays using fluorescence anisotropy to study the interaction of recombinant TR and GATA2-Zf with regulatory elements present in the TSH beta promoter. Surprisingly, we observed that ligand (T3) weakens TR binding to a negative regulatory element (NRE) present in the TSH beta promoter. We also show that TR may interact with GATA2-Zf in the absence of ligand, but T3 is crucial for increasing the affinity of this complex for different GATA response elements (GATA-REs). Importantly, these results indicate that TR complex formation enhances DNA binding of the TR-GATA2 in a ligand-dependent manner. Conclusions: Our findings extend previous results obtained in vivo, further improving our understanding of how liganded nuclear receptors down-regulate gene transcription, with the cooperative binding of transcription factors to DNA forming the core of this process.

Gene Expression Complex Networks: Synthesis, Identification, and Analysis

Thanks to recent advances in molecular biology, allied to an ever increasing amount of experimental data, the functional state of thousands of genes can now be extracted simultaneously by using methods such as cDNA microarrays and RNA-Seq. Particularly important related investigations are the modeling and identification of gene regulatory networks from expression data sets. Such a knowledge is fundamental for many applications, such as disease treatment, therapeutic intervention strategies and drugs design, as well as for planning high-throughput new experiments. Methods have been developed for gene networks modeling and identification from expression profiles. However, an important open problem regards how to validate such approaches and its results. This work presents an objective approach for validation of gene network modeling and identification which comprises the following three main aspects: (1) Artificial Gene Networks (AGNs) model generation through theoretical models of complex networks, which is used to simulate temporal expression data; (2) a computational method for gene network identification from the simulated data, which is founded on a feature selection approach where a target gene is fixed and the expression profile is observed for all other genes in order to identify a relevant subset of predictors; and (3) validation of the identified AGN-based network through comparison with the original network. The proposed framework allows several types of AGNs to be generated and used in order to simulate temporal expression data. The results of the network identification method can then be compared to the original network in order to estimate its properties and accuracy. Some of the most important theoretical models of complex networks have been assessed: the uniformly-random Erdos-Renyi (ER), the small-world Watts-Strogatz (WS), the scale-free Barabasi-Albert (BA), and geographical networks (GG). The experimental results indicate that the inference method was sensitive to average degree k variation, decreasing its network recovery rate with the increase of k. The signal size was important for the inference method to get better accuracy in the network identification rate, presenting very good results with small expression profiles. However, the adopted inference method was not sensible to recognize distinct structures of interaction among genes, presenting a similar behavior when applied to different network topologies. In summary, the proposed framework, though simple, was adequate for the validation of the inferred networks by identifying some properties of the evaluated method, which can be extended to other inference methods.

Feature selection environment for genomic applications

Background: Feature selection is a pattern recognition approach to choose important variables according to some criteria in order to distinguish or explain certain phenomena (i.e., for dimensionality reduction). There are many genomic and proteomic applications that rely on feature selection to answer questions such as selecting signature genes which are informative about some biological state, e. g., normal tissues and several types of cancer; or inferring a prediction network among elements such as genes, proteins and external stimuli. In these applications, a recurrent problem is the lack of samples to perform an adequate estimate of the joint probabilities between element states. A myriad of feature selection algorithms and criterion functions have been proposed, although it is difficult to point the best solution for each application. Results: The intent of this work is to provide an open-source multiplataform graphical environment for bioinformatics problems, which supports many feature selection algorithms, criterion functions and graphic visualization tools such as scatterplots, parallel coordinates and graphs. A feature selection approach for growing genetic networks from seed genes ( targets or predictors) is also implemented in the system. Conclusion: The proposed feature selection environment allows data analysis using several algorithms, criterion functions and graphic visualization tools. Our experiments have shown the software effectiveness in two distinct types of biological problems. Besides, the environment can be used in different pattern recognition applications, although the main concern regards bioinformatics tasks.

Characterization of global transcription profile of normal and HPV-immortalized keratinocytes and their response to TNF treatment

Background: Persistent infection by high risk HPV types (e.g. HPV-16, -18, -31, and -45) is the main risk factor for development of cervical intraepithelial neoplasia and cervical cancer. Tumor necrosis factor (TNF) is a key mediator of epithelial cell inflammatory response and exerts a potent cytostatic effect on normal or HPV16, but not on HPV18 immortalized keratinocytes. Moreover, several cervical carcinoma-derived cell lines are resistant to TNF anti-proliferative effect suggesting that the acquisition of TNF-resistance may constitute an important step in HPV-mediated carcinogenesis. In the present study, we compared the gene expression profiles of normal and HPV16 or 18 immortalized human keratinocytes before and after treatment with TNF for 3 or 60 hours. Methods: In this study, we determined the transcriptional changes 3 and 60 hours after TNF treatment of normal, HPV16 and HPV18 immortalized keratinocytes by microarray analysis. The expression pattern of two genes observed by microarray was confirmed by Northern Blot. NF-kappa B activation was also determined by electrophoretic mobility shift assay (EMSA) using specific oligonucleotides and nuclear protein extracts. Results: We observed the differential expression of a common set of genes in two TNF-sensitive cell lines that differs from those modulated in TNF-resistant ones. This information was used to define genes whose differential expression could be associated with the differential response to TNF, such as: KLK7 (kallikrein 7), SOD2 (superoxide dismutase 2), 100P (S100 calcium binding protein P), PI3 (protease inhibitor 3, skin-derived), CSTA (cystatin A), RARRES1 (retinoic acid receptor responder 1), and LXN (latexin). The differential expression of the KLK7 and SOD2 transcripts was confirmed by Northern blot. Moreover, we observed that SOD2 expression correlates with the differential NF-kappa B activation exhibited by TNF-sensitive and TNF-resistant cells. Conclusion: This is the first in depth analysis of the differential effect of TNF on normal and HPV16 or HPV18 immortalized keratinocytes. Our findings may be useful for the identification of genes involved in TNF resistance acquisition and candidate genes which deregulated expression may be associated with cervical disease establishment and/or progression.

Functional microarray analysis suggests repressed cell-cell signaling and cell survival-related modules inhibit progression of head and neck squamous cell carcinoma

Background: Cancer shows a great diversity in its clinical behavior which cannot be easily predicted using the currently available clinical or pathological markers. The identification of pathways associated with lymph node metastasis (N+) and recurrent head and neck squamous cell carcinoma (HNSCC) may increase our understanding of the complex biology of this disease. Methods: Tumor samples were obtained from untreated HNSCC patients undergoing surgery. Patients were classified according to pathologic lymph node status (positive or negative) or tumor recurrence (recurrent or non-recurrent tumor) after treatment (surgery with neck dissection followed by radiotherapy). Using microarray gene expression, we screened tumor samples according to modules comprised by genes in the same pathway or functional category. Results: The most frequent alterations were the repression of modules in negative lymph node (N0) and in non-recurrent tumors rather than induction of modules in N+ or in recurrent tumors. N0 tumors showed repression of modules that contain cell survival genes and in non-recurrent tumors cell-cell signaling and extracellular region modules were repressed. Conclusions: The repression of modules that contain cell survival genes in N0 tumors reinforces the important role that apoptosis plays in the regulation of metastasis. In addition, because tumor samples used here were not microdissected, tumor gene expression data are represented together with the stroma, which may reveal signaling between the microenvironment and tumor cells. For instance, in non-recurrent tumors, extracellular region module was repressed, indicating that the stroma and tumor cells may have fewer interactions, which disable metastasis development. Finally, the genes highlighted in our analysis can be implicated in more than one pathway or characteristic, suggesting that therapeutic approaches to prevent tumor progression should target more than one gene or pathway, specially apoptosis and interactions between tumor cells and the stroma.

Identification of protein-coding and non-coding RNA expression profiles in CD34(+) and in stromal cells in refractory anemia with ringed sideroblasts

Background: Myelodysplastic syndromes (MDS) are a group of clonal hematological disorders characterized by ineffective hematopoiesis with morphological evidence of marrow cell dysplasia resulting in peripheral blood cytopenia. Microarray technology has permitted a refined high-throughput mapping of the transcriptional activity in the human genome. Non-coding RNAs (ncRNAs) transcribed from intronic regions of genes are involved in a number of processes related to post-transcriptional control of gene expression, and in the regulation of exon-skipping and intron retention. Characterization of ncRNAs in progenitor cells and stromal cells of MDS patients could be strategic for understanding gene expression regulation in this disease. Methods: In this study, gene expression profiles of CD34(+) cells of 4 patients with MDS of refractory anemia with ringed sideroblasts (RARS) subgroup and stromal cells of 3 patients with MDS-RARS were compared with healthy individuals using 44 k combined intron-exon oligoarrays, which included probes for exons of protein-coding genes, and for non-coding RNAs transcribed from intronic regions in either the sense or antisense strands. Real-time RT-PCR was performed to confirm the expression levels of selected transcripts. Results: In CD34(+) cells of MDS-RARS patients, 216 genes were significantly differentially expressed (q-value <= 0.01) in comparison to healthy individuals, of which 65 (30%) were non-coding transcripts. In stromal cells of MDS-RARS, 12 genes were significantly differentially expressed (q-value <= 0.05) in comparison to healthy individuals, of which 3 (25%) were non-coding transcripts. Conclusions: These results demonstrated, for the first time, the differential ncRNA expression profile between MDS-RARS and healthy individuals, in CD34(+) cells and stromal cells, suggesting that ncRNAs may play an important role during the development of myelodysplastic syndromes.

Interference with Hemozoin Formation Represents an Important Mechanism of Schistosomicidal Action of Antimalarial Quinoline Methanols

Background: The parasitic trematode Schistosoma mansoni is one of the major causative agents of human schistosomiasis, which afflicts 200 million people worldwide. Praziquantel remains the main drug used for schistosomiasis treatment, and reliance on the single therapy has been prompting the search for new therapeutic compounds against this disease. Our group has demonstrated that heme crystallization into hemozoin (Hz) within the S. mansoni gut is a major heme detoxification route with lipid droplets involved in this process and acting as a potential chemotherapeutical target. In the present work, we investigated the effects of three antimalarial compounds, quinine (QN), quinidine (QND) and quinacrine (QCR) in a murine schistosomiasis model by using a combination of biochemical, cell biology and molecular biology approaches. Methodology/Principal Findings: Treatment of S. mansoni-infected female Swiss mice with daily intraperitoneal injections of QN, and QND (75 mg/kg/day) from the 11(th) to 17(th) day after infection caused significant decreases in worm burden (39%-61%) and egg production (42%-98%). Hz formation was significantly inhibited (40%-65%) in female worms recovered from QN- and QND-treated mice and correlated with reduction in the female worm burden. We also observed that QN treatment promoted remarkable ultrastructural changes in male and female worms, particularly in the gut epithelium and reduced the granulomatous reaction to parasite eggs trapped in the liver. Microarray gene expression analysis indicated that QN treatment increased the expression of transcripts related to musculature, protein synthesis and repair mechanisms. Conclusions: The overall significant reduction in several disease burden parameters by the antimalarial quinoline methanols indicates that interference with Hz formation in S. mansoni represents an important mechanism of schistosomicidal action of these compounds and points out the heme crystallization process as a valid chemotherapeutic target to treat schistosomiasis.

Haplotype distribution of five nuclear genes based on network genealogies and Bayesian inference indicates that Trypanosoma cruzi hybrid strains are polyphyletic

Chagas disease is still a major public health problem in Latin America. Its causative agent, Trypanosoma cruzi, can be typed into three major groups, T. cruzi I, T. cruzi II and hybrids. These groups each have specific genetic characteristics and epidemiological distributions. Several highly virulent strains are found in the hybrid group; their origin is still a matter of debate. The null hypothesis is that the hybrids are of polyphyletic origin, evolving independently from various hybridization events. The alternative hypothesis is that all extant hybrid strains originated from a single hybridization event. We sequenced both alleles of genes encoding EF-1 alpha, actin and SSU rDNA of 26 T. cruzi strains and DHFR-TS and TR of 12 strains. This information was used for network genealogy analysis and Bayesian phylogenies. We found T. cruzi I and T. cruzi II to be monophyletic and that all hybrids had different combinations of T. cruzi I and T. cruzi II haplotypes plus hybrid-specific haplotypes. Bootstrap values (networks) and posterior probabilities (Bayesian phylogenies) of clades supporting the monophyly of hybrids were far below the 95% confidence interval, indicating that the hybrid group is polyphyletic. We hypothesize that T. cruzi I and T. cruzi II are two different species and that the hybrids are extant representatives of independent events of genome hybridization, which sporadically have sufficient fitness to impact on the epidemiology of Chagas disease.

Intercellular transport of epidermis-expressed MADS domain transcription factors and their effect on plant morphology and floral transition

P>During the lifetime of an angiosperm plant various important processes such as floral transition, specification of floral organ identity and floral determinacy, are controlled by members of the MADS domain transcription factor family. To investigate the possible non-cell-autonomous function of MADS domain proteins, we expressed GFP-tagged clones of AGAMOUS (AG), APETALA3 (AP3), PISTILLATA (PI) and SEPALLATA3 (SEP3) under the control of the MERISTEMLAYER1 promoter in Arabidopsis thaliana plants. Morphological analyses revealed that epidermal overexpression was sufficient for homeotic changes in floral organs, but that it did not result in early flowering or terminal flower phenotypes that are associated with constitutive overexpression of these proteins. Localisations of the tagged proteins in these plants were analysed with confocal laser scanning microscopy in leaf tissue, inflorescence meristems and floral meristems. We demonstrated that only AG is able to move via secondary plasmodesmata from the epidermal cell layer to the subepidermal cell layer in the floral meristem and to a lesser extent in the inflorescence meristem. To study the homeotic effects in more detail, the capacity of trafficking AG to complement the ag mutant phenotype was compared with the capacity of the non-inwards-moving AP3 protein to complement the ap3 mutant phenotype. While epidermal expression of AG gave full complementation, AP3 appeared not to be able to drive all homeotic functions from the epidermis, perhaps reflecting the difference in mobility of these proteins.

The Influence of Different Land Uses on the Structure of Archaeal Communities in Amazonian Anthrosols Based on 16S rRNA and amoA Genes

Soil from the Amazonian region is usually regarded as unsuitable for agriculture because of its low organic matter content and low pH; however, this region also contains extremely rich soil, the Terra Preta Anthrosol. A diverse archaeal community usually inhabits acidic soils, such as those found in the Amazon. Therefore, we hypothesized that this community should be sensitive to changes in the environment. Here, the archaeal community composition of Terra Preta and adjacent soil was examined in four different sites in the Brazilian Amazon under different anthropic activities. The canonical correspondence analysis of terminal restriction fragment length polymorphisms has shown that the archaeal community structure was mostly influenced by soil attributes that differentiate the Terra Preta from the adjacent soil (i.e., pH, sulfur, and organic matter). Archaeal 16S rRNA gene clone libraries indicated that the two most abundant genera in both soils were Candidatus nitrosphaera and Canditatus nitrosocaldus. An ammonia monoxygenase gene (amoA) clone library analysis indicated that, within each site, there was no significant difference between the clone libraries of Terra Preta and adjacent soils. However, these clone libraries indicated there were significant differences between sites. Quantitative PCR has shown that Terra Preta soils subjected to agriculture displayed a higher number of amoA gene copy numbers than in adjacent soils. On the other hand, soils that were not subjected to agriculture did not display significant differences on amoA gene copy numbers between Terra Preta and adjacent soils. Taken together, our findings indicate that the overall archaeal community structure in these Amazonian soils is determined by the soil type and the current land use.

Searching for Moniliophthora perniciosa pathogenicity genes

The basidiomycete Moniliophthora perniciosa is the causal agent of witches` broom disease of Theobroma cacao (cacao). Pathogenesis mechanisms of this hemibiotrophic fungus are largely unknown. An approach to identify putative pathogenicity genes is searching for sequences induced in mycelia grown under in vitro conditions. Using this approach, genes from M. perniciosa induced under limiting nitrogen and light were identified from a cDNA library enriched by suppression subtractive hybridization as potential putative pathogenicity genes. From the 159 identified unique sequences, 59 were annotated and classified by gene ontology. Two sequences were categorized as ""Defence genes, Virulence, and Cell response"" presumably coding for allergenic proteins, whose homologues from other fungi are inducers of animal or plant defences. Differential gene expression was evaluated by quantitative amplification of reversed transcripts (RT-qPCR) of the putative identified genes coding for the two allergenic proteins (Aspf13 and 88KD), and for the enzymes Arylsulfatase (AS); Aryl-Alcohol Oxidase; Aldo-Keto Reductase (AK); Cytochrome P450 (P450); Phenylalanine Ammonia-Lyase; and Peroxidase from mycelia grown under contrasting N concentrations. All genes were validated for differential expression, except for the putative Peroxidase. The same eight genes were analysed for expression in susceptible plants inoculated with M. perniciosa, and six were induced during the early asymptomatic stage of the disease. In infected host tissues, transcripts of 88KD and AS were found more abundant at the biotrophic phase, while those from Aspf13, AK, PAL, and P450 accumulated at the necrotrophic phase, enabling to suggest that mycelia transition from biotrophic to necrotrophic might occur earlier than currently considered. These sequences appeared to be virulence life-style genes, which encode factors or enzymes that enable invasion, colonization or intracellular survival, or manipulate host factors to benefit the pathogen`s own survival in the hostile environment. (C) 2010 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.