Quantitative integration of geophysical and hydraulic data : from the local towards the regional scale range

Des progrès significatifs ont été réalisés dans le domaine de l'intégration quantitative des données géophysique et hydrologique l'échelle locale. Cependant, l'extension à de plus grandes échelles des approches correspondantes constitue encore un défi majeur. Il est néanmoins extrêmement important de relever ce défi pour développer des modèles fiables de flux des eaux souterraines et de transport de contaminant. Pour résoudre ce problème, j'ai développé une technique d'intégration des données hydrogéophysiques basée sur une procédure bayésienne de simulation séquentielle en deux étapes. Cette procédure vise des problèmes à plus grande échelle. L'objectif est de simuler la distribution d'un paramètre hydraulique cible à partir, d'une part, de mesures d'un paramètre géophysique pertinent qui couvrent l'espace de manière exhaustive, mais avec une faible résolution (spatiale) et, d'autre part, de mesures locales de très haute résolution des mêmes paramètres géophysique et hydraulique. Pour cela, mon algorithme lie dans un premier temps les données géophysiques de faible et de haute résolution à travers une procédure de réduction déchelle. Les données géophysiques régionales réduites sont ensuite reliées au champ du paramètre hydraulique à haute résolution. J'illustre d'abord l'application de cette nouvelle approche dintégration des données à une base de données synthétiques réaliste. Celle-ci est constituée de mesures de conductivité hydraulique et électrique de haute résolution réalisées dans les mêmes forages ainsi que destimations des conductivités électriques obtenues à partir de mesures de tomographic de résistivité électrique (ERT) sur l'ensemble de l'espace. Ces dernières mesures ont une faible résolution spatiale. La viabilité globale de cette méthode est testée en effectuant les simulations de flux et de transport au travers du modèle original du champ de conductivité hydraulique ainsi que du modèle simulé. Les simulations sont alors comparées. Les résultats obtenus indiquent que la procédure dintégration des données proposée permet d'obtenir des estimations de la conductivité en adéquation avec la structure à grande échelle ainsi que des predictions fiables des caractéristiques de transports sur des distances de moyenne à grande échelle. Les résultats correspondant au scénario de terrain indiquent que l'approche d'intégration des données nouvellement mise au point est capable d'appréhender correctement les hétérogénéitées à petite échelle aussi bien que les tendances à gande échelle du champ hydraulique prévalent. Les résultats montrent également une flexibilté remarquable et une robustesse de cette nouvelle approche dintégration des données. De ce fait, elle est susceptible d'être appliquée à un large éventail de données géophysiques et hydrologiques, à toutes les gammes déchelles. Dans la deuxième partie de ma thèse, j'évalue en détail la viabilité du réechantillonnage geostatique séquentiel comme mécanisme de proposition pour les méthodes Markov Chain Monte Carlo (MCMC) appliquées à des probmes inverses géophysiques et hydrologiques de grande dimension . L'objectif est de permettre une quantification plus précise et plus réaliste des incertitudes associées aux modèles obtenus. En considérant une série dexemples de tomographic radar puits à puits, j'étudie deux classes de stratégies de rééchantillonnage spatial en considérant leur habilité à générer efficacement et précisément des réalisations de la distribution postérieure bayésienne. Les résultats obtenus montrent que, malgré sa popularité, le réechantillonnage séquentiel est plutôt inefficace à générer des échantillons postérieurs indépendants pour des études de cas synthétiques réalistes, notamment pour le cas assez communs et importants où il existe de fortes corrélations spatiales entre le modèle et les paramètres. Pour résoudre ce problème, j'ai développé un nouvelle approche de perturbation basée sur une déformation progressive. Cette approche est flexible en ce qui concerne le nombre de paramètres du modèle et lintensité de la perturbation. Par rapport au rééchantillonage séquentiel, cette nouvelle approche s'avère être très efficace pour diminuer le nombre requis d'itérations pour générer des échantillons indépendants à partir de la distribution postérieure bayésienne. - Significant progress has been made with regard to the quantitative integration of geophysical and hydrological data at the local scale. However, extending corresponding approaches beyond the local scale still represents a major challenge, yet is critically important for the development of reliable groundwater flow and contaminant transport models. To address this issue, I have developed a hydrogeophysical data integration technique based on a two-step Bayesian sequential simulation procedure that is specifically targeted towards larger-scale problems. The objective is to simulate the distribution of a target hydraulic parameter based on spatially exhaustive, but poorly resolved, measurements of a pertinent geophysical parameter and locally highly resolved, but spatially sparse, measurements of the considered geophysical and hydraulic parameters. To this end, my algorithm links the low- and high-resolution geophysical data via a downscaling procedure before relating the downscaled regional-scale geophysical data to the high-resolution hydraulic parameter field. I first illustrate the application of this novel data integration approach to a realistic synthetic database consisting of collocated high-resolution borehole measurements of the hydraulic and electrical conductivities and spatially exhaustive, low-resolution electrical conductivity estimates obtained from electrical resistivity tomography (ERT). The overall viability of this method is tested and verified by performing and comparing flow and transport simulations through the original and simulated hydraulic conductivity fields. The corresponding results indicate that the proposed data integration procedure does indeed allow for obtaining faithful estimates of the larger-scale hydraulic conductivity structure and reliable predictions of the transport characteristics over medium- to regional-scale distances. The approach is then applied to a corresponding field scenario consisting of collocated high- resolution measurements of the electrical conductivity, as measured using a cone penetrometer testing (CPT) system, and the hydraulic conductivity, as estimated from electromagnetic flowmeter and slug test measurements, in combination with spatially exhaustive low-resolution electrical conductivity estimates obtained from surface-based electrical resistivity tomography (ERT). The corresponding results indicate that the newly developed data integration approach is indeed capable of adequately capturing both the small-scale heterogeneity as well as the larger-scale trend of the prevailing hydraulic conductivity field. The results also indicate that this novel data integration approach is remarkably flexible and robust and hence can be expected to be applicable to a wide range of geophysical and hydrological data at all scale ranges. In the second part of my thesis, I evaluate in detail the viability of sequential geostatistical resampling as a proposal mechanism for Markov Chain Monte Carlo (MCMC) methods applied to high-dimensional geophysical and hydrological inverse problems in order to allow for a more accurate and realistic quantification of the uncertainty associated with the thus inferred models. Focusing on a series of pertinent crosshole georadar tomographic examples, I investigated two classes of geostatistical resampling strategies with regard to their ability to efficiently and accurately generate independent realizations from the Bayesian posterior distribution. The corresponding results indicate that, despite its popularity, sequential resampling is rather inefficient at drawing independent posterior samples for realistic synthetic case studies, notably for the practically common and important scenario of pronounced spatial correlation between model parameters. To address this issue, I have developed a new gradual-deformation-based perturbation approach, which is flexible with regard to the number of model parameters as well as the perturbation strength. Compared to sequential resampling, this newly proposed approach was proven to be highly effective in decreasing the number of iterations required for drawing independent samples from the Bayesian posterior distribution.

A novel role for proline- and acid-rich basic region leucine zipper (PAR bZIP) proteins in the transcriptional regulation of a BH3-only proapoptotic gene.

Proline- and acid-rich (PAR) basic region leucine zipper (bZIP) proteins thyrotroph embryonic factor (TEF), D-site-binding protein (DBP), and hepatic leukemia factor have been involved in neurotransmitter homeostasis and amino acid metabolism. Here we demonstrate a novel role for these proteins in the transcriptional control of a BH3-only gene. PAR bZIP proteins are able to transactivate the promoter of bcl-gS. This promoter is particularly responsive to TEF activation and is silenced by NFIL3, a repressor that shares the consensus binding site with PAR bZIP proteins. Consistently, transfection of TEF induces the expression of endogenous bcl-gS in cancer cells, and this induction is independent of p53. A naturally occurring variant of DBP (tDBP), lacking the transactivation domain, has been identified and shown to impede the formation of active TEF dimers in a competitive manner and to reduce the TEF-dependent induction of bcl-gS. Of note, treatment of cancer cells with etoposide induces TEF activation and promotes the expression of bcl-gS. Furthermore, blockade of bcl-gS or TEF expression by a small interfering RNA strategy or transfection with tDBP significantly reduces the etoposide-mediated apoptotic cell death. These findings represent the first described role for PAR bZIP proteins in the regulation of a gene involved in the execution of apoptosis.

A novel tool to analyze MRI recurrence patterns in glioblastoma.

At least 10% of glioblastoma relapses occur at distant and even contralateral locations. This disseminated growth limits surgical intervention and contributes to neurological morbidity. Preclinical data pointed toward a role for temozolomide (TMZ) in reducing radiotherapy-induced glioma cell invasiveness. Our objective was to develop and validate a new analysis tool of MRI data to examine the clinical recurrence pattern of glioblastomas. MRIcro software was used to map the location and extent of initial preoperative and recurrent tumors on MRI of 63 patients in the European Organisation for Research and Treatment of Cancer (EORTC) 26981/22981/National Cancer Institute of Canada (NCIC) CE.3 study into the same stereotaxic space. This allowed us to examine changes of site and distance between the initial and the recurrent tumor on the group level. Thirty of the 63 patients were treated using radiotherapy, while the other patients completed a radiotherapy-plus-TMZ treatment. Baseline characteristics (median age, KPS) and outcome data (progression-free survival, overall survival) of the patients included in this analysis resemble those of the general study cohort. The patient groups did not differ in the promoter methylation status of methyl guanine methyltransferase (MGMT). Overall frequency of distant recurrences was 20%. Analysis of recurrence patterns revealed no difference between the groups in the size of the recurrent tumor or in the differential effect on the distance of the recurrences from the preoperative tumor location. The data show the feasibility of groupwise recurrence pattern analysis. An effect of TMZ treatment on the recurrence pattern in the EORTC 26981/22981/NCIC CE.3 study could not be demonstrated.

Regulation of the DNA-binding and transcriptional activities of Xenopus laevis NFI-X by a novel C-terminal domain.

The nuclear factor I (NFI) family consists of sequence-specific DNA-binding proteins that activate both transcription and adenovirus DNA replication. We have characterized three new members of the NFI family that belong to the Xenopus laevis NFI-X subtype and differ in their C-termini. We show that these polypeptides can activate transcription in HeLa and Drosophila Schneider line 2 cells, using an activation domain that is subdivided into adjacent variable and subtype-specific domains each having independent activation properties in chimeric proteins. Together, these two domains constitute the full NFI-X transactivation potential. In addition, we find that the X. laevis NFI-X proteins are capable of activating adenovirus DNA replication through their conserved N-terminal DNA-binding domains. Surprisingly, their in vitro DNA-binding activities are specifically inhibited by a novel repressor domain contained within the C-terminal part, while the dimerization and replication functions per se are not affected. However, inhibition of DNA-binding activity in vitro is relieved within the cell, as transcriptional activation occurs irrespective of the presence of the repressor domain. Moreover, the region comprising the repressor domain participates in transactivation. Mechanisms that may allow the relief of DNA-binding inhibition in vivo and trigger transcriptional activation are discussed.

Avian erythrocytes have functional mitochondria, opening novel perspectives for birds as animal models in the study of ageing.

BACKGROUND: In contrast to mammalian erythrocytes, which have lost their nucleus and mitochondria during maturation, the erythrocytes of almost all other vertebrate species are nucleated throughout their lifespan. Little research has been done however to test for the presence and functionality of mitochondria in these cells, especially for birds. Here, we investigated those two points in erythrocytes of one common avian model: the zebra finch (Taeniopygia guttata). RESULTS: Transmission electron microscopy showed the presence of mitochondria in erythrocytes of this small passerine bird, especially after removal of haemoglobin interferences. High-resolution respirometry revealed increased or decreased rates of oxygen consumption by erythrocytes in response to the addition of respiratory chain substrates or inhibitors, respectively. Fluorometric assays confirmed the production of mitochondrial superoxide by avian erythrocytes. Interestingly, measurements of plasmatic oxidative markers indicated lower oxidative stress in blood of the zebra finch compared to a size-matched mammalian model, the mouse. CONCLUSIONS: Altogether, those findings demonstrate that avian erythrocytes possess functional mitochondria in terms of respiratory activities and reactive oxygen species (ROS) production. Interestingly, since blood oxidative stress was lower for our avian model compared to a size-matched mammalian, our results also challenge the idea that mitochondrial ROS production could have been one actor leading to this loss during the course of evolution. Opportunities to assess mitochondrial functioning in avian erythrocytes open new perspectives in the use of birds as models for longitudinal studies of ageing via lifelong blood sampling of the same subjects.

Where will conflicts between alien and rare species occur after climate and land-use change? A test with a novel combined modelling approach

Protecting native biodiversity against alien invasive species requires powerful methods to anticipate these invasions and to protect native species assumed to be at risk. Here, we describe how species distribution models (SDMs) can be used to identify areas predicted as suitable for rare native species and also predicted as highly susceptible to invasion by alien species, at present and under future climate and land-use scenarios. To assess the condition and dynamics of such conflicts, we developed a combined predictive modelling (CPM) approach, which predicts species distributions by combining two SDMs fitted using subsets of predictors classified as acting at either regional or local scales. We illustrate the CPM approach for an alien invader and a rare species associated to similar habitats in northwest Portugal. Combined models predict a wider variety of potential species responses, providing more informative projections of species distributions and future dynamics than traditional, non-combined models. They also provide more informative insight regarding current and future rare-invasive conflict areas. For our studied species, conflict areas of highest conservation relevance are predicted to decrease over the next decade, supporting previous reports that some invasive species may contract their geographic range and impact due to climate change. More generally, our results highlight the more informative character of the combined approach to address practical issues in conservation and management programs, especially those aimed at mitigating the impact of invasive plants, land-use and climate changes in sensitive regions

The novel hydroxylamine derivative NG-094 suppresses polyglutamine protein toxicity in Caenorhabditis elegans.

Aggregation-prone polyglutamine (polyQ) expansion proteins cause several neurodegenerative disorders, including Huntington disease. The pharmacological activation of cellular stress responses could be a new strategy to combat protein conformational diseases. Hydroxylamine derivatives act as co-inducers of heat-shock proteins (HSPs) and can enhance HSP expression in diseased cells, without significant adverse effects. Here, we used Caenorhabditis elegans expressing polyQ expansions with 35 glutamines fused to the yellow fluorescent protein (Q35-YFP) in body wall muscle cells as a model system to investigate the effects of treatment with a novel hydroxylamine derivative, NG-094, on the progression of polyQ diseases. NG-094 significantly ameliorated polyQ-mediated animal paralysis, reduced the number of Q35-YFP aggregates and delayed polyQ-dependent acceleration of aging. Micromolar concentrations of NG-094 in animal tissues with only marginal effects on the nematode fitness sufficed to confer protection against polyQ proteotoxicity, even when the drug was administered after disease onset. NG-094 did not reduce insulin/insulin-like growth factor 1-like signaling, but conferred cytoprotection by a mechanism involving the heat-shock transcription factor HSF-1 that potentiated the expression of stress-inducible HSPs. NG-094 is thus a promising candidate for tests on mammalian models of polyQ and other protein conformational diseases.

Novel relationship between Vegf expression and retinal pigmented epithelium permeability following light damage in the mouse retina

Purpose:In the retina, the balance between pro- and anti-angiogenic factors is critical for angiogenesis control but is also involved in cell survival and maintenance. For instance, the anti-angiogenic factor PEDF is neuroprotective for photoreceptors (PRs) in models of retinal degeneration. We previously reported upregulation of VEGF (24h to 48h post lesion) in the light-damage (LD) model. Furthermore, systemic delivery of PEDF, as well as lentiviral gene transfer of an anti-VEGF antibody rescue PRs from cell death. Studies in vitro show that VEGF induces retinal endothelial cells apoptosis via the alteration of the Akt1/p38 MAPK signalling pathway under hypoxic conditions. Thus, in this study, we investigate the effect of high levels of VEGF on retinal pigmented epithelium (RPE) permeability and molecular targets expression after light-induced PR degeneration. Methods:To characterize the action of VEGF in the retina during the course of LD, we exposed adult Balb/c mice to 5'000 lux for 1h, and we collected neural retinas and eye-cups (containing RPE) at different time points after the LD. We analysed protein expression by Elisa and Western blotting. In order to study RPE cell permeability after the LD we stained β-catenin on flat mounted RPE. Results:In the neural retina, preliminary results indicate that high levels of VEGF induce a significant upregulation of VEGF receptor 2, whereas VEGF receptor 1 expression is decreased. Concomitantly with VEGF upregulation, LD increases the Src phosphorylation between 24h to 48h. Furthermore, we observe that β-catenin translocates to the cytoplasm of RPE cells between 24h to 36h after the lesion, indicating an increase on the RPE permeability, which could contribute indirectly to the deleterious effect of VEGF observed during light-induced PR apoptosis. Conclusions:This study further involves VEGF in LD and highlights the prime importance of angiogenic factor balance for PR survival. Our results suggest that PR apoptosis is augmented by RPE cell permeability, which may induce high level of VEGF and could be deleterious. The specific action of RPE permeability on PR survival and the role of Src in the retina are under investigation.

A comprehensive characterization of genome-wide copy number aberrations in colorectal cancer reveals novel oncogenes and patterns of alterations.

To develop a comprehensive overview of copy number aberrations (CNAs) in stage-II/III colorectal cancer (CRC), we characterized 302 tumors from the PETACC-3 clinical trial. Microsatellite-stable (MSS) samples (n = 269) had 66 minimal common CNA regions, with frequent gains on 20 q (72.5%), 7 (41.8%), 8 q (33.1%) and 13 q (51.0%) and losses on 18 (58.6%), 4 q (26%) and 21 q (21.6%). MSS tumors have significantly more CNAs than microsatellite-instable (MSI) tumors: within the MSI tumors a novel deletion of the tumor suppressor WWOX at 16 q23.1 was identified (p<0.01). Focal aberrations identified by the GISTIC method confirmed amplifications of oncogenes including EGFR, ERBB2, CCND1, MET, and MYC, and deletions of tumor suppressors including TP53, APC, and SMAD4, and gene expression was highly concordant with copy number aberration for these genes. Novel amplicons included putative oncogenes such as WNK1 and HNF4A, which also showed high concordance between copy number and expression. Survival analysis associated a specific patient segment featured by chromosome 20 q gains to an improved overall survival, which might be due to higher expression of genes such as EEF1B2 and PTK6. The CNA clustering also grouped tumors characterized by a poor prognosis BRAF-mutant-like signature derived from mRNA data from this cohort. We further revealed non-random correlation between CNAs among unlinked loci, including positive correlation between 20 q gain and 8 q gain, and 20 q gain and chromosome 18 loss, consistent with co-selection of these CNAs. These results reinforce the non-random nature of somatic CNAs in stage-II/III CRC and highlight loci and genes that may play an important role in driving the development and outcome of this disease.

Ampakine CX546 bolsters energetic response of astrocytes: a novel target for cognitive-enhancing drugs acting as alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptor modulators.

Glutamate was previously shown to enhance aerobic glycolysis i.e. increase glucose utilization and lactate production with no change in oxygen levels, in mouse cortical astrocytes by a mechanism involving glutamate uptake. It is reported here that a similar response is produced in both hippocampal and cerebellar astrocytes. Application of the cognitive-enhancing drug CX546 promoted further enhancement of glucose utilization by astrocytes from each brain area following glutamate exposure. alpha-Amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors represent the purported molecular target of cognitive-enhancing drugs such as CX546, and the presence of AMPA receptor subunits GluR1-4 was evidenced in astrocytes from all three regions by immunocytochemistry. AMPA itself did not stimulate aerobic glycolysis, but in the presence of CX546, a strong enhancement of glucose utilization and lactate production was obtained in cortical, hippocampal and cerebellar astrocytes. The effect of CX546 was concentration-dependent, with an EC(50) of 93.2 microm in cortical astrocytes. AMPA-induced glucose utilization in the presence of CX546 was prevented by the AMPA receptor antagonist 6-cyano-7-nitroquinoxaline-2,3-dione (CNQX) and the negative modulator GYKI 52466. In addition, the metabolic effect of CX546 in the presence of AMPA was mimicked by the AMPA receptor modulator cyclothiazide. Our data suggest that astrocyte energetics represents a novel target for cognitive-enhancing drugs acting as AMPA receptor modulators.

Nineteen additional unpredicted transcripts from human chromosome 21.

The identification of all human chromosome 21 (HC21) genes is a necessary step in understanding the molecular pathogenesis of trisomy 21 (Down syndrome). The first analysis of the sequence of 21q included 127 previously characterized genes and predicted an additional 98 novel anonymous genes. Recently we evaluated the quality of this annotation by characterizing a set of HC21 open reading frames (C21orfs) identified by mapping spliced expressed sequence tags (ESTs) and predicted genes (PREDs), identified only in silico. This study underscored the limitations of in silico-only gene prediction, as many PREDs were incorrectly predicted. To refine the HC21 annotation, we have developed a reliable algorithm to extract and stringently map sequences that contain bona fide 3' transcript ends to the genome. We then created a specific 21q graphical display allowing an integrated view of the data that incorporates new ESTs as well as features such as CpG islands, repeats, and gene predictions. Using these tools we identified 27 new putative genes. To validate these, we sequenced previously cloned cDNAs and carried out RT-PCR, 5'- and 3'-RACE procedures, and comparative mapping. These approaches substantiated 19 new transcripts, thus increasing the HC21 gene count by 9.5%. These transcripts were likely not previously identified because they are small and encode small proteins. We also identified four transcriptional units that are spliced but contain no obvious open reading frame. The HC21 data presented here further emphasize that current gene prediction algorithms miss a substantial number of transcripts that nevertheless can be identified using a combination of experimental approaches and multiple refined algorithms.

Ultrasound-guided suprascapular nerve block, description of a novel supraclavicular approach.

BACKGROUND AND OBJECTIVES: The suprascapular nerve (SSN) block is frequently performed for different shoulder pain conditions and for perioperative and postoperative pain control after shoulder surgery. Blind and image-guided techniques have been described, all of which target the nerve within the supraspinous fossa or at the suprascapular notch. This classic target point is not always ideal when ultrasound (US) is used because it is located deep under the muscles, and hence the nerve is not always visible. Blocking the nerve in the supraclavicular region, where it passes underneath the omohyoid muscle, could be an attractive alternative. METHODS: In the first step, 60 volunteers were scanned with US, both in the supraclavicular and the classic target area. The visibility of the SSN in both regions was compared. In the second step, 20 needles were placed into or immediately next to the SSN in the supraclavicular region of 10 cadavers. The accuracy of needle placement was determined by injection of dye and following dissection. RESULTS: In the supraclavicular region of volunteers, the nerve was identified in 81% of examinations (95% confidence interval [CI], 74%-88%) and located at a median depth of 8 mm (interquartile range, 6-9 mm). Near the suprascapular notch (supraspinous fossa), the nerve was unambiguously identified in 36% of examinations (95% CI, 28%-44%) (P < 0.001) and located at a median depth of 35 mm (interquartile range, 31-38 mm; P < 0.001). In the cadaver investigation, the rate of correct needle placement of the supraclavicular approach was 95% (95% CI, 86%-100%). CONCLUSIONS: Visualization of the SSN with US is better in the supraclavicular region as compared with the supraspinous fossa. The anatomic dissections confirmed that our novel supraclavicular SSN block technique is accurate.