Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution

  1. Korin Sahinyan
  2. Darren M Blackburn
  3. Marie-Michelle Simon
  4. Felicia Lazure
  5. Tony Kwan
  6. Guillaume Bourque
  7. Vahab D Soleimani  Is a corresponding author
  1. McGill University, Canada
  2. McGill University and Genome Quebec Innovation Centre, Canada
  3. McGill University and Génome Québec Innovation Centre, Canada

Abstract

Myofibers are the main components of skeletal muscle, which is the largest tissue in the body. Myofibers are highly adaptive and can be altered under different biological and disease conditions. Therefore, transcriptional and epigenetic studies on myofibers are crucial to discover how chromatin alterations occur in the skeletal muscle under different conditions. However, due to the heterogenous nature of skeletal muscle, studying myofibers in isolation proves to be a challenging task. Single cell sequencing has permitted the study of the epigenome of isolated myonuclei. While this provides sequencing with high dimensionality, the sequencing depth is lacking, which makes comparisons between different biological conditions difficult. Here we report the first implementation of single myofiber ATAC-Seq, which allows for the sequencing of an individual myofiber at a depth sufficient for peak calling and for comparative analysis of chromatin accessibility under various physiological and disease conditions. Application of this technique revealed significant differences in chromatin accessibility between resting and regenerating myofibers, as well as between myofibers from a mouse model of Duchenne Muscular Dystrophy (mdx) and wild type (WT) counterparts. This technique can lead to a wide application in the identification of chromatin regulatory elements and epigenetic mechanisms in muscle fibers during development and in muscle-wasting diseases.

Data availability

The data discussed in this study have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession numbers GSE173676 and GSE171534

The following data sets were generated

Article and author information

Author details

  1. Korin Sahinyan

    Department of Human Genetics, McGill University, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
  2. Darren M Blackburn

    Department of Human Genetics, McGill University, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
  3. Marie-Michelle Simon

    McGill University and Genome Quebec Innovation Centre, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
  4. Felicia Lazure

    Department of Human Genetics, McGill University, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
  5. Tony Kwan

    McGill University and Génome Québec Innovation Centre, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
  6. Guillaume Bourque

    Department of Human Genetics, McGill University, Montreal, Canada
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3933-9656
  7. Vahab D Soleimani

    Department of Human Genetics, McGill University, Montreal, Canada
    For correspondence
    vahab.soleimani@mcgill.ca
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2154-4894

Funding

Natural Sciences and Engineering Research Council of Canada

  • Vahab D Soleimani

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. YM Dennis Lo, The Chinese University of Hong Kong, Hong Kong

Ethics

Animal experimentation: All procedures that were performed on animals were approved by the McGill University Animal Care Committee (UACC), protocol #7512

Version history

  1. Preprint posted: June 15, 2021 (view preprint)
  2. Received: August 4, 2021
  3. Accepted: February 9, 2022
  4. Accepted Manuscript published: February 21, 2022 (version 1)
  5. Version of Record published: March 7, 2022 (version 2)

Copyright

© 2022, Sahinyan et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 5,399
    views
  • 489
    downloads
  • 5
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Korin Sahinyan
  2. Darren M Blackburn
  3. Marie-Michelle Simon
  4. Felicia Lazure
  5. Tony Kwan
  6. Guillaume Bourque
  7. Vahab D Soleimani
(2022)
Application of ATAC-Seq for genome-wide analysis of the chromatin state at single myofiber resolution
eLife 11:e72792.
https://doi.org/10.7554/eLife.72792

Share this article

https://doi.org/10.7554/eLife.72792

Further reading

    1. Genetics and Genomics
    2. Immunology and Inflammation
    Jean-David Larouche, Céline M Laumont ... Claude Perreault
    Research Article

    Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

    1. Genetics and Genomics
    Pianpian Zhao, Zhifeng Sheng ... Hou-Feng Zheng
    Research Article

    The ‘diabetic bone paradox’ suggested that type 2 diabetes (T2D) patients would have higher areal bone mineral density (BMD) but higher fracture risk than individuals without T2D. In this study, we found that the genetically predicted T2D was associated with higher BMD and lower risk of fracture in both weighted genetic risk score (wGRS) and two-sample Mendelian randomization (MR) analyses. We also identified ten genomic loci shared between T2D and fracture, with the top signal at SNP rs4580892 in the intron of gene RSPO3. And the higher expression in adipose subcutaneous and higher protein level in plasma of RSPO3 were associated with increased risk of T2D, but decreased risk of fracture. In the prospective study, T2D was observed to be associated with higher risk of fracture, but BMI mediated 30.2% of the protective effect. However, when stratified by the T2D-related risk factors for fracture, we observed that the effect of T2D on the risk of fracture decreased when the number of T2D-related risk factors decreased, and the association became non-significant if the T2D patients carried none of the risk factors. In conclusion, the genetically determined T2D might not be associated with higher risk of fracture. And the shared genetic architecture between T2D and fracture suggested a top signal around RSPO3 gene. The observed effect size of T2D on fracture risk decreased if the T2D-related risk factors could be eliminated. Therefore, it is important to manage the complications of T2D to prevent the risk of fracture.