|
|
GEO help: Mouse over screen elements for information. |
|
Status |
Public on May 23, 2021 |
Title |
O2_PreTx |
Sample type |
SRA |
|
|
Source name |
Hematopoietic stem cells
|
Organism |
Mus musculus |
Characteristics |
age at procedure: 24 months old transposition day: 1 Sex: Male
|
Extracted molecule |
genomic DNA |
Extraction protocol |
Hematopoietic stem cells were isolated from the bone marrow of mice through a standardized process of bone marrow extraction, creation of a single cell suspension, and flow cytometry based cell sorting. Cells were lysed and transposed following the Corces et al. Nature Methods 2017 Omni-ATAC seq protocol. Transposed DNA was purifed using a Zymo 5 Clean and Concentrator kit and stored at -20 C until library preparation. Libraries were quality controlled using Qubit and Agilent Bioanalyzer. ATAC-seq libraries were made with Nextera sequencing adapters and NEBnext Master Mix as in Corces et al. Nature Methods 2017. Libraries were purified with double-sided AMPure XP beads.
|
|
|
Library strategy |
ATAC-seq |
Library source |
genomic |
Library selection |
other |
Instrument model |
Illumina HiSeq 4000 |
|
|
Data processing |
Sequencing quality control was performed using FastQC and MultiQC, which revealed Nextera sequencing adapter contamination that was removed using the Trim Galore Cutadapt wrapper in paired end mode, followed by confirmation of adapter removal with FastQC. For pathway enrichment analysis, we annotated promoter proximal peaks within ± 1000bp of the transcription start site (TSS) using TxDb.Mmusculus.UCSC.mm10.knownGene. Rsubread was used to align samples to the mm10 genome. Rsamtools was used to sort and index output BAM files. BAM files were quality controlled by looking at the distribution of reads aligning to chromosomes and the distribution of insert sizes. We next utilized Genrich (Harvard FAS Informatics) to call peaks in ATAC seq mode using a maximum q value threshold of 0.05 for peak calling, removing PCR duplicates and excluding reads that mapped to the mitochondrial genome and the Y chromosome. NarrowPeak files generated by Genrich were quality controlled for distribution of reads mapping to mm10 genomic regions using the R CHIPseeker package and Tx.Db.mmusculus.UCSC.mm10.knowngene annotation. For differential ATAC-seq, a GRanges consensus peak object was generated using ChIPQC48 containing non-redundant open regions present in any sample, excluding known mm10 ATAC-seq blacklisted regions obtained from Encode (ENCFF547MET). Nucleosome free regions (<100bp fragment length) from sample BAM files were counted using FeatureCounts against the consensus peak annotation. Sample counts were imported into a dds object using DESeq2, and exploratory data analysis was performed using rlog normalized count data and PCA. For differential peak accessibility calling between groups, we used DESeq2 using the results function and the DESeq2 estimator. Genome_build: We used ReactomePA to test for Reactome pathway enrichment in peaks with increased relative peak accessibility in a given comparison using a log2 fold change threshold of >= 0.05 and a pvalue cutoff of 0.05.
|
|
|
Submission date |
Jun 02, 2020 |
Last update date |
May 23, 2021 |
Contact name |
Paul Vincent Dellorusso |
E-mail(s) |
[email protected]
|
Phone |
5166552491
|
Organization name |
Columbia University Irving Medical Center
|
Department |
Genetics and Development
|
Lab |
Passegue Lab
|
Street address |
701 West 168th Street, 724 Hammer Health Sciences Center
|
City |
New York |
State/province |
NY |
ZIP/Postal code |
10032 |
Country |
USA |
|
|
Platform ID |
GPL21103 |
Series (1) |
GSE151665 |
Young and Old Pre- and Post-Transplantation ATAC Seq |
|
Relations |
BioSample |
SAMN15084600 |
SRA |
SRX8455210 |
Supplementary file |
Size |
Download |
File type/resource |
GSM4587843_O2-PreTx.narrowPeak.gz |
780.9 Kb |
(ftp)(http) |
NARROWPEAK |
SRA Run Selector |
Raw data are available in SRA |
Processed data provided as supplementary file |
|
|
|
|
|