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 (This dataset does not contain this module)

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

Soverview

Mechanism analysis of drug resistance-related DEGs for each cell type in one patient or cell lines (This dataset does not contain this module)

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: PMID34715028

Dataset summary for PMID34715028

Datast informationDatasetPMID34715028
PMID34715028
Raw data ID/linkEGAD00001008166
OrganismHomo sapiens
Sourcepatients
TissueTumor tissue
Cancer type level1Kidney cancer
Cancer type level2Clear cell renal cell carcinoma (ccRCC)
Regimennivolumab
Drug typeImmunotherapy
Sample sizepost (resistant 1, sensitive 1)
Cell number27332
Extract protocol10x genomics
Data processingCellRanger 5.0.0
Public datePublic on Oct 28, 2021
DescriptionThis dataset has 2 patients with post-treatment samples.

Top

Drug summary for PMID34715028

Drug nameDrugBank IDDrug typeTargets nameTargets uniprot ID
"Nivolumab"

DB09035

biotechPDCD1; CD274Q15116; Q9NZQ7

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.)

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.)

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

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.)

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.
Gene symbolGene idavg_log2FCp_valp_val_adjpct.1pct.2Cell typeTimepoint

BTG1

ENSG00000133639-0.3276990.00e+000.00e+000.9780.993CD8+ T cellspost

RAB11FIP1

ENSG00000156675-0.330410.00e+000.00e+000.1580.329CD8+ T cellspost

ELL2

ENSG00000118985-0.3318550.00e+000.00e+000.1290.317CD8+ T cellspost

TNFRSF9

ENSG00000049249-0.3340470.00e+000.00e+000.2860.445CD8+ T cellspost

JMY

ENSG00000152409-0.3354850.00e+000.00e+000.1120.297CD8+ T cellspost

NR3C1

ENSG00000113580-0.3367320.00e+000.00e+000.4080.548CD8+ T cellspost

RALGAPA1

ENSG00000174373-0.339440.00e+000.00e+000.2150.396CD8+ T cellspost

AZIN1

ENSG00000155096-0.341330.00e+000.00e+000.1570.343CD8+ T cellspost

NAA50

ENSG00000121579-0.3418120.00e+000.00e+000.170.356CD8+ T cellspost

ELF1

ENSG00000120690-0.343690.00e+000.00e+000.5120.654CD8+ T cellspost

CALM2

ENSG00000143933-0.343820.00e+000.00e+000.6650.779CD8+ T cellspost

FYN

ENSG00000010810-0.3444870.00e+000.00e+000.5310.703CD8+ T cellspost

GNG2

ENSG00000186469-0.3473440.00e+000.00e+000.4520.588CD8+ T cellspost

NFE2L2

ENSG00000116044-0.3503840.00e+000.00e+000.190.381CD8+ T cellspost

EIF4A2

ENSG00000156976-0.3514940.00e+000.00e+000.5490.704CD8+ T cellspost

OAZ1

ENSG00000104904-0.3523540.00e+000.00e+000.8220.907CD8+ T cellspost

DNAJB9

ENSG00000128590-0.3524570.00e+000.00e+000.0750.288CD8+ T cellspost

DTHD1

ENSG00000197057-0.3528540.00e+000.00e+000.2780.464CD8+ T cellspost

ISG20

ENSG00000172183-0.3556110.00e+000.00e+000.3310.476CD8+ T cellspost

PLK3

ENSG00000173846-0.3558250.00e+000.00e+000.0660.254CD8+ T cellspost

RELB

ENSG00000104856-0.3561380.00e+000.00e+000.0810.296CD8+ T cellspost

RSRC2

ENSG00000111011-0.3598780.00e+000.00e+000.3180.509CD8+ T cellspost

EMD

ENSG00000102119-0.3609170.00e+000.00e+000.2940.497CD8+ T cellspost

ARIH1

ENSG00000166233-0.3624070.00e+000.00e+000.2380.429CD8+ T cellspost

PTMA

ENSG00000187514-0.3627550.00e+000.00e+000.9890.999CD8+ T cellspost

MX1

ENSG00000157601-0.3628450.00e+000.00e+000.1320.284CD8+ T cellspost

SRSF5

ENSG00000100650-0.3632340.00e+000.00e+000.6340.771CD8+ T cellspost

PNPLA8

ENSG00000135241-0.3637760.00e+000.00e+000.1020.305CD8+ T cellspost

ARID4B

ENSG00000054267-0.3658280.00e+000.00e+000.4740.639CD8+ T cellspost

GTF2B

ENSG00000137947-0.3673920.00e+000.00e+000.1410.359CD8+ T cellspost

IL21R

ENSG00000103522-0.3686660.00e+000.00e+000.1620.362CD8+ T cellspost

HLA-C

ENSG00000206435-0.3687410.00e+000.00e+000.9991CD8+ T cellspost

SFPQ

ENSG00000116560-0.3717920.00e+000.00e+000.5170.679CD8+ T cellspost

PECAM1

ENSG00000261371-0.374720.00e+000.00e+000.1970.428CD8+ T cellspost

STIP1

ENSG00000168439-0.3824340.00e+000.00e+000.1880.385CD8+ T cellspost

ARF4

ENSG00000168374-0.3832610.00e+000.00e+000.2050.425CD8+ T cellspost

FAM118A

ENSG00000100376-0.3838230.00e+000.00e+000.0940.314CD8+ T cellspost

HERPUD1

ENSG00000051108-0.3875430.00e+000.00e+000.3150.504CD8+ T cellspost

MCOLN2

ENSG00000153898-0.3896510.00e+000.00e+000.1390.319CD8+ T cellspost

XBP1

ENSG00000100219-0.3944930.00e+000.00e+000.2010.401CD8+ T cellspost

PLIN2

ENSG00000147872-0.3958940.00e+000.00e+000.1050.276CD8+ T cellspost

LITAF

ENSG00000189067-0.396250.00e+000.00e+000.3520.528CD8+ T cellspost

ZFP36L1

ENSG00000185650-0.3963530.00e+000.00e+000.6390.753CD8+ T cellspost

NSMCE3

ENSG00000282130-0.3972810.00e+000.00e+000.1630.384CD8+ T cellspost

CLK1

ENSG00000013441-0.3979970.00e+000.00e+000.3360.549CD8+ T cellspost

TGFB1

ENSG00000105329-0.3981770.00e+000.00e+000.4730.55CD8+ T cellspost

WHRN

ENSG00000095397-0.4000190.00e+000.00e+000.1030.323CD8+ T cellspost

BPGM

ENSG00000172331-0.4004890.00e+000.00e+000.0840.268CD8+ T cellspost

G3BP2

ENSG00000138757-0.4032880.00e+000.00e+000.3020.51CD8+ T cellspost

PDE4B

ENSG00000184588-0.4036080.00e+000.00e+000.2620.455CD8+ T cellspost
Page: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

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.

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.)
Enrichment results for 1 known drug resistance mechanisms in TME cells

boxplot
MechnismPvalFDRSignatureGenesetOverlapBackgroundCell typeEnriched Genes
Regulation by the Disease Microenvironment8.40e-038.40e-035171861023467CD4+ T cellsBCL2,CDKN1B,NFKB1,TGFB1,CXCR4,FYN,SOCS3,STAT3,VIM,TSPYL2
Regulation by the Disease Microenvironment4.80e-054.80e-055001861423467CD8+ T cellsBCL2,CDKN1B,NFKB1,TGFB1,TGFBR2,CXCR4,FYN,GATA3,RHOB,HNRNPA2B1,SOCS3,STAT3,VIM,TSPYL2

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.)
Up-regulated HallmarkUp-regulated KEGGUp-regulated GO BP
boxplotboxplotboxplot

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.)
Down-regulated HallmarkDown-regulated KEGGDown-regulated GO BP
boxplotboxplotboxplot

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.)

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.)

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.)

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.)
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.This dataset does not contain malignant cells, nor does it predict any miRNA that corresponds to the top10 DEGs.

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.)

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.)
Motifs and TFs regulating up-regulated DEGs in resistant group of the post-treatment samples
motifmotifNESAUCTF_highConfCell_type
motifswissregulon__hs__GMEB24.760.091GMEB2 (directAnnotation). CD4+ T cells
motiftransfac_pro__M064083.910.0802ZSCAN30 (directAnnotation). CD4+ T cells
motifswissregulon__hs__ETV63.850.0794ETV6 (directAnnotation). CD4+ T cells
motiftransfac_pro__M006553.780.0785ETV4 (directAnnotation). CD4+ T cells
motiftransfac_pro__M038443.720.0778ELK4 (directAnnotation). CD4+ T cells
motiftaipale_tf_pairs__ERF_HOXA3_RSCGGAWNNNNNNNYMATTA_CAP_repr3.650.0769ERF; HOXA3 (directAnnotation). CD4+ T cells
motiftaipale_tf_pairs__TEAD4_CEBPD_NTTRCGYAANNNNNNRGWATGY_CAP3.650.0768CEBPD; TEAD4 (directAnnotation). CD4+ T cells
motiftaipale_tf_pairs__ELF2_NATGCGGAAGTR_HT3.640.0767ELF2 (directAnnotation). CD4+ T cells
motifjaspar__MA1509.13.60.0762IRF6 (directAnnotation). CD4+ T cells
motifjaspar__MA1939.13.50.075SREBF2 (directAnnotation). CD4+ T cells
Page: 1 2 3 4 5 6 7

Motifs and TFs regulating down-regulated DEGs in resistant group of the post-treatment samples
motifmotifNESAUCTF_highConfCell_type
motifjaspar__MA0604.18.190.12ATF1 (inferredBy_Orthology). CD4+ T cells
motiftaipale_cyt_meth__CREB1_NRTGACGTN_eDBD_repr6.460.103CREB1 (directAnnotation). CD4+ T cells
motiftransfac_pro__M018615.610.0939ATF1 (directAnnotation). CD4+ T cells
motifmetacluster_111.34.980.0874ATF1; ATF3; CREB1; CREB1; CREB1; CREB1; CREB1; CREM; CREM; CREM; CREM; E4F1; E4F1 (directAnnotation). CREB1; CREB1; CREM; JDP2 (inferredBy_Orthology). CD4+ T cells
motifdbtfbs__GMEB1_K562_ENCSR928KOR_merged_N14.860.0861GMEB1 (directAnnotation). CD4+ T cells
motiftransfac_pro__M008014.770.0852ATF1; ATF2; ATF3; ATF4; ATF7; CREB1; CREM (directAnnotation). ATF7-NPFF (inferredBy_Orthology). CD4+ T cells
motifdbtfbs__GMEB2_HepG2_ENCSR745VSQ_merged_N14.770.0852GMEB2 (directAnnotation). CD4+ T cells
motifmetacluster_111.44.730.0848ATF1; ATF1; ATF1; ATF1; ATF1; ATF1; ATF2; ATF2; ATF2; ATF2; ATF2; ATF2; ATF2; ATF2; ATF2; ATF2; ATF3; ATF3; ATF3; ATF3; ATF3; ATF4; ATF4; ATF5; ATF6; ATF6B; ATF7; ATF7; ATF7; BATF3; BATF3; CREB1; CREB1; CREB1; CREB1; CREB1; CREB1; CREB1; CREB1; CREB1; CREB1; CREB3; CREB3L4; CREB3L4; CREB5; CREB5; CREB5; CREM; CREM; CREM; CREM; CREM; FOSB; FOSL1; FOSL1; FOSL2; JDP2; JDP2; JDP2; JUN; JUNB; JUNB; JUNB; JUND; JUND; JUND; JUND; JUND; XBP1; XBP1; XBP1; ZNF629 (directAnnotation). ATF1; ATF1; ATF1; ATF2; ATF2; ATF3; ATF6; ATF6B; ATF7; ATF7; ATF7; ATF7-NPFF; ATF7-NPFF; ATF7-NPFF; ATF7-NPFF; BATF; BATF2; BATF3; CREB1; CREB1; CREB3; CREB3L1; CREB3L2; CREB3L3; CREB3L4; CREM; FOS; FOSB; FOSL1; FOSL2; JDP2; JDP2 (inferredBy_Orthology). CD4+ T cells
motifhocomoco__GMEB2_HUMAN.H11MO.0.D4.460.0821GMEB2 (directAnnotation). CD4+ T cells
motiftransfac_public__M000394.430.0817CREB1 (directAnnotation). CD4+ T cells
Page: 1 2 3 4 5 6 7 8 9 10 11

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.)

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

ETS1

cisbp__M02645upCD4+ T cellspost

PRDM1

metacluster_187.2upCD4+ T cellspost

ATF3

metacluster_111.3downCD4+ T cellspost

CREM

metacluster_111.3downCD4+ T cellspost

ATF3

transfac_pro__M00801downCD4+ T cellspost

CREM

transfac_pro__M00801downCD4+ T cellspost

ATF3

metacluster_111.4downCD4+ T cellspost

CREM

metacluster_111.4downCD4+ T cellspost

FOSB

metacluster_111.4downCD4+ T cellspost

FOSL2

metacluster_111.4downCD4+ T cellspost
Page: 1 2




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."