Symbolic Formulae for Linear Mixed Models

A statistical model is a mathematical representation of an often simplified or idealised data-generating process. In this paper, we focus on a particular type of statistical model, called linear mixed models (LMMs), that is widely used in many disciplines e.g.~agriculture, ecology, econometrics, psychology. Mixed models, also commonly known as multi-level, nested, hierarchical or panel data models, incorporate a combination of fixed and random effects, with LMMs being a special case. The inclusion of random effects in particular gives LMMs considerable flexibility in accounting for many types of complex correlated structures often found in data. This flexibility, however, has given rise to a number of ways by which an end-user can specify the precise form of the LMM that they wish to fit in statistical software. In this paper, we review the software design for specification of the LMM (and its special case, the linear model), focusing in particular on the use of high-level symbolic model formulae and two popular but contrasting R-packages in lme4 and asreml

Taxonomic variations in the gut microbiome of gout patients with and without tophi might have a functional impact on urate metabolism

Objective: To evaluate the taxonomic composition of the gut microbiome in gout patients with and without tophi formation, and predict bacterial functions that might have an impact on urate metabolism. Methods: Hypervariable V3–V4 regions of the bacterial 16S rRNA gene from fecal samples of gout patients with and without tophi (n=33 and n=25, respectively) were sequenced and compared to fecal samples from 53 healthy controls. We explored predictive functional profles using bioinformatics in order to identify diferences in taxonomy and metabolic pathways. Results: We identifed a microbiome characterized by the lowest richness and a higher abundance of Phascolarctobacterium, Bacteroides, Akkermansia, and Ruminococcus_gnavus_group genera in patients with gout without tophi when compared to controls. The Proteobacteria phylum and the Escherichia-Shigella genus were more abundant in patients with tophaceous gout than in controls. Fold change analysis detected nine genera enriched in healthy controls compared to gout groups (Bifdobacterium, Butyricicoccus, Oscillobacter, Ruminococcaceae_UCG_010, Lachnospiraceae_ND2007_group, Haemophilus, Ruminococcus_1, Clostridium_sensu_stricto_1, and Ruminococcaceae_ UGC_013). We found that the core microbiota of both gout groups shared Bacteroides caccae, Bacteroides stercoris ATCC 43183, and Bacteroides coprocola DSM 17136. These bacteria might perform functions linked to one-carbon metabo‑ lism, nucleotide binding, amino acid biosynthesis, and purine biosynthesis. Finally, we observed diferences in key bacterial enzymes involved in urate synthesis, degradation, and elimination. Conclusion: Our fndings revealed that taxonomic variations in the gut microbiome of gout patients with and with‑ out tophi might have a functional impact on urate metabolism. Keywords: Gout, Gut microbiota, Uric acid metabolis

The Endogenous Th17 Response in NO<inf>2</inf>-Promoted Allergic Airway Disease Is Dispensable for Airway Hyperresponsiveness and Distinct from Th17 Adoptive Transfer

Severe, glucocorticoid-resistant asthma comprises 5-7% of patients with asthma. IL-17 is a biomarker of severe asthma, and the adoptive transfer of Th17 cells in mice is sufficient to induce glucocorticoid-resistant allergic airway disease. Nitrogen dioxide (NO2) is an environmental toxin that correlates with asthma severity, exacerbation, and risk of adverse outcomes. Mice that are allergically sensitized to the antigen ovalbumin by exposure to NO2 exhibit a mixed Th2/Th17 adaptive immune response and eosinophil and neutrophil recruitment to the airway following antigen challenge, a phenotype reminiscent of severe clinical asthma. Because IL-1 receptor (IL-1R) signaling is critical in the generation of the Th17 response in vivo, we hypothesized that the IL-1R/Th17 axis contributes to pulmonary inflammation and airway hyperresponsiveness (AHR) in NO2-promoted allergic airway disease and manifests in glucocorticoid-resistant cytokine production. IL-17A neutralization at the time of antigen challenge or genetic deficiency in IL-1R resulted in decreased neutrophil recruitment to the airway following antigen challenge but did not protect against the development of AHR. Instead, IL-1R-/- mice developed exacerbated AHR compared to WT mice. Lung cells from NO2-allergically inflamed mice that were treated in vitro with dexamethasone (Dex) during antigen restimulation exhibited reduced Th17 cytokine production, whereas Th17 cytokine production by lung cells from recipient mice of in vitro Th17-polarized OTII T-cells was resistant to Dex. These results demonstrate that the IL-1R/Th17 axis does not contribute to AHR development in NO2-promoted allergic airway disease, that Th17 adoptive transfer does not necessarily reflect an endogenously-generated Th17 response, and that functions of Th17 responses are contingent on the experimental conditions in which they are generated. © 2013 Martin et al

Serogroup C meningococci in Italy in the era of conjugate menC vaccination

<p>Abstract</p> <p>Background</p> <p>To assess changes in the pattern of Invasive Meningococcal Disease (IMD) in Italy after the introduction of conjugate menC vaccine in the National Vaccine Plan 2005–2007 and to provide information for developing timely and appropriate public health interventions, analyses of microbiological features of isolates and clinical characteristics of patients have been carried out. In Italy, the number of serogroup C meningococci fell progressively following the introduction of the MenC conjugate vaccine, recommended by the Italian Ministry of Health but implemented according to different regional strategies.</p> <p>Methods</p> <p>IMD cases from January 2005 through July 2008 reported to the National Meningococcal Surveillance System were considered for this study. Serogrouping and sero/subtyping were performed on 179 serogroup C strains received at the National Reference Laboratory of the Istituto Superiore di Sanità. Antibiotic susceptibility testing was possible for 157 isolates. MLST (Multilocus sequence typing), <it>por</it>A VRs (Variable Region) typing, PFGE (Pulsed Field Gel Electrophoresis), VNTR (Variable Number Tandem Repeats) analyses were performed on all C:2a and C:2b meningococci (n = 147), following standard procedures.</p> <p>Results</p> <p>In 2005 and 2008, IMD showed an incidence of 0.5 and 0.3 per 100,000 inhabitants, respectively. While the incidence due to serogroup B remained stable, IMD incidence due to serogroup C has decreased since 2006. In particular, the decrease was significant among infants. C:2a and C:2b were the main serotypes, all C:2a strains belonged to ST-11 clonal complex and all C:2b to ST-8/A4. Clinical manifestations and outcome of infections underlined more severe disease caused by C:2a isolates. Two clusters due to C:2a/ST-11 meningococci were reported in the North of Italy in December 2007 and July 2008, respectively, with a high rate of septicaemia and fatal outcome.</p> <p>Conclusion</p> <p>Public health surveillance of serogroup C invasive meningococcal disease and microbiological/molecular characterization of the isolates requires particular attention, since the hyper-invasive ST-11 predominantly affected adolescents and young adults for whom meningococcal vaccination was not recommended in the 2005–2007 National Vaccine Plan.</p

Drug utilization and cost in a Medicaid population: A simulation study of community vs. mail order pharmacy

<p>Abstract</p> <p>Background</p> <p>Outpatient drugs are dispensed through both community and mail order pharmacies. There is no empirical evidence that substitution of community pharmacy with mail order reduces overall drug expenditures. The need for evaluating the potential effects on utilization and costs of the possible extension of mail order services in Medicaid provides the rationale for conducting this study. This study compares drug utilization and drug product cost in community vs. mail order pharmacy dispensing services in a Medicaid population.</p> <p>Methods</p> <p>This study is a retrospective cohort study comparing utilization and cost patterns in community vs. mail order pharmacy. A simulation model was employed to assess drug utilization and cost in mail order pharmacy using community pharmacy claim data. The model assumed that courses of drug therapy (CDT) in mail order pharmacy would have utilization patterns similar to those found in community pharmacy. A 95% confidence interval surrounding changes in average utilization and average cost were estimated using bootstrap analysis. A sensitivity analysis was performed by varying drug selection criteria and supply, fill point, and medication possession ratio (MPR). Sub-analyses were performed to address differences between mail order and community pharmacy related to therapeutic class and dual-eligible patients.</p> <p>Data for the study derived from pharmacy claims database of Ohio Medicaid State program for the period January 2000-September 2004. Drug claims were aggregated to obtain a set of CDTs representing unique patient IDs and unique drug products. Drug product cost estimates excluded dispensing fees and were used to estimate the cost reduction required in mail order to become cost neutral in comparison with community pharmacy.</p> <p>Results</p> <p>The baseline model revealed that the use of mail order vs. community pharmacy would result in a 5.5% increase in drug utilization and a 5.4% cost reduction required in mail order to become cost neutral. Results from Ohio Medicaid drugs for chronic use revealed a 5.1% increase in utilization and a 4.9% cost reduction required to become cost neutral in comparison with community pharmacy.</p> <p>Conclusion</p> <p>The results of the simulation model indicate that mail order pharmacy increases drug utilization and can also increase drug product cost if the cost per unit is not reduced accordingly. Prior consideration should be given to the patient population, day-supply, disease, therapy, and insurance characteristics to ensure the appropriate use of mail order pharmacy services.</p

Assessment of the Antibiofilm Performance of Chitosan-Based Surfaces in Marine Environments

18 pages, 6 figures.-- This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) licenseMarine biofouling is a natural process often associated with biofilm formation on submerged surfaces, creating a massive economic and ecological burden. Although several antifouling paints have been used to prevent biofouling, growing ecological concerns emphasize the need to develop new and environmentally friendly antifouling approaches such as bio-based coatings. Chitosan (CS) is a natural polymer that has been widely used due to its outstanding biological properties, including non-toxicity and antimicrobial activity. This work aims to produce and characterize poly (lactic acid) (PLA)-CS surfaces with CS of different molecular weight (Mw) at different concentrations for application in marine paints. Loligo opalescens pens, a waste from the fishery industry, were used as a CS source. The antimicrobial activity of the CS and CS-functionalized surfaces was assessed against Cobetia marina, a model proteobacterium for marine biofouling. Results demonstrate that CS targets the bacterial cell membrane, and PLA-CS surfaces were able to reduce the number of culturable cells up to 68% compared to control, with this activity dependent on CS Mw. The antifouling performance was corroborated by Optical Coherence Tomography since PLA-CS surfaces reduced the biofilm thickness by up to 36%, as well as the percentage and size of biofilm empty spaces. Overall, CS coatings showed to be a promising approach to reducing biofouling in marine environments mimicked in this work, contributing to the valorization of fishing waste and encouraging further research on this topic.This work was financially supported by: LA/P/0045/2020 (ALiCE), UIDB/00511/2020 and UIDP/00511/2020 (LEPABE) were funded by national funds through the FCT/MCTES (PIDDAC); project HealthyWaters (NORTE-01-0145-FEDER-000069) was supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF); CVMAR + I—Industrial Innovation and Marine Biotechnology Valorization” project was funded by INTERREG V Espanha Portugal (POCTEP) (0302_CVMAR_I_1_P). The research was also supported by the SurfSAFE project funded by the European Union’s Horizon 2020 research and innovation programme under grant agreement No. 952471. L.C.G. and A.I.B. thank the Portuguese Foundation for Science and Technology (FCT) for the financial support of their work contracts through the Scientific Employment Stimulus— Individual Call (ref. CEECIND/01700/2017 and 2020.03447.CEECIND, respectively). R.T.-S. acknowledges the receipt of a junior researcher fellowship from the Project PTDC/CTM-COM/4844/2020 (NanoCAT) supported by national funds through the FCT/MCTES (PIDDAC). M.L. and M.J.R. also acknowledge FCT for their PhD grants (2022.11196.BD and SFRH/BD/140080/2018, respectively). J.A.V. and J.V. also thank Xunta de Galicia by Xunta de Galicia (Grupos de Potencial Crecimiento, IN607B 2018/2019) for the financial support. Support from the EURO-MIC COST Action (CA20130) is acknowledgePeer reviewe

Regional differences in mitochondrial DNA methylation in human post-mortem brain tissue

Background: DNA methylation is an important epigenetic mechanism involved in gene regulation, with alterations in DNA methylation in the nuclear genome being linked to numerous complex diseases. Mitochondrial DNA methylation is a phenomenon that is receiving ever-increasing interest, particularly in diseases characterized by mitochondrial dysfunction; however, most studies have been limited to the investigation of specific target regions. Analyses spanning the entire mitochondrial genome have been limited, potentially due to the amount of input DNA required. Further, mitochondrial genetic studies have been previously confounded by nuclear-mitochondrial pseudogenes. Methylated DNA Immunoprecipitation Sequencing is a technique widely used to profile DNA methylation across the nuclear genome; however, reads mapped to mitochondrial DNA are often discarded. Here, we have developed an approach to control for nuclear-mitochondrial pseudogenes within Methylated DNA Immunoprecipitation Sequencing data. We highlight the utility of this approach in identifying differences in mitochondrial DNA methylation across regions of the human brain and pre-mortem blood. Results: We were able to correlate mitochondrial DNA methylation patterns between the cortex, cerebellum and blood. We identified 74 nominally significant differentially methylated regions (p < 0.05) in the mitochondrial genome, between anatomically separate cortical regions and the cerebellum in matched samples (N = 3 matched donors). Further analysis identified eight significant differentially methylated regions between the total cortex and cerebellum after correcting for multiple testing. Using unsupervised hierarchical clustering analysis of the mitochondrial DNA methylome, we were able to identify tissue-specific patterns of mitochondrial DNA methylation between blood, cerebellum and cortex. Conclusions: Our study represents a comprehensive analysis of the mitochondrial methylome using pre-existing Methylated DNA Immunoprecipitation Sequencing data to identify brain region-specific patterns of mitochondrial DNA methylation

Nuclear translocation of haeme oxygenase-1 is associated to prostate cancer

The role of oxidative stress in prostate cancer has been increasingly recognised. Acute and chronic inflammations generate reactive oxygen species that result in damage to cellular structures. Haeme oxygenase-1 (HO-1) has cytoprotective effects against oxidative damage. We hypothesise that modulation of HO-1 expression may be involved in the process of prostate carcinogenesis and prostate cancer progression. We thus studied HO-1 expression and localisation in 85 samples of organ-confined primary prostate cancer obtained via radical prostatectomy (Gleason grades 4–9) and in 39 specimens of benign prostatic hyperplasia (BPH). We assessed HO-1 expression by immunohistochemical staining. No significant difference was observed in the cytoplasmic positive reactivity among tumours (84%), non-neoplastic surrounding parenchyma (89%), or BPH samples (87%) (P=0.53). Haeme oxygenase-1 immunostaining was detected in the nuclei of prostate cancer cells in 55 of 85 (65%) patients but less often in non-neoplastic surrounding parenchyma (30 of 85, 35%) or in BPH (9 of 39, 23%) (P<0.0001). Immunocytochemical and western blot analysis showed HO-1 only in the cytoplasmic compartment of PC3 and LNCaP prostate cancer cell lines. Treatment with hemin, a well-known specific inducer of HO-1, led to clear nuclear localisation of HO-1 in both cell lines and highly induced HO-1 expression in both cellular compartments. These findings have demonstrated, for the first time, that HO-1 expression and nuclear localisation can define a new subgroup of prostate cancer primary tumours and that the modulation of HO-1 expression and its nuclear translocation could represent new avenues for therapy

Observation of associated near-side and away-side long-range correlations in √sNN=5.02 TeV proton-lead collisions with the ATLAS detector

Two-particle correlations in relative azimuthal angle (Δϕ) and pseudorapidity (Δη) are measured in √sNN=5.02  TeV p+Pb collisions using the ATLAS detector at the LHC. The measurements are performed using approximately 1  μb-1 of data as a function of transverse momentum (pT) and the transverse energy (ΣETPb) summed over 3.1<η<4.9 in the direction of the Pb beam. The correlation function, constructed from charged particles, exhibits a long-range (2<|Δη|<5) “near-side” (Δϕ∼0) correlation that grows rapidly with increasing ΣETPb. A long-range “away-side” (Δϕ∼π) correlation, obtained by subtracting the expected contributions from recoiling dijets and other sources estimated using events with small ΣETPb, is found to match the near-side correlation in magnitude, shape (in Δη and Δϕ) and ΣETPb dependence. The resultant Δϕ correlation is approximately symmetric about π/2, and is consistent with a dominant cos⁡2Δϕ modulation for all ΣETPb ranges and particle pT

Measurement of the cross-section of high transverse momentum vector bosons reconstructed as single jets and studies of jet substructure in pp collisions at √s = 7 TeV with the ATLAS detector

This paper presents a measurement of the cross-section for high transverse momentum W and Z bosons produced in pp collisions and decaying to all-hadronic final states. The data used in the analysis were recorded by the ATLAS detector at the CERN Large Hadron Collider at a centre-of-mass energy of √s = 7 TeV;{\rm Te}{\rm V}$and correspond to an integrated luminosity of$4.6\;{\rm f}{{{\rm b}}^{-1}}$. The measurement is performed by reconstructing the boosted W or Z bosons in single jets. The reconstructed jet mass is used to identify the W and Z bosons, and a jet substructure method based on energy cluster information in the jet centre-of-mass frame is used to suppress the large multi-jet background. The cross-section for events with a hadronically decaying W or Z boson, with transverse momentum${{p}_{{\rm T}}}\gt 320\;{\rm Ge}{\rm V}$and pseudorapidity$|\eta |\lt 1.9$, is measured to be${{\sigma }_{W+Z}}=8.5\pm 1.7$ pb and is compared to next-to-leading-order calculations. The selected events are further used to study jet grooming techniques