In Vivo Target Gene Activation via CRISPR/Cas9-Mediated Trans-epigenetic Modulation.

Current genome-editing systems generally rely on inducing DNA double-strand breaks (DSBs). This may limit their utility in clinical therapies, as unwanted mutations caused by DSBs can have deleterious effects. CRISPR/Cas9 system has recently been repurposed to enable target gene activation, allowing regulation of endogenous gene expression without creating DSBs. However, in vivo implementation of this gain-of-function system has proven difficult. Here, we report a robust system for in vivo activation of endogenous target genes through trans-epigenetic remodeling. The system relies on recruitment of Cas9 and transcriptional activation complexes to target loci by modified single guide RNAs. As proof-of-concept, we used this technology to treat mouse models of diabetes, muscular dystrophy, and acute kidney disease. Results demonstrate that CRISPR/Cas9-mediated target gene activation can be achieved in vivo, leading to measurable phenotypes and amelioration of disease symptoms. This establishes new avenues for developing targeted epigenetic therapies against human diseases. VIDEO ABSTRACT.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project.

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Current Status of and Perspectives on the Application of Marmosets in Neurobiology.

The common marmoset (Callithrix jacchus), a small New World primate, is receiving substantial attention in the neuroscience and biomedical science fields because its anatomical features, functional and behavioral characteristics, and reproductive features and its amenability to available genetic modification technologies make it an attractive experimental subject. In this review, I outline the progress of marmoset neuroscience research and summarize both the current status (opportunities and limitations) of and the future perspectives on the application of marmosets in neuroscience and disease modeling.

Mitochondrial abnormalities contribute to muscle weakness in a Dnajb6 deficient zebrafish model.

Mutations in DNAJB6 are a well-established cause of limb girdle muscular dystrophy type D1 (LGMD D1). Patients with LGMD D1 develop progressive muscle weakness with histology showing fibre damage, autophagic vacuoles, and aggregates. Whilst there are many reports of LGMD D1 patients, the role of DNAJB6 in the muscle is still unclear. In this study, we developed a loss of function zebrafish model in order to investigate the role of Dnajb6. Using a double dnajb6a and dnajb6b mutant model, we show that loss of Dnajb6 leads to a late onset muscle weakness. Interestingly, we find that adult fish lacking Dnajb6 do not have autophagy or myofibril defects, however, they do show mitochondrial changes and damage. This study demonstrates that loss of Dnajb6 causes mitochondrial defects and suggests that this contributes to muscle weakness in LGMD D1. These findings expand our knowledge of the role of Dnajb6 in the muscle and provides a model to screen novel therapies for LGMD D1.

Genetic disruption of mammalian endoplasmic reticulum-associated protein degradation: Human phenotypes and animal and cellular disease models.

Endoplasmic reticulum-associated protein degradation (ERAD) is a stringent quality control mechanism through which misfolded, unassembled and some native proteins are targeted for degradation to maintain appropriate cellular and organelle homeostasis. Several in vitro and in vivo ERAD-related studies have provided mechanistic insights into ERAD pathway activation and its consequent events; however, a majority of these have investigated the effect of ERAD substrates and their consequent diseases affecting the degradation process. In this review, we present all reported human single-gene disorders caused by genetic variation in genes that encode ERAD components rather than their substrates. Additionally, after extensive literature survey, we present various genetically manipulated higher cellular and mammalian animal models that lack specific components involved in various stages of the ERAD pathway.

Distinctive In Vitro Phenotypes in iPSC-Derived Neurons From Patients With Gain- and Loss-of-Function SCN2A Developmental and Epileptic Encephalopathy.

SCN2A encodes NaV1.2, an excitatory neuron voltage-gated sodium channel and a major monogenic cause of neurodevelopmental disorders, including developmental and epileptic encephalopathies (DEE) and autism. Clinical presentation and pharmocosensitivity vary with the nature of SCN2A variant dysfunction and can be divided into gain-of-function (GoF) cases with pre- or peri-natal seizures and loss-of-function (LoF) patients typically having infantile spasms after 6 months of age. We established and assessed patient induced pluripotent stem cell (iPSC) - derived neuronal models for two recurrent SCN2A DEE variants with GoF R1882Q and LoF R853Q associated with early- and late-onset DEE, respectively. Two male patient-derived iPSC isogenic pairs were differentiated using Neurogenin-2 overexpression yielding populations of cortical-like glutamatergic neurons. Functional properties were assessed using patch clamp and multielectrode array recordings and transcriptomic profiles obtained with total mRNA sequencing after 2-4 weeks in culture. At 3 weeks of differentiation, increased neuronal activity at cellular and network levels was observed for R1882Q iPSC-derived neurons. In contrast, R853Q neurons showed only subtle changes in excitability after 4 weeks and an overall reduced network activity after 7 weeks in vitro. Consistent with the reported efficacy in some GoF SCN2A patients, phenytoin (sodium channel blocker) reduced the excitability of neurons to the control levels in R1882Q neuronal cultures. Transcriptomic alterations in neurons were detected for each variant and convergent pathways suggested potential shared mechanisms underlying SCN2A DEE. In summary, patient iPSC-derived neuronal models of SCN2A GoF and LoF pathogenic variants causing DEE show specific functional and transcriptomic in vitro phenotypes.

The 36th International Mammalian Genome Conference: A scientific gathering under the cherry blossoms in Tsukuba.

The 36th International Mammalian Genome Conference (IMGC) was held in a hybrid format at the Tsukuba International Congress Center in Tsukuba, Ibaraki, Japan, for 4 days from March 28 to 31, 2023. This international conference on functional genomics of mouse, human, and other mammalian species attracted 246 participants in total, of which 129 were from outside Japan, including Europe, the United States and Asia, and 117 participants were from Japan. The conference included three technical workshops, keynote lectures by domestic researchers, commemorative lectures for the conference awards, 57 oral presentations, and 97 poster presentations. The event was a great success. Topics included the establishment and analysis of disease models using genetically engineered or spontaneous mutant mice, systems genetic analysis using mouse strains such as wild-derived mice and recombinant inbred mouse strains, infectious diseases, immunology, and epigenetics. In addition, as a joint program, a two-day RIKEN Symposium was held, and active discussions continued over the four-day period. Also, there was a trainee symposium, in which young researchers were encouraged to participate, and excellent papers were selected as oral presentations in the main session.

Recent advances regarding the potential roles of invariant natural killer T cells in cardiovascular diseases with immunological and inflammatory backgrounds.

Invariant natural killer T (iNKT) cells, which bear αß-type T-cell antigen-receptors (TCRs), recognize glycolipid antigens in a cluster of differentiation 1d (CD1d)-restricted manner. Regarding these cells, the unique modes of thymic selection and maturation elucidate innateness, irrespective of them also being members of the adaptive immune system as a T-cell. iNKT cells develop and differentiate into NKT1 [interferon γ (IFN-γ)-producing], NKT2 [interleukin 4 (IL-4)/IL-13-producing], or NKT17 (IL-17-producing) subsets in the thymus. After egress, NKT10 (IL-10-producing), follicular helper NKT (NKTfh; IL-21-producing), and regulatory NKT (NKTreg) subsets emerge following stimulation in the periphery. Moreover, iNKT cells have been shown to possess several physiological or pathological roles. iNKT cells exhibit dual alleviating or aggravating roles in experimentally induced immune and/or inflammatory diseases in mice. These findings indicate that the modulation of iNKT cells can be employed for therapeutic use or prevention of human diseases. In this review, we discuss the potential roles of iNKT cells in the development of immune/inflammatory diseases of the cardiovascular system, with emphasis on atherosclerosis, aortic aneurysms, and cardiac remodeling.

Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision.

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

Atrial proteomic profiling reveals a switch towards profibrotic gene expression program in CREM-IbΔC-X mice with persistent atrial fibrillation.

BACKGROUND: Overexpression of the CREM (cAMP response element-binding modulator) isoform CREM-IbΔC-X in transgenic mice (CREM-Tg) causes the age-dependent development of spontaneous AF. PURPOSE: To identify key proteome signatures and biological processes accompanying the development of persistent AF through integrated proteomics and bioinformatics analysis. METHODS: Atrial tissue samples from three CREM-Tg mice and three wild-type littermates were subjected to unbiased mass spectrometry-based quantitative proteomics, differential expression and pathway enrichment analysis, and protein-protein interaction (PPI) network analysis. RESULTS: A total of 98 differentially expressed proteins were identified. Gene ontology analysis revealed enrichment for biological processes regulating actin cytoskeleton organization and extracellular matrix (ECM) dynamics. Changes in ITGAV, FBLN5, and LCP1 were identified as being relevant to atrial fibrosis and structural based on expression changes, co-expression patterns, and PPI network analysis. Comparative analysis with previously published datasets revealed a shift in protein expression patterns from ion-channel and metabolic regulators in young CREM-Tg mice to profibrotic remodeling factors in older CREM-Tg mice. Furthermore, older CREM-Tg mice exhibited protein expression patterns reminiscent of those seen in humans with persistent AF. CONCLUSIONS: This study uncovered distinct temporal changes in atrial protein expression patterns with age in CREM-Tg mice consistent with the progressive evolution of AF. Future studies into the role of the key differentially abundant proteins identified in this study in AF progression may open new therapeutic avenues to control atrial fibrosis and substrate development in AF.

Expanded in vivo substrate profile of the yeast N-terminal acetyltransferase NatC.

N-terminal acetylation is a conserved protein modification among eukaryotes. The yeast Saccharomyces cerevisiae is a valuable model system for studying this modification. The bulk of protein N-terminal acetylation in S. cerevisiae is catalyzed by the N-terminal acetyltransferases NatA, NatB, and NatC. Thus far, proteome-wide identification of the in vivo protein substrates of yeast NatA and NatB has been performed by N-terminomics. Here, we used S. cerevisiae deleted for the NatC catalytic subunit Naa30 and identified 57 yeast NatC substrates by N-terminal combined fractional diagonal chromatography analysis. Interestingly, in addition to the canonical N-termini starting with ML, MI, MF, and MW, yeast NatC substrates also included MY, MK, MM, MA, MV, and MS. However, for some of these substrate types, such as MY, MK, MV, and MS, we also uncovered (residual) non-NatC NAT activity, most likely due to the previously established redundancy between yeast NatC and NatE/Naa50. Thus, we have revealed a complex interplay between different NATs in targeting methionine-starting N-termini in yeast. Furthermore, our results showed that ectopic expression of human NAA30 rescued known NatC phenotypes in naa30Δ yeast, as well as partially restored the yeast NatC Nt-acetylome. Thus, we demonstrate an evolutionary conservation of NatC from yeast to human thereby underpinning future disease models to study pathogenic NAA30 variants. Overall, this work offers increased biochemical and functional insights into NatC-mediated N-terminal acetylation and provides a basis for future work to pinpoint the specific molecular mechanisms that link the lack of NatC-mediated N-terminal acetylation to phenotypes of NatC deletion yeast.

Rare and de novo coding variants in chromodomain genes in Chiari I malformation.

Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10-10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10-10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10-6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10-9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

A patient-based medaka alg2 mutant as a model for hypo-N-glycosylation.

Defects in the evolutionarily conserved protein-glycosylation machinery during embryonic development are often fatal. Consequently, congenital disorders of glycosylation (CDG) in human are rare. We modelled a putative hypomorphic mutation described in an alpha-1,3/1,6-mannosyltransferase (ALG2) index patient (ALG2-CDG) to address the developmental consequences in the teleost medaka (Oryzias latipes). We observed specific, multisystemic, late-onset phenotypes, closely resembling the patient's syndrome, prominently in the facial skeleton and in neuronal tissue. Molecularly, we detected reduced levels of N-glycans in medaka and in the patient's fibroblasts. This hypo-N-glycosylation prominently affected protein abundance. Proteins of the basic glycosylation and glycoprotein-processing machinery were over-represented in a compensatory response, highlighting the regulatory topology of the network. Proteins of the retinal phototransduction machinery, conversely, were massively under-represented in the alg2 model. These deficiencies relate to a specific failure to maintain rod photoreceptors, resulting in retinitis pigmentosa characterized by the progressive loss of these photoreceptors. Our work has explored only the tip of the iceberg of N-glycosylation-sensitive proteins, the function of which specifically impacts on cells, tissues and organs. Taking advantage of the well-described human mutation has allowed the complex interplay of N-glycosylated proteins and their contribution to development and disease to be addressed.

A new model of axon degeneration in the mouse optic nerve using repeat intraocular pressure challenge.

We characterize a new experimental model for inducing retinal ganglion cell (RGC) dysfunction and degeneration in mice. C57BL/6J mice were subjected to two acute periods of intraocular pressure (IOP) elevation (50 mmHg for 30 min) by cannulation of the anterior chamber. We used full-field electroretinography and visual evoked potentials (VEPs) to measure subsequent changes in retina and optic nerve function, and histochemical techniques to assess RGC survival and optic nerve structure. In 12 month old mice, a single IOP challenge caused loss and subsequent recovery of RGC function over the following 28 days with minimal cell death and no observed axonal damage. A second identical IOP challenge resulted in persistent RGC dysfunction and significant (36%) loss of RGC somas. This was accompanied by a 16.7% delay in the latency and a 27.6% decrease in the amplitude of the VEP. Severe axonal damage was seen histologically with enlargement of axons, myelin disruption, reduced axon density, and the presence of glial scarring. In contrast, younger 3 month old mice when exposed to a single or repeat IOP challenge showed quicker RGC functional recovery after a single challenge and full functional recovery after a repeat challenge with no detectable optic nerve dysfunction. These data demonstrate a highly reproducible and minimally invasive method for inducing RGC degeneration and axonal damage in mice. Resilience of the optic nerve to damage is highly dependent on animal age. The time-defined nature of functional versus structural loss seen in this model stands to facilitate investigation of neuroglial responses in the retina after IOP injury and the associated evaluation of neuroprotective treatment strategies. Further, the model may be used to investigate the impact of aging and the cellular switch between neurorecovery and neurodegeneration.

More than a small adult brain: Lessons from chemotherapy-induced cognitive impairment for modelling paediatric brain disorders.

Childhood is recognised as a period of immense physical and emotional development, and this, in part, is driven by underlying neurophysiological transformations. These neurodevelopmental processes are unique to the paediatric brain and are facilitated by augmented rates of neuroplasticity and expanded neural stem cell populations within neurogenic niches. However, given the immaturity of the developing central nervous system, innate protective mechanisms such as neuroimmune and antioxidant responses are functionally naïve which results in periods of heightened sensitivity to neurotoxic insult. This is highly relevant in the context of paediatric cancer, and in particular, the neurocognitive symptoms associated with treatment, such as surgery, radio- and chemotherapy. The vulnerability of the developing brain may increase susceptibility to damage and persistent symptomology, aligning with reports of more severe neurocognitive dysfunction in children compared to adults. It is therefore surprising, given this intensified neurocognitive burden, that most of the pre-clinical, mechanistic research focuses exclusively on adult populations and extrapolates findings to paediatric cohorts. Given this dearth of age-specific research, throughout this review we will draw comparisons with neurodevelopmental disorders which share comparable pathways to cancer treatment related side-effects. Furthermore, we will examine the unique nuances of the paediatric brain along with the somatic systems which influence neurological function. In doing so, we will highlight the importance of developing in vitro and in vivo paediatric disease models to produce age-specific discovery and clinically translatable research.

PRIMCAT-CRC: a data-driven whole disease simulation model for predicting the number of colorectal cancer patients by stage and treatment line in Australia.

OBJECTIVES: Effective healthcare planning, resource allocation, and budgeting require accurate predictions of the number of patients needing treatment at specific cancer stages and treatment lines. The PRedicting the population health economic IMpact of current and new Cancer Treatments Colorectal Cancer (PRIMCAT-CRC) simulation model was developed to meet this requirement for all CRC stages and relevant molecular profiles in Australia. METHODS: Real-world data was used to estimate treatment utilisation and time-to-event distributions. This populated a discrete-event simulation, projecting the number of patients receiving treatment across all disease stages and treatment lines for CRC and forecasting the number of patients likely to utilise future treatments. Illustrative analyses were undertaken, estimating treatments across disease stages and treatment lines over a 5-year period (2022-2026). We demonstrated the model's applicability through a case study introducing pembrolizumab as a first-line treatment for mismatch-repair deficient stage IV. RESULTS: Clinical registry data from 7,163 patients informed the model. The model forecasts 15,738 incident and 2,821 prevalent cases requiring treatment in 2022, rising to 15,921 and 2,871 respectively by 2026. Projections show that over 2022-2026, there will be a total of 116,752 treatments initiated, with 43% intended for stage IV disease. The introduction of pembrolizumab is projected for 706 patients annually, totalling 3,530 individuals starting treatment with pembrolizumab over the forecasted period, without significantly altering downstream utilisation of subsequent treatments. CONCLUSIONS: PRIMCAT-CRC is a versatile tool that can be used to estimate the eligible patient populations for novel cancer therapies, thereby reducing uncertainty for policymakers in decisions to publicly reimburse new treatments.

Identification of tumour-infiltrating myeloid subsets associated with overall survival in lung squamous cell carcinoma.

Lung squamous cell carcinoma (LUSC) is a primary subtype of lung cancer with limited therapeutic options and poor prognosis, and tumour-infiltrating myeloid cells (TIMs) are key regulators of LUSC. However, the correlation between the abundance of TIM subtypes and clinical outcomes of LUSC remains unexplored. This study aimed to develop and validate a prognostic model for low- and high-risk patients with LUSC based on myeloid cell microenvironments. TIM markers in the tumoural (T) and stromal (S) regions were quantified using immunohistochemistry for 502 LUSC patients. L1-penalized Cox regression was used to develop a myeloid survival score (MSS) model based on the training cohort, followed by validation in distinct cohorts from multiple centres. RNA sequencing and immunostaining were used to examine the mechanisms of myeloid cells in LUSC progression and predict potential drug targets and therapeutic agents. Of the 12 myeloid markers, CD163T, CD163S, and S100A12T were highly associated with overall survival (OS) in LUSC patients. The MSS of the three myeloid signatures accurately categorized LUSC patients into risk categories, with an observable difference in OS between the training and validation cohorts. Tumours with high MSS were associated with enhanced antioxidative ability and hedgehog signalling and a shift to a more pro-tumorigenic microenvironment, accompanied by a reduced tumour cell immunogenicity and increased CD8+ T cell exhaustion patterns. Additionally, in high-risk patients, potential drug targets and compounds regulating hedgehog signalling were identified. Our study provides the first prognostic myeloid signature for LUSC, which may help advance precision medicine. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Determining the toxicological effects of indoor air pollution on both a healthy and an inflammatory-comprised model of the alveolar epithelial barrier in vitro.

Exposure to indoor air pollutants (IAP) has increased recently, with people spending more time indoors (i.e. homes, offices, schools and transportation). Increased exposures of IAP on a healthy population are poorly understood, and those with allergic respiratory conditions even less so. The objective of this study, therefore, was to implement a well-characterised in vitro model of the human alveolar epithelial barrier (A549 + PMA differentiated THP-1 incubated with and without IL-13, IL-5 and IL-4) to determine the effects of a standardised indoor particulate (NIST 2583) on both a healthy lung model and one modelling a type-II (stimulated with IL-13, IL-5 and IL-4) inflammatory response (such as asthma).Using concentrations from the literature, and an environmentally appropriate exposure we investigated 232, 464 and 608ng/cm2 of NIST 2583 respectively. Membrane integrity (blue dextran), viability (trypan blue), genotoxicity (micronucleus (Mn) assay) and (pro-)/(anti-)inflammatory effects (IL-6, IL-8, IL-33, IL-10) were then assessed 24 h post exposure to both models. Models were exposed using a physiologically relevant aerosolisation method (VitroCell Cloud 12 exposure system).No changes in Mn frequency or membrane integrity in either model were noted when exposed to any of the tested concentrations of NIST 2583. A significant decrease (p < 0.05) in cell viability at the highest concentration was observed in the healthy model. Whilst cell viability in the "inflamed" model was decreased at the lower concentrations (significantly (p < 0.05) after 464ng/cm2). A significant reduction (p < 0.05) in IL-10 and a significant increase in IL-33 was seen after 24 h exposure to NIST 2583 (464, 608ng/cm2) in the "inflamed" model.Collectively, the results indicate the potential for IAP to cause the onset of a type II response as well as exacerbating pre-existing allergic conditions. Furthermore, the data imposes the importance of considering unhealthy individuals when investigating the potential health effects of IAP. It also highlights that even in a healthy population these particles have the potential to induce this type II response and initiate an immune response following exposure to IAP.

Vaccination for communicable endemic diseases: optimal allocation of initial and booster vaccine doses.

For some communicable endemic diseases (e.g., influenza, COVID-19), vaccination is an effective means of preventing the spread of infection and reducing mortality, but must be augmented over time with vaccine booster doses. We consider the problem of optimally allocating a limited supply of vaccines over time between different subgroups of a population and between initial versus booster vaccine doses, allowing for multiple booster doses. We first consider an SIS model with interacting population groups and four different objectives: those of minimizing cumulative infections, deaths, life years lost, or quality-adjusted life years lost due to death. We solve the problem sequentially: for each time period, we approximate the system dynamics using Taylor series expansions, and reduce the problem to a piecewise linear convex optimization problem for which we derive intuitive closed-form solutions. We then extend the analysis to the case of an SEIS model. In both cases vaccines are allocated to groups based on their priority order until the vaccine supply is exhausted. Numerical simulations show that our analytical solutions achieve results that are close to optimal with objective function values significantly better than would be obtained using simple allocation rules such as allocation proportional to population group size. In addition to being accurate and interpretable, the solutions are easy to implement in practice. Interpretable models are particularly important in public health decision making.

Impact assessment of self-medication on COVID-19 prevalence in Gauteng, South Africa, using an age-structured disease transmission modelling framework.

OBJECTIVE: To assess the impact of self-medication on the transmission dynamics of COVID-19 across different age groups, examine the interplay of vaccination and self-medication in disease spread, and identify the age group most prone to self-medication. METHODS: We developed an age-structured compartmentalized epidemiological model to track the early dynamics of COVID-19. Age-structured data from the Government of Gauteng, encompassing the reported cumulative number of cases and daily confirmed cases, were used to calibrate the model through a Markov Chain Monte Carlo (MCMC) framework. Subsequently, uncertainty and sensitivity analyses were conducted on the model parameters. RESULTS: We found that self-medication is predominant among the age group 15-64 (74.52%), followed by the age group 0-14 (34.02%), and then the age group 65+ (11.41%). The mean values of the basic reproduction number, the size of the first epidemic peak (the highest magnitude of the disease), and the time of the first epidemic peak (when the first highest magnitude occurs) are 4.16499, 241,715 cases, and 190.376 days, respectively. Moreover, we observed that self-medication among individuals aged 15-64 results in the highest spreading rate of COVID-19 at the onset of the outbreak and has the greatest impact on the first epidemic peak and its timing. CONCLUSION: Studies aiming to understand the dynamics of diseases in areas prone to self-medication should account for this practice. There is a need for a campaign against COVID-19-related self-medication, specifically targeting the active population (ages 15-64).