ATP regulates the differentiation of mammalian skeletal muscle by activation of a P2X5 receptor on satellite cells

ATP is well known for its role as an intracellular energy source. However, there is increasing awareness of its role as an extracellular messenger molecule (Burnstock, 1997). Although evidence for the presence of receptors for extracellular ATP on skeletal myoblasts was first published in 1983 (Kolb and Wakelam), their physiological function has remained unclear. In this paper we used primary cultures of rat skeletal muscle satellite cells to investigate the role of purinergic signaling in muscle formation. Using immunocytochemistry, RT-PCR, and electrophysiology, we demonstrate that the ionotropic P2X5 receptor is present on satellite cells and that activation of a P2X receptor inhibits proliferation, stimulates expression of markers of muscle cell differentiation, including myogenin, p21, and myosin heavy chain, and increases the rate of myotube formation. Furthermore, we demonstrate that ATP application results in a significant and rapid increase in the phosphorylation of MAPKs, particularly p38, and that inhibition of p38 activity can prevent the effect of ATP on cell number. These results not only demonstrate the existence of a novel regulator of skeletal muscle differentiation, namely ATP, but also a new role for ionotropic P2X receptors in the control of cell fate

Whole genome expression as a quantitative trait

Abstract Surprisingly, whole genome analyses of complex human neurological and psychiatric disorders have revealed that many genetic risk factors are likely to influence gene expression rather than alter protein sequences. Previous analyses of neurological diseases have shown that genetic variability in gene expression levels of deposited proteins influence disease risk. With this background, we have embarked on a comprehensive project to determine the effects of common genetic variability on whole genome gene expression

ensemblQueryR: fast, flexible and high-throughput querying of Ensembl LD API endpoints in R

We present ensemblQueryR, a package providing an R interface to the Ensembl REST API that facilitates flexible, fast, user-friendly and R workflow integrable querying of Ensembl REST API linkage disequilibrium (LD) endpoints, optimised for high-throughput querying. ensemblQueryR achieves this through functions that are intuitive and amenable to custom code integration, use of familiar R object types as inputs and outputs, code optimisation and optional parallelisation functionality. For each LD endpoint, ensemblQueryR provides two functions, permitting both single-query and multi-query modes of operation. The multi-query functions are optimised for large query sizes and provide optional parallelisation to leverage available computational resources and minimise processing time. We demonstrate that ensemblQueryR has improved performance in terms of random access memory (RAM) usage and speed, delivering a 10-fold speed increase over analogous software whilst using a third of the RAM. Finally, ensemblQueryR is near-agnostic to operating system and computational architecture through availability of Docker and singularity images, making this tool widely accessible to the scientific community

ensemblQueryR: fast, flexible and high-throughput querying of Ensembl LD API endpoints in R

We present ensemblQueryR, an R package for querying Ensembl linkage disequilibrium (LD) endpoints. This package is flexible, fast and user-friendly, and optimised for high-throughput querying. ensemblQueryR uses functions that are intuitive and amenable to custom code integration, familiar R object types as inputs and outputs as well as providing parallelisation functionality. For each Ensembl LD endpoint, ensemblQueryR provides two functions, permitting both single- and multi-query modes of operation. The multi-query functions are optimised for large query sizes and provide optional parallelisation to leverage available computational resources and minimise processing time. We demonstrate improved computational performance of ensemblQueryR over an exisiting tool in terms of random access memory (RAM) usage and speed, delivering a 10-fold speed increase whilst using a third of the RAM. Finally, ensemblQueryR is near-agnostic to operating system and computational architecture through Docker and singularity images, making this tool widely accessible to the scientific community

ggtranscript: an R package for the visualization and interpretation of transcript isoforms using ggplot2

MOTIVATION: The advent of long-read sequencing technologies has increased demand for the visualisation and interpretation of transcripts. However, tools that perform such visualizations remain inflexible and lack the ability to easily identify differences between transcript structures. Here, we introduce ggtranscript, an R package that provides a fast and flexible method to visualize and compare transcripts. As a ggplot2 extension, ggtranscript inherits the functionality and familiarity of ggplot2 making it easy to use. AVAILABILITY: ggtranscript is an R package available at https://github.com/dzhang32/ggtranscript (DOI: https://doi.org/10.5281/zenodo.6374061) via an open-source MIT license. Further is available at https://dzhang32.github.io/ggtranscript/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online

Regional genetic correlations highlight relationships between neurodegenerative disease loci and the immune system

Neurodegenerative diseases, including Alzheimer’s and Parkinson’s disease, are devastating complex diseases resulting in physical and psychological burdens on patients and their families. There have been important efforts to understand their genetic basis leading to the identification of disease risk-associated loci involved in several molecular mechanisms, including immune-related pathways. Regional, in contrast to genome-wide, genetic correlations between pairs of immune and neurodegenerative traits have not been comprehensively explored, but could uncover additional immune-mediated risk-associated loci. Here, we systematically assess the role of the immune system in five neurodegenerative diseases by estimating regional genetic correlations between these diseases and immune-cell-derived single-cell expression quantitative trait loci (sc-eQTLs). We also investigate correlations between diseases and protein levels. We observe significant (FDR < 0.01) correlations between sc-eQTLs and neurodegenerative diseases across 151 unique genes, spanning both the innate and adaptive immune systems, across most diseases tested. With Parkinson’s, for instance, RAB7L1 in CD4+ naïve T cells is positively correlated and KANSL1-AS1 is negatively correlated across all adaptive immune cell types. Follow-up colocalization highlight candidate causal risk genes. The outcomes of this study will improve our understanding of the immune component of neurodegeneration, which can warrant repurposing of existing immunotherapies to slow disease progression

Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies

We are building an open-access database of regional human brain expression designed to allow the genome-wide assessment of genetic variability on expression. Array and RNA sequencing technologies make assessment of genome-wide expression possible. Human brain tissue is a challenging source for this work because it can only be obtained several and variable hours post-mortem and after varying agonal states. These variables alter RNA integrity in a complex manner. In this report, we assess the effect of post-mortem delay, agonal state and age on gene expression, and the utility of pH and RNA integrity number as predictors of gene expression as measured on 1266 Affymetrix Exon Arrays. We assessed the accuracy of the array data using QuantiGene, as an independent non-PCR-based method. These quality control parameters will allow database users to assess data accuracy. We report that within the parameters of this study post-mortem delay, agonal state and age have little impact on array quality, array data are robust to variable RNA integrity, and brain pH has only a small effect on array performance. QuantiGene gave very similar expression profiles as array data. This study is the first step in our initiative to make human, regional brain expression freely available

Leveraging omic features with F3UTER enables identification of unannotated 3'UTRs for synaptic genes

There is growing evidence for the importance of 3' untranslated region (3'UTR) dependent regulatory processes. However, our current human 3'UTR catalogue is incomplete. Here, we develop a machine learning-based framework, leveraging both genomic and tissue-specific transcriptomic features to predict previously unannotated 3'UTRs. We identify unannotated 3'UTRs associated with 1,563 genes across 39 human tissues, with the greatest abundance found in the brain. These unannotated 3'UTRs are significantly enriched for RNA binding protein (RBP) motifs and exhibit high human lineage-specificity. We find that brain-specific unannotated 3'UTRs are enriched for the binding motifs of important neuronal RBPs such as TARDBP and RBFOX1, and their associated genes are involved in synaptic function. Our data is shared through an online resource F3UTER ( https://astx.shinyapps.io/F3UTER/ ). Overall, our data improves 3'UTR annotation and provides additional insights into the mRNA-RBP interactome in the human brain, with implications for our understanding of neurological and neurodevelopmental diseases

Genome‐Wide Association Study of Pericardial Fat Area in 28 161 UK Biobank Participants

BACKGROUND: Pericardial adipose tissue (PAT) is the visceral adipose tissue compartment surrounding the heart. Experimental and observational research has suggested that greater PAT deposition might mediate cardiovascular disease, independent of general or subcutaneous adiposity. We characterize the genetic architecture of adiposity‐adjusted PAT and identify causal associations between PAT and adverse cardiac magnetic resonance imaging measures of cardiac structure and function in 28 161 UK Biobank participants. METHODS AND RESULTS: The PAT phenotype was extracted from cardiac magnetic resonance images using an automated image analysis tool previously developed and validated in this cohort. A genome‐wide association study was performed with PAT area set as the phenotype, adjusting for age, sex, and other measures of obesity. Functional mapping and Bayesian colocalization were used to understand the biologic role of identified variants. Mendelian randomization analysis was used to examine potential causal links between genetically determined PAT and cardiac magnetic resonance–derived measures of left ventricular structure and function. We discovered 12 genome‐wide significant variants, with 2 independent sentinel variants (rs6428792, P =4.20×10 −9 and rs11992444, P =1.30×10 −12 ) at 2 distinct genomic loci, that were mapped to 3 potentially causal genes: T‐box transcription factor 15 ( TBX15 ), tryptophanyl tRNA synthetase 2, mitochondrial ( WARS2 ) and early B‐cell factor‐2 ( EBF2 ) through functional annotation. Bayesian colocalization additionally suggested a role of RP4‐712E4.1. Genetically predicted differences in adiposity‐adjusted PAT were causally associated with adverse left ventricular remodeling. CONCLUSIONS: This study provides insights into the genetic architecture determining differential PAT deposition, identifies causal links with left structural and functional parameters, and provides novel data about the pathophysiological importance of adiposity distribution