Population structure of three Psammodromus species in the Iberian Peninsula.

The knowledge of a species' population structure is essential for the development of adequate conservation actions as well as for the understanding of its evolution. The population structure is unknown in all species of the Genus Psammodromus, including the Western Sand Racer (Psammodromus occidentalis; a recently described species), the Edward's Sand Racer (P. edwardsianus) and the Spanish Sand Racer (P. hispanicus). In this article, the genetic variability and population structure of Psammodromus edwardsianus, P. hispanicus, and P. occidentalis were studied in the Iberian Peninsula covering their natural geographic distribution. Mitochondrial DNA showed genetically different units in all species with higher genetic variability in their southern populations (latitudinal variation). Genetic differentiation was different among species and contrasted to those of species with similar characteristics. Our results therefore highlight the importance of species-specific studies analysing population structure.

Impact of model complexity on cross-temporal transferability in Maxent species distribution models: An assessment using paleobotanical data

Maximum entropy modeling (Maxent) is a widely used algorithm for predicting species distributions across space and time. Properly assessing the uncertainty in such predictions is non-trivial and requires validation with independent datasets. Notably, model complexity (number of model parameters) remains a major concern in relation to overfitting and, hence, transferability of Maxent models. An emerging approach is to validate the cross-temporal transferability of model predictions using paleoecological data. In this study, we assess the effect of model complexity on the performance of Maxent projections across time using two European plant species (Alnus giutinosa (L.) Gaertn. and Corylus avellana L) with an extensive late Quaternary fossil record in Spain as a study case. We fit 110 models with different levels of complexity under present time and tested model performance using AUC (area under the receiver operating characteristic curve) and AlCc (corrected Akaike Information Criterion) through the standard procedure of randomly partitioning current occurrence data. We then compared these results to an independent validation by projecting the models to mid-Holocene (6000 years before present) climatic conditions in Spain to assess their ability to predict fossil pollen presence-absence and abundance. We find that calibrating Maxent models with default settings result in the generation of overly complex models. While model performance increased with model complexity when predicting current distributions, it was higher with intermediate complexity when predicting mid-Holocene distributions. Hence, models of intermediate complexity resulted in the best trade-off to predict species distributions across time. Reliable temporal model transferability is especially relevant for forecasting species distributions under future climate change. Consequently, species-specific model tuning should be used to find the best modeling settings to control for complexity, notably with paleoecological data to independently validate model projections. For cross-temporal projections of species distributions for which paleoecological data is not available, models of intermediate complexity should be selected.

Life span and reproductive cost explain interspecific variation in the optimal onset of reproduction.

Fitness can be profoundly influenced by the age at first reproduction (AFR), but to date the AFR-fitness relationship only has been investigated intraspecifically. Here, we investigated the relationship between AFR and average lifetime reproductive success (LRS) across 34 bird species. We assessed differences in the deviation of the Optimal AFR (i.e., the species-specific AFR associated with the highest LRS) from the age at sexual maturity, considering potential effects of life history as well as social and ecological factors. Most individuals adopted the species-specific Optimal AFR and both the mean and Optimal AFR of species correlated positively with life span. Interspecific deviations of the Optimal AFR were associated with indices reflecting a change in LRS or survival as a function of AFR: a delayed AFR was beneficial in species where early AFR was associated with a decrease in subsequent survival or reproductive output. Overall, our results suggest that a delayed onset of reproduction beyond maturity is an optimal strategy explained by a long life span and costs of early reproduction. By providing the first empirical confirmations of key predictions of life-history theory across species, this study contributes to a better understanding of life-history evolution.

Predicted decrease in global climate suitability masks regional complexity of invasive fruit fly species response to climate change

Climate change affects the rate of insect invasions as well as the abundance, distribution and impacts of such invasions on a global scale. Among the principal analytical approaches to predicting and understanding future impacts of biological invasions are Species Distribution Models (SDMs), typically in the form of correlative Ecological Niche Models (ENMs). An underlying assumption of ENMs is that species-environment relationships remain preserved during extrapolations in space and time, although this is widely criticised. The semi-mechanistic modelling platform, CLIMEX, employs a top-down approach using species ecophysiological traits and is able to avoid some of the issues of extrapolation, making it highly applicable to investigating biological invasions in the context of climate change. The tephritid fruit flies (Diptera: Tephritidae) comprise some of the most successful invasive species and serious economic pests around the world. Here we project 12 tephritid species CLIMEX models into future climate scenarios to examine overall patterns of climate suitability and forecast potential distributional changes for this group. We further compare the aggregate response of the group against species-specific responses. We then consider additional drivers of biological invasions to examine how invasion potential is influenced by climate, fruit production and trade indices. Considering the group of tephritid species examined here, climate change is predicted to decrease global climate suitability and to shift the cumulative distribution poleward. However, when examining species-level patterns, the predominant directionality of range shifts for 11 of the 12 species is eastward. Most notably, management will need to consider regional changes in fruit fly species invasion potential where high fruit production, trade indices and predicted distributions of these flies overlap.

Coevolution of restriction factors and retroviruses

Les virus exploitent la machinerie cellulaire de l'hôte pour se répliquer. Ils doivent s'adapter pour infecter la cellule hôte de manière optimale tout en échappant à la vigilance du système de défense de l'hôte. Ainsi l'hôte et les virus se livrent à de constantes batailles évolutives. Mon travail de thèse a porté sur l'étude des signatures évolutives de facteurs de l'hôte agissant comme des 'facteurs de restriction' en bloquant la réplication rétrovirale chez les primates. Plus spécifiquement, mon travail a visé à utiliser des données évolutives pour renseigner les analyses fonctionnelles et la biologie. Nous avons étudié le facteur anti-VIH-1 nommé TRIM5a (i) chez les prosimiens pour mieux comprendre son rôle dans le contrôle d'un lentivirus endogène, (ii) dans son activité contre d'autres anciennes infections représentées par des rétrovirus endogènes humains et (iii) en tant que protéine capable de générer des mutants de la capside. Premièrement nous nous sommes intéressés à TRIM5a chez deux espèces de lémuriens dont Microcebus murinus qui porte le lentivirus endogène PSIV dans son génome depuis plusieurs millions d'années,. Nous avons observé que TRIM5a chez M. murinus a un spectre d'activité antivirale réduit à l'opposé de TRIM5a chez le Lemur catta - non porteur du PSIV endogène - qui bloque une large variété de rétrovirus dont le PSIV. De ce fait TRIM5a aurait pu contribuer à protéger certaines espèces de lémuriens vis-à-vis d'anciennes infections par le PSIV. A l'inverse du PSIV, des virus dérivés des rétrovirus endogènes humains HERV-K and HERV-H se sont révélés largement résistants à l'inhibition par TRIM5a. Ces données illustrent une absence de protection par TRIM5a face à d'autres anciennes infections rétrovirales. Puis, pour évaluer l'impact de la protéine TRIM5a humaine sur le VIH-1, nous avons testé l'effet de mutations des résidues sous sélection positive dans la capside du VIH-1 sur l'inhibition par TRIM5a. Nos résultats montrent que TRIM5a ne jouerait pas un rôle significatif dans l'évolution de la capside du VIH-1. Enfin notre travail a porté sur le facteur anti-VIH-1 SAMHD1 récemment découvert, que nous avons séquencé chez 25 espèces de primates. L'analyse évolutive des sites sous sélection positive et des expériences fonctionnelles ont permis d'identifier le domaine de SAMHD1 interagissant avec la protéine lentivirale Vpx. De même que d'autres protéines virales contrecarrent les facteurs de restriction en les menant à la dégradation, nous avons observé que Vpx induit la dégradation de SAMHD1 de manière spécifique à l'espèce. Ces découvertes contribuent à comprendre comment les facteurs de restriction et les virus co-évoluent pour se neutraliser l'un l'autre. - Viruses hijack the host cellular machinery to replicate. They adapt to infect optimally host cells while escaping host defense systems. Viruses and the host coevolve in an evolutionary struggle. My thesis work has been devoted to study the evolutionary signatures of host factors acting as restriction factors that block retroviral replication in primates. Specifically, my work aimed at using evolutionary data to inform functional analyses and biology. We studied the anti-HIV-1 factor TRIM5a (i) in prosimians to better understand its possible role in the control of an endogenous lentivirus, (ii) in its activity against other ancient infections - as represented by HERVs, and (iii) as a protein capable of generating escape mutants in the viral capsid. First, my work focused on two lemur species, one of which was the gray mouse lemur that carries the endogenous lentivirus PSIV integrated in its genome for several million years. TRIM5a from gray mouse lemur exhibited a limited antiviral spectrum as opposed to TRIM5a from ring-tailed lemur - not a host of PSIV - that is able to block diverse retroviruses notably PSIV. These results support the possible contribution of TRIM5a in protecting lemur species from ancient infection by PSIV. In contrast, chimeric viruses derived from two human endogenous retroviruses were broadly resistant to TRIM5a-mediated restriction, suggesting TRIM5a lack of activity against other types of ancient infections. To evaluate the recent impact of human TRIM5a on HIV-1 evolution, we tested whether variants at positively selected sites in the HIV-1 capsid affected the ability of human TRIM5a alleles to restrict HIV-1. Our results indicate that TRIM5a does not play a significant role in the evolution of HIV1 capsid. At last, our work concentrated on the newly discovered anti-HIV-1 restriction factor SAMHD1. We determined its coding sequence in a panel of 25 species of primates. Evolutionary analyses of positively selected sites in SAMHD1 domains and functional assays identified the domain of SAMHD1 interacting with the lentiviral protein Vpx. Similar to other viral countermeasures targeting cellular restriction factors, Vpx was responsible of the degradation of SAMHD1 orthologs in a species-specific manner. These findings contributed to understanding how restriction factors and viruses evolve to counteract each other.

21st century climate change threatens mountain flora unequally across Europe

Continental-scale assessments of 21st century global impacts of climate change on biodiversity have forecasted range contractions for many species. These coarse resolution studies are, however, of limited relevance for projecting risks to biodiversity in mountain systems, where pronounced microclimatic variation could allow species to persist locally, and are ill-suited for assessment of species-specific threat in particular regions. Here, we assess the impacts of climate change on 2632 plant species across all major European mountain ranges, using high-resolution (ca. 100 m) species samples and data expressing four future climate scenarios. Projected habitat loss is greater for species distributed at higher elevations; depending on the climate scenario, we find 36-55% of alpine species, 31-51% of subalpine species and 19-46% of montane species lose more than 80% of their suitable habitat by 2070-2100. While our high-resolution analyses consistently indicate marked levels of threat to cold-adapted mountain florae across Europe, they also reveal unequal distribution of this threat across the various mountain ranges. Impacts on florae from regions projected to undergo increased warming accompanied by decreased precipitation, such as the Pyrenees and the Eastern Austrian Alps, will likely be greater than on florae in regions where the increase in temperature is less pronounced and rainfall increases concomitantly, such as in the Norwegian Scandes and the Scottish Highlands. This suggests that change in precipitation, not only warming, plays an important role in determining the potential impacts of climate change on vegetation.

Threats from climate change to terrestrial vertebrate hotspots in europe.

We identified hotspots of terrestrial vertebrate species diversity in Europe and adjacent islands. Moreover, we assessed the extent to which by the end of the 21(st) century such hotspots will be exposed to average monthly temperature and precipitation patterns which can be regarded as extreme if compared to the climate experienced during 1950-2000. In particular, we considered the entire European sub-continent plus Turkey and a total of 1149 species of terrestrial vertebrates. For each species, we developed species-specific expert-based distribution models (validated against field data) which we used to calculate species richness maps for mammals, breeding birds, amphibians, and reptiles. Considering four global circulation model outputs and three emission scenarios, we generated an index of risk of exposure to extreme climates, and we used a bivariate local Moran's I to identify the areas with a significant association between hotspots of diversity and high risk of exposure to extreme climates. Our results outline that the Mediterranean basin represents both an important hotspot for biodiversity and especially for threatened species for all taxa. In particular, the Iberian and Italian peninsulas host particularly high species richness as measured over all groups, while the eastern Mediterranean basin is particularly rich in amphibians and reptiles; the islands (both Macaronesian and Mediterranean) host the highest richness of threatened species for all taxa occurs. Our results suggest that the main hotspots of biodiversity for terrestrial vertebrates may be extensively influenced by the climate change projected to occur over the coming decades, especially in the Mediterranean bioregion, posing serious concerns for biodiversity conservation.

Modelling the distribution of bats in relation to landscape structure in a temperate mountain environment

1. Landscape modification is often considered the principal cause of population decline in many bat species. Thus, schemes for bat conservation rely heavily on knowledge about species-landscape relationships. So far, however, few studies have quantified the possible influence of landscape structure on large-scale spatial patterns in bat communities. 2. This study presents quantitative models that use landscape structure to predict (i) spatial patterns in overall community composition and (ii) individual species' distributions through canonical correspondence analysis and generalized linear models, respectively. A geographical information system (GIS) was then used to draw up maps of (i) overall community patterns and (ii) distribution of potential species' habitats. These models relied on field data from the Swiss Jura mountains. 3. Fight descriptors of landscape structure accounted for 30% of the variation in bat community composition. For some species, more than 60% of the variance in distribution could be explained by landscape structure. Elevation, forest or woodland cover, lakes and suburbs, were the most frequent predictors. 4. This study shows that community composition in bats is related to landscape structure through species-specific relationships to resources. Due to their nocturnal activities and the difficulties of remote identification, a comprehensive bat census is rarely possible, and we suggest that predictive modelling of the type described here provides an indispensable conservation tool.

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Interactions of tumor necrosis factor (TNF) and TNF receptor family members in the mouse and human.

Ligands of the tumor necrosis factor superfamily (TNFSF) (4-1BBL, APRIL, BAFF, CD27L, CD30L, CD40L, EDA1, EDA2, FasL, GITRL, LIGHT, lymphotoxin alpha, lymphotoxin alphabeta, OX40L, RANKL, TL1A, TNF, TWEAK, and TRAIL) bind members of the TNF receptor superfamily (TNFRSF). A comprehensive survey of ligand-receptor interactions was performed using a flow cytometry-based assay. All ligands engaged between one and five receptors, whereas most receptors only bound one to three ligands. The receptors DR6, RELT, TROY, NGFR, and mouse TNFRH3 did not interact with any of the known TNFSF ligands, suggesting that they either bind other types of ligands, function in a ligand-independent manner, or bind ligands that remain to be identified. The study revealed that ligand-receptor pairs are either cross-reactive between human and mouse (e.g. Tweak/Fn14, RANK/RANKL), strictly species-specific (GITR/GITRL), or partially species-specific (e.g. OX40/OX40L, CD40/CD40L). Interestingly, the receptor binding patterns of lymphotoxin alpha and alphabeta are redundant in the human but not in the mouse system. Ligand oligomerization allowed detection of weak interactions, such as that of human TNF with mouse TNFR2. In addition, mouse APRIL exists as two different splice variants differing by a single amino acid. Although human APRIL does not interact with BAFF-R, the shorter variant of mouse APRIL exhibits weak but detectable binding to mouse BAFF-R.

Role of NOX enzymes in phagocyte response and signaling upon Chlamydia infection

Les membres de l'ordre des Chlamydiales peuvent infecter un choix étendu d'animaux, insectes, et protistes. Comme toutes bactéries intracellulaires obligatoires, les Chlamydiales ont besoin d'une cellule hôte pour se répliquer. Chaque fois qu'une cellule est infectée une lutte commence entre les mécanismes de défense de la cellule et l'arsenal de facteurs de virulence de la bactérie. Dans cette thèse nous nous sommes intéressés à déterminer le rôle de deux mécanismes de l'immunité innée de l'hôte. En premier, nous avons étudié les NADPH oxidases, une source de molécules superoxydantes (MSO). Leur rôle dans la restriction de la réplication de Waddlia chondrophila et Estrella iausannensis a été étudié dans l'organisme modèle Dictyostelium discoideum et les macrophages humains. Différentes protéines Nox étaient nécessaires pour contrôler la réplication de W. chondrophila ou E. Iausannensis. De plus, nous avons déterminé que parmi les Chlamydiales, cinq espèces possédaient une catalase. Cette enzyme peut dégrader l'eau oxygénée, une MSO. L'activité de la catalase a été démontrée in vitro et dans les corps élémentaires. Avant de pouvoir étudier le rôle de NOX2 dans des macrophages infectés avec E. Iausannensis, nous avons dû établir la capacité de la bactérie à se répliquer clans les macrophages avec son trafic intracellulaire. Le deuxième mécanisme d'immunité innée que nous avons étudié est l'autophagie. Dans les cellules infectées l'autophagie permet de digérer les bactéries envahissantes. Deux protéines de la voie autophagique (Atg1 et Atg8) jouent un rôle dans la restriction de la croissance de W. chondrophila dans D. discoideum. D'avantage d'études sur l'immunité innée et les bactéries apparentés aux Chlamydia sont indispensables, car les réponses paraissent être spécifiques pour chaque espèce. - Members of the Chlamydiales order are able to infect a large variety of animals, insects, and protists. These obligate intracellular bacteria require a host cell for replication. Each time a cell is infected a struggle begins between the virulence arsenal of the bacteria and the defense mechanisms activated by the host. Each bacterial species will exhibit a selection of virulence factors that will allow it to overcome the defense of the host in some species, but not others. In this thesis we were interested in dissecting the role of two host innate immunity mechanisms. First we determined the role of NADPH oxidases, a source of reactive oxygen species (ROS), in restricting replication of Waddlia chondrophila and EstreHa lausannensis in the model organism Dictyostelium discoideum and human macrophages. Different Nox proteins were required to restrict growth of W. chondrophila and E. lausannensis. Additionally, we determined that five Chlamydia- related bacterial species encode for catalase, an enzyme that is able to degrade hydrogen peroxide, a ROS. The activity of the catalase was demonstrated in vitro and in elementary bodies. To study the role of NOX2 in macrophages for E. lausannensis we first had to determine the ability of E. lausannensis to grow in macrophages. Besides demonstrating its replication we also determined the intracellular trafficking of E. lausannensis. The second innate immunity mechanism studied was autophagy. Through autophagy bacteria can be targeted to degradation. Atg1 and Atg8, two autophagic proteins appeared restrict W. chondrophila replication in D. discoideum. More studies on innate immunity and Chlamydia-related bacteria are required. It appears that the responses to innate immunity are species specific and it will be difficult to generalize data obtained for W. chondrophila to the Chlamydiales order.

Application of detection probabilities to the design of amphibian monitoring programs in temporary ponds

Failure to detect a species in an area where it is present is a major source of error in biological surveys. We assessed whether it is possible to optimize single-visit biological monitoring surveys of highly dynamic freshwater ecosystems by framing them a priori within a particular period of time. Alternatively, we also searched for the optimal number of visits and when they should be conducted. We developed single-species occupancy models to estimate the monthly probability of detection of pond-breeding amphibians during a four-year monitoring program. Our results revealed that detection probability was species-specific and changed among sampling visits within a breeding season and also among breeding seasons. Thereby, the optimization of biological surveys with minimal survey effort (a single visit) is not feasible as it proves impossible to select a priori an adequate sampling period that remains robust across years. Alternatively, a two-survey combination at the beginning of the sampling season yielded optimal results and constituted an acceptable compromise between sampling efficacy and survey effort. Our study provides evidence of the variability and uncertainty that likely affects the efficacy of monitoring surveys, highlighting the need of repeated sampling in both ecological studies and conservation management.

Regulation of horizontal gene transfer of the catabolic island ICEclc

ICEclc is a mobile genetic element found in two copies on the chromosome of the bacterium Pseudomonas knackmussii B13. ICEclc harbors genes encoding metabolic pathways for the degradation of chlorocatechols (CLC) and 2-aminophenol (2AP). At low frequencies, ICEclc excises from the chromosome, closes into a circular DNA molecule which can transfer to another bacterium via conjugation. Once in the recipient cell, ICEclc can reintegrate into the chromosome by site-specific recombination. This thesis aimed at identifying the regulatory network underlying the decisions for ICEclc horizontal transfer (HGT). The first chapter is an introduction on integrative and conjugative elements (ICEs) more in general, of which ICEclc is one example. In particular I emphasized the current knowledge of regulation and conjugation machineries of the different classes of ICE. In the second chapter, I describe a transcriptional analysis using microarrays and other experiments to understand expression of ICEclc in exponential and stationary phase. By overlaying transcriptomic profiles with Northern hybridizations and RT- PCR data, we established a transcription map for the entire core region of ICEclc, a region assumed to encode the ICE conjugation process. We also demonstrated how transcription of the ICEclc core is maximal in stationary phase, which correlates to expression of reporter genes fused to key ICEclc promoters. In the third chapter, I present a transcriptome analysis of ICEclc in a variety of different host species, in order to explore whether there are species-specific differences. In the fourth chapter, I focus on the role of a curious ICEclc-encoded TetR-type transcriptional repressor. We find that this gene, which we name mfsR, not only controls its own expression but that of a set of genes for a putative multi-drug efflux pump (mfsABC) as well. By using a combination of biochemical and molecular biology techniques, I could show that MfsR specifically binds to operator boxes in two ICEclc promoters (PmfsR and PmfsA), inhibiting the transcription of both the mfsR and mfsABC-orf38184 operons. Although we could not detect a clear phenotype of an mfsABC deletion, we discuss the implications of pump gene reorganizations in ICEclc and close relatives. In the fifth chapter, we find that mfsR not only controls its own expression and that of the mfsABC operon, but is also indirectly controlling ICEclc transfer. Using gene deletions, microarrays, transfer assays and microscopy-based reporter fusions, we demonstrate that mfsR actually controls a small operon of three regulatory genes. The last gene of this mfsR operon, orf17162, encodes a LysR-type activator that when deleted strongly impairs ICEclc transfer. Interestingly, deletion of mfsR leads to transfer competence in almost all cells, thereby overruling the bistability process in the wild-type. In the final sixth chapter, I discuss the relevance of the present thesis and the resulting perspectives for future studies.

Biochemical analysis of female mice urine with reference to endocrine function: a key tool for estrus detection.

Species-specific chemical signals released through urine, sweat, saliva and feces are involved in communication between animals. Urinary biochemical constituents along with pheromones may contribute to variation across reproductive cycles and facilitate to estrus detection. Hence, the present study was designed to analyze such biochemical profiles, such as proteins, carbohydrates, lipids, fatty acids, in response with steroid hormones such as estradiol and progesterone. The experimental groups were normal, prepubertal, ovariectomized, and ovariectomized with estrogentreated female mice. In normal mice, the protein and lipid concentrations in urine were significantly higher in proestrus and estrus phases and the quantity of fatty acids was also comparatively higher in estrus. Furthermore, certain fatty acids, namely tridecanoic, palmitic and oleic acids, were present during proestrus and estrus phases, but were exclusively absent in ovariectomized mice. However, the carbohydrate level was equally maintained throughout the four phases of estrous cycle. For successful communication, higher concentrations of protein and specific fatty acids in estrus are directly involved. The significant increase in estradiol at estrus and progesterone at metestrus seems to be of greater importance in the expression pattern of biochemical constituents and may play a notable role in estrous cycle regulation. Thus, we conclude that the variations observed in the concentration of the biochemical constituents depend on the phase of the reproductive cycle as well as hormonal status of animals. The appearance of protein and specific fatty acids during estrus phase raises the possibility to use these as a urinary indicators for estrus detection.

The molecular signature of the stroma response in prostate cancer-induced osteoblastic bone metastasis highlights expansion of hematopoietic and prostate epithelial stem cell niches.

The reciprocal interaction between cancer cells and the tissue-specific stroma is critical for primary and metastatic tumor growth progression. Prostate cancer cells colonize preferentially bone (osteotropism), where they alter the physiological balance between osteoblast-mediated bone formation and osteoclast-mediated bone resorption, and elicit prevalently an osteoblastic response (osteoinduction). The molecular cues provided by osteoblasts for the survival and growth of bone metastatic prostate cancer cells are largely unknown. We exploited the sufficient divergence between human and mouse RNA sequences together with redefinition of highly species-specific gene arrays by computer-aided and experimental exclusion of cross-hybridizing oligonucleotide probes. This strategy allowed the dissection of the stroma (mouse) from the cancer cell (human) transcriptome in bone metastasis xenograft models of human osteoinductive prostate cancer cells (VCaP and C4-2B). As a result, we generated the osteoblastic bone metastasis-associated stroma transcriptome (OB-BMST). Subtraction of genes shared by inflammation, wound healing and desmoplastic responses, and by the tissue type-independent stroma responses to a variety of non-osteotropic and osteotropic primary cancers generated a curated gene signature ("Core" OB-BMST) putatively representing the bone marrow/bone-specific stroma response to prostate cancer-induced, osteoblastic bone metastasis. The expression pattern of three representative Core OB-BMST genes (PTN, EPHA3 and FSCN1) seems to confirm the bone specificity of this response. A robust induction of genes involved in osteogenesis and angiogenesis dominates both the OB-BMST and Core OB-BMST. This translates in an amplification of hematopoietic and, remarkably, prostate epithelial stem cell niche components that may function as a self-reinforcing bone metastatic niche providing a growth support specific for osteoinductive prostate cancer cells. The induction of this combinatorial stem cell niche is a novel mechanism that may also explain cancer cell osteotropism and local interference with hematopoiesis (myelophthisis). Accordingly, these stem cell niche components may represent innovative therapeutic targets and/or serum biomarkers in osteoblastic bone metastasis.