Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Testing differentiation in diploid populations.

We examine the power of different exact tests of differentiation for diploid populations. Since there is not necessarily random mating within populations, the appropriate hypothesis to construct exact tests is that of independent sampling of genotypes. There are two categories of tests, FST-estimator tests and goodness of fit tests. In this latter category, we distinguish "allelic statistics", which account for the nature of alleles within genotypes, from "genotypic statistics" that do not. We show that the power of FST-estimator tests and of allelic goodness of fit tests are similar when sampling is balanced, and higher than the power of genotypic goodness of fit tests. When sampling is unbalanced, the most powerful tests are shown to belong to the allelic goodness of fit group.

Binding free energy differences in a TCR-peptide-MHC complex induced by a peptide mutation: a simulation analysis.

Recognition by the T-cell receptor (TCR) of immunogenic peptides presented by class I major histocompatibility complexes (MHCs) is the determining event in the specific cellular immune response against virus-infected cells or tumor cells. It is of great interest, therefore, to elucidate the molecular principles upon which the selectivity of a TCR is based. These principles can in turn be used to design therapeutic approaches, such as peptide-based immunotherapies of cancer. In this study, free energy simulation methods are used to analyze the binding free energy difference of a particular TCR (A6) for a wild-type peptide (Tax) and a mutant peptide (Tax P6A), both presented in HLA A2. The computed free energy difference is 2.9 kcal/mol, in good agreement with the experimental value. This makes possible the use of the simulation results for obtaining an understanding of the origin of the free energy difference which was not available from the experimental results. A free energy component analysis makes possible the decomposition of the free energy difference between the binding of the wild-type and mutant peptide into its components. Of particular interest is the fact that better solvation of the mutant peptide when bound to the MHC molecule is an important contribution to the greater affinity of the TCR for the latter. The results make possible identification of the residues of the TCR which are important for the selectivity. This provides an understanding of the molecular principles that govern the recognition. The possibility of using free energy simulations in designing peptide derivatives for cancer immunotherapy is briefly discussed.

Diagnostics VIH tardifs en 2009 en Suisse: motivation a un dépistage systématique du VIH [Late diagnoses of HIV infection in 2009 in Switzerland: towards a universal HIV testing?]

Individuals who are unaware of their HIV infection constitute a fragile and a critical population both from a personal and from a community-based point of view: if the infection is diagnosed too late, the outcome can be complicated by AIDS-defining diseases, whose prognosis may remain unfavourable even after the initiation of a potent antiretroviral therapy. These patients contribute to the virus spreading into the community, owing to a high viral load. It is now necessary to recognise the limits of a risk behaviour and disease-driven HIV screening policy. Since 2006, the American guidelines recommend a routine HIV testing for all patients age 13 to 64 years unless they specifically refuse the test (opt-out). In Switzerland, the recommendations remain so far risk-centred.

Preoperative upper gastrointestinal testing can help predicting long-term outcome after gastric banding for morbid obesity.

BACKGROUND: Gastric banding (GB) is one of the most popular bariatric procedures for morbid obesity. Apart from causing weight loss by alimentary restriction, it can interfere with functions of the esophagus and upper stomach. The aim of this study was to evaluate if the results of extensive preoperative upper GI testing were correlated with long-term outcome and complications after GB. METHODS: Using a prospectively maintained computerized database including all the patients undergoing bariatric operations in both our hospitals, we performed a retrospective analysis of the patients who underwent complete upper gastrointestinal (GI) testing (endoscopy, pH monitoring, and manometry) before GB. RESULTS: One hundred thirty-four patients underwent complete testing before GB. Abnormal pH monitoring (increased total reflux time, increased diurnal reflux time, increased number of reflux episodes) predicted the development of complications and especially pouch dilatation and food intolerance. The mean De Meester score was higher among patients who developed complications than in the remaining ones (25.4 vs 17.7, P=0.03). High lower esophageal sphincter pressure also predicted progressive long-term food intolerance. Endoscopic findings were not predictive of the long-term outcome. CONCLUSIONS: There is some association between the function of the upper digestive tract and long-term complications after gastric banding. Abnormal pH monitoring predicts overall long-term complications, especially food intolerance with or without reflux, and pouch dilatation, and a high lower esophageal sphincter pressure predicts long-term food intolerance. Extended upper gastrointestinal testing with endoscopy, 24-h pH monitoring, and esophageal manometry is probably worthwhile in selecting patients for gastric banding.

Five-year monitoring of a gay-friendly voluntary counselling and testing facility in Switzerland: who got tested and why?

ABSTRACT: BACKGROUND: An increase in new HIV cases among men who have sex with men (MSM) has been reported in Switzerland since 2001. A rapid result HIV testing for MSM through voluntary counselling and testing (VCT) facility ("Checkpoint") was opened in Geneva in 2005. This gay-friendly facility, the first to open in Switzerland, provides testing for sexually transmitted infections (STI) and rapid result HIV testing and counselling. Our objective was to analyze Checkpoint's activity over its first five years of activity and its ability to attract at-risk MSM. METHODS: We used routine data collected anonymously about the facility activity (number of clients, number of tests, and test results) and about the characteristics of the clientele (sociodemographic data, sexual risk behaviour, and reasons for testing) from 2005 to 2009. RESULTS: The yearly number of HIV tests performed increased from 249 in 2005 to 561 in 2009. The annual proportion of positive tests among tests performed varied between 2% and 3%. Among MSM clients, the median annual number of anal intercourse (AI) partners was three. Roughly 30% of all MSM clients had at least one unprotected anal intercourse (UAI) experience in the previous 12 months with a partner of different/unknown HIV status.The main reason for testing in 2007, 2008, and 2009 was "sexual risk exposure" (~40%), followed by "routine" testing (~30%) and "condom stopping in the beginning of a new steady relationship" (~10%). Clients who came to the facility after a sexual risk exposure, compared to clients who came for "routine testing" or "condom stopping" reasons, had the highest number of AI partners in the previous 12 months, were more likely to have had UAI with a partner of different/unknown HIV status in the previous 12 months (respectively 57.3%, 12.5%, 23.5%), more likely to have had an STI diagnosed in the past (41.6%, 32.2%, 22.9%), and more likely to report recent feelings of sadness or depression (42.6%; 32.8%, 18.5%). CONCLUSION: Many of Checkpoint's clients reported elevated sexual risk exposure and risk factors, and the annual proportion of new HIV cases in the facility is stable. This VCT facility attracts the intended population and appears to be a useful tool contributing to the fight against the HIV epidemic among MSM in Switzerland.

Simulation of stromatolite growth using diffusion-limited aggregation and a discrete model of biogeochemical exchanges in microbial mat

Resume : Mieux comprendre les stromatolithes et les tapis microbiens est un sujet important en biogéosciences puisque cela aide à l'étude des premières formes de vie sur Terre, a mieux cerner l'écologie des communautés microbiennes et la contribution des microorganismes a la biominéralisation, et même à poser certains fondements dans les recherches en exobiologie. D'autre part, la modélisation est un outil puissant utilisé dans les sciences naturelles pour appréhender différents phénomènes de façon théorique. Les modèles sont généralement construits sur un système d'équations différentielles et les résultats sont obtenus en résolvant ce système. Les logiciels disponibles pour implémenter les modèles incluent les logiciels mathématiques et les logiciels généraux de simulation. L'objectif principal de cette thèse est de développer des modèles et des logiciels pour aider a comprendre, via la simulation, le fonctionnement des stromatolithes et des tapis microbiens. Ces logiciels ont été développés en C++ en ne partant d'aucun pré-requis de façon a privilégier performance et flexibilité maximales. Cette démarche permet de construire des modèles bien plus spécifiques et plus appropriés aux phénomènes a modéliser. Premièrement, nous avons étudié la croissance et la morphologie des stromatolithes. Nous avons construit un modèle tridimensionnel fondé sur l'agrégation par diffusion limitée. Le modèle a été implémenté en deux applications C++: un moteur de simulation capable d'exécuter un batch de simulations et de produire des fichiers de résultats, et un outil de visualisation qui permet d'analyser les résultats en trois dimensions. Après avoir vérifié que ce modèle peut en effet reproduire la croissance et la morphologie de plusieurs types de stromatolithes, nous avons introduit un processus de sédimentation comme facteur externe. Ceci nous a mené a des résultats intéressants, et permis de soutenir l'hypothèse que la morphologie des stromatolithes pourrait être le résultat de facteurs externes autant que de facteurs internes. Ceci est important car la classification des stromatolithes est généralement fondée sur leur morphologie, imposant que la forme d'un stromatolithe est dépendante de facteurs internes uniquement (c'est-à-dire les tapis microbiens). Les résultats avancés dans ce mémoire contredisent donc ces assertions communément admises. Ensuite, nous avons décidé de mener des recherches plus en profondeur sur les aspects fonctionnels des tapis microbiens. Nous avons construit un modèle bidimensionnel de réaction-diffusion fondé sur la simulation discrète. Ce modèle a été implémenté dans une application C++ qui permet de paramétrer et exécuter des simulations. Nous avons ensuite pu comparer les résultats de simulation avec des données du monde réel et vérifier que le modèle peut en effet imiter le comportement de certains tapis microbiens. Ainsi, nous avons pu émettre et vérifier des hypothèses sur le fonctionnement de certains tapis microbiens pour nous aider à mieux en comprendre certains aspects, comme la dynamique des éléments, en particulier le soufre et l'oxygène. En conclusion, ce travail a abouti à l'écriture de logiciels dédiés à la simulation de tapis microbiens d'un point de vue tant morphologique que fonctionnel, suivant deux approches différentes, l'une holistique, l'autre plus analytique. Ces logiciels sont gratuits et diffusés sous licence GPL (General Public License). Abstract : Better understanding of stromatolites and microbial mats is an important topic in biogeosciences as it helps studying the early forms of life on Earth, provides clues re- garding the ecology of microbial ecosystems and their contribution to biomineralization, and gives basis to a new science, exobiology. On the other hand, modelling is a powerful tool used in natural sciences for the theoretical approach of various phenomena. Models are usually built on a system of differential equations and results are obtained by solving that system. Available software to implement models includes mathematical solvers and general simulation software. The main objective of this thesis is to develop models and software able to help to understand the functioning of stromatolites and microbial mats. Software was developed in C++ from scratch for maximum performance and flexibility. This allows to build models much more specific to a phenomenon rather than general software. First, we studied stromatolite growth and morphology. We built a three-dimensional model based on diffusion-limited aggregation. The model was implemented in two C++ applications: a simulator engine, which can run a batch of simulations and produce result files, and a Visualization tool, which allows results to be analysed in three dimensions. After verifying that our model can indeed reproduce the growth and morphology of several types of stromatolites, we introduced a sedimentation process as an external factor. This lead to interesting results, and allowed to emit the hypothesis that stromatolite morphology may be the result of external factors as much as internal factors. This is important as stromatolite classification is usually based on their morphology, imposing that a stromatolite shape is dependant on internal factors only (i.e. the microbial mat). This statement is contradicted by our findings, Second, we decided to investigate deeper the functioning of microbial mats, We built a two-dimensional reaction-diffusion model based on discrete simulation, The model was implemented in a C++ application that allows setting and running simulations. We could then compare simulation results with real world data and verify that our model can indeed mimic the behaviour of some microbial mats. Thus, we have proposed and verified hypotheses regarding microbial mats functioning in order to help to better understand them, e.g. the cycle of some elements such as oxygen or sulfur. ln conclusion, this PhD provides a simulation software, dealing with two different approaches. This software is free and available under a GPL licence.

Pulmonary embolism severity index and troponin testing for the selection of low-risk patients with acute symptomatic pulmonary embolism.

Summary Background: The combination of the Pulmonary Embolism Severity Index (PESI) and troponin testing could help physicians identify appropriate patients with acute pulmonary embolism (PE) for early hospital discharge. Methods: This prospective cohort study included a total of 567 patients from a single center registry with objectively confirmed acute symptomatic PE. On the basis of the PESI, each patient was classified into 1 of 5 classes (I to V). At the time of hospital admission, patients had troponin I (cTnI) levels measured. The endpoint of the study was all-cause mortality within 30 days after diagnosis. We calculated the mortality rates in 4 patient groups: group 1: PESI class I-II plus cTnI <0.1 ng mL(-1); group 2: PESI classes III-V plus cTnI <0.1 ng mL(-1); group 3: PESI classes I-II plus cTnI >/= 0.1 ng mL(-1); and group 4: PESI classes III-V plus cTnI >/= 0.1 ng mL(-1). Results: The study cohort had a 30-day mortality of 10% (95% confidence interval [CI], 7.6 to 12.5%). Mortality rates in the 4 groups were 1.3%, 14.2%, 0% and 15.4%, respectively. Compared to non-elevated cTnl, the low-risk PESI had a higher negative predictive value (NPV) (98.9% vs 90.8%) and negative likelihood ratio (NLR) (0.1 vs 0.9) for predicting mortality. The addition of non-elevated cTnI to low-risk PESI did not improve the NPV or the NLR compared to either test alone. Conclusions: Compared to cTnl testing, PESI classification more accurately identified patients with PE who are at low risk of all-cause death within 30-days of presentation.

Assessing Patient's belief on HIV testing before elective surgery: place for improvements in Switzerland

Background and aim of the study: In Switzerland no HIV test is performed without the patient's consent based on a Voluntary Counseling and Testing policy (VCT). We hypothesized that a substantial proportion of patients going through an elective surgery falsely believed that an HIV test was performed on a routine basis and that the lack of transmission of result was interpreted as being HIV negative. Material and method: All patients with elective orthopedic surgery during 2007 were contacted by phone in 2008. A structured questionnaire assessed their belief about routine preoperative blood analysis (glycemia, coagulation capacity, HIV serology and cholesterol) as well as result awareness and interpretation. Variables included age and gender. Analysis were conducted using the software JMP 6.0.3. Results: 1123 patients were included. 130 (12%) were excluded (i.e. unreachable, unable to communicate on the phone, not operated). 993 completed the survey (89%). Median age was 51 (16-79). 50% were female. 376 (38%) patients thought they had an HIV test performed before surgery but none of them had one. 298 (79%) interpreted the absence of result as a negative HIV test. A predictive factor to believe an HIV test had been done was an age below 50 years old (45% vs 33% for 16-49 years old and 50-79 years old respectively, p <0.001). No difference was observed between genders. Conclusion: In Switzerland, nearly 40% of the patients falsely thought an HIV test had been performed on a routine basis before surgery and were erroneously reassured about their HIV status. These results should either improve the information given to the patient regarding preoperative exams, or motivate public health policy to consider HIV opt-out screening, as patients are already expecting it.