Human-based modelling and simulations are becoming ubiquitous in biomedical science due to their ability to augment experimental and clinical investigations. Cardiac electrophysiology is one of the most advanced areas, with cardiac modelling and simulation being considered for virtual testing of pharmacological therapies and medical devices. Current models present inconsistencies with experimental data, which limit further progress. In this study, we present the design, development, calibration and independent validation of a human-based ventricular model (ToR-ORd) for simulations of electrophysiology and excitation-contraction coupling, from ionic to whole-organ dynamics, including the electrocardiogram. Validation based on substantial multiscale simulations supports the credibility of the ToR-ORd model under healthy and key disease conditions, as well as drug blockade. In addition, the process uncovers new theoretical insights into the biophysical properties of the L-type calcium current, which are critical for sodium and calcium dynamics. These insights enable the reformulation of L-type calcium current, as well as replacement of the hERG current model.

Progressive tissue remodeling after myocardial infarction (MI) promotes cardiac arrhythmias. This process is well studied in young animals, but little is known about pro-arrhythmic changes in aged animals. Senescent cells accumulate with age and accelerate age-associated diseases. Senescent cells interfere with cardiac function and outcome post-MI with age, but studies have not been performed in larger animals, and the mechanisms are unknown. Specifically, age-associated changes in timecourse of senescence and related changes in inflammation and fibrosis are not well understood. Additionally, the cellular and systemic role of senescence and its inflammatory milieu in influencing arrhythmogenesis with age is not clear, particularly in large animal models with cardiac electrophysiology more similar to humans than previously studied animal models. Here, we investigated the role of senescence in regulating inflammation, fibrosis, and arrhythmogenesis in young and aged infarcted rabbits. Aged rabbits exhibited increased peri-procedural mortality and arrhythmogenic electrophysiological remodeling at the infarct border zone (IBZ) compared to young rabbits. Studies of the aged infarct zone revealed persistent myofibroblast senescence and increased inflammatory signaling over a 12-week timecourse. Senescent IBZ myofibroblasts in aged rabbits appear to be coupled to myocytes, and our computational modeling showed that senescent myofibroblast-cardiomyocyte coupling prolongs action potential duration (APD) and facilitates conduction block permissive of arrhythmias. Aged infarcted human ventricles show levels of senescence consistent with aged rabbits, and senescent myofibroblasts also couple to IBZ myocytes. Our findings suggest that therapeutic interventions targeting senescent cells may mitigate arrhythmias post-MI with age.

Inter-individual variation in gene expression has been shown to be heritable and is often associated with differences in disease susceptibility between individuals. Many studies focused on mapping associations between genetic and gene regulatory variation, yet much less attention has been paid to the evolutionary processes that shape the observed differences in gene regulation between individuals in humans or any other primate. To begin addressing this gap, we performed a comparative analysis of gene expression variability and expression quantitative trait loci (eQTLs) in humans and chimpanzees, using gene expression data from primary heart samples. We found that expression variability in both species is often determined by non-genetic sources, such as cell-type heterogeneity. However, we also provide evidence that inter-individual variation in gene regulation can be genetically controlled, and that the degree of such variability is generally conserved in humans and chimpanzees. In particular, we found a significant overlap of orthologous genes associated with eQTLs in both species. We conclude that gene expression variability in humans and chimpanzees often evolves under similar evolutionary pressures.

After cardiac ischaemia, a prolonged decrease of coronary microvascular perfusion often occurs even after flow is restored in an upstream artery. This 'no-reflow' phenomenon worsens patient prognosis. In the brain, after stroke, a similar post-ischaemic 'no-reflow' has been attributed to capillary constriction by contractile pericytes. We now show that occlusion of a rat coronary artery, followed by reperfusion, blocks 40% of cardiac capillaries and halves perfused blood volume within the affected region. Capillary blockages colocalised strongly with pericytes, where capillary diameter was reduced by 37%. The pericyte relaxant adenosine increased capillary diameter by 21% at pericyte somata, decreased capillary block by 25% and increased perfusion volume by 57%. Thus, cardiac pericytes constrict coronary capillaries and reduce microvascular blood flow after ischaemia, despite re-opening of the culprit artery. Cardiac pericytes are therefore a novel therapeutic target in ischaemic heart disease.

To generate energy efficiently, the cell is uniquely challenged to co-ordinate the abundance of electron transport chain protein subunits expressed from both nuclear and mitochondrial genomes. How an effective stoichiometry of this many constituent subunits is co-ordinated post-transcriptionally remains poorly understood. Here we show that Cerox1, an unusually abundant cytoplasmic long noncoding RNA (lncRNA), modulates the levels of mitochondrial complex I subunit transcripts in a manner that requires binding to microRNA-488-3p. Increased abundance of Cerox1 cooperatively elevates complex I subunit protein abundance and enzymatic activity, decreases reactive oxygen species production, and protects against the complex I inhibitor rotenone. Cerox1 function is conserved across placental mammals: human and mouse orthologues effectively modulate complex I enzymatic activity in mouse and human cells, respectively. Cerox1 is the first lncRNA demonstrated, to our knowledge, to regulate mitochondrial oxidative phosphorylation and, with miR-488-3p, represent novel targets for the modulation of complex I activity.

Myocardial infarction (MI) is accompanied by severe energy deprivation and extensive epigenetic changes. However, how energy metabolism and chromatin modifications are interlinked during MI and heart repair has been poorly explored. Here, we examined the effect of different carbon sources that are involved in the major metabolic pathways of acetyl-CoA synthesis on myocardial infarction and found that elevation of acetyl-CoA by sodium octanoate (8C) significantly improved heart function in ischemia reperfusion (I/R) rats. Mechanistically, 8C reduced I/R injury by promoting histone acetylation which in turn activated the expression of antioxidant genes and inhibited cardiomyocyte (CM) apoptosis. Furthermore, we elucidated that 8C-promoted histone acetylation and heart repair were carried out by metabolic enzyme medium-chain acyl-CoA dehydrogenase (MCAD) and histone acetyltransferase Kat2a, suggesting that 8C dramatically improves cardiac function mainly through metabolic acetyl-CoA-mediated histone acetylation. Therefore, our study uncovers an interlinked metabolic/epigenetic network comprising 8C, acetyl-CoA, MCAD, and Kat2a to combat heart injury.

Chromatin accessibility is modulated in a variety of ways to create open and closed chromatin states, both of which are critical for eukaryotic gene regulation. At the single molecule level, how accessibility is regulated of the chromatin fiber composed of canonical or variant nucleosomes is a fundamental question in the field. Here, we developed a single-molecule tracking method where we could analyze thousands of canonical H3 and centromeric variant nucleosomes imaged by high-speed atomic force microscopy. This approach allowed us to investigate how changes in nucleosome dynamics in vitro inform us about transcriptional potential in vivo. By high-speed atomic force microscopy, we tracked chromatin dynamics in real time and determined the mean square displacement and diffusion constant for the variant centromeric CENP-A nucleosome. Furthermore, we found that an essential kinetochore protein CENP-C reduces the diffusion constant and mobility of centromeric nucleosomes along the chromatin fiber. We subsequently interrogated how CENP-C modulates CENP-A chromatin dynamics in vivo. Overexpressing CENP-C resulted in reduced centromeric transcription and impaired loading of new CENP-A molecules. From these data, we speculate that factors altering nucleosome mobility in vitro, also correspondingly alter transcription in vivo. Subsequently, we propose a model in which variant nucleosomes encode their own diffusion kinetics and mobility, and where binding partners can suppress or enhance nucleosome mobility.

We recently found that nucleosomes directly block access of CRISPR/Cas9 to DNA (Horlbeck et al., 2016). Here, we build on this observation with a comprehensive algorithm that incorporates chromatin, position, and sequence features to accurately predict highly effective single guide RNAs (sgRNAs) for targeting nuclease-dead Cas9-mediated transcriptional repression (CRISPRi) and activation (CRISPRa). We use this algorithm to design next-generation genome-scale CRISPRi and CRISPRa libraries targeting human and mouse genomes. A CRISPRi screen for essential genes in K562 cells demonstrates that the large majority of sgRNAs are highly active. We also find CRISPRi does not exhibit any detectable non-specific toxicity recently observed with CRISPR nuclease approaches. Precision-recall analysis shows that we detect over 90% of essential genes with minimal false positives using a compact 5 sgRNA/gene library. Our results establish CRISPRi and CRISPRa as premier tools for loss- or gain-of-function studies and provide a general strategy for identifying Cas9 target sites.

Cancer survivors suffer from progressive frailty, multimorbidity, and premature morbidity. We hypothesise that therapy-induced senescence and senescence progression via bystander effects are significant causes of this premature ageing phenotype. Accordingly, the study addresses the question whether a short anti-senescence intervention is able to block progression of radiation-induced frailty and disability in a pre-clinical setting. Male mice were sublethally irradiated at 5 months of age and treated (or not) with either a senolytic drug (Navitoclax or dasatinib + quercetin) for 10 days or with the senostatic metformin for 10 weeks. Follow-up was for 1 year. Treatments commencing within a month after irradiation effectively reduced frailty progression (p<0.05) and improved muscle (p<0.01) and liver (p<0.05) function as well as short-term memory (p<0.05) until advanced age with no need for repeated interventions. Senolytic interventions that started late, after radiation-induced premature frailty was manifest, still had beneficial effects on frailty (p<0.05) and short-term memory (p<0.05). Metformin was similarly effective as senolytics. At therapeutically achievable concentrations, metformin acted as a senostatic neither via inhibition of mitochondrial complex I, nor via improvement of mitophagy or mitochondrial function, but by reducing non-mitochondrial reactive oxygen species production via NADPH oxidase 4 inhibition in senescent cells. Our study suggests that the progression of adverse long-term health and quality-of-life effects of radiation exposure, as experienced by cancer survivors, might be rescued by short-term adjuvant anti-senescence interventions.

White adipose tissue (WAT) remodeling is dictated by coordinated interactions between adipocytes and resident stromal-vascular cells; however, the functional heterogeneity of adipose stromal cells has remained unresolved. We combined single-cell RNA-sequencing and FACS to identify and isolate functionally distinct subpopulations of PDGFRβ+ stromal cells within visceral WAT of adult mice. LY6C- CD9- PDGFRβ+ cells represent highly adipogenic visceral adipocyte precursor cells ('APCs'), whereas LY6C+ PDGFRβ+ cells represent fibro-inflammatory progenitors ('FIPs'). FIPs lack adipogenic capacity, display pro-fibrogenic/pro-inflammatory phenotypes, and can exert an anti-adipogenic effect on APCs. The pro-inflammatory phenotype of PDGFRβ+ cells is regulated, at least in part, by NR4A nuclear receptors. These data highlight the functional heterogeneity of visceral WAT perivascular cells, and provide insight into potential cell-cell interactions impacting adipogenesis and inflammation. These improved strategies to isolate FIPs and APCs from visceral WAT will facilitate the study of physiological WAT remodeling and mechanisms leading to metabolic dysfunction.

Respiratory chemoreceptors regulate breathing in response to changes in tissue CO2/H+. Blood flow is a fundamental determinant of tissue CO2/H+, yet little is known regarding how regulation of vascular tone in chemoreceptor regions contributes to respiratory behavior. Previously, we showed in rat that CO2/H+-vasoconstriction in the retrotrapezoid nucleus (RTN) supports chemoreception by a purinergic-dependent mechanism (Hawkins et al., 2017). Here, we show in mice that CO2/H+ dilates arterioles in other chemoreceptor regions, thus demonstrating CO2/H+ vascular reactivity in the RTN is unique. We also identify P2Y2 receptors in RTN smooth muscle cells as the substrate responsible for this response. Specifically, pharmacological blockade or genetic deletion of P2Y2 from smooth muscle cells blunted the ventilatory response to CO2, and re-expression of P2Y2 receptors only in RTN smooth muscle cells fully rescued the CO2/H+ chemoreflex. These results identify P2Y2 receptors in RTN smooth muscle cells as requisite determinants of respiratory chemoreception.

Klebsiella pneumoniae is a respiratory, blood, liver, and bladder pathogen of significant clinical concern. We show that the adaptor protein, SKAP2, is required for protection against K. pneumoniae (ATCC 43816) pulmonary infections. Skap2-/- mice had 100-fold higher bacterial burden when compared to wild-type and burden was controlled by SKAP2 expression in innate immune cells. Skap2-/- neutrophils and monocytes were present in infected lungs, and the neutrophils degranulated normally in response to K. pneumoniae infection in mice; however, K. pneumoniae-stimulated reactive oxygen species (ROS) production in vitro was abolished. K. pneumoniae-induced neutrophil ROS response required the activity of SFKs, Syk, Btk, PLCγ2, and PKC. The loss of SKAP2 significantly hindered the K. pneumoniae-induced phosphorylation of SFKs, Syk, and Pyk2 implicating SKAP2 as proximal to their activation in pathogen-signaling pathways. In conclusion, SKAP2-dependent signaling in neutrophils is essential for K. pneumoniae-activated ROS production and for promoting bacterial clearance during infection.

Translational readthrough gives rise to low abundance proteins with C-terminal extensions beyond the stop codon. To identify functional translational readthrough, we estimated the readthrough propensity (RTP) of all stop codon contexts of the human genome by a new regression model in silico, identified a nucleotide consensus motif for high RTP by using this model, and analyzed all readthrough extensions in silico with a new predictor for peroxisomal targeting signal type 1 (PTS1). Lactate dehydrogenase B (LDHB) showed the highest combined RTP and PTS1 probability. Experimentally we show that at least 1.6% of the total cellular LDHB is targeted to the peroxisome by a conserved hidden PTS1. The readthrough-extended lactate dehydrogenase subunit LDHBx can also co-import LDHA, the other LDH subunit, into peroxisomes. Peroxisomal LDH is conserved in mammals and likely contributes to redox equivalent regeneration in peroxisomes.

Currently there is great interest in targeting mitochondrial oxidative phosphorylation (OXPHOS) in cancer. However, notwithstanding the targeting of mutant dehydrogenases, nearly all hopeful 'mito-therapeutics' cannot discriminate cancerous from non-cancerous OXPHOS and thus suffer from a limited therapeutic index. Using acute myeloid leukemia (AML) as a model, herein, we leveraged an in-house diagnostic biochemical workflow to identify 'actionable' bioenergetic vulnerabilities intrinsic to cancerous mitochondria. Consistent with prior reports, AML growth and proliferation was associated with a hyper-metabolic phenotype which included increases in basal and maximal respiration. However, despite having nearly 2-fold more mitochondria per cell, clonally expanding hematopoietic stem cells, leukemic blasts, as well as chemoresistant AML were all consistently hallmarked by intrinsic OXPHOS limitations. Remarkably, by performing experiments across a physiological span of ATP free energy, we provide direct evidence that leukemic mitochondria are particularly poised to consume ATP. Relevant to AML biology, acute restoration of oxidative ATP synthesis proved highly cytotoxic to leukemic blasts, suggesting that active OXPHOS repression supports aggressive disease dissemination in AML. Together, these findings argue against ATP being the primary output of leukemic mitochondria and provide proof-of-principle that restoring, rather than disrupting, OXPHOS may represent an untapped therapeutic avenue for combatting hematological malignancy and chemoresistance.