The cellular mechanisms of neuronal swelling underlying cytotoxic edema.

Cytotoxic brain edema triggered by neuronal swelling is the chief cause of mortality following brain trauma and cerebral infarct. Using fluorescence lifetime imaging to analyze contributions of intracellular ionic changes in brain slices, we find that intense Na(+) entry triggers a secondary increase in intracellular Cl(-) that is required for neuronal swelling and death. Pharmacological and siRNA-mediated knockdown screening identified the ion exchanger SLC26A11 unexpectedly acting as a voltage-gated Cl(-) channel that is activated upon neuronal depolarization to membrane potentials lower than -20 mV. Blockade of SLC26A11 activity attenuates both neuronal swelling and cell death. Therefore cytotoxic neuronal edema occurs when sufficient Na(+) influx and depolarization is followed by Cl(-) entry via SLC26A11. The resultant NaCl accumulation causes subsequent neuronal swelling leading to neuronal death. These findings shed light on unique elements of volume control in excitable cells and lay the ground for the development of specific treatments for brain edema.

Natural antisense transcript of Period2, Per2AS, regulates the amplitude of the mouse circadian clock.

In mammals, a set of core clock genes form transcription-translation feedback loops to generate circadian oscillations. We and others recently identified a novel transcript at the Period2 (Per2) locus that is transcribed from the antisense strand of Per2 This transcript, Per2AS, is expressed rhythmically and antiphasic to Per2 mRNA, leading to our hypothesis that Per2AS and Per2 mutually inhibit each other's expression and form a double negative feedback loop. By perturbing the expression of Per2AS, we found that Per2AS transcription, but not transcript, represses Per2 However, Per2 does not repress Per2AS, as Per2 knockdown led to a decrease in the Per2AS level, indicating that Per2AS forms a single negative feedback loop with Per2 and maintains the level of Per2 within the oscillatory range. Per2AS also regulates the amplitude of the circadian clock, and this function cannot be solely explained through its interaction with Per2, as Per2 knockdown does not recapitulate the phenotypes of Per2AS perturbation. Overall, our data indicate that Per2AS is an important regulatory molecule in the mammalian circadian clock machinery. Our work also supports the idea that antisense transcripts of core clock genes constitute a common feature of circadian clocks, as they are found in other organisms.

The oscillation of mitotic kinase governs cell cycle latches in mammalian cells.

The mammalian cell cycle alternates between two phases - S-G2-M with high levels of A- and B-type cyclins (CycA and CycB, respectively) bound to cyclin-dependent kinases (CDKs), and G1 with persistent degradation of CycA and CycB by an activated anaphase promoting complex/cyclosome (APC/C) bound to Cdh1 (also known as FZR1 in mammals; denoted APC/C:Cdh1). Because CDKs phosphorylate and inactivate Cdh1, these two phases are mutually exclusive. This 'toggle switch' is flipped from G1 to S by cyclin-E bound to a CDK (CycE:CDK), which is not degraded by APC/C:Cdh1, and from M to G1 by Cdc20-bound APC/C (APC/C:Cdc20), which is not inactivated by CycA:CDK or CycB:CDK. After flipping the switch, cyclin E is degraded and APC/C:Cdc20 is inactivated. Combining mathematical modelling with single-cell timelapse imaging, we show that dysregulation of CycB:CDK disrupts strict alternation of the G1-S and M-G1 switches. Inhibition of CycB:CDK results in Cdc20-independent Cdh1 'endocycles', and sustained activity of CycB:CDK drives Cdh1-independent Cdc20 endocycles. Our model provides a mechanistic explanation for how whole-genome doubling can arise, a common event in tumorigenesis that can drive tumour evolution.

A dynamical model of growth and maturation in Drosophila.

The decision to stop growing and mature into an adult is a critical point in development that determines adult body size, impacting multiple aspects of an adult's biology. In many animals, growth cessation is a consequence of hormone release that appears to be tied to the attainment of a particular body size or condition. Nevertheless, the size-sensing mechanism animals use to initiate hormone synthesis is poorly understood. Here, we develop a simple mathematical model of growth cessation in Drosophila melanogaster, which is ostensibly triggered by the attainment of a critical weight (CW) early in the last instar. Attainment of CW is correlated with the synthesis of the steroid hormone ecdysone, which causes a larva to stop growing, pupate, and metamorphose into the adult form. Our model suggests that, contrary to expectation, the size-sensing mechanism that initiates metamorphosis occurs before the larva reaches CW; that is, the critical-weight phenomenon is a downstream consequence of an earlier size-dependent developmental decision, not a decision point itself. Further, this size-sensing mechanism does not require a direct assessment of body size but emerges from the interactions between body size, ecdysone, and nutritional signaling. Because many aspects of our model are evolutionarily conserved among all animals, the model may provide a general framework for understanding how animals commit to maturing from their juvenile to adult form.

Feedback in the ß-catenin destruction complex imparts bistability and cellular memory.

Wnt ligands are considered classical morphogens, for which the strength of the cellular response is proportional to the concentration of the ligand. Herein, we show an emergent property of bistability arising from feedback among the Wnt destruction complex proteins that target the key transcriptional co-activator ß-catenin for degradation. Using biochemical reconstitution, we identified positive feedback between the scaffold protein Axin and the kinase glycogen synthase kinase 3 (GSK3). Theoretical modeling of this feedback between Axin and GSK3 suggested that the activity of the destruction complex exhibits bistable behavior. We experimentally confirmed these predictions by demonstrating that cellular cytoplasmic ß-catenin concentrations exhibit an "all-or-none" response with sustained memory (hysteresis) of the signaling input. This bistable behavior was transformed into a graded response and memory was lost through inhibition of GSK3. These findings provide a mechanism for establishing decisive, switch-like cellular response and memory upon Wnt pathway stimulation.

Rapid, high-throughput, cost-effective whole-genome sequencing of SARS-CoV-2 using a condensed library preparation of the Illumina DNA Prep kit.

The ongoing COVID-19 pandemic necessitates cost-effective, high-throughput, and timely whole-genome sequencing (WGS) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viruses for outbreak investigations, identifying variants of concern (VoC), characterizing vaccine breakthrough infections, and public health surveillance. In addition, the enormous demand for WGS on supply chains and the resulting shortages of laboratory supplies necessitated the use of low-reagent and low-consumable methods. Here, we report an optimized library preparation method (the BCCDC cutdown method) that can be used in a high-throughput scenario, where one technologist can perform 576 library preparations (6 plates of 96 samples) over the course of one 8-hour shift. The same protocol can also be used in a rapid turnaround time scenario, from primary samples (up to 96 samples) to loading on a sequencer in an 8-hour shift. This new method uses Freed et al.'s 1,200 bp primer sets (Biol Methods Protoc 5:bpaa014, 2020, https://doi.org/10.1093/biomethods/bpaa014) and a modified and condensed Illumina DNA Prep workflow (Illumina, CA, USA). Compared to the original protocol, the application of this new method using hundreds of clinical specimens demonstrated equivalent results to the full-length DNA Prep workflow at 45% of the cost, 15% of consumables required (such as pipet tips), 25% of manual hands-on time, and 15% of on-instrument time if performing on a liquid handler, with no compromise in sequence quality. Results demonstrate that this new method is a rapid, simple, cost-effective, and high-quality SARS-CoV-2 WGS protocol. IMPORTANCE: Sequencing has played an invaluable role in the response to the COVID-19 pandemic. Ongoing work in this area, however, demands optimization of laboratory workflow to increase sequencing capacity, improve turnaround time, and reduce cost without compromising sequence quality. This report describes an optimized DNA library preparation method for improved whole-genome sequencing of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pathogen. The workflow advantages summarized here include significant time, cost, and consumable savings, which suggest that this new method is an efficient, scalable, and pragmatic alternative for SARS-CoV-2 whole-genome sequencing.

Clinical Severity of Severe Acute Respiratory Syndrome Coronavirus 2 Omicron Variant Relative to Delta in British Columbia, Canada: A Retrospective Analysis of Whole-Genome Sequenced Cases.

BACKGROUND: In late 2021, the Omicron severe acute respiratory syndrome coronavirus 2 variant emerged and rapidly replaced Delta as the dominant variant. The increased transmissibility of Omicron led to surges in case rates and hospitalizations; however, the true severity of the variant remained unclear. We aimed to provide robust estimates of Omicron severity relative to Delta. METHODS: This retrospective cohort study was conducted with data from the British Columbia COVID-19 Cohort, a large provincial surveillance platform with linkage to administrative datasets. To capture the time of cocirculation with Omicron and Delta, December 2021 was chosen as the study period. Whole-genome sequencing was used to determine Omicron and Delta variants. To assess the severity (hospitalization, intensive care unit [ICU] admission, length of stay), we conducted adjusted Cox proportional hazard models, weighted by inverse probability of treatment weights (IPTW). RESULTS: The cohort was composed of 13 128 individuals (7729 Omicron and 5399 Delta). There were 419 coronavirus disease 2019 hospitalizations, with 118 (22%) among people diagnosed with Omicron (crude rate = 1.5% Omicron, 5.6% Delta). In multivariable IPTW analysis, Omicron was associated with a 50% lower risk of hospitalization compared with Delta (adjusted hazard ratio [aHR] = 0.50, 95% confidence interval [CI] = 0.43 to 0.59), a 73% lower risk of ICU admission (aHR = 0.27, 95% CI = 0.19 to 0.38), and a 5-day shorter hospital stay (aß = -5.03, 95% CI = -8.01 to -2.05). CONCLUSIONS: Our analysis supports findings from other studies that have demonstrated lower risk of severe outcomes in Omicron-infected individuals relative to Delta.

Emergence of SARS-CoV-2 Delta Variant and Effect of Nonpharmaceutical Interventions, British Columbia, Canada.

In British Columbia, Canada, initial growth of the SARS-CoV-2 Delta variant was slower than that reported in other jurisdictions. Delta became the dominant variant (>50% prevalence) within ≈7-13 weeks of first detection in regions within the United Kingdom and United States. In British Columbia, it remained at <10% of weekly incident COVID-19 cases for 13 weeks after first detection on March 21, 2021, eventually reaching dominance after 17 weeks. We describe the growth of Delta variant cases in British Columbia during March 1-June 30, 2021, and apply retrospective counterfactual modeling to examine factors for the initially low COVID-19 case rate after Delta introduction, such as vaccination coverage and nonpharmaceutical interventions. Growth of COVID-19 cases in the first 3 months after Delta emergence was likely limited in British Columbia because additional nonpharmaceutical interventions were implemented to reduce levels of contact at the end of March 2021, soon after variant emergence.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Clinical severity of Omicron subvariants BA.1, BA.2, and BA.5 in a population-based cohort study in British Columbia, Canada.

The SARS-CoV-2 variant Omicron emerged in late 2021. In British Columbia (BC), Canada, and globally, three genetically distinct subvariants of Omicron, BA.1, BA.2, and BA.5, emerged and became dominant successively within an 8-month period. SARS-CoV-2 subvariants continue to circulate in the population, acquiring new mutations that have the potential to alter infectivity, immunity, and disease severity. Here, we report a propensity-matched severity analysis from residents of BC over the course of the Omicron wave, including 39,237 individuals infected with BA.1, BA.2, or BA.5 based on paired high-quality sequence data and linked to comprehensive clinical outcomes data between December 23, 2021 and August 31, 2022. Relative to BA.1, BA.2 cases were associated with a 15% and 28% lower risk of hospitalization and intensive care unit (ICU) admission (aHRhospital = 1.17; 95% confidence interval [CI] = 1.096-1.252; aHRICU = 1.368; 95% CI = 1.152-1.624), whereas BA.5 infections were associated with an 18% higher risk of hospitalization (aHRhospital = 1.18; 95% CI = 1.133-1.224) after accounting for age, sex, comorbidities, vaccination status, geography, and social determinants of health. Phylogenetic analysis revealed no specific subclades associated with more severe clinical outcomes for any Omicron subvariant. In summary, BA.1, BA.2, and BA.5 subvariants were associated with differences in clinical severity, emphasizing how variant-specific monitoring programs remain critical components of patient and population-level public health responses as the pandemic continues.

Mathematical analysis of robustness of oscillations in models of the mammalian circadian clock.

Circadian rhythms in a wide range of organisms are mediated by molecular mechanisms based on transcription-translation feedback. In this paper, we use bifurcation theory to explore mathematical models of genetic oscillators, based on Kim & Forger's interpretation of the circadian clock in mammals. At the core of their models is a negative feedback loop whereby PER proteins (PER1 and PER2) bind to and inhibit their transcriptional activator, BMAL1. For oscillations to occur, the dissociation constant of the PER:BMAL1 complex, [Formula: see text], must be ≤ 0.04 nM, which is orders of magnitude smaller than a reasonable expectation of 1-10 nM for this protein complex. We relax this constraint by two modifications to Kim & Forger's 'single negative feedback' (SNF) model: first, by introducing a multistep reaction chain for posttranscriptional modifications of Per mRNA and posttranslational phosphorylations of PER, and second, by replacing the first-order rate law for degradation of PER in the nucleus by a Michaelis-Menten rate law. These modifications increase the maximum allowable [Formula: see text] to ~2 nM. In a third modification, we consider an alternative rate law for gene transcription to resolve an unrealistically large rate of Per2 transcription at very low concentrations of BMAL1. Additionally, we studied extensions of the SNF model to include a second negative feedback loop (involving REV-ERB) and a supplementary positive feedback loop (involving ROR). Contrary to Kim & Forger's observations of these extended models, we find that, with our modifications, the supplementary positive feedback loop makes the oscillations more robust than observed in the models with one or two negative feedback loops. However, all three models are similarly robust when accounting for circadian rhythms (~24 h period) with [Formula: see text] ≥ 1 nM. Our results provide testable predictions for future experimental studies.

From the Belousov-Zhabotinsky reaction to biochemical clocks, traveling waves and cell cycle regulation.

In the last 20 years, a growing army of systems biologists has employed quantitative experimental methods and theoretical tools of data analysis and mathematical modeling to unravel the molecular details of biological control systems with novel studies of biochemical clocks, cellular decision-making, and signaling networks in time and space. Few people know that one of the roots of this new paradigm in cell biology can be traced to a serendipitous discovery by an obscure Russian biochemist, Boris Belousov, who was studying the oxidation of citric acid. The story is told here from an historical perspective, tracing its meandering path through glycolytic oscillations, cAMP signaling, and frog egg development. The connections among these diverse themes are drawn out by simple mathematical models (nonlinear differential equations) that share common structures and properties.

Comparative Single-Dose mRNA and ChAdOx1 Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2, Including Variants of Concern: Test-Negative Design, British Columbia, Canada.

BACKGROUND: In British Columbia, Canada, most adults 50-69 years old became eligible for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) vaccine in April 2021, with chimpanzee adenoviral vectored vaccine (ChAdOx1) restricted to ≥55-year-olds and second doses deferred ≥6 weeks to optimize single-dose coverage. METHODS: Among adults 50-69 years old, single-dose messenger RNA (mRNA) and ChAdOx1 vaccine effectiveness (VE) against SARS-CoV-2 infection and hospitalization, including variant-specific, was assessed by test-negative design between 4 April and 2 October 2021. RESULTS: Single-dose VE included 11 861 cases and 99 544 controls. Median of postvaccination follow-up was 32 days (interquartile range, 15-52 days). Alpha, Gamma, and Delta variants comprised 23%, 18%, and 56%, respectively, of genetically characterized viruses. At 21-55 days postvaccination, single-dose mRNA and ChAdOx1 VE (95% confidence interval [CI]) was 74% (71%-76%) and 59% (53%-65%) against any infection and 86% (80%-90%) and 94% (85%-97%) against hospitalization, respectively. VE (95% CI) was similar against Alpha and Gamma infections for mRNA (80% [76%-84%] and 80% [75%-84%], respectively) and ChAdOx1 (69% [60%-76%] and 66% [56%-73%], respectively). mRNA VE was lower at 63% (95% CI, 56%-69%) against Delta but 85% (95% CI, 71%-92%) against Delta-associated hospitalization (nonestimable for ChAdOx1). CONCLUSIONS: A single mRNA or ChAdOx1 vaccine dose gave important protection against SARS-CoV-2, including early variants of concern. ChAdOx1 VE was lower against infection, but 1 dose of either vaccine reduced the hospitalization risk by >85% to at least 8 weeks postvaccination. Findings inform program options, including longer dosing intervals.

Single-dose mRNA Vaccine Effectiveness Against Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), Including Alpha and Gamma Variants: A Test-negative Design in Adults 70 Years and Older in British Columbia, Canada.

BACKGROUND: Randomized-controlled trials of messenger RNA (mRNA) vaccine protection against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) included relatively few elderly participants. We assess single-dose mRNA vaccine effectiveness (VE) in adults ≥ 70 years old in British Columbia, Canada, where second doses were deferred by up to 16 weeks and where a spring 2021 wave uniquely included codominant circulation of Alpha (B.1.1.7) and Gamma (P.1) variants of concern (VOC). METHODS: Analyses included community-dwelling adults ≥ 70 years old with specimen collection between 4 April (epidemiological week 14) and 1 May (week 17) 2021. Adjusted VE was estimated by test-negative design. Cases were reverse-transcription polymerase chain reaction (RT-PCR) test-positive for SARS-CoV-2, and controls were test-negative. Vaccine status was defined by receipt of a single-dose ≥ 21 days before specimen collection, but a range of intervals was assessed. Variant-specific VE was estimated against viruses genetically characterized as Alpha, Gamma or non-VOC lineages. RESULTS: VE analyses included 16 993 specimens: 1226 (7%) test-positive cases and 15 767 test-negative controls. Of 1131 (92%) genetically characterized viruses, 509 (45%), 314 (28%), and 276 (24%) were Alpha, Gamma, and non-VOC lineages, respectively. At 0-13 days postvaccination, VE was negligible at 14% (95% confidence interval [CI], 0-26) but increased from 43% (95% CI, 30-53) at 14-20 days to 75% (95% CI, 63-83) at 35-41 days postvaccination. VE at ≥ 21 days postvaccination was 65% (95% CI, 58-71) overall: 72% (95% CI, 58-81), 67% (95% CI, 57-75), and 61% (95% CI, 45-72) for non-VOC, Alpha, and Gamma variants, respectively. CONCLUSIONS: A single dose of mRNA vaccine reduced the risk of SARS-CoV-2 by about two-thirds in adults ≥ 70 years old, with protection only minimally reduced against Alpha and Gamma variants.

Two-Dose Severe Acute Respiratory Syndrome Coronavirus 2 Vaccine Effectiveness With Mixed Schedules and Extended Dosing Intervals: Test-Negative Design Studies From British Columbia and Quebec, Canada.

BACKGROUND: The Canadian coronavirus disease 2019 (COVID-19) immunization strategy deferred second doses and allowed mixed schedules. We compared 2-dose vaccine effectiveness (VE) by vaccine type (mRNA and/or ChAdOx1), interval between doses, and time since second dose in 2 of Canada's larger provinces. METHODS: Two-dose VE against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection or hospitalization among adults ≥18 years, including due to Alpha, Gamma, and Delta variants of concern (VOCs), was assessed ≥14 days postvaccination by test-negative design studies separately conducted in British Columbia and Quebec, Canada, between 30 May and 27 November (epi-weeks 22-47) 2021. RESULTS: In both provinces, all homologous or heterologous mRNA and/or ChAdOx1 2-dose schedules were associated with ≥90% reduction in SARS-CoV-2 hospitalization risk for ≥7 months. With slight decline from a peak of >90%, VE against infection was ≥80% for ≥6 months following homologous mRNA vaccination, lower by ∼10% when both doses were ChAdOx1 but comparably high following heterologous ChAdOx1 + mRNA receipt. Findings were similar by age group, sex, and VOC. VE was significantly higher with longer 7-8-week versus manufacturer-specified 3-4-week intervals between mRNA doses. CONCLUSIONS: Two doses of any mRNA and/or ChAdOx1 combination gave substantial and sustained protection against SARS-CoV-2 hospitalization, spanning Delta-dominant circulation. ChAdOx1 VE against infection was improved by heterologous mRNA series completion. A 7-8-week interval between first and second doses improved mRNA VE and may be the optimal schedule outside periods of intense epidemic surge. Findings support interchangeability and extended intervals between SARS-CoV-2 vaccine doses, with potential global implications for low-coverage areas and, going forward, for children.

Cov2clusters: genomic clustering of SARS-CoV-2 sequences.

BACKGROUND: The COVID-19 pandemic remains a global public health concern. Advances in sequencing technologies has allowed for high numbers of SARS-CoV-2 whole genome sequence (WGS) data and rapid sharing of sequences through global repositories to enable almost real-time genomic analysis of the pathogen. WGS data has been used previously to group genetically similar viral pathogens to reveal evidence of transmission, including methods that identify distinct clusters on a phylogenetic tree. Identifying clusters of linked cases can aid in the regional surveillance and management of the disease. In this study, we present a novel method for producing stable genomic clusters of SARS-CoV-2 cases, cov2clusters, and compare the accuracy and stability of our approach to previous methods used for phylogenetic clustering using real-world SARS-CoV-2 sequence data obtained from British Columbia, Canada. RESULTS: We found that cov2clusters produced more stable clusters than previously used phylogenetic clustering methods when adding sequence data through time, mimicking an increase in sequence data through the pandemic. Our method also showed high accuracy when predicting epidemiologically informed clusters from sequence data. CONCLUSIONS: Our new approach allows for the identification of stable clusters of SARS-CoV-2 from WGS data. Producing high-resolution SARS-CoV-2 clusters from sequence data alone can a challenge and, where possible, both genomic and epidemiological data should be used in combination.

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.

Understanding virtual patients efficiently and rigorously by combining machine learning with dynamical modelling.

Individual biological organisms are characterized by daunting heterogeneity, which precludes describing or understanding populations of 'patients' with a single mathematical model. Recently, the field of quantitative systems pharmacology (QSP) has adopted the notion of virtual patients (VPs) to cope with this challenge. A typical population of VPs represents the behavior of a heterogeneous patient population with a distribution of parameter values over a mathematical model of fixed structure. Though this notion of VPs is a powerful tool to describe patients' heterogeneity, the analysis and understanding of these VPs present new challenges to systems pharmacologists. Here, using a model of the hypothalamic-pituitary-adrenal axis, we show that an integrated pipeline that combines machine learning (ML) and bifurcation analysis can be used to effectively and efficiently analyse the behaviors observed in populations of VPs. Compared with local sensitivity analyses, ML allows us to capture and analyse the contributions of simultaneous changes of multiple model parameters. Following up with bifurcation analysis, we are able to provide rigorous mechanistic insight regarding the influences of ML-identified parameters on the dynamical system's behaviors. In this work, we illustrate the utility of this pipeline and suggest that its wider adoption will facilitate the use of VPs in the practice of systems pharmacology.

Nucleation of stem cell domains in a bistable activator-inhibitor model of the shoot apical meristem.

Shoot apical meristems (SAMs) give rise to all above-ground tissues of a plant. Expansion of meristematic tissue is derived from the growth and division of stem cells that reside in a central zone of the SAM. This reservoir of stem cells is maintained by expression of a transcription factor WUSCHEL that is responsible for the development of stem cells in the central zone. WUSCHEL expression is self-activating and downregulated by a signaling pathway initiated by CLAVATA proteins, which are upregulated by WUSCHEL. This classic activator-inhibitor network can generate localized patterns of WUSCHEL activity by a Turing instability provided certain constraints on reaction rates and diffusion constants of WUSCHEL and CLAVATA are satisfied, and most existing mathematical models of nucleation and confinement of stem cells in the SAM rely on Turing's mechanism. However, Turing patterns have certain properties that are inconsistent with observed patterns of stem cell differentiation in the SAM. As an alternative mechanism, we propose a model for stem cell confinement based on a bistable-switch in WUSCHEL-CLAVATA interactions. We study the bistable-switch mechanism for pattern formation in a spatially continuous domain and in a discrete cellularized tissue in the presence of a non-uniform field of a rapidly diffusing hormone. By comparing domain formation by Turing and bistable-switch mechanisms in these contexts, we show that bistable switching provides a superior account of nucleation and confinement of the stem cell domain under reasonable assumptions on reaction rates and diffusion constants.

Rapid Increase in SARS-CoV-2 P.1 Lineage Leading to Codominance with B.1.1.7 Lineage, British Columbia, Canada, January-April 2021.

Several severe acute respiratory syndrome coronavirus 2 variants of concern (VOCs) emerged in late 2020; lineage B.1.1.7 initially dominated globally. However, lineages B.1.351 and P.1 represent potentially greater risk for transmission and immune escape. In British Columbia, Canada, B.1.1.7 and B.1.351 were first identified in December 2020 and P.1 in February 2021. We combined quantitative PCR and whole-genome sequencing to assess relative contribution of VOCs in nearly 67,000 infections during the first 16 weeks of 2021 in British Columbia. B.1.1.7 accounted for <10% of screened or sequenced specimens early on, increasing to >50% by week 8. P.1 accounted for <10% until week 10, increased rapidly to peak at week 12, and by week 13 codominated within 10% of rates of B.1.1.7. B.1.351 was a minority throughout. This rapid expansion of P.1 but suppression of B.1.351 expands our understanding of population-level VOC patterns and might provide clues to fitness determinants for emerging VOCs.