Simulations in evolution. II. Relative fitness and the propagation of mutants.

In Neo-Darwinism, variation and natural selection are the two evolutionary mechanisms which propel biological evolution. Our previous article presented a histogram model [1] consisting in populations of individuals whose number changed under the influence of variation and/or fitness, the total population remaining constant. Individuals are classified into bins, and the content of each bin is calculated generation after generation by an Excel spreadsheet. Here, we apply the histogram model to a stable population with fitness F(1)=1.00 in which one or two fitter mutants emerge. In a first scenario, a single mutant emerged in the population whose fitness was greater than 1.00. The simulations ended when the original population was reduced to a single individual. The histogram model was validated by excellent agreement between its predictions and those of a classical continuous function (Eqn. 1) which predicts the number of generations needed for a favorable mutation to spread throughout a population. But in contrast to Eqn. 1, our histogram model is adaptable to more complex scenarios, as demonstrated here. In the second and third scenarios, the original population was present at time zero together with two mutants which differed from the original population by two higher and distinct fitness values. In the fourth scenario, the large original population was present at time zero together with one fitter mutant. After a number of generations, when the mutant offspring had multiplied, a second mutant was introduced whose fitness was even greater. The histogram model also allows Shannon entropy (SE) to be monitored continuously as the information content of the total population decreases or increases. The results of these simulations illustrate, in a graphically didactic manner, the influence of natural selection, operating through relative fitness, in the emergence and dominance of a fitter mutant.

Development of CBCT-based prostate setup correction strategies and impact of rectal distension.

BACKGROUND: Cone-beam computed tomography (CBCT) image-guided radiotherapy (IGRT) systems are widely used tools to verify and correct the target position before each fraction, allowing to maximize treatment accuracy and precision. In this study, we evaluate automatic three-dimensional intensity-based rigid registration (RR) methods for prostate setup correction using CBCT scans and study the impact of rectal distension on registration quality. METHODS: We retrospectively analyzed 115 CBCT scans of 10 prostate patients. CT-to-CBCT registration was performed using (a) global RR, (b) bony RR, or (c) bony RR refined by a local prostate RR using the CT clinical target volume (CTV) expanded with 1-to-20-mm varying margins. After propagation of the manual CT contours, automatic CBCT contours were generated. For evaluation, a radiation oncologist manually delineated the CTV on the CBCT scans. The propagated and manual CBCT contours were compared using the Dice similarity and a measure based on the bidirectional local distance (BLD). We also conducted a blind visual assessment of the quality of the propagated segmentations. Moreover, we automatically quantified rectal distension between the CT and CBCT scans without using the manual CBCT contours and we investigated its correlation with the registration failures. To improve the registration quality, the air in the rectum was replaced with soft tissue using a filter. The results with and without filtering were compared. RESULTS: The statistical analysis of the Dice coefficients and the BLD values resulted in highly significant differences (p<10(-6)) for the 5-mm and 8-mm local RRs vs the global, bony and 1-mm local RRs. The 8-mm local RR provided the best compromise between accuracy and robustness (Dice median of 0.814 and 97% of success with filtering the air in the rectum). We observed that all failures were due to high rectal distension. Moreover, the visual assessment confirmed the superiority of the 8-mm local RR over the bony RR. CONCLUSION: The most successful CT-to-CBCT RR method proved to be the 8-mm local RR. We have shown the correlation between its registration failures and rectal distension. Furthermore, we have provided a simple (easily applicable in routine) and automatic method to quantify rectal distension and to predict registration failure using only the manual CT contours.

The tectonically controlled emplacement of a vertically sheeted gabbro-pyroxenite intrusion: Feeder-zone of an ocean-island volcano (Fuerteventura, Canary Islands)

The Miocene PX1 gabbro-pyroxenite pluton, Fuerteventura, Canary Islands, is a 3.5 x 5.5 km shallow-level intrusion (0.15-0.2 GPa and 1100-1120 degrees C), interpreted as the feeder-zone to an ocean-island volcano. It displays a vertical magmatic banding expressed in five 50 to 100 metre-wide NNE-SSW trending alkaline gabbro sequences alternating with pyroxenites. This emplacement geometry was controlled by brittle to ductile shear zones, generated by a regional E-W extensional tectonic setting that affected Fuerteventura during the Miocene. At a smaller scale, the PX1 gabbro and pyroxenite bands consist of metre-thick differentiation units, which suggest emplacement by periodic injection of magma pulses as vertical dykes that amalgamated, similarly to a sub-volcanic sheeted dyke complex. Individual dykes underwent internal differentiation following a solidification front parallel to the dyke edges. This solidification front may have been favoured by a significant lateral/horizontal thermal gradient, expressed by the vertical banding in the gabbros, the fractionation asymmetry within individual dykes and the migmatisation of the wall rocks. Pyroxenitic layers result from the fractionation and accumulation of clinopyroxene +/- olivine +/- plagioclase crystals from a mildly alkaline basaltic liquid. They are interpreted as truncated differentiation sequences, from which residual melts were extracted at various stages of their chemical evolution by subsequent dyke intrusions, either next to or within the crystallising unit. Compaction and squeezing of the crystal mush is ascribed to the incoming and inflating magma pulses. The expelled interstitial liquid was likely collected and erupted along with the magma flowing through the newly injected dykes. Clinopyroxene mineral orientation - as evidenced by EBSD and micro X-ray tomography investigations - displays a marked pure-shear component, supporting the interpretation of the role of compaction in the generation of the pyroxenites. Conversely, gabbro sequences underwent minor melt extraction and are believed to represent crystallised coalesced magma batches emplaced at lower rates at the end of eruptive cycles. Clinopyroxene orientations in gabbros record a simple shear component suggesting syn-magmatic deformation parallel to observed NNE-SSW trending shear zones induced by the regional tensional stress field. This emplacement model implies a crystallisation time of 1 to 5 years for individual dykes, consistent with PX1 emplacement over less than 0.5 My. A minimum amount of approximately 150 km(3) of magma is needed to generate the pluton, part of it having been erupted through the Central Volcanic Centre of Fuerteventura. If the regional extensional tectonic regime controls the PX1 feeder-zone initiation and overall geometry, rates and volumes of magma depend on other, source-related factors. High injection rates are likely to induce intrusion growth rates larger than could be accommodated by the regional extension. In this case, dyke intrusion by propagation of a weak tip, combined with the inability of magma to circulate through previously emplaced and crystallised dykes could result in an increase of non-lithostatic pressure on previously emplaced mushy dyke walls; thus generating strong pure-shear compaction within the pluton feeder-zone and interstitial melt expulsion. These compaction-dominated processes are recorded by the cumulitic pyroxenite bands. (C) 2010 Elsevier B.V. All rights reserved.

Expression analysis of the AtPHO1 gene family in Arabidopsis thaliana and the involvement of AtPHO1 in the regulation of stomatal movements

Résumé Le transfert du phosphate des racines vers les feuilles s'effectue par la voie du xylème. Il a été précédemment démontré que la protéine AtPHO1 était indispensable au transfert du phosphate dans les vaisseaux du xylème des racines chez la plante modèle Arabidopsis thaliana. Le séquençage et l'annotation du génome d'Arabidopsis ont permis d'identifier dix séquences présentant un niveau de similarité significatif avec le gène AtPHO1 et constituant une nouvelle famille de gène appelé la famille de AtPHO1. Basée sur une étude moléculaire et génétique, cette thèse apporte des éléments de réponse pour déterminer le rôle des membres de ia famille de AtPHO1 chez Arabidopsis, inconnue à ce jour. Dans un premier temps, une analyse bioinformatique des séquences protéiques des membres de la famille de AtPHO1 a révélé la présence dans leur région N-terminale d'un domaine nommé SPX. Ce dernier est conservé parmi de nombreuses protéines impliquées dans l'homéostasie du phosphate chez la levure, renforçant ainsi l'hypothèse que les membres de la famille de AtPHO1 auraient comme AtPHO1 un rôle dans l'équilibre du phosphate dans la plante. En parallèle, la localisation tissulaire de l'expression des gènes AtPHO dans Arabidopsis a été identifiée par l'analyse de plantes transgéniques exprimant le gène rapporteur uidA sous le contrôle des promoteurs respectifs des gènes AtPHO. Un profil d'expression de chaque gène AtPHO au cours du développement de la plante a été obtenu. Une expression prédominante au niveau des tissus vasculaires des racines, des feuilles, des tiges et des fleurs a été observée, suggérant que les gènes AtPHO pourraient avoir des fonctions redondantes au niveau du transfert de phosphate dans le cylindre vasculaire de ces différents organes. Toutefois, plusieurs régions promotrices des gènes AtPHO contrôlent également un profil d'expression GUS non-vasculaire, indiquant un rôle putatif des gènes AtPHO dans l'acquisition ou le recyclage de phosphate dans la plante. Dans un deuxième temps, l'analyse de l'expression des gènes AtPHO durant une carence en phosphate a établi que seule l'expression des gènes AtPHO1, AtPHO1; H1 et AtPHO1; H10 est régulée par cette carence. Une étude approfondie de leur expression en réponse à des traitements affectant l'homéostasie du phosphate dans la plante a ensuite démontré leur régulation par différentes voies de signalisation. Ensuite, une analyse détaillée de la régulation de l'expression du gène AtPHO1; H1O dans des feuilles d'Arabidopsis blessées ou déshydratées a révélé que ce gène constitue le premìer gène marqueur d'une nouvelle voie de signalisation induite par l'OPDA, pas par le JA et dépendante de la protéine COI1. Ces résultats démontrent pour la première fois que l'OPDA et le JA peuvent activer différents gènes via des voies de signalisation dépendantes de COI1. Enfin, cette thèse révèle l'identification d'un nouveau rôle de la protéine AtPHO1 dans la régulation de l'action de l'ABA au cours des processus de fermeture stomatique et de germination des graines chez Arabidopsis. Bien que les fonctions exactes des protéines AtPHO restent à être déterminées, ce travail de thèse suggère leur implication dans la propagation de différents signaux dans la plante via la modulation du potentiel membranaire et/ou l'affectation de la composition en ions des cellules comme le font de nombreux transporteurs ou régulateur du transport d'ions. Summary Phosphate is transferred from the roots to the shoot via the xylem. The requirement for AtPHO1 protein to transfer phosphate to the xylem vessels of the root has been previously demonstrated in Arabidopsis thaliana. The sequencing and the annotation of the Arabidopsis genome had allowed the identification of ten sequences that show a significant level of similarity with the AtPHO1 gene. These 10 genes, of unknown functions, constitute a new gene family called the AtPHO1 gene family. Based on a molecular and genetics study, this thesis reveals some information needed to understand the role of the AtPHO1 family members in the plant Arabidopsis. First, a bioinformatics study revealed that the AtPHO sequences contained, in the N-terminal hydrophilic region, a motif called SPX and conserved among multiple proteins involved in phosphate homeostasis in yeast. This finding reinforces the hypothesis that all AtPHO1 family members have, as AtPHO1, a role in phosphate homeostasis. In parallel, we identified the pattern of expression of AtPHO genes in Arabidopsis via analysis of transgenic plants expressing the uidA reporter gene under the control of respective AtPHO promoter regions. The results exhibit a predominant expression of AtPHO genes in vascular tissues of all organs of the plant, implying that these AtPHO genes could have redundant functions in the transfer of phosphate to the vascular cylinder of various organs. The GUS expression pattern for several AtPHO promoter regions was also detected in non-vascular tissue indicating a broad role of AtPHO genes in the acquisition or in the recycling of phosphate in the plant. In a second step, the analysis of the expression of AtPHO genes during phosphate starvation established that only the expression of the AtPHO1, AtPHO1; H1 and AtPHO1; H10 genes were regulated by Pi starvation. Interestingly, different signalling pathways appeared to regulate these three genes during various treatments affecting Pi homeostasis in the plant. The third chapter presents a detailed analysis of the signalling pathways regulating the expression of the AtPHO1; H10 gene in Arabidopsis leaves during wound and dehydrated stresses. Surprisingly, the expression of AtPHO1; H10 was found to be regulated by OPDA (the precursor of JA) but not by JA itself and via the COI1 protein (the central regulator of the JA signalling pathway). These results demonstrated for the first time that OPDA and JA could activate distinct genes via COI1-dependent pathways. Finally, this thesis presents the identification of a novel role of the AtPHO1 protein in the regulation of ABA action in Arabidopsis guard cells and during seed germination. Although the exact role and function of AtPHO1 still need to be determined, these last findings suggest that AtPHO1 and by extension other AtPHO proteins could mediate the propagation of various signals in the plant by modulating the membrane potential and/or by affecting cellular ion composition, as it is the case for many ion transporters or regulators of ion transport.

Comparison of porcine and human coagulation by thrombelastometry.

Introduction: Although the pig is a standard model for the evaluation of various diseases in humans, including coagulopathy, it is not clear whether results in animals can be extrapolated to man.Materials and methods: In 75 anesthetized pigs, we assessed reagent-supported thrombelastometry (ExTEM (R)), platelet-blocked thrombelastometry (FibTEM (R)), and aprotinin thrombelastometry (ApTEM (R)). Results were compared to values from 13 anesthetized humans.Results (median, 95% CI): ExTEM (R) : While clot strength was comparable in pigs (66 mm, 65-67 mm) and in humans (64 mm, 60-68 mm; NS), clotting time in animals was longer (pigs 64 s, 62-66 s; humans 55 s, 49-71 s; P<0.05) and clot formation time shorter (pigs 52 s, 49-54 s; humans 83 s, 67-98 s, P<0.001). The clot lysis index at 30 minutes was lower in animals (96.9%, 95.1-97.3%) than in humans (99.5%, 98.6-99.9%; P<0.001). ApTEM (R) showed no hyperfibrinolysis in animals. Modification of the anesthesia protocol in animals resulted in significant ExTEM (R) changes. FibTEM (R) : Complete platelet inhibition yielded significantly higher platelet contribution to clot strength in pigs (79%, 76-81%) than in humans (73%, 71-77%; P<0.05), whereas fibrinogen contribution to clot strength was higher in humans (27%, 24-29%) than in animals (21%, 19-24%; P<0.05).Conclusions: Maximum clot firmness is comparable in human and porcine blood. However, clot lysis, platelet and fibrinogen contribution to clot strength, as well as initiation and propagation of clotting, are considerably different between pigs and humans. In addition, anesthesic drugs seem to influence thrombelastometry in animals. Accordingly, coagulation abnormalities in pigs subjected to diseases may not necessarily represent the coagulation profile in sick patients. (C) 2011 Elsevier Ltd. All rights reserved.

An activation-independent role of transcription factors in insulator function.

Chromatin insulators are defined as transcriptionally neutral elements that prevent negative or positive influence from extending across chromatin to a promoter. Here we show that yeast subtelomeric anti-silencing regions behave as boundaries to telomere-driven silencing and also allow discontinuous propagation of silent chromatin. These two facets of insulator activity, boundary and silencing discontinuity, can be recapitulated by tethering various transcription activation domains to tandem sites on DNA. Importantly, we show that these insulator activities do not involve direct transcriptional activation of the reporter promoter. These findings predict that certain promoters behave as insulators and partition genomes in functionally independent domains.

Technical note: Implementation of a geodatabase of published and unpublished data on the catastrophic Vaiont landslide

On 9 October 1963 a catastrophic landslide suddenly occurred on the southern slope of the Vaiont dam reservoir. A mass of approximately 270 million m3 collapsed into the reservoir generating a wave that overtopped the dam and hit the town of Longarone and other villages nearby. Several investigations and interpretations of the slope collapse have been carried out during the last 45 years, however, a comprehensive explanation of both the triggering and the dynamics of the phenomenon has yet to be provided. In order to re-evaluate the currently existing information on the slide, an electronic bibliographic database and an ESRI-geodatabase have been developed. The chronology of the collected documentation showed that most of the studies for re-evaluating the failure mechanisms were conducted in the last decade, as a consequence of knowledge, methods and techniques recently acquired. The current contents of the geodatabase will improve definition of the structural setting that influenced the slide and led to the the propagation of the displaced rock mass. The objectives, structure and contents of the e-bibliography and Geodatabase are indicated, together with a brief description on the possible use of the alphanumeric and spatial contents of the databases.

Brain network analysis of patients with Temporal Lobe Epilepsy

Patients with Temporal Lobe Epilepsy (TLE) suffer from widespread subtle white matter abnormalities and abnormal functional connectivity extending beyond the affected lobe, as revealed by Diffusion Tensor MR Imaging, volumetric and functional MRI studies. Diffusion Spectrum Imaging (DSI) is a diffusion imaging technique with high angular resolution for improving the mapping of white matter pathways. In this study, we used DSI, connectivity matrices and topological measures to investigate how the alteration in structural connectivity influences whole brain structural networks. Eleven patients with right-sided TLE and hippocampal sclerosis and 18 controls underwent our DSI protocol at 3T. The cortical and subcortical grey matters were parcellated into 86 regions of interest and the connectivity between every region pair was estimated using global tractography and a connectivity matrix (the adjacency matrix of the structural network). We then compared the networks of patients and controls using topological measures. In patients, we found a higher characteristic path length and a lower clustering coefficient compared to controls. Local measures at node level of the clustering and efficiency showed a significant difference after a multiple comparison correction (Bonferroni). These significant nodes were located within as well outside the temporal lobe, and the localisation of most of them was consistent with regions known to be part of epileptic networks in TLE. Our results show altered connectivity patterns that are concordant with the mapping of functional epileptic networks in patients with TLE. Further studies are needed to establish the relevance of these findings for the propagation of epileptic activity, cognitive deficits in medial TLE and outcome of epilepsy surgery in individual patients.

Analysis of a boundary protein delimiting chromatin domains in yeast

SUMMARY The expression state of a eukaryotic gene depends in part on its location in the chromosome. This position effect results from the organization of eukaryotic genomes into discrete functional domains, defined by local differences in chromatin structure. The expression of genes within each domain appears to be defined and maintained by the concerted action of regulatory elements such as promoters, enhancers, silencers and locus control regions. Individual domains may be bordered by boundary elements that separate regions of permissive and silent chromatin. When located next to chromosomal elements such as telomeres, genes can be subjected to epigenetic silencing. In yeast, this is mediated by the propagation of the SIR proteins from telomeres towards more centromeric regions. Particular transcription factors can protect downstream genes from silencing when tethered between the gene and the telomere, and they may thus act as chromatin domain boundaries. Here we have studied one of these transcription factors, CTF-1, that binds directly histone H3. A deletion mutagenesis localized the barrier activity to CTF-1 histone-binding domain. A saturating point mutagenesis of this domain identified several amino-acid substitutions that similarly inhibited the boundary and histone-binding activities. Chromatin immunoprecipitation experiments indicated that the barrier protein efficiently prevents the spreading of SIR proteins, and that it separates domains of hypoacetylated and hyperacetylated histones. Together, these results suggest a mechanism by which proteins such as CTF-1 may interact directly with histone H3 to prevent the propagation of a silent chromatin structure, thereby defining boundaries of permissive and silent chromatin domains. RESUME L'expression des gènes eucaryotes dépend en partie de leur localisation sur les chromosomes. Cet effet de position résulte de l'organisation des génomes eucaryotes en domaines fonctionnels, définis par des changements locaux au niveau de la structure de la chromatine. Dans chacun de ces domaines, l'expression des gènes est définie et maintenue par l'action concertée de différents éléments régulateurs tels que les promoteurs, les amplificateurs, les silenceurs et les locus control régions. Ces domaines peuvent être entourés par des éléments barrière, séparant les régions de chromatine répressive des régions permissive pour l'expression des gènes. Lorsqu'ils se situent à proximité d'éléments chromosomiques comme les telomères, les gènes peuvent être réprimés de manière épigénétique. Chez la levure, cette répression est établie par la propagation des protéines SIR depuis les télomères vers les régions centromériques. Certains facteurs de transcription peuvent empêcher la répression d'un gène, lorsqu'ils sont placés entre ce gène et le télomère. Nous avons étudié un de ces facteurs, CTF-1, qui a la particularité de lier directement l'histone H3. La délétion de certaines parties de CTF-1 a permis de déterminer que la région responsable de l'activité barrière correspond au domaine d'interaction avec H3. Plusieurs mutations points effectuées dans ce domaine inhibent à la fois l'activité barrière et la capacité de lier H3. Des expériences d'immuno-précipitation de la chromatine indiquent que la protéine barrière CTF-1 prévient efficacement la propagation des protéines SIR et sépare des domaines contenant des histones hypo-acétylées de ceux constitués d'histones hyper-acétylées. Ces résultats suggèrent que CTF-1 interagit directement avec l'histone H3 pour empêcher la propagation de la chromatine répressive, délimitant ainsi des domaines de chromatine permissive et des domaines de chromatine silencieuse.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Characterization of insulators and barrier elements for gene therapy

Gene transfer that relies on integrating vectors often suffers from epigenetic or regulatory effects that influence the expression of the therapeutic gene and=or of cellular genes located near the vector integration site in the chromosome. Insulator elements act to block gene activation by enhancers, while chromatin domain boundary or barrier sequences prevent gene-silencing effects. At present, the modes of action of insulator and barriers are poorly understood, and their use in the context of gene therapies remains to be documented. Using combinations of reporter genes coding for indicator fluorescent proteins, we constructed assay systems that allow the quantification of the insulator or barrier activities of genetic elements in individual cells. This presentation will illustrate how these assay systems were used to identify short DNA elements that insulate nearby genes from activation by viral vector elements, and=or that block the propagation of a silent chromatin structure that leads to gene silencing. We will show that some barrier elements do not merely block repressive effects, but that they can act to stabilize and sustain transgene expression. We will illustrate that this may be beneficial when transgenes are introduced into stem or precursor cells using non-viral vectors, where later differentiation may lead to the silencing of the therapeutic gene. We will show that these elements can be used to maintain efficient transgene expression upon the differentiation of murine precursor cells towards myofibers, in a model of cell therapy for muscle dystrophies.

Rate and processes of river network rearrangement during incipient faulting: The case of the Cahabon river, Guatemala

Deeply incised river networks are generally regarded as robust features that are not easily modified by erosion or tectonics. Although the reorganization of deeply incised drainage systems has been documented, the corresponding importance with regard to the overall landscape evolution of mountain ranges and the factors that permit such reorganizations are poorly understood. To address this problem, we have explored the rapid drainage reorganization that affected the Cahabon River in Guatemala during the Quaternary. Sediment-provenance analysis, field mapping, and electrical resistivity tomography (ERT) imaging are used to reconstruct the geometry of the valley before the river was captured. Dating of the abandoned valley sediments by the Be-10-Al-26 burial method and geomagnetic polarity analysis allow us to determine the age of the capture events and then to quantify several processes, such as the rate of tectonic deformation of the paleovalley, the rate of propagation of post-capture drainage reversal, and the rate at which canyons that formed at the capture sites have propagated along the paleovalley. Transtensional faulting started 1 to 3 million years ago, produced ground tilting and ground faulting along the Cahabon River, and thus generated differential uplift rate of 0.3 +/- 0.1 up to 0.7 +/- 0.4 mm . y(-1) along the river's course. The river responded to faulting by incising the areas of relative uplift and depositing a few tens of meters of sediment above the areas of relative subsidence. Then, the river experienced two captures and one avulsion between 700 ky and 100 ky. The captures breached high-standing ridges that separate the Cahabon River from its captors. Captures occurred at specific points where ridges are made permeable by fault damage zones and/or soluble rocks. Groundwater flow from the Cahabon River down to its captors likely increased the erosive power of the captors thus promoting focused erosion of the ridges. Valley-fill formation and capture occurred in close temporal succession, suggesting a genetic link between the two. We suggest that the aquifers accumulated within the valley-fills, increased the head along the subterraneous system connecting the Cahabon River to its captors, and promoted their development. Upon capture, the breached valley experienced widespread drainage reversal toward the capture sites. We attribute the generalized reversal to combined effects of groundwater sapping in the valley-fill, axial drainage obstruction by lateral fans, and tectonic tilting. Drainage reversal increased the size of the captured areas by a factor of 4 to 6. At the capture sites, 500 m deep canyons have been incised into the bedrock and are propagating upstream at a rate of 3 to 11 mm . y(-1) deepening at a rate of 0.7 to 1 5 mm . y(-1). At this rate, 1 to 2 million years will be necessary for headward erosion to completely erase the topographic expression of the paleovalley. It is concluded that the rapid reorganization of this drainage system was made possible by the way the river adjusted to the new tectonic strain field, which involved transient sedimentation along the river's course. If the river had escaped its early reorganization and had been given the time necessary to reach a new dynamic equilibrium, then the transient conditions that promoted capture would have vanished and its vulnerability to capture would have been strongly reduced.

Regulation of Nav1.7 and Nav1.8 voltage-gated sodium channels by Nedd4-2 and β-subunits and its implication in neuropathic pain

In mammals, the presence of excitable cells in muscles, heart and nervous system is crucial and allows fast conduction of numerous biological information over long distances through the generation of action potentials (AP). Voltage-gated sodium channels (Navs) are key players in the generation and propagation of AP as they are responsible for the rising phase of the AP. Navs are heteromeric proteins composed of a large pore-forming a-subunit (Nav) and smaller ß-auxiliary subunits. There are ten genes encoding for Navl.l to Nav1.9 and NaX channels, each possessing its own specific biophysical properties. The excitable cells express differential combinations of Navs isoforms, generating a distinct electrophysiological signature. Noteworthy, only when anchored at the membrane are Navs functional and are participating in sodium conductance. In addition to the intrinsic properties of Navs, numerous regulatory proteins influence the sodium current. Some proteins will enhance stabilization of membrane Navs while others will favour internalization. Maintaining equilibrium between the two is of crucial importance for controlling cellular excitability. The E3 ubiquitin ligase Nedd4-2 is a well-characterized enzyme that negatively regulates the turnover of many membrane proteins including Navs. On the other hand, ß-subunits are known since long to stabilize Navs membrane anchoring. Peripheral neuropathic pain is a disabling condition resulting from nerve injury. It is characterized by the dysregulation of Navs expressed in dorsal root ganglion (DRG) sensory neurons as highlighted in different animal models of neuropathic pain. Among Navs, Nav1.7 and Nav1.8 are abundantly and specifically expressed in DRG sensory neurons and have been recurrently incriminated in nociception and neuropathic pain development. Using the spared nerve injury (SNI) experimental model of neuropathic pain in mice, I observed a specific reduction of Nedd4-2 in DRG sensory neurons. This decrease subsequently led to an upregulation of Nav1.7 and Nav1.8 protein and current, in the axon and the DRG neurons, respectively, and was sufficient to generate neuropathic pain-associated hyperexcitability. Knocking out Nedd4-2 specifically in nociceptive neurons led to the same increase of Nav1.7 and Nav1.8 concomitantly with an increased thermal sensitivity in mice. Conversely, rescuing Nedd4-2 downregulation using viral vector transfer attenuated neuropathic pain mechanical hypersensitivity. This study demonstrates the significant role of Nedd4-2 in regulating cellular excitability in vivo and its involvement in neuropathic pain development. The role of ß-subunits in neuropathic pain was already demonstrated in our research group. Because of their stabilization role, the increase of ßl, ß2 and ß3 subunits in DRGs after SNI led to increased Navs anchored at the membrane. Here, I report a novel mechanism of regulation of a-subunits by ß- subunits in vitro; ßl and ß3-subunits modulate the glycosylation pattern of Nav1.7, which might account for stabilization of its membrane expression. This opens new perspectives for investigation Navs state of glycosylation in ß-subunits dependent diseases, such as in neuropathic pain. - Chez les mammifères, la présence de cellules excitables dans les muscles, le coeur et le système nerveux est cruciale; elle permet la conduction rapide de nombreuses informations sur de longues distances grâce à la génération de potentiels d'action (PA). Les canaux sodiques voltage-dépendants (Navs) sont des participants importants dans la génération et la propagation des PA car ils sont responsables de la phase initiale de dépolarisation du PA. Les Navs sont des protéines hétéromériques composées d'une grande sous-unité a (formant le pore du canal) et de petites sous-unités ß accompagnatrices. Il existe dix gènes qui codent pour les canaux sodiques, du Nav 1.1 au Nav 1.9 ainsi que NaX, chacun possédant des propriétés biophysiques spécifiques. Les cellules excitables expriment différentes combinaisons des différents isoformes de Navs, qui engendrent une signature électrophysiologique distincte. Les Navs ne sont fonctionnels et ne participent à la conductibilité du Na+, que s'ils sont ancrés à la membrane plasmique. En plus des propriétés intrinsèques des Navs, de nombreuses protéines régulatrices influencent également le courant sodique. Certaines protéines vont favoriser l'ancrage et la stabilisation des Navs exprimés à la membrane, alors que d'autres vont plutôt favoriser leur internalisation. Maintenir l'équilibre des deux processus est crucial pour contrôler l'excitabilité cellulaire. Dans ce contexte, Nedd4-2, de la famille des E3 ubiquitin ligase, est une enzyme bien caractérisée qui régule l'internalisation de nombreuses protéines, notamment celle des Navs. Inversement, les sous-unités ß sont connues depuis longtemps pour stabiliser l'ancrage des Navs à la membrane. La douleur neuropathique périphérique est une condition débilitante résultant d'une atteinte à un nerf. Elle est caractérisée par la dérégulation des Navs exprimés dans les neurones sensoriels du ganglion spinal (DRG). Ceci a été démontré à de multiples occasions dans divers modèles animaux de douleur neuropathique. Parmi les Navs, Nav1.7 et Nav1.8 sont abondamment et spécifiquement exprimés dans les neurones sensoriels des DRG et ont été impliqués de façon récurrente dans le développement de la douleur neuropathique. En utilisant le modèle animal de douleur neuropathique d'épargne du nerf sural (spared nerve injury, SNI) chez la souris, j'ai observé une réduction spécifique des Nedd4-2 dans les neurones sensoriels du DRG. Cette diminution avait pour conséquence l'augmentation de l'expression des protéines et des courants de Nav 1.7 et Nav 1.8, respectivement dans l'axone et les neurones du DRG, et était donc suffisante pour créer l'hyperexcitabilité associée à la douleur neuropathique. L'invalidation pour le gène codant pour Nedd4-2 dans une lignée de souris génétiquement modifiées a conduit à de similaires augmentations de Nav1.7 et Nav1.8, parallèlement à une augmentation à la sensibilité thermique. A l'opposé, rétablir une expression normale de Nedd4-2 en utilisant un vecteur viral a eu pour effet de contrecarrer le développement de l'hypersensibilité mécanique lié à ce modèle de douleur neuropathique. Cette étude démontre le rôle important de Nedd4-2 dans la régulation de l'excitabilité cellulaire in vivo et son implication dans le développement des douleurs neuropathiques. Le rôle des sous-unités ß dans les douleurs neuropathiques a déjà été démontré dans notre groupe de recherche. A cause de leur rôle stabilisateur, l'augmentation des sous-unités ßl, ß2 et ß3 dans les DRG après SNI, conduit à une augmentation des Navs ancrés à la membrane. Dans mon travail de thèse, j'ai observé un nouveau mécanisme de régulation des sous-unités a par les sous-unités ß in vitro. Les sous-unités ßl et ß3 régulent l'état de glycosylation du canal Nav1.7, et stabilisent son expression membranaire. Ceci ouvre de nouvelles perspectives dans l'investigation de l'état de glycosylation des Navs dans des maladies impliquant les sous-unités ß, notamment les douleurs neuropathiques.

Influence of bedrock structures on the spatial pattern of erosional landforms in small alpine catchments

Structural settings and lithological characteristics are traditionally assumed to influence the development of erosional landforms, such as gully networks and rock couloirs, in steep mountain rock basins. The structural control of erosion of two small alpine catchments of distinctive rock types is evaluated by comparing the correspondences between the orientations of their gullies and rock couloirs with (1) the sliding orientations of potential slope failures mechanisms, and (2) the orientation of the maximum joint frequency, this latter being considered as the direction exploited primarily by erosion and mass wasting processes. These characteristic orientations can be interpreted as structural weaknesses contributing to the initiation and propagation of erosion. The morphostructural analysis was performed using digital elevation models and field observations. The catchment comprised of magmatic intrusive rocks shows a clear structural control, mostly expressed through potential wedges failure. Such joint configurations have a particular geometry that encourages the development of gullies in hard rock, e.g. through enhanced gravitational and hydrological erosional processes. In the catchment underlain by sedimentary rocks, penetrative joints that act as structural weaknesses seem to be exploited by gullies and rock couloirs. However, the lithological setting and bedding configuration prominently control the development of erosional landforms, and influence not only the local pattern of geomorphic features, but the general morphology of the catchment. The orientations of the maximum joint frequency are clearly associated with the gully network, suggesting that its development is governed by anisotropy in rock strength. These two catchments are typical of bedrock-dominated basins prone to intense processes of debris supply. This study suggests a quantitative approach for describing the relationship between bedrock jointing and geomorphic features geometry. Incorporation of bedrock structure can be relevant when studying processes governing the transfer of clastic material, for the assessment of sediment yields and in landforms evolution models.

Role of glutamate in neuron-glia metabolic coupling.

The coupling between synaptic activity and glucose utilization (neurometabolic coupling) is a central physiologic principle of brain function that has provided the basis for 2-deoxyglucose-based functional imaging with positron emission tomography. Approximately 10 y ago we provided experimental evidence that indicated a central role of glutamate signaling on astrocytes in neurometabolic coupling. The basic mechanism in neurometabolic coupling is the glutamate-stimulated aerobic glycolysis in astrocytes, such that the sodium-coupled reuptake of glutamate by astrocytes and the ensuing activation of the Na(+)-K(+) ATPase triggers glucose uptake and its glycolytic processing, which results in the release of lactate from astrocytes. Lactate can then contribute to the activity-dependent fueling of the neuronal energy demands associated with synaptic transmission. Analyses of this coupling have been extended in vivo and have defined the methods of coupling for inhibitory neurotransmission as well as its spatial extent in relation to the propagation of metabolic signals within the astrocytic syncytium. On the basis of a large body of experimental evidence, we proposed an operational model, "the astrocyte-neuron lactate shuttle." A series of results obtained by independent laboratories have provided further support for this model. This body of evidence provides a molecular and cellular basis for interpreting data that are obtained with functional brain imaging studies.