The metabolic ER stress sensor IRE1α suppresses alternative activation of macrophages and impairs energy expenditure in obesity.

Obesity is associated with metabolic inflammation and endoplasmic reticulum (ER) stress, both of which promote metabolic disease progression. Adipose tissue macrophages (ATMs) are key players orchestrating metabolic inflammation, and ER stress enhances macrophage activation. However, whether ER stress pathways underlie ATM regulation of energy homeostasis remains unclear. Here, we identified inositol-requiring enzyme 1α (IRE1α) as a critical switch governing M1-M2 macrophage polarization and energy balance. Myeloid-specific IRE1α abrogation in Ern1f/f; Lyz2-Cre mice largely reversed high-fat diet (HFD)-induced M1-M2 imbalance in white adipose tissue (WAT) and blocked HFD-induced obesity, insulin resistance, hyperlipidemia and hepatic steatosis. Brown adipose tissue (BAT) activity, WAT browning and energy expenditure were significantly higher in Ern1f/f; Lyz2-Cre mice. Furthermore, IRE1α ablation augmented M2 polarization of macrophages in a cell-autonomous manner. Thus, IRE1α senses protein unfolding and metabolic and immunological states, and consequently guides ATM polarization. The macrophage IRE1α pathway drives obesity and metabolic syndrome through impairing BAT activity and WAT browning.

Safe and just Earth system boundaries.

The stability and resilience of the Earth system and human well-being are inseparably linked1-3, yet their interdependencies are generally under-recognized; consequently, they are often treated independently4,5. Here, we use modelling and literature assessment to quantify safe and just Earth system boundaries (ESBs) for climate, the biosphere, water and nutrient cycles, and aerosols at global and subglobal scales. We propose ESBs for maintaining the resilience and stability of the Earth system (safe ESBs) and minimizing exposure to significant harm to humans from Earth system change (a necessary but not sufficient condition for justice)4. The stricter of the safe or just boundaries sets the integrated safe and just ESB. Our findings show that justice considerations constrain the integrated ESBs more than safety considerations for climate and atmospheric aerosol loading. Seven of eight globally quantified safe and just ESBs and at least two regional safe and just ESBs in over half of global land area are already exceeded. We propose that our assessment provides a quantitative foundation for safeguarding the global commons for all people now and into the future.

The vulnerability of aging states: A survival analysis across premodern societies.

How states and great powers rise and fall is an intriguing enigma of human history. Are there any patterns? Do polities become more vulnerable over time as they age? We analyze longevity in hundreds of premodern states using survival analysis to help provide initial insights into these questions. This approach is commonly used to study the risk of death in biological organisms or failure in mechanical systems. The results reveal that the risk of state termination increased steeply over approximately the first two centuries after formation and stabilized thereafter. This provides the first quantitative support for the hypothesis that the resilience of political states decreases over time. Potential mechanisms that could drive such declining resilience include environmental degradation, increasing complexity, growing inequality, and extractive institutions. While the cases are from premodern times, such dynamics and drivers of vulnerability may remain relevant today.

Long-distance facilitation of coastal ecosystem structure and resilience.

Biotic interactions that hierarchically organize ecosystems by driving ecological and evolutionary processes across spatial scales are ubiquitous in our biosphere. Biotic interactions have been extensively studied at local and global scales, but how long-distance, cross-ecosystem interactions at intermediate landscape scales influence the structure, function, and resilience of ecological systems remains poorly understood. We used remote sensing, modeling, and field data to test the hypothesis that the long-distance impact of an invasive species dramatically affects one of the largest tidal flat ecosystems in East Asia. We found that the invasion of exotic cordgrass Spartina alterniflora can produce long-distance effects on native species up to 10 km away, driving decadal coastal ecosystem transitions. The invasive cordgrass at low elevations facilitated the expansion of the native reed Phragmites australis at high elevations, leading to the massive loss and reduced resilience of the iconic Suaeda salsa "Red Beach" marshes at intermediate elevations, largely as a consequence of reduced soil salinity across the landscape. Our results illustrate the complex role that long-distance interactions can play in shaping landscape structure and ecosystem resilience and in bridging the gap between local and global biotic interactions.

Climate Endgame: Exploring catastrophic climate change scenarios.

Prudent risk management requires consideration of bad-to-worst-case scenarios. Yet, for climate change, such potential futures are poorly understood. Could anthropogenic climate change result in worldwide societal collapse or even eventual human extinction? At present, this is a dangerously underexplored topic. Yet there are ample reasons to suspect that climate change could result in a global catastrophe. Analyzing the mechanisms for these extreme consequences could help galvanize action, improve resilience, and inform policy, including emergency responses. We outline current knowledge about the likelihood of extreme climate change, discuss why understanding bad-to-worst cases is vital, articulate reasons for concern about catastrophic outcomes, define key terms, and put forward a research agenda. The proposed agenda covers four main questions: 1) What is the potential for climate change to drive mass extinction events? 2) What are the mechanisms that could result in human mass mortality and morbidity? 3) What are human societies' vulnerabilities to climate-triggered risk cascades, such as from conflict, political instability, and systemic financial risk? 4) How can these multiple strands of evidence-together with other global dangers-be usefully synthesized into an "integrated catastrophe assessment"? It is time for the scientific community to grapple with the challenge of better understanding catastrophic climate change.

The DEAD-box helicase RCF1 plays roles in miRNA biogenesis and RNA splicing in Arabidopsis.

RCF1 is a highly conserved DEAD-box RNA helicase found in yeast, plants, and mammals. Studies about the functions of RCF1 in plants are limited. Here, we uncovered the functions of RCF1 in Arabidopsis thaliana as a player in pri-miRNA processing and splicing, as well as in pre-mRNA splicing. A mutant with miRNA biogenesis defects was isolated, and the defect was traced to a recessive point mutation in RCF1 (rcf1-4). We show that RCF1 promotes D-body formation and facilitates the interaction between pri-miRNAs and HYL1. Finally, we show that intron-containing pri-miRNAs and pre-mRNAs exhibit a global splicing defect in rcf1-4. Together, this work uncovers roles for RCF1 in miRNA biogenesis and RNA splicing in Arabidopsis.

Non-nociceptive roles of opioids in the CNS: opioids' effects on neurogenesis, learning, memory and affect.

Mortality due to opioid use has grown to the point where, for the first time in history, opioid-related deaths exceed those caused by car accidents in many states in the United States. Changes in the prescribing of opioids for pain and the illicit use of fentanyl (and derivatives) have contributed to the current epidemic. Less known is the impact of opioids on hippocampal neurogenesis, the functional manipulation of which may improve the deleterious effects of opioid use. We provide new insights into how the dysregulation of neurogenesis by opioids can modify learning and affect, mood and emotions, processes that have been well accepted to motivate addictive behaviours.

Electroacupuncture Ameliorates Depression-Like Behaviors Comorbid to Chronic Neuropathic Pain via Tet1-Mediated Restoration of Adult Neurogenesis.

Although electroacupuncture (EA) stimulation is a widely used therapy for chronic pain and comorbid psychiatric disorders, its long-term effects on chronic neuropathic pain-induced depression and the underlying mechanisms remain elusive. In the present study, we found that EA stimulation was able to restore adult neurogenesis in the ventral dentate gyrus (DG), by both increasing neuronal differentiation and restoring the normal morphology of newborn dendrites, in mice with spared nerve injury surgery. By ablating the Nestin+ neural stem cells (NSCs) via diphtheria toxin fragment A expression, we further proved that neurogenesis in the ventral DG was crucial to the long-term, but not the immediate antidepressant effect of EA, nor was it associated with nociception. Furthermore, we found that the restoration of neurogenesis was dependent on Tet1-mediated epigenetic modification upon EA treatment. Tet1 could bind to the promoter of the Prox1 gene, thus catalyzing its demethylation and facilitating its expression, which finally contributed to the restoration of neurogenesis and amelioration of depression-like behaviors induced by chronic neuropathic pain. Thus, we conclude that EA stimulation restores inhibited Tet1 expression in hippocampal NSCs of mice with chronic neuropathic pain, and increased Tet1 expression ameliorates hypermethylation of Prox1 and restores normal adult neurogenesis in the ventral DG, which contributes to the long-term antidepressant effect of EA.

Bioprinted dermis with human adipose tissue-derived microvascular fragments promotes wound healing.

Tissue-engineered skin is an effective material for treating large skin defects in a clinical setting. However, its use is limited owing to vascular complications. Human adipose tissue-derived microvascular fragments (HaMVFs) are vascularized units that form vascular networks by rapid reassembly. In this study, we designed a vascularized bionic skin tissue using a three-dimensional (3D) bioprinter of HaMVFs and human fibroblasts encapsulated in a hybrid hydrogel composed of GelMA, HAMA, and fibrinogen. Tissues incorporating HaMVFs showed good in vitro vascularization and mechanical properties after UV crosslinking and thrombin exposure. Thus, the tissue could be sutured appropriately to the wound. In vivo, the vascularized 3D bioprinted skin promoted epidermal regeneration, collagen maturation in the dermal tissue, and vascularization of the skin tissue to accelerate wound healing. Overall, vascularized 3D bioprinted skin with HaMVFs is an effective material for treating skin defects and may be clinically applicable to reduce the necrosis rate of skin grafts.

Strong and consistent effects of waterbird composition on HPAI H5 occurrences across Europe.

Since 2014, highly pathogenic avian influenza (HPAI) H5 viruses of clade 2.3.4.4 have been dominating the outbreaks across Europe, causing massive deaths among poultry and wild birds. However, the factors shaping these broad-scale outbreak patterns, especially those related to waterbird community composition, remain unclear. In particular, we do not know whether these risk factors differ from those of other H5 clades. Addressing this knowledge gap is important for predicting and preventing future HPAI outbreaks. Using extensive waterbird survey datasets from about 6883 sites, we here explored the effect of waterbird community composition on HPAI H5Nx (clade 2.3.4.4) spatial patterns in the 2016/2017 and 2020/2021 epidemics in Europe, and compared it with the 2005/2006 HPAI H5N1 (clade 2.2) epidemic. We showed that HPAI H5 occurrences in wild birds in the three epidemics were strongly associated with very similar waterbird community attributes, which suggested that, in nature, similar interspecific transmission processes operate between the HPAI H5 subtypes or clades. Importantly, community phylogenetic diversity consistently showed a negative association with H5 occurrence in all three epidemics, suggesting a dilution effect of phylogenetic diversity. In contrast, waterbird community variables showed much weaker associations with HPAI H5Nx occurrence in poultry. Our results demonstrate that models based on previous epidemics can predict future HPAI H5 patterns in wild birds, implying that it is important to include waterbird community factors in future HPAI studies to predict outbreaks and improve surveillance activities.

Fish microplastic ingestion may induce tipping points of aquatic ecosystems.

Microplastics can be ingested by a wide range of aquatic animals. Extensive studies have demonstrated that microplastic ingestion-albeit often not lethal-can affect a range of species life-history traits. However, it remains unclear how the sublethal effects of microplastics on individual levels scale up to influence ecosystem-level dynamics through cascading trophic interactions. Here we employ a well-studied, empirically fed three-species trophic chain model, which was parameterized to mimic a common type of aquatic ecosystems to examine how microplastic ingestion by fish on an intermediate trophic level can produce cascading effects on the species at both upper and lower trophic levels. We show that gradually increasing microplastics in the ingested substances of planktivorous fish may cause population structure effects such as skewed size distributions (i.e. reduced average body length vs. increased maximal body size), and induce abrupt declines in fish biomass and reproduction. Our model analysis demonstrates that these abrupt changes correspond to an ecosystem-level tipping point, crossing which difficult-to-reverse ecosystem degradation can happen. Importantly, microplastic pollution may interact with other anthropogenic stressors to reduce safe operating space of aquatic ecosystems. Our work contributes to better understanding complex effects of microplastic pollution and anticipating tipping points of aquatic ecosystems in a changing world. It also calls attention to an emerging threat that novel microplastic contaminants may lead to unexpected and abrupt degradation of aquatic ecosystems, and invites systematic studies on the ecosystem-level consequences of microplastic exposure.

Electroacupuncture alleviates the relapse of pain-related aversive memory by activating KOR and inhibiting GABAergic neurons in the insular cortex.

Pain-related aversive memory is common in chronic pain patients. Electroacupuncture has been demonstrated to block pain-related aversive memory. The insular cortex is a key region closely related to aversive behaviors. In our study, a potential mechanism underlying the effect of electroacupuncture treatment on pain-related aversive memory behaviors relative to the insular cortex was investigated. Our study used the chemogenetic method, pharmacological method, electroacupuncture intervention, and behavioral detection. Our study showed that both inhibition of gamma-aminobutyric acidergic neurons and activation of the kappa opioid receptor in the insular cortex blocked the pain-related aversive memory behaviors induced by 2 crossover injections of carrageenan in mice; conversely, both the activation of gamma-aminobutyric acidergic neurons and inhibition of kappa opioid receptor in the insular cortex play similar roles in inducing pain-related aversive memory behaviors following 2 crossover injections of carrageenan. In addition, activation of gamma-aminobutyric acidergic neurons in the insular cortex reversed the effect of kappa opioid receptor activation in the insular cortex. Moreover, electroacupuncture effectively blocked pain-related aversive memory behaviors in model mice, which was reversed by both activation of gamma-aminobutyric acidergic neurons and inhibition of kappa opioid receptor in the insular cortex. The effect of electroacupuncture on blocking pain-related aversive memory behaviors may be related to the activation of the kappa opioid receptor and inhibition of gamma-aminobutyric acidergic neurons in the insular cortex.

Effect of deep neuromuscular block on the quality of early recovery after sleeve gastrectomy in obese patients: a randomized controlled trial.

BACKGROUND: Deep neuromuscular block (NMB) has been shown to improve surgical conditions and alleviate post-operative pain in bariatric surgery compared with moderate NMB. We hypothesized that deep NMB could also improve the quality of early recovery after laparoscopic sleeve gastrectomy (LSG). METHODS: Eighty patients were randomized to receive either deep (post-tetanic count 1-3) or moderate (train-of-four count 1-3) NMB. The QoR-15 questionnaire was used to evaluate the quality of early recovery at 1 day before surgery (T0), 24 and 48 h after surgery (T2, T3). Additionally, we recorded diaphragm excursion (DE), postoperative pain, surgical condition, cumulative dose of analgesics, time of first flatus and ambulation, post-operative nausea and vomiting, time of tracheal tube removal and hospitalization time. MAIN RESULTS: The quality of recovery was significantly better 24 h after surgery in patients who received a deep versus moderate block (114.4 ± 12.9 versus 102.1 ± 18.1). Diaphragm excursion was significantly greater in the deep NMB group when patients performed maximal inspiration at T2 and T3 (P < 0.05). Patients who underwent deep NMB reported lower visceral pain scores 40 min after surgery; additionally, these patients experienced lower pain during movement at T3 (P < 0.05). Optimal surgical conditions were rated in 87.5% and 64.6% of all measurements during deep and moderate NMB respectively (P < 0.001). The time to tracheal tube removal was significantly longer in the deep NMB group (P = 0.001). There were no differences in other outcomes. CONCLUSION: In obese patients receiving deep NMB during LSG, we observed improved QoR-15 scores, greater diaphragmatic excursions, improved surgical conditions, and visceral pain scores were lower. More evidence is needed to determine the effects of deep NMB on these outcomes. TRIAL REGISTRATION: ChiCTR2200065919. Date of retrospectively registered: 18/11/2022.

Global inequality remotely sensed.

Economic inequality is notoriously difficult to quantify as reliable data on household incomes are missing for most of the world. Here, we show that a proxy for inequality based on remotely sensed nighttime light data may help fill this gap. Individual households cannot be remotely sensed. However, as households tend to segregate into richer and poorer neighborhoods, the correlation between light emission and economic thriving shown in earlier studies suggests that spatial variance of remotely sensed light per person might carry a signal of economic inequality. To test this hypothesis, we quantified Gini coefficients of the spatial variation in average nighttime light emitted per person. We found a significant relationship between the resulting light-based inequality indicator and existing estimates of net income inequality. This correlation between light-based Gini coefficients and traditional estimates exists not only across countries, but also on a smaller spatial scale comparing the 50 states within the United States. The remotely sensed character makes it possible to produce high-resolution global maps of estimated inequality. The inequality proxy is entirely independent from traditional estimates as it is based on observed light emission rather than self-reported household incomes. Both are imperfect estimates of true inequality. However, their independent nature implies that the light-based proxy could be used to constrain uncertainty in traditional estimates. More importantly, the light-based Gini maps may provide an estimate of inequality where previously no data were available at all.

Evaluation of Low-Cost CO2 Sensors Using Reference Instruments and Standard Gases for Indoor Use.

CO2 monitoring is important for carbon emission evaluation. Low-cost and medium-precision sensors (LCSs) have become an exploratory direction for CO2 observation under complex emission conditions in cities. Here, we used a calibration method that improved the accuracy of SenseAir K30 CO2 sensors from ±30 ppm to 0.7-4.0 ppm for a CO2-monitoring instrument named the SENSE-IAP, which has been used in several cities, such as in Beijing, Jinan, Fuzhou, Hangzhou, and Wuhan, in China since 2017. We conducted monthly to yearly synchronous observations using the SENSE-IAP along with reference instruments (Picarro) and standard gas to evaluate the performance of the LCSs for indoor use with relatively stable environments. The results show that the precision and accuracy of the SENSE-IAP compared to the standard gases were rather good in relatively stable indoor environments, with the short-term (daily scale) biases ranging from -0.9 to 0.2 ppm, the root mean square errors (RMSE) ranging from 0.7 to 1.6 ppm, the long-term (monthly scale) bias ranging from -1.6 to 0.5 ppm, and the RMSE ranging from 1.3 to 3.2 ppm. The accuracy of the synchronous observations with Picarro was in the same magnitude, with an RMSE of 2.0-3.0 ppm. According to our evaluation, standard instruments or reliable standard gases can be used as a reference to improve the accuracy of the SENSE-IAP. If calibrated daily using standard gases, the bias of the SENSE-IAP can be maintained within 1.0 ppm. If the standard gases are hard to access frequently, we recommend a calibration frequency of at least three months to maintain an accuracy within 3 ppm.

Parallel analysis of RNA ends reveals global microRNA-mediated target RNA cleavage in maize.

MicroRNAs (miRNAs) are endogenous 20-24-nucleotide non-coding RNAs that play important regulatory roles in many biological processes in eukaryotes. miRNAs modulate the expression of target genes at the post-transcriptional level by transcript cleavage or translational inhibition. The identification of miRNA target genes has been extensively investigated in Arabidopsis and rice, but an in-depth global analysis of miRNA-mediated target regulation is still lacking in maize. Here, we report a transcriptome-wide identification of miRNA targets by analyzing parallel analysis of RNA ends (PARE) datasets derived from nine different tissues at five developmental stages of the maize (Zea mays L.) B73 cultivar. In total, 246 targets corresponding to 60 miRNAs from 25 families were identified, including transcription factors and other genes. In addition, PARE analysis revealed that miRNAs guide specific target transcript cleavage in a tissue-preferential manner. Primary transcripts of MIR159c and MIR169e were found to be cleaved by mature miR159 and miR169, respectively, indicating a negative-feedback regulatory mechanism in miRNA biogenesis. Moreover, several miRNA-target gene pairs involved in seed germination were identified and experimentally validated. Our PARE analyses generated a wide and detailed miRNA-target interaction atlas, which provides a valuable resource for investigating the roles of miRNAs and their targets in maize.

Functional traits explain waterbirds' host status, subtype richness, and community-level infection risk for avian influenza.

Species functional traits can influence pathogen transmission processes, and consequently affect species' host status, pathogen diversity, and community-level infection risk. We here investigated, for 143 European waterbird species, effects of functional traits on host status and pathogen diversity (subtype richness) for avian influenza virus at species level. We then explored the association between functional diversity and HPAI H5Nx occurrence at the community level for 2016/17 and 2021/22 epidemics in Europe. We found that both host status and subtype richness were shaped by several traits, such as diet guild and dispersal ability, and that the community-weighted means of these traits were also correlated with community-level risk of H5Nx occurrence. Moreover, functional divergence was negatively associated with H5Nx occurrence, indicating that functional diversity can reduce infection risk. Our findings highlight the value of integrating trait-based ecology into the framework of diversity-disease relationship, and provide new insights for HPAI prediction and prevention.

Identification of genetic and immune signatures for the recurrence of HER2-positive breast cancer after trastuzumab-based treatment.

PURPOSE: To determine the genetic and immune features associated with the recurrence of human epidermal growth factor receptor2-positive (HER2 +) breast cancer (BC) after trastuzumab-based treatment. METHODS: A retrospective cohort study of 48 patients who received trastuzumab-based treatment was divided into recurrent and non-recurrent groups according to clinical follow-up. Baseline samples from all 48 patients were analyzed for genetic variation, HLA allele type, gene expression, and immune features, which were linked to HER2 + BC recurrence. Statistics included logistic regression models, Kaplan-Meier plots, and Univariate Cox proportional hazards models. RESULTS: Compared with the non-recurrent group, the extracellular matrix-related pathway and 3 Hallmark gene sets were enriched in the recurrent group. The infiltration levels of immature B cells and activated B cells were significantly increased in the non-recurrent group, which correlated remarkably with improved overall survival (OS) in two other published gene expression datasets, including TCGA and METABRIC. In the TCGA cohort (n = 275), activated B cells (HR 0.23, 95%CI 0.13-0.43, p < 0.0001), and immature B cells (HR 0.26, 95%CI 0.12-0.59, p < 0.0001). In the METABRIC cohort (n = 236), activated B cells (HR 0.60, 95%CI 0.43-0.83, p = 0.002), and immature B cells (HR 0.65, 95%CI 0.47-0.91, p = 0.011). Cox regression suggested that immature B cells and activated B cells were protective factors for outcome OS. CONCLUSIONS: Aberrant activation of multiple pathways and low baseline tumor-infiltrating B cells are related to HER2 + BC trastuzumab-based recurrence, which primarily affects the antitumor activity of trastuzumab.

The Agonist of Adenosine A1 Receptor Induced Desensitization of delta Opioid receptor-mediated Raf-1/MEK/ERK Signaling by Feedback Phosphorylation of Raf-1-Ser289/296/301.

Our previous study found that activation of adenosine A1 receptor (A1R) induced phosphorylation of delta opioid receptor (DOR) and desensitization of its downstream signaling molecules, cAMP and Akt. To further investigate the effect of A1R agonist on DOR signaling and the underlying mechanism, we examined the effect of A1R activation upon binding of its agonist N6-cyclohexyl-adenosine (CHA) on DOR-mediated Raf-1/MEK/ERK activation, and found that prolonged CHA exposure resulted in downregulation of DOR-mediated Raf-1/MEK/ERK signaling pathway. CHA-treatment time dependently attenuated Raf-1-Ser338 phosphorylation induced by [D-Pen2,5] enkephalin (DPDPE), a specific agonist of DOR, and further caused downregulation of the Raf-1/MEK/ERK signaling pathway activated by DOR agonist. Moreover, CHA exposure time-dependently induced the phosphorylation of Raf-1-Ser289/296/301, the inhibitory phosphorylation sites that were regulated by negative feedback, thereby inhibiting activation of the MEK/ERK pathway, and this effect could be blocked by MEK inhibitor U0126. Finally, we proved that the heterologous desensitization of the Raf-1/MEK/ERK cascade was essential in the regulation of anti-nociceptive effect of DOR agonists by confirming that such effect was inhibited by pretreatment of CHA. Therefore, we conclude that the activation of A1R inhibits DOR-mediated MAPK signaling pathway via heterologous desensitization of the Raf-1/MEK/ERK cascade, which is a result of ERK-mediated Raf-1-Ser289/296/301 phosphorylation mediated by activation of A1R.

Avalanches during epithelial tissue growth; Uniform Growth and a drosophila eye disc model.

Epithelial tissues constitute an exotic type of active matter with non-linear properties reminiscent of amorphous materials. In the context of a proliferating epithelium, modeled by the quasistatic vertex model, we identify novel discrete tissue scale rearrangements, i.e. cellular rearrangement avalanches, which are a form of collective cell movement. During the avalanches, the vast majority of cells retain their neighbors, and the resulting cellular trajectories are radial in the periphery, a vortex in the core. After the onset of these avalanches, the epithelial area grows discontinuously. The avalanches are found to be stochastic, and their strength is correlated with the density of cells in the tissue. Overall, avalanches redistribute accumulated local spatial pressure along the tissue. Furthermore, the distribution of avalanche magnitudes is found to obey a power law, with an exponent consistent with sheer induced avalanches in amorphous materials. To understand the role of avalanches in organ development, we simulate epithelial growth of the Drosophila eye disc during the third instar using a computational model, which includes both chemical and mechanistic signaling. During the third instar, the morphogenetic furrow (MF), a ~10 cell wide wave of apical area constriction propagates through the epithelium. These simulations are used to understand the details of the growth process, the effect of the MF on the growth dynamics on the tissue scale, and to make predictions for experimental observations. The avalanches are found to depend on the strength of the apical constriction of cells in the MF, with a stronger apical constriction leading to less frequent and more pronounced avalanches. The results herein highlight the dependence of simulated tissue growth dynamics on relaxation timescales, and serve as a guide for in vitro experiments.