As a complementary approach to identify variants located in genomic regions with regulatory activity in plasma cells, we analyzed previously published31 (link) chromatin immunoprecipitation sequencing (ChIP-seq) data for the H3K4me3 histone modification38 (link). Briefly, L363 cells (DSMZ) were cross-linked with 1% paraformaldehyde (ThermoFisher, #28908). DNA was sonicated into 200–400 bp fragments (Bioruptor Pico Sonication System, Diagenode, Belgium). For pull-down, we used 1–10 μg of H3K4me3 antibody (Millipore, #04-745). Fragments were de-cross-linked and purified (Zymogen, #D5205). ChIP-seq libraries were prepared using the ThruPLEX DNA-seq Kit (Rubicon Genomics, #R400406) and sequenced on Illumina HiSeq 2500 sequencer (paired-end; 2 × 125 cycles). De-multiplexing and generation of FASTQ files was performed using bcl2fastq v.1.8. FastQC (v0.11.5)39 (link) was used to assess read quality low-quality bases were removed using Trimmomatic (v.0.36)40 (link),41 (link) prior to alignment. using Bowtie2 (v.2.3.0)41 (link). Coverage in 50 bp over the SOHL2 region was calculated with the GenomicAlignments and GenomicRanges R-packages42 (coverage and binnedAverage functions) and scaled to Counts-per-million (CPM) relative to the total number of reads per library.
Free full text: Click here