Comparative genomics of the odorant-binding and chemosensory protein gene families across the arthropoda: origin and evolutionary history of the chemosensory system

Chemoreception is a biological process essential for the survival of animals, as it allows the recognition of important volatile cues for the detection of food, egg-laying substrates, mates or predators, among other purposes. Furthermore, its role in pheromone detection may contribute to evolutionary processes such as reproductive isolation and speciation. This key role in several vital biological processes makes chemoreception a particularly interesting system for studying the role of natural selection in molecular adaptation. Two major gene families are involved in the perireceptor events of the chemosensory system: the odorant-binding protein (OBP) and chemosensory protein (CSP) families. Here, we have conducted an exhaustive comparative genomic analysis of these gene families in twenty Arthropoda species. We show that the evolution of the OBP and CSP gene families is highly dynamic, with a high number of gains and losses of genes, pseudogenes and independent origins of subfamilies. Taken together, our data clearly support the birth-and-death model for the evolution of these gene families with an overall high gene-turnover rate. Moreover, we show that the genome organization of the two families is significantly more clustered than expected by chance and, more important, that this pattern appears to be actively maintained across the Drosophila phylogeny. Finally, we suggest the homologous nature of the OBP and CSP gene families, dating back their MRCA (most recent common ancestor) to 380¿420 Mya, and we propose a scenario for the origin and diversification of these families.

Ecological and historical drivers of diversification in the fly genus Chiastocheta Pokorny.

Coevolution is among the main forces shaping the biodiversity on Earth. In Eurasia, one of the best-known plant-insect interactions showing highly coevolved features involves the fly genus Chiastocheta and its host-plant Trollius. Although this system has been widely studied from an ecological point of view, the phylogenetic relationships and biogeographic history of the flies have remained little investigated. In this integrative study, we aim to test the monophyly of the five Chiastocheta eco-morphological groups, defined by Pellmyr in 1992, by inferring a mitochondrial phylogeny. We further apply a new approach to assess the effect of (i) different molecular substitution rates and (ii) phylogenetic uncertainty on the inference of the spatio-temporal evolution of the group. From a taxonomic point of view, we demonstrate that only two of Pellmyr's groups (rotundiventris and dentifera) are phylogenetically supported, the other species appearing para- or polyphyletic. We also identify the position of C. lophota, which was not included in previous surveys. From a spatio-temporal perspective, we show that the genus arose during the Pliocene in Europe. Our results also indicate that at least four large-scale dispersal events are required to explain the current distribution of Chiastocheta. Moreover, each dispersal to or from Asia is associated with a host-shift and seems to correspond to an increase in speciation rates. Finally, we highlight the correlation between diversification and climatic fluctuations, which indicate that the cycles of global cooling over the last million years had an influence on the radiation of the group.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

Characterising avian haemosporidian communities : ecological factors, parasite persistence and co-infection

Les parasites jouent un rôle clef dans l'évolution des comportements et des traits d'histoire de vie de leurs hôtes. Le parasitisme s'avère parfois dévastateur à l'échelle de population d'hôtes, et peut également altérer certains traits associés à la valeur sélective d'un individu infecté, tels que son succès reproducteur ou encore son taux de mortalité. La coévolution hôte/parasite, qui représente l'une des forces sélectives les plus puissantes dans l'évolution des organismes, peut également conduire les partenaires de l'association parasitaire à s'adapter localement à des environnements hétérogènes. Cette thèse porte sur l'étude de parasites aviaires, du genre Plasmodium, Haemopro- teus et Leucocytozoon (Haemosporidae), naturellement associés à différentes populations de mésanges charbonnières (Parus major) et d'hirondelles des fenêtres (Delichon ur- bicum). Dans un premier temps, nous avons cherché à déterminer comment se distribuent ces parasites au sein de différentes populations hôtes et si ces communautés de parasites sont structurées. Par la suite, la principale question à laquelle nous voulions répondre était de savoir comment ces parasites, et notamment après coexistence de plusieurs lignées génétiques d'Haemosporidae au sein dun même-individu (i.e. co-infection), affectent la physiologie et le succès de reproducteur des hôtes. Nos résultats suggèrent que la distribution des Haemosporidae est principalement gouvernée par la présence d'insectes vecteurs et que la persistance de l'infection chez les hôtes varie en fonction du genre d'Haemosporidae (Chapitre 1-2). Par ailleurs, nous avons trouvé que des lignées de parasite génétiquement distinctes peuvent avoir des effets contrastés sur leurs hôtes. Par exemple, les hôtes exhibent des différences de parasitémie marquées en fonction des lignées de parasites responsable de l'infection. De plus, le succès reproducteur ainsi que la charge parasitaire des mésanges infectées par Plasmodium ou Haemoproteus n'étaient pas affecté par l'infection simultanée avec Leucocytozoon (Chapitre 2-3). Dans le Chapitre 4, j'ai examiné la capacité immunitaire de mésanges charbonnières infectées par des hémosporidies. Les résultats n'ont pas été concluant, et je suggère fortement une réévaluation de ceux-ci dans de futures études. Les mésanges charbonnières ne semblent pas signaler leur statut infectieux par la coloration de leur plumage (Chapitre 5); toutefois, la coloration noire des plumes reflète l'état de stress oxydatif des mésanges, qui dépend lui-même de l'infection parasitaire. La coloration verte pourrait également indiquer la qualité des soins paxentaux délivrés par les mésanges adultes femelles à leurs petits, comme le suggère la corrélation que nous avons observée entre la masse des jeunes d'une nichée et la coloration de leur mère. Les hirondelles capturées en Algérie souffrent plus de l'infection que celles échantillon¬nées en Europe (Chapitre 6). Les similitudes observées entre les communautés de par¬asites affectant les populations européennes et celles des populations nord-africaines suggèrent que la transmission des parasites a lieu lors de la migration vers le sud. A l'instar de nos observations sur les mésanges dans les chapitres 2 et 3, les hirondelles co-infectées ne montrent pas d'altérations de leur condition physique. Cette thèse démontre qu'il existe, au sein des populations de mésanges charbonnières, des interactions antagonistes entre, d'une part, les parasites et leurs hôtes et d'autre part, entre différent parasites. Le résultat de ces interactions antagonistes varie en fonction des espèces et de la zone géographique considérée. Nous avons démontré que les interactions ne suivent pas toujours la théorie, puisque la coevolution qui, en suivant le concept de la virulence, devrait augmenter la charge parasitaire et diminuer la condition physique des hôtes, ne montre pourtant pas d'impact négatif sur les populations de mésanges. Nous pouvons maintenant concentrer nos efforts à la caractérisation des interactions antagonistes. De plus, grâce aux avancées des méthodes moléculaires, nous pouvons suivre et étudier en détails comment ces interactions se manifestent et quels sont leurs effets sur la condition physique des hôtes. - Parasites are key in shaping various behavioural and life-history traits of their hosts. The influence of parasitism on host populations varies from slight to devastating and might influence such parameters as mortality rates or reproductive success. Host-parasite coevolution is one of the most powerful selective forces in evolution and can lead to local adaptation of parasites and hosts in spatially structured environments. In this thesis, I studied haemosporidian parasites in different populations of great tits (Parus major) and house martins (Delichon urbicum). Firstly, I wanted to determine how parasites are distributed and if parasite communities are structured. The main question I wanted to address hereafter was how parasites, and specifically infection with multiple genera of parasites (i.e. co-infection) influenced host physiology and reproductive success. I found that parasite distribution is environmentally driven and could therefore be closely linked to vector prevalence; and that the stability of parasite infection over time is genus-dependent (Chapter 1 - 2). I further found that different haemosporidian lineages might interact differently with their hosts as parasitaemia was strongly lineage-specific and that the presence of Leucocytozoon parasites showed no correlation to Plasmodium or Haemoproteus parasitaemia, nor to great tit reproductive success (Chapter 2-3). In Chapter 4 I examined immune capacity of haemosporidian-infected great tits. The results proved inconclusive, and I strongly suggest re-evaluation hereof in future work. Great tits do not appear to signal parasite infection through plumage colouration (Chapter 5); however, infection did have a link to oxidative stress resistance which is strongly signalled through the black breast stripe, with darker males being more resistant and darker females less resistant. Females might incur different costs associated with darker stripes. This would allow reversal of signaling function. Green colouration could also serve as a cue for female provisioning quality as indicated by the strong correlation between colouration and chick body mass. Breeding house martins caught in Algeria suffer greater haemosporidian infection than European populations (Chapter 6). Similar parasite communities in European and North-African populations suggest transmission of parasites may occur during southward migration. Similarly to what was observed in great tits in Chapter 2 and 3, no relationship was found between parasite co-infection and Swiss house martin body condition. This thesis demonstrates that host-parasite and inter-parasite antagonistic interac¬tions exist in great tit populations. How these interactions play out is species dependent and varies geographically. I have demonstrated that interactions do not always follow the theory, as co-infection - which under the concept of virulence should increase parasitaemia and decrease body condition - showed no negative impact on great tit populations. We can now concentrate our efforts on characterising these antagonistic interactions, and with the advance in molecular methods, track and investigate how these interactions play out and what the effect on host fitness is.

The ecological significance of arbuscular mycorrhizal fungal effects on clonal reproduction in plants

The population ecology of clonal plants depends on the number and distribution of ramets formed during growth. Variation in clonal reproduction has previously been explained by variation in effects of abiotic resource heterogeneity and by plant genotypic variation. Different co-occurring species of the mutualistic arbuscular mycorrhizal fungi (AMF) have been shown to differentially alter growth traits of Prunella vulgaris which we hypothesize would lead to changes in clonal reproduction. Two experiments were carried out to test whether different co-occurring mycorrhizal fungi significantly influence clonal reproduction of P. vulgaris whether this effect also occurs when P. vulgaris is growing in an artificial plant community and how the effects compare with plant genotype effects on clonal growth of P. vulgaris. In the first experiment the number of ramets of P. vulgaris growing in a plant community of simulated calcareous grassland was significantly affected by inoculation with different mycorrhizal fungi. The number of ramets produced by P. vulgaris differed by a factor of up to 1.8 with different mycorrhizal fungi. The fungal effects on the number of new ramets were independent of their effects on the biomass of P. vulgaris. In a second experiment 17 different genotypes of P. vulgaris were inoculated with different mycorrhizal fungi. There were significant main effects of genotypes and mycorrhizal fungi on clonal reproduction of P. vulgaris. The effect of different mycorrhizal fungi contributed more than the effect of plant genotype to variation in size and ramet production. However mean stolon length and spacer length which determine the spatial arrangement of ramets were only significantly affected by plant genotype. There were no mycorrhizal fungal X plant genotype interactions on clonal growth of P. vulgaris indicating that there is no obvious evidence that selection pressures would favor further coevolution between P. vulgaris and mycorrhizal fungal species. In natural communities plants can be colonized by several different AMF at the same time. The effect of the mixed AMF treatment on the growth and clonal reproduction of P. vulgaris could not be predicted from the responses of the plants to the single AMF To what extent however the patterns of colonization by different AMF differ among plants in a natural community is unknown. Since the effects of AMF on growth and clonal reproduction occur on a population of P. vulgaris in a microcosm plant community and because the effects are also as great as those caused by plant genotypic variation we conclude that the effects are strong enough to potentially affect population size and variation of clonal plants in communities.

Comparative Genomics Suggests that the Fungal Pathogen Pneumocystis Is an Obligate Parasite Scavenging Amino Acids from Its Host's Lungs.

Pneumocystis jirovecii is a fungus causing severe pneumonia in immuno-compromised patients. Progress in understanding its pathogenicity and epidemiology has been hampered by the lack of a long-term in vitro culture method. Obligate parasitism of this pathogen has been suggested on the basis of various features but remains controversial. We analysed the 7.0 Mb draft genome sequence of the closely related species Pneumocystis carinii infecting rats, which is a well established experimental model of the disease. We predicted 8'085 (redundant) peptides and 14.9% of them were mapped onto the KEGG biochemical pathways. The proteome of the closely related yeast Schizosaccharomyces pombe was used as a control for the annotation procedure (4'974 genes, 14.1% mapped). About two thirds of the mapped peptides of each organism (65.7% and 73.2%, respectively) corresponded to crucial enzymes for the basal metabolism and standard cellular processes. However, the proportion of P. carinii genes relative to those of S. pombe was significantly smaller for the "amino acid metabolism" category of pathways than for all other categories taken together (40 versus 114 against 278 versus 427, P<0.002). Importantly, we identified in P. carinii only 2 enzymes specifically dedicated to the synthesis of the 20 standard amino acids. By contrast all the 54 enzymes dedicated to this synthesis reported in the KEGG atlas for S. pombe were detected upon reannotation of S. pombe proteome (2 versus 54 against 278 versus 427, P<0.0001). This finding strongly suggests that species of the genus Pneumocystis are scavenging amino acids from their host's lung environment. Consequently, they would have no form able to live independently from another organism, and these parasites would be obligate in addition to being opportunistic. These findings have implications for the management of patients susceptible to P. jirovecii infection given that the only source of infection would be other humans.

Ecological-economic optimization of biodiversity conservation under climate change

Substantial investment in climate change research has led to dire predictions of the impacts and risks to biodiversity. The Intergovernmental Panel on Climate Change fourth assessment report(1) cites 28,586 studies demonstrating significant biological changes in terrestrial systems(2). Already high extinction rates, driven primarily by habitat loss, are predicted to increase under climate change(3-6). Yet there is little specific advice or precedent in the literature to guide climate adaptation investment for conserving biodiversity within realistic economic constraints(7). Here we present a systematic ecological and economic analysis of a climate adaptation problem in one of the world's most species-rich and threatened ecosystems: the South African fynbos. We discover a counterintuitive optimal investment strategy that switches twice between options as the available adaptation budget increases. We demonstrate that optimal investment is nonlinearly dependent on available resources, making the choice of how much to invest as important as determining where to invest and what actions to take. Our study emphasizes the importance of a sound analytical framework for prioritizing adaptation investments(4). Integrating ecological predictions in an economic decision framework will help support complex choices between adaptation options under severe uncertainty. Our prioritization method can be applied at any scale to minimize species loss and to evaluate the robustness of decisions to uncertainty about key assumptions.

Case studies and mathematical models of ecological speciation. 1. Cichlids in a crater lake.

A recent study of a pair of sympatric species of cichlids in Lake Apoyo in Nicaragua is viewed as providing probably one of the most convincing examples of sympatric speciation to date. Here, we describe and study a stochastic, individual-based, explicit genetic model tailored for this cichlid system. Our results show that relatively rapid (<20,000 generations) colonization of a new ecological niche and (sympatric or parapatric) speciation via local adaptation and divergence in habitat and mating preferences are theoretically plausible if: (i) the number of loci underlying the traits controlling local adaptation, and habitat and mating preferences is small; (ii) the strength of selection for local adaptation is intermediate; (iii) the carrying capacity of the population is intermediate; and (iv) the effects of the loci influencing nonrandom mating are strong. We discuss patterns and timescales of ecological speciation identified by our model, and we highlight important parameters and features that need to be studied empirically to provide information that can be used to improve the biological realism and power of mathematical models of ecological speciation.

Behavioural types and ecological effects in a natural population of the cooperative cichlid Neolamprologus pulcher

The ecological relevance of behavioural syndromes is little studied in cooperative breeding systems where it is assumed that the behavioural type might influence individual decisions on helping and dispersal (e.g. shy, nonaggressive and nonexplorative individuals remain philopatric and helpful, whereas bold, aggressive, explorative individuals compete for vacancies outside their group and disperse). We measured the behavioural type of 19 subordinates in the cooperatively breeding cichlid fish Neolamprologus pulcher in their natural environment by quantifying six behavioural traits up to four times ('trials') in three different contexts, by presenting them with a conspecific intruder, a predator or nothing inside a tube. We found only moderate within-context repeatability (intraclass correlation coefficients) of the focal individual's behaviour, except for attacking either the conspecific or the predator inside the tube. The focal individual's attack rate of the tube was also positively affected by its group size. Averaging traits per context removed the between-trial variation, and consequently the across-context repeatability was very high for all six traits, except for territory maintenance. Trait values depended significantly on the context, except for territory defence. Consequently, individuals could be classified into different behavioural types based on their reaction towards the tube, but surprisingly, and opposite to laboratory studies in this species, ranging propensity and territory maintenance were not included in this behavioural syndrome. We suggest that more studies are needed to compare standardized focal personality tests (e.g. exploration propensity) with actual behaviour observed in nature (e.g. ranging and dispersal).