ABCtoolbox: a versatile toolkit for approximate Bayesian computations.

BACKGROUND: The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. RESULTS: Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. CONCLUSION: ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Two-dimensional probabilistic inversion of plane-wave electromagnetic data: Methodology, model constraints and joint inversion with electrical resistivity data

Probabilistic inversion methods based on Markov chain Monte Carlo (MCMC) simulation are well suited to quantify parameter and model uncertainty of nonlinear inverse problems. Yet, application of such methods to CPU-intensive forward models can be a daunting task, particularly if the parameter space is high dimensional. Here, we present a 2-D pixel-based MCMC inversion of plane-wave electromagnetic (EM) data. Using synthetic data, we investigate how model parameter uncertainty depends on model structure constraints using different norms of the likelihood function and the model constraints, and study the added benefits of joint inversion of EM and electrical resistivity tomography (ERT) data. Our results demonstrate that model structure constraints are necessary to stabilize the MCMC inversion results of a highly discretized model. These constraints decrease model parameter uncertainty and facilitate model interpretation. A drawback is that these constraints may lead to posterior distributions that do not fully include the true underlying model, because some of its features exhibit a low sensitivity to the EM data, and hence are difficult to resolve. This problem can be partly mitigated if the plane-wave EM data is augmented with ERT observations. The hierarchical Bayesian inverse formulation introduced and used herein is able to successfully recover the probabilistic properties of the measurement data errors and a model regularization weight. Application of the proposed inversion methodology to field data from an aquifer demonstrates that the posterior mean model realization is very similar to that derived from a deterministic inversion with similar model constraints.

Heritability of renal function in hypertensive families of African descent in the Seychelles (Indian Ocean).

BACKGROUND: We estimated the heritability of three measures of glomerular filtration rate (GFR) in hypertensive families of African descent in the Seychelles (Indian Ocean). METHODS: Families with at least two hypertensive siblings and an average of two normotensive siblings were identified through a national hypertension register. Using the ASSOC program in SAGE (Statistical Analysis in Genetic Epidemiology), the age- and gender-adjusted narrow sense heritability of GFR was estimated by maximum likelihood assuming multivariate normality after power transformation. ASSOC can calculate the additive polygenic component of the variance of a trait from pedigree data in the presence of other familial correlations. The effects of body mass index (BMI), blood pressure, natriuresis, along with sodium to potassium ratio in urine and diabetes, were also tested as covariates. RESULTS: Inulin clearance, 24-hour creatinine clearance, and GFR based on the Cockcroft-Gault formula were available for 348 persons from 66 pedigrees. The age- and gender-adjusted correlations (+/- SE) were 0.51 (+/- 0.04) between inulin clearance and creatinine clearance, 0.53 (+/- 0.04) between inulin clearance and Cockcroft-Gault formula and 0.66 (+/- 0.03) between creatinine clearance and Cockcroft-Gault formula. The age- and gender-adjusted heritabilities (+/- SE) of GFR were 0.41 (+/- 0.10) for inulin clearance, 0.52 (+/- 0.13) for creatinine clearance, and 0.82 (+/- 0.09) for Cockcroft-Gault formula. Adjustment for BMI slightly lowered the correlations and heritabilities for all measurements whereas adjustment for blood pressure had virtually no effect. CONCLUSION: The significant heritability estimates of GFR in our sample of families of African descent confirm the familial aggregation of this trait and justify further analyses aimed at discovering genetic determinants of GFR.

Adaptively weighted maximum likelihood estimation of discrete distributions

SummaryDiscrete data arise in various research fields, typically when the observations are count data.I propose a robust and efficient parametric procedure for estimation of discrete distributions. The estimation is done in two phases. First, a very robust, but possibly inefficient, estimate of the model parameters is computed and used to indentify outliers. Then the outliers are either removed from the sample or given low weights, and a weighted maximum likelihood estimate (WML) is computed.The weights are determined via an adaptive process such that if the data follow the model, then asymptotically no observation is downweighted.I prove that the final estimator inherits the breakdown point of the initial one, and that its influence function at the model is the same as the influence function of the maximum likelihood estimator, which strongly suggests that it is asymptotically fully efficient.The initial estimator is a minimum disparity estimator (MDE). MDEs can be shown to have full asymptotic efficiency, and some MDEs have very high breakdown points and very low bias under contamination. Several initial estimators are considered, and the performances of the WMLs based on each of them are studied.It results that in a great variety of situations the WML substantially improves the initial estimator, both in terms of finite sample mean square error and in terms of bias under contamination. Besides, the performances of the WML are rather stable under a change of the MDE even if the MDEs have very different behaviors.Two examples of application of the WML to real data are considered. In both of them, the necessity for a robust estimator is clear: the maximum likelihood estimator is badly corrupted by the presence of a few outliers.This procedure is particularly natural in the discrete distribution setting, but could be extended to the continuous case, for which a possible procedure is sketched.RésuméLes données discrètes sont présentes dans différents domaines de recherche, en particulier lorsque les observations sont des comptages.Je propose une méthode paramétrique robuste et efficace pour l'estimation de distributions discrètes. L'estimation est faite en deux phases. Tout d'abord, un estimateur très robuste des paramètres du modèle est calculé, et utilisé pour la détection des données aberrantes (outliers). Cet estimateur n'est pas nécessairement efficace. Ensuite, soit les outliers sont retirés de l'échantillon, soit des faibles poids leur sont attribués, et un estimateur du maximum de vraisemblance pondéré (WML) est calculé.Les poids sont déterminés via un processus adaptif, tel qu'asymptotiquement, si les données suivent le modèle, aucune observation n'est dépondérée.Je prouve que le point de rupture de l'estimateur final est au moins aussi élevé que celui de l'estimateur initial, et que sa fonction d'influence au modèle est la même que celle du maximum de vraisemblance, ce qui suggère que cet estimateur est pleinement efficace asymptotiquement.L'estimateur initial est un estimateur de disparité minimale (MDE). Les MDE sont asymptotiquement pleinement efficaces, et certains d'entre eux ont un point de rupture très élevé et un très faible biais sous contamination. J'étudie les performances du WML basé sur différents MDEs.Le résultat est que dans une grande variété de situations le WML améliore largement les performances de l'estimateur initial, autant en terme du carré moyen de l'erreur que du biais sous contamination. De plus, les performances du WML restent assez stables lorsqu'on change l'estimateur initial, même si les différents MDEs ont des comportements très différents.Je considère deux exemples d'application du WML à des données réelles, où la nécessité d'un estimateur robuste est manifeste : l'estimateur du maximum de vraisemblance est fortement corrompu par la présence de quelques outliers.La méthode proposée est particulièrement naturelle dans le cadre des distributions discrètes, mais pourrait être étendue au cas continu.

Better-quasi-order : ideals and spaces

This thesis deals with combinatorics, order theory and descriptive set theory. The first contribution is to the theory of well-quasi-orders (wqo) and better-quasi-orders (bqo). The main result is the proof of a conjecture made by Maurice Pouzet in 1978 his thèse d'état which states that any wqo whose ideal completion remainder is bqo is actually bqo. Our proof relies on new results with both a combinatorial and a topological flavour concerning maps from a front into a compact metric space. The second contribution is of a more applied nature and deals with topological spaces. We define a quasi-order on the subsets of every second countable To topological space in a way that generalises the Wadge quasi-order on the Baire space, while extending its nice properties to virtually all these topological spaces. The Wadge quasi-order of reducibility by continuous functions is wqo on Borei subsets of the Baire space, this quasi-order is however far less satisfactory for other important topological spaces such as the real line, as Hertling, Ikegami and Schlicht notably observed. Some authors have therefore studied reducibility with respect to some classes of discontinuous functions to remedy this situation. We propose instead to keep continuity but to weaken the notion of function to that of relation. Using the notion of admissible representation studied in Type-2 theory of effectivity, we define the quasi-order of reducibility by relatively continuous relations. We show that this quasi-order both refines the classical hierarchies of complexity and is wqo on the Borei subsets of virtually every second countable To space - including every (quasi-)Polish space. -- Cette thèse se situe dans les domaines de la combinatoire, de la théorie des ordres et de la théorie descriptive. La première contribution concerne la théorie des bons quasi-ordres (wqo) et des meilleurs quasi-ordres (bqo). Le résultat principal est la preuve d'une conjecture, énoncée par Pouzet en 1978 dans sa thèse d'état, qui établit que tout wqo dont l'ensemble des idéaux non principaux ordonnés par inclusion forme un bqo est alors lui-même un bqo. La preuve repose sur de nouveaux résultats, qui allient la combinatoire et la topologie, au sujet des fonctions d'un front vers un espace métrique compact. La seconde contribution de cette thèse traite de la complexité topologique dans le cadre des espaces To à base dénombrable. Dans le cas de l'espace de Baire, le quasi-ordre de Wadge est un wqo sur les sous-ensembles Boréliens qui a suscité énormément d'intérêt. Cependant cette relation de réduction par fonctions continues s'avère bien moins satisfaisante pour d'autres espaces d'importance tels que la droite réelle, comme l'ont fait notamment remarquer Hertling, Schlicht et Ikegami. Nous proposons de conserver la continuité et d'affaiblir la notion de fonction pour celle de relation. Pour ce faire, nous utilisons la notion de représentation admissible étudiée en « Type-2 theory of effectivity » initiée par Weihrauch. Nous introduisons alors le quasi-ordre de réduction par relations relativement continues et montrons que celui-ci à la fois raffine les hiérarchies classiques de complexité topologique et forme un wqo sur les sous-ensembles Boréliens de chaque espace quasi-Polonais.

Effect on renal function of tenofovir co-administered with either efavirenz or boosted lopinavir or boosted atazanavir: the Swiss HIV Cohort Study

Reduced re'nal function has been reported with tenofovir disoproxil fumarate (TDF). It is not clear whether TDF co-administered with a boosted protease inhibitor (PI) leads to a greater decline in renal function than TDF co-administered with a non-nucleoside reverse transcriptase inhibitor (NNRTI).Methods: We selected ail antiretroviral therapy-naive patients in the Swiss HIV Cohort Study (SHCS) with calibrated or corrected serum creatinine measurements starting antiretroviral therapy with TDF and either efavirenz (EFV) or the ritonavir-boosted PIs, lopinavir (LPV/r) or atazanavir (ATV/r). As a measure of renal function, we used the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation to estimate the glomerular filtration rate (eGFR). We calculated the difference in eGFR over time between two therapies using a marginal model for repeated measures. In weighted analyses, observations were weighted by the product of their point of treatment and censoring weights to adjust for differences both in the sort of patients starting each therapy and in the sort of patients remaining on each therapy over time.Results: By March 2011, 940 patients with at least one creatinine measurement on a first therapy with either TDF and EFV (n=484), TDF and LPVlr (n=269) or TDF and ATV/r (n=187) had been followed for a median of 1. 7, 1.2 and 1.3 years, respectively. Table 1 shows the difference in average estimated GFR (eGFR) over time since starting cART for two marginal models. The first model was not adjusted for potential confounders; the second mode! used weights to adjust for confounders. The results suggest a greater decline in renal function during the first 6 months if TDF is used with a PI rather than with an NNRTI, but no further difference between these therapies after the first 6 months. TDF and ATV/r may lead to a greater decline in the first 6 months than TDF and LPVlr.Conclusions: TDF co-administered with a boosted PI leads to a greater de cline in renal function over the first 6 months of therapy than TDF co-administered with an NNRTI; this decline may be worse with ATV/r than with LPV/r.

Heritability, determinants and reference values of renal length: a family-based population study.

OBJECTIVES: In this population-based study, reference values were generated for renal length, and the heritability and factors associated with kidney length were assessed. METHODS: Anthropometric parameters and renal ultrasound measurements were assessed in randomly selected nuclear families of European ancestry (Switzerland). The adjusted narrow sense heritability of kidney size parameters was estimated by maximum likelihood assuming multivariate normality after power transformation. Gender-specific reference centiles were generated for renal length according to body height in the subset of non-diabetic non-obese participants with normal renal function. RESULTS: We included 374 men and 419 women (mean ± SD, age 47 ± 18 and 48 ± 17 years, BMI 26.2 ± 4 and 24.5 ± 5 kg/m(2), respectively) from 205 families. Renal length was 11.4 ± 0.8 cm in men and 10.7 ± 0.8 cm in women; there was no difference between right and left renal length. Body height, weight and estimated glomerular filtration rate (eGFR) were positively associated with renal length, kidney function negatively, age quadratically, whereas gender and hypertension were not. The adjusted heritability estimates of renal length and volume were 47.3 ± 8.5 % and 45.5 ± 8.8 %, respectively (P < 0.001). CONCLUSION: The significant heritability of renal length and volume highlights the familial aggregation of this trait, independently of age and body size. Population-based references for renal length provide a useful guide for clinicians. KEY POINTS: • Renal length and volume are heritable traits, independent of age and size. • Based on a European population, gender-specific reference values/percentiles are provided for renal length. • Renal length correlates positively with body length and weight. • There was no difference between right and left renal lengths in this study. • This negates general teaching that the left kidney is larger and longer.

Altered NKG2D function in NK cells induced by chronic exposure to NKG2D ligand-expressing tumor cells.

NKG2D is an activation receptor that allows natural killer (NK) cells to detect diseased host cells. The engagement of NKG2D with corresponding ligand results in surface modulation of the receptor and reduced function upon subsequent receptor engagement. However, it is not clear whether in addition to modulation the NKG2D receptor complex and/or its signaling capacity is preserved. We show here that the prolonged encounter with tumor cell-bound, but not soluble, ligand can completely uncouple the NKG2D receptor from the intracellular mobilization of calcium and the exertion of cell-mediated cytolysis. However, cytolytic effector function is intact since NKG2D ligand-exposed NK cells can be activated via the Ly49D receptor. While NKG2D-dependent cytotoxicity is impaired, prolonged ligand exposure results in constitutive interferon gamma (IFNgamma) production, suggesting sustained signaling. The functional changes are associated with a reduced presence of the relevant signal transducing adaptors DNAX-activating protein of 10 kDa (DAP-10) and killer cell activating receptor-associated protein/DNAX-activating protein of 12 kDa (KARAP/DAP-12). That is likely the consequence of constitutive NKG2D engagement and signaling, since NKG2D function and adaptor expression is restored to normal when the stimulating tumor cells are removed. Thus, the chronic exposure to tumor cells expressing NKG2D ligand alters NKG2D signaling and may facilitate the evasion of tumor cells from NK cell reactions.

Differential effects of two anti-apo-cytochrome c-specific monoclonal antibodies on the function of apo-cytochrome c-specific murine T cell clones.

A murine monoclonal antibody (SJL 2-4) specific for the antigen apo-cytochrome c was shown to inhibit both antigen-induced proliferation and lymphokine secretion by an apo-cytochrome c-specific BALB/c helper T cell clone. The inhibition was specific because additional apo-cytochrome c-specific T cell clones were not inhibited by the same monoclonal antibody. Time course studies of the inhibition indicated that the initial 8 hr of contact between T cell clones and antigen-presenting cells were critical for activation of the T cell clones. Inhibition of T cell functions by antigen-specific antibodies appeared to correlate with the antibody-antigen binding constant because a second monoclonal antibody (Cyt-1-59), with identical specificity but with a lower affinity constant for apo-cytochrome c, had very little inhibitory effect on the proliferation or lymphokine secretion of apo-cytochrome c-specific T cell clones.

Structure-function relationships of the alpha1b-adrenergic receptor.

The alpha1b-adrenergic receptor (AR) is a member of the large superfamily of seven transmembrane domain (TMD) G protein-coupled receptors (GPCR). Combining site-directed mutagenesis of the alpha1b-AR with computational simulations of receptor dynamics, we have explored the conformational changes underlying the process of receptor activation, i.e. the transition between the inactive and active states. Our findings suggest that the structural constraint stabilizing the alpha1b-AR in the inactive form is a network of H-bonding interactions amongst conserved residues forming a polar pocket and R143 of the DRY sequence at the end of TMDIII. We have recently reported that point mutations of D142, of the DRY sequence and of A293 in the distal portion of the third intracellular loop resulted in ligand-independent (constitutive) activation of the alpha1b-AR. These constitutively activating mutations could induce perturbations resulting in the shift of R143 out of the polar pocket. The main role of R143 may be to mediate receptor activation by triggering the exposure of several basic amino acids of the intracellular loops towards the G protein. Our investigation has been extended also to the biochemical events involved in the desensitization process of alpha1b-AR. Our results indicate that immediately following agonist-induced activation, the alpha1b-AR can undergo rapid agonist-induced phosphorylation and desensitization. Different members of the G protein coupled receptor kinase family can play a role in agonist-induced regulation of the alpha1b-AR. In addition, constitutively active alpha1b-AR mutants display different phosphorylation and internalization features. The future goal is to further elucidate the molecular mechanism underlying the complex equilibrium between activation and inactivation of the alpha1b-AR and its regulation by pharmacological substances. These findings can help to elucidate the mechanism of action of various agents displaying properties of agonists or inverse agonists at the adrenergic system.

Connexin-36 contributes to control function of insulin-producing cells.

Connexin-36 (Cx36) is a gap junction protein expressed by the insulin-producing beta-cells. We investigated the contribution of this protein in normal beta-cell function by using a viral gene transfer approach to alter Cx36 content in the insulin-producing line of INS-1E cells and rat pancreatic islets. Transcripts for Cx43, Cx45, and Cx36 were detected by reverse transcriptase-PCR in freshly isolated pancreatic islets, whereas only a transcript for Cx36 was detected in INS-1E cells. After infection with a sense viral vector, which induced de novo Cx36 expression in the Cx-defective HeLa cells we used to control the transgene expression, Western blot, immunofluorescence, and freeze-fracture analysis showed a large increase of Cx36 within INS-1E cell membranes. In contrast, after infection with an antisense vector, Cx36 content was decreased by 80%. Glucose-induced insulin release and insulin content were decreased, whether infected INS-1E cells expressed Cx36 levels that were largely higher or lower than those observed in wild-type control cells. In both cases, basal insulin secretion was unaffected. Comparable observations on basal secretion and insulin content were made in freshly isolated rat pancreatic islets. The data indicate that large changes in Cx36 alter insulin content and, at least in INS-1E cells, also affect glucose-induced insulin release.

In vitro engineering of human stratified urothelium: analysis of its morphology and function.

PURPOSE: Gastric or intestinal patches, commonly used for reconstructive cystoplasty, may induce severe metabolic complications. The use of bladder tissues reconstructed in vitro could avoid these complications. We compared cellular differentiation and permeability characteristics of human native with in vitro cultured stratified urothelium. MATERIALS AND METHODS: Human stratified urothelium was induced in vitro. Morphology was studied with light and electron microscopy and expression of key cellular proteins was assessed using immunohistochemistry. Permeability coefficients were determined by measuring water, urea, ammonia and proton fluxes across the urothelium. RESULTS: As in native urothelium the stratified urothelial construct consisted of basal membrane and basal, intermediate and superficial cell layers. The apical membrane of superficial cells formed villi and glycocalices, and tight junctions and desmosomes were developed. Immunohistochemistry showed similarities and differences in the expression of cytokeratins, integrin and cellular adhesion proteins. In the cultured urothelium cytokeratin 20 and integrin subunits alpha6 and beta4 were absent, and symplekin was expressed diffusely in all layers. Uroplakins were clearly expressed in the superficial umbrella cells of the urothelial constructs, however, they were also present in intermediate and basal cells. Symplekin and uroplakins were expressed only in the superficial cells of native bladder tissue. The urothelial constructs showed excellent viability, and functionally their permeabilities for water, urea and ammonia were no different from those measured in native human urothelium. Proton permeability was even lower in the constructs compared to that of native urothelium. CONCLUSIONS: Although the in vitro cultured human stratified urothelium did not show complete terminal differentiation of its superficial cells, it retained the same barrier characteristics against the principal urine components. These results indicate that such in vitro cultured urothelium, after being grown on a compliant degradable support or in coculture with smooth muscle cells, is suitable for reconstructive cystoplasty.

Immune regulation and function of the NOD-like receptors NLRC5 and NLRP3

AbstractThe vertebrate immune system is composed of the innate and the adaptive branches. Innate immune cells represent the first line of defense and detect pathogens through pattern recognition receptors (PRRs), detecting evolutionary conserved pathogen- and danger- associated molecular patterns. Engagement of these receptors initiates the inflammatory response, but also instructs antigen-specific adaptive immune cells. NOD-like receptors (NLRs) are an important group of PRRs, leading to the production of inflammatory mediators and favoring antigen presentation to Τ lymphocytes through the regulation of major histocompatibility complex (MHC) molecules.In this work we focused our attention on selected NOD-like receptors (NLRs) and their role at the interface between innate and adaptive immunity. First, we describe a new regulatory mechanism controlling IL-1 production. Our results indicate that type I interferons (IFNs) block NLRP1 and NLRP3 inflammasome activity and interfere with LPS-driven proIL-Ια and -β induction. As type I IFNs are produced upon viral infections, these anti-inflammatory effects of type I IFN could be relevant in the context of superinfections, but could also help explaining the efficacy of IFN-β in multiple sclerosis treatment.The second project addresses the role of a novel NLR family member, called NLRC5. The function of this NLR is still matter of debate, as it has been proposed as both an inhibitor and an activator of different inflammatory pathways. We found that the expression of this protein is restricted to immune cells and is positively regulated by IFNs. We generated Nlrc5-deficient mice and found that this NLR plays an essential role in Τ, NKT and, NK lymphocytes, in which it drives the expression of MHC class I molecules. Accordingly, we could show that CD8+ Τ cell-mediated killing of target lymphocytes lacking NLRC5 is strongly impaired. Moreover, NLRC5 expression was found to be low in many lymphoid- derived tumor cell lines, a mechanism that could be exploited by tumors to escape immunosurveillance.Finally, we found NLRC5 to be involved in the production of IL-10 by CD4+ Τ cells, as Nlrc5- deficient Τ lymphocytes produced less of this cytokine upon TCR triggering. In line with these observations, Mrc5-deficient CD4+ Τ cells expanded more than control cells when transferred into lymphopenic hosts and led to a more rapid appearance of colitis symptoms. Therefore, our work gives novel insights on the function of NLRC5 by using knockout mice, and strongly supports the idea that NLRs direct not only innate, but also adaptive immune responses.