DRMref Logo

Home

Download

Statistics

Landscape

Help

Contact

Center for Computational Systems Medicine
leaf

Dataset summary

leaf

Drug summary

leaf

Umap of single cell types and conditions (resistant and sensitive)

leaf

Comparison of cell composition between the resistant and sensitive groups (This dataset does not contain this module)

leaf

Comparison of intra-tumor heterogeneity (ITH) and epithelial-mesenchymal transition (EMT) scores of malignant cells between the resistant and sensitive groups

leaf

Difference of cell-cell interactions between the resistant and sensitive groups

Soverview

Differentially expressed genes (DEGs) between the resistant and sensitive groups for each cell type

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in the pre-treatment samples (This dataset does not contain this module)

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in the post-treatment samples (This dataset does not contain this module)

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in one patient or cell lines

Soverview

MicroRNAs (miRNAs) regulating drug resistance-related DEGs for Malignant cells

Soverview

Motifs and transcription factors (TFs) regulating drug resistance-related DEGs for each cell type

     •

Differentially expressed genes act as transcription factors

Dataset: GSE164551

Dataset summary for GSE164551

Datast informationDatasetGSE164551
PMID33558511
Raw data ID/linkhttps://doi.org/10.7910/DVN/1RKYQ8
OrganismHomo sapiens
Sourcepatients
TissueBone marrow aspirate
Cancer type level1Multiple myeloma
Cancer type level2Refractory multiple myeloma (MM)
RegimenCAR-T
Drug typeImmunotherapy
Sample sizeresistant 1, sensitive 1
Cell number7736
Extract protocol10x genomics
Data processingCellRanger 3.1.0
Public datePublic on Jan 11, 2021
DescriptionThis dataset has 1 patient with a sensitive pre-treatment sample and a resistant post-treatment sample.

Top

Drug summary for GSE164551

Drug nameDrugBank IDDrug typeTargets nameTargets uniprot ID
"CAR-T" is not included in the drug list.

Top

Umap of single cell types and conditions (resistant and sensitive)

check buttonUmap of cell types and conditions (resistant and sensitive). If this dataset has both pre-treatment and post-treatment samples, the cell types were annotated together. If this dataset is a cell line dataset, the UMAP of cell clusters and the cell ratio of drug-resistant and sensitive groups within each cluster will also be shown. (If the image exists, the user can click on it to enlarge it in a new window.)
Umap of cell typesboxplotUmap of conditionsboxplot

Top

Comparison of cell composition between the resistant and sensitive groups

check buttonIn the pre-treatment samples, the cell composition of the resistant and sensitive groups were compared.
* represents significant difference.
(If the image exists, the user can click on it to enlarge it in a new window.)

check buttonIn the post-treatment samples, the cell composition of the resistant and sensitive groups were compared.
* represents significant difference.
(If the image exists, the user can click on it to enlarge it in a new window.)

Top

Comparison of ITH and EMT scores of malignant cells between the resistant and sensitive groups

check buttonComparison of ITH (intra-tumoral heterogeneity) and EMT (epithelial-mesenchymal transition) scores of malignant cells between the resistant and sensitive groups. ITH represents intratumoral heterogeneity, EMT represents epithelial mesenchymal transition, pre represents pre-treatment, post represents post-treatment, post2 represents longer post-treatment. (If the image exists, the user can click on it to enlarge it in a new window.)
Comparison of ITH scoresboxplotComparison of EMT scoresboxplot

Top

Difference of cell-cell interactions between the resistant and sensitive groups

check buttonIn the pre-treatment samples, the difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonIn the post-treatment samples, the difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The difference in cell-cell interactions between the resistant and sensitive groups was compared. The numbers indicate the count of cell-cell interactions that were either higher or lower in the resistant group compared to the sensitive group. The dot plots of significant ligand-receptor pairs in the resistant group and sensitive group are subsequently displayed. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

resistant vs. sensitive

boxplot
Dot plot of significant ligand-receptor pairs in the resistant groupDot plot of significant ligand-receptor pairs in the sensitive group
boxplotboxplot

Top

Differentially expressed genes (DEGs) between the resistant and sensitive groups for each cell type

check buttonIn the pre-treatment samples, the differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.

check buttonIn the post-treatment samples, the differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The following are differentially expressed genes (|avg_log2FC|>0.25 & p_val_adj<0.05) between the resistant and sensitive groups for each cell type.
Gene symbolGene idavg_log2FCp_valp_val_adjpct.1pct.2Cell typeTimepoint

PSMB6

ENSG00000142507-0.2967193.92e-216.17e-170.0750.262ErythrocytesNA

RBM8A

ENSG00000265241-0.2973671.25e-161.97e-120.1080.277ErythrocytesNA

E2F4

ENSG00000205250-0.300857.98e-171.26e-120.2140.416ErythrocytesNA

ELOB

ENSG00000103363-0.3023381.47e-102.31e-060.4480.549ErythrocytesNA

NDUFA3

ENSG00000276061-0.3028056.40e-221.01e-170.0850.282ErythrocytesNA

NAP1L4

ENSG00000273562-0.3030488.90e-161.40e-110.2120.398ErythrocytesNA

TFR2

ENSG00000106327-0.3031882.06e-153.24e-110.1540.327ErythrocytesNA

BSDC1

ENSG00000160058-0.3058052.40e-123.79e-080.2430.408ErythrocytesNA

ST13

ENSG00000100380-0.3059277.39e-121.17e-070.2650.411ErythrocytesNA

ZNF451

ENSG00000112200-0.3064151.21e-191.90e-150.4280.683ErythrocytesNA

CLK1

ENSG00000013441-0.3068243.24e-125.11e-080.1880.343ErythrocytesNA

ARPC3

ENSG00000111229-0.3068366.18e-179.74e-130.1320.31ErythrocytesNA

TOMM20

ENSG00000173726-0.3070641.10e-181.74e-140.090.267ErythrocytesNA

USP14

ENSG00000101557-0.3070672.21e-183.48e-140.1010.283ErythrocytesNA

RHCE

ENSG00000188672-0.3083661.18e-131.86e-090.2170.389ErythrocytesNA

GPCPD1

ENSG00000125772-0.3096373.81e-096.01e-050.2990.441ErythrocytesNA

SLC38A5

ENSG00000017483-0.3118092.80e-184.41e-140.1570.35ErythrocytesNA

HIPK1

ENSG00000163349-0.3118585.48e-148.64e-100.170.334ErythrocytesNA

AP2M1

ENSG00000161203-0.3120981.45e-162.29e-120.2010.386ErythrocytesNA

NDUFA1

ENSG00000125356-0.3129791.15e-181.82e-140.1130.301ErythrocytesNA

YWHAH

ENSG00000128245-0.3134335.49e-258.66e-210.0490.25ErythrocytesNA

MIS18BP1

ENSG00000129534-0.3152174.18e-206.58e-160.0750.261ErythrocytesNA

HLA-C

ENSG00000206435-0.316033.43e-225.41e-180.1830.424ErythrocytesNA

COX8A

ENSG00000176340-0.318124.67e-207.36e-160.0820.265ErythrocytesNA

SUZ12

ENSG00000178691-0.3186958.29e-151.31e-100.1080.267ErythrocytesNA

SLC2A1

ENSG00000117394-0.3193653.96e-176.24e-130.5150.763ErythrocytesNA

RBM3

ENSG00000102317-0.320053.66e-185.78e-140.1030.285ErythrocytesNA

CD47

ENSG00000196776-0.3216819.96e-161.57e-110.1370.312ErythrocytesNA

SOX6

ENSG00000110693-0.3217133.12e-144.92e-100.2340.411ErythrocytesNA

WSB1

ENSG00000109046-0.3218922.47e-193.90e-150.0820.262ErythrocytesNA

SERF2

ENSG00000140264-0.3220122.79e-084.40e-040.7960.787ErythrocytesNA

PSMD6

ENSG00000163636-0.3220512.11e-173.33e-130.1370.317ErythrocytesNA

EIF2S2

ENSG00000125977-0.3225393.11e-184.90e-140.1030.281ErythrocytesNA

GCLM

ENSG00000023909-0.3243472.98e-204.69e-160.1080.304ErythrocytesNA

ANP32E

ENSG00000143401-0.326221.40e-162.20e-120.1130.288ErythrocytesNA

TCP1

ENSG00000120438-0.328177.40e-211.17e-160.070.253ErythrocytesNA

ARGLU1

ENSG00000134884-0.3283094.91e-197.74e-150.0850.263ErythrocytesNA

ENY2

ENSG00000120533-0.3290455.00e-227.87e-180.0780.273ErythrocytesNA

ATP6V0A1

ENSG00000033627-0.3292041.55e-162.44e-120.1240.299ErythrocytesNA

SRSF5

ENSG00000100650-0.3300951.97e-173.11e-130.1650.358ErythrocytesNA

TFDP2

ENSG00000114126-0.3303553.65e-165.76e-120.2450.433ErythrocytesNA

PABPN1

ENSG00000100836-0.3311152.77e-244.36e-200.0570.257ErythrocytesNA

EDF1

ENSG00000107223-0.3335411.82e-242.86e-200.0540.251ErythrocytesNA

HBD

ENSG00000223609-0.3338646.18e-079.73e-030.9360.909ErythrocytesNA

VAPA

ENSG00000101558-0.3341652.31e-173.64e-130.1810.377ErythrocytesNA

FRYL

ENSG00000075539-0.3342543.65e-145.75e-100.250.43ErythrocytesNA

NDUFC1

ENSG00000109390-0.3352953.95e-246.23e-200.070.275ErythrocytesNA

HLA-B

ENSG00000206450-0.3367181.26e-201.99e-160.1750.412ErythrocytesNA

UGCG

ENSG00000148154-0.3367944.50e-177.10e-130.1230.296ErythrocytesNA

PNISR

ENSG00000132424-0.3390271.57e-172.48e-130.1470.335ErythrocytesNA
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in the pre-treatment samples

check buttonEnrichment results for 5 known drug-resistant mechanisms in malignant cells and 1 known mechanism in TME cells.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group of pre-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group of pre-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window.Complete files can be downloaded from the download section.)

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in the post-treatment samples

check buttonEnrichment results for 5 known drug resistance mechanisms in malignant cells and 1 known mechanism in TME cells.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group of post-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group of post-treatment samples. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)

Top

Mechanism analysis of drug resistance-related DEGs for each cell type in one patient or cell lines

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. Enrichment results for 5 known drug resistance mechanisms in malignant cells and 1 known mechanism in TME cells are shown here.The "Signature" column represents the differentially expressed genes (DEGs) used for enrichment. The "Geneset" column represents the genes included in this mechanism. The "Overlap" column represents the number of genes in the "Geneset" that are enriched in the "Signature". The "Enriched Genes" column represents the genes in the "Geneset" that are enriched in the "Signature". (If the image exists, the user can click on it to enlarge it in a new window.)
Enrichment results for 5 known drug resistance mechanisms in malignant cells
boxplot
MechnismPvalFDRSignatureGenesetOverlapBackgroundCell typeEnriched Genes
Unusual Activation of Pro-survival Pathway5.90e-022.90e-0112781223467Malignant cellsEGR1,CXCL8
Epigenetic Alteration of DNA, RNA or Protein2.10e-015.20e-0112455123467Malignant cellsCXCL8
Aberration of the Drug's Therapeutic Target1.00e+001.00e+001290023467Malignant cells
Drug Inactivation by Structure Modification1.00e+001.00e+001227023467Malignant cells
Irregularity in Drug Uptake and Drug Efflux1.00e+001.00e+001237023467Malignant cells

Enrichment results for 1 known drug resistance mechanisms in TME cells

boxplot
MechnismPvalFDRSignatureGenesetOverlapBackgroundCell typeEnriched Genes
Regulation by the Disease Microenvironment1.00e+001.00e+003186023467B cellsNA
Regulation by the Disease Microenvironment6.60e-036.60e-03148186523467CD4+ T cellsBRAF,TGFB1,SOCS3,VIM,ZEB1
Regulation by the Disease Microenvironment1.80e-011.80e-01276186423467CD8+ T cellsNFKB1,TGFB1,CXCR4,VIM
Regulation by the Disease Microenvironment5.10e-015.10e-01467186423467ErythrocytesBCL2L1,MTDH,HNRNPA2B1,MALAT1
Regulation by the Disease Microenvironment9.20e-039.20e-035241861023467Mono_MacroBRAF,CTNNB1,TGFB1,VEGFA,ZEB2,IL1B,HNRNPA2B1,SOCS3,VIM,MALAT1
Regulation by the Disease Microenvironment4.80e-014.80e-0181186123467NK cellsVIM

check buttonHallmark, KEGG and GOBP enrichment results for up-regulated genes in resistant group. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)
Up-regulated HallmarkUp-regulated KEGGUp-regulated GO BP

boxplot

boxplot

boxplot

check buttonHallmark, KEGG and GOBP enrichment results for down-regulated genes in resistant group. Only show the top 50 pathways with fdr<0.05. If certain cell types are not displayed, it means that there are no significant pathways based on DEGs in those cell types. (If the image exists, the user can click on it to enlarge it in a new window. Complete files can be downloaded from the download section.)
Down-regulated HallmarkDown-regulated KEGGDown-regulated GO BP

boxplot

boxplot

boxplot

Top

MicroRNAs (miRNAs) regulating drug resistance-related DEGs for malignant cells

check buttonIn the pre-treatment samples, the miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for malignant cells are shown here. If the gene does not have any miRNAs with a score greater than 80, it will not be shown. Here only the names of miRNAs with a score greater than 90 are shown, but some names will not be shown because of too many overlaps. (Complete files containing all datasets and cell types can be downloaded from the download section.)

check buttonIn the post-treatment samples, the miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for Malignant cells are shown here. If the gene does not have any miRNAs with a score greater than 80, it will not be shown. Here only the names of miRNAs with a score greater than 90 are shown, but some names will not be shown because of too many overlaps. (Complete files containing all datasets and cell types can be downloaded from the download section.)

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. The miRNAs with a score greater than 80 that regulate the top10 abs(avg_log2FC) drug resistance-related DEGs for Malignant cells are shown here. Here only the names of miRNAs with a score greater than 80 are shown, but some names will not be shown because of too many overlaps. If the gene does not have any miRNAs with a score greater than 90, it will not be shown. (Complete files containing all datasets and cell types can be downloaded from the download section.)
Regulating up-regulated DEGs in resistant groupRegulating down-regulated DEGs in resistant group
This dataset does not contain malignant cells, nor does it predict any miRNA that corresponds to the top10 DEGs.boxplot

Top

Motifs and transcription factors (TFs) regulating drug resistance-related DEGs for each cell type

check buttonIn the pre-treatment samples, the motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)

check buttonIn the post-treatment samples, the motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)

check buttonThis dataset only has one patient or is a cell line dataset, the resistant sample is the post-treatment sample and the sensitive sample is the pre-treatment sample. Motifs and TFs that regulate drug resistance-related DEGs for each cell type are shown here. (Complete files can be downloaded from the download section.)
Motifs and TFs regulating up-regulated DEGs in resistant group
motifmotifNESAUCTF_highConfCell_type
motiftransfac_pro__M028484.050.193EGR1 (inferredBy_Orthology). CD4+ T cells
motiftaipale_tf_pairs__ETV5_HES7_NNCACGTGNNNNCCGGAANN_CAP4.020.192ETV5; HES7 (directAnnotation). CD4+ T cells
motiftransfac_pro__M074753.810.185TBX5 (inferredBy_Orthology). CD4+ T cells
motifcisbp__M007753.790.184ZBTB1 (inferredBy_Orthology). CD4+ T cells
motifmetacluster_84.23.760.183RARA (directAnnotation). CD4+ T cells
motiftaipale_tf_pairs__E2F1_ELK1_SGCGCNNNNNNNNNNCGGAAGN_CAP_repr3.740.182E2F1; ELK1 (directAnnotation). CD4+ T cells
motifmetacluster_141.53.710.181NR3C1; THAP1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY1; YY2; YY2; ZNF597 (directAnnotation). YY1; YY1; YY1; YY2; ZFP42 (inferredBy_Orthology). CD4+ T cells
motifjaspar__MA1650.13.680.18ZBTB14 (directAnnotation). CD4+ T cells
motiftransfac_pro__M020903.620.178E2F4 (directAnnotation). CD4+ T cells
motiftransfac_pro__M080053.560.176TCFL5 (inferredBy_Orthology). CD4+ T cells
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28

Motifs and TFs regulating down-regulated DEGs in resistant group
motifmotifNESAUCTF_highConfCell_type
motiftransfac_pro__M012527.710.61E2F6 (directAnnotation). B cells
motiftfdimers__MD003107.020.558IKZF1; OTX1; OTX2 (directAnnotation). B cells
motifmetacluster_76.16.650.53ZNF345; ZNF345 (directAnnotation). MYB; MYBL1 (inferredBy_Orthology). B cells
motifmetacluster_156.26.560.523ATF3; ATF4; ATF4; ATF4; ATF4; ATF4; ATF4; CEBPG; CEBPG; DDIT3; MYC (directAnnotation). ATF4; ATF4; ATF4; CEBPG; DDIT3 (inferredBy_Orthology). B cells
motiftransfac_pro__M055316.50.518IKZF5 (inferredBy_Orthology). B cells
motiftfdimers__MD001036.440.514ZNF333 (directAnnotation). B cells
motifcisbp__M008386.020.482FOXP1 (directAnnotation). B cells
motifmetacluster_130.15.990.48ZNF713; ZNF713 (directAnnotation). HSF1; HSF2; HSF4; HSF5; HSFX1; HSFX2; HSFX3; HSFX4; HSFY1; HSFY2 (inferredBy_Orthology). B cells
motifmetacluster_139.95.890.473TFE3 (inferredBy_Orthology). B cells
motiftfdimers__MD002985.880.471CEBPA; CEBPB; CEBPD; CEBPE; CEBPG; FOXL1 (directAnnotation). B cells
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57

Top

Differentially expressed genes act as transcription factors

check buttonThis table shows the differentially expressed genes that act as transcription factors. (Complete files can be downloaded from the download section.)
TFmotifExpressionCell typeTimepoint

ZEB1

tfdimers__MD00527upCD4+ T cellsNA

MECP2

metacluster_135.10upCD8+ T cellsNA

E2F4

metacluster_45.6downErythrocytesNA

KLF13

hocomoco__KLF13_HUMAN.H11MO.0.DdownErythrocytesNA

YBX1

metacluster_79.33downErythrocytesNA

BACH1

metacluster_22.37upMono/MacroNA

BACH1

metacluster_157.2upMono/MacroNA

BACH1

metacluster_30.4upMono/MacroNA

EGR1

cisbp__M01864upMono/MacroNA

EGR1

metacluster_131.7upMono/MacroNA
Page: 1




1"Sun X, Zhang Y, Li H, Zhou Y, Shi S, Chen Z, He X, Zhang H, Li F, Yin J, Mou M, Wang Y, Qiu Y, Zhu F. DRESIS: the first comprehensive landscape of drug resistance information. Nucleic Acids Res. 2023 Jan 6;51(D1):D1263-D1275. doi: 10.1093/nar/gkac812. PMID: 36243960; PMCID: PMC9825618."
2"Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, Sajed T, Johnson D, Li C, Sayeeda Z, Assempour N, Iynkkaran I, Liu Y, Maciejewski A, Gale N, Wilson A, Chin L, Cummings R, Le D, Pon A, Knox C, Wilson M. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018 Jan 4;46(D1):D1074-D1082. doi: 10.1093/nar/gkx1037. PMID: 29126136; PMCID: PMC5753335."