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

Point-of-care serological assays for delayed SARS-CoV-2 case identification among health-care workers in the UK: a prospective multicentre cohort study

BACKGROUND: Health-care workers constitute a high-risk population for acquisition of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Capacity for acute diagnosis via PCR testing was limited for individuals with mild to moderate SARS-CoV-2 infection in the early phase of the COVID-19 pandemic and a substantial proportion of health-care workers with suspected infection were not tested. We aimed to investigate the performance of point-of-care and laboratory serology assays and their utility in late case identification, and to estimate SARS-CoV-2 seroprevalence. METHODS: We did a prospective multicentre cohort study between April 8 and June 12, 2020, in two phases. Symptomatic health-care workers with mild to moderate symptoms were eligible to participate 14 days after onset of COVID-19 symptoms, as per the Public Health England (PHE) case definition. Health-care workers were recruited to the asymptomatic cohort if they had not developed PHE-defined COVID-19 symptoms since Dec 1, 2019. In phase 1, two point-of-care lateral flow serological assays, the Onsite CTK Biotech COVID-19 split IgG/IgM Rapid Test (CTK Bitotech, Poway, CA, USA) and the Encode SARS-CoV-2 split IgM/IgG One Step Rapid Test Device (Zhuhai Encode Medical Engineering, Zhuhai, China), were evaluated for performance against a laboratory immunoassay (EDI Novel Coronavirus COVID-19 IgG ELISA kit [Epitope Diagnostics, San Diego, CA, USA]) in 300 samples from health-care workers and 100 pre-COVID-19 negative control samples. In phase 2 (n=6440), serosurveillance was done among 1299 (93·4%) of 1391 health-care workers reporting symptoms, and in a subset of asymptomatic health-care workers (405 [8·0%] of 5049). FINDINGS: There was variation in test performance between the lateral flow serological assays; however, the Encode assay displayed reasonable IgG sensitivity (127 of 136; 93·4% [95% CI 87·8-96·9]) and specificity (99 of 100; 99·0% [94·6-100·0]) among PCR-proven cases and good agreement (282 of 300; 94·0% [91·3-96·7]) with the laboratory immunoassay. By contrast, the Onsite assay had reduced sensitivity (120 of 136; 88·2% [95% CI 81·6-93·1]) and specificity (94 of 100; 94·0% [87·4-97·8]) and agreement (254 of 300; 84·7% [80·6-88·7]). Five (7%) of 70 PCR-positive cases were negative across all assays. Late changes in lateral flow serological assay bands were recorded in 74 (9·3%) of 800 cassettes (35 [8·8%] of 400 Encode assays; 39 [9·8%] of 400 Onsite assays), but only seven (all Onsite assays) of these changes were concordant with the laboratory immunoassay. In phase 2, seroprevalence among the workforce was estimated to be 10·6% (95% CI 7·6-13·6) in asymptomatic health-care workers and 44·7% (42·0-47·4) in symptomatic health-care workers. Seroprevalence across the entire workforce was estimated at 18·0% (95% CI 17·0-18·9). INTERPRETATION: Although a good positive predictive value was observed with both lateral flow serological assays and ELISA, this agreement only occurred if the pre-test probability was modified by a strict clinical case definition. Late development of lateral flow serological assay bands would preclude postal strategies and potentially home testing. Identification of false-negative results among health-care workers across all assays suggest caution in interpretation of IgG results at this stage; for now, testing is perhaps best delivered in a clinical setting, supported by government advice about physical distancing. FUNDING: None

Structural similarity assessment for drug sensitivity prediction in cancer

<p>Abstract</p> <p>Background</p> <p>The ability to predict drug sensitivity in cancer is one of the exciting promises of pharmacogenomic research. Several groups have demonstrated the ability to predict drug sensitivity by integrating chemo-sensitivity data and associated gene expression measurements from large anti-cancer drug screens such as NCI-60. The general approach is based on comparing gene expression measurements from sensitive and resistant cancer cell lines and deriving drug sensitivity profiles consisting of lists of genes whose expression is predictive of response to a drug. Importantly, it has been shown that such profiles are generic and can be applied to cancer cell lines that are not part of the anti-cancer screen. However, one limitation is that the profiles can not be generated for untested drugs (i.e., drugs that are not part of an anti-cancer drug screen). In this work, we propose using an existing drug sensitivity profile for drug A as a substitute for an untested drug B given high structural similarities between drugs A and B.</p> <p>Results</p> <p>We first show that structural similarity between pairs of compounds in the NCI-60 dataset highly correlates with the similarity between their activities across the cancer cell lines. This result shows that structurally similar drugs can be expected to have a similar effect on cancer cell lines. We next set out to test our hypothesis that we can use existing drug sensitivity profiles as substitute profiles for untested drugs. In a cross-validation experiment, we found that the use of substitute profiles is possible without a significant loss of prediction accuracy if the substitute profile was generated from a compound with high structural similarity to the untested compound.</p> <p>Conclusion</p> <p>Anti-cancer drug screens are a valuable resource for generating omics-based drug sensitivity profiles. We show that it is possible to extend the usefulness of existing screens to untested drugs by deriving substitute sensitivity profiles from structurally similar drugs part of the screen.</p

A lipidomic approach to identify cold-induced changes in Arabidopsis membrane lipid composition

Lipidomic analysis using electrospray ionization triple quadrupole mass spectrometry can be employed to monitor lipid changes that occur during cold and freezing stress of plants. Here we describe the analysis of Arabidopsis thaliana polar glycerolipids with normal and oxidized acyl chains, sampled during cold and freezing treatments. Mass spectral data are processed using the online capabilities of LipidomeDB Data Calculation Environment

Carbonic Anhydrase 5 Regulates Acid-Base Homeostasis in Zebrafish

The regulation of the acid-base balance in cells is essential for proper cellular homeostasis. Disturbed acid-base balance directly affects cellular physiology, which often results in various pathological conditions. In every living organism, the protein family of carbonic anhydrases regulate a broad variety of homeostatic processes. Here we describe the identification, mapping and cloning of a zebrafish carbonic anhydrase 5 (ca5) mutation, collapse of fins (cof), which causes initially a collapse of the medial fins followed by necrosis and rapid degeneration of the embryo. These phenotypical characteristics can be mimicked in wild-type embryos by acetazolamide treatment, suggesting that CA5 activity in zebrafish is essential for a proper development. In addition we show that CA5 regulates acid-base balance during embryonic development, since lowering the pH can compensate for the loss of CA5 activity. Identification of selective modulators of CA5 activity could have a major impact on the development of new therapeutics involved in the treatment of a variety of disorders

RNA polymerase II stalling promotes nucleosome occlusion and pTEFb recruitment to drive immortalization by Epstein-Barr virus

Epstein-Barr virus (EBV) immortalizes resting B-cells and is a key etiologic agent in the development of numerous cancers. The essential EBV-encoded protein EBNA 2 activates the viral C promoter (Cp) producing a message of ~120 kb that is differentially spliced to encode all EBNAs required for immortalization. We have previously shown that EBNA 2-activated transcription is dependent on the activity of the RNA polymerase II (pol II) C-terminal domain (CTD) kinase pTEFb (CDK9/cyclin T1). We now demonstrate that Cp, in contrast to two shorter EBNA 2-activated viral genes (LMP 1 and 2A), displays high levels of promoter-proximally stalled pol II despite being constitutively active. Consistent with pol II stalling, we detect considerable pausing complex (NELF/DSIF) association with Cp. Significantly, we observe substantial Cp-specific pTEFb recruitment that stimulates high-level pol II CTD serine 2 phosphorylation at distal regions (up to +75 kb), promoting elongation. We reveal that Cp-specific pol II accumulation is directed by DNA sequences unfavourable for nucleosome assembly that increase TBP access and pol II recruitment. Stalled pol II then maintains Cp nucleosome depletion. Our data indicate that pTEFb is recruited to Cp by the bromodomain protein Brd4, with polymerase stalling facilitating stable association of pTEFb. The Brd4 inhibitor JQ1 and the pTEFb inhibitors DRB and Flavopiridol significantly reduce Cp, but not LMP1 transcript production indicating that Brd4 and pTEFb are required for Cp transcription. Taken together our data indicate that pol II stalling at Cp promotes transcription of essential immortalizing genes during EBV infection by (i) preventing promoter-proximal nucleosome assembly and ii) necessitating the recruitment of pTEFb thereby maintaining serine 2 CTD phosphorylation at distal regions

Hv-CBF2A overexpression in barley accelerates COR gene transcript accumulation and acquisition of freezing tolerance during cold acclimation

Abstract C-Repeat Binding Factors (CBFs) are DNAbinding transcriptional activators of gene pathways imparting freezing tolerance. Poaceae contain three CBF subfamilies, two of which, HvCBF3/CBFIII and HvCBF4/CBFIV, are unique to this taxon. To gain mechanistic insight into HvCBF4/CBFIV CBFs we overexpressed Hv-CBF2A in spring barley (Hordeum vulgare) cultivar ‘Golden Promise’. The Hv-CBF2A overexpressing lines exhibited stunted growth, poor yield, and greater freezing tolerance compared to non-transformed ‘Golden Promise’. Differences in freezing tolerance were apparent only upon cold acclimation. During cold acclimation freezing tolerance of the Hv-CBF2A overexpressing lines increased more rapidly than that of ‘Golden Promise’ and paralleled the freezing tolerance of the winter hardy barley ‘Dicktoo’. Transcript levels of candidate CBF target genes, COR14B and DHN5 were increased in the overexpressor lines at warm temperatures, and at cold temperatures they accumulated to much higher levels in the Hv-CBF2A overexpressors than in ‘Golden Promise’. Hv-CBF2A overexpression also increased transcript levels of other CBF genes at FROST RESISTANCE-H2-H2 (FR-H2) possessing CRT/DRE sites in their upstream regions, the most notable of which was CBF12. CBF12 transcript levels exhibited a relatively constant incremental increase above levels in ‘Golden Promise’ both at warm and cold. These data indicate that Hv-CBF2A activates target genes at warm temperatures and that transcript accumulation for some of these targets is greatly enhanced by cold temperatures

Expression of the proteolysis-inducing factor core peptide mRNA is upregulated in both tumour and adjacent normal tissue in gastro-oesophageal malignancy

Gastro-oesophageal cancer is associated with a high incidence of cachexia. Proteolysis-inducing factor (PIF) has been identified as a possible cachectic factor and studies suggest that PIF is produced exclusively by tumour cells. We investigated PIF core peptide (PIF-CP) mRNA expression in tumour and benign tissue from patients with gastro-oesophageal cancer and in gastro-oesophageal biopsies for healthy volunteers. Tumour tissue and adjacent benign tissue were collected from patients with gastric and oesophageal cancer (n=46) and from benign tissue only in healthy controls (n=11). Expression of PIF-CP mRNA was quantified by real-time PCR. Clinical and pathological information along with nutritional status was collected prospectively. In the cancer patients, PIF-CP mRNA was detected in 27 (59%) tumour samples and 31 (67%) adjacent benign tissue samples. Four (36%) gastro-oesophageal biopsies from healthy controls also expressed PIF-CP mRNA. Expression was higher in tumour tissue (P=0.031) and benign tissue (P=0.022) from cancer patients compared with healthy controls. In the cancer patients, tumour and adjacent benign tissue PIF-CP mRNA concentrations were correlated with each other (P<0.0001, r=0.73) but did not correlate with weight loss or prognosis. Although PIF-CP mRNA expression is upregulated in both tumour and adjacent normal tissue in gastro-oesophageal malignancy, expression does not relate to prognosis or cachexia. Post-translational modification of PIF may be a key step in determining the biological role of PIF in the patient with advanced cancer and cachexia

Capsular profiling of the Cronobacter genus and the association of specific Cronobacter sakazakii and C. malonaticus capsule types with neonatal meningitis and necrotizing enterocolitis

Background: Cronobacter sakazakii and C. malonaticus can cause serious diseases especially in infants where they are associated with rare but fatal neonatal infections such as meningitis and necrotising enterocolitis. Methods: This study used 104 whole genome sequenced strains, covering all seven species in the genus, to analyse capsule associated clusters of genes involved in the biosynthesis of the O-antigen, colanic acid, bacterial cellulose, enterobacterial common antigen (ECA), and a previously uncharacterised K-antigen. Results: Phylogeny of the gnd and galF genes flanking the O-antigen region enabled the defining of 38 subgroups which are potential serotypes. Two variants of the colanic acid synthesis gene cluster (CA1 and CA2) were found which differed with the absence of galE in CA2. Cellulose (bcs genes) were present in all species, but were absent in C. sakazakii sequence type (ST) 13 and clonal complex (CC) 100 strains. The ECA locus was found in all strains. The K-antigen capsular polysaccharide Region 1 (kpsEDCS) and Region 3 (kpsMT) genes were found in all Cronobacter strains. The highly variable Region 2 genes were assigned to 2 homology groups (K1 and K2). C. sakazakii and C. malonaticus isolates with capsular type [K2:CA2:Cell+] were associated with neonatal meningitis and necrotizing enterocolitis. Other capsular types were less associated with clinical infections. Conclusion: This study proposes a new capsular typing scheme which identifies a possible important virulence trait associated with severe neonatal infections. The various capsular polysaccharide structures warrant further investigation as they could be relevant to macrophage survival, desiccation resistance, environmental survival, and biofilm formation in the hospital environment, including neonatal enteral feeding tubes