Synchronization of EEG: bivariate and multivariate measures.

Synchronization behavior of electroencephalographic (EEG) signals is important for decoding information processing in the human brain. Modern multichannel EEG allows a transition from traditional measurements of synchronization in pairs of EEG signals to whole-brain synchronization maps. The latter can be based on bivariate measures (BM) via averaging over pair-wise values or, alternatively, on multivariate measures (MM), which directly ascribe a single value to the synchronization in a group. In order to compare BM versus MM, we applied nine different estimators to simulated multivariate time series with known parameters and to real EEGs.We found widespread correlations between BM and MM, which were almost frequency-independent for all the measures except coherence. The analysis of the behavior of synchronization measures in simulated settings with variable coupling strength, connection probability, and parameter mismatch showed that some of them, including S-estimator, S-Renyi, omega, and coherence, aremore sensitive to linear interdependences,while others, like mutual information and phase locking value, are more responsive to nonlinear effects. Onemust consider these properties together with the fact thatMM are computationally less expensive and, therefore, more efficient for the large-scale data sets than BM while choosing a synchronization measure for EEG analysis.

Assessment of left coronary artery beat-to-beat repositioning in order to propose an optimal acquisition window for coronary MRA

Introduction: Coronary magnetic resonance angiography (MRA) is a medical imaging technique that involves collecting data from consecutive heartbeats, always at the same time in the cardiac cycle, in order to minimize heart motion artifacts. This technique relies on the assumption that coronary arteries always follow the same trajectory from heartbeat to heartbeat. Until now, choosing the acquisition window in the cardiac cycle was based exclusively on the position of minimal coronary motion. The goal of this study was to test the hypothesis that there are time intervals during the cardiac cycle when coronary beat-to-beat repositioning is optimal. The repositioning uncertainty values in these time intervals were then compared with the intervals of low coronary motion in order to propose an optimal acquisition window for coronary MRA. Methods: Cine breath-hold x-ray angiograms with synchronous ECG were collected from 11 patients who underwent elective routine diagnostic coronarography. Twenty-three bifurcations of the left coronary artery were selected as markers to evaluate repositioning uncertainty and velocity during cardiac cycle. Each bifurcation was tracked by two observers, with the help of a user-assisted algorithm implemented in Matlab (The Mathworks, Natick, MA, USA) that compared the trajectories of the markers coming from consecutive heartbeats and computed the coronary repositioning uncertainty with steps of 50ms until 650ms after the R-wave. Repositioning uncertainty was defined as the diameter of the smallest circle encompassing the points to be compared at the same time after the R-wave. Student's t-tests with a false discovery rate (FDR, q=0.1) correction for multiple comparison were applied to see whether coronary repositioning and velocity vary statistically during cardiac cycle. Bland-Altman plots and linear regression were used to assess intra- and inter-observer agreement. Results: The analysis of left coronary artery beat-to-beat repositioning uncertainty shows a tendency to have better repositioning in mid systole (less than 0.84±0.58mm) and mid diastole (less than 0.89±0.6mm) than in the rest of the cardiac cycle (highest value at 50ms=1.35±0.64mm). According to Student's t-tests with FDR correction for multiple comparison (q=0.1), two intervals, in mid systole (150-200ms) and mid diastole (550-600ms), provide statistically better repositioning in comparison with the early systole and the early diastole. Coronary velocity analysis reveals that left coronary artery moves more slowly in end systole (14.35±11.35mm/s at 225ms) and mid diastole (11.78±11.62mm/s at 625ms) than in the rest of the cardiac cycle (highest value at 25ms: 55.96±22.34mm/s). This was confirmed by Student's t-tests with FDR correction for multiple comparison (q=0.1, FDR-corrected p-value=0.054): coronary velocity values at 225, 575 and 625ms are not much different between them but they are statistically inferior to all others. Bland-Altman plots and linear regression show that intra-observer agreement (y=0.97x+0.02 with R²=0.93 at 150ms) is better than inter-observer (y=0.8x+0.11 with R²=0.67 at 150ms). Discussion: The present study has demonstrated that there are two time intervals in the cardiac cycle, one in mid systole and one in mid diastole, where left coronary artery repositioning uncertainty reaches points of local minima. It has also been calculated that the velocity is the lowest in end systole and mid diastole. Since systole is less influenced by heart rate variability than diastole, it was finally proposed to test an acquisition window between 150 and 200ms after the R-wave.

Interhemispheric coupling between the posterior sylvian regions impacts successful auditory temporal order judgment

Introduction: Accurate registration of the relative timing between the occurrence of sensory events on a sub-second time scale is crucial for both sensory-motor and cognitive functions (Mauk and Buonomano, 2004; Habib, 2000). Support for this assumption comes notably from evidence that temporal processing impairments are implicated in a range of neurological and psychiatric conditions (e.g. Buhusi & Meck, 2005). For instance, deficits in fast auditory temporal integration have been regularly put forward as resulting in phonologic discrimination impairments at the basis of speech comprehension deficits characterizing e.g. dyslexia (Habib, 2000). At least two aspects of the brain mechanisms of temporal order judgment remain unknown. First, it is unknown when during the course of stimulus processing a temporal ,,stamp‟ is established to guide TOJ perception. Second, the extent of interplay between the cerebral hemispheres in engendering accurate TOJ performance is unresolved Methods: We investigated the spatiotemporal brain dynamics of auditory temporal order judgment (aTOJ) using electrical neuroimaging analyses of auditory evoked potentials (AEPs) recorded while participants completed a near-threshold task requiring spatial discrimination of left-right and right-left sound sequences. Results: AEPs to sound pairs modulated topographically as a function of aTOJ accuracy over the 39-77ms post-stimulus period, indicating the engagement of distinct configurations of brain networks during early auditory processing stages. Source estimations revealed that accurate and inaccurate performance were linked to bilateral posterior sylvian regions activity (PSR). However, activity within left, but not right, PSR predicted behavioral performance suggesting that left PSR activity during early encoding phases of pairs of auditory spatial stimuli appears critical for the perception of their order of occurrence. Correlation analyses of source estimations further revealed that activity between left and right PSR was significantly correlated in the inaccurate but not accurate condition, indicating that aTOJ accuracy depends on the functional de-coupling between homotopic PSR areas. Conclusions: These results support a model of temporal order processing wherein behaviorally relevant temporal information - i.e. a temporal 'stamp'- is extracted within the early stages of cortical processes within left PSR but critically modulated by inputs from right PSR. We discuss our results with regard to current models of temporal of temporal order processing, namely gating and latency mechanisms.

Extremal properties of certain risk models

Cette thèse s'intéresse à étudier les propriétés extrémales de certains modèles de risque d'intérêt dans diverses applications de l'assurance, de la finance et des statistiques. Cette thèse se développe selon deux axes principaux, à savoir: Dans la première partie, nous nous concentrons sur deux modèles de risques univariés, c'est-à- dire, un modèle de risque de déflation et un modèle de risque de réassurance. Nous étudions le développement des queues de distribution sous certaines conditions des risques commun¬s. Les principaux résultats sont ainsi illustrés par des exemples typiques et des simulations numériques. Enfin, les résultats sont appliqués aux domaines des assurances, par exemple, les approximations de Value-at-Risk, d'espérance conditionnelle unilatérale etc. La deuxième partie de cette thèse est consacrée à trois modèles à deux variables: Le premier modèle concerne la censure à deux variables des événements extrême. Pour ce modèle, nous proposons tout d'abord une classe d'estimateurs pour les coefficients de dépendance et la probabilité des queues de distributions. Ces estimateurs sont flexibles en raison d'un paramètre de réglage. Leurs distributions asymptotiques sont obtenues sous certaines condi¬tions lentes bivariées de second ordre. Ensuite, nous donnons quelques exemples et présentons une petite étude de simulations de Monte Carlo, suivie par une application sur un ensemble de données réelles d'assurance. L'objectif de notre deuxième modèle de risque à deux variables est l'étude de coefficients de dépendance des queues de distributions obliques et asymétriques à deux variables. Ces distri¬butions obliques et asymétriques sont largement utiles dans les applications statistiques. Elles sont générées principalement par le mélange moyenne-variance de lois normales et le mélange de lois normales asymétriques d'échelles, qui distinguent la structure de dépendance de queue comme indiqué par nos principaux résultats. Le troisième modèle de risque à deux variables concerne le rapprochement des maxima de séries triangulaires elliptiques obliques. Les résultats théoriques sont fondés sur certaines hypothèses concernant le périmètre aléatoire sous-jacent des queues de distributions. -- This thesis aims to investigate the extremal properties of certain risk models of interest in vari¬ous applications from insurance, finance and statistics. This thesis develops along two principal lines, namely: In the first part, we focus on two univariate risk models, i.e., deflated risk and reinsurance risk models. Therein we investigate their tail expansions under certain tail conditions of the common risks. Our main results are illustrated by some typical examples and numerical simu¬lations as well. Finally, the findings are formulated into some applications in insurance fields, for instance, the approximations of Value-at-Risk, conditional tail expectations etc. The second part of this thesis is devoted to the following three bivariate models: The first model is concerned with bivariate censoring of extreme events. For this model, we first propose a class of estimators for both tail dependence coefficient and tail probability. These estimators are flexible due to a tuning parameter and their asymptotic distributions are obtained under some second order bivariate slowly varying conditions of the model. Then, we give some examples and present a small Monte Carlo simulation study followed by an application on a real-data set from insurance. The objective of our second bivariate risk model is the investigation of tail dependence coefficient of bivariate skew slash distributions. Such skew slash distributions are extensively useful in statistical applications and they are generated mainly by normal mean-variance mixture and scaled skew-normal mixture, which distinguish the tail dependence structure as shown by our principle results. The third bivariate risk model is concerned with the approximation of the component-wise maxima of skew elliptical triangular arrays. The theoretical results are based on certain tail assumptions on the underlying random radius.

Addition of inspiratory resistance increases the amplitude of the slow component of O2 uptake kinetics.

The contribution of respiratory muscle work to the development of the O(2) consumption (Vo(2)) slow component is a point of controversy because it has been shown that the increased ventilation in hypoxia is not associated with a concomitant increase in Vo(2) slow component. The first purpose of this study was thus to test the hypothesis of a direct relationship between respiratory muscle work and Vo(2) slow component by manipulating inspiratory resistance. Because the conditions for a Vo(2) slow component specific to respiratory muscle can be reached during intense exercise, the second purpose was to determine whether respiratory muscles behave like limb muscles during heavy exercise. Ten trained subjects performed two 8-min constant-load heavy cycling exercises with and without a threshold valve in random order. Vo(2) was measured breath by breath by using a fast gas exchange analyzer, and the Vo(2) response was modeled after removal of the cardiodynamic phase by using two monoexponential functions. As anticipated, when total work was slightly increased with loaded inspiratory resistance, slight increases in base Vo(2), the primary phase amplitude, and peak Vo(2) were noted (14.2%, P < 0.01; 3.5%, P > 0.05; and 8.3%, P < 0.01, respectively). The bootstrap method revealed small coefficients of variation for the model parameter, including the slow-component amplitude and delay (15 and 19%, respectively), indicating an accurate determination for this critical parameter. The amplitude of the Vo(2) slow component displayed a 27% increase from 8.1 +/- 3.6 to 10.3 +/- 3.4 ml. min(-1). kg(-1) (P < 0.01) with the addition of inspiratory resistance. Taken together, this increase and the lack of any differences in minute volume and ventilatory parameters between the two experimental conditions suggest the occurrence of a Vo(2) slow component specific to the respiratory muscles in loaded condition.

Total hemoglobin mass--a new parameter to detect blood doping?

PURPOSE: All kinds of blood manipulations aim to increase the total hemoglobin mass (tHb-mass). To establish tHb-mass as an effective screening parameter for detecting blood doping, the knowledge of its normal variation over time is necessary. The aim of the present study, therefore, was to determine the intraindividual variance of tHb-mass in elite athletes during a training year emphasizing off, training, and race seasons at sea level. METHODS: tHb-mass and hemoglobin concentration ([Hb]) were determined in 24 endurance athletes five times during a year and were compared with a control group (n = 6). An analysis of covariance was used to test the effects of training phases, age, gender, competition level, body mass, and training volume. Three error models, based on 1) a total percentage error of measurement, 2) the combination of a typical percentage error (TE) of analytical origin with an absolute SD of biological origin, and 3) between-subject and within-subject variance components as obtained by an analysis of variance, were tested. RESULTS: In addition to the expected influence of performance status, the main results were that the effects of training volume (P = 0.20) and training phases (P = 0.81) on tHb-mass were not significant. We found that within-subject variations mainly have an analytical origin (TE approximately 1.4%) and a very small SD (7.5 g) of biological origin. CONCLUSION: tHb-mass shows very low individual oscillations during a training year (<6%), and these oscillations are below the expected changes in tHb-mass due to Herythropoetin (EPO) application or blood infusion (approximately 10%). The high stability of tHb-mass over a period of 1 year suggests that it should be included in an athlete's biological passport and analyzed by recently developed probabilistic inference techniques that define subject-based reference ranges.

Discovery of catalases in members of the Chlamydiales order.

Catalase is an important virulence factor for survival in macrophages and other phagocytic cells. In Chlamydiaceae, no catalase had been described so far. With the sequencing and annotation of the full genomes of Chlamydia-related bacteria, the presence of different catalase-encoding genes has been documented. However, their distribution in the Chlamydiales order and the functionality of these catalases remain unknown. Phylogeny of chlamydial catalases was inferred using MrBayes, maximum likelihood, and maximum parsimony algorithms, allowing the description of three clade 3 and two clade 2 catalases. Only monofunctional catalases were found (no catalase-peroxidase or Mn-catalase). All presented a conserved catalytic domain and tertiary structure. Enzymatic activity of cloned chlamydial catalases was assessed by measuring hydrogen peroxide degradation. The catalases are enzymatically active with different efficiencies. The catalase of Parachlamydia acanthamoebae is the least efficient of all (its catalytic activity was 2 logs lower than that of Pseudomonas aeruginosa). Based on the phylogenetic analysis, we hypothesize that an ancestral class 2 catalase probably was present in the common ancestor of all current Chlamydiales but was retained only in Criblamydia sequanensis and Neochlamydia hartmannellae. The catalases of class 3, present in Estrella lausannensis and Parachlamydia acanthamoebae, probably were acquired by lateral gene transfer from Rhizobiales, whereas for Waddlia chondrophila they likely originated from Legionellales or Actinomycetales. The acquisition of catalases on several occasions in the Chlamydiales suggests the importance of this enzyme for the bacteria in their host environment.

Urea nitrogen, creatinine, and uric acid levels in postmortem serum, vitreous humor, and pericardial fluid.

Urea nitrogen, creatinine, and uric acid are relatively stable in postmortem serum and may, therefore, be used for diagnostic purposes when chronic kidney disease and end-stage renal failure are investigated as causes of death. Nevertheless, uncertainties remain in defining the best alternative to postmortem serum for the identification and assessment of significantly decreased kidney function. In this study, we investigated urea nitrogen, creatinine, and uric acid levels in postmortem serum, pericardial fluid, and vitreous humor in a series of medico-legal cases (500 autopsies) with various causes of death. No postmortem interval-related differences were observed in any of the investigated fluids for any analyzed parameter, confirming the biochemical stability of all compounds after death. Data analysis failed to reveal statistically significant differences between postmortem serum and pericardial fluid urea nitrogen, creatinine, and uric acid concentrations. Conversely, statistically significant differences were observed in all analyzed biomarkers between postmortem serum and vitreous humor levels, with lower concentrations of all markers measured in vitreous. The results of this study suggest that, in order to estimate as accurately as possible blood analyte concentrations at the time of death, pericardial fluid should be preferred to vitreous humor.

Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.

Rôle des cellules souches épithéliales de la cornée et de leur micro-environnement dans le traitement des destructions de la surface oculaire [Darwin or Lamarck? Understanding the ocular surface and its normal or abnormal differentiation in order to cure ocular surface destruction with corneal opacification].

According to the World Health Organization, 5.1% of blindnesses or visual impairments are related to corneal opacification. Cornea is a transparent tissue placed in front of the color of the eye. Its transparency is mandatory for vision. The ocular surface is a functional unit including the cornea and all the elements involved in maintaining its transparency i.e., the eyelids, the conjunctiva, the lymphoid tissue of the conjunctiva, the limbus, the lacrymal glands and the tear film. The destruction of the ocular surface is a disease caused by : traumatisms, infections, chronic inflammations, cancers, toxics, unknown causes or congenital abnormalities. The treatment of the ocular surface destruction requires a global strategy including all the elements that are involved in its physiology. The microenvironnement of the ocular surface must first be restored, i.e., the lids, the conjunctiva, the limbus and the structures that secrete the different layers of the tear film. In a second step, the transparency of the cornea can be reconstructed. A corneal graft performed in a healthy ocular surface microenvironnement will have a better survival rate. To achieve these goals, a thorough understanding of the renewal of the epitheliums and the role of the epithelial stem cells are mandatory.

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.

A hierarchical approach to model parameter optimization for developmental systems.

In the context of Systems Biology, computer simulations of gene regulatory networks provide a powerful tool to validate hypotheses and to explore possible system behaviors. Nevertheless, modeling a system poses some challenges of its own: especially the step of model calibration is often difficult due to insufficient data. For example when considering developmental systems, mostly qualitative data describing the developmental trajectory is available while common calibration techniques rely on high-resolution quantitative data. Focusing on the calibration of differential equation models for developmental systems, this study investigates different approaches to utilize the available data to overcome these difficulties. More specifically, the fact that developmental processes are hierarchically organized is exploited to increase convergence rates of the calibration process as well as to save computation time. Using a gene regulatory network model for stem cell homeostasis in Arabidopsis thaliana the performance of the different investigated approaches is evaluated, documenting considerable gains provided by the proposed hierarchical approach.

Biogeochemical evolution of a marine shore tailings deposit during bioremediation

Summary : Mining activities produce enormous amounts of waste material known as tailings which are composed of fine to medium size particles. These tailings often contain sulfides, which oxidation can lead to acid and metal contamination of water; therefore they need to be remediated. In this work a tailings bioremediation approach was investigated by an interdisciplinary study including geochemistry, mineralogy and microbiology. The aim of the work was to study the effect of the implementation of wetland above oxidizing tailings on the hydrogeology and the biogeochemical element cycles, and to assess the system evolution over time. To reach these goals, biogeochemical processes occurring in a marine shore tailings deposit were investigated. The studied tailings deposit is located at the Bahìa de Ite, Pacific Ocean, southern Peru, where between 1940 and 1996 the tailings were discharged from the two porphyry copper mines Cuajone and Toquepala. After the end of deposition, a remediation approach was initiated in 1997 with a wetland implementation above the oxidizing tailings. Around 90% of the tailings deposits (total 16 km2) were thus remediated, except the central delta area and some areas close to the shoreline. The multi-stable isotope study showed that the tailings were saturated with fresh water in spite of the marine setting, due to the high hydraulic gradient resulting from the wetland implementation. Submarine groundwater discharge (SGD) was the major source of SO4 2-, C1-, Na+, Fe2+, and Mn2+ input into the tailings at the original shelf-seawater interface. The geochemical study (aquatic geochemistry and X-Ray diffraction (XRD) and sequential extractions from the solid fraction) showed that iron and sulfur oxidation were the main processes in the non-remediated tailings, which showed a top a low-pH oxidation zone with strong accumulation of efflorescent salts at the surface due to capillary upward transport of heavy metals (Fe, Cu, Zn, Mn, Cd, Co, and Ni) in the arid climate. The study showed also that the implementation of the wetland resulted in very low concentrations of heavy metals in solution (mainly under the detection limit) due to the near neutral pH and more reducing conditions (100-150 mV). The heavy metals, which were taken from solution, precipitated as hydroxides and sulfides or were bound to organic matter. The bacterial community composition analysis by Terminal Restriction Fragment Length Polymorphism (T-RFLP) and cloning and sequencing of 16S rRNA genes combined with a detailed statistical analysis revealed a high correlation between the bacterial distribution and the geochemical variables. Acidophilic autotrophic oxidizing bacteria were dominating the oxidizing tailings, whereas neutrophilic and heterotrophic reducing bacteria were driving the biogeochemical processes in the remediated tailings below the wetland. At the subsurface of the remediated tailings, an iron cycling was highlighted with oxidation and reduction processes due to micro-aerophilic niches provided by the plant rhizosphere in this overall reducing environment. The in situ bioremediation experiment showed that the main parameter to take into account for the effectiveness was the water table and chemistry which controls the system. The constructed remediation cells were more efficient and rapid in metal removal when saturation conditions were available. This study showed that the bioremediation by wetland implementation could be an effective and rapid treatment for some sulfidic mine tailings deposits. However, the water saturation of the tailings has to be managed on a long-term basis in order to guarantee stability. Résumé : L'activité minière produit d'énormes quantités de déchets géologiques connus sous le nom de « tailings » composées de particules de taille fine à moyenne. Ces déchets contiennent souvent des sulfures dont l'oxydation conduit à la formation d'effluents acides contaminés en métaux, d'où la nécessité d'effectuer une remédiation des sites de stockage concernés. Le but de ce travail est dans un premier temps d'étudier l'effet de la bio-remédiation d'un dépôt de tailings oxydés sur l'hydrogéologie du système et les cycles biogéochimiques des éléments et en second lieu, d'évaluer l'évolution du processus de remédiation dans le temps. Le site étudié dans ce travail est situé dans la Bahía de Ite, au sud du Pérou, au bord de l'Océan Pacifique. Les déchets miniers en question sont déposés dans un environnement marin. De 1940 à 1996, les déchets de deux mines de porphyre cuprifère - Cuajone et Toquepala - ont été acheminés sur le site via la rivière Locumba. En 1997, une première remédiation a été initiée avec la construction d'une zone humide sur les tailings. Depuis, environ 90% de la surface du dépôt (16 km2) a été traité, les parties restantes étant la zone centrale du delta du Locumba et certaines zones proches de la plage. Malgré la proximité de l'océan, les études isotopiques menées dans le cadre de ce travail ont montré que les tailings étaient saturés en eau douce. Cette saturation est due à la pression hydraulique résultant de la mise en place des zones humides. Un écoulement d'eau souterrain sous-marin a été à détecté à l'interface entre les résidus et l'ancien fond marin. En raison de la géologie locale, il constitue une source d'entrée de SO4 2-, Cl-, Na+, FeZ+, et Mn2+ dans le système. L'analyse de la géochimie aquatique, la Diffraction aux Rayons X (XRD) et l'extraction séquentielle ont montré que l'oxydation du fer et .des sulfures est le principal processus se produisant dans les déchets non remédiés. Ceci a entraîné le développement d'une zone d'oxydation à pH bas induisant une forte accumulation des sels efflorescents, conséquence de la migration capillaire des métaux lourds (Fe, Cu, Zn, Mn, Cd, Co et Ni) de la solution vers la surface dans ce climat aride. Cette étude a montré également que la construction de la zone humide a eu comme résultats une précipitation des métaux dans des phases minérales en raison du pH neutre et des conditions réductrices (100-150mV). Les métaux lourds ont précipité sous la forme d'hydroxydes et de sulfures ou sont adsorbés à la matière organique. L'analyse de la composition de la communauté bactérienne à l'aide la technique T-RFLP (Terminal Restriction Fragment Length Polymorphism) et par le clonage/séquençage des gènes de l'ARNr 16S a été combinée à une statistique détaillée. Cette dernière a révélé une forte corrélation entre la distribution de bactéries spécifiques et la géochimie : Les bactéries autotrophes acidophiles dominent dans les déchets oxydés non remédiés, tandis que des bactéries hétérotrophes neutrophiles ont mené les processus microbiens dans les déchets remédiés sous la zone humide. Sous la surface de la zone humide, nos analyses ont également mis en évidence un cycle du fer par des processus d'oxydoréduction rendus possibles par la présence de niches micro-aérées par la rhizosphère dans cet environnement réducteur. L'expérience de bio-remédiation in situ a montré que les paramètres clés qui contrôlent l'efficacité du traitement sont le niveau de la nappe aquifère et la chimie de l'eau. Les cellules de remédiation se sont montrées plus efficaces et plus rapides lorsque le système a pu être saturé en eau. Finalement, cette étude a montré que la bio-remédiation de déchets miniers par la construction de zones humides est un moyen de traitement efficace, rapide et peu coûteux. Cependant, la saturation en eau du système doit être gérée sur le long terme afin de garantir la stabilité de l'ensemble du système.