Genome reconstruction and comparative genomics of pneumocystis jirovecii

Pneumocystis jirovecii is a fungus belonging to a basal lineage of the Ascomycotina, the Taphrinomycotina subphylum. It is a parasite specific to humans that dwells primarily in the lung and can cause severe pneumonia in individuals with debilitated immune system. Despite its clinical importance, many aspects of its biology remain poorly understood, at least in part because of the lack of a continuous in vitro cultivation system. The present thesis consists in the genome reconstruction and comparative genomics of P. jirovecii. It is made of three parts: (i) the de novo sequencing of P. jirovecii genome starting from a single broncho- alveolar lavage fluid of a single patient (ii) the de novo sequencing of the genome of the plant pathogen Taphrina deformans, a fungus closely related to P. jirovecii, and (iii) the genome scale comparison of P. jirovecii to other Taphrinomycotina members. Enrichment in P. jirovecii cells by immuno-precipitation, whole DNA random amplification, two complementary high throughput DNA sequencing methods, and in silico sorting and assembly of sequences were used for the de novo reconstruction of P. jirovecii genome from the microbiota of a single clinical specimen. An iterative ad hoc pipeline as well as numerical simulations was used to recover P. jirovecii sequences while purging out contaminants and assembly or amplification chimeras. This strategy produced a 8.1 Mb assembly, which encodes 3,898 genes. Homology searches, mapping on biochemical pathways atlases, and manual validations revealed that this genome lacks (i) most of the enzymes dedicated to the amino acids biosyntheses, and (ii) most virulence factors observed in other fungi, e.g. the glyoxylate shunt pathway and specific peptidases involved in the degradation of the host cell membrane. The same analyses applied to the available genomic sequences from Pneumocystis carinii the species infecting rats and Pneumocystis murina the species infecting mice revealed the same deficiencies. The genome sequencing of T. deformans yielded a 13 Mb assembly, which encodes 5,735 genes. T. deformans possesses enzymes involved plant cell wall degradation, secondary metabolism, the glyoxylate cycle, detoxification, sterol biosynthesis, as well as the biosyntheses of plant hormones such as abscisic acid or indole-3-acetic acid. T. deformans also harbors gene subsets that have counterparts in plant saprophytes or pathogens, which is consistent with its alternate saprophytic and pathogenic lifestyles. Mating genes were also identified. The homothallism of this fungus suggests a mating-type switching mechanism. Comparative analyses indicated that 81% of P. jirovecii genes are shared with eight other Taphrinomycotina members, including T. deformans, P. carinii and P. murina. These genes are mostly involved in housekeeping activities. The genes specific to the Pneumocystis genus represent 8%, and are involved in RNA metabolism and signaling. The signaling is known to be crucial for interaction of Pneumocystis spp with their environment. Eleven percent are unique to P. jirovecii and encode mostly proteins of unknown function. These genes in conjunction with other ones (e.g. the major surface glycoproteins) might govern the interaction of P. jirovecii with its human host cells, and potentially be responsible of the host specificity. P. jirovecii exhibits a reduced genome in size with a low GC content, and most probably scavenges vital compounds such as amino acids and cholesterol from human lungs. Consistently, its genome encodes a large set of transporters (ca. 22% of its genes), which may play a pivotal role in the acquisition of these compounds. All these features are generally observed in obligate parasite of various kingdoms (bacteria, protozoa, fungi). Moreover, epidemiological studies failed to evidence a free-living form of the fungus and Pneumocystis spp were shown to co-evolved with their hosts. Given also the lack of virulence factors, our observations strongly suggest that P. jirovecii is an obligate parasite specialized in the colonization of human lungs, and which causes disease only in individuals with compromised immune system. The same conclusion is most likely true for all other Pneumocystis spp in their respective mammalian host. - Pneumocystis jirovecii est un champignon appartenant à ine branche basale des Ascomycotina, le sous-embranchement des Taphrinomycotina. C'est un parasite spécifique aux humains qui réside principalement dans les poumons, et qui peut causer des pneumonies sévères chez des individus ayant un système immunitaire déficient. En dépit de son importance clinique, de nombreux aspects de sa biologie demeurent,largement méconnus, au moins en partie à cause de l'absence d'un système de culture in vitro continu. Cette thèse traite de la reconstruction du génome et de la génomique comparative de P. jirovecii. Elle comporte trois parties: (i) le séquençage de novo du génome de P. jirovecii à partir d'un lavage broncho-alvéolaire provenant d'un seul patient, (ii) le séquençage de novo du génome d'un champignon pathogène de plante Taphrina deformans qui est phylogénétiquement proche de P. jirovecii, et (iii) la comparaison du génome de P. jirovecii à celui d'autres membres du sous-embranchement des Taphrinomycotina. Un enrichissement en cellules de P. jirovecii par immuno-précipitation, une amplification aléatoire des molécules d'ADN, deux méthodes complémentaires de séquençage à haut débit, un tri in silico et un assemblage des séquences ont été utilisés pour reconstruire de novo le génome de P. jirovecii à partir du microbiote d'un seul échantillon clinique. Un pipeline spécifique ainsi que des simulations numériques ont été utilisés pour récupérer les séquences de P. jirovecii tout en éliminant les séquences contaminants et les chimères d'amplification ou d'assemblage. Cette stratégie a produit un assemblage de 8.1 Mb, qui contient 3898 gènes. Les recherches d'homologies, de cartographie des voies métaboliques et des validations manuelles ont révélé que ce génome est dépourvu (i) de la plupart des enzymes dédiées à la biosynthèse des acides aminés, et (ii) de la plupart des facteurs de virulence observés chez d'autres champignons, par exemple, le cycle du glyoxylate ainsi que des peptidases spécifiques impliquées dans la dégradation de la membrane de la cellule hôte. Les analyses appliquées aux données génomiques disponibles de Pneumocystis carinii, l'espèce infectant les rats, et de Pneumocystis murina, l'espèce infectant les souris, ont révélé les mêmes déficiences. Le séquençage du génome de T. deformans a généré un assemblage de 13.3 Mb qui contient 5735 gènes. T. deformans possède les gènes codant pour les enzymes impliquées dans la dégradation des parois cellulaires des plantes, le métabolisme secondaire, le cycle du glyoxylate, la détoxification, la biosynthèse des stérols ainsi que la biosynthèse d'hormones de plantes telles que l'acide abscissique ou l'acide indole 3-acétique. T. deformans possède également des sous-ensembles de gènes présents exclusivement chez des saprophytes ou des pathogènes de plantes, ce qui est consistent avec son mode de vie alternatif saprophyte et pathogène. Des gènes impliqués dans la conjugaison ont été identifiés. L'homothallisme de ce champignon suggère mécanisme de permutation du type conjuguant. Les analyses comparatives ont démontré que 81% des gènes de P. jirovecii sont présent chez les autres membres du sous-embranchement des Taphrinomycotina. Ces gènes sont essentiellement impliqués dans le métabolisme basai. Les gènes spécifiques au genre Pneumocystis représentent 8%, et sont impliqués dans le métabolisme de l'ARN et la signalisation. La signalisation est connue pour être cruciale pour l'interaction des espèces de Pneumocystis avec leur environnement. Les gènes propres à P. jirovecii représentent 11% et codent en majorité pour des protéines dont la fonction est inconnue. Ces gènes en conjonction avec d'autres (par exemple, les glycoprotéines de surface), pourraient être déterminants dans l'interaction de P. jirovecii avec les cellules de l'hôte humain, et être potentiellement responsable de la spécificité d'hôte. P. jirovecii possède un génome de taille réduite à faible pourcentage en GC et récupère très probablement des composés vitaux comme les acides aminés et le cholestérol à partir des poumons humains. De manière consistante, son génome code pour de nombreux transporteurs (22% de ses gènes), qui pourraient jouer un rôle essentiel dans l'acquisition de ces composés. Ces caractéristiques sont généralement observées chez les parasites obligatoires de plusieurs règnes (bactéries, protozoaires, champignons). De plus, les études épidémiologiques n'ont pas réussi à prouver l'existence d'ime forme vivant librement du champignon. Etant donné également l'absence de facteurs de virulence, nos observations suggèrent que P. jirovecii est un parasite obligatoire spécialisé dans la colonisation des poumons humains, ne causant une maladie que chez des individus ayant un système immunitaire compromis. La même conclusion est très probablement applicable à toutes les autres espèces de Pneumocystis dans leur hôte mammifère respectif.

ExPASy: SIB bioinformatics resource portal.

ExPASy (http://www.expasy.org) has worldwide reputation as one of the main bioinformatics resources for proteomics. It has now evolved, becoming an extensible and integrative portal accessing many scientific resources, databases and software tools in different areas of life sciences. Scientists can henceforth access seamlessly a wide range of resources in many different domains, such as proteomics, genomics, phylogeny/evolution, systems biology, population genetics, transcriptomics, etc. The individual resources (databases, web-based and downloadable software tools) are hosted in a 'decentralized' way by different groups of the SIB Swiss Institute of Bioinformatics and partner institutions. Specifically, a single web portal provides a common entry point to a wide range of resources developed and operated by different SIB groups and external institutions. The portal features a search function across 'selected' resources. Additionally, the availability and usage of resources are monitored. The portal is aimed for both expert users and people who are not familiar with a specific domain in life sciences. The new web interface provides, in particular, visual guidance for newcomers to ExPASy.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

An ontology to clarify homology-related concepts.

Although homology is a fundamental concept in biology and is one of the shared channels of communication universal to all biology, it is difficult to find a consensus definition. Indeed, the interpretations of homology have changed as biology has progressed. New terms, such as paramorphism, have been introduced into the literature with mixed success. In addition, different research fields operate with different definitions of homology, for example the mechanistic usage of evo-devo is not strictly historical and would not be acceptable in cladistics. This makes a global understanding of homology complex, whereas the integration of evolutionary concepts into bioinformatics and genomics is increasingly important. We propose an ontology organizing homology and related concepts and hope this solution will also facilitate the integration and sharing of knowledge among the community.

Constraints on the late thermotectonic evolution of the Western Alps - Evidence for episodic rapid uplift

Forty-two new apatite and zircon fission track ages are presented for samples from the Western Alps in southern Switzerland, northern Italy, and southeastern France. Measured ages plotted against assumed closure temperatures yield cooling patterns for the final cooling, uplift, and exhumation of the Western Alps. Similar fission track zircon ages in the Penninic Gran Paradiso massif, Dent Blanche nappe, Sesia-Lanzo Zone, and Ivrea Zone indicate cooling of all four units to approximately 225-degrees-C by 33 Ma. Differences in apatite ages reveal differential cooling of the four blocks between 33 Ma and the present. In the Sesia-Lanzo Zone, similarity of apatite ages regardless of elevation, together with near-volcanic confined fission track length patterns suggest rapid cooling and uplift at approximately 25 Ma compared with slow cooling of other Western Alps units around 12 Ma. Uplift is thus not continuous but episodic, often over a short time interval beyond the resolution of other methods. Such episodes of uplift, as revealed here in the Sesia-Lanzo Zone, may be the rule rather than the exception.

Duplication and concerted evolution in a master sex determiner under balancing selection.

The transformer (tra) gene is a key regulator in the signalling hierarchy controlling all aspects of somatic sexual differentiation in Drosophila and other insects. Here, we show that six of the seven sequenced ants have two copies of tra. Surprisingly, the two paralogues are always more similar within species than among species. Comparative sequence analyses indicate that this pattern is owing to the ongoing concerted evolution after an ancestral duplication rather than independent duplications in each of the six species. In particular, there was strong support for inter-locus recombination between the paralogues of the ant Atta cephalotes. In the five species where the location of paralogues is known, they are adjacent to each other in four cases and separated by only few genes in the fifth case. Because there have been extensive genomic rearrangements in these lineages, this suggests selection acting to conserve their synteny. In three species, we also find a signature of positive selection in one of the paralogues. In three bee species where information is available, the tra gene is also duplicated, the copies are adjacent and in at least one species there was recombination between paralogues. These results suggest that concerted evolution plays an adaptive role in the evolution of this gene family.

Trophic specialization influences the rate of environmental niche evolution in damselfishes (Pomacentridae).

The rate of environmental niche evolution describes the capability of species to explore the available environmental space and is known to vary among species owing to lineage-specific factors. Trophic specialization is a main force driving species evolution and is responsible for classical examples of adaptive radiations in fishes. We investigate the effect of trophic specialization on the rate of environmental niche evolution in the damselfish, Pomacentridae, which is an important family of tropical reef fishes. First, phylogenetic niche conservatism is not detected in the family using a standard test of phylogenetic signal, and we demonstrate that the environmental niches of damselfishes that differ in trophic specialization are not equivalent while they still overlap at their mean values. Second, we estimate the relative rates of niche evolution on the phylogenetic tree and show the heterogeneity among rates of environmental niche evolution of the three trophic groups. We suggest that behavioural characteristics related to trophic specialization can constrain the evolution of the environmental niche and lead to conserved niches in specialist lineages. Our results show the extent of influence of several traits on the evolution of the environmental niche and shed new light on the evolution of damselfishes, which is a key lineage in current efforts to conserve biodiversity in coral reefs.

Can fisheries-induced evolution shift reference points for fisheries management?

Biological reference points are important tools for fisheries management. Reference points are not static, but may change when a population's environment or the population itself changes. Fisheries-induced evolution is one mechanism that can alter population characteristics, leading to "shifting" reference points by modifying the underlying biological processes or by changing the perception of a fishery system. The former causes changes in "true" reference points, whereas the latter is caused by changes in the yardsticks used to quantify a system's status. Unaccounted shifts of either kind imply that reference points gradually lose their intended meaning. This can lead to increased precaution, which is safe, but potentially costly. Shifts can also occur in more perilous directions, such that actual risks are greater than anticipated. Our qualitative analysis suggests that all commonly used reference points are susceptible to shifting through fisheries-induced evolution, including the limit and "precautionary" reference points for spawning-stock biomass, Blim and Bpa, and the target reference point for fishing mortality, F0.1. Our findings call for increased awareness of fisheries-induced changes and highlight the value of always basing reference points on adequately updated information, to capture all changes in the biological processes that drive fish population dynamics.

Alpine tectonic evolution of a Jurassic subduction-accretionary complex: Deformation, kinematics and 40Ar/39Ar age constraints on the Mesozoic low-grade schists of the Circum-Rhodope Belt in the eastern Rhodope-Thrace region, Bulgaria-Greece

Deformation of the Circum-Rhodope Belt Mesozoic (Middle Triassic to earliest Lower Cretaceous) low-grade schists underneath an arc-related ophiolitic magmatic suite and associated sedimentary successions in the eastern Rhodope-Thrace region occurred as a two-episode tectonic process: (i) Late Jurassic deformation of arc to margin units resulting from the eastern Rhodope-Evros arc-Rhodope terrane continental margin collision and accretion to that margin, and (ii) Middle Eocene deformation related to the Tertiary crustal extension and final collision resulting in the closure of the Vardar ocean south of the Rhodope terrane. The first deformational event D-1 is expressed by Late Jurassic NW-N vergent fold generations and the main and subsidiary planar-linear structures. Although overprinting, these structural elements depict uniform bulk north-directed thrust kinematics and are geometrically compatible with the increments of progressive deformation that develops in same greenschist-facies metamorphic grade. It followed the Early-Middle Jurassic magmatic evolution of the eastern Rhodope-Evros arc established on the upper plate of the southward subducting Maliac-Meliata oceanic lithosphere that established the Vardar Ocean in a supra-subduction back-arc setting. This first event resulted in the thrust-related tectonic emplacement of the Mesozoic schists in a supra-crustal level onto the Rhodope continental margin. This Late Jurassic-Early Cretaceous tectonic event related to N-vergent Balkan orogeny is well-constrained by geochronological data and traced at a regional-scale within distinct units of the Carpatho-Balkan Belt. Following subduction reversal towards the north whereby the Vardar Ocean was subducted beneath the Rhodope margin by latest Cretaceous times, the low-grade schists aquired a new position in the upper plate, and hence, the Mesozoic schists are lacking the Cretaceous S-directed tectono-metamorphic episode whose effects are widespread in the underlying high-grade basement. The subduction of the remnant Vardar Ocean located behind the colliding arc since the middle Cretaceous was responsible for its ultimate closure, Early Tertiary collision with the Pelagonian block and extension in the region caused the extensional collapse related to the second deformational event D-2. This extensional episode was experienced passively by the Mesozoic schists located in the hanging wall of the extensional detachments in Eocene times. It resulted in NE-SW oriented open folds representing corrugation antiforms of the extensional detachment surfaces, brittle faulting and burial history beneath thick Eocene sediments as indicated by 42.1-39.7 Ma Ar-40/Ar-39 mica plateau ages obtained in the study. The results provide structural constraints for the involvement components of Jurassic paleo-subduction zone in a Late Jurassic arc-continental margin collisional history that contributed to accretion-related crustal growth of the Rhodope terrane. (C) 2011 Elsevier Ltd. All rights reserved.

Constraints on palaeodrainage evolution induced by uplift and exhumation on the southern flank of the Zagros-Iranian Plateau

Foreland sedimentary rocks from the northern Fars region of Iran contain a record of deformation associated with the Cenozoic collision between Arabia and Eurasia that resulted in formation of the Zagros orogen. The timing of the deformation associated with this event is poorly known. To address this we conducted a study of Miocene foreland sedimentary rocks (19.7-14.8 Ma) of the Chahar-Makan syncline using clast composition, clay mineralogy and low-temperature fission-track dating. The results showed that most of the sedimentary rocks were sourced from ophiolitic rocks. Detrital apatite fission-track (AFT) age signatures of Miocene sedimentary rocks record exhumation in the hanging wall of the Main Zagros Thrust and confirm that the change from underthrusting of the stretched Arabian margin to widespread crustal thickening and deformation in the Zagros region is no younger than 19.7 Ma. A transition from Late Oligocene to Mesozoic-Eocene AFT detrital age signatures between 19.7-16.6 Ma and 16.6-13.8 Ma is interpreted to reflect a possible rearrangement of palaeodrainage distribution that resulted from folding and expansion-uplift of the Zagros-Iranian Plateau region.

Functional MRI evaluation of liver tumour response after radiofrequency: short- and mid-term evolution of diffusion parameters

Purpose: To evaluate the short- and mid-term evolutions of the apparent diffusion coefficient of lesions treated with RF, in order to determine if the ADC can be used as a marker of tumour response. Methods and Materials: Twenty patients were treated for a liver malignancy with RF and were examined on a 1.5 T/3.0 T machine with T2, gadolinium-enhanced T1 and diffusion sequences: before treatment (< 1 month), just after treatment (< 1 month) and midterm (3-6 months). The ADC was measured in the whole lesion and in the area with the most restricted diffusion (MRDA). The ROI size was also measured on the diffusion map. The Pearson/ANOVA tests were used. Results: All patients were successfully treated with complete disappearance of CE. The lesional size on T2 showed a negative evolution in time (p < 0.002). The ADC in the whole lesion showed a bell-shaped evolution (increasing just after RF, then decreasing, p = 0.02). The ROI size on the diffusion map followed a similar course (p = 0.01). For the MRDA, such evolutions were also found, but they were not significant. There was a negative correlation between CE and the ADC (p < 0.02) and between the lesional size on T2 and ADC (p = 0.03) in the whole lesion. There were also positive correlations between the ROI size and ADC (p = 0.0008) and between CE and the size on T2 (p = 0.0002). The ADC in MRDA showed some non-significant correlations with other variables. Conclusion: The lesions successfully treated with RF have a clear and predictable evolution in terms of T2 size, CE and ADC.

Domain shuffling in a sensor protein contributed to the evolution of insect pathogenicity in plant-beneficial Pseudomonas protegens.

Pseudomonas protegens is a biocontrol rhizobacterium with a plant-beneficial and an insect pathogenic lifestyle, but it is not understood how the organism switches between the two states. Here, we focus on understanding the function and possible evolution of a molecular sensor that enables P. protegens to detect the insect environment and produce a potent insecticidal toxin specifically during insect infection but not on roots. By using quantitative single cell microscopy and mutant analysis, we provide evidence that the sensor histidine kinase FitF is a key regulator of insecticidal toxin production. Our experimental data and bioinformatic analyses indicate that FitF shares a sensing domain with DctB, a histidine kinase regulating carbon uptake in Proteobacteria. This suggested that FitF has acquired its specificity through domain shuffling from a common ancestor. We constructed a chimeric DctB-FitF protein and showed that it is indeed functional in regulating toxin expression in P. protegens. The shuffling event and subsequent adaptive modifications of the recruited sensor domain were critical for the microorganism to express its potent insect toxin in the observed host-specific manner. Inhibition of the FitF sensor during root colonization could explain the mechanism by which P. protegens differentiates between the plant and insect host. Our study establishes FitF of P. protegens as a prime model for molecular evolution of sensor proteins and bacterial pathogenicity.

The influence of learning on evolution: a mathematical framework.

The Baldwin effect can be observed if phenotypic learning influences the evolutionary fitness of individuals, which can in turn accelerate or decelerate evolutionary change. Evidence for both learning-induced acceleration and deceleration can be found in the literature. Although the results for both outcomes were supported by specific mathematical or simulation models, no general predictions have been achieved so far. Here we propose a general framework to predict whether evolution benefits from learning or not. It is formulated in terms of the gain function, which quantifies the proportional change of fitness due to learning depending on the genotype value. With an inductive proof we show that a positive gain-function derivative implies that learning accelerates evolution, and a negative one implies deceleration under the condition that the population is distributed on a monotonic part of the fitness landscape. We show that the gain-function framework explains the results of several specific simulation models. We also use the gain-function framework to shed some light on the results of a recent biological experiment with fruit flies.