Cell freezing protocol suitable for ATAC-Seq on motor neurons derived from human induced pluripotent stem cells

Genetic Associations With an Amyotrophic Lateral Sclerosis Reversal Phenotype

BACKGROUND AND OBJECTIVES

The term "ALS Reversal" describes patients who initially meet diagnostic criteria for amyotrophic lateral sclerosis (ALS) or had clinical features most consistent with progressive muscular atrophy (PMA) but subsequently demonstrated substantial and sustained clinical improvement. The objective of this genome-wide association study (GWAS) was to identify correlates of this unusual clinical phenotype.

METHODS

Participants were recruited from a previously created database of individuals with the ALS Reversal phenotype. Whole-genome sequencing (WGS) data were compared with ethnicity-matched patients with typically progressive ALS enrolled through the CReATe Consortium's Phenotype-Genotype-Biomarker (PGB) study. These results were replicated using an independent ethnically matched WGS data set from Target ALS. Significant results were further explored with available databases of genetic regulatory markers and expression quantitative trait loci (eQTL) analysis.

RESULTS

WGS from 22 participants with documented ALS Reversals was compared with the PGB primary cohort (n = 103) and the Target ALS validation cohort (n = 140). Two genetic loci met predefined criteria for statistical significance (two-sided permutation p ≤ 0.01) and remained plausible after fine-mapping. The lead single nucleotide variant (SNV) from the first locus was rs4242007 (primary cohort GWAS OR = 12.0, 95% CI 4.1 to 34.6), which is in an IGFBP7 intron and is in near-perfect linkage disequilibrium with a SNV in the IGFBP7 promoter region. Both SNVs are associated with decreased frontal cortex IGFBP7 expression in eQTL data sets. Notably, 3 Reversals, but none of the typically progressive individuals (n = 243), were homozygous for rs4242007. The importance of the second locus, located near GRIP1, is uncertain given the absence of an associated effect on nearby gene transcription.

DISCUSSION

We found a significant association between the Reversal phenotype and an IGFBP7 noncoding SNV that is associated with IGFBP7 expression. This is biologically relevant as IGFBP7 is a reported inhibitor of the insulin growth factor-1 (IGF-1) receptor that activates the possibly neuroprotective IGF-1 signaling pathway. This finding is limited by small sample size but suggests that there may be merit in further exploration of IGF-1 pathway signaling as a therapeutic mechanism for ALS.

TRIAL REGISTRATION INFORMATION

This study was registered with ClinicalTrials.gov (NCT03464903) on March 14, 2018. The first participant was enrolled on June 22, 2018.

Tissue characterization at an enhanced resolution across spatial omics platforms with deep generative model

Recent advances in spatial omics have expanded the spectrum of profiled molecular categories beyond transcriptomics. However, many of these technologies are constrained by limited spatial resolution, hindering our ability to deeply characterize intricate tissue architectures. Existing computational methods primarily focus on the resolution enhancement of transcriptomics data, lacking the adaptability to address the emerging spatial omics technologies that profile various omics types. Here, we introduce soScope, a unified generative framework designed to enhance data quality and spatial resolution for molecular profiles obtained from diverse spatial technologies. soScope aggregates multimodal tissue information from omics, spatial relations and images, and jointly infers omics profiles at enhanced resolutions with omics-specific modeling through distribution priors. With comprehensive evaluations on diverse spatial omics platforms, including Visium, Xenium, spatial-CUT&Tag, and slide-DNA/RNA-seq, soScope improves performances in identifying biologically meaningful intestine and kidney architectures, revealing embryonic heart structure that cannot be resolved at the original resolution and correcting sample and technical biases arising from sequencing and sample processing. Furthermore, soScope extends to spatial multiomics technology spatial-CITE-seq and spatial ATAC-RNA-seq, leveraging cross-omics reference for simultaneous multiomics enhancement. soScope provides a versatile tool to improve the utilization of continually expanding spatial omics technologies and resources.

Disease related changes in ATAC-seq of iPSC-derived motor neuron lines from ALS patients and controls

Amyotrophic Lateral Sclerosis (ALS), like many other neurodegenerative diseases, is highly heritable, but with only a small fraction of cases explained by monogenic disease alleles. To better understand sporadic ALS, we report epigenomic profiles, as measured by ATAC-seq, of motor neuron cultures derived from a diverse group of 380 ALS patients and 80 healthy controls. We find that chromatin accessibility is heavily influenced by sex, the iPSC cell type of origin, ancestry, and the inherent variance arising from sequencing. Once these covariates are corrected for, we are able to identify ALS-specific signals in the data. Additionally, we find that the ATAC-seq data is able to predict ALS disease progression rates with similar accuracy to methods based on biomarkers and clinical status. These results suggest that iPSC-derived motor neurons recapitulate important disease-relevant epigenomic changes.

Vocal learning-associated convergent evolution in mammalian proteins and regulatory elements

Vocal production learning ("vocal learning") is a convergently evolved trait in vertebrates. To identify brain genomic elements associated with mammalian vocal learning, we integrated genomic, anatomical, and neurophysiological data from the Egyptian fruit bat (Rousettus aegyptiacus) with analyses of the genomes of 215 placental mammals. First, we identified a set of proteins evolving more slowly in vocal learners. Then, we discovered a vocal motor cortical region in the Egyptian fruit bat, an emergent vocal learner, and leveraged that knowledge to identify active cis-regulatory elements in the motor cortex of vocal learners. Machine learning methods applied to motor cortex open chromatin revealed 50 enhancers robustly associated with vocal learning whose activity tended to be lower in vocal learners. Our research implicates convergent losses of motor cortex regulatory elements in mammalian vocal learning evolution.

Cell-type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell-type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell-type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell-type-specific cis-regulatory changes. We find that cell-type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell-type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Improved quality metrics for association and reproducibility in chromatin accessibility data using mutual information

BACKGROUND

Correlation metrics are widely utilized in genomics analysis and often implemented with little regard to assumptions of normality, homoscedasticity, and independence of values. This is especially true when comparing values between replicated sequencing experiments that probe chromatin accessibility, such as assays for transposase-accessible chromatin via sequencing (ATAC-seq). Such data can possess several regions across the human genome with little to no sequencing depth and are thus non-normal with a large portion of zero values. Despite distributed use in the epigenomics field, few studies have evaluated and benchmarked how correlation and association statistics behave across ATAC-seq experiments with known differences or the effects of removing specific outliers from the data. Here, we developed a computational simulation of ATAC-seq data to elucidate the behavior of correlation statistics and to compare their accuracy under set conditions of reproducibility.

RESULTS

Using these simulations, we monitored the behavior of several correlation statistics, including the Pearson's R and Spearman's [Formula: see text] coefficients as well as Kendall's [Formula: see text] and Top-Down correlation. We also test the behavior of association measures, including the coefficient of determination R[Formula: see text], Kendall's W, and normalized mutual information. Our experiments reveal an insensitivity of most statistics, including Spearman's [Formula: see text], Kendall's [Formula: see text], and Kendall's W, to increasing differences between simulated ATAC-seq replicates. The removal of co-zeros (regions lacking mapped sequenced reads) between simulated experiments greatly improves the estimates of correlation and association. After removing co-zeros, the R[Formula: see text] coefficient and normalized mutual information display the best performance, having a closer one-to-one relationship with the known portion of shared, enhanced loci between simulated replicates. When comparing values between experimental ATAC-seq data using a random forest model, mutual information best predicts ATAC-seq replicate relationships.

CONCLUSIONS

Collectively, this study demonstrates how measures of correlation and association can behave in epigenomics experiments. We provide improved strategies for quantifying relationships in these increasingly prevalent and important chromatin accessibility assays.

Cell type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell type-specific cis-regulatory changes. We find that cell type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Integrative profiling of gene expression and chromatin accessibility elucidates specific transcriptional networks in porcine neutrophils

Neutrophils are vital components of the immune system for limiting the invasion and proliferation of pathogens in the body. Surprisingly, the functional annotation of porcine neutrophils is still limited. The transcriptomic and epigenetic assessment of porcine neutrophils from healthy pigs was performed by bulk RNA sequencing and transposase accessible chromatin sequencing (ATAC-seq). First, we sequenced and compared the transcriptome of porcine neutrophils with eight other immune cell transcriptomes to identify a neutrophil-enriched gene list within a detected neutrophil co-expression module. Second, we used ATAC-seq analysis to report for the first time the genome-wide chromatin accessible regions of porcine neutrophils. A combined analysis using both transcriptomic and chromatin accessibility data further defined the neutrophil co-expression network controlled by transcription factors likely important for neutrophil lineage commitment and function. We identified chromatin accessible regions around promoters of neutrophil-specific genes that were predicted to be bound by neutrophil-specific transcription factors. Additionally, published DNA methylation data from porcine immune cells including neutrophils were used to link low DNA methylation patterns to accessible chromatin regions and genes with highly enriched expression in porcine neutrophils. In summary, our data provides the first integrative analysis of the accessible chromatin regions and transcriptional status of porcine neutrophils, contributing to the Functional Annotation of Animal Genomes (FAANG) project, and demonstrates the utility of chromatin accessible regions to identify and enrich our understanding of transcriptional networks in a cell type such as neutrophils.

Parallel recovery of chromatin accessibility and gene expression dynamics from frozen human regulatory T cells

Epigenetic features such as DNA accessibility dictate transcriptional regulation in a cell type- and cell state- specific manner, and mapping this in health vs. disease in clinically relevant material is opening the door to new mechanistic insights and new targets for therapy. Assay for Transposase Accessible Chromatin Sequencing (ATAC-seq) allows chromatin accessibility profiling from low cell input, making it tractable on rare cell populations, such as regulatory T (Treg) cells. However, little is known about the compatibility of the assay with cryopreserved rare cell populations. Here we demonstrate the robustness of an ATAC-seq protocol comparing primary Treg cells recovered from fresh or cryopreserved PBMC samples, in the steady state and in response to stimulation. We extend this method to explore the feasibility of conducting simultaneous quantitation of chromatin accessibility and transcriptome from a single aliquot of 50,000 cryopreserved Treg cells. Profiling of chromatin accessibility and gene expression in parallel within the same pool of cells controls for cellular heterogeneity and is particularly beneficial when constrained by limited input material. Overall, we observed a high correlation of accessibility patterns and transcription factor dynamics between fresh and cryopreserved samples. Furthermore, highly similar transcriptomic profiles were obtained from whole cells and from the supernatants recovered from ATAC-seq reactions. We highlight the feasibility of applying these techniques to profile the epigenomic landscape of cells recovered from cryopreservation biorepositories.

Robust Generation of Ready-to-Use Cryopreserved Motor Neurons from Human Pluripotent Stem Cells for Disease Modeling

Human pluripotent stem cell (hPSC)-derived motor neurons (MNs) act as models for motor neuron diseases (MNDs), such as amyotrophic lateral sclerosis (ALS) or spinal muscular atrophy. However, the MN differentiation efficiency and viability following cryopreservation require further development for application in large-scale studies and drug screening. Here, we developed a robust protocol to convert hPSCs into MN cryopreservation stocks (hPSCs were converted into >92% motor neural progenitors and >91% MNs). Near-mature MNs were cryopreserved at a high thawing survival rate and 89% MN marker expression on day 32. Moreover, these MNs exhibited classical electrophysiological properties and neuromuscular junction (NMJ) formation ability within only 4−6 days after thawing. To apply this platform as an MND model, MN stocks were generated from SOD1G85R, SOD1G85G isogenic control, and sporadic ALS hPSC lines. The thawed ALS MNs expressed ALS-specific cytopathies, including SOD1 protein aggregation and TDP-43 redistribution. Thus, a stable and robust protocol was developed to generate ready-to-use cryopreserved MNs without further neuronal maturation processes for application in MND mechanistic studies, NMJ model establishment, and large-scale drug screening.

Spatial-CUT&Tag: Spatially resolved chromatin modification profiling at the cellular level

Spatial omics emerged as a new frontier of biological and biomedical research. Here, we present spatial-CUT&Tag for spatially resolved genome-wide profiling of histone modifications by combining in situ CUT&Tag chemistry, microfluidic deterministic barcoding, and next-generation sequencing. Spatially resolved chromatin states in mouse embryos revealed tissue-type-specific epigenetic regulations in concordance with ENCODE references and provide spatial information at tissue scale. Spatial-CUT&Tag revealed epigenetic control of the cortical layer development and spatial patterning of cell types determined by histone modification in mouse brain. Single-cell epigenomes can be derived in situ by identifying 20-micrometer pixels containing only one nucleus using immunofluorescence imaging. Spatial chromatin modification profiling in tissue may offer new opportunities to study epigenetic regulation, cell function, and fate decision in normal physiology and pathogenesis.

Answer ALS, a large-scale resource for sporadic and familial ALS combining clinical and multi-omics data from induced pluripotent cell lines

Answer ALS is a biological and clinical resource of patient-derived, induced pluripotent stem (iPS) cell lines, multi-omic data derived from iPS neurons and longitudinal clinical and smartphone data from over 1,000 patients with ALS. This resource provides population-level biological and clinical data that may be employed to identify clinical-molecular-biochemical subtypes of amyotrophic lateral sclerosis (ALS). A unique smartphone-based system was employed to collect deep clinical data, including fine motor activity, speech, breathing and linguistics/cognition. The iPS spinal neurons were blood derived from each patient and these cells underwent multi-omic analytics including whole-genome sequencing, RNA transcriptomics, ATAC-sequencing and proteomics. The intent of these data is for the generation of integrated clinical and biological signatures using bioinformatics, statistics and computational biology to establish patterns that may lead to a better understanding of the underlying mechanisms of disease, including subgroup identification. A web portal for open-source sharing of all data was developed for widespread community-based data analytics.

Cell type-specific chromatin accessibility analysis in the mouse and human brain

The Assay for Transposase Accessible Chromatin by sequencing (ATAC-seq) is becoming popular in the neuroscience field where chromatin regulation is thought to be involved in neurodevelopment, activity-dependent gene regulation, hormonal and environmental responses, and pathophysiology of neuropsychiatric disorders. The advantages of using ATAC-seq include a small amount of material needed, fast protocol, and the ability to capture a range of gene regulatory elements with a single assay. With increasing interest in chromatin research, it is an imperative to have feasible, reliable assays that are compatible with a range of neuroscience study designs. Here we tested three protocols for neuronal chromatin accessibility analysis, including a varying brain tissue freezing method followed by fluorescence-activated nuclei sorting (FANS) and ATAC-seq. Our study shows that the cryopreservation method impacts the number of open chromatin regions identified from frozen brain tissue using ATAC-seq. However, we show that all protocols generate consistent and robust data and enable the identification of functional regulatory elements in neuronal cells. Our study implies that the broad biological interpretation of chromatin accessibility data is not significantly affected by the freezing condition. We also reveal additional challenges of doing chromatin analysis on post-mortem human brain tissue. Overall, ATAC-seq coupled with FANS is a powerful method to capture cell-type-specific chromatin accessibility information in mouse and human brain. Our study provides alternative brain preservation methods that generate high-quality ATAC-seq data while fitting in different study designs, and further encourages the use of this method to uncover the role of epigenetic (dys)regulation in the brain.

An integrated multi-omic analysis of iPSC-derived motor neurons from C9ORF72 ALS patients

Neurodegenerative diseases are challenging for systems biology because of the lack of reliable animal models or patient samples at early disease stages. Induced pluripotent stem cells (iPSCs) could address these challenges. We investigated DNA, RNA, epigenetics, and proteins in iPSC-derived motor neurons from patients with ALS carrying hexanucleotide expansions in C9ORF72. Using integrative computational methods combining all omics datasets, we identified novel and known dysregulated pathways. We used a C9ORF72 Drosophila model to distinguish pathways contributing to disease phenotypes from compensatory ones and confirmed alterations in some pathways in postmortem spinal cord tissue of patients with ALS. A different differentiation protocol was used to derive a separate set of C9ORF72 and control motor neurons. Many individual -omics differed by protocol, but some core dysregulated pathways were consistent. This strategy of analyzing patient-specific neurons provides disease-related outcomes with small numbers of heterogeneous lines and reduces variation from single-omics to elucidate network-based signatures.

Successful ATAC-Seq From Snap-Frozen Equine Tissues

An assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) has become an increasingly popular method to assess genome-wide chromatin accessibility in isolated nuclei from fresh tissues. However, many biobanks contain only snap-frozen tissue samples. While ATAC-seq has been applied to frozen brain tissues in human, its applicability in a wide variety of tissues in horse remains unclear. The Functional Annotation of Animal Genome (FAANG) project is an international collaboration aimed to provide high quality functional annotation of animal genomes. The equine FAANG initiative has generated a biobank of over 80 tissues from two reference female animals and experiments to begin to characterize tissue specificity of genome function for prioritized tissues have been performed. Due to the logistics of tissue collection and storage, extracting nuclei from a large number of tissues for ATAC-seq at the time of collection is not always practical. To assess the feasibility of using stored frozen tissues for ATAC-seq and to provide a guideline for the equine FAANG project, we compared ATAC-seq results from nuclei isolated from frozen tissue to cryopreserved nuclei (CN) isolated at the time of tissue harvest in liver, a highly cellular homogenous tissue, and lamina, a relatively acellular tissue unique to the horse. We identified 20,000-33,000 accessible chromatin regions in lamina and 22-61,000 in liver, with consistently more peaks identified using CN isolated at time of tissue collection. Our results suggest that frozen tissues are an acceptable substitute when CN are not available. For more challenging tissues such as lamina, nuclei extraction at the time of tissue collection is still preferred for optimal results. Therefore, tissue type and accessibility to intact nuclei should be considered when designing ATAC-seq experiments.

Downregulation by CNNM2 of ATP5MD expression in the 10q24.32 schizophrenia-associated locus involved in impaired ATP production and neurodevelopment

Genome-wide association studies (GWAS) have accelerated the discovery of numerous genetic variants associated with schizophrenia. However, most risk variants show a small effect size (odds ratio (OR) <1.2), suggesting that more functional risk variants remain to be identified. Here, we employed region-based multi-marker analysis of genomic annotation (MAGMA) to identify additional risk loci containing variants with large OR value from Psychiatry Genomics Consortium (PGC2) schizophrenia GWAS data and then employed summary-data-based mendelian randomization (SMR) to prioritize schizophrenia susceptibility genes. The top-ranked susceptibility gene ATP5MD, encoding an ATP synthase membrane subunit, is observed to be downregulated in schizophrenia by the risk allele of CNNM2-rs1926032 in the schizophrenia-associated 10q24.32 locus. The Atp5md knockout (KO) in mice was associated with abnormal startle reflex and gait, and ATP5MD knockdown (KD) in human induced pluripotent stem cell-derived neurons disrupted the neural development and mitochondrial respiration and ATP production. Moreover, CNNM2-rs1926032 KO could induce downregulation of ATP5MD expression and disruptions of mitochondrial respiration and ATP production. This study constitutes an important mechanistic component that links schizophrenia-associated CNNM2 regions to disruption in energy adenosine system modulation and neuronal function by long-distance chromatin domain downregulation of ATP5MD. This pathogenic mechanism provides therapeutic implications for schizophrenia.

Isolation of genetically manipulated neural progenitors and immature neurons from embryonic mouse neocortex by FACS

The embryonic mammalian neocortex includes neural progenitors and neurons at various stages of differentiation. The regulatory mechanisms underlying multiple aspects of neocortical development—including cell division, neuronal fate commitment, neuronal migration, and neuronal differentiation—have been explored using in utero electroporation and virus infection. Here, we describe a protocol for investigation of the effects of genetic manipulation on neural development through direct isolation of neural progenitors and neurons from the mouse embryonic neocortex by fluorescence-activated cell sorting.

For complete details on the use and execution of this protocol, please refer to Tsuboi et al. (2018) and Sakai et al. (2019).

Direct isolation of neural progenitors and neurons from the mouse embryonic neocortex

Purification of neural cells at various stages of differentiation using wild-type mice

Protocol enables biochemical analyses of neural cells genetically manipulated in vivo

Applicable for transcriptomic and epigenomic analyses

Chromatin accessibility in neuropsychiatric disorders

Epigenetic mechanisms have recently emerged as critical regulators of brain function in health and disease. By controlling the accessibility and the expression of specific genes, these pathways can mediate transient and long-lasting changes in neuronal function in both physiological and pathological contexts. Due to the extreme complexity of the epigenetic regulatory landscape, emerging methods that directly assay chromatin accessibility are of particular interest. Here, I review recent insights gained on open and closed chromatin states in the brain, with emphasis on neuropsychiatric disorders. These advances generated an invaluable wealth of information that can help us better understand gene regulation in the brain and its impairments that contribute to the development of disease.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Profiling chromatin regulatory landscape: insights into the development of ChIP-seq and ATAC-seq

Chromatin regulatory landscape plays a critical role in many disease processes and embryo development. Epigenome sequencing technologies such as chromatin immunoprecipitation sequencing (ChIP-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) have enabled us to dissect the pan-genomic regulatory landscape of cells and tissues in both time and space dimensions by detecting specific chromatin state and its corresponding transcription factors. Pioneered by the advancement of chromatin immunoprecipitation-chip (ChIP-chip) technology, abundant epigenome profiling technologies have become available such as ChIP-seq, DNase I hypersensitive site sequencing (DNase-seq), ATAC-seq and so on. The advent of single-cell sequencing has revolutionized the next-generation sequencing, applications in single-cell epigenetics are enriched rapidly. Epigenome sequencing technologies have evolved from low-throughput to high-throughput and from bulk sample to the single-cell scope, which unprecedentedly benefits scientists to interpret life from different angles. In this review, after briefly introducing the background knowledge of epigenome biology, we discuss the development of epigenome sequencing technologies, especially ChIP-seq & ATAC-seq and their current applications in scientific research. Finally, we provide insights into future applications and challenges.

Systematic alteration of ATAC-seq for profiling open chromatin in cryopreserved nuclei preparations from livestock tissues

The use of Assay for Transposase-Accessible Chromatin (ATAC-seq) to profile chromatin accessibility has surged over the past years, but its applicability to tissues has been very limited. With the intent of preserving nuclear architecture during long-term storage, cryopreserved nuclei preparations from chicken lung were used to optimize ATAC-seq. Sequencing data were compared with existing DNase-seq, ChIP-seq, and RNA-seq data to evaluate library quality, ultimately resulting in a modified ATAC-seq method capable of generating high quality chromatin accessibility data from cryopreserved nuclei preparations. Using this method, nucleosome-free regions (NFR) identified in chicken lung overlapped half of DNase-I hypersensitive sites, coincided with active histone modifications, and specifically marked actively expressed genes. Notably, sequencing only the subnucleosomal fraction dramatically improved signal, while separation of subnucleosomal reads post-sequencing did not improve signal or peak calling. The broader applicability of this modified ATAC-seq technique was tested using cryopreserved nuclei preparations from pig tissues, resulting in NFR that were highly consistent among biological replicates. Furthermore, tissue-specific NFR were enriched for binding motifs of transcription factors related to tissue-specific functions, and marked genes functionally enriched for tissue-specific processes. Overall, these results provide insights into the optimization of ATAC-seq and a platform for profiling open chromatin in animal tissues.

Aberrant Development Corrected in Adult-Onset Huntington's Disease iPSC-Derived Neuronal Cultures via WNT Signaling Modulation

Aberrant neuronal development and the persistence of mitotic cellular populations have been implicated in a multitude of neurological disorders, including Huntington's disease (HD). However, the mechanism underlying this potential pathology remains unclear. We used a modified protocol to differentiate induced pluripotent stem cells (iPSCs) from HD patients and unaffected controls into neuronal cultures enriched for medium spiny neurons, the cell type most affected in HD. We performed single-cell and bulk transcriptomic and epigenomic analyses and demonstrated that a persistent cyclin D1+ neural stem cell (NSC) population is observed selectively in adult-onset HD iPSCs during differentiation. Treatment with a WNT inhibitor abrogates this NSC population while preserving neurons. Taken together, our findings identify a mechanism that may promote aberrant neurodevelopment and adult neurogenesis in adult-onset HD striatal neurons with the potential for therapeutic compensation.

Detect accessible chromatin using ATAC-sequencing, from principle to applications

Background

Chromatin accessibility is crucial for gene expression regulation in specific cells and in multiple biological processes. Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) is an effective way to reveal chromatin accessibility at a genome-wide level. Through ATAC-seq, produced reads from a small number of cells reflect accessible regions that correspond to nucleosome positioning and transcription factor binding sites, due to probing hyperactive Tn5 transposase to DNA sequence.

Conclusion

In this review, we summarize both principle and features of ATAC-seq, highlight its applications in basic and clinical research. ATAC-seq has generated comprehensive chromatin accessible maps, and is becoming a powerful tool to understand dynamic gene expression regulation in stem cells, early embryos and tumors.

Hair-follicle-associated pluripotent stem cells derived from cryopreserved intact human hair follicles sustain multilineage differentiation potential

The bulge area of the hair follicle contains hair-follicle-associated pluripotent (HAP) stem cells. Here, we present effective cryopreservation procedures of the human hair follicle that preserve the differentiation potential of HAP stem cells. Whole hair follicles isolated from human scalp were cryopreserved by a slow-rate cooling medium and stored in liquid nitrogen. A careful thawing method was used to collect the upper parts of the human hair follicles which were cultured for four weeks in a Dulbecco’s Modified Eagle’s Medium with fetal bovine serum (FBS). Proliferating hair follicle cells were then shifted to DMEM/Ham’s Nutrient Mixture F-12 medium without FBS and allowed to grow for one week. These proliferating cells were able to produce HAP stem cell colonies with multilineage differentiation capacity. They produced keratinocytes, smooth muscle cells, cardiac muscle cells, neurons and glial cells. Interestingly, these cryopreserved hair follicles produced pluripotent HAP stem cell colonies similar to fresh follicles. These findings suggest that the cryopreserved whole human hair follicle preserves the ability to produce HAP stem cells, which will enable any individual to preserve a bank of these stem cells for personalized regenerative medicine.

An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

High Quality ATAC-Seq Data Recovered from Cryopreserved Breast Cell Lines and Tissue

DNA accessibility to transcription regulators varies between cells and modulates gene expression patterns. Several “open” chromatin profiling methods that provide valuable insight into the activity of these regulatory regions have been developed. However, their application to clinical samples has been limited despite the discovery that the Analysis of Transposase-Accessible Chromatin followed by sequencing (ATAC-seq) method can be performed using fewer cells than other techniques. Obtaining fresh rather than stored samples and a lack of adequate optimization and quality controls are major barriers to ATAC’s clinical implementation. Here, we describe an optimized ATAC protocol in which we varied nuclear preparation conditions and transposase concentrations and applied rigorous quality control measures before testing fresh, flash frozen, and cryopreserved breast cells and tissue. We obtained high quality data from small cell number. Furthermore, the genomic distribution of sequencing reads, their enrichment at transcription start sites, and transcription factor footprint analyses were similar between cryopreserved and fresh samples. This updated method is applicable to clinical samples, including cells from fine needle aspiration and tissues obtained via core needle biopsy or surgery. Chromatin accessibility analysis using patient samples will greatly expand the range of translational research and personalized medicine by identification of clinically-relevant epigenetic features.

Prenatal Alcohol Exposure: Profiling Developmental DNA Methylation Patterns in Central and Peripheral Tissues

Background: Prenatal alcohol exposure (PAE) can alter the development of neurobiological systems, leading to lasting neuroendocrine, neuroimmune, and neurobehavioral deficits. Although the etiology of this reprogramming remains unknown, emerging evidence suggests DNA methylation as a potential mediator and biomarker for the effects of PAE due to its responsiveness to environmental cues and relative stability over time. Here, we utilized a rat model of PAE to examine the DNA methylation profiles of rat hypothalami and leukocytes at four time points during early development to assess the genome-wide impact of PAE on the epigenome and identify potential biomarkers of PAE. Our model of PAE resulted in blood alcohol levels of ~80-150 mg/dl throughout the equivalent of the first two trimesters of human pregnancy. Hypothalami were analyzed on postnatal days (P) 1, 8, 15, 22 and leukocytes at P22 to compare central and peripheral markers. Genome-wide DNA methylation analysis was performed by methylated DNA immunoprecipitation followed by next-generation sequencing. Results: PAE resulted in lasting changes to DNA methylation profiles across all four ages, with 118 differentially methylated regions (DMRs) displaying persistent alterations across the developmental period at a false-discovery rate (FDR) < 0.05. In addition, 299 DMRs showed the same direction of change in the hypothalamus and leukocytes of P22 pups at an FDR < 0.05, with some genes overlapping with the developmental profile findings. The majority of these DMRs were located in intergenic regions, which contained several computationally-predicted transcription factor binding sites. Differentially methylated genes were generally involved in immune function, epigenetic remodeling, metabolism, and hormonal signaling, as determined by gene ontology analyses. Conclusions: Persistent DNA methylation changes in the hypothalamus may be associated with the long-term physiological and neurobehavioral alterations in observed in PAE. Furthermore, correlations between epigenetic alterations in peripheral tissues and those in the brain will provide a foundation for the development of biomarkers of fetal alcohol spectrum disorder (FASD). Finally, findings from studies of PAE provide important insight into the etiology of neurodevelopmental and mental health disorders, as they share numerous phenotypes and comorbidities.

Chromatin accessibility landscape of articular knee cartilage reveals aberrant enhancer regulation in osteoarthritis

Osteoarthritis (OA) is a common joint disorder with increasing impact in an aging society. While genetic and transcriptomic analyses have revealed some genes and non-coding loci associated to OA, the pathogenesis remains incompletely understood. Chromatin profiling, which provides insight into gene regulation, has not been reported in OA mainly due to technical difficulties. Here, we employed Assay for Transposase-Accessible Chromatin with high throughput sequencing (ATAC-seq) to map the accessible chromatin landscape in articular knee cartilage of OA patients. We identified 109,215 accessible chromatin regions for cartilages, of which 71% were annotated as enhancers. By overlaying them with genetic and DNA methylation data, we have determined potential OA-relevant enhancers and their putative target genes. Furthermore, through integration with RNA-seq data, we characterized genes that are altered both at epigenomic and transcriptomic levels in OA. These genes are enriched in pathways regulating ossification and mesenchymal stem cell (MSC) differentiation. Consistently, the differentially accessible regions in OA are enriched for MSC-specific enhancers and motifs of transcription factor families involved in osteoblast differentiation. In conclusion, we demonstrate how direct chromatin profiling of clinical tissues can provide comprehensive epigenetic information for a disease and suggest candidate genes and enhancers of translational potential.

ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data

BACKGROUND

ATAC-seq (Assays for Transposase-Accessible Chromatin using sequencing) is a recently developed technique for genome-wide analysis of chromatin accessibility. Compared to earlier methods for assaying chromatin accessibility, ATAC-seq is faster and easier to perform, does not require cross-linking, has higher signal to noise ratio, and can be performed on small cell numbers. However, to ensure a successful ATAC-seq experiment, step-by-step quality assurance processes, including both wet lab quality control and in silico quality assessment, are essential. While several tools have been developed or adopted for assessing read quality, identifying nucleosome occupancy and accessible regions from ATAC-seq data, none of the tools provide a comprehensive set of functionalities for preprocessing and quality assessment of aligned ATAC-seq datasets.

RESULTS

We have developed a Bioconductor package, ATACseqQC, for easily generating various diagnostic plots to help researchers quickly assess the quality of their ATAC-seq data. In addition, this package contains functions to preprocess aligned ATAC-seq data for subsequent peak calling. Here we demonstrate the utilities of our package using 25 publicly available ATAC-seq datasets from four studies. We also provide guidelines on what the diagnostic plots should look like for an ideal ATAC-seq dataset.

CONCLUSIONS

This software package has been used successfully for preprocessing and assessing several in-house and public ATAC-seq datasets. Diagnostic plots generated by this package will facilitate the quality assessment of ATAC-seq data, and help researchers to evaluate their own ATAC-seq experiments as well as select high-quality ATAC-seq datasets from public repositories such as GEO to avoid generating hypotheses or drawing conclusions from low-quality ATAC-seq experiments. The software, source code, and documentation are freely available as a Bioconductor package at https://bioconductor.org/packages/release/bioc/html/ATACseqQC.html .

Assay for Transposase Accessible Chromatin (ATAC-Seq) to Chart the Open Chromatin Landscape of Human Pancreatic Islets

The regulatory mechanisms that ensure an accurate control of gene transcription are central to cellular function, development and disease. Such mechanisms rely largely on noncoding regulatory sequences that allow the establishment and maintenance of cell identity and tissue-specific cellular functions.The study of chromatin structure and nucleosome positioning allowed revealing transcription factor accessible genomic sites with regulatory potential, facilitating the comprehension of tissue-specific cis-regulatory networks. Recently a new technique coupled with high-throughput sequencing named Assay for Transposase Accessible Chromatin (ATAC-seq) emerged as an efficient method to chart open chromatin genome wide. The application of such technique to different cell types allowed unmasking tissue-specific regulatory elements and characterizing cis-regulatory networks. Herein we describe the implementation of the ATAC-seq method to human pancreatic islets, a tissue playing a central role in the control of glucose metabolism.

Integrated Analysis of Whole-Genome ChIP-Seq and RNA-Seq Data of Primary Head and Neck Tumor Samples Associates HPV Integration Sites with Open Chromatin Marks

Chromatin alterations mediate mutations and gene expression changes in cancer. Chromatin immunoprecipitation followed by sequencing (ChIP-Seq) has been utilized to study genome-wide chromatin structure in human cancer cell lines, yet numerous technical challenges limit comparable analyses in primary tumors. Here we have developed a new whole-genome analytic pipeline to optimize ChIP-Seq protocols on patient-derived xenografts from human papillomavirus-related (HPV⁺) head and neck squamous cell carcinoma (HNSCC) samples. We further associated chromatin aberrations with gene expression changes from a larger cohort of the tumor and normal samples with RNA-Seq data. We detect differential histone enrichment associated with tumor-specific gene expression variation, sites of HPV integration in the human genome, and HPV-associated histone enrichment sites upstream of cancer driver genes, which play central roles in cancer-associated pathways. These comprehensive analyses enable unprecedented characterization of the complex network of molecular changes resulting from chromatin alterations that drive HPV-related tumorigenesis. Cancer Res; 77(23); 6538-50. ©2017 AACR.

Epigenetics of Atherosclerosis: Emerging Mechanisms and Methods

Atherosclerosis is a vascular pathology characterized by inflammation and plaque build-up within arterial vessel walls. Vessel occlusion, often occurring after plaque rupture, can result in myocardial and cerebral infarction. Epigenetic changes are increasingly being associated with atherosclerosis and are of interest from both therapeutic and biomarker perspectives. Emerging genomic approaches that profile DNA methylation, chromatin accessibility, post-translational histone modifications, transcription factor binding, and RNA expression in low or single cell populations are poised to enhance our spatiotemporal understanding of atherogenesis. Here, we review recent therapeutically relevant epigenetic discoveries and emerging technologies that may generate new opportunities for atherosclerosis research.

Single-cell transcriptome conservation in cryopreserved cells and tissues

A variety of single-cell RNA preparation procedures have been described. So far, protocols require fresh material, which hinders complex study designs. We describe a sample preservation method that maintains transcripts in viable single cells, allowing one to disconnect time and place of sampling from subsequent processing steps. We sequence single-cell transcriptomes from >1000 fresh and cryopreserved cells using 3'-end and full-length RNA preparation methods. Our results confirm that the conservation process did not alter transcriptional profiles. This substantially broadens the scope of applications in single-cell transcriptomics and could lead to a paradigm shift in future study designs.

Genome-wide epigenomic profiling for biomarker discovery

A myriad of diseases is caused or characterized by alteration of epigenetic patterns, including changes in DNA methylation, post-translational histone modifications, or chromatin structure. These changes of the epigenome represent a highly interesting layer of information for disease stratification and for personalized medicine. Traditionally, epigenomic profiling required large amounts of cells, which are rarely available with clinical samples. Also, the cellular heterogeneity complicates analysis when profiling clinical samples for unbiased genome-wide biomarker discovery. Recent years saw great progress in miniaturization of genome-wide epigenomic profiling, enabling large-scale epigenetic biomarker screens for disease diagnosis, prognosis, and stratification on patient-derived samples. All main genome-wide profiling technologies have now been scaled down and/or are compatible with single-cell readout, including: (i) Bisulfite sequencing to determine DNA methylation at base-pair resolution, (ii) ChIP-Seq to identify protein binding sites on the genome, (iii) DNaseI-Seq/ATAC-Seq to profile open chromatin, and (iv) 4C-Seq and HiC-Seq to determine the spatial organization of chromosomes. In this review we provide an overview of current genome-wide epigenomic profiling technologies and main technological advances that allowed miniaturization of these assays down to single-cell level. For each of these technologies we evaluate their application for future biomarker discovery. We will focus on (i) compatibility of these technologies with methods used for clinical sample preservation, including methods used by biobanks that store large numbers of patient samples, and (ii) automation of these technologies for robust sample preparation and increased throughput.