Samples
In order to determine the potential utility of targeting BRD4 in the treatment of AML, we analyzed the expression of BRD4 based on public RNA-seq data in AML samples. The standardized gene expression data was used to assess the prognostic significance of BRD4 and the correlation between the two, and the overall survival rate of BRD4 and AML patients was also analyzed.
Cell culture
The human leukemia cell lines (NB4, Kasumi-1, HL-60, MV4-11 and K562) and mouse leukemia cell line (P388-D1) were all from the Chinese Academy of Sciences Cell Bank, they were all verified by short tandem repeat analysis in 2019 and 2020. Cells were cultured in RPMI1640 medium containing 10% fetal bovine serum (Termo Fisher Scientifc), and penicillin and streptomycin (Millipore, Billerica, MA, USA), 37℃, 5% CO2, cultured in a humidified incubator, and routinely tested for mycoplasma.
Cell viability determination
GNE-987 is dissolved in 100% DMSO, the stock solution has a concentration of 10 mmol, and it is placed in a refrigerator at -80℃. AML cells were planted in a 96-well cell culture plate with a cell density of 2×104 in each well, and were treated with GNE-987 with different concentration gradients. The primary leukemia cells were separated from the bone marrow of children by Ficoll-Hypake centrifugation, and then planted in a 96-well plate with a density of 1×105 cells in the culture medium. Cells treated with 0.05% dimethyl sulfoxide (DMSO) in complete medium without GNE-987 were used as controls. After 24 hours of drug treatment, according to the manufacturer's instructions, the cell viability was determined by the Cell Counting Kit-8 (CCK-8) assay (Dojindo Molecular Technologies, Tokyo, Japan). Each concentration is in triplicate, repeated in at least three independent experiments. Graph Prism software 8.3.0(GraphPad Software Inc., San Diego, CA, USA) was used to calculate the half maximum inhibitory concentration (IC50) of GNE-987.
Soft agar clone formation analysis
Prepare 1.2% and 0.7% agarose, autoclave and place in a 55℃ water bath, prepare 2×RPMI1640 medium containing 20% FBS, 2× penicillin and streptomycin, and filter out with a 0.2 micron filter bacteria. Lower layer gel: 1.2% agarose gel was mixed with 2× medium 1:1, added to a 6-well plate, 1.5ml per well, and solidified at room temperature. Cell count: AML cells treated with GNE-987 at different concentration gradients were washed with PBS, mixed and diluted with new medium and adjusted to 5×103/ ml with 100ul of cells suspension in each well. Upper glue: Mix 0.7% agarose gel with 2× medium 1:1, add 100ul cell suspension, after mixing wells, add 1.5ml to each well. Put it in a 37°C, CO2 incubator, add culture medium every three days, and harvest the cells after about 3 weeks. Count and compare the number of cells treated with different concentrations of GNE-987, and finally analyze to calculate the rate of monoclonal formation.
Preparation and infection of lentivirus
Short hairpin RNA (shRNA) targeting VHL (GGAGCCTAGTCAAGCCTGAGA CATCCGTTGATGTGCAATGCG) constructed in the pLKO.1 lentiviral vector. The CDS region of the VHL gene was searched in pubmed, synthesized and constructed into the PLVX-EF1a-puro vector. Short hairpin RNA (shRNA) targeting LYL1 (Table 1) were constructed in the pLKO.1-puro lentiviral vector (IGE BIOTECHNOLOGY LTD, Guangzhou, China). When preparing lentivirus, envelope plasmid and packaging plasmid were purchased from Addgene (pMD2.G: #12,259; psPAX2: #12,260; Cambridge, MA, USA). Co-transfect pMD2.G, psPAX2 and the transfer plasmid with polyethyleneimine into 293FT cells (linear MW 25,000 Da, 5mg/ml, pH 7.0) (cat. No. 23966-1; Polysciences, Warrington, PA, USA) according to the manufacturer’s instructions. After 6h, completely replace the medium with fresh medium. Collect the virus supernatant 48h after transfection and filter it with a 0.22μm filter. Then, AML cells were infected with lentivirus for 24 hours in the presence of 10 μg/mL polyene (Sigma-Aldrich). Screen stable cell lines with puromycin (Sigma-Aldrich).
RNA preparation and real-time PCR expression analysis
Total RNA was extracted from cell pellets using TRIzol®reagent (Invitrogen, CA, USA), according to the manufacturer’s protocol. For cDNA synthesis, 1µg of total RNA was converted to cDNA using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, CA, USA). Quantitative real-time PCR analysis was carried out using LightCycler® 480 SYBR Green I Master mix (cat. No. 04707516001; Roche, Penzberg, Germany) with a LightCycler 480 Real Time System (Roche), according to the manufacturer’s protocol. mRNA expression levels were calculated using the Ct method with glyceraldehyde 3-phosphate dehydrogenase (GAPDH) expression as an internal reference. Primer sequences are listed in Table 2.
Cellcycleanalysis
24 hours after adding different concentrations of GNE-987 to the AML cell line, the cell line was trypsinized, washed and fixed in 70% ethanol at 4℃ overnight. Then wash the cells with cold phosphate buffered saline (PBS) and resuspend them in 0.5 ml of PI/RNase staining fermentation broth (cat. No. 550825; BD Pharmingen™, San Diego, CA, USA), then incubate at room temperature for 15 min. Flow cytometry was performed using Beckman Gallios™ Flow Cytometer (Beckman, Krefeld, Germany), and the cell cycle was analyzed using MultiCycle AV DNA analysis software (Verity Software House, Topsham, ME, USA).
Cell apoptosis analysis
Add different concentrations of GNE-987 to the cell line, collect it after 24h, centrifuge at 2000rpm, 3min, discard the supernatant, wash once with cold PBS, centrifuge at 4000rpm, 3min, take the supernatant, suspend it in 1× binding buffer, use fluorescein isothiocyanate (FITC)- AnnexinV apoptosis kit and PI solution staining (cat. No. 556420; BD Biosciences, Franklin Lakes, NJ, USA), follow the manufacturer's instructions. Cell counting method was used to analyze cell apoptosis (Beckman Gallios™ Flow Cytometer; Beckman).
Western blotting analysis
Use the following antibodies for Western blotting analysis, BRD2 (cat.No.5848 s; 1:1000; Cell Signaling Technology, Boston, MA, USA), BRD3 (cat. No. 11859-1-AP; 1:1000; Proteintech, Chicago, IL, USA), BRD4 (cat. No. 13440 s; 1:1000; Cell Signaling Technology), VHL (cat. No. 68547 s; 1:1000; Cell Signaling Technology), LYL1(cat. No. sc-374164; 1:1000; Santa Cruz Biotechnology), and PARP (cat. No. 9542; 1:1000; Cell Signaling Technology) with glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (cat. No. MA3374; 1:1000; Millipore) as a reference protein. Peroxidase-conjugated Afniure goat anti-rabbit IgG(H+L) IgG(H+L) (cat. 111-035-003; 1:5000) and goat anti-mouse IgG(H+L) (cat. No. 115-035-003; 1:5000) purchased from Jackson ImmunoResearch Laboratories, Inc. (West Grove, PA, USA). Define the role of the proteasome, MG132 (cat. No. 474787, Sigma-Aldrich, St. Louis, MO, USA) inhibits the proteasome activity. After 24h treatment with different concentrations of GNE-987, the cells were collected, and the BRD2, BRD3, BRD4, PARP, VHL, LYL1 and GAPDH proteins were determined by Western blotting analysis.
Study the anti-tumor effect of GNE-987 in Vivo
All experimental animal procedures in this study were approved and licensed by the Animal Care and Use Committee of Children's Hospital of Soochow University (CAMSU-AP#: JP-2018-1). SPF grade BALB/c mice were from Linghang Biotechnology Co., Ltd. (Shanghai, China). Five-week-old female mice (n=5 in each group) were injected with 3×105 P388-D1 cells via the tail vein. 2 days after injection of cells, each mouse was injected with luciferase into the abdominal cavity and immediately anesthetized with isoflurane gas. Then, each group of mice was imaged using the NightOWL in-vivo Imaging System (BERTHOLD, Germany). After the tumor fluorescence signal appeared (Day 2), mice in the experimental group were injected intraperitoneally with 0.5mg/KgGNE-987; in the control group, mice were injected intraperitoneally with GNE-987 (5%®HS15), once a day, 9 times in total (Day2~Day10). Continue to use the NightOWL in-vivo Imaging System to image each group of mice on the 4th, 7th, and 10th days. The mice were weighed daily to observe the color and mobility of the fur. The liver, spleen, kidneys and intestines of the experimental group and control group mice were collected, and the organ size was observed and weighed. Each organ specimen was subjected to immunohistochemistry and HE staining. Primary antibody BRD4 (cat: No. 13440 s; 1:1000; Cell Signaling Technology), cleaved-caspase 3 (cat: No. GB11009-1, 1:300, Servicebio, Boston, MA, USA) and Ki67 (cat: No. ab15580, 1:300, Abcam, Cambridge, UK), has been used according to the manufacturer’s recommendations.
RNA sequencing and data processing
RNA-seq was carried out according to the protocols suggested by Novogene Bioinformatics Technology Co., Ltd. (Beijing, China). First, the total RNA was reverse transcribed into cDNA to construct a library, and then the cDNA library was sequenced. Filter the original reads and map the clean reads according to HISAT. Then calculate the gene expression level (as the fragment mapped per million reads per kilobase exon model). Using DESeq2 analysis, differentially expressed genes were identified (P<0.05 and fold change>2 or fold change<0.5). For enrichment analysis, differentially expressed genes were analyzed using the GSEA software (UC San Diego and Broad Institute,).
Chromatin immunoprecipitation sequencing (ChIP-seq)
3-5×107 cells were cross-linked with 1% formaldehyde for 10 minutes, and neutralized with 1.25M glycine at room temperature for 5 minutes. Use Bioruptor (Diagenode, Liège, Belgium) to collect, lyse and sonicate fixed cells. Sonicated chromatin was incubated with anti-histone H3 (acetyl K27) antibody (cat. No. ab4729; Abcam, Cambridge, UK) overnight at 4°C.DNA was eluted and purified using a QIAquick PCR purification kit (cat. No. 208106; Qiagen, Hilden, Germany). The samples are sequenced on the novaseq 6000 platform (Novogene Bioinformatics Technology Co., Ltd. Beijing, China) and a BGISEQ 2000 platform (Beijing Genomics Institute (Shenzhen, China)). Raw data of ChIP-Seq H3K27ac analysis was aligned to the reference genome (UCSC hg38) using Bowtie2 (v 2.3.5) [17], with alignment parameters -p 4 -q -x. Peaks were identified using MACS2 (v2.0.9) [18], with parameters -g hs -n test -B -q 0.01. The bedgraph files generated by MACS2 were converted to bigwig files using the UCSC bedGraphToBigWig tool, and then bigwig files were visualized by Integrative Genomics Viewer (IGV) [19]. Super-enhancers were then identified using the ROSE (Rank Order of Super Enhancers) method [20, 21], with parameters -s 12500 -t 2000.
Statistical analysis
All experiments were carried out independently at least three times. Statistical analysis was carried out using SPSS software version 21.0 (IBM Corporation, NY, USA). Survival analysis was performed by Kaplan-Meier estimates with log-rank tests. Student's t test was used to compare the percentage of apoptosis, BRD4 mRNA level and cell viability. Normally distributed measurement data are expressed as mean± standard deviation. Use the t test to compare the differences between the two groups. Non-normally distributed data are expressed in quartiles (usually the median plus the range) and are compared using the Mann-Whitney U test. All experiments are two-tailed, P<0.05 is statistically significant.