Discrete and continuous mathematical models of sharp-fronted collective cell migration and invasion

Mathematical models describing the spatial spreading and invasion of populations of biological cells are often developed in a continuum modelling framework using reaction-diffusion equations. While continuum models based on linear diffusion are routinely employed and known to capture key experimental observations, linear diffusion fails to predict well-defined sharp fronts that are often observed experimentally. This observation has motivated the use of nonlinear degenerate diffusion, however these nonlinear models and the associated parameters lack a clear biological motivation and interpretation. Here we take a different approach by developing a stochastic discrete lattice-based model incorporating biologically-inspired mechanisms and then deriving the reaction-diffusion continuum limit. Inspired by experimental observations, agents in the simulation deposit extracellular material, that we call a substrate, locally onto the lattice, and the motility of agents is taken to be proportional to the substrate density. Discrete simulations that mimic a two--dimensional circular barrier assay illustrate how the discrete model supports both smooth and sharp-fronted density profiles depending on the rate of substrate deposition. Coarse-graining the discrete model leads to a novel partial differential equation (PDE) model whose solution accurately approximates averaged data from the discrete model. The new discrete model and PDE approximation provides a simple, biologically motivated framework for modelling the spreading, growth and invasion of cell populations with well-defined sharp frontsComment: 47 Pages, 8 Figure

Housing Arrangement and Location Determine the Likelihood of Housing Loss Due to Wildfire

Surging wildfires across the globe are contributing to escalating residential losses and have major social, economic, and ecological consequences. The highest losses in the U.S. occur in southern California, where nearly 1000 homes per year have been destroyed by wildfires since 2000. Wildfire risk reduction efforts focus primarily on fuel reduction and, to a lesser degree, on house characteristics and homeowner responsibility. However, the extent to which land use planning could alleviate wildfire risk has been largely missing from the debate despite large numbers of homes being placed in the most hazardous parts of the landscape. Our goal was to examine how housing location and arrangement affects the likelihood that a home will be lost when a wildfire occurs. We developed an extensive geographic dataset of structure locations, including more than 5500 structures that were destroyed or damaged by wildfire since 2001, and identified the main contributors to property loss in two extensive, fire-prone regions in southern California. The arrangement and location of structures strongly affected their susceptibility to wildfire, with property loss most likely at low to intermediate structure densities and in areas with a history of frequent fire. Rates of structure loss were higher when structures were surrounded by wildland vegetation, but were generally higher in herbaceous fuel types than in higher fuel-volume woody types. Empirically based maps developed using housing pattern and location performed better in distinguishing hazardous from non-hazardous areas than maps based on fuel distribution. The strong importance of housing arrangement and location indicate that land use planning may be a critical tool for reducing fire risk, but it will require reliable delineations of the most hazardous locations

Identifying uncertainties in Arctic climate change projections

Wide ranging climate changes are expected in the Arctic by the end of the 21st century, but projections of the size of these changes vary widely across current global climate models. This variation represents a large source of uncertainty in our understanding of the evolution of Arctic climate. Here we systematically quantify and assess the model uncertainty in Arctic climate changes in two CO2 doubling experiments: a multimodel ensemble (CMIP3) and an ensemble constructed using a single model (HadCM3) with multiple parameter perturbations (THC-QUMP). These two ensembles allow us to assess the contribution that both structural and parameter variations across models make to the total uncertainty and to begin to attribute sources of uncertainty in projected changes. We find that parameter uncertainty is an major source of uncertainty in certain aspects of Arctic climate. But also that uncertainties in the mean climate state in the 20th century, most notably in the northward Atlantic ocean heat transport and Arctic sea ice volume, are a significant source of uncertainty for projections of future Arctic change. We suggest that better observational constraints on these quantities will lead to significant improvements in the precision of projections of future Arctic climate change

Scalable production of large quantities of defect-free few-layer graphene by shear exfoliation in liquids

To progress from the laboratory to commercial applications, it will be necessary to develop industrially scalable methods to produce large quantities of defect-free graphene. Here we show that high-shear mixing of graphite in suitable stabilizing liquids results in large-scale exfoliation to give dispersions of graphene nanosheets. X-ray photoelectron spectroscopy and Raman spectroscopy show the exfoliated flakes to be unoxidized and free of basal-plane defects. We have developed a simple model that shows exfoliation to occur once the local shear rate exceeds 10(4) s(-1). By fully characterizing the scaling behaviour of the graphene production rate, we show that exfoliation can be achieved in liquid volumes from hundreds of millilitres up to hundreds of litres and beyond. The graphene produced by this method performs well in applications from composites to conductive coatings. This method can be applied to exfoliate BN, MoS2 and a range of other layered crystals

The impact of viral mutations on recognition by SARS-CoV-2 specific T cells.

We identify amino acid variants within dominant SARS-CoV-2 T cell epitopes by interrogating global sequence data. Several variants within nucleocapsid and ORF3a epitopes have arisen independently in multiple lineages and result in loss of recognition by epitope-specific T cells assessed by IFN-γ and cytotoxic killing assays. Complete loss of T cell responsiveness was seen due to Q213K in the A∗01:01-restricted CD8+ ORF3a epitope FTSDYYQLY207-215; due to P13L, P13S, and P13T in the B∗27:05-restricted CD8+ nucleocapsid epitope QRNAPRITF9-17; and due to T362I and P365S in the A∗03:01/A∗11:01-restricted CD8+ nucleocapsid epitope KTFPPTEPK361-369. CD8+ T cell lines unable to recognize variant epitopes have diverse T cell receptor repertoires. These data demonstrate the potential for T cell evasion and highlight the need for ongoing surveillance for variants capable of escaping T cell as well as humoral immunity.This work is supported by the UK Medical Research Council (MRC); Chinese Academy of Medical Sciences(CAMS) Innovation Fund for Medical Sciences (CIFMS), China; National Institute for Health Research (NIHR)Oxford Biomedical Research Centre, and UK Researchand Innovation (UKRI)/NIHR through the UK Coro-navirus Immunology Consortium (UK-CIC). Sequencing of SARS-CoV-2 samples and collation of data wasundertaken by the COG-UK CONSORTIUM. COG-UK is supported by funding from the Medical ResearchCouncil (MRC) part of UK Research & Innovation (UKRI),the National Institute of Health Research (NIHR),and Genome Research Limited, operating as the Wellcome Sanger Institute. T.I.d.S. is supported by a Well-come Trust Intermediate Clinical Fellowship (110058/Z/15/Z). L.T. is supported by the Wellcome Trust(grant number 205228/Z/16/Z) and by theUniversity of Liverpool Centre for Excellence in Infectious DiseaseResearch (CEIDR). S.D. is funded by an NIHR GlobalResearch Professorship (NIHR300791). L.T. and S.C.M.are also supported by the U.S. Food and Drug Administration Medical Countermeasures Initiative contract75F40120C00085 and the National Institute for Health Research Health Protection Research Unit (HPRU) inEmerging and Zoonotic Infections (NIHR200907) at University of Liverpool inpartnership with Public HealthEngland (PHE), in collaboration with Liverpool School of Tropical Medicine and the University of Oxford.L.T. is based at the University of Liverpool. M.D.P. is funded by the NIHR Sheffield Biomedical ResearchCentre (BRC – IS-BRC-1215-20017). ISARIC4C is supported by the MRC (grant no MC_PC_19059). J.C.K.is a Wellcome Investigator (WT204969/Z/16/Z) and supported by NIHR Oxford Biomedical Research Centreand CIFMS. The views expressed are those of the authors and not necessarily those of the NIHR or MRC

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio

Basic science232. Certolizumab pegol prevents pro-inflammatory alterations in endothelial cell function

Background: Cardiovascular disease is a major comorbidity of rheumatoid arthritis (RA) and a leading cause of death. Chronic systemic inflammation involving tumour necrosis factor alpha (TNF) could contribute to endothelial activation and atherogenesis. A number of anti-TNF therapies are in current use for the treatment of RA, including certolizumab pegol (CZP), (Cimzia ®; UCB, Belgium). Anti-TNF therapy has been associated with reduced clinical cardiovascular disease risk and ameliorated vascular function in RA patients. However, the specific effects of TNF inhibitors on endothelial cell function are largely unknown. Our aim was to investigate the mechanisms underpinning CZP effects on TNF-activated human endothelial cells. Methods: Human aortic endothelial cells (HAoECs) were cultured in vitro and exposed to a) TNF alone, b) TNF plus CZP, or c) neither agent. Microarray analysis was used to examine the transcriptional profile of cells treated for 6 hrs and quantitative polymerase chain reaction (qPCR) analysed gene expression at 1, 3, 6 and 24 hrs. NF-κB localization and IκB degradation were investigated using immunocytochemistry, high content analysis and western blotting. Flow cytometry was conducted to detect microparticle release from HAoECs. Results: Transcriptional profiling revealed that while TNF alone had strong effects on endothelial gene expression, TNF and CZP in combination produced a global gene expression pattern similar to untreated control. The two most highly up-regulated genes in response to TNF treatment were adhesion molecules E-selectin and VCAM-1 (q 0.2 compared to control; p > 0.05 compared to TNF alone). The NF-κB pathway was confirmed as a downstream target of TNF-induced HAoEC activation, via nuclear translocation of NF-κB and degradation of IκB, effects which were abolished by treatment with CZP. In addition, flow cytometry detected an increased production of endothelial microparticles in TNF-activated HAoECs, which was prevented by treatment with CZP. Conclusions: We have found at a cellular level that a clinically available TNF inhibitor, CZP reduces the expression of adhesion molecule expression, and prevents TNF-induced activation of the NF-κB pathway. Furthermore, CZP prevents the production of microparticles by activated endothelial cells. This could be central to the prevention of inflammatory environments underlying these conditions and measurement of microparticles has potential as a novel prognostic marker for future cardiovascular events in this patient group. Disclosure statement: Y.A. received a research grant from UCB. I.B. received a research grant from UCB. S.H. received a research grant from UCB. All other authors have declared no conflicts of interes