Marcin Grynberg

Comparative genomic analyses of the human fungal pathogens Coccidioides and their relatives.

While most Ascomycetes tend to associate principally with plants, the dimorphic fungi Coccidioides immitis and Coccidioides posadasii are primary pathogens of immunocompetent mammals, including humans. Infection results from environmental exposure to Coccidiodies, which is believed to grow as a soil saprophyte in arid deserts. To investigate hypotheses about the life history and evolution of Coccidioides, the genomes of several Onygenales, including C. immitis and C. posadasii; a close, nonpathogenic relative, Uncinocarpus reesii; and a more diverged pathogenic fungus, Histoplasma capsulatum, were sequenced and compared with those of 13 more distantly related Ascomycetes. This analysis identified increases and decreases in gene family size associated with a host/substrate shift from plants to animals in the Onygenales. In addition, comparison among Onygenales genomes revealed evolutionary changes in Coccidioides that may underlie its infectious phenotype, the identification of which may facilitate improved treatment and prevention of coccidioidomycosis. Overall, the results suggest that Coccidioides species are not soil saprophytes, but that they have evolved to remain associated with their dead animal hosts in soil, and that Coccidioides metabolism genes, membrane-related proteins, and putatively antigenic compounds have evolved in response to interaction with an animal host.

The retinitis pigmentosa GTPase regulator (RPGR)- interacting protein: Subserving RPGR function and participating in disk morphogenesis

Retinitis pigmentosa is a photoreceptor degenerative disease leading to blindness in adulthood. Leber congenital amaurosis (LCA) describes a more severe condition with visual deficit in early childhood. Defects in the retinitis pigmentosa GTPase regulator (RPGR) and an RPGR-interacting protein (RPGRIP) are known causes of retinitis pigmentosa and LCA, respectively. Both proteins localize in the photoreceptor connecting cilium (CC), a thin bridge linking the cell body and the light-sensing outer segment. We show that RPGR is absent in the CC of photoreceptors lacking RPGRIP, but not vice versa. Mice lacking RPGRIP elaborate grossly oversized outer segment disks resembling a cytochalasin D-induced defect and have a more severe disease than mice lacking RPGR. Mice lacking both proteins are phenotypically indistinguishable from mice lacking RPGRIP alone. In vitro, RPGRIP forms homodimer and elongated filaments via interactions involving its coiled-coil and C-terminal domains. We conclude that RPGRIP is a stable polymer in the CC where it tethers RPGR and that RPGR depends on RPGRIP for subcellular localization and normal function. Our data suggest that RPGRIP is also required for disk morphogenesis, putatively by regulating actin cytoskeleton dynamics. The latter hypothesis may be consistent with a distant homology between the C-terminal domain of RPGRIP and an actin-fragmin kinase, predicted by fold recognition algorithms. A defect in RPGRIP encompasses loss of both functions, hence the more severe clinical manifestation as LCA.

Comparative Genomic Analysis of Human Fungal Pathogens Causing Paracoccidioidomycosis

Paracoccidioides is a fungal pathogen and the cause of paracoccidioidomycosis, a health-threatening human systemic mycosis endemic to Latin America. Infection by Paracoccidioides, a dimorphic fungus in the order Onygenales, is coupled with a thermally regulated transition from a soil-dwelling filamentous form to a yeast-like pathogenic form. To better understand the genetic basis of growth and pathogenicity in Paracoccidioides, we sequenced the genomes of two strains of Paracoccidioides brasiliensis (Pb03 and Pb18) and one strain of Paracoccidioides lutzii (Pb01). These genomes range in size from 29.1 Mb to 32.9 Mb and encode 7,610 to 8,130 genes. To enable genetic studies, we mapped 94% of the P. brasiliensis Pb18 assembly onto five chromosomes. We characterized gene family content across Onygenales and related fungi, and within Paracoccidioides we found expansions of the fungal-specific kinase family FunK1. Additionally, the Onygenales have lost many genes involved in carbohydrate metabolism and fewer genes involved in protein metabolism, resulting in a higher ratio of proteases to carbohydrate active enzymes in the Onygenales than their relatives. To determine if gene content correlated with growth on different substrates, we screened the non-pathogenic onygenale Uncinocarpus reesii, which has orthologs for 91% of Paracoccidioides metabolic genes, for growth on 190 carbon sources. U. reesii showed growth on a limited range of carbohydrates, primarily basic plant sugars and cell wall components; this suggests that Onygenales, including dimorphic fungi, can degrade cellulosic plant material in the soil. In addition, U. reesii grew on gelatin and a wide range of dipeptides and amino acids, indicating a preference for proteinaceous growth substrates over carbohydrates, which may enable these fungi to also degrade animal biomass. These capabilities for degrading plant and animal substrates suggest a duality in lifestyle that could enable pathogenic species of Onygenales to transfer from soil to animal hosts.

Recurrent Horizontal Transfer of Bacterial Toxin Genes to Eukaryotes

In this work, we report likely recurrent horizontal (lateral) gene transfer events of genes encoding pore-forming toxins of the aerolysin family between species belonging to different kingdoms of life. Clustering based on pairwise similarity and phylogenetic analysis revealed several distinct aerolysin sequence groups, each containing proteins from multiple kingdoms of life. These results strongly support at least six independent transfer events between distantly related phyla in the evolutionary history of one protein family and discount selective retention of ancestral genes as a plausible explanation for this patchy phylogenetic distribution. We discuss the possible roles of these proteins and show evidence for a convergent new function in two extant species. We hypothesize that certain gene families are more likely to be maintained following horizontal gene transfer from commensal or pathogenic organism to its host if they 1) can function alone; and 2) are immediately beneficial for the ecology of the organism, as in the case of pore-forming toxins which can be utilized in multicellular organisms for defense and predation.

Extending the Aerolysin Family: From Bacteria to Vertebrates

A number of bacterial virulence factors have been observed to adopt structures similar to that of aerolysin, the principal toxin of Aeromonas species. However, a comprehensive description of architecture and structure of the aerolysin-like superfamily has not been determined. In this study, we define a more compact aerolysin-like domain – or aerolysin fold – and show that this domain is far more widely spread than anticipated since it can be found throughout kingdoms. The aerolysin-fold could be found in very diverse domain and functional contexts, although a toxic function could often be assigned. Due to this diversity, the borders of the superfamily could not be set on a sequence level. As a border-defining member, we therefore chose pXO2-60 – a protein from the pathogenic pXO2 plasmid of Bacillus anthracis. This fascinating protein, which harbors a unique ubiquitin-like fold domain at the C-terminus of the aerolysin-domain, nicely illustrates the diversity of the superfamily. Its putative role in the virulence of B. anthracis and its three dimensional model are discussed.

LTR retrotransposons in fungi.

Transposable elements with long terminal direct repeats (LTR TEs) are one of the best studied groups of mobile elements. They are ubiquitous elements present in almost all eukaryotic genomes. Their number and state of conservation can be a highlight of genome dynamics. We searched all published fungal genomes for LTR-containing retrotransposons, including both complete, functional elements and remnant copies. We identified a total of over 66,000 elements, all of which belong to the Ty1/Copia or Ty3/Gypsy superfamilies. Most of the detected Gypsy elements represent Chromoviridae, i.e. they carry a chromodomain in the pol ORF. We analyzed our data from a genome-ecology perspective, looking at the abundance of various types of LTR TEs in individual genomes and at the highest-copy element from each genome. The TE content is very variable among the analyzed genomes. Some genomes are very scarce in LTR TEs (<50 elements), others demonstrate huge expansions (>8000 elements). The data shows that transposon expansions in fungi usually involve an increase both in the copy number of individual elements and in the number of element types. The majority of the highest-copy TEs from all genomes are Ty3/Gypsy transposons. Phylogenetic analysis of these elements suggests that TE expansions have appeared independently of each other, in distant genomes and at different taxonomical levels. We also analyzed the evolutionary relationships between protein domains encoded by the transposon pol ORF and we found that the protease is the fastest evolving domain whereas reverse transcriptase and RNase H evolve much slower and in correlation with each other.

Independent Subtilases Expansions in Fungi Associated With Animals

Many socially important fungi encode an elevated number of subtilisin-like serine proteases, which have been shown to be involved in fungal mutualisms with grasses and in parasitism of insects, nematodes, plants, other fungi, and mammalian skin. These proteins have endopeptidase activities and constitute a significant part of fungal secretomes. Here, we use comparative genomics to investigate the relationship between the quality and quantity of serine proteases and the ability of fungi to cause disease in invertebrate and vertebrate animals. Our screen of previously unexamined fungi allowed us to annotate and identify nearly 1000 subtilisin-containing proteins and to describe six new categories of serine proteases. Architectures of predicted proteases reveal novel combinations of subtilisin domains with other, co-occurring domains. Phylogenetic analysis of the most common clade of fungal proteases, proteinase K, showed that gene family size changed independently in fungi, pathogenic to invertebrates (Hypocreales) and vertebrates (Onygenales). Interestingly, simultaneous expansions in the S8 and S53 families of subtilases in a single fungal species are rare. Our analysis finds that closely related systemic human pathogens may not show the same gene family expansions, and that related pathogens and nonpathogens may show the same type of gene family expansion. Therefore, the number of proteases does not appear to relate to pathogenicity. Instead, we hypothesize that the number of fungal serine proteases in a species is related to the use of the animal as a food source, whether it is dead or alive.

Adaptative potential of the Lactococcus lactis IL594 strain encoded in its 7 plasmids.

The extrachromosomal gene pool plays a significant role both in evolution and in the environmental adaptation of bacteria. The L. lactis subsp. lactis IL594 strain contains seven plasmids, named pIL1 to pIL7, and is the parental strain of the plasmid-free L. lactis IL1403, which is one of the best characterized lactococcal strains of LAB. Complete nucleotide sequences of pIL1 (6,382 bp), pIL2 (8,277 bp), pIL3 (19,244 bp), pIL4 (48,979), pIL5 (23,395), pIL6 (28,435 bp) and pIL7 (28,546) were established and deposited in the generally accessible database (GeneBank). Nine highly homologous repB-containing replicons, belonging to the lactococcal theta-type replicons, have been identified on the seven plasmids. Moreover, a putative region involved in conjugative plasmid mobilization was found on four plasmids, through identification of the presence of mob genes and/or oriT sequences. Detailed bioinformatic analysis of the plasmid nucleotide sequences provided new insight into the repertoire of plasmid-encoded functions in L. lactis, and indicated that plasmid genes from IL594 strain can be important for L. lactis adaptation to specific environmental conditions (e.g. genes coding for proteins involved in DNA repair or cold shock response) as well as for technological processes (e.g. genes encoding citrate and lactose utilization, oligopeptide transport, restriction-modification system). Moreover, global gene analysis indicated cooperation between plasmid- and chromosome-encoded metabolic pathways.

A novel protein kinase-like domain in a selenoprotein, widespread in the tree of life.

Selenoproteins serve important functions in many organisms, usually providing essential oxidoreductase enzymatic activity, often for defense against toxic xenobiotic substances. Most eukaryotic genomes possess a small number of these proteins, usually not more than 20. Selenoproteins belong to various structural classes, often related to oxidoreductase function, yet a few of them are completely uncharacterised. Here, the structural and functional prediction for the uncharacterised selenoprotein O (SELO) is presented. Using bioinformatics tools, we predict that SELO protein adopts a three-dimensional fold similar to protein kinases. Furthermore, we argue that despite the lack of conservation of the “classic” catalytic aspartate residue of the archetypical His-Arg-Asp motif, SELO kinases might have retained catalytic phosphotransferase activity, albeit with an atypical active site. Lastly, the role of the selenocysteine residue is considered and the possibility of an oxidoreductase-regulated kinase function for SELO is discussed. The novel kinase prediction is discussed in the context of functional data on SELO orthologues in model organisms, FMP40 a.k.a.YPL222W (yeast), and ydiU (bacteria). Expression data from bacteria and yeast suggest a role in oxidative stress response. Analysis of genomic neighbourhoods of SELO homologues in the three domains of life points toward a role in regulation of ABC transport, in oxidative stress response, or in basic metabolism regulation. Among bacteria possessing SELO homologues, there is a significant over-representation of aquatic organisms, also of aerobic ones. The selenocysteine residue in SELO proteins occurs only in few members of this protein family, including proteins from Metazoa, and few small eukaryotes (Ostreococcus, stramenopiles). It is also demonstrated that enterobacterial mchC proteins involved in maturation of bactericidal antibiotics, microcins, form a distant subfamily of the SELO proteins. The new protein structural domain, with a putative kinase function assigned, expands the known kinome and deserves experimental determination of its biological role within the cell-signaling network.

Sensor Domains Encoded in Bacillus anthracis Virulence Plasmids Prevent Sporulation by Hijacking a Sporulation Sensor Histidine Kinase

Anthrax toxin and capsule, determinants for successful infection by Bacillus anthracis, are encoded on the virulence plasmids pXO1 and pXO2, respectively. Each of these plasmids also encodes proteins that are highly homologous to the signal sensor domain of a chromosomally encoded major sporulation sensor histidine kinase (BA2291) in this organism. B. anthracis Sterne overexpressing the plasmid pXO2-61-encoded signal sensor domain exhibited a significant decrease in sporulation that was suppressed by the deletion of the BA2291 gene. Expression of the sensor domains from the pXO1-118 and pXO2-61 genes in Bacillus subtilis strains carrying the B. anthracis sporulation sensor kinase BA2291 gene resulted in BA2291-dependent inhibition of sporulation. These results indicate that sporulation sensor kinase BA2291 is converted from an activator to an inhibitor of sporulation in its native host by the virulence plasmid-encoded signal sensor domains. We speculate that activation of these signal sensor domains contributes to the initiation of B. anthracis sporulation in the bloodstream of its infected host, a salient characteristic in the virulence of this organism, and provides an additional role for the virulence plasmids in anthrax pathogenesis.