Bile Duct Injuries During Urgent Cholecystectomy at a Safety Net Teaching Hospital: Attending Experience and Time of Day May Matter.

Background: Bile duct injury (BDI) is one of the most severe complications during cholecystectomy. Early identification of risk factors for BDI may permit risk reduction strategies and inform patient consent.Objective: This study aimed to define patient, provider, and systemic factors associated with BDI; BDI incidence; and short-term outcomes of BDI after urgent cholecystectomy.Methods: Patients who underwent urgent cholecystectomy for acute cholecystitis were retrospectively screened (2020-2022). All patients who sustained BDI were included without exclusions. Demographics, clinical data, and outcomes were collected and compared with descriptive statistics.Results: During the study period, BDI occurred in 4 (0.5%) of 728 patients who underwent urgent cholecystectomy for acute cholecystitis. Most BDI cases (75%) took place overnight or during the weekend. The attending surgeon was almost exclusively (75%) in their first year of practice. BDI was recognized during index operation in 2 cases (50%). Hepatobiliary surgery performed the bile duct repair in all 4 cases. Two complications occurred (50%). All patients were followed by hepatobiliary surgery in the outpatient setting and returned to their baseline level of function within 2 months of hospital discharge.Conclusion: Most BDI occurred in procedures attended by first-year faculty during after hours cholecystectomies, suggesting a role for increased proctorship in early career attendings in addition to in-hours cholecystectomy for acute cholecystitis. The timely return to baseline function experienced by these patients emphasizes the favorable outcomes associated with early recognition of BDI and involvement of hepatobiliary surgery. Further examination with multicenter evaluation would be beneficial to validate these study findings.

Real-Time Cell Cycle Imaging in a 3D Cell Culture Model of Melanoma, Quantitative Analysis, Optical Clearing, and Mathematical Modeling.

Aberrant cell cycle progression is a hallmark of solid tumors. Therefore, cell cycle analysis is an invaluable technique to study cancer cell biology. However, cell cycle progression has been most commonly assessed by methods that are limited to temporal snapshots or that lack spatial information. In this chapter, we describe a technique that allows spatiotemporal real-time tracking of cell cycle progression of individual cells in a multicellular context. The power of this system lies in the use of 3D melanoma spheroids generated from melanoma cells engineered with the fluorescent ubiquitination-based cell cycle indicator (FUCCI). This technique, combined with mathematical modeling, allows us to gain further and more detailed insight into several relevant aspects of solid cancer cell biology, such as tumor growth, proliferation, invasion, and drug sensitivity.

Modelling count data with partial differential equation models in biology.

Partial differential equation (PDE) models are often used to study biological phenomena involving movement-birth-death processes, including ecological population dynamics and the invasion of populations of biological cells. Count data, by definition, is non-negative, and count data relating to biological populations is often bounded above by some carrying capacity that arises through biological competition for space or nutrients. Parameter estimation, parameter identifiability, and making model predictions usually involves working with a measurement error model that explicitly relating experimental measurements with the solution of a mathematical model. In many biological applications, a typical approach is to assume the data are normally distributed about the solution of the mathematical model. Despite the widespread use of the standard additive Gaussian measurement error model, the assumptions inherent in this approach are rarely explicitly considered or compared with other options. Here, we interpret scratch assay data, involving migration, proliferation and delays in a population of cancer cells using a reaction-diffusion PDE model. We consider relating experimental measurements to the PDE solution using a standard additive Gaussian measurement error model alongside a comparison to a more biologically realistic binomial measurement error model. While estimates of model parameters are relatively insensitive to the choice of measurement error model, model predictions for data realisations are very sensitive. The standard additive Gaussian measurement error model leads to biologically inconsistent predictions, such as negative counts and counts that exceed the carrying capacity across a relatively large spatial region within the experiment. Furthermore, the standard additive Gaussian measurement error model requires estimating an additional parameter compared to the binomial measurement error model. In contrast, the binomial measurement error model leads to biologically plausible predictions and is simpler to implement. We provide open source Julia software on GitHub to replicate all calculations in this work, and we explain how to generalise our approach to deal with coupled PDE models with several dependent variables through a multinomial measurement error model, as well as pointing out other potential generalisations by linking our work with established practices in the field of generalised linear models.

Prophylactic Prednisolone Promotes AAV5 Hepatocyte Transduction Through the Novel Mechanism of AAV5 Coreceptor Platelet-Derived Growth Factor Receptor Alpha Upregulation and Innate Immune Suppression.

Adeno-associated virus (AAV) vectors are used to deliver therapeutic transgenes, but host immune responses may interfere with transduction and transgene expression. We evaluated prophylactic corticosteroid treatment on AAV5-mediated expression in liver tissue. Wild-type C57BL/6 mice received 6 × 1013 vg/kg AAV5-HLP-hA1AT, an AAV5 vector carrying a human α1-antitrypsin (hA1AT) gene with a hepatocyte-specific promoter. Mice received 4 weeks of daily 2 mg/kg prednisolone or water starting day -1 or 0 before vector dosing. Mice that received prophylactic corticosteroids had significantly higher serum hA1AT protein than mice that did not, starting at 6 weeks and persisting to the study end at 12 weeks, potentially through a decrease in the number of low responders. RNAseq and proteomic analyses investigating mechanisms mediating the improvement of transgene expression found that prophylactic corticosteroid treatment upregulated the AAV5 coreceptor platelet-derived growth factor receptor alpha (PDGFRα) on hepatocytes and downregulated its competitive ligand PDGFα, thus increasing the uptake of AAV5 vectors. Evidently, prophylactic corticosteroid treatment also suppressed acute immune responses to AAV. Together, these mechanisms resulted in increased uptake and preservation of the transgene, allowing more vector genomes to be available to assemble into stable, full-length structures mediating long-term transgene expression. Prophylactic corticosteroids represent a potential actionable strategy to improve AAV5-mediated transgene expression and decrease intersubject variability.

Formation and Growth of Co-Culture Tumour Spheroids: New Compartment-Based Mathematical Models and Experiments.

Co-culture tumour spheroid experiments are routinely performed to investigate cancer progression and test anti-cancer therapies. Therefore, methods to quantitatively characterise and interpret co-culture spheroid growth are of great interest. However, co-culture spheroid growth is complex. Multiple biological processes occur on overlapping timescales and different cell types within the spheroid may have different characteristics, such as differing proliferation rates or responses to nutrient availability. At present there is no standard, widely-accepted mathematical model of such complex spatio-temporal growth processes. Typical approaches to analyse these experiments focus on the late-time temporal evolution of spheroid size and overlook early-time spheroid formation, spheroid structure and geometry. Here, using a range of ordinary differential equation-based mathematical models and parameter estimation, we interpret new co-culture experimental data. We provide new biological insights about spheroid formation, growth, and structure. As part of this analysis we connect Greenspan's seminal mathematical model to co-culture data for the first time. Furthermore, we generalise a class of compartment-based spheroid mathematical models that have previously been restricted to one population so they can be applied to multiple populations. As special cases of the general model, we explore multiple natural two population extensions to Greenspan's seminal model and reveal biological mechanisms that can describe the internal dynamics of growing co-culture spheroids and those that cannot. This mathematical and statistical modelling-based framework is well-suited to analyse spheroids grown with multiple different cell types and the new class of mathematical models provide opportunities for further mathematical and biological insights.

Distinct Epithelial-Innate Immune Cell Transcriptional Circuits Underlie Airway Hyperresponsiveness in Asthma.

Rationale: Indirect airway hyperresponsiveness (AHR) is a highly specific feature of asthma, but the underlying mechanisms responsible for driving indirect AHR remain incompletely understood. Objectives: To identify differences in gene expression in epithelial brushings obtained from individuals with asthma who were characterized for indirect AHR in the form of exercise-induced bronchoconstriction (EIB). Methods: RNA-sequencing analysis was performed on epithelial brushings obtained from individuals with asthma with EIB (n = 11) and without EIB (n = 9). Differentially expressed genes (DEGs) between the groups were correlated with measures of airway physiology, sputum inflammatory markers, and airway wall immunopathology. On the basis of these relationships, we examined the effects of primary airway epithelial cells (AECs) and specific epithelial cell-derived cytokines on both mast cells (MCs) and eosinophils (EOS). Measurements and Main Results: We identified 120 DEGs in individuals with and without EIB. Network analyses suggested critical roles for IL-33-, IL-18-, and IFN-γ-related signaling among these DEGs. IL1RL1 expression was positively correlated with the density of MCs in the epithelial compartment, and IL1RL1, IL18R1, and IFNG were positively correlated with the density of intraepithelial EOS. Subsequent ex vivo modeling demonstrated that AECs promote sustained type 2 (T2) inflammation in MCs and enhance IL-33-induced T2 gene expression. Furthermore, EOS increase the expression of IFNG and IL13 in response to both IL-18 and IL-33 as well as exposure to AECs. Conclusions: Circuits involving epithelial interactions with MCs and EOS are closely associated with indirect AHR. Ex vivo modeling indicates that epithelial-dependent regulation of these innate cells may be critical in indirect AHR and modulating T2 and non-T2 inflammation in asthma.

Growth and adaptation mechanisms of tumour spheroids with time-dependent oxygen availability.

Tumours are subject to external environmental variability. However, in vitro tumour spheroid experiments, used to understand cancer progression and develop cancer therapies, have been routinely performed for the past fifty years in constant external environments. Furthermore, spheroids are typically grown in ambient atmospheric oxygen (normoxia), whereas most in vivo tumours exist in hypoxic environments. Therefore, there are clear discrepancies between in vitro and in vivo conditions. We explore these discrepancies by combining tools from experimental biology, mathematical modelling, and statistical uncertainty quantification. Focusing on oxygen variability to develop our framework, we reveal key biological mechanisms governing tumour spheroid growth. Growing spheroids in time-dependent conditions, we identify and quantify novel biological adaptation mechanisms, including unexpected necrotic core removal, and transient reversal of the tumour spheroid growth phases.

Rhinovirus infection of the airway epithelium enhances mast cell immune responses via epithelial-derived interferons.

BACKGROUND: Mast cells (MCs) within the airway epithelium in asthma are closely related to airway dysfunction, but cross talk between airway epithelial cells (AECs) and MCs in asthma remains incompletely understood. Human rhinovirus (RV) infections are key triggers for asthma progression, and AECs from individuals with asthma may have dysregulated antiviral responses. OBJECTIVE: We utilized primary AECs in an ex vivo coculture model system to examine cross talk between AECs and MCs after epithelial rhinovirus infection. METHODS: Primary AECs were obtained from 11 children with asthma and 10 healthy children, differentiated at air-liquid interface, and cultured in the presence of laboratory of allergic diseases 2 (LAD2) MCs. AECs were infected with rhinovirus serogroup A 16 (RV16) for 48 hours. RNA isolated from both AECs and MCs underwent RNA sequencing. Direct effects of epithelial-derived interferons on LAD2 MCs were examined by real-time quantitative PCR. RESULTS: MCs increased expression of proinflammatory and antiviral genes in AECs. AECs demonstrated a robust antiviral response after RV16 infection that resulted in significant changes in MC gene expression, including upregulation of genes involved in antiviral responses, leukocyte activation, and type 2 inflammation. Subsequent ex vivo modeling demonstrated that IFN-ß induces MC type 2 gene expression. The effects of AEC donor phenotype were small relative to the effects of viral infection and the presence of MCs. CONCLUSIONS: There is significant cross talk between AECs and MCs, which are present in the epithelium in asthma. Epithelial-derived interferons not only play a role in viral suppression but also further alter MC immune responses including specific type 2 genes.

Computationally efficient framework for diagnosing, understanding and predicting biphasic population growth.

Throughout the life sciences, biological populations undergo multiple phases of growth, often referred to as biphasic growth for the commonly encountered situation involving two phases. Biphasic population growth occurs over a massive range of spatial and temporal scales, ranging from microscopic growth of tumours over several days, to decades-long regrowth of corals in coral reefs that can extend for hundreds of kilometres. Different mathematical models and statistical methods are used to diagnose, understand and predict biphasic growth. Common approaches can lead to inaccurate predictions of future growth that may result in inappropriate management and intervention strategies being implemented. Here, we develop a very general computationally efficient framework, based on profile likelihood analysis, for diagnosing, understanding and predicting biphasic population growth. The two key components of the framework are as follows: (i) an efficient method to form approximate confidence intervals for the change point of the growth dynamics and model parameters and (ii) parameter-wise profile predictions that systematically reveal the influence of individual model parameters on predictions. To illustrate our framework we explore real-world case studies across the life sciences.

Vector genome loss and epigenetic modifications mediate decline in transgene expression of AAV5 vectors produced in mammalian and insect cells.

Recombinant adeno-associated virus (rAAV) vectors are often produced in HEK293 or Spodoptera frugiperda (Sf)-based cell lines. We compared expression profiles of "oversized" (∼5,000 bp) and "standard-sized" (4,600 bp) rAAV5-human α1-antitrypsin (rAAV5-hA1AT) vectors manufactured in HEK293 or Sf cells and investigated molecular mechanisms mediating expression decline. C57BL/6 mice received 6 × 1013 vg/kg of vector, and blood and liver samples were collected through week 57. For all vectors, peak expression (weeks 12-24) declined by 50% to week 57. For Sf- and HEK293-produced oversized vectors, serum hA1AT was initially comparable, but in weeks 12-57, Sf vectors provided significantly higher expression. For HEK293 oversized vectors, liver genomes decreased continuously through week 57 and significantly correlated with A1AT protein. In RNA-sequencing analysis, HEK293 vector-treated mice had significantly higher inflammatory responses in liver at 12 weeks compared with Sf vector- and vehicle-treated mice. Thus, HEK293 vector genome loss led to decreased transgene protein. For Sf-produced vectors, genomes did not decrease from peak expression. Instead, vector genome accessibility significantly decreased from peak to week 57 and correlated with transgene RNA. Vector DNA interactions with active histone marks (H3K27ac/H3K4me3) were significantly reduced from peak to week 57, suggesting that epigenetic regulation impacts transgene expression of Sf-produced vectors.

Tumor Cell Extrinsic Synaptogyrin 3 Expression as a Diagnostic and Prognostic Biomarker in Head and Neck Cancer.

Over 70% of oropharyngeal head and neck squamous cell carcinoma (HNSC) cases in the United States are positive for human papillomavirus (HPV) yet biomarkers for stratifying oropharyngeal head and neck squamous cell carcinoma (HNSC) patient risk are limited. We used immunogenomics to identify differentially expressed genes in immune cells of HPV(+) and HPV(-) squamous carcinomas. Candidate genes were tested in clinical specimens using both quantitative RT-PCR and IHC and validated by IHC using the Carolina Head and Neck Cancer Study (CHANCE) tissue microarray of HNSC cases. We performed multiplex immunofluorescent staining to confirm expression within the immune cells of HPV(+) tumors, receiver operating characteristic (ROC) curve analyses, and assessed survival outcomes. The neuronal gene Synaptogyrin-3 (SYNGR3) is robustly expressed in immune cells of HPV(+) squamous cancers. Multiplex immunostaining and single cell RNA-seq analyses confirmed SYNGR3 expression in T cells, but also unexpectedly in B cells of HPV(+) tumors. ROC curve analyses revealed that combining SYNGR3 and p16 provides more sensitivity and specificity for HPV detection compared to p16 IHC alone. SYNGR3-high HNSC patients have significantly better prognosis with five-year OS and DSS rates of 60% and 71%, respectively. Moreover, combining p16 localization and SYNGR3 expression can further risk stratify HPV(+) patients such that high cytoplasmic, low nuclear p16 do significantly worse (Hazard Ratio, 8.6; P = 0.032) compared to patients with high cytoplasmic, high nuclear p16. SYNGR3 expression in T and B cells is associated with HPV status and enhanced survival outcomes of HNSC patients.

Patient experience before and after treatment with idecabtagene vicleucel (ide-cel, bb2121): qualitative analysis of patient interviews in the KarMMa trial.

OBJECTIVE: To understand the experience of patients with relapsed and refractory multiple myeloma (RRMM) receiving idecabtagene vicleucel (ide-cel), a B-cell maturation antigen-directed chimeric antigen receptor T cell therapy, in the pivotal, phase 2 KarMMa trial. METHODS: Optional semi-structured interviews before leukapheresis (pre-treatment) captured expectations and after ide-cel infusion (1, 2, and 3 months post-treatment), assessed treatment experience, ide-cel advantages/disadvantages, and health and well-being. In a mixed-method analysis, treatment experiences were categorized by clinical response status, health and well-being, and self-reported recovery after infusion. RESULTS: Pre-treatment interviews indicated unmet treatment needs. In post-treatment interviews, most patients reported the positives of ide-cel outweighed negatives (69%, n = 27/39). Most common advantages of ide-cel were efficacy (18-64%), favorable side-effect profile (46-68%), and recovery time (13-18%); most common disadvantages were related to side effects (13-20%). When analyzed by clinical response, patients most often had stringent complete or very good partial response and improved health and well-being with mild or severe recovery from the infusion (27/58, 47%). Most patients with minimal clinical response reported mild infusion recovery (5/6, 83%). CONCLUSIONS: Patient interviews before ide-cel treatment showed unmet needs in triple-class exposed RRMM. Post-treatment experiences generally favored ide-cel versus previously received treatments.

TP53, CDKN2A/P16, and NFE2L2/NRF2 regulate the incidence of pure- and combined-small cell lung cancer in mice.

Studies have shown that Nrf2E79Q/+ is one of the most common mutations found in human tumors. To elucidate how this genetic change contributes to lung cancer, we compared lung tumor development in a genetically-engineered mouse model (GEMM) with dual Trp53/p16 loss, the most common mutations found in human lung tumors, in the presence or absence of Nrf2E79Q/+. Trp53/p16-deficient mice developed combined-small cell lung cancer (C-SCLC), a mixture of pure-SCLC (P-SCLC) and large cell neuroendocrine carcinoma. Mice possessing the LSL-Nrf2E79Q mutation showed no difference in the incidence or latency of C-SCLC compared with Nrf2+/+ mice. However, these tumors did not express NRF2 despite Cre-induced recombination of the LSL-Nrf2E79Q allele. Trp53/p16-deficient mice also developed P-SCLC, where activation of the NRF2E79Q mutation associated with a higher incidence of this tumor type. All C-SCLCs and P-SCLCs were positive for NE-markers, NKX1-2 (a lung cancer marker) and negative for P63 (a squamous cell marker), while only P-SCLC expressed NRF2 by immunohistochemistry. Analysis of a consensus NRF2 pathway signature in human NE+-lung tumors showed variable activation of NRF2 signaling. Our study characterizes the first GEMM that develops C-SCLC, a poorly-studied human cancer and implicates a role for NRF2 activation in SCLC development.

Activation of the CREB Coactivator CRTC2 by Aberrant Mitogen Signaling promotes oncogenic functions in HPV16 positive head and neck cancer.

Head and neck squamous cell carcinoma (HNSCC) is the 6th most common cancer worldwide and incidence rates are continuing to rise globally. Patients often present with locally advanced disease and a staggering 50% chance of relapse following treatment. Aberrant activation of adaptive response signaling pathways, such as the cAMP/PKA pathway, induce an array of genes associated with known cancer pathways that promote tumorigenesis and drug resistance. We identified the cAMP Regulated Transcription Coactivator 2 (CRTC2) to be overexpressed and constitutively activated in HNSCCs and this confers poor prognosis. CRTCs are regulated through their subcellular localization and we show that CRTC2 is exclusively nuclear in HPV(+) HNSCC, thus constitutively active, due to non-canonical Mitogen-Activated Kinase Kinase 1 (MEKK1)-mediated activation via a MEKK1-p38 signaling axis. Loss-of-function and pharmacologic inhibition experiments decreased CRTC2/CREB transcriptional activity by reducing nuclear CRTC2 via nuclear import inhibition and/or by eviction of CRTC2 from the nucleus. This shift in localization was associated with decreased proliferation, migration, and invasion. Our results suggest that small molecules that inhibit nuclear CRTC2 and p38 activity may provide therapeutic benefit to patients with HPV(+) HNSCC.

A stochastic mathematical model of 4D tumour spheroids with real-time fluorescent cell cycle labelling.

In vitro tumour spheroids have been used to study avascular tumour growth and drug design for over 50 years. Tumour spheroids exhibit heterogeneity within the growing population that is thought to be related to spatial and temporal differences in nutrient availability. The recent development of real-time fluorescent cell cycle imaging allows us to identify the position and cell cycle status of individual cells within the growing spheroid, giving rise to the notion of a four-dimensional (4D) tumour spheroid. We develop the first stochastic individual-based model (IBM) of a 4D tumour spheroid and show that IBM simulation data compares well with experimental data using a primary human melanoma cell line. The IBM provides quantitative information about nutrient availability within the spheroid, which is important because it is difficult to measure these data experimentally.

Location of eosinophils in the airway wall is critical for specific features of airway hyperresponsiveness and T2 inflammation in asthma.

BACKGROUND: Eosinophils are implicated as effector cells in asthma, but the functional implications of the precise location of eosinophils in the airway wall is poorly understood. We aimed to quantify eosinophils in the different compartments of the airway wall and associate these findings with clinical features of asthma and markers of airway inflammation. METHODS: In this cross-sectional study, we utilised design-based stereology to accurately partition the numerical density of eosinophils in both the epithelial compartment and the subepithelial space (airway wall area below the basal lamina including the submucosa) in individuals with and without asthma and related these findings to airway hyperresponsiveness (AHR) and features of airway inflammation. RESULTS: Intraepithelial eosinophils were linked to the presence of asthma and endogenous AHR, the type that is most specific for asthma. In contrast, both intraepithelial and subepithelial eosinophils were associated with type 2 (T2) inflammation, with the strongest association between IL5 expression and intraepithelial eosinophils. Eosinophil infiltration of the airway wall was linked to a specific mast cell phenotype that has been described in asthma. We found that interleukin (IL)-33 and IL-5 additively increased cysteinyl leukotriene (CysLT) production by eosinophils and that the CysLT LTC4 along with IL-33 increased IL13 expression in mast cells and altered their protease profile. CONCLUSIONS: We conclude that intraepithelial eosinophils are associated with endogenous AHR and T2 inflammation and may interact with intraepithelial mast cells via CysLTs to regulate airway inflammation.