Predicting the deleterious effects of mutation load in fragmented populations.

Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure (F(ST)) and within-deme heterozygosity (H(S)). A similar part of variance in median time to extinction was explained by a combination of local population size (N) and heterozygosity (H(S)). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.

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.

In vivo gene transfer into the ocular ciliary muscle mediated by ultrasound and microbubbles.

This study aimed to assess application of ultrasound (US) combined with microbubbles (MB) to transfect the ciliary muscle of rat eyes. Reporter DNA plasmids encoding for Gaussia luciferase, β-galactosidase or the green fluorescent protein (GFP), alone or mixed with 50% Artison MB, were injected into the ciliary muscle, with or without US exposure (US set at 1 MHz, 2 W/cm(2), 50% duty cycle for 2 min). Luciferase activity was measured in ocular fluids at 7 and 30 days after sonoporation. At 1 week, the US+MB treatment showed a significant increase in luminescence compared with control eyes, injected with plasmid only, with or without MB (×2.6), and, reporter proteins were localized in the ciliary muscle by histochemical analysis. At 1 month, a significant decrease in luciferase activity was observed in all groups. A rise in lens and ciliary muscle temperature was measured during the procedure but did not result in any observable or microscopic damages at 1 and 8 days. The feasibility to transfer gene into the ciliary muscle by US and MB suggests that sonoporation may allow intraocular production of proteins for the treatment of inflammatory, angiogenic and/or degenerative retinal diseases.

Markers for sepsis diagnosis in the forensic setting: state of the art.

Reliable diagnoses of sepsis remain challenging in forensic pathology routine despite improved methods of sample collection and extensive biochemical and immunohistochemical investigations. Macroscopic findings may be elusive and have an infectious or non-infectious origin. Blood culture results can be difficult to interpret due to postmortem contamination or bacterial translocation. Lastly, peripheral and cardiac blood may be unavailable during autopsy. Procalcitonin, C-reactive protein, and interleukin-6 can be measured in biological fluids collected during autopsy and may be used as in clinical practice for diagnostic purposes. However, concentrations of these parameters may be increased due to etiologies other than bacterial infections, indicating that a combination of biomarkers could more effectively discriminate non-infectious from infectious inflammations. In this article, we propose a review of the literature pertaining to the diagnostic performance of classical and novel biomarkers of inflammation and bacterial infection in the forensic setting.

Early diagenesis of bone and tooth apatite in fluvial and marine settings: Constraints from combined oxygen isotope, nitrogen and REE analysis

Fossil bones and teeth of Late Pleistocene terrestrial mammals from Rhine River gravels (RS) and the North Sea (NS), that have been exposed to chemically and isotopically distinct diagenetic fluids (fresh water versus seawater), were investigated to study the effects of early diagenesis on biogenic apatite. Changes in phosphate oxygen isotopic composition (delta O-18(PO4)), nitrogen content (wt.% N) and rare earth element (REE) concentrations were measured along profiles within bones that have not been completely fossilized, and in skeletal tissues (bone, dentine, enamel) with different susceptibilities to diagenetic alteration. Early diagenetic changes of elemental and isotopic compositions of apatite in fossil bone are related to the loss of the stabilizing collagen matrix. The REE concentration is negatively correlated with the nitrogen content, and therefore the amount of collagen provides a sensitive proxy for early diagenetic alteration. REE patterns of RS and NS bones indicate initial fossilization in a fresh water fluid with similar REE compositions. Bones from both settings have nearly collagen-free, REE-, U-, F- and Sr-enriched altered outer rims, while the collagen-bearing bone compacta in the central part often display early diagenetic pyrite void-fillings. However, NS bones exposed to Holocene seawater have outer rim delta O-18(PO4) values that are 1.1 to 2.6 parts per thousand higher compared to the central part of the same bones (delta O-18(PO4) = 18.2 +/- 0.9 parts per thousand, n = 19). Surprisingly, even the collagen-rich bone compacta with low REE contents and apatite crystallinity seems altered, as NS tooth enamel (delta O-18(PO4) =15.0 +/- 0.3 parts per thousand, n=4) has about 3%. lower delta O-18(PO4) values, values that are also similar to those of enamel from RS teeth. Therefore, REE concentration, N content and apatite crystallinity are in this case only poor proxies for the alteration of delta O-18(PO4) values. Seawater exposure of a few years up to 8 kyr can change the delta O-18(PO4) values of the bone apatite by > 3 parts per thousand. Therefore, bones fossilized in marine settings must be treated with caution for palaeoclimatic reconstructions. However, enamel seems to preserve pristine delta O-18(PO4) values on this time scale. Using species-specific calibrations for modern mammals, a mean delta O-18(H2O) value can be reconstructed for Late Pleistocene mammalian drinking water of around -9.2 +/- 0.5 parts per thousand, which is similar to that of Late Pleistocene groundwater from central Europe. (c) 2008 Elsevier B.V. All rights reserved.

Origin of CO2 and carbonate veins in mantle-derived xenoliths in the Pannonian Basin

The origin and evolution of CO2 inclusions and calcite veins in peridotite xenoliths of the Pannonian Basin, Hungary, were investigated by means of petrographic investigation and stable isotope analyses. The fluid inclusions recovered in paragenetic olivine and clinopyroxene belong to distinct populations: type A (texturally early) inclusions with regular shapes (often with negative crystal forms) forming intragranular trails, type B (texturally late) inclusions defining randomly oriented trails that reach grain boundaries Type B inclusions are often associated with silicate melt (type C) inclusions Stable carbon isotope compositions in inclusion-hosted CO2 were obtained by vacuum crushing followed by conventional dual inlet as well as continuous flow mass spectrometry in order to eliminate possible lab artifacts. Olivines, clino- and orthopyroxenes of the host peridotite have oxygen isotope compositions from 5.3 to 6.0 parts per thousand (relative to V-SMOW). without any relationship with xenolith texture. Some of the xenoliths contained calcite in various forms veins and infillings in silicate globules in veins, secondary carbonate veins filling cracks and metasomatic veins with diffuse margins The former two carbonate types have delta C-13 values around -13 parts per thousand (relative to V-PDB) and low Sr contents (<05 wt %), whereas the third type,veins with high-temperature metasomatic features have a delta C-13 value of -5 0 parts per thousand and high Sr contents up to 34 wt.% In spite of the mantle-like delta C-13 value and the unusually high Sr content typical for mantle-derived carbonate, trace element compositions have proven a crustal origin. This observation supports the conclusions of earlier studies that the carbonate melt droplets found on peridotite xenoliths in the alkaline basalts represent mobilized sedimentary carbonate. The large delta C-13 range and the C-12-enrichment in the carbonates can be attributed to devolanlization of the migrating carbonate or infiltration of surficial fluids containing C-12-rich dissolved carbon Carbon isotope compositions of inclusion-hosted CO2 range from -17 8 to -4.8 parts per thousand (relative to V-PDB) with no relation to the amount of CO2 released by vacuum crushing. Low-delta C-13 values measured by stepwise heating under vacuum suggest that the carbon component is pristine and not related to surficial contamination, and that primary mantle fluids with delta C-13 values around -5 parts per thousand were at least partly preserved in the xenoliths Tectonic reworking and heating by the basaltic magma resulted in partial CO2 release and local C-13-depletion. (C) 2010 Elsevier B V All rights reserved

The unsolved puzzle of neuropathogenesis in glutaric aciduria type I.

Glutaric aciduria type I (GA-I) is a cerebral organic aciduria caused by deficiency of glutaryl-Co-A dehydrogenase (GCDH). GCDH deficiency leads to accumulation of glutaric acid (GA) and 3-hydroxyglutaric acid (3-OHGA), two metabolites that are believed to be neurotoxic, in brain and body fluids. The disorder usually becomes clinically manifest during a catabolic state (e.g. intercurrent illness) with an acute encephalopathic crisis that results in striatal necrosis and in a permanent dystonic-dyskinetic movement disorder. The results of numerous in vitro and in vivo studies have pointed to three main mechanisms involved in the metabolite-mediated neuronal damage: excitotoxicity, impairment of energy metabolism and oxidative stress. There is evidence that during a metabolic crisis GA and its metabolites are produced endogenously in the CNS and accumulate because of limiting transport mechanisms across the blood-brain barrier. Despite extensive experimental work, the relative contribution of the proposed pathogenic mechanisms remains unclear and specific therapeutic approaches have yet to be developed. Here, we review the experimental evidence and try to delineate possible pathogenetic models and approaches for future studies.

Lack of MRI neurohypophyseal bright signal in a child with congenital nephrogenic diabetes insipidus

Congenital nephrogenic diabetes insipidus (CNDI) is a rare disease characterized by the inability of the kidney to respond to arginine vasopressin (AVP). The absence of the neurohypophyseal 'bright signal' on T1 sequence magnetic resonance imaging (MRI) is considered as an argument in favour of the diagnosis of central diabetes insipidus (CDI). This observation is challenged as we hereby present a case of a child diagnosed with CNDI and who did not present MRI pituitary bright signal. A 6-month-old male presented with failure to thrive, polyuria and polydypsia. Family history revealed that the mother, 35 years of age, had been presenting polydypsia and polyuria, and she was investigated at the age of 6 years with no concluding diagnosis. The patient's physical exam showed a weight of 5215 g (−3 DS) and clinical signs of dehydration. The patient's plasma sodium level was 155 mmol/L, osmolality 305 mOsm/kg and urine osmolality 150 mOsm/kg. Brain MRI showed in T1 sequences the absence of the posterior pituitary bright signal suggesting the diagnosis of CDI (Figure 1). The child was treated with synthetic AVP analogue 1-desamino-8-d-arginine vasopressin (DDAVP) without improvement, which led to the consideration of CNDI. The diagnosis was confirmed by an elevated serum level of AVP of 214 pmol/L (reference value ≤4.34 pmol/L) and by genetic analysis demonstrating and T106C mutation in the V2R (X-linked CNDI). The child was treated with thiazide diuretic and increased fluids with restricted sodium intake. This resulted in catch-up growth and improved neurological development. A follow-up MRI was performed 6 months after the start of therapy with the same technique. At that time, the child's weight had improved to 9310 g (−1.5 DS) corresponding to a gain of 22 g per day, and he did not present any clinical signs of dehydration and had a normal plasma level of sodium (140 mmol/L). MRI showed that the bright signal was still absent.

Fluid migration through thrust faults in the Helvetic nappes (Western Swiss Alps)

The delta(18)O, delta(13)C and Sr-87/Sr-86 values of calcite and organic matter were measured for carbonate mylonites from numerous thrusts in the Helvetic Alps. Carbonate mylonites in most of the thrusts retain essentially unaltered protolith delta(18)O and delta(13)C values, consistent with there having been little to no advection of isotopically distinct fluid through these faults. Only carbonate mylonites from the basal thrusts of the largest nappes have delta(18)O and/or delta(13)C values that differ from those of their protoliths. The zone of isotopic alteration/exchange is confined to c. 10 to 20 meters of these fault contacts. We propose the fluids that migrated through these faults contained variable amounts of organically derived carbon and radiogenic strontium, and were probably derived from dewatering of the sedimentary rocks and prograde metamorphic reactions in the nappes' root zones. Apart from the basal thrusts of the largest nappes that behaved as narrow, laterally extensive conduits for fluids, there is little isotopic evidence that large quantities of fluids passed through most of the carbonate-hosted thrusts in the Helvetic Alps.

Spontaneous aggregation and global polar ordering in squirmer suspensions

We have developed numerical simulations of three dimensional suspensions of active particles to characterize the capabilities of the hydrodynamic stresses induced by active swimmers to promote global order and emergent structures in active suspensions. We have considered squirmer suspensions embedded in a fluid modeled under a Lattice Boltzmann scheme. We have found that active stresses play a central role to decorrelate the collective motion of squirmers and that contractile squirmers develop significant aggregates.

The release of antidiuretic hormone is appropriate in response to hypovolemia and/or sodium administration in children with severe head injury: a trial of lactated Ringer's solution versus hypertonic saline.

We conducted an open, randomized, and prospective study to determine the effect of hypertonic saline on the secretion of antidiuretic hormone (ADH) and aldosterone in children with severe head injury (Glasgow coma scale <8). Thirty-one consecutive patients at a level III pediatric intensive care unit at a children's hospital received either lactated Ringer's solution (Ringer's group, n = 16) or hypertonic saline (Hypertonic Saline group, n = 15) over a 3-day period. Serum ADH levels were significantly larger in the Hypertonic Saline group as compared with the Ringer's group (P = 0.001; analysis of variance) and were correlated to sodium intake (Ringer's group: r = 0.39, R(2) = 0.15, P = 0.02; Hypertonic Saline group: r = 0.42, R(2) = 0.18, P = 0.02) and volume of fluids given IV (Ringer's group: r = 0.38, R(2) = 0.15, P = 0.02; Hypertonic Saline group: r = 0.32, R(2) = 0.1, P = not significant). Correlation of ADH to plasma osmolality was significant if plasma osmolality was >280 mOsm/kg (r = 0.5, R(2) = 0.25, P = 0.06), indicating an osmotic threshold for ADH release. Serum aldosterone levels were larger on the first day than during Days 2 and 3 in both groups and inversely correlated to serum sodium levels only in the Ringer's group (r = -0.55, R(2) = 0.3, P < 0.001). This group received a significantly larger fluid volume on Day 1 (P = 0.05, Mann-Whitney U-test) than did patients in the Hypertonic Saline group, indicating hypovolemia during the first day. Head-injured children have appropriate levels of ADH. They may be hypovolemic during the first day of treatment, especially if they receive lactated Ringer's solution. IMPLICATIONS: In head-injured patients, we recommend fluid restriction to avoid inappropriate secretion of antidiuretic hormone. In a prospective, randomized, and controlled study in 31 children, we were able to show that the antidiuretic hormone levels are appropriate in response to hypovolemia, sodium load, or both.

Gravity-driven slug motion in capillary tubes

The velocity of a liquid slug falling in a capillary tube is lower than predicted for Poiseuille flow due to presence of menisci, whose shapes are determined by the complex interplay of capillary, viscous, and gravitational forces. Due to the presence of menisci, a capillary pressure proportional to surface curvature acts on the slug and streamlines are bent close to the interface, resulting in enhanced viscous dissipation at the wedges. To determine the origin of drag-force increase relative to Poiseuille flow, we compute the force resultant acting on the slug by integrating Navier-Stokes equations over the liquid volume. Invoking relationships from differential geometry we demonstrate that the additional drag is due to viscous forces only and that no capillary drag of hydrodynamic origin exists (i.e., due to hydrodynamic deformation of the interface). Requiring that the force resultant is zero, we derive scaling laws for the steady velocity in the limit of small capillary numbers by estimating the leading order viscous dissipation in the different regions of the slug (i.e., the unperturbed Poiseuille-like bulk, the static menisci close to the tube axis and the dynamic regions close to the contact lines). Considering both partial and complete wetting, we find that the relationship between dimensionless velocity and weight is, in general, nonlinear. Whereas the relationship obtained for complete-wetting conditions is found in agreement with the experimental data of Bico and Quere [J. Bico and D. Quere, J. Colloid Interface Sci. 243, 262 (2001)], the scaling law under partial-wetting conditions is validated by numerical simulations performed with the Volume of Fluid method. The simulated steady velocities agree with the behavior predicted by the theoretical scaling laws in presence and in absence of static contact angle hysteresis. The numerical simulations suggest that wedge-flow dissipation alone cannot account for the entire additional drag and that the non-Poiseuille dissipation in the static menisci (not considered in previous studies) has to be considered for large contact angles.

Post-mortem β-hydroxybutyrate determination in synovial fluid.

β-hydroxybutyrate concentrations were determined in blood and synovial fluid in a series of medico-legal cases including hypothermia fatalities, individuals found dead in a cold environment and non-hypothermia cases with various, non-traumatic causes of death. Hypothermia was considered to be the cause of death according to circumstantial elements indicating exposure to cold, autopsy findings, biochemical investigation results and exclusion of other causes of death. The intention of this study was to characterize β-hydroxybutyrate distribution in synovial fluid and assess its usefulness for the postmortem diagnosis of antemortem abnormalities in blood β-hydroxybutyrate levels. Unenhanced CT scans, autopsies, histology, neuropathology, toxicology, and biochemistry were systematically performed. Within the limited number of subjects included in the study, the results indicate that abnormalities in antemortem β-hydroxybutyrate blood levels, as may be observed in hypothermia fatalities, are reflected in postmortem synovial fluid values. These preliminary findings notwithstanding, synovial fluid analysis to determine β-hydroxybutyrate is unlikely to be generally applied due to the more invasive collection technique it requires and could be limited to special cases in which biological fluids systematically collected upon autopsy are unavailable.

Quantify this! Report on a round table discussion on quantitative mass spectrometry in proteomics.

Following the success of the first round table in 2001, the Swiss Proteomic Society has organized two additional specific events during its last two meetings: a proteomic application exercise in 2002 and a round table in 2003. Such events have as their main objective to bring together, around a challenging topic in mass spectrometry, two groups of specialists, those who develop and commercialize mass spectrometry equipment and software, and expert MS users for peptidomics and proteomics studies. The first round table (Geneva, 2001) entitled "Challenges in Mass Spectrometry" was supported by brief oral presentations that stressed critical questions in the field of MS development or applications (Stöcklin and Binz, Proteomics 2002, 2, 825-827). Topics such as (i) direct analysis of complex biological samples, (ii) status and perspectives for MS investigations of noncovalent peptide-ligant interactions; (iii) is it more appropriate to have complementary instruments rather than a universal equipment, (iv) standardization and improvement of the MS signals for protein identification, (v) what would be the new generation of equipment and finally (vi) how to keep hardware and software adapted to MS up-to-date and accessible to all. For the SPS'02 meeting (Lausanne, 2002), a full session alternative event "Proteomic Application Exercise" was proposed. Two different samples were prepared and sent to the different participants: 100 micro g of snake venom (a complex mixture of peptides and proteins) and 10-20 micro g of almost pure recombinant polypeptide derived from the shrimp Penaeus vannamei carrying an heterogeneous post-translational modification (PTM). Among the 15 participants that received the samples blind, eight returned results and most of them were asked to present their results emphasizing the strategy, the manpower and the instrumentation used during the congress (Binz et. al., Proteomics 2003, 3, 1562-1566). It appeared that for the snake venom extract, the quality of the results was not particularly dependant on the strategy used, as all approaches allowed Lication of identification of a certain number of protein families. The genus of the snake was identified in most cases, but the species was ambiguous. Surprisingly, the precise identification of the recombinant almost pure polypeptides appeared to be much more complicated than expected as only one group reported the full sequence. Finally the SPS'03 meeting reported here included a round table on the difficult and challenging task of "Quantification by Mass Spectrometry", a discussion sustained by four selected oral presentations on the use of stable isotopes, electrospray ionization versus matrix-assisted laser desorption/ionization approaches to quantify peptides and proteins in biological fluids, the handling of differential two-dimensional liquid chromatography tandem mass spectrometry data resulting from high throughput experiments, and the quantitative analysis of PTMs. During these three events at the SPS meetings, the impressive quality and quantity of exchanges between the developers and providers of mass spectrometry equipment and software, expert users and the audience, were a key element for the success of these fruitful events and will have definitively paved the way for future round tables and challenging exercises at SPS meetings.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.