New perspectives in the use of ink evidence in forensic science : part III : operational applications and evaluation

The research reported in this series of article aimed at (1) automating the search of questioned ink specimens in ink reference collections and (2) at evaluating the strength of ink evidence in a transparent and balanced manner. These aims require that ink samples are analysed in an accurate and reproducible way and that they are compared in an objective and automated way. This latter requirement is due to the large number of comparisons that are necessary in both scenarios. A research programme was designed to (a) develop a standard methodology for analysing ink samples in a reproducible way, (b) comparing automatically and objectively ink samples and (c) evaluate the proposed methodology in forensic contexts. This report focuses on the last of the three stages of the research programme. The calibration and acquisition process and the mathematical comparison algorithms were described in previous papers [C. Neumann, P. Margot, New perspectives in the use of ink evidence in forensic science-Part I: Development of a quality assurance process for forensic ink analysis by HPTLC, Forensic Sci. Int. 185 (2009) 29-37; C. Neumann, P. Margot, New perspectives in the use of ink evidence in forensic science-Part II: Development and testing of mathematical algorithms for the automatic comparison of ink samples analysed by HPTLC, Forensic Sci. Int. 185 (2009) 38-50]. In this paper, the benefits and challenges of the proposed concepts are tested in two forensic contexts: (1) ink identification and (2) ink evidential value assessment. The results show that different algorithms are better suited for different tasks. This research shows that it is possible to build digital ink libraries using the most commonly used ink analytical technique, i.e. high-performance thin layer chromatography, despite its reputation of lacking reproducibility. More importantly, it is possible to assign evidential value to ink evidence in a transparent way using a probabilistic model. It is therefore possible to move away from the traditional subjective approach, which is entirely based on experts' opinion, and which is usually not very informative. While there is room for the improvement, this report demonstrates the significant gains obtained over the traditional subjective approach for the search of ink specimens in ink databases, and the interpretation of their evidential value.

Niche, distribution and global change : modeling insights into biogeography and conservation biology Olivier Broennimann

RESUME Les évidences montrant que les changements globaux affectent la biodiversité s'accumulent. Les facteurs les plus influant dans ce processus sont les changements et destructions d'habitat, l'expansion des espèces envahissantes et l'impact des changements climatiques. Une évaluation pertinente de la réponse des espèces face à ces changements est essentielle pour proposer des mesures permettant de réduire le déclin actuel de la biodiversité. La modélisation de la répartition d'espèces basée sur la niche (NBM) est l'un des rares outils permettant cette évaluation. Néanmoins, leur application dans le contexte des changements globaux repose sur des hypothèses restrictives et demande une interprétation critique. Ce travail présente une série d'études de cas investiguant les possibilités et limitations de cette approche pour prédire l'impact des changements globaux. Deux études traitant des menaces sur les espèces rares et en danger d'extinction sont présentées. Les caractéristiques éco-géographiques de 118 plantes avec un haut degré de priorité de conservation sont revues. La prévalence des types de rareté sont analysées en relation avec leur risque d'extinction UICN. La revue souligne l'importance de la conservation à l'échelle régionale. Une évaluation de la rareté à échelle globale peut être trompeuse pour certaine espèces car elle ne tient pas en compte des différents degrés de rareté que présente une espèce à différentes échelles spatiales. La deuxième étude test une approche pour améliorer l'échantillonnage d'espèces rares en incluant des phases itératives de modélisation et d'échantillonnage sur le terrain. L'application de l'approche en biologie de la conservation (illustrée ici par le cas du chardon bleu, Eryngium alpinum), permettrait de réduire le temps et les coûts d'échantillonnage. Deux études sur l'impact des changements climatiques sur la faune et la flore africaine sont présentées. La première étude évalue la sensibilité de 227 mammifères africains face aux climatiques d'ici 2050. Elle montre qu'un nombre important d'espèces pourrait être bientôt en danger d'extinction et que les parcs nationaux africains (principalement ceux situé en milieux xériques) pourraient ne pas remplir leur mandat de protection de la biodiversité dans le futur. La seconde étude modélise l'aire de répartition en 2050 de 975 espèces de plantes endémiques du sud de l'Afrique. L'étude propose l'inclusion de méthodes améliorant la prédiction des risques liés aux changements climatiques. Elle propose également une méthode pour estimer a priori la sensibilité d'une espèce aux changements climatiques à partir de ses propriétés écologiques et des caractéristiques de son aire de répartition. Trois études illustrent l'utilisation des modèles dans l'étude des invasions biologiques. Une première étude relate l'expansion de la laitue sáuvage (Lactuca serriola) vers le nord de l'Europe en lien avec les changements du climat depuis 250 ans. La deuxième étude analyse le potentiel d'invasion de la centaurée tachetée (Centaures maculosa), une mauvaise herbe importée en Amérique du nord vers 1890. L'étude apporte la preuve qu'une espèce envahissante peut occuper une niche climatique différente après introduction sur un autre continent. Les modèles basés sur l'aire native prédisent de manière incorrecte l'entier de l'aire envahie mais permettent de prévoir les aires d'introductions potentielles. Une méthode alternative, incluant la calibration du modèle à partir des deux aires où l'espèce est présente, est proposée pour améliorer les prédictions de l'invasion en Amérique du nord. Je présente finalement une revue de la littérature sur la dynamique de la niche écologique dans le temps et l'espace. Elle synthétise les récents développements théoriques concernant le conservatisme de la niche et propose des solutions pour améliorer la pertinence des prédictions d'impact des changements climatiques et des invasions biologiques. SUMMARY Evidences are accumulating that biodiversity is facing the effects of global change. The most influential drivers of change in ecosystems are land-use change, alien species invasions and climate change impacts. Accurate projections of species' responses to these changes are needed to propose mitigation measures to slow down the on-going erosion of biodiversity. Niche-based models (NBM) currently represent one of the only tools for such projections. However, their application in the context of global changes relies on restrictive assumptions, calling for cautious interpretations. In this thesis I aim to assess the effectiveness and shortcomings of niche-based models for the study of global change impacts on biodiversity through the investigation of specific, unsolved limitations and suggestion of new approaches. Two studies investigating threats to rare and endangered plants are presented. I review the ecogeographic characteristic of 118 endangered plants with high conservation priority in Switzerland. The prevalence of rarity types among plant species is analyzed in relation to IUCN extinction risks. The review underlines the importance of regional vs. global conservation and shows that a global assessment of rarity might be misleading for some species because it can fail to account for different degrees of rarity at a variety of spatial scales. The second study tests a modeling framework including iterative steps of modeling and field surveys to improve the sampling of rare species. The approach is illustrated with a rare alpine plant, Eryngium alpinum and shows promise for complementing conservation practices and reducing sampling costs. Two studies illustrate the impacts of climate change on African taxa. The first one assesses the sensitivity of 277 mammals at African scale to climate change by 2050 in terms of species richness and turnover. It shows that a substantial number of species could be critically endangered in the future. National parks situated in xeric ecosystems are not expected to meet their mandate of protecting current species diversity in the future. The second study model the distribution in 2050 of 975 endemic plant species in southern Africa. The study proposes the inclusion of new methodological insights improving the accuracy and ecological realism of predictions of global changes studies. It also investigates the possibility to estimate a priori the sensitivity of a species to climate change from the geographical distribution and ecological proprieties of the species. Three studies illustrate the application of NBM in the study of biological invasions. The first one investigates the Northwards expansion of Lactuca serriola L. in Europe during the last 250 years in relation with climate changes. In the last two decades, the species could not track climate change due to non climatic influences. A second study analyses the potential invasion extent of spotted knapweed, a European weed first introduced into North America in the 1890s. The study provides one of the first empirical evidence that an invasive species can occupy climatically distinct niche spaces following its introduction into a new area. Models fail to predict the current full extent of the invasion, but correctly predict areas of introduction. An alternative approach, involving the calibration of models with pooled data from both ranges, is proposed to improve predictions of the extent of invasion on models based solely on the native range. I finally present a review on the dynamic nature of ecological niches in space and time. It synthesizes the recent theoretical developments to the niche conservatism issues and proposes solutions to improve confidence in NBM predictions of the impacts of climate change and species invasions on species distributions.

Accuracy of four cardiovascular risk scores in a Swiss population-based cohort

Purpose: To assess the global cardiovascular (CV) risk of an individual, several scores have been developed. However, their accuracy and comparability need to be evaluated in populations others from which they were derived. The aim of this study was to compare the predictive accuracy of 4 CV risk scores using data of a large population-based cohort. Methods: Prospective cohort study including 4980 participants (2698 women, mean age± SD: 52.7±10.8 years) in Lausanne, Switzerland followed for an average of 5.5 years (range 0.2 - 8.5). Two end points were assessed: 1) coronary heart disease (CHD), and 2) CV diseases (CVD). Four risk scores were compared: original and recalibrated Framingham coronary heart disease scores (1998 and 2001); original PROCAM score (2002) and its recalibrated version for Switzerland (IAS-AGLA); Reynolds risk score. Discrimination was assessed using Harrell's C statistics, model fitness using Akaike's information criterion (AIC) and calibration using pseudo Hosmer-Lemeshow test. The sensitivity, specificity and corresponding 95% confidence intervals were assessed for each risk score using the highest risk category ([20+ % at 10 years) as the "positive" test. Results: Recalibrated and original 1998 and original 2001 Framingham scores show better discrimination (>0.720) and model fitness (low AIC) for CHD and CVD. All 4 scores are correctly calibrated (Chi2<20). The recalibrated Framingham 1998 score has the best sensitivities, 37.8% and 40.4%, for CHD and CVD, respectively. All scores present specificities >90%. Framingham 1998, PROCAM and IAS-AGLA scores include the greatest proportion of subjects (>200) in the high risk category whereas recalibrated Framingham 2001 and Reynolds include <=44 subjects. Conclusion: In this cohort, we see variations of accuracy between risk scores, the original Framingham 2001 score demonstrating the best compromise between its accuracy and its limited selection of subjects in the highest risk category. We advocate that national guidelines, based on independently validated data, take into account calibrated CV risk scores for their respective countries.

A new method for the determination of plutonium and americium using high pressure microwave digestion and alpha-spectrometry or ICP-SMS

Plutonium and americium are radionuclides particularly difficult to measure in environmental samples because they are alpha-emitters and therefore necessitate a careful separation before any measurement, either using radiometric methods or ICP-SMS. Recent developments in extraction chromatography resins such as Eichrom (R) TRU and TEVA have resolved many of the analytical problems but drawbacks such as low recovery and spectral interferences still occasionally occur. Here, we report on the use of the new Eichrom (R) DGA resin in association with TEVA resin and high pressure microwave acid leaching for the sequential determination of plutonium and americium in environmental samples. The method results in average recoveries of 83 +/- 15% for plutonium and 73 +/- 22% for americium (n = 60), and a less than 10% deviation from reference values of four IAEA reference materials and three samples from intercomparisons exercises. The method is also suitable for measuring Pu-239 in water samples at the mu Bq/l level, if ICP-SMS is used for the measurement.

Does the physical disector method provide an accurate estimation of sensory neuron number in rat dorsal root ganglia?

The physical disector is a method of choice for estimating unbiased neuron numbers; nevertheless, calibration is needed to evaluate each counting method. The validity of this method can be assessed by comparing the estimated cell number with the true number determined by a direct counting method in serial sections. We reconstructed a 1/5 of rat lumbar dorsal root ganglia taken from two experimental conditions. From each ganglion, images of 200 adjacent semi-thin sections were used to reconstruct a volumetric dataset (stack of voxels). On these stacks the number of sensory neurons was estimated and counted respectively by physical disector and direct counting methods. Also, using the coordinates of nuclei from the direct counting, we simulate, by a Matlab program, disector pairs separated by increasing distances in a ganglion model. The comparison between the results of these approaches clearly demonstrates that the physical disector method provides a valid and reliable estimate of the number of sensory neurons only when the distance between the consecutive disector pairs is 60 microm or smaller. In these conditions the size of error between the results of physical disector and direct counting does not exceed 6%. In contrast when the distance between two pairs is larger than 60 microm (70-200 microm) the size of error increases rapidly to 27%. We conclude that the physical dissector method provides a reliable estimate of the number of rat sensory neurons only when the separating distance between the consecutive dissector pairs is no larger than 60 microm.

Fast quantification of ten psychotropic drugs and metabolites in human plasma by ultra-high performance liquid chromatography tandem mass spectrometry for therapeutic drug monitoring.

A sensitive and selective ultra-high performance liquid chromatography (UHPLC) tandem mass spectrometry (MS/MS) method was developed for the fast quantification of ten psychotropic drugs and metabolites in human plasma for the needs of our laboratory (amisulpride, asenapine, desmethyl-mirtazapine, iloperidone, mirtazapine, norquetiapine, olanzapine, paliperidone, quetiapine and risperidone). Stable isotope-labeled internal standards were used for all analytes, to compensate for the global method variability, including extraction and ionization variations. Sample preparation was performed by generic protein precipitation with acetonitrile. Chromatographic separation was achieved in less than 3.0min on an Acquity UPLC BEH Shield RP18 column (2.1mm×50mm; 1.7μm), using a gradient elution of 10mM ammonium formate buffer pH 3.0 and acetonitrile at a flow rate of 0.4ml/min. The compounds were quantified on a tandem quadrupole mass spectrometer operating in positive electrospray ionization mode, using multiple reaction monitoring. The method was fully validated according to the latest recommendations of international guidelines. Eight point calibration curves were used to cover a large concentration range 0.5-200ng/ml for asenapine, desmethyl-mirtazapine, iloperidone, mirtazapine, olanzapine, paliperidone and risperidone, and 1-1500ng/ml for amisulpride, norquetiapine and quetiapine. Good quantitative performances were achieved in terms of trueness (93.1-111.2%), repeatability (1.3-8.6%) and intermediate precision (1.8-11.5%). Internal standard-normalized matrix effects ranged between 95 and 105%, with a variability never exceeding 6%. The accuracy profiles (total error) were included in the acceptance limits of ±30% for biological samples. This method is therefore suitable for both therapeutic drug monitoring and pharmacokinetic studies.

Phylogenetic relationships and divergence times among dormice (Rodentia, Gliridae) based on three nuclear genes

We examined phylogenetic relationships among six species representing three subfamilies, Glirinae, Graphiurinae and Leithiinae with sequences from three nuclear protein-coding genes (apolipoprotein B, APOB; interphotoreceptor retinoid-binding protein, IRBP; recombination-activating gene 1, RAG1). Phylogenetic trees reconstructed from maximum-parsimony (MP), maximum-likelihood (ML) and Bayesian-inference (BI) analyses showed the monophyly of Glirinae (Glis and Glirulus) and Leithiinae (Dryomys, Eliomys and Muscardinus) with strong support, although the branch length maintaining this relationship was very short, implying rapid diversification among the three subfamilies. Divergence time estimates were calculated from ML (local clock model) and Bayesian-dating method using a calibration point of 25 Myr (million years) ago for the divergence between Glis and Glirulus, and 55 Myr ago for the split between lineages of Gliridae and Sciuridae on the basis of fossil records. The results showed that each lineage of Graphiuros, Glis, Glirulus and Muscardinus dates from the Late Oligocene to the Early Miocene period, which is mostly in agreement with fossil records. Taking into account that warm climate harbouring a glirid-favoured forest dominated from Europe to Asia during this period, it is considered that this warm environment triggered the prosperity of the glirid species through the rapid diversification. Glirulus japonicas is suggested to be a relict of this ancient diversification during the warm period.

Stable isotopes in lake Geneva carbonate sediments and molluscs: Review and new data

New isotopic results on bulk carbonate and mollusc (gastropods and bivalves) samples from Lake Geneva (Switzerland), spanning the period from the Oldest Dryas to the present day, are compared with pre-existing stable isotope data. According to preliminary calibration of modern samples, Lake Geneva endogenic calcite precipitates at or near oxygen isotopic equilibrium with ambient water, confirming the potential of this large lake to record paleoenvironmental and paleoclimatic changes. The onset of endogenic calcite precipitation at the beginning of the Allerod biozone is clearly indicated by the oxygen isotopic signature of bulk carbonate. A large change in delta(13)C values occurs during the Preboreal. This carbon shift is likely to be due to a change in bioproductivity and/or to a `'catchment effect'', the contribution of biogenic CO2 from the catchment area to the dissolved inorganic carbon reservoir of the lake water becoming significant only during the Preboreal. Gastropods are confirmed as valuable for studies of changes in paleotemperature and in paleowater isotopic composition, despite the presence of a vital effect. Mineralogical evidence indicates an increased detrital influence upon sedimentation since the Subboreal time period. On the other hand, stable isotope measurements of Subatlantic carbonate sediments show values comparable to those of pure endogenic calcite and of gastropods (taking into account the vital effect). This apparent disagreement still remains difficult to explain.

Determination of picogram levels of midazolam, and 1- and 4-hydroxymidazolam in human plasma by gas chromatography-negative chemical ionization-mass spectrometry.

Midazolam is a widely accepted probe for phenotyping cytochrome P4503A. A gas chromatography-mass spectrometry (GC-MS)-negative chemical ionization method is presented which allows measuring very low levels of midazolam (MID), 1-OH midazolam (1OHMID) and 4-OH midazolam (4OHMID), in plasma, after derivatization with the reagent N-tert-butyldimethylsilyl-N-methyltrifluoroacetamide. The standard curves were linear over a working range of 20 pg/ml to 5 ng/ml for the three compounds, with the mean coefficients of correlation of the calibration curves (n = 6) being 0.999 for MID and 1OHMID, and 1.0 for 4OHMID. The mean recoveries measured at 100 pg/ml, 500 pg/ml, and 2 ng/ml, ranged from 76 to 87% for MID, from 76 to 99% for 1OHMID, from 68 to 84% for 4OHMID, and from 82 to 109% for N-ethyloxazepam (internal standard). Intra- (n = 7) and inter-day (n = 8) coefficients of variation determined at three concentrations ranged from 1 to 8% for MID, from 2 to 13% for 1OHMID and from 1 to 14% for 4OHMID. The percent theoretical concentrations (accuracy) were within +/-8% for MID and 1OHMID, within +/-9% for 4OHMID at 500 pg/ml and 2 ng/ml, and within +/-28% for 4OHMID at 100 pg/ml. The limits of quantitation were found to be 10 pg/ml for the three compounds. This method can be used for phenotyping cytochrome P4503A in humans following the administration of a very low oral dose of midazolam (75 microg), without central nervous system side-effects.

Determination of meropenem in plasma and filtrate-dialysate from patients under continuous veno-venous haemodiafiltration by SPE-LC.

Meropenem, a carbapenem antibiotic displaying a broad spectrum of antibacterial activity, is administered in Medical Intensive Care Unit to critically ill patients undergoing continuous veno-venous haemodiafiltration (CVVHDF). However, there are limited data available to substantial rational dosing decisions in this condition. In an attempt to refine our knowledge and propose a rationally designed dosage regimen, we have developed a HPLC method to determine meropenem after solid-phase extraction (SPE) of plasma and dialysate fluids obtained from patients under CVVHDF. The assay comprises the simultaneous measurement of meropenem's open-ring metabolite UK-1a, whose fate has never been studied in CVVHDF patients. The clean-up procedure involved a SPE on C18 cartridge. Matrix components were eliminated with phosphate buffer pH 7.4 followed by 15:85 MeOH-phosphate buffer pH 7.4. Meropenem and UK-1a were subsequently desorbed with MeOH. The eluates were evaporated under nitrogen at room temperature (RT) and reconstituted in phosphate buffer pH 7.4. Separation was performed at RT on a Nucleosil 100-5 microm C18 AB cartridge column (125 x 4 mm I.D.) equipped with a guard column (8 x 4 mm I.D.) with UV-DAD detection set at 208 nm. The mobile phase was 1 ml min(-1), using a step-wise gradient elution program: %MeOH/0.005 M tetrabutylammonium chloride pH 7.4; 10/90-50/50 in 27 min. Over the range of 5-100 microg ml(-1), the regression coefficient of the calibration curves (plasma and dialysate) were >0.998. The absolute extraction recoveries of meropenem and UK-1a in plasma and filtrate-dialysate were stable and ranged from 88-93 to 72-77% for meropenem, and from 95-104 to 75-82% for UK-1a. In plasma and filtrate-dialysate, respectively, the mean intra-assay precision was 4.1 and 2.6% for meropenem and 4.2 and 3.7% for UK-1a. The inter-assay variability was 2.8 and 3.6% for meropenem and 2.3 and 2.8% for UK-1a. The accuracy was satisfactory for both meropenem and UK-1a with deviation never exceeding 9.0% of the nominal concentrations. The stability of meropenem, studied in biological samples left at RT and at +4 degrees C, was satisfactory with < 5% degradation after 1.5 h in blood but reached 22% in filtrate-dialysate samples stored at RT for 8 h, precluding accurate measurements of meropenem excreted unchanged in the filtrate-dialysate left at RT during the CVVHDF procedure. The method reported here enables accurate measurements of meropenem in critically ill patients under CVVHDF, making dosage individualisation possible in such patients. The levels of the metabolite UK-1a encountered in this population of patients were higher than those observed in healthy volunteers but was similar to those observed in patients with renal impairment under hemodialysis.

A LC-tandem MS assay for the simultaneous measurement of new antiretroviral agents: Raltegravir, maraviroc, darunavir, and etravirine.

Raltegravir (RAL), maraviroc (MVC), darunavir (DRV), and etravirine (ETV) are new antiretroviral agents with significant potential for drug interactions. This work describes a sensitive and accurate liquid chromatography-tandem mass spectrometry (LC-MS/MS) method for the determination of plasma drug levels. Single-step extraction of RAL, MVC, DRV, ETV and RTV from plasma (100 microl) is performed by protein precipitation using 600 microl of acetonitrile, after the addition of 100 microl darunavir-d(9) (DRV-d(9)) at 1000 ng/ml in MeOH/H(2)O 50/50 as internal standard (I.S.). The mixture is vortexed, sonicated for 10 min, vortex-mixed again and centrifuged. An aliquot of supernatant (150 microl) is diluted 1:1 with a mixture of 20 mM ammonium acetate/MeOH 40/60 and 10 microl is injected onto a 2.1 x 50 mm Waters Atlantis-dC18 3 microm analytical column. Chromatographic separations are performed using a gradient program with 2 mM ammonium acetate containing 0.1% formic acid and acetonitrile with 0.1% formic acid. Analytes quantification is performed by electrospray ionisation-triple quadrupole mass spectrometry using the selected reaction monitoring detection in the positive mode. The method has been validated over the clinically relevant concentrations ranging from 12.5 to 5000 ng/ml, 2.5 to 1000 ng/ml, 25 to 10,000 ng/ml, 10 to 4000 ng/ml, and 5 to 2000 ng/ml for RAL, MRV, DRV, ETV and RTV, respectively. The extraction recovery for all antiretroviral drugs is always above 91%. The method is precise, with mean inter-day CV% within 5.1-9.8%, and accurate (range of inter-day deviation from nominal values -3.3 to +5.1%). In addition our method enables the simultaneous assessment of raltegravir-glucuronide. This is the first analytical method allowing the simultaneous assay of antiretroviral agents targeted to four different steps of HIV replication. The proposed method is suitable for the Therapeutic Drug Monitoring Service of these new regimen combinations administered as salvage therapy to patients having experienced treatment failure, and for whom exposure, tolerance and adherence assessments are critical.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

A validated assay by liquid chromatography-tandem mass spectrometry for the simultaneous quantification of elvitegravir and rilpivirine in HIV positive patients.

Because of the large variability in the pharmacokinetics of anti-HIV drugs, therapeutic drug monitoring in patients may contribute to optimize the overall efficacy and safety of antiretroviral therapy. An LC-MS/MS method for the simultaneous assay in plasma of the novel antiretroviral agents rilpivirine (RPV) and elvitegravir (EVG) has been developed to that endeavor. Plasma samples (100 μL) extraction is performed by protein precipitation with acetonitrile, and the supernatant is subsequently diluted 1:1 with 20-mM ammonium acetate/MeOH 50:50. After reverse-phase chromatography, quantification of RPV and EVG, using matrix-matched calibration samples, is performed by electrospray ionization-triple quadrupole mass spectrometry by selected reaction monitoring detection using the positive mode. The stable isotopic-labeled compounds RPV-(13) C6 and EVG-D6 were used as internal standards. The method was validated according to FDA recommendations, including assessment of extraction yield, matrix effects variability (<6.4%), as well as EVG and RPV short and long-term stability in plasma. Calibration curves were validated over the clinically relevant concentrations ranging from 5 to 2500 ng/ml for RPV and from 50 to 5000 ng/ml for EVG. The method is precise (inter-day CV%: 3-6.3%) and accurate (3.8-7.2%). Plasma samples were found to be stable (<15%) in all considered conditions (RT/48 h, +4°C/48 h, -20°C/3 months and 60°C/1 h). Selected metabolite profiles analysis in patients' samples revealed the presence of EVG glucuronide, that was well separated from parent EVG, allowing to exclude potential interferences through the in-source dissociation of glucuronide to parent drug. This new, rapid and robust LCMS/MS assay for the simultaneous quantification of plasma concentrations of these two major new anti-HIV drugs EVG and RPV offers an efficient analytical tool for clinical pharmacokinetics studies and routine therapeutic drug monitoring service. Copyright © 2013 John Wiley & Sons, Ltd.

Sensitive determination of D-lactic acid and L-lactic acid in urine by high-performance liquid chromatography-tandem mass spectrometry.

D-lactic acid in urine originates mainly from bacterial production in the intestinal tract. Increased D-lactate excretion as observed in patients affected by short bowel syndrome or necrotizing enterocolitis reflects D-lactic overproduction. Therefore, there is a need for a reliable and sensitive method able to detect D-lactic acid even at subclinical elevation levels. A new and highly sensitive method for the simultaneous determination of L- and D-lactic acid by a two-step procedure has been developed. This method is based on the concentration of lactic acid enantiomers from urine by supported liquid extraction followed by high-performance liquid chromatography-tandem mass spectrometry. The separation was achieved by the use of an Astec Chirobiotic? R chiral column under isocratic conditions. The calibration curves were linear over the ranges of 2-400 and 0.5-100 µmol/L respectively for L- and D-lactic acid. The limit of detection of D-lactic acid was 0.125 µmol/L and its limit of quantification was 0.5 µmol/L. The overall accuracy and precision were well within 10% of the nominal values. The developed method is suitable for production of reference values in children and could be applied for accurate routine analysis.

Fluid-rock interaction during late stages of the Alpine exhumation

Résumé de la thèseLa fracturation des roches au cours de phases compressives ou extensives est un souvent évoquée pour expliquer la circulation de fluide au sein des roches cristallines. Dans le cadre de cette thèse, la circulation des fluides lors de l'exhumation tardive des Alpes a été étudiée en utilisant deux approches différentes: analyses structurales de la déformation fragile d'une part et analyses géochimiques des roches et des minéraux (isotopes stables, datations U/Pb, thermochronologie (U-Th)/He) d'autre part. Cette approche combinée a permis de mieux comprendre l'interaction existante entre les fluides métamorphiques et les fluides météoriques, ainsi que leur interaction avec les roches encaissantes. Le travail a été effectué dans la zone Pennique du Valais suisse.La première partie était focalisée sur la déformation fragile, le but étant de définir les différents types de déformations existantes et de déterminer l'âge relatif des différentes familles de failles. Dans la région d'étude, quatre domaines ont été distingués. Chacun d'eux comportent deux types de structures fragiles, certaines sont minéralisées alors que d'autre non. Au sein de chaque domaine, la direction principale des structures minéralisées correspond à l'orientation des accidents tectoniques majeurs de la région (Aosta- Ranzola Line au Sud, Rhône Line au Nord et Simplon Fault Zone à l'Est), alors que les structures non- minéralisées montrent des orientations plus variables. Ainsi, le premier type de structure est interprété comme résultant d'une dislocation tectonique alors que le deuxième type de structure résulterait d'une dislocation gravitaire locale. Il n'est néanmoins pas possible de classer chronologiquement la formation de ces deux types de structure ni d'attribuer un âge relatif aux changements d'orientation des contraintes majeures.La deuxième étude a été effectuée dans la région de la zone de faille du Simplon. Dans cette zone, la composition isotopique des minéraux ayant cristallisé à l'intérieur des fractures tardives permet de distinguer différents types de circulation de fluide. Les valeurs δ180 du quartz de la roche encaissante ainsi que ceux des veines tardives du bloque inférieur de la faille sont comparables. Ces valeurs indiquent un rééquilibrage et un tamponnage isotopique des fluides tardifs au contact de la roche encaissante lors de la fracturation de cette dernière et de la cristallisation des veines tardives. La même situation est observée dans la partie nord du bloque supérieur ainsi que dans sa partie sud. Ceci n'est néanmoins pas le cas pour la partie centrale du bloque supérieur où les valeurs isotopiques des minéraux dans les veines tardives sont approximativement 3 %o plus basses (avec des valeurs extrêmes négatifs), indiquant une contribution d'eau météorique aux fluides circulant dans les veines. Ces données suggèrent qu'une infiltration d'eau météorique a pu avoir lieu dans le bloque supérieur, où la fracturation des roches est plus intensive car le déplacement relatif le long de la faille y fut plus important, et la température maximale du métamorphisme plus basse. La troisième contribution traite de la géo-thermochronologie de la zone de contact entre la klippe de la Dent Blanche et la nappe de Tsaté. De petits zircons euhédraux ont été trouvés dans un plan de faille minéralisé (parallèle à la Faille du Rhône, voir première partie de l'étude), riche en hématite et quartz, de la zone d'étude. Les analyses U/Pb donnent des âges radiométriques autour de 270 - 280 Ma aux zircons extraits de la minéralisation ainsi que ceux extraits de la roche encaissante, ce qui correspond à l'âge de la nappe de la Dent Blanche et non celui de la nappe du Tsaté qui est elle-même classiquement interprétée comme une ophiolite Jurassique de l'Océan Liguro-Piémontais. Ces données suggèrent que les zircons contenus dans la veine ont été hérités de la roche encaissante. Les résultats (U-Th)/He indiquent un âge de refroidissement différent pour la roche encaissante (25.5 ± 2.0 Ma) que celui de la minéralisation (17.7 ±1.4 Ma). Le thermomètre isotopique quartz-hématite indique une température d'équilibre, et donc de mise en place de la minéralisation, d'environ 170 °C, température très proche de la température de -180 °C de fermeture du zircon pour le système (U-Th)/He. Ceci suggère que l'âge de refroidissement des zircons de la minéralisation correspond aussi à l'âge de formation de la faille.Thesis abstractFluid circulation in fractured rocks is a common process in geology, and it is generally the consequence of faulting and fracturing during both tectonic compression and extension. This thesis is focused on fluid circulation during late stages of the Alpine exhumation. After a structural analysis of the late brittle deformation of the studied samples, several analytical methods (stable isotope investigations, U/Pb radiometric dating, (U-Th)/He thermochronology) have been applied to understand the interaction of metamorphic and meteoric fluids with one another as well as with the host rock. This thesis is articulated around three study directions. All studies were conducted in the Penninic Zone of the Valais, Switzerland. The first study deals with late, brittle deformation and focuses on the different deformation styles and on the relative age of the different families of fractures. In order to do this, late brittle structures observed in four different domains have been subdivided as a function of the existence (or not) and type of mineralization. Comparisons between mineralized and non-mineralized strike directions for all four domains show that mineralized structures follow the strike orientation of major tectonic movements indicated in the Penninic Zone of the Valais (Aosta-Ranzola Line to the S, Rhône Line to the Ν and Simplon Fault Zone to the E), whereas non-mineralized fractures have a more variable strike orientation. This difference could be interpreted as indicative of tectonic-related faulting (mineralized structures) vs. local, collapse-related faulting (non-mineralized fractures), but it is not strong enough to indicate a relative age of the late brittle structures, and/or a change in the orientation of the strain field in post-Miocene times. The second studied area is focused on the Simplon Fault Zone (SFZ). Stable isotope analyses of minerals filling these late fractures indicate that there are two different fluid circulation systems in the footwall and hanging wall of the SFZ. In the footwall, δ180 values of quartz from both the host rock and the late veins range from +10 %o to +12 %o. This is consistent with buffering of circulating fluids by the host rock during fracturing and vein precipitation. In the hanging wall, δΙ80 values for quartz crystals from the host rock and the late veins are similar in both the northern and southern parts of the detachment that are both affected by the same degree of metamorphism (greenschist to the Ν and amphibolite to the S). This is not the case in the central part of the SFZ, where there is a jump from amphibolite facies in the footwall to greenschist facies in the hanging wall. δ,80 values for quartz from the hanging wall late veins are approximately 3.0 %o lower (down to negative values in some cases) than the values observed in the footwall These data suggest that infiltration of meteoric water may have occurred in the most fractured parts of the hanging wall, where relative displacement on the SFZ was the greatest and the peak temperature lower. In the less fractured footwall the δ180 values reflect a host rock-buffered system.The third study is focused on geo-thermochronology at the contact between the Dent Blanche nappe and the Tsaté nappe where small, euhedral zircons were found in a hematite- and quartz-rich mineralization on a late normal fault plane parallel to the Rhône Line (see first part of the study). U/Pb analysis indicates that the zircons - both in the late mineralization and in the host rock - have absolute radiometric ages clustering around 270 - 280 Ma, which is the accepted age for intrusive rocks from the Austroalpine Dent Blanche units but not for the Tsaté nappe. The latter is classically interpreted as an ophiolitic remnant of the Jurassic Liguro-Piemontais Ocean. U/Pb analyses suggest that zircons in late mineralization are all inherited from the host rock; however, results of (U-Th)/He analyses indicate that cooling ages for the host rocks are different to the cooling ages for the zircons in late mineralization. Indeed, the calculated cooling age for the Arolla gneiss is 25.5 ± 2.0 Ma, whilst the cooling age for the associated mineralized fault plane is 17.7 ±1.4 Ma. Oxygen stable isotope fractionation between quartz and hematite in the same late mineralization corresponds to temperatures of about 170 °C. The proximity of the calculated emplacement temperature for the mineralization and the lower accepted closure temperature for zircon in the (U-Th)/He system (-180 °C) imply that the age of 17.7 ± 1.4 Ma can also be interpreted as the formation age of this late brittle fault.Résumé grand publicLa circulation des fluides dans les roches fracturées est typique de nombreux processus géologiques, et très souvent est la conséquence de la fracturation des roches. Cette thèse aborde la question de la circulation des fluides pendant les dernières phases du soulèvement des Alpes. Après une analyse structurale de la fracturation directement sur le terrain, plusieurs méthodes géochimiques ont été appliquées pour comprendre l'interaction entre les différents fluides circulants, et avec leur propre roche mère. L'étude, concentrée sur trois directions principales, a été conduite dans la zone Pennique du Valais suisse. La première partie traite de la déformation cassante dans le secteur cité. L'analyse détaillée des fractures a permis de les subdiviser en structures minéralisées et non-minéralisées, sur quatre domaines différents. La comparaison entre les directions des structures minéralisées et non-minéralisées a permis de montrer que les premières suivent l'orientation des accidents tectoniques majeurs de la région, alors que les structures non- minéralisées ont une orientation plus variable. Cette différence pourrait être interprétée comme indication d'une dislocation tectonique (structures minéralisées) contre une dislocation gravitaire locale (structures non-minéralisées), mais elle n'est pas assez forte pour indiquer un âge relatif des structures tardives et/ou un changement de l'orientation des contraintes après -20 Ma vers le présent.A partir de ces observations, la deuxième étude est concentrée dans la région de la faille du Simplon. Les analyses géochimiques sur les minéraux remplissant les structures tardives indiquent qu'il y a deux différents systèmes de circulation des fluides dans les deux parties (toit et mur) de la faille. Dans le mur, les valeurs isotopiques des minéraux cristallisés à partir d'un fluide tardif sont les mêmes de ceux de la roche mère, donc il y a eu rééquilibration chimique entre fluide et roche pendant la fracturation de cette dernière et la précipitation des minéraux. Dans le toit, les valeurs isotopiques dans la roche mère et dans les minéraux des veines tardives sont comparables dans les parties Ν et S de la faille, où les roches du toit et du mur ont atteint une température maximale - pendant phase prograde de la formation des Alpes - comparable. Au contraire, dans la partie centrale, où le mur a atteint des températures maximales plus élevées par rapport au toit, les valeurs géochimiques des minéralisations tardives du toit sont parfois plus basses que les valeurs observées dans le mur. Ces données suggèrent que l'infiltration de l'eau de surface aurait pu se produire dans la partie plus fracturée du toit, où le déplacement relatif le long de la faille était majeur et les températures maximales mineures. Au contraire, les données géochimiques du mur de la partie centrale indiquent un système isotopique équilibré par la roche mère.La troisième partie de ce travail se base sur l'étude géochimique intégrée des isotopes stables d'Oxygène et radioactifs du Plomb, Uranium, Thorium et Hélium, auprès d'une faille normale minéralisée et des roches de la région à cheval entre deux nappes, la nappe de la Dent Blanche et la nappe de Tsaté. Ici, des petits zircons ont été trouvés dans la minéralisation citée, riche en hématite et quartz. L'analyse radiométrique Uranium/Plomb a montré que les zircons dans la minéralisation et dans les roches autour ont des âges comparables (autour 280 Ma). Cela signifie que les zircons dans la minéralisation tardive ont été hérités de la roche mère pendant la fracturation et la circulation des fluides tardives. De l'autre coté, les résultats des analyses Uranium-Thorium/Hélium indiquent que les âges de refroidissement pour les roches mères sont différents comparés aux âges de refroidissement pour les zircons dans la minéralisation tardive: ces derniers sont plus jeunes d'environ 8 Ma (autour 25 Ma et autour 17 Ma respectivement). Les analyses des isotopes de l'oxygène sur quartz et hématite dans la même minéralisation donnent une température de mise en place de cette dernière d'environ 170° C. La température de fermeture du système chimique des zircons dans le système (Uranium-Thorium)/Hélium est d'environ 180 °C: la proximité de ces deux températures implique que l'âge de refroidissement de la minéralisation tardive peut également être interprété comme âge de formation de la faille.