Species distributions and interactions in a changing world : from understanding to forecasting

Les écosystèmes fournissent de nombreuses ressources et services écologiques qui sont utiles à la population humaine. La biodiversité est une composante essentielle des écosystèmes et maintient de nombreux services. Afin d'assurer la permanence des services écosystémiques, des mesures doivent être prises pour conserver la biodiversité. Dans ce but, l'acquisition d'informations détaillées sur la distribution de la biodiversité dans l'espace est essentielle. Les modèles de distribution d'espèces (SDMs) sont des modèles empiriques qui mettent en lien des observations de terrain (présences ou absences d'une espèce) avec des descripteurs de l'environnement, selon des courbes de réponses statistiques qui décrive la niche réalisée des espèces. Ces modèles fournissent des projections spatiales indiquant les lieux les plus favorables pour les espèces considérées. Le principal objectif de cette thèse est de fournir des projections plus réalistes de la distribution des espèces et des communautés en montagne pour le climat présent et futur en considérant non-seulement des variables abiotiques mais aussi biotiques. Les régions de montagne et l'écosystème alpin sont très sensibles aux changements globaux et en même temps assurent de nombreux services écosystémiques. Cette thèse est séparée en trois parties : (i) fournir une meilleure compréhension du rôle des interactions biotiques dans la distribution des espèces et l'assemblage des communautés en montagne (ouest des Alpes Suisses), (ii) permettre le développement d'une nouvelle approche pour modéliser la distribution spatiale de la biodiversité, (iii) fournir des projections plus réalistes de la distribution future des espèces ainsi que de la composition des communautés. En me focalisant sur les papillons, bourdons et plantes vasculaires, j'ai détecté des interactions biotiques importantes qui lient les espèces entre elles. J'ai également identifié la signature du filtre de l'environnement sur les communautés en haute altitude confirmant l'utilité des SDMs pour reproduire ce type de processus. A partir de ces études, j'ai contribué à l'amélioration méthodologique des SDMs dans le but de prédire les communautés en incluant les interactions biotiques et également les processus non-déterministes par une approche probabiliste. Cette approche permet de prédire non-seulement la distribution d'espèces individuelles, mais également celle de communautés dans leur entier en empilant les projections (S-SDMs). Finalement, j'ai utilisé cet outil pour prédire la distribution d'espèces et de communautés dans le passé et le futur. En particulier, j'ai modélisé la migration post-glaciaire de Trollius europaeus qui est à l'origine de la structure génétique intra-spécifique chez cette espèce et évalué les risques de perte face au changement climatique. Finalement, j'ai simulé la distribution des communautés de bourdons pour le 21e siècle afin d'évaluer les changements probables dans ce groupe important de pollinisateurs. La diversité fonctionnelle des bourdons va être altérée par la perte d'espèces spécialistes de haute altitude et ceci va influencer la pollinisation des plantes en haute altitude. - Ecosystems provide a multitude of resources and ecological services, which are useful to human. Biodiversity is an essential component of those ecosystems and guarantee many services. To assure the permanence of ecosystem services for future generation, measure should be applied to conserve biodiversity. For this purpose, the acquisition of detailed information on how biodiversity implicated in ecosystem function is distributed in space is essential. Species distribution models (SDMs) are empirical models relating field observations to environmental predictors based on statistically-derived response surfaces that fit the realized niche. These models result in spatial predictions indicating locations of the most suitable environment for the species and may potentially be applied to predict composition of communities and their functional properties. The main objective of this thesis was to provide more accurate projections of species and communities distribution under current and future climate in mountains by considering not solely abiotic but also biotic drivers of species distribution. Mountain areas and alpine ecosystems are considered as particularly sensitive to global changes and are also sources of essential ecosystem services. This thesis had three main goals: (i) a better ecological understanding of biotic interactions and how they shape the distribution of species and communities, (ii) the development of a novel approach to the spatial modeling of biodiversity, that can account for biotic interactions, and (iii) ecologically more realistic projections of future species distributions, of future composition and structure of communities. Focusing on butterfly and bumblebees in interaction with the vegetation, I detected important biotic interactions for species distribution and community composition of both plant and insects along environmental gradients. I identified the signature of environmental filtering processes at high elevation confirming the suitability of SDMs for reproducing patterns of filtering. Using those case-studies, I improved SDMs by incorporating biotic interaction and accounting for non-deterministic processes and uncertainty using a probabilistic based approach. I used improved modeling to forecast the distribution of species through the past and future climate changes. SDMs hindcasting allowed a better understanding of the spatial range dynamic of Trollius europaeus in Europe at the origin of the species intra-specific genetic diversity and identified the risk of loss of this genetic diversity caused by climate change. By simulating the future distribution of all bumblebee species in the western Swiss Alps under nine climate change scenarios for the 21st century, I found that the functional diversity of this pollinator guild will be largely affected by climate change through the loss of high elevation specialists. In turn, this will have important consequences on alpine plant pollination.

Differential control of Notch1 gene transcription by Klf4 and Sp3 transcription factors in normal versus cancer-derived keratinocytes.

In specific cell types like keratinocytes, Notch signaling plays an important pro-differentiation and tumor suppressing function, with down-modulation of the Notch1 gene being associated with cancer development. Besides being controlled by p53, little else is known on regulation of Notch1 gene expression in this context. We report here that transcription of this gene is driven by a TATA-less "sharp peak" promoter and that the minimal functional region of this promoter, which extends from the -342 bp position to the initiation codon, is differentially active in normal versus cancer cells. This GC rich region lacks p53 binding sites, but binds Klf4 and Sp3. This finding is likely to be of biological significance, as Klf4 and, to a lesser extent, Sp3 are up-regulated in a number of cancer cells where Notch1 expression is down-modulated, and Klf4 over-expression in normal cells is sufficient to down-modulate Notch1 gene transcription. The combined knock-down of Klf4 and Sp3 was necessary for the reverse effect of increasing Notch1 transcription, consistent with the two factors exerting an overlapping repressor function through their binding to the Notch1 promoter.

Extracellular protease and phospholipase C are controlled by the global regulatory gene gacA in the biocontrol strain Pseudomonas fluorescens CHA0.

Pseudomonas fluorescens strain CHA0 protects plants from various root diseases. Antibiotic metabolites synthesized by this strain play an important role in disease suppression; their production is mediated by the global activator gene gacA. Here we show by complementation that the gacA gene is also essential for the expression of two extracellular enzymes in P. fluorescens CHA0: phospholipase C and a 47-kDa metalloprotease. In contrast, the production of another exoenzyme, lipase, is not regulated by the gacA gene. Protease, phospholipase and antibiotics of P. fluorescens are all known to be optimally produced at the end of exponential growth; thus, the gacA gene appears to be a general stationary-phase regulator.

In vitro functional correction of Hermansky-Pudlak Syndrome type-1 by lentiviral-mediated gene transfer.

Hermansky-Pudlak syndrome (HPS) is a genetic disorder characterized by oculocutaneous albinism, bleeding tendency and susceptibility to pulmonary fibrosis. No curative therapy is available. Genetic correction directed to the lungs, bone marrow and/or gastro-intestinal tract might provide alternative forms of treatment for the diseases multi-systemic complications. We demonstrate that lentiviral-mediated gene transfer corrects the expression and function of the HPS1 gene in patient dermal melanocytes, which opens the way to development of gene therapy for HPS.

Effects of the social environment on the survival and fungal resistance of ant brood

The phenotype of social animals can be influenced by genetic, maternal and environmental effects, which include social interactions during development. In social insects, the social environment and genetic origin of brood can each influence a whole suite of traits, from individual size to caste differentiation. Here, we investigate to which degree the social environment during development affects the survival and fungal resistance of ant brood of known maternal origin. We manipulated one component of the social environment, the worker/brood ratio, of brood originating from single queens of Formica selysi. We monitored the survival of brood and measured the head size and ability to resist the entomopathogenic fungus Beauveria bassiana of the resulting callow workers. The worker/brood ratio and origin of eggs affected the survival and maturation time of the brood and the size of the resulting callow workers. The survival of the callow workers varied greatly according to their origin, both in controls and when challenged with B. bassiana. However, there was no interaction between the fungal challenge and either the worker/brood ratio or origin of eggs, suggesting that these factors did not affect parasite resistance in the conditions tested. Overall, the social conditions during brood rearing and the origin of eggs had a strong impact on brood traits that are important for fitness. We detected a surprisingly large amount of variation among queens in the survival of their brood reared in standard queenless conditions, which calls for further studies on genetic, maternal and social effects influencing brood development in the social insects.

Suppression of short interfering RNA-mediated gene silencing by the structural proteins of hepatitis C virus.

Viruses have evolved strategies to overcome the antiviral effects of the host at different levels. Besides specific defence mechanisms, the host responds to viral infection via the interferon pathway and also by RNA interference (RNAi). However, several viruses have been identified that suppress RNAi. We addressed the question of whether hepatitis C virus (HCV) suppresses RNAi, using cell lines constitutively expressing green fluorescent protein (GFP) and inducibly expressing HCV proteins. It was found that short interfering RNA-mediated GFP gene silencing was inhibited when the entire HCV polyprotein was expressed. Further studies showed that HCV structural proteins, and in particular envelope protein 2 (E2), were responsible for this inhibition. Co-precipitation assays demonstrated that E2 bound to Argonaute-2 (Ago-2), a member of the RNA-induced silencing complex, RISC. Thus, HCV E2 that interacts with Ago-2 is able to suppress RNAi.

A peptide encoded by the human MAGE3 gene and presented by HLA-B44 induces cytolytic T lymphocytes that recognize tumor cells expressing MAGE3.

The human MAGE3 gene is expressed in a significant proportion of tumors of various histological types, but is silent in normal adult tissues other than testis and placenta. Antigens encoded by MAGE3 may therefore be useful targets for specific antitumor immunization. Two antigenic peptides encoded by the MAGE3 gene have been reported previously. One is presented to cytolytic T lymphocytes (CTL) by HLA-A1, the other by HLA-A2 molecules. Here we show that MAGE3 also codes for a peptide that is presented to CTL by HLA-B44. MAGE3 peptides containing the HLA-B44 peptide binding motif were synthesized. Peptide MEVDPIGHLY, which showed the strongest binding to HLA-B44, was used to stimulate blood T lymphocytes from normal HLA-B44 donors. CTL clones were obtained that recognized not only HLA-B44 cells sensitized with the peptide, but also HLA-B44 tumor cell lines expressing MAGE3. The proportion of metastatic melanomas expressing the MAGE3/HLA-B44 antigen should amount to approximately 17% in the Caucasian population, since 24% of individuals carry the HLA-B44 allele and 76% of these tumors express MAGE3.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Exogenous mRNA delivery and bioavailability in gene transfer mediated by piggyBac transposition.

BACKGROUND: Up to now, the different uptake pathways and the subsequent intracellular trafficking of plasmid DNA have been largely explored. By contrast, the mode of internalization and the intracellular routing of an exogenous mRNA in transfected cells are poorly investigated and remain to be elucidated. The bioavailability of internalized mRNA depends on its intracellular routing and its potential accumulation in dynamic sorting sites for storage: stress granules and processing bodies. This question is of particular significance when a secure transposon-based system able to integrate a therapeutic transgene into the genome is used. Transposon vectors usually require two components: a plasmid DNA, carrying the gene of interest, and a source of transposase allowing the integration of the transgene. The principal drawback is the lasting presence of the transposase, which could remobilize the transgene once it has been inserted. Our study focused on the pharmacokinetics of the transposition process mediated by the piggyBac transposase mRNA transfection. Exogenous mRNA internalization and trafficking were investigated towards a better apprehension and fine control of the piggyBac transposase bioavailability. RESULTS: The mRNA prototype designed in this study provides a very narrow expression window of transposase, which allows high efficiency transposition with no cytotoxicity. Our data reveal that exogenous transposase mRNA enters cells by clathrin and caveolae-mediated endocytosis, before finishing in late endosomes 3 h after transfection. At this point, the mRNA is dissociated from its carrier and localized in stress granules, but not in cytoplasmic processing bodies. Some weaker signals have been observed in stress granules at 18 h and 48 h without causing prolonged production of the transposase. So, we designed an mRNA that is efficiently translated with a peak of transposase production 18 h post-transfection without additional release of the molecule. This confines the integration of the transgene in a very small time window. CONCLUSION: Our results shed light on processes of exogenous mRNA trafficking, which are crucial to estimate the mRNA bioavailability, and increase the biosafety of transgene integration mediated by transposition. This approach provides a new way for limiting the transgene copy in the genome and their remobilization by mRNA engineering and trafficking.

The mif gene is transcriptionally regulated by glucose in insulin-secreting cells.

Macrophage migration inhibitory factor (MIF) is an important regulator of glucose homeostasis. In pancreatic beta-cells, MIF expression is regulated by glucose and its secretion potentiates the glucose-induced insulin secretion. The molecular mechanisms by which glucose mediates its effect on MIF expression are not elucidated. Herein, we report that incubating the differentiated insulin-secreting cell line INS-1 in high glucose concentration increases MIF transcriptional activity as well as the reporter gene activity driven by the -1033 to +63 bp fragment of the MIF promoter. A minimal region located between -187 and -98 bp of this promoter sequence contributes both to basal activity and glucose-responsiveness of the gene. Within this promoter region, two cis-binding sequences were identified by mobility shift assays and footprinting experiments. Both cis-elements interact with nuclear proteins expressed specifically in insulin-secreting cells. In conclusion, we identified a minimal region of the MIF promoter which contributes to the glucose stimulation of the mif gene in insulin-secreting cells.

BAFF production by antigen-presenting cells provides T cell co-stimulation.

The B cell-activating factor from the tumor necrosis factor family (BAFF) is an important regulator of B cell immunity. Recently, we demonstrated that recombinant BAFF also provides a co-stimulatory signal to T cells. Here, we studied expression of BAFF in peripheral blood leukocytes and correlated this expression with BAFF T cell co-stimulatory function. BAFF is produced by antigen-presenting cells (APC). Blood dendritic cells (DC) as well as DC differentiated in vitro from monocytes or CD34+ stem cells express BAFF mRNA. Exposure to bacterial products further up-regulates BAFF production in these cells. A low level of BAFF transcription, up-regulated upon TCR stimulation, was also detected in T cells. Functionally, blockade of endogenous BAFF produced by APC and, to a lesser extent, by T cells inhibits T cell activation. Altogether, this indicates that BAFF may regulate T cell immunity during APC-T cell interactions and as an autocrine factor once T cells have detached from the APC.

Human synthetic lethal inference as potential anti-cancer target gene detection

Background: Two genes are called synthetic lethal (SL) if mutation of either alone is not lethal, but mutation of both leads to death or a significant decrease in organism's fitness. The detection of SL gene pairs constitutes a promising alternative for anti-cancer therapy. As cancer cells exhibit a large number of mutations, the identification of these mutated genes' SL partners may provide specific anti-cancer drug candidates, with minor perturbations to the healthy cells. Since existent SL data is mainly restricted to yeast screenings, the road towards human SL candidates is limited to inference methods. Results: In the present work, we use phylogenetic analysis and database manipulation (BioGRID for interactions, Ensembl and NCBI for homology, Gene Ontology for GO attributes) in order to reconstruct the phylogenetically-inferred SL gene network for human. In addition, available data on cancer mutated genes (COSMIC and Cancer Gene Census databases) as well as on existent approved drugs (DrugBank database) supports our selection of cancer-therapy candidates.Conclusions: Our work provides a complementary alternative to the current methods for drug discovering and gene target identification in anti-cancer research. Novel SL screening analysis and the use of highly curated databases would contribute to improve the results of this methodology.

Characterization and evolution of the novel gene family FAM90A in primates originated by multiple duplication and rearrangement events

Genomic plasticity of human chromosome 8p23.1 region is highly influenced by two groups of complex segmental duplications (SDs), termed REPD and REPP, that mediate different kinds of rearrangements. Part of the difficulty to explain the wide range of phenotypes associated with 8p23.1 rearrangements is that REPP and REPD are not yet well characterized, probably due to their polymorphic status. Here, we describe a novel primate-specific gene family, named FAM90A (family with sequence similarity 90), found within these SDs. According to the current human reference sequence assembly, the FAM90A family includes 24 members along 8p23.1 region plus a single member on chromosome 12p13.31, showing copy number variation (CNV) between individuals. These genes can be classified into subfamilies I and II, which differ in their upstream and 5′-untranslated region sequences, but both share the same open reading frame and are ubiquitously expressed. Sequence analysis and comparative fluorescence in situ hybridization studies showed that FAM90A subfamily II suffered a big expansion in the hominoid lineage, whereas subfamily I members were likely generated sometime around the divergence of orangutan and African great apes by a fusion process. In addition, the analysis of the Ka/Ks ratios provides evidence of functional constraint of some FAM90A genes in all species. The characterization of the FAM90A gene family contributes to a better understanding of the structural polymorphism of the human 8p23.1 region and constitutes a good example of how SDs, CNVs and rearrangements within themselves can promote the formation of new gene sequences with potential functional consequences.

T-Cell Receptors Binding Orientation over Peptide/MHC Class I Is Driven by Long-Range Interactions.

Crystallographic data about T-Cell Receptor - peptide - major histocompatibility complex class I (TCRpMHC) interaction have revealed extremely diverse TCR binding modes triggering antigen recognition. Understanding the molecular basis that governs TCR orientation over pMHC is still a considerable challenge. We present a simplified rigid approach applied on all non-redundant TCRpMHC crystal structures available. The CHARMM force field in combination with the FACTS implicit solvation model is used to study the role of long-distance interactions between the TCR and pMHC. We demonstrate that the sum of the coulomb interactions and the electrostatic solvation energies is sufficient to identify two orientations corresponding to energetic minima at 0° and 180° from the native orientation. Interestingly, these results are shown to be robust upon small structural variations of the TCR such as changes induced by Molecular Dynamics simulations, suggesting that shape complementarity is not required to obtain a reliable signal. Accurate energy minima are also identified by confronting unbound TCR crystal structures to pMHC. Furthermore, we decompose the electrostatic energy into residue contributions to estimate their role in the overall orientation. Results show that most of the driving force leading to the formation of the complex is defined by CDR1,2/MHC interactions. This long-distance contribution appears to be independent from the binding process itself, since it is reliably identified without considering neither short-range energy terms nor CDR induced fit upon binding. Ultimately, we present an attempt to predict the TCR/pMHC binding mode for a TCR structure obtained by homology modeling. The simplicity of the approach and the absence of any fitted parameters make it also easily applicable to other types of macromolecular protein complexes.

Characterization of human gene expression changes after olive oil ingestion: an exploratory approach

Olive oil consumption is protective against risk factors for cardiovascular and cancer diseases. A nutrigenomic approach was performed to assess whether changes in gene expression could occur in human peripheral blood mononuclear cells after oli ve oil ingestion at postprandial state. Six healthy male volunteers ingested, at fasting state, 50 ml of olive oil. Prior to intervention a 1-week washout period with a controlled diet and sunflower oil as the only source of fat was followed. During the 3 days before and on the intervention day, a very low-phenolic compound diet was followed. At baseline (0 h) and at post-ingestion (6 h), total RNA was isolated and gene expression (29,082 genes) was evaluated by microarray. From microarray data, nutrient-gene interactions were observed in genes related to metabolism, cellular processes, cancer, and atherosclerosis (e.g. USP48 by 2.16; OGT by 1.68-fold change) and associated processes such as inflammation (e.g. AKAP13 by 2.30; IL-10 by 1.66-fold change) and DNA damage (e.g. DCLRE1C by 1.47; POLK by 1.44- fold change). When results obtained by microarray were verified by qRT-PCR in nine genes, full concordance was achieved only in the case of up-regulated genes. Changes were observed at a real-life dose of olive oil, as it is daily consumed in some Mediterranean areas. Our results support the hypothesis that postprandial protective changes related to olive oil consumption could be mediated through gene expression changes.