Spatial distribution of economic activities: an empirical approach usingself-organizing maps

The aim of this paper is to analyse the colocation patterns of industries and firms. We study the spatial distribution of firms from different industries at a microgeographic level and from this identify the main reasons for this locational behaviour. The empirical application uses data from Mercantile Registers of Spanish firms (manufacturers and services). Inter-sectorial linkages are shown using self-organizing maps. Key words: clusters, microgeographic data, self-organizing maps, firm location JEL classification: R10, R12, R34

Statistical tests for scaling in the inter-event times of earthquakes in California

We explore in depth the validity of a recently proposed scaling law for earthquake inter-event time distributions in the case of the Southern California, using the waveform cross-correlation catalog of Shearer et al. Two statistical tests are used: on the one hand, the standard two-sample Kolmogorov-Smirnov test is in agreement with the scaling of the distributions. On the other hand, the one-sample Kolmogorov-Smirnov statistic complemented with Monte Carlo simulation of the inter-event times, as done by Clauset et al., supports the validity of the gamma distribution as a simple model of the scaling function appearing on the scaling law, for rescaled inter-event times above 0.01, except for the largest data set (magnitude greater than 2). A discussion of these results is provided.

Comprehensive analysis of RET common and rare variant patientsin a series of Spanish Hirschsprung patients confirms a synergistic effect of both kinds of events

Background. RET is the major gene associated to Hirschsprung disease (HSCR) with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In the present study, we have performed a comprehensive study of our HSCR series evaluating the involvement of both RET rare variants (RVs) and common variants (CVs) in the context of the disease. Methods. RET mutational screening was performed by dHPLC and direct sequencing for the identification of RVs. In addition Taqman technology was applied for the genotyping of 3 RET CVs previously associated to HSCR, including a variant lying in an enhancer domain within RET intron 1 (rs2435357). Statistical analyses were performed using the SPSS v.17.0 to analyze the distribution of the variants. Results. Our results confirm the strongest association to HSCR for the "enhancer" variant, and demonstrate a significantly higher impact of it in male versus female patients. Integration of the RET RVs and CVs analysis showed that in 91.66% of cases with both kinds of mutational events, the enhancer allele is in trans with the allele bearing the RET RV. Conclusions. A gender effect exists on both the transmission and distribution of rare coding and common HSCR causing mutations. In addition, these RET CVs and RVs seem to act in a synergistic way leading to HSCR phenotype.

Evidence for classification of c.1852_1853AA>GC in MLH1 as a neutral variant for Lynch syndrome

Background: Lynch syndrome (LS) is an autosomal dominant inherited cancer syndrome characterized by early onset cancers of the colorectum, endometrium and other tumours. A significant proportion of DNA variants in LS patients are unclassified. Reports on the pathogenicity of the c.1852_1853AA>GC (p.Lys618Ala) variant of the MLH1 gene are conflicting. In this study, we provide new evidence indicating that this variant has no significant implications for LS.Methods: The following approach was used to assess the clinical significance of the p.Lys618Ala variant: frequency in a control population, case-control comparison, co-occurrence of the p.Lys618Ala variant with a pathogenic mutation, co-segregation with the disease and microsatellite instability in tumours from carriers of the variant. We genotyped p.Lys618Ala in 1034 individuals (373 sporadic colorectal cancer [CRC] patients, 250 index subjects from families suspected of having LS [revised Bethesda guidelines] and 411 controls). Three well-characterized LS families that fulfilled the Amsterdam II Criteria and consisted of members with the p.Lys618Ala variant were included to assess co-occurrence and co-segregation. A subset of colorectal tumour DNA samples from 17 patients carrying the p.Lys618Ala variant was screened for microsatellite instability using five mononucleotide markers.Results: Twenty-seven individuals were heterozygous for the p.Lys618Ala variant; nine had sporadic CRC (2.41%), seven were suspected of having hereditary CRC (2.8%) and 11 were controls (2.68%). There were no significant associations in the case-control and case-case studies. The p.Lys618Ala variant was co-existent with pathogenic mutations in two unrelated LS families. In one family, the allele distribution of the pathogenic and unclassified variant was in trans, in the other family the pathogenic variant was detected in the MSH6 gene and only the deleterious variant co-segregated with the disease in both families. Only two positive cases of microsatellite instability (2/17, 11.8%) were detected in tumours from p.Lys618Ala carriers, indicating that this variant does not play a role in functional inactivation of MLH1 in CRC patients.Conclusions: The p.Lys618Ala variant should be considered a neutral variant for LS. These findings have implications for the clinical management of CRC probands and their relatives.

Analysis of the multi-copy gene family FAM90A as a copy number variant in different ethnic backgrounds

Copy number variants contribute extensively to inter-individual genomic differences, but little is known about their inter-population variability and diversity. In a previous study (Bosch et al., 2007; 16:2572-2582), we reported that the primate-specific gene family FAM90A, which accounts for as many as 25 members in the human reference assembly, has expanded the number of FAM90A clusters across the hominoid lineage. Here we examined the copy number variability of FAM90A genes in 260 HapMap samples of European, African, and Asian ancestry, and showed significant inter-population differences (p<0.0001). Based on the recent study of Stranger et al. (2007; 315:848-853), we also explored the correlation between copy number variability and expression levels of the FAM90A gene family. Despite the high genomic variability, we found a low correlation between FAM90A copy number and expression levels, which could be due to the action of independent trans-acting factors. Our results show that FAM90A is highly variable in copy number between individuals and between populations. However, this variability has little impact on gene expression levels, thus highlighting the importance of genomic variability for genes located in regions containing segmental duplications.

Bilinear models for spatio-temporal point distribution analysis: application to extrapolation of left ventricular, biventricular and whole heart cardiac dynamics

In this work we describe the usage of bilinear statistical models as a means of factoring the shape variability into two components attributed to inter-subject variation and to the intrinsic dynamics of the human heart. We show that it is feasible to reconstruct the shape of the heart at discrete points in the cardiac cycle. Provided we are given a small number of shape instances representing the same heart atdifferent points in the same cycle, we can use the bilinearmodel to establish this. Using a temporal and a spatial alignment step in the preprocessing of the shapes, around half of the reconstruction errors were on the order of the axial image resolution of 2 mm, and over 90% was within 3.5 mm. From this, weconclude that the dynamics were indeed separated from theinter-subject variability in our dataset.

Topographical body fat distribution links to amino acid and lipid metabolism in healthy non-obese women.

Visceral adiposity is increasingly recognized as a key condition for the development of obesity related disorders, with the ratio between visceral adipose tissue (VAT) and subcutaneous adipose tissue (SAT) reported as the best correlate of cardiometabolic risk. In this study, using a cohort of 40 obese females (age: 25-45 y, BMI: 28-40 kg/m(2)) under healthy clinical conditions and monitored over a 2 weeks period we examined the relationships between different body composition parameters, estimates of visceral adiposity and blood/urine metabolic profiles. Metabonomics and lipidomics analysis of blood plasma and urine were employed in combination with in vivo quantitation of body composition and abdominal fat distribution using iDXA and computerized tomography. Of the various visceral fat estimates, VAT/SAT and VAT/total abdominal fat ratios exhibited significant associations with regio-specific body lean and fat composition. The integration of these visceral fat estimates with metabolic profiles of blood and urine described a distinct amino acid, diacyl and ether phospholipid phenotype in women with higher visceral fat. Metabolites important in predicting visceral fat adiposity as assessed by Random forest analysis highlighted 7 most robust markers, including tyrosine, glutamine, PC-O 44∶6, PC-O 44∶4, PC-O 42∶4, PC-O 40∶4, and PC-O 40∶3 lipid species. Unexpectedly, the visceral fat associated inflammatory profiles were shown to be highly influenced by inter-days and between-subject variations. Nevertheless, the visceral fat associated amino acid and lipid signature is proposed to be further validated for future patient stratification and cardiometabolic health diagnostics.

Distribution of neuronal and metabolic markers in the human auditory cortex

SUMMARY The human auditory cortex, located on the supratemporal plane of the temporal lobe, is divided in a primary auditory area and several non-primary areas surrounding it. These different areas show anatomical and functional differences. Many studies have focussed on auditory areas in non-human primates, using investigation techniques such as electrophysiological recordings, tracing of neural connections, or immunohistochemical and histochemical staining. Some of these studies have suggested parallel and hierarchical organization of the cortical auditory areas as well as subcortical auditory relays. In humans, only few studies have investigated these regions immunohistochemically, but activation and lesion studies speak in favour of parallel and hierarchical organization, very similar to that of non-human primates. Calcium-binding proteins and metabolic markers were used to investigate possible correlates of hierarchical and parallel organization in man. Calcium-binding proteins, parvalbumin, calretinin and calbindin, modulate the concentration of intracellular free calcium ions and were found in distinct subpopulations of GABAergic neurons in non-human primates species. In our study, their distribution showed several differences between auditory areas: the primary auditory area was darkly stained for both parvalbumin and calbindin, and their expression rapidly decreased while moving away from the primary area. This staining pattern suggests a hierarchical organization of the areas, in which the more darkly stained areas could correspond to an earlier integration level and the areas showing light staining may correspond to higher level integration areas. Parallel organization of primary and non-primary auditory areas was suggested by the complementarity, within a given area, between parvalbumin and calbindin expression across layers. To investigate the possible differences in the energetic metabolism of the cortical auditory areas, several metabolic markers were used: cytochrome oxidase and LDH1 were used as oxidative metabolism markers and LDH5 was used as glycolytic metabolism marker. The results obtained show a difference in the expression of enzymes involved in oxidative metabolism between areas. In the primary auditory area the oxidative metabolism markers were maximally expressed in layer IV. In contrast, higher order areas showed maximal staining in supragranular layers. The expression of LDH5 varied in patches, but did not differ between the different hierarchical auditory areas. The distribution of the two LDH enzymes isoforms also provides information about cellular aspects of metabolic organization, since neurons expressed the LDH1 isoform whereas astrocytes express primarily LDH5, but some astrocytes also contained the LDH1 isoform. This cellular distribution pattern supports the hypothesis of the existence of an astrocyte-neuron lactate shuttle, previously suggested in rodent studies, and in particular of lactate transfer from astrocytes, which produce lactate from the glucose obtained from the circulation, to neurons that use lactate as energy substrate. In conclusion, the hypothesis of parallel and hierarchical organization of the auditory areas can be supported by CaBPs, cytochrome oxidase and LDH1 distribution. Moreover, the two LDHs cellular distribution pattern support the hypothesis of an astrocyte-neuron lactate shuttle in human cortex.

The HSR on chromosome 1 of the house mouse, Mus domesticus: distribution and frequency in Switzerland.

A total of 357 house mice (Mus domesticus) from 83 localities uniformly distributed throughout Switzerland were screened for the presence of a homogenously staining region (HSR) on chromosome 1. Altogether 47 mice from 11 localities were HSR/+ or HSR/HSR. One sample of 11 individuals all had an HSR/HSR karyotype. Almost all mice with the variant were collected from the Rhone valley (HSR frequency: 61%) and Val Bregaglia (HSR frequency: 81%). For samples from most of the area of Switzerland, the HSR was absent. There was no strong association between the geographic distribution of the HSR and the areas of occurrence of metacentrics. However, at Chiggiogna the HSR was found on Rb (1.3). Possible explanations for the HSR polymorphism are discussed.

Time-lapse and probabilistic inversion strategies for plane-wave electromagnetic data.

L'utilisation efficace des systèmes géothermaux, la séquestration du CO2 pour limiter le changement climatique et la prévention de l'intrusion d'eau salée dans les aquifères costaux ne sont que quelques exemples qui démontrent notre besoin en technologies nouvelles pour suivre l'évolution des processus souterrains à partir de la surface. Un défi majeur est d'assurer la caractérisation et l'optimisation des performances de ces technologies à différentes échelles spatiales et temporelles. Les méthodes électromagnétiques (EM) d'ondes planes sont sensibles à la conductivité électrique du sous-sol et, par conséquent, à la conductivité électrique des fluides saturant la roche, à la présence de fractures connectées, à la température et aux matériaux géologiques. Ces méthodes sont régies par des équations valides sur de larges gammes de fréquences, permettant détudier de manières analogues des processus allant de quelques mètres sous la surface jusqu'à plusieurs kilomètres de profondeur. Néanmoins, ces méthodes sont soumises à une perte de résolution avec la profondeur à cause des propriétés diffusives du champ électromagnétique. Pour cette raison, l'estimation des modèles du sous-sol par ces méthodes doit prendre en compte des informations a priori afin de contraindre les modèles autant que possible et de permettre la quantification des incertitudes de ces modèles de façon appropriée. Dans la présente thèse, je développe des approches permettant la caractérisation statique et dynamique du sous-sol à l'aide d'ondes EM planes. Dans une première partie, je présente une approche déterministe permettant de réaliser des inversions répétées dans le temps (time-lapse) de données d'ondes EM planes en deux dimensions. Cette stratégie est basée sur l'incorporation dans l'algorithme d'informations a priori en fonction des changements du modèle de conductivité électrique attendus. Ceci est réalisé en intégrant une régularisation stochastique et des contraintes flexibles par rapport à la gamme des changements attendus en utilisant les multiplicateurs de Lagrange. J'utilise des normes différentes de la norme l2 pour contraindre la structure du modèle et obtenir des transitions abruptes entre les régions du model qui subissent des changements dans le temps et celles qui n'en subissent pas. Aussi, j'incorpore une stratégie afin d'éliminer les erreurs systématiques de données time-lapse. Ce travail a mis en évidence l'amélioration de la caractérisation des changements temporels par rapport aux approches classiques qui réalisent des inversions indépendantes à chaque pas de temps et comparent les modèles. Dans la seconde partie de cette thèse, j'adopte un formalisme bayésien et je teste la possibilité de quantifier les incertitudes sur les paramètres du modèle dans l'inversion d'ondes EM planes. Pour ce faire, je présente une stratégie d'inversion probabiliste basée sur des pixels à deux dimensions pour des inversions de données d'ondes EM planes et de tomographies de résistivité électrique (ERT) séparées et jointes. Je compare les incertitudes des paramètres du modèle en considérant différents types d'information a priori sur la structure du modèle et différentes fonctions de vraisemblance pour décrire les erreurs sur les données. Les résultats indiquent que la régularisation du modèle est nécessaire lorsqu'on a à faire à un large nombre de paramètres car cela permet d'accélérer la convergence des chaînes et d'obtenir des modèles plus réalistes. Cependent, ces contraintes mènent à des incertitudes d'estimations plus faibles, ce qui implique des distributions a posteriori qui ne contiennent pas le vrai modèledans les régions ou` la méthode présente une sensibilité limitée. Cette situation peut être améliorée en combinant des méthodes d'ondes EM planes avec d'autres méthodes complémentaires telles que l'ERT. De plus, je montre que le poids de régularisation des paramètres et l'écart-type des erreurs sur les données peuvent être retrouvés par une inversion probabiliste. Finalement, j'évalue la possibilité de caractériser une distribution tridimensionnelle d'un panache de traceur salin injecté dans le sous-sol en réalisant une inversion probabiliste time-lapse tridimensionnelle d'ondes EM planes. Etant donné que les inversions probabilistes sont très coûteuses en temps de calcul lorsque l'espace des paramètres présente une grande dimension, je propose une stratégie de réduction du modèle ou` les coefficients de décomposition des moments de Legendre du panache de traceur injecté ainsi que sa position sont estimés. Pour ce faire, un modèle de résistivité de base est nécessaire. Il peut être obtenu avant l'expérience time-lapse. Un test synthétique montre que la méthodologie marche bien quand le modèle de résistivité de base est caractérisé correctement. Cette méthodologie est aussi appliquée à un test de trac¸age par injection d'une solution saline et d'acides réalisé dans un système géothermal en Australie, puis comparée à une inversion time-lapse tridimensionnelle réalisée selon une approche déterministe. L'inversion probabiliste permet de mieux contraindre le panache du traceur salin gr^ace à la grande quantité d'informations a priori incluse dans l'algorithme. Néanmoins, les changements de conductivités nécessaires pour expliquer les changements observés dans les données sont plus grands que ce qu'expliquent notre connaissance actuelle des phénomenès physiques. Ce problème peut être lié à la qualité limitée du modèle de résistivité de base utilisé, indiquant ainsi que des efforts plus grands devront être fournis dans le futur pour obtenir des modèles de base de bonne qualité avant de réaliser des expériences dynamiques. Les études décrites dans cette thèse montrent que les méthodes d'ondes EM planes sont très utiles pour caractériser et suivre les variations temporelles du sous-sol sur de larges échelles. Les présentes approches améliorent l'évaluation des modèles obtenus, autant en termes d'incorporation d'informations a priori, qu'en termes de quantification d'incertitudes a posteriori. De plus, les stratégies développées peuvent être appliquées à d'autres méthodes géophysiques, et offrent une grande flexibilité pour l'incorporation d'informations additionnelles lorsqu'elles sont disponibles. -- The efficient use of geothermal systems, the sequestration of CO2 to mitigate climate change, and the prevention of seawater intrusion in coastal aquifers are only some examples that demonstrate the need for novel technologies to monitor subsurface processes from the surface. A main challenge is to assure optimal performance of such technologies at different temporal and spatial scales. Plane-wave electromagnetic (EM) methods are sensitive to subsurface electrical conductivity and consequently to fluid conductivity, fracture connectivity, temperature, and rock mineralogy. These methods have governing equations that are the same over a large range of frequencies, thus allowing to study in an analogous manner processes on scales ranging from few meters close to the surface down to several hundreds of kilometers depth. Unfortunately, they suffer from a significant resolution loss with depth due to the diffusive nature of the electromagnetic fields. Therefore, estimations of subsurface models that use these methods should incorporate a priori information to better constrain the models, and provide appropriate measures of model uncertainty. During my thesis, I have developed approaches to improve the static and dynamic characterization of the subsurface with plane-wave EM methods. In the first part of this thesis, I present a two-dimensional deterministic approach to perform time-lapse inversion of plane-wave EM data. The strategy is based on the incorporation of prior information into the inversion algorithm regarding the expected temporal changes in electrical conductivity. This is done by incorporating a flexible stochastic regularization and constraints regarding the expected ranges of the changes by using Lagrange multipliers. I use non-l2 norms to penalize the model update in order to obtain sharp transitions between regions that experience temporal changes and regions that do not. I also incorporate a time-lapse differencing strategy to remove systematic errors in the time-lapse inversion. This work presents improvements in the characterization of temporal changes with respect to the classical approach of performing separate inversions and computing differences between the models. In the second part of this thesis, I adopt a Bayesian framework and use Markov chain Monte Carlo (MCMC) simulations to quantify model parameter uncertainty in plane-wave EM inversion. For this purpose, I present a two-dimensional pixel-based probabilistic inversion strategy for separate and joint inversions of plane-wave EM and electrical resistivity tomography (ERT) data. I compare the uncertainties of the model parameters when considering different types of prior information on the model structure and different likelihood functions to describe the data errors. The results indicate that model regularization is necessary when dealing with a large number of model parameters because it helps to accelerate the convergence of the chains and leads to more realistic models. These constraints also lead to smaller uncertainty estimates, which imply posterior distributions that do not include the true underlying model in regions where the method has limited sensitivity. This situation can be improved by combining planewave EM methods with complimentary geophysical methods such as ERT. In addition, I show that an appropriate regularization weight and the standard deviation of the data errors can be retrieved by the MCMC inversion. Finally, I evaluate the possibility of characterizing the three-dimensional distribution of an injected water plume by performing three-dimensional time-lapse MCMC inversion of planewave EM data. Since MCMC inversion involves a significant computational burden in high parameter dimensions, I propose a model reduction strategy where the coefficients of a Legendre moment decomposition of the injected water plume and its location are estimated. For this purpose, a base resistivity model is needed which is obtained prior to the time-lapse experiment. A synthetic test shows that the methodology works well when the base resistivity model is correctly characterized. The methodology is also applied to an injection experiment performed in a geothermal system in Australia, and compared to a three-dimensional time-lapse inversion performed within a deterministic framework. The MCMC inversion better constrains the water plumes due to the larger amount of prior information that is included in the algorithm. The conductivity changes needed to explain the time-lapse data are much larger than what is physically possible based on present day understandings. This issue may be related to the base resistivity model used, therefore indicating that more efforts should be given to obtain high-quality base models prior to dynamic experiments. The studies described herein give clear evidence that plane-wave EM methods are useful to characterize and monitor the subsurface at a wide range of scales. The presented approaches contribute to an improved appraisal of the obtained models, both in terms of the incorporation of prior information in the algorithms and the posterior uncertainty quantification. In addition, the developed strategies can be applied to other geophysical methods, and offer great flexibility to incorporate additional information when available.

Role of β-TrCP1, variant and homologue in HIV-1 life cycle

Summary The CD4 molecule plays a key role in AIDS pathogenesis, it is required for entry of the virus into permissive cells and its subsequent down-modulation of the cell surface is a hallmark of HN-1 infected cells. The virus encodes no less than three proteins that participate in this process: Nef, Vpu and Env. Vpu protein interacts with CD4 within the endoplasmic reticulum of infected cells, where it targets CD4 for degradation through the interaction with a cellular protein named ß-TrCP1. This F-box protein functions as the substrate recognition subunit of the SCF ß-Trcr E3 ubiquitin ligase, which normally induce the ubiquitination and subsequent degradation of various proteins such as ß-catenin and IxBa. Mammals possess a homologue of ß-TrCP1, HOS, also named ß-TrCP2 which has a cytoplasmic subcellular distribution. Structural analysis of the ligand-binding domain of both homologues shows striking surface similarities. Both F-box proteins have a redundant role in a number of cellular processes; however the potential role of ß-TrCP2 in HIV-1 infected cells has not been evaluated. In the present study, we assessed the existence of génetic variants of BRTC, encoding ß-TrCP1, and evaluated whether these variants would affect CD4 down-modulation. Additionally, we determined whether ß-TrCP2 shares with its homologue structural and functional properties that would allow it to bind Vpu, modulate CD4 expression, and thus participate in HN-1 pathogenesis. We identified a single nucleotide polymorphism present in the human population with an allelic frequency of 0.03 that leads to the substitution of alanine 507 by a serine. However, we showed by transient transfection in HeLa CD4+ cells that this variant behaves as ß-TrCP1 with respect to CD4 down-modulation. We established transient expression systems in HeLa CD4+ cells to test whether ß-TrCP2 is implicated in Vpu-mediated CD4 down-modulation. We show by coimmunoprecipitation experiments that ß-TrCP2 binds Vpu and is able to induce CD4 down-modulation as efficiently as ß-TrCP1. In two different cell lines, HeLa CD4+ and Jurkat, Vpu-mediated CD4 down-modulation could not be completely reversed through the silencing of endogenous ß-TrCP 1 or ß-TrCP2 individually, but required both genes to be silenced simultaneously. We evaluated the role of ß-TrCP1 and ß-TrCP2 in HIV-1 life cycle using silencing prior to actual viral infection. Both ß-TrCP1 and ß-TrCP2 contributed to CD4 down-modulation during aone-cycle viral infection iri Ghost cells. In addition, the combined silencing of both homologues in the absence of env and nef reversed CD4 down-modulation, showing that ß-TrCP 1 and ß-TrCP2 represent the main and additive effectors of HIV-1 encoded Vpu. In addition, we showed that silencing of ß-TrCPI but not ß-TrCP2 induced a decrease of HIV-1 LTR-driven expression. In a transient transfection system with Tat and a LTR luciferase reporter, both homologues modulated LTR-driven expression. The present study revealed that ß-TrCP2 represents a novel protein participating in HIV-1 cycle and complete comprehension of the complex interplay occurring between the two F-Box will improve our understanding of HIV-1 infection. Résumé La molécule CD4 joue un rôle clef dans la pathogenèse du SIDA ; elle est requise pour l'entrée du virus dans les cellules permissives et la diminution de sa concentration au niveau de la surface cellulaire est une importante caractéristique des cellules infectées par le VIH-1. Le virus encode pas moins de trois protéines qui participent à ce processus Nef, Vpu et Env. La protéine Vpu lie CD4 au niveau du réticulum endoplasmique et induit sa dégradation en interagissant avec une protéine cellulaire nommée ß-TrCP 1. Cette protéine de type F-Box est une sous unité du complexe ubiquitine-ligase E3 SCFß-TrCP. Elle permet la reconnaissance du substrat par le complexe qui induit l'ubiquitination et la subséquente dégradation de diverses protéines cellulaires comme la ß-catenin ou IκBα. Les mammifères possèdent un homologue à ß-TrCP1appelé ß-TrCP2 (ou HOS). L'analyse comparative du domaine permettant la reconnaissance des substrats des deux homologues montre de frappantes similarités. Le rôle de ß-TrCP2 dans le cycle viral du VIH-1 n'a pas encore été évalué. Lors de cette étude, nous avons recherché l'existence de variants génétique de BTRC (codant pour ß-TrCP1) et nous avons évalué si ces variants pourraient affecter la dégradation des molécules CD4 induite par le virus. Nous avons ainsi identifié un polymorphisme présent dans la population humaine avec une fréquence allélique de 0.03 qui consiste en une substitution de l'alanine 507 par une sérine. Nous avons cependant montré par transfection dans des cellules HeLa CD4+ que ce variant se comporte comme ß-TrCP 1 en ce qui concerne la modulation de CD4. De plus, nous avons déterminé si ß-TrCP2 partageait avec son homologue des propriétés structurelles et fonctionnelles qui lui permettraient de lier Vpu, moduler la concentration de CD4 et ainsi prendre part à la pathogenèse du SIDA. Pour ce faire, nous avons établi un système d'expression temporaire dans des cellules HeLa CD4+. Par co-immunoprécipitation, nous avons montré que ß-TrCP2 lie Vpu et est capable d'induire la dégradation de CD4 aussi efficacement que ß-TrCP1. Dans deux différentes lignées cellulaires, HeLa CD4+ et Jurkat, la dégradation de CD4 n'a pu être complètement inhibée par le silencing individuel de ß-TrCP 1 ou ß-TrCP2, mais nécessitait le silencing simultané des 2 gènes. Nous avons évalué le rôle des deux homologues dans le cycle viral du VIH-1 en infectant des cellules Ghost avec le virus après avoir effectué un silencing des deux protéines. Nous avons ainsi montré que ß-TrCP 1 et ß-TrCP2 contribuent de manière additive à la dégradation de CD4 induite par une infection du VIH-1. Le silencing combiné des deux homologues inhiba complètement cette dégradation en l'absence de env et nef, prouvant qu'aucune autre voie ne participe à ce processus: En outre, nous avons montré que le silencing de ß-TrCP 1 mais pas celui de ß-TrCP2 induisait une diminution de l'expression virale sous contrôle du LTR. Nous n'avons cependant pas été en mesure de reconstituer cet effet en exprimant Tat et un gène reporteur sous contrôle du LTR dans des cellules HeLa CD4+. Le présent travail révèle que ß-TrCP2 représente une nouvelle protéine participant dans le cycle viral du VIH-1. Une complète compréhension de l'effet de chacun des deux homologues sur le cycle viral permettra d'améliorer notre compréhension de l'infection par le VIH-1.

Monoclonal antibodies identify a CEA crossreacting antigen of 95 kD (NCA-95) distinct in antigenicity and tissue distribution from the previously described NCA of 55 kD.

During the selection of monoclonal antibodies (MAb) raised against purified carcinoembryonic antigen (CEA), two MAbs were identified which immunoprecipitated a glycoprotein of 95 kD present both in perchloric acid extracts of normal lung and on the surface of normal granulocytes. This antigen was distinct from the previously reported normal glycoprotein crossreacting with CEA (NCA) which had a molecular weight of 55 kD. The difference between the smaller and the larger crossreacting antigens termed NCA-55 and NCA-95, respectively, was demonstrated by SDS-polyacrylamide gel electrophoresis, by elution from Sephadex-G200 and by selective binding to a series of anti-CEA MAb. Out of six MAb which all bound CEA purified from colon carcinoma, three did not react with these two crossreacting antigens, one bound only NCA-95, one reacted only with NCA-55 and one reacted with both NCA-55 and NCA-95. Immunoadsorbent purified preparations of 125I labelled NCA-95 and NCA-55 were found useful for the screening of new anti-CEA MAb. In addition, when tested on frozen sections of colon carcinoma, normal spleen, normal lung and pancreas, each type of MAb gave a clearly different pattern of reactivity. The three anti-CEA MAb which did not bind any of the crossreacting antigens stained only the colon carcinoma cells; the MAb binding to either one of the two types of NCA gave a similar pattern of reactivity both on colon carcinoma cells and on granulocytes. However, on normal lung and pancreas, the MAb binding NCA-55 stained granulocytes as well as bronchiolar and alveolar epithelial cells in lung and inter- and intra-lobular duct epithelial cells in pancreas, whereas the MAb binding only NCA-95 stained only the granulocytes. Thus, the newly identified NCA-95 appears to differ from NCA-55 not only in terms of molecular size and antigenicity but also by the fact that in normal lung and pancreas it is found in granulocytes but not in epithelial cells.

Ambulatory assessment of 3D ground reaction force using plantar pressure distribution.

This study aimed to use the plantar pressure insole for estimating the three-dimensional ground reaction force (GRF) as well as the frictional torque (T(F)) during walking. Eleven subjects, six healthy and five patients with ankle disease participated in the study while wearing pressure insoles during several walking trials on a force-plate. The plantar pressure distribution was analyzed and 10 principal components of 24 regional pressure values with the stance time percentage (STP) were considered for GRF and T(F) estimation. Both linear and non-linear approximators were used for estimating the GRF and T(F) based on two learning strategies using intra-subject and inter-subjects data. The RMS error and the correlation coefficient between the approximators and the actual patterns obtained from force-plate were calculated. Our results showed better performance for non-linear approximation especially when the STP was considered as input. The least errors were observed for vertical force (4%) and anterior-posterior force (7.3%), while the medial-lateral force (11.3%) and frictional torque (14.7%) had higher errors. The result obtained for the patients showed higher error; nevertheless, when the data of the same patient were used for learning, the results were improved and in general slight differences with healthy subjects were observed. In conclusion, this study showed that ambulatory pressure insole with data normalization, an optimal choice of inputs and a well-trained nonlinear mapping function can estimate efficiently the three-dimensional ground reaction force and frictional torque in consecutive gait cycle without requiring a force-plate.

Are niche-based species distribution models transferable in space?

Aim To assess the geographical transferability of niche-based species distribution models fitted with two modelling techniques. Location Two distinct geographical study areas in Switzerland and Austria, in the subalpine and alpine belts. Methods Generalized linear and generalized additive models (GLM and GAM) with a binomial probability distribution and a logit link were fitted for 54 plant species, based on topoclimatic predictor variables. These models were then evaluated quantitatively and used for spatially explicit predictions within (internal evaluation and prediction) and between (external evaluation and prediction) the two regions. Comparisons of evaluations and spatial predictions between regions and models were conducted in order to test if species and methods meet the criteria of full transferability. By full transferability, we mean that: (1) the internal evaluation of models fitted in region A and B must be similar; (2) a model fitted in region A must at least retain a comparable external evaluation when projected into region B, and vice-versa; and (3) internal and external spatial predictions have to match within both regions. Results The measures of model fit are, on average, 24% higher for GAMs than for GLMs in both regions. However, the differences between internal and external evaluations (AUC coefficient) are also higher for GAMs than for GLMs (a difference of 30% for models fitted in Switzerland and 54% for models fitted in Austria). Transferability, as measured with the AUC evaluation, fails for 68% of the species in Switzerland and 55% in Austria for GLMs (respectively for 67% and 53% of the species for GAMs). For both GAMs and GLMs, the agreement between internal and external predictions is rather weak on average (Kulczynski's coefficient in the range 0.3-0.4), but varies widely among individual species. The dominant pattern is an asymmetrical transferability between the two study regions (a mean decrease of 20% for the AUC coefficient when the models are transferred from Switzerland and 13% when they are transferred from Austria). Main conclusions The large inter-specific variability observed among the 54 study species underlines the need to consider more than a few species to test properly the transferability of species distribution models. The pronounced asymmetry in transferability between the two study regions may be due to peculiarities of these regions, such as differences in the ranges of environmental predictors or the varied impact of land-use history, or to species-specific reasons like differential phenotypic plasticity, existence of ecotypes or varied dependence on biotic interactions that are not properly incorporated into niche-based models. The lower variation between internal and external evaluation of GLMs compared to GAMs further suggests that overfitting may reduce transferability. Overall, a limited geographical transferability calls for caution when projecting niche-based models for assessing the fate of species in future environments.

Conservation of the Eurasian Lynx (Lynx lynx) in a fragmented landscape - habitat models, dispersal and potential distribution

Abstract: The expansion of a recovering population - whether re-introduced or spontaneously returning - is shaped by (i) biological (intrinsic) factors such as the land tenure system or dispersal, (ii) the distribution and availability of resources (e.g. prey), (iii) habitat and landscape features, and (iv) human attitudes and activities. In order to develop efficient conservation and recovery strategies, we need to understand all these factors and to predict the potential distribution and explore ways to reach it. An increased number of lynx in the north-western Swiss Alps in the nineties lead to a new controversy about the return of this cat. When the large carnivores were given legal protection in many European countries, most organizations and individuals promoting their protection did not foresee the consequences. Management plans describing how to handle conflicts with large predators are needed to find a balance between "overabundance" and extinction. Wildlife and conservation biologists need to evaluate the various threats confronting populations so that adequate management decisions can be taken. I developed a GIS probability model for the lynx, based on habitat information and radio-telemetry data from the Swiss Jura Mountains, in order to predict the potential distribution of the lynx in this mountain range, which is presently only partly occupied by lynx. Three of the 18 variables tested for each square kilometre describing land use, vegetation, and topography, qualified to predict the probability of lynx presence. The resulting map was evaluated with data from dispersing subadult lynx. Young lynx that were not able to establish home ranges in what was identified as good lynx habitat did not survive their first year of independence, whereas the only one that died in good lynx habitat was illegally killed. Radio-telemetry fixes are often used as input data to calibrate habitat models. Radio-telemetry is the only way to gather accurate and unbiased data on habitat use of elusive larger terrestrial mammals. However, it is time consuming and expensive, and can therefore only be applied in limited areas. Habitat models extrapolated over large areas can in turn be problematic, as habitat characteristics and availability may change from one area to the other. I analysed the predictive power of Ecological Niche Factor Analysis (ENFA) in Switzerland with the lynx as focal species. According to my results, the optimal sampling strategy to predict species distribution in an Alpine area lacking available data would be to pool presence cells from contrasted regions (Jura Mountains, Alps), whereas in regions with a low ecological variance (Jura Mountains), only local presence cells should be used for the calibration of the model. Dispersal influences the dynamics and persistence of populations, the distribution and abundance of species, and gives the communities and ecosystems their characteristic texture in space and time. Between 1988 and 2001, the spatio-temporal behaviour of subadult Eurasian lynx in two re-introduced populations in Switzerland was studied, based on 39 juvenile lynx of which 24 were radio-tagged to understand the factors influencing dispersal. Subadults become independent from their mothers at the age of 8-11 months. No sex bias neither in the dispersal rate nor in the distance moved was detected. Lynx are conservative dispersers, compared to bear and wolf, and settled within or close to known lynx occurrences. Dispersal distances reached in the high lynx density population - shorter than those reported in other Eurasian lynx studies - are limited by habitat restriction hindering connections with neighbouring metapopulations. I postulated that high lynx density would lead to an expansion of the population and validated my predictions with data from the north-western Swiss Alps where about 1995 a strong increase in lynx abundance took place. The general hypothesis that high population density will foster the expansion of the population was not confirmed. This has consequences for the re-introduction and recovery of carnivores in a fragmented landscape. To establish a strong source population in one place might not be an optimal strategy. Rather, population nuclei should be founded in several neighbouring patches. Exchange between established neighbouring subpopulations will later on take place, as adult lynx show a higher propensity to cross barriers than subadults. To estimate the potential population size of the lynx in the Jura Mountains and to assess possible corridors between this population and adjacent areas, I adapted a habitat probability model for lynx distribution in the Jura Mountains with new environmental data and extrapolated it over the entire mountain range. The model predicts a breeding population ranging from 74-101 individuals and from 51-79 individuals when continuous habitat patches < 50 km2 are disregarded. The Jura Mountains could once be part of a metapopulation, as potential corridors exist to the adjoining areas (Alps, Vosges Mountains, and Black Forest). Monitoring of the population size, spatial expansion, and the genetic surveillance in the Jura Mountains must be continued, as the status of the population is still critical. ENFA was used to predict the potential distribution of lynx in the Alps. The resulting model divided the Alps into 37 suitable habitat patches ranging from 50 to 18,711 km2, covering a total area of about 93,600 km2. When using the range of lynx densities found in field studies in Switzerland, the Alps could host a population of 961 to 1,827 residents. The results of the cost-distance analysis revealed that all patches were within the reach of dispersing lynx, as the connection costs were in the range of dispersal cost of radio-tagged subadult lynx moving through unfavorable habitat. Thus, the whole Alps could once be considered as a metapopulation. But experience suggests that only few disperser will cross unsuitable areas and barriers. This low migration rate may seldom allow the spontaneous foundation of new populations in unsettled areas. As an alternative to natural dispersal, artificial transfer of individuals across the barriers should be considered. Wildlife biologists can play a crucial role in developing adaptive management experiments to help managers learning by trial. The case of the lynx in Switzerland is a good example of a fruitful cooperation between wildlife biologists, managers, decision makers and politician in an adaptive management process. This cooperation resulted in a Lynx Management Plan which was implemented in 2000 and updated in 2004 to give the cantons directives on how to handle lynx-related problems. This plan was put into practice e.g. in regard to translocation of lynx into unsettled areas. Résumé: L'expansion d'une population en phase de recolonisation, qu'elle soit issue de réintroductions ou d'un retour naturel dépend 1) de facteurs biologiques tels que le système social et le mode de dispersion, 2) de la distribution et la disponibilité des ressources (proies), 3) de l'habitat et des éléments du paysage, 4) de l'acceptation de l'espèce par la population locale et des activités humaines. Afin de pouvoir développer des stratégies efficaces de conservation et de favoriser la recolonisation, chacun de ces facteurs doit être pris en compte. En plus, la distribution potentielle de l'espèce doit pouvoir être déterminée et enfin, toutes les possibilités pour atteindre les objectifs, examinées. La phase de haute densité que la population de lynx a connue dans les années nonante dans le nord-ouest des Alpes suisses a donné lieu à une controverse assez vive. La protection du lynx dans de nombreux pays européens, promue par différentes organisations, a entraîné des conséquences inattendues; ces dernières montrent que tout plan de gestion doit impérativement indiquer des pistes quant à la manière de gérer les conflits, tout en trouvant un équilibre entre l'extinction et la surabondance de l'espèce. Les biologistes de la conservation et de la faune sauvage doivent pour cela évaluer les différents risques encourus par les populations de lynx, afin de pouvoir rapidement prendre les meilleuresmdécisions de gestion. Un modèle d'habitat pour le lynx, basé sur des caractéristiques de l'habitat et des données radio télémétriques collectées dans la chaîne du Jura, a été élaboré afin de prédire la distribution potentielle dans cette région, qui n'est que partiellement occupée par l'espèce. Trois des 18 variables testées, décrivant pour chaque kilomètre carré l'utilisation du sol, la végétation ainsi que la topographie, ont été retenues pour déterminer la probabilité de présence du lynx. La carte qui en résulte a été comparée aux données télémétriques de lynx subadultes en phase de dispersion. Les jeunes qui n'ont pas pu établir leur domaine vital dans l'habitat favorable prédit par le modèle n'ont pas survécu leur première année d'indépendance alors que le seul individu qui est mort dans l'habitat favorable a été braconné. Les données radio-télémétriques sont souvent utilisées pour l'étalonnage de modèles d'habitat. C'est un des seuls moyens à disposition qui permette de récolter des données non biaisées et précises sur l'occupation de l'habitat par des mammifères terrestres aux moeurs discrètes. Mais ces méthodes de- mandent un important investissement en moyens financiers et en temps et peuvent, de ce fait, n'être appliquées qu'à des zones limitées. Les modèles d'habitat sont ainsi souvent extrapolés à de grandes surfaces malgré le risque d'imprécision, qui résulte des variations des caractéristiques et de la disponibilité de l'habitat d'une zone à l'autre. Le pouvoir de prédiction de l'Analyse Ecologique de la Niche (AEN) dans les zones où les données de présence n'ont pas été prises en compte dans le calibrage du modèle a été analysée dans le cas du lynx en Suisse. D'après les résultats obtenus, la meilleure mé- thode pour prédire la distribution du lynx dans une zone alpine dépourvue d'indices de présence est de combiner des données provenant de régions contrastées (Alpes, Jura). Par contre, seules les données sur la présence locale de l'espèce doivent être utilisées pour les zones présentant une faible variance écologique tel que le Jura. La dispersion influence la dynamique et la stabilité des populations, la distribution et l'abondance des espèces et détermine les caractéristiques spatiales et temporelles des communautés vivantes et des écosystèmes. Entre 1988 et 2001, le comportement spatio-temporel de lynx eurasiens subadultes de deux populations réintroduites en Suisse a été étudié, basé sur le suivi de 39 individus juvéniles dont 24 étaient munis d'un collier émetteur, afin de déterminer les facteurs qui influencent la dispersion. Les subadultes se sont séparés de leur mère à l'âge de 8 à 11 mois. Le sexe n'a pas eu d'influence sur le nombre d'individus ayant dispersés et la distance parcourue au cours de la dispersion. Comparé à l'ours et au loup, le lynx reste très modéré dans ses mouvements de dispersion. Tous les individus ayant dispersés se sont établis à proximité ou dans des zones déjà occupées par des lynx. Les distances parcourues lors de la dispersion ont été plus courtes pour la population en phase de haute densité que celles relevées par les autres études de dispersion du lynx eurasien. Les zones d'habitat peu favorables et les barrières qui interrompent la connectivité entre les populations sont les principales entraves aux déplacements, lors de la dispersion. Dans un premier temps, nous avons fait l'hypothèse que les phases de haute densité favorisaient l'expansion des populations. Mais cette hypothèse a été infirmée par les résultats issus du suivi des lynx réalisé dans le nord-ouest des Alpes, où la population connaissait une phase de haute densité depuis 1995. Ce constat est important pour la conservation d'une population de carnivores dans un habitat fragmenté. Ainsi, instaurer une forte population source à un seul endroit n'est pas forcément la stratégie la plus judicieuse. Il est préférable d'établir des noyaux de populations dans des régions voisines où l'habitat est favorable. Des échanges entre des populations avoisinantes pourront avoir lieu par la suite car les lynx adultes sont plus enclins à franchir les barrières qui entravent leurs déplacements que les individus subadultes. Afin d'estimer la taille de la population de lynx dans le Jura et de déterminer les corridors potentiels entre cette région et les zones avoisinantes, un modèle d'habitat a été utilisé, basé sur un nouveau jeu de variables environnementales et extrapolé à l'ensemble du Jura. Le modèle prédit une population reproductrice de 74 à 101 individus et de 51 à 79 individus lorsque les surfaces d'habitat d'un seul tenant de moins de 50 km2 sont soustraites. Comme des corridors potentiels existent effectivement entre le Jura et les régions avoisinantes (Alpes, Vosges, et Forêt Noire), le Jura pourrait faire partie à l'avenir d'une métapopulation, lorsque les zones avoisinantes seront colonisées par l'espèce. La surveillance de la taille de la population, de son expansion spatiale et de sa structure génétique doit être maintenue car le statut de cette population est encore critique. L'AEN a également été utilisée pour prédire l'habitat favorable du lynx dans les Alpes. Le modèle qui en résulte divise les Alpes en 37 sous-unités d'habitat favorable dont la surface varie de 50 à 18'711 km2, pour une superficie totale de 93'600 km2. En utilisant le spectre des densités observées dans les études radio-télémétriques effectuées en Suisse, les Alpes pourraient accueillir une population de lynx résidents variant de 961 à 1'827 individus. Les résultats des analyses de connectivité montrent que les sous-unités d'habitat favorable se situent à des distances telles que le coût de la dispersion pour l'espèce est admissible. L'ensemble des Alpes pourrait donc un jour former une métapopulation. Mais l'expérience montre que très peu d'individus traverseront des habitats peu favorables et des barrières au cours de leur dispersion. Ce faible taux de migration rendra difficile toute nouvelle implantation de populations dans des zones inoccupées. Une solution alternative existe cependant : transférer artificiellement des individus d'une zone à l'autre. Les biologistes spécialistes de la faune sauvage peuvent jouer un rôle important et complémentaire pour les gestionnaires de la faune, en les aidant à mener des expériences de gestion par essai. Le cas du lynx en Suisse est un bel exemple d'une collaboration fructueuse entre biologistes de la faune sauvage, gestionnaires, organes décisionnaires et politiciens. Cette coopération a permis l'élaboration du Concept Lynx Suisse qui est entré en vigueur en 2000 et remis à jour en 2004. Ce plan donne des directives aux cantons pour appréhender la problématique du lynx. Il y a déjà eu des applications concrètes sur le terrain, notamment par des translocations d'individus dans des zones encore inoccupées.