Somatic Gain of KRAS Function in the Endothelium Is Sufficient to Cause Vascular Malformations That Require MEK but Not PI3K Signaling
Jason E Fish, Carlos Perfecto Flores Suarez, Emilie Boudreau, Alexander M Herman, Manuel Cantu Gutierrez, Dakota Gustafson, Peter V DiStefano, Meng Cui, Zhiqi Chen, Karen Berman De Ruiz, Taylor S Schexnayder, Christopher S Ward, Ivan Radovanovic, Joshua D Wythe, Jason E Fish, Carlos Perfecto Flores Suarez, Emilie Boudreau, Alexander M Herman, Manuel Cantu Gutierrez, Dakota Gustafson, Peter V DiStefano, Meng Cui, Zhiqi Chen, Karen Berman De Ruiz, Taylor S Schexnayder, Christopher S Ward, Ivan Radovanovic, Joshua D Wythe
Abstract
Rationale: We previously identified somatic activating mutations in the KRAS (Kirsten rat sarcoma viral oncogene homologue) gene in the endothelium of the majority of human sporadic brain arteriovenous malformations; a disorder characterized by direct connections between arteries and veins. However, whether this genetic abnormality alone is sufficient for lesion formation, as well as how active KRAS signaling contributes to arteriovenous malformations, remains unknown.
Objective: To establish the first in vivo models of somatic KRAS gain of function in the endothelium in both mice and zebrafish to directly observe the phenotypic consequences of constitutive KRAS activity at a cellular level in vivo, and to test potential therapeutic interventions for arteriovenous malformations.
Methods and results: Using both postnatal and adult mice, as well as embryonic zebrafish, we demonstrate that endothelial-specific gain of function mutations in Kras (G12D or G12V) are sufficient to induce brain arteriovenous malformations. Active KRAS signaling leads to altered endothelial cell morphogenesis and increased cell size, ectopic sprouting, expanded vessel lumen diameter, and direct connections between arteries and veins. Furthermore, we show that these lesions are not associated with altered endothelial growth dynamics or a lack of proper arteriovenous identity but instead seem to feature exuberant angiogenic signaling. Finally, we demonstrate that KRAS-dependent arteriovenous malformations in zebrafish are refractory to inhibition of the downstream effector PI3K but instead require active MEK (mitogen-activated protein kinase kinase 1) signaling.
Conclusions: We demonstrate that active KRAS expression in the endothelium is sufficient for brain arteriovenous malformations, even in the setting of uninjured adult vasculature. Furthermore, the finding that KRAS-dependent lesions are reversible in zebrafish suggests that MEK inhibition may represent a promising therapeutic treatment for arteriovenous malformation patients. Graphical Abstract: A graphical abstract is available for this article.
Keywords: brain; cell size; endothelium, vascular; models, animal; vascular disease.
Conflict of interest statement
None.
Figures
References
- Solomon RA, Connolly ES., Jr. Arteriovenous malformations of the brain. N Engl J Med 20173761859–1866doi: 10.1056/NEJMra1607407
- Yun JH, Kwon DH, Lee EJ, Lee DH, Ahn JS, Kwun BD. New nidus formation adjacent to the target site of an arteriovenous malformation treated by Gamma Knife surgery. J Neurosurg 2012117suppl120–125doi: 10.3171/2012.8.GKS12994
- Nikolaev SI, Vetiska S, Bonilla X, Boudreau E, Jauhiainen S, Rezai Jahromi B, Khyzha N, DiStefano PV, Suutarinen S, Kiehl TR, et al. Somatic activating KRAS mutations in arteriovenous malformations of the brain. N Engl J Med 2018378250–261doi: 10.1056/NEJMoa1709449
- Oka M, Kushamae M, Aoki T, Yamaguchi T, Kitazato K, Abekura Y, Kawamata T, Mizutani T, Miyamoto S, Takagi Y. KRAS G12D or G12V mutation in human brain arteriovenous malformations. World Neurosurg 2019126e1365–e1373doi: 10.1016/j.wneu.2019.03.105
- Priemer DS, Vortmeyer AO, Zhang S, Chang HY, Curless KL, Cheng L. Activating KRAS mutations in arteriovenous malformations of the brain: frequency and clinicopathologic correlation. Hum Pathol 20198933–39doi: 10.1016/j.humpath.2019.04.004
- Hong T, Yan Y, Li J, Radovanovic I, Ma X, Shao YW, Yu J, Ma Y, Zhang P, Ling F, et al. High prevalence of KRAS/BRAF somatic mutations in brain and spinal cord arteriovenous malformations. Brain 201914223–34doi: 10.1093/brain/awy307
- Al-Olabi L, Polubothu S, Dowsett K, Andrews KA, Stadnik P, Joseph AP, Knox R, Pittman A, Clark G, Baird W, et al. Mosaic RAS/MAPK variants cause sporadic vascular malformations which respond to targeted therapy. J Clin Invest 20181281496–1508doi: 10.1172/JCI98589
- Serebriiskii IG, Connelly C, Frampton G, Newberg J, Cooke M, Miller V, Ali S, Ross JS, Handorf E, Arora S, et al. Comprehensive characterization of RAS mutations in colon and rectal cancers in old and young patients. Nat Commun. 2019;10:3722. doi: 10.1038/s41467-019-11530-0.
- Lapinski PE, Doosti A, Salato V, North P, Burrows PE, King PD. Somatic second hit mutation of RASA1 in vascular endothelial cells in capillary malformation-arteriovenous malformation. Eur J Med Genet 20186111–16doi: 10.1016/j.ejmg.2017.10.004
- Limaye N, Wouters V, Uebelhoer M, Tuominen M, Wirkkala R, Mulliken JB, Eklund L, Boon LM, Vikkula M. Somatic mutations in angiopoietin receptor gene TEK cause solitary and multiple sporadic venous malformations. Nat Genet 200941118–124doi: 10.1038/ng.272
- Limaye N, Kangas J, Mendola A, Godfraind C, Schlögel MJ, Helaers R, Eklund L, Boon LM, Vikkula M. Somatic activating PIK3CA mutations cause venous malformation. Am J Hum Genet 201597914–921doi: 10.1016/j.ajhg.2015.11.011
- Shirley MD, Tang H, Gallione CJ, Baugher JD, Frelin LP, Cohen B, North PE, Marchuk DA, Comi AM, Pevsner J. Sturge-Weber syndrome and port-wine stains caused by somatic mutation in GNAQ. N Engl J Med 20133681971–1979doi: 10.1056/NEJMoa1213507
- Huang L, Couto JA, Pinto A, Alexandrescu S, Madsen JR, Greene AK, Sahin M, Bischoff J. Somatic GNAQ mutation is enriched in brain endothelial cells in Sturge-Weber syndrome. Pediatr Neurol 20176759–63doi: 10.1016/j.pediatrneurol.2016.10.010
- Hammer J, Seront E, Duez S, Dupont S, Van Damme A, Schmitz S, Hoyoux C, Chopinet C, Clapuyt P, Hammer F, et al. Sirolimus is efficacious in treatment for extensive and/or complex slow-flow vascular malformations: a monocentric prospective phase II study. Orphanet J Rare Dis. 2018;13:191. doi: 10.1186/s13023-018-0934-z.
- Li X, Cai Y, Goines J, Pastura P, Brichta L, Lane A, Le Cras TD, Boscolo E. Ponatinib combined with rapamycin causes regression of murine venous malformation. Arterioscler Thromb Vasc Biol 201939496–512doi: 10.1161/ATVBAHA.118.312315
- Boscolo E, Limaye N, Huang L, Kang KT, Soblet J, Uebelhoer M, Mendola A, Natynki M, Seront E, Dupont S, et al. Rapamycin improves TIE2-mutated venous malformation in murine model and human subjects. J Clin Invest 20151253491–3504doi: 10.1172/JCI76004
- Hong CC, Peterson QP, Hong JY, Peterson RT. Artery/vein specification is governed by opposing phosphatidylinositol-3 kinase and MAP kinase/ERK signaling. Curr Biol 2006161366–1372doi: 10.1016/j.cub.2006.05.046
- Deng Y, Larrivée B, Zhuang ZW, Atri D, Moraes F, Prahst C, Eichmann A, Simons M. Endothelial RAF1/ERK activation regulates arterial morphogenesis. Blood 20131213988–96, S1doi: 10.1182/blood-2012-12-474601
- Wythe JD, Dang LT, Devine WP, Boudreau E, Artap ST, He D, Schachterle W, Stainier DY, Oettgen P, Black BL, et al. ETS factors regulate Vegf-dependent arterial specification. Dev Cell 20132645–58doi: 10.1016/j.devcel.2013.06.007
- Shin M, Beane TJ, Quillien A, Male I, Zhu LJ, Lawson ND. Vegfa signals through ERK to promote angiogenesis, but not artery differentiation. Development 20161433796–3805doi: 10.1242/dev.137919
- Fish JE, Cantu Gutierrez M, Dang LT, Khyzha N, Chen Z, Veitch S, Cheng HS, Khor M, Antounians L, Njock MS, et al. Dynamic regulation of VEGF-inducible genes by an ERK/ERG/p300 transcriptional network. Development 20171442428–2444doi: 10.1242/dev.146050
- Yaeger R, Corcoran RB. Targeting alterations in the RAF-MEK pathway. Cancer Discov 20199329–341doi: 10.1158/-18-1321
- Wang Y, Nakayama M, Pitulescu ME, Schmidt TS, Bochenek ML, Sakakibara A, Adams S, Davy A, Deutsch U, Lüthi U, et al. Ephrin-B2 controls VEGF-induced angiogenesis and lymphangiogenesis. Nature 2010465483–486doi: 10.1038/nature09002
- Jackson EL, Willis N, Mercer K, Bronson RT, Crowley D, Montoya R, Jacks T, Tuveson DA. Analysis of lung tumor initiation and progression using conditional expression of oncogenic K-ras. Genes Dev 2001153243–3248doi: 10.1101/gad.943001
- Murphy PA, Lam MT, Wu X, Kim TN, Vartanian SM, Bollen AW, Carlson TR, Wang RA. Endothelial Notch4 signaling induces hallmarks of brain arteriovenous malformations in mice. Proc Natl Acad Sci U S A 200810510901–10906doi: 10.1073/pnas.0802743105
- Ridder DA, Lang MF, Salinin S, Röderer JP, Struss M, Maser-Gluth C, Schwaninger M. TAK1 in brain endothelial cells mediates fever and lethargy. J Exp Med 20112082615–2623doi: 10.1084/jem.20110398
- Mahmoud M, Allinson KR, Zhai Z, Oakenfull R, Ghandi P, Adams RH, Fruttiger M, Arthur HM. Pathogenesis of arteriovenous malformations in the absence of endoglin. Circ Res 20101061425–1433doi: 10.1161/CIRCRESAHA.109.211037
- Tual-Chalot S, Mahmoud M, Allinson KR, Redgrave RE, Zhai Z, Oh SP, Fruttiger M, Arthur HM. Endothelial depletion of Acvrl1 in mice leads to arteriovenous malformations associated with reduced endoglin expression. PLoS One. 2014;9:e98646. doi: 10.1371/journal.pone.0098646.
- Nielsen CM, Huang L, Murphy PA, Lawton MT, Wang RA. Mouse Models of Cerebral Arteriovenous Malformation. Stroke 201647293–300doi: 10.1161/STROKEAHA.115.002869
- Park SO, Wankhede M, Lee YJ, Choi EJ, Fliess N, Choe SW, Oh SH, Walter G, Raizada MK, Sorg BS, et al. Real-time imaging of de novo arteriovenous malformation in a mouse model of hereditary hemorrhagic telangiectasia. J Clin Invest 20091193487–3496doi: 10.1172/JCI39482
- Chen W, Sun Z, Han Z, Jun K, Camus M, Wankhede M, Mao L, Arnold T, Young WL, Su H. De novo cerebrovascular malformation in the adult mouse after endothelial Alk1 deletion and angiogenic stimulation. Stroke 201445900–902doi: 10.1161/STROKEAHA.113.003655
- Leblanc GG, Golanov E, Awad IA, Young WL; Biology of Vascular Malformations of the Brain NINDS Workshop Collaborators Biology of vascular malformations of the brain. Stroke 200940e694–e702doi: 10.1161/STROKEAHA.109.563692
- Lawton MT, Rutledge WC, Kim H, Stapf C, Whitehead KJ, Li DY, Krings T, terBrugge K, Kondziolka D, Morgan MK, et al. Brain arteriovenous malformations. Nat Rev Dis Primers. 2015;1:15008. doi: 10.1038/nrdp.2015.8.
- Xu B, Wu YQ, Huey M, Arthur HM, Marchuk DA, Hashimoto T, Young WL, Yang GY. Vascular endothelial growth factor induces abnormal microvasculature in the endoglin heterozygous mouse brain. J Cereb Blood Flow Metab 200424237–244doi: 10.1097/01.WCB.0000107730.66603.51
- Hao Q, Zhu Y, Su H, Shen F, Yang GY, Kim H, Young WL. VEGF induces more severe cerebrovascular dysplasia in endoglin than in Alk1 mice. Transl Stroke Res 20101197–201doi: 10.1007/s12975-010-0020-x
- Choi EJ, Walker EJ, Shen F, Oh SP, Arthur HM, Young WL, Su H. Minimal homozygous endothelial deletion of Eng with VEGF stimulation is sufficient to cause cerebrovascular dysplasia in the adult mouse. Cerebrovasc Dis 201233540–547doi: 10.1159/000337762
- Sugden WW, Meissner R, Aegerter-Wilmsen T, Tsaryk R, Leonard EV, Bussmann J, Hamm MJ, Herzog W, Jin Y, Jakobsson L, et al. Endoglin controls blood vessel diameter through endothelial cell shape changes in response to haemodynamic cues. Nat Cell Biol 201719653–665doi: 10.1038/ncb3528
- Rochon ER, Menon PG, Roman BL. Alk1 controls arterial endothelial cell migration in lumenized vessels. Development 20161432593–2602doi: 10.1242/dev.135392
- Eerola I, Boon LM, Mulliken JB, Burrows PE, Dompmartin A, Watanabe S, Vanwijck R, Vikkula M. Capillary malformation-arteriovenous malformation, a new clinical and genetic disorder caused by RASA1 mutations. Am J Hum Genet 2003731240–1249doi: 10.1086/379793
- Revencu N, Boon LM, Mulliken JB, Enjolras O, Cordisco MR, Burrows PE, Clapuyt P, Hammer F, Dubois J, Baselga E, et al. Parkes Weber syndrome, vein of Galen aneurysmal malformation, and other fast-flow vascular anomalies are caused by RASA1 mutations. Hum Mutat 200829959–965doi: 10.1002/humu.20746
- Chen D, Teng JM, North PE, Lapinski PE, King PD. RASA1-dependent cellular export of collagen IV controls blood and lymphatic vascular development. J Clin Invest 20191293545–3561doi: 10.1172/JCI124917
- Henkemeyer M, Rossi DJ, Holmyard DP, Puri MC, Mbamalu G, Harpal K, Shih TS, Jacks T, Pawson T. Vascular system defects and neuronal apoptosis in mice lacking ras GTPase-activating protein. Nature 1995377695–701doi: 10.1038/377695a0
- Lapinski PE, Kwon S, Lubeck BA, Wilkinson JE, Srinivasan RS, Sevick-Muraca E, King PD. RASA1 maintains the lymphatic vasculature in a quiescent functional state in mice. J Clin Invest 2012122733–747doi: 10.1172/JCI46116
- Lapinski PE, Lubeck BA, Chen D, Doosti A, Zawieja SD, Davis MJ, King PD. RASA1 regulates the function of lymphatic vessel valves in mice. J Clin Invest 20171272569–2585doi: 10.1172/JCI89607
- Lubeck BA, Lapinski PE, Bauler TJ, Oliver JA, Hughes ED, Saunders TL, King PD. Blood vascular abnormalities in Rasa1(R780Q) knockin mice: implications for the pathogenesis of capillary malformation-arteriovenous malformation. Am J Pathol 20141843163–3169doi: 10.1016/j.ajpath.2014.08.018
- Barclay SF, Inman KW, Luks VL, McIntyre JB, Al-Ibraheemi A, Church AJ, Perez-Atayde AR, Mangray S, Jeng M, Kreimer SR, et al. A somatic activating NRAS variant associated with kaposiform lymphangiomatosis. Genet Med 2019211517–1524doi: 10.1038/s41436-018-0390-0
- Manevitz-Mendelson E, Leichner GS, Barel O, Davidi-Avrahami I, Ziv-Strasser L, Eyal E, Pessach I, Rimon U, Barzilai A, Hirshberg A, et al. Somatic NRAS mutation in patient with generalized lymphatic anomaly. Angiogenesis 201821287–298doi: 10.1007/s10456-018-9595-8
- Konczyk DJ, Goss JA, Smits PJ, Huang AY, Al-Ibraheemi A, Sudduth CL, Warman ML, Greene AK. Arteriovenous malformation associated with a HRAS mutation. Hum Genet 20191381419–1421doi: 10.1007/s00439-019-02072-y
- Li QF, Decker-Rockefeller B, Bajaj A, Pumiglia K. Activation of Ras in the vascular endothelium induces brain vascular malformations and hemorrhagic stroke. Cell Rep 2018242869–2882doi: 10.1016/j.celrep.2018.08.025
- Martinelli E, Morgillo F, Troiani T, Ciardiello F. Cancer resistance to therapies against the EGFR-RAS-RAF pathway: the role of MEK. Cancer Treat Rev 20175361–69doi: 10.1016/j.ctrv.2016.12.001
- Merchant M, Moffat J, Schaefer G, Chan J, Wang X, Orr C, Cheng J, Hunsaker T, Shao L, Wang SJ, et al. Combined MEK and ERK inhibition overcomes therapy-mediated pathway reactivation in RAS mutant tumors. PLoS One. 2017;12:e0185862. doi: 10.1371/journal.pone.0185862.
- Jessen WJ, Miller SJ, Jousma E, Wu J, Rizvi TA, Brundage ME, Eaves D, Widemann B, Kim MO, Dombi E, et al. MEK inhibition exhibits efficacy in human and mouse neurofibromatosis tumors. J Clin Invest 2013123340–347doi: 10.1172/JCI60578
- Wang MT, Fer N, Galeas J, Collisson EA, Kim SE, Sharib J, McCormick F. Blockade of leukemia inhibitory factor as a therapeutic approach to KRAS driven pancreatic cancer. Nat Commun. 2019;10:3055. doi: 10.1038/s41467-019-11044-9.
- Tetsu O, McCormick F. ETS-targeted therapy: can it substitute for MEK inhibitors? Clin Transl Med. 2017;6:16. doi: 10.1186/s40169-017-0147-4.
- Luo J, Solimini NL, Elledge SJ. Principles of cancer therapy: oncogene and non-oncogene addiction. Cell 2009136823–837doi: 10.1016/j.cell.2009.02.024
- Li D, March ME, Gutierrez-Uzquiza A, Kao C, Seiler C, Pinto E, Matsuoka LS, Battig MR, Bhoj EJ, Wenger TL, et al. ARAF recurrent mutation causes central conducting lymphatic anomaly treatable with a MEK inhibitor. Nat Med 2019251116–1122doi: 10.1038/s41591-019-0479-2
- Triana P, Dore M, Cerezo VN, Cervantes M, Sánchez AV, Ferrero MM, González MD, Lopez-Gutierrez JC. Sirolimus in the treatment of vascular anomalies. Eur J Pediatr Surg 20172786–90doi: 10.1055/s-0036-1593383
- Lekwuttikarn R, Lim YH, Admani S, Choate KA, Teng JMC. Genotype-guided medical treatment of an arteriovenous malformation in a child. JAMA Dermatol 2019155256–257doi: 10.1001/jamadermatol.2018.4653
- Villefranc JA, Amigo J, Lawson ND. Gateway compatible vectors for analysis of gene function in the zebrafish. Dev Dyn 20072363077–3087doi: 10.1002/dvdy.21354
- Zhou Y, Cashman TJ, Nevis KR, Obregon P, Carney SA, Liu Y, Gu A, Mosimann C, Sondalle S, Peterson RE, et al. Latent TGF-β binding protein 3 identifies a second heart field in zebrafish. Nature 2011474645–648doi: 10.1038/nature10094
- Nagai T, Ibata K, Park ES, Kubota M, Mikoshiba K, Miyawaki A. A variant of yellow fluorescent protein with fast and efficient maturation for cell-biological applications. Nat Biotechnol 20022087–90doi: 10.1038/nbt0102-87
- Zhao L, Borikova AL, Ben-Yair R, Guner-Ataman B, MacRae CA, Lee RT, Burns CG, Burns CE. Notch signaling regulates cardiomyocyte proliferation during zebrafish heart regeneration. Proc Natl Acad Sci U S A 20141111403–1408doi: 10.1073/pnas.1311705111
- Kwan KM, Fujimoto E, Grabher C, Mangum BD, Hardy ME, Campbell DS, Parant JM, Yost HJ, Kanki JP, Chien CB. The Tol2kit: a multisite gateway-based construction kit for Tol2 transposon transgenesis constructs. Dev Dyn 20072363088–3099doi: 10.1002/dvdy.21343
- Birnbaum RY, Clowney EJ, Agamy O, Kim MJ, Zhao J, Yamanaka T, Pappalardo Z, Clarke SL, Wenger AM, Nguyen L, et al. Coding exons function as tissue-specific enhancers of nearby genes. Genome Res 2012221059–1068doi: 10.1101/gr.133546.111
- Mo A, Mukamel EA, Davis FP, Luo C, Henry GL, Picard S, Urich MA, Nery JR, Sejnowski TJ, Lister R, et al. Epigenomic signatures of neuronal diversity in the mammalian brain. Neuron 2015861369–1384doi: 10.1016/j.neuron.2015.05.018
- Liu J, Krautzberger AM, Sui SH, Hofmann OM, Chen Y, Baetscher M, Grgic I, Kumar S, Humphreys BD, Humphreys B, et al. Cell-specific translational profiling in acute kidney injury. J Clin Invest 20141241242–1254doi: 10.1172/JCI72126
- Madisen L, Zwingman TA, Sunkin SM, Oh SW, Zariwala HA, Gu H, Ng LL, Palmiter RD, Hawrylycz MJ, Jones AR, et al. A robust and high-throughput Cre reporting and characterization system for the whole mouse brain. Nat Neurosci 201013133–140doi: 10.1038/nn.2467
- Hayashi S, McMahon AP. Efficient recombination in diverse tissues by a tamoxifen-inducible form of Cre: a tool for temporally regulated gene activation/inactivation in the mouse. Dev Biol 2002244305–318doi: 10.1006/dbio.2002.0597
- Zudaire E, Gambardella L, Kurcz C, Vermeren S. A computational tool for quantitative analysis of vascular networks. PLoS One. 2011;6:e27385. doi: 10.1371/journal.pone.0027385.
- Respress JL, Wehrens XH. Transthoracic echocardiography in mice. J Vis Exp. 2010:1738.
- Jin SW, Beis D, Mitchell T, Chen JN, Stainier DY. Cellular and molecular analyses of vascular tube and lumen formation in zebrafish. Development 20051325199–5209doi: 10.1242/dev.02087
- Traver D, Paw BH, Poss KD, Penberthy WT, Lin S, Zon LI. Transplantation and in vivo imaging of multilineage engraftment in zebrafish bloodless mutants. Nat Immunol 200341238–1246doi: 10.1038/ni1007
- Proulx K, Lu A, Sumanas S. Cranial vasculature in zebrafish forms by angioblast cluster-derived angiogenesis. Dev Biol 201034834–46doi: 10.1016/j.ydbio.2010.08.036
- Parsons MJ, Pisharath H, Yusuff S, Moore JC, Siekmann AF, Lawson N, Leach SD. Notch-responsive cells initiate the secondary transition in larval zebrafish pancreas. Mech Dev 2009126898–912doi: 10.1016/j.mod.2009.07.002
- Vanhollebeke B, Stone OA, Bostaille N, Cho C, Zhou Y, Maquet E, Gauquier A, Cabochette P, Fukuhara S, Mochizuki N, et al. Tip cell-specific requirement for an atypical Gpr124- and Reck-dependent Wnt/beta-catenin pathway during brain angiogenesis. Elife. 2015;4:e06489.
- Roman BL, Pham VN, Lawson ND, Kulik M, Childs S, Lekven AC, Garrity DM, Moon RT, Fishman MC, Lechleider RJ, et al. Disruption of acvrl1 increases endothelial cell number in zebrafish cranial vessels. Development 20021293009–3019
- Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods 201714417–419doi: 10.1038/nmeth.4197
- Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. 2015;4:1521. doi: 10.12688/f1000research.7563.2.
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. doi: 10.1186/s13059-014-0550-8.
- Blighe K, Rana S, Lewis M. EnhancedVolcano: Publication-ready volcano plots with enhanced colouring and labeling. R package version 152. . 2020.
- Kolde R. pheatmap: Pretty Heatmaps. R package version 1012. . 2019.
- Ge SX, Jung D, Yao R. ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics 2020362628–2629doi: 10.1093/bioinformatics/btz931
Source: PubMed