Using simulation studies to evaluate statistical methods

Simulation studies are computer experiments that involve creating data by pseudorandom sampling. The key strength of simulation studies is the ability to understand the behaviour of statistical methods because some 'truth' (usually some parameter/s of interest) is known from the process of generating the data. This allows us to consider properties of methods, such as bias. While widely used, simulation studies are often poorly designed, analysed and reported. This tutorial outlines the rationale for using simulation studies and offers guidance for design, execution, analysis, reporting and presentation. In particular, this tutorial provides: a structured approach for planning and reporting simulation studies, which involves defining aims, data-generating mechanisms, estimands, methods and performance measures ('ADEMP'); coherent terminology for simulation studies; guidance on coding simulation studies; a critical discussion of key performance measures and their estimation; guidance on structuring tabular and graphical presentation of results; and new graphical presentations. With a view to describing recent practice, we review 100 articles taken from Volume 34 of Statistics in Medicine that included at least one simulation study and identify areas for improvement.Comment: 31 pages, 9 figures (2 in appendix), 8 tables (1 in appendix

Bim and Bmf synergize to induce apoptosis in Neisseria gonorrhoeae infection

Abstract: Bcl-2 family proteins including the pro-apoptotic BH3-only proteins are central regulators of apoptotic cell death. Here we show by a focused siRNA miniscreen that the synergistic action of the BH3-only proteins Bim and Bmf is required for apoptosis induced by infection with Neisseria gonorrhoeae (Ngo). While Bim and Bmf were associated with the cytoskeleton of healthy cells, they both were released upon Ngo infection. Loss of Bim and Bmf from the cytoskeleton fraction required the activation of Jun-N-terminal kinase-1 (JNK-1), which in turn depended on Rac-1. Depletion and inhibition of Rac-1, JNK-1, Bim, or Bmf prevented the activation of Bak and Bax and the subsequent activation of caspases. Apoptosis could be reconstituted in Bim-depleted and Bmf-depleted cells by additional silencing of antiapoptotic Mcl-1 and Bcl-XL, respectively. Our data indicate a synergistic role for both cytoskeletal-associated BH3-only proteins, Bim, and Bmf, in an apoptotic pathway leading to the clearance of Ngo-infected cells. Author Summary: A variety of physiological death signals, as well as pathological insults, trigger apoptosis, a genetically programmed form of cell death. Pathogens often induce host cell apoptosis to establish a successful infection. Neisseria gonorrhoeae (Ngo), the etiological agent of the sexually transmitted disease gonorrhoea, is a highly adapted obligate human-specific pathogen and has been shown to induce apoptosis in infected cells. Here we unveil the molecular mechanisms leading to apoptosis of infected cells. We show that Ngo-mediated apoptosis requires a special subset of proapoptotic proteins from the group of BH3-only proteins. BH3-only proteins act as stress sensors to translate toxic environmental signals to the initiation of apoptosis. In a siRNA-based miniscreen, we found Bim and Bmf, BH3-only proteins associated with the cytoskeleton, necessary to induce host cell apoptosis upon infection. Bim and Bmf inactivated different inhibitors of apoptosis and thereby induced cell death in response to infection. Our data unveil a novel pathway of infection-induced apoptosis that enhances our understanding of the mechanism by which BH3-only proteins control apoptotic cell death

The Adaptor Molecule Nck Localizes the WAVE Complex to Promote Actin Polymerization during CEACAM3-Mediated Phagocytosis of Bacteria

Background: CEACAM3 is a granulocyte receptor mediating the opsonin-independent recognition and phagocytosis of human-restricted CEACAM-binding bacteria. CEACAM3 function depends on an intracellular immunoreceptor tyrosine-based activation motif (ITAM)-like sequence that is tyrosine phosphorylated by Src family kinases upon receptor engagement. The phosphorylated ITAM-like sequence triggers GTP-loading of Rac by directly associating with the guanine nucleotide exchange factor (GEF) Vav. Rac stimulation in turn is critical for actin cytoskeleton rearrangements that generate lamellipodial protrusions and lead to bacterial uptake. Principal Findings: In our present study we provide biochemical and microscopic evidence that the adaptor proteins Nck1 and Nck2, but not CrkL, Grb2 or SLP-76, bind to tyrosine phosphorylated CEACAM3. The association is phosphorylation-dependent and requires the Nck SH2 domain. Overexpression of the isolated Nck1 SH2 domain, RNAi-mediated knock-down of Nck1, or genetic deletion of Nck1 and Nck2 interfere with CEACAM3-mediated bacterial internalization and with the formation of lamellipodial protrusions. Nck is constitutively associated with WAVE2 and directs the actin nucleation promoting WAVE complex to tyrosine phosphorylated CEACAM3. In turn, dominant-negative WAVE2 as well as shRNA-mediated knock-down of WAVE2 or the WAVE-complex component Nap1 reduce internalization of bacteria. Conclusions: Our results provide novel mechanistic insight into CEACAM3-initiated phagocytosis. We suggest that the CEACAM3 ITAM-like sequence is optimized to co-ordinate a minimal set of cellular factors needed to efficiently trigger actin-based lamellipodial protrusions and rapid pathogen engulfment

Factors associated with infant mortality in Nepal: a comparative analysis of Nepal demographic and health surveys (NDHS) 2006 and 2011

Background: Infant mortality is one of the priority public health issues in developing countries like Nepal. The infant mortality rate (IMR) was 48 and 46 per 1000 live births for the year 2006 and 2011, respectively, a slight reduction during the 5 years’ period. A comprehensive analysis that has identified and compared key factors associated with infant mortality is limited in Nepal, and, therefore, this study aims to fill the gap. Methods: Datasets from Nepal Demographic and Health Surveys (NDHS) 2006 and 2011 were used to identify and compare the major factors associated with infant mortality. Both surveys used multistage stratified cluster sampling techniques. A total of 8707 and 10,826 households were interviewed in 2006 and 2011, with more than 99% response rate in both studies. The survival information of singleton live-born infants born 5 years preceding the two surveys were extracted from the ‘childbirth’ dataset. Multiple logistic regression analysis using a hierarchical modelling approach with the backward elimination method was conducted. Complex Samples Analysis was used to adjust for unequal selection probability due to the multistage stratified cluster-sampling procedure used in both NDHS.Results: Based on NDHS 2006, ecological region, succeeding birth interval, breastfeeding status and type of delivery assistance were found to be significant predictors of infant mortality. Infants born in hilly region (AOR = 0.43, p = 0.013) and with professional assistance (AOR = 0.27, p = 0.039) had a lower risk of mortality. On the other hand, infants with succeeding birth interval less than 24 months (AOR = 6.66, p = 0.001) and those who were never breastfed (AOR = 1.62, p = 0.044) had a higher risk of mortality. Based on NDHS 2011, birth interval (preceding and succeeding) and baby’s size at birth were identified to be significantly associated with infant mortality. Infants born with preceding birth interval (AOR = 1.94, p = 0.022) or succeeding birth interval (AOR = 3.22, p = 0.002) shorter than 24 months had higher odds of mortality while those born with a very large or larger than average size had significantly lowered odds (AOR = 0.17, p = 0.008) of mortality. Conclusion: IMR and associated risk factors differ between NDHS 2006 and 2011 except ‘succeeding birth interval’ which attained significant status in the both study periods. This study identified the ecological region, birth interval, delivery assistant type, baby’s birth size and breastfeeding status as significant predictors of infant mortality

Assessment of the In Vivo Toxicity of Gold Nanoparticles

The environmental impact of nanoparticles is evident; however, their toxicity due to their nanosize is rarely discussed. Gold nanoparticles (GNPs) may serve as a promising model to address the size-dependent biological response to nanoparticles because they show good biocompatibility and their size can be controlled with great precision during their chemical synthesis. Naked GNPs ranging from 3 to 100 nm were injected intraperitoneally into BALB/C mice at a dose of 8 mg/kg/week. GNPs of 3, 5, 50, and 100 nm did not show harmful effects; however, GNPs ranging from 8 to 37 nm induced severe sickness in mice. Mice injected with GNPs in this range showed fatigue, loss of appetite, change of fur color, and weight loss. Starting from day 14, mice in this group exhibited a camel-like back and crooked spine. The majority of mice in these groups died within 21 days. Injection of 5 and 3 nm GNPs, however, did not induce sickness or lethality in mice. Pathological examination of the major organs of the mice in the diseased groups indicated an increase of Kupffer cells in the liver, loss of structural integrity in the lungs, and diffusion of white pulp in the spleen. The pathological abnormality was associated with the presence of gold particles at the diseased sites, which were verified by ex vivo Coherent anti-Stoke Raman scattering microscopy. Modifying the surface of the GNPs by incorporating immunogenic peptides ameliorated their toxicity. This reduction in the toxicity is associated with an increase in the ability to induce antibody response. The toxicity of GNPs may be a fundamental determinant of the environmental toxicity of nanoparticles

Tetanus toxin Hc fragment induces the formation of ceramide platforms and protects neuronal cells against oxidative stress

Tetanus toxin (TeTx) is the protein, synthesized by the anaerobic bacteria Clostridium tetani, which causes tetanus disease. TeTx gains entry into target cells by means of its interaction with lipid rafts, which are membrane domains enriched in sphingomyelin and cholesterol. However, the exact mechanism of host membrane binding remains to be fully established. In the present study we used the recombinant carboxyl terminal fragment from TeTx (Hc-TeTx), the domain responsible for target neuron binding, showing that Hc-TeTx induces a moderate but rapid and sustained increase in the ceramide/sphingomyelin ratio in primary cultures of cerebellar granule neurons and in NGF-differentiated PC12 cells, as well as induces the formation of ceramide platforms in the plasma membrane. The mentioned increase is due to the promotion of neutral sphingomyelinase activity and not to the de novo synthesis, since GW4869, a specific neutral sphingomyelinase inhibitor, prevents neutral sphingomyelinase activity increase and formation of ceramide platforms. Moreover, neutral sphingomyelinase inhibition with GW4869 prevents Hc-TeTx-triggered signaling (Akt phosphorylation), as well as the protective effect of Hc-TeTx on PC12 cells subjected to oxidative stress, while siRNA directed against nSM2 prevents protection by Hc-TeTx of NSC-34 cells against oxidative insult. Finally, neutral sphingomyelinase activity seems not to be related with the internalization of Hc-TeTx into PC12 cells. Thus, the presented data shed light on the mechanisms triggered by TeTx after membrane binding, which could be related with the events leading to the neuroprotective action exerted by the Hc-TeTx fragment

Allosteric Regulation of Fibronectin/α5β1 Interaction by Fibronectin-Binding MSCRAMMs

Citation: Liang, X. W., Garcia, B. L., Visai, L., Prabhakaran, S., Meenan, N. A. G., Potts, J. R., . . . Hook, M. (2016). Allosteric Regulation of Fibronectin/alpha(5)beta(1) Interaction by Fibronectin-Binding MSCRAMMs. Plos One, 11(7), 17. doi:10.1371/journal.pone.0159118Adherence ofmicrobes to host tissues is a hallmark of infectious disease and is often mediated by a class of adhesins termed MSCRAMMs (Microbial Surface Components Recognizing Adhesive Matrix Molecules). Numerous pathogens express MSCRAMMs that specifically bind the heterodimeric human glycoprotein fibronectin (Fn). In addition to roles in adhesion, Fn-binding MSCRAMMs exploit physiological Fn functions. For example, several pathogens can invade host cells by a mechanism whereby MSCRAMM-bound Fn bridges interaction with alpha(5)beta(1) integrin. Here, we investigate two Fn-binding MSCRAMMs, FnBPA (Staphylococcus aureus) and BBK32 (Borrelia burgdorferi) to probe structure-activity relationships of MSCRAMM-induced Fn/alpha(5)beta(1) integrin activation. Circular dichroism, fluorescence resonance energy transfer, and dynamic light scattering techniques uncover a conformational rearrangement of Fn involving domains distant from the MSCRAMM binding site. Surface plasmon resonance experiments demonstrate a significant enhancement of Fn/alpha(5)beta(1) integrin affinity in the presence of FnBPA or BBK32. Detailed kinetic analysis of these interactions reveal that this change in affinity can be attributed solely to an increase in the initial Fn/alpha(5)beta(1) on-rate and that this rate-enhancement is dependent on high-affinity Fn-binding by MSCRAMMs. These data implicate MSCRAMM-induced perturbation of specific intramolecular contacts within the Fn heterodimer resulting in activation by exposing previously cryptic alpha(5)beta(1) interaction motifs. By correlating structural changes in Fn to a direct measurement of increased Fn/alpha(5)beta(1) affinity, this work significantly advances our understanding of the structural basis for the modulation of integrin function by Fn-binding MSCRAMMs

The Na(+)–H(+ )exchanger-1 induces cytoskeletal changes involving reciprocal RhoA and Rac1 signaling, resulting in motility and invasion in MDA-MB-435 cells

INTRODUCTION: An increasing body of evidence shows that the tumour microenvironment is essential in driving neoplastic progression. The low serum component of this microenvironment stimulates motility/invasion in human breast cancer cells via activation of the Na(+)–H(+ )exchanger (NHE) isoform 1, but the signal transduction systems that underlie this process are still poorly understood. We undertook the present study to elucidate the role and pattern of regulation by the Rho GTPases of this serum deprivation-dependent activation of both NHE1 and subsequent invasive characteristics, such as pseudopodia and invadiopodia protrusion, directed cell motility and penetration of normal tissues. METHODS: The present study was performed in a well characterized human mammary epithelial cell line representing late stage metastatic progression, MDA-MB-435. The activity of RhoA and Rac1 was modified using their dominant negative and constitutively active mutants and the activity of NHE1, cell motility/invasion, F-actin content and cell shape were measured. RESULTS: We show for the first time that serum deprivation induces NHE1-dependent morphological and cytoskeletal changes in metastatic cells via a reciprocal interaction of RhoA and Rac1, resulting in increased chemotaxis and invasion. Deprivation changed cell shape by reducing the amount of F-actin and inducing the formation of leading edge pseudopodia. Serum deprivation inhibited RhoA activity and stimulated Rac1 activity. Rac1 and RhoA were antagonistic regulators of both basal and stimulated tumour cell NHE1 activity. The regulation of NHE1 activity by RhoA and Rac1 in both conditions was mediated by an alteration in intracellular proton affinity of the exchanger. Interestingly, the role of each of these G-proteins was reversed during serum deprivation; basal NHE1 activity was regulated positively by RhoA and negatively by Rac1, whereas RhoA negatively and Rac1 positively directed the stimulation of NHE1 during serum deprivation. Importantly, the same pattern of RhoA and Rac1 regulation found for NHE1 activity was observed in both basal and serum deprivation dependent increases in motility, invasion and actin cytoskeletal organization. CONCLUSION: Our findings suggest that the reported antagonistic roles of RhoA and Rac1 in cell motility/invasion and cytoskeletal organization may be due, in part, to their concerted action on NHE1 activity as a convergence point

Novel Association Strategy with Copy Number Variation for Identifying New Risk Loci of Human Diseases

Copy number variations (CNV) are important causal genetic variations for human disease; however, the lack of a statistical model has impeded the systematic testing of CNVs associated with disease in large-scale cohort.Here, we developed a novel integrated strategy to test CNV-association in genome-wide case-control studies. We converted the single-nucleotide polymorphism (SNP) signal to copy number states using a well-trained hidden Markov model. We mapped the susceptible CNV-loci through SNP site-specific testing to cope with the physiological complexity of CNVs. We also ensured the credibility of the associated CNVs through further window-based CNV-pattern clustering. Genome-wide data with seven diseases were used to test our strategy and, in total, we identified 36 new susceptible loci that are associated with CNVs for the seven diseases: 5 with bipolar disorder, 4 with coronary artery disease, 1 with Crohn's disease, 7 with hypertension, 9 with rheumatoid arthritis, 7 with type 1 diabetes and 3 with type 2 diabetes. Fifteen of these identified loci were validated through genotype-association and physiological function from previous studies, which provide further confidence for our results. Notably, the genes associated with bipolar disorder converged in the phosphoinositide/calcium signaling, a well-known affected pathway in bipolar disorder, which further supports that CNVs have impact on bipolar disorder.Our results demonstrated the effectiveness and robustness of our CNV-association analysis and provided an alternative avenue for discovering new associated loci of human diseases