Comparative secretome analysis of Striga and Cuscuta species identifies candidate virulence factors for two evolutionarily independent parasitic plant lineages.
Cuscuta
Striga
Host-parasite interaction
Parasitic plants
Secretomes
Transcriptomics
Virulence factor
Journal
BMC plant biology
ISSN: 1471-2229
Titre abrégé: BMC Plant Biol
Pays: England
ID NLM: 100967807
Informations de publication
Date de publication:
06 Apr 2024
06 Apr 2024
Historique:
received:
21
07
2023
accepted:
20
03
2024
medline:
7
4
2024
pubmed:
7
4
2024
entrez:
6
4
2024
Statut:
epublish
Résumé
Many parasitic plants of the genera Striga and Cuscuta inflict huge agricultural damage worldwide. To form and maintain a connection with a host plant, parasitic plants deploy virulence factors (VFs) that interact with host biology. They possess a secretome that represents the complement of proteins secreted from cells and like other plant parasites such as fungi, bacteria or nematodes, some secreted proteins represent VFs crucial to successful host colonisation. Understanding the genome-wide complement of putative secreted proteins from parasitic plants, and their expression during host invasion, will advance understanding of virulence mechanisms used by parasitic plants to suppress/evade host immune responses and to establish and maintain a parasite-host interaction. We conducted a comparative analysis of the secretomes of root (Striga spp.) and shoot (Cuscuta spp.) parasitic plants, to enable prediction of candidate VFs. Using orthogroup clustering and protein domain analyses we identified gene families/functional annotations common to both Striga and Cuscuta species that were not present in their closest non-parasitic relatives (e.g. strictosidine synthase like enzymes), or specific to either the Striga or Cuscuta secretomes. For example, Striga secretomes were strongly associated with 'PAR1' protein domains. These were rare in the Cuscuta secretomes but an abundance of 'GMC oxidoreductase' domains were found, that were not present in the Striga secretomes. We then conducted transcriptional profiling of genes encoding putatively secreted proteins for the most agriculturally damaging root parasitic weed of cereals, S. hermonthica. A significant portion of the Striga-specific secretome set was differentially expressed during parasitism, which we probed further to identify genes following a 'wave-like' expression pattern peaking in the early penetration stage of infection. We identified 39 genes encoding putative VFs with functions such as cell wall modification, immune suppression, protease, kinase, or peroxidase activities, that are excellent candidates for future functional studies. Our study represents a comprehensive secretome analysis among parasitic plants and revealed both similarities and differences in candidate VFs between Striga and Cuscuta species. This knowledge is crucial for the development of new management strategies and delaying the evolution of virulence in parasitic weeds.
Sections du résumé
BACKGROUND
BACKGROUND
Many parasitic plants of the genera Striga and Cuscuta inflict huge agricultural damage worldwide. To form and maintain a connection with a host plant, parasitic plants deploy virulence factors (VFs) that interact with host biology. They possess a secretome that represents the complement of proteins secreted from cells and like other plant parasites such as fungi, bacteria or nematodes, some secreted proteins represent VFs crucial to successful host colonisation. Understanding the genome-wide complement of putative secreted proteins from parasitic plants, and their expression during host invasion, will advance understanding of virulence mechanisms used by parasitic plants to suppress/evade host immune responses and to establish and maintain a parasite-host interaction.
RESULTS
RESULTS
We conducted a comparative analysis of the secretomes of root (Striga spp.) and shoot (Cuscuta spp.) parasitic plants, to enable prediction of candidate VFs. Using orthogroup clustering and protein domain analyses we identified gene families/functional annotations common to both Striga and Cuscuta species that were not present in their closest non-parasitic relatives (e.g. strictosidine synthase like enzymes), or specific to either the Striga or Cuscuta secretomes. For example, Striga secretomes were strongly associated with 'PAR1' protein domains. These were rare in the Cuscuta secretomes but an abundance of 'GMC oxidoreductase' domains were found, that were not present in the Striga secretomes. We then conducted transcriptional profiling of genes encoding putatively secreted proteins for the most agriculturally damaging root parasitic weed of cereals, S. hermonthica. A significant portion of the Striga-specific secretome set was differentially expressed during parasitism, which we probed further to identify genes following a 'wave-like' expression pattern peaking in the early penetration stage of infection. We identified 39 genes encoding putative VFs with functions such as cell wall modification, immune suppression, protease, kinase, or peroxidase activities, that are excellent candidates for future functional studies.
CONCLUSIONS
CONCLUSIONS
Our study represents a comprehensive secretome analysis among parasitic plants and revealed both similarities and differences in candidate VFs between Striga and Cuscuta species. This knowledge is crucial for the development of new management strategies and delaying the evolution of virulence in parasitic weeds.
Identifiants
pubmed: 38582844
doi: 10.1186/s12870-024-04935-7
pii: 10.1186/s12870-024-04935-7
doi:
Types de publication
Journal Article
Langues
eng
Sous-ensembles de citation
IM
Pagination
251Subventions
Organisme : Biotechnology and Biological Sciences Research Council
ID : BB/J011703/1
Pays : United Kingdom
Organisme : Biotechnology and Biological Sciences Research Council
ID : BB/J011703/1
Pays : United Kingdom
Informations de copyright
© 2024. The Author(s).
Références
Westwood JH, Yoder JI, Timko MP. The evolution of parasitism in plants. Trends Plant Sci. 2010;15:227–35.
pubmed: 20153240
doi: 10.1016/j.tplants.2010.01.004
Parker C, Riches CR. Parasitic Weeds of the World: Biology and Control. CAB International. 1993.
Dorr I. How striga parasitizes its host: a TEM and SEM Study. Annuals of Botany. 1997;79:463–72.
doi: 10.1006/anbo.1996.0385
Spallek T, Mutuku M, Shirasu K. The genus Striga: a witch profile: profile of striga, the witchweed. Mol Plant Pathol. 2013;14:861–9.
pubmed: 23841683
pmcid: 6638688
doi: 10.1111/mpp.12058
Vaughn KC. Conversion of the searching hyphae of dodder into xylic and phloic hyphae: a cytochemical and immunocytochemical investigation. Int J Plant Sci. 2006;167:1099–114.
doi: 10.1086/507872
Runo S, Kuria EK. Habits of a highly successful cereal killer Striga. PLOS Pathogens. 2018;14:1006731.
doi: 10.1371/journal.ppat.1006731
Masanga J, Mwangi BN, Kibet W, Sagero P, Wamalwa M, Oduor R, et al. Physiological and ecological warnings that dodders pose an exigent threat to farmlands in Eastern Africa. Plant Physiol. 2021;185:1457–67.
pubmed: 33661304
pmcid: 8355486
doi: 10.1093/plphys/kiab034
Yoshida S, Cui S, Ichihashi Y, Shirasu K. The haustorium, a specialized invasive organ in parasitic plants. Annu Rev Plant Biol. 2016;67:643–67.
pubmed: 27128469
doi: 10.1146/annurev-arplant-043015-111702
Bleischwitz M, Albert M, Fuchsbauer H-L, Kaldenhoff R. Significance of Cuscutain, a cysteine protease from Cuscuta reflexa, in host-parasite interactions. BMC Plant Biol. 2010;10:227.
pubmed: 20964874
pmcid: 3017850
doi: 10.1186/1471-2229-10-227
Spallek T, Melnyk CW, Wakatake T, Zhang J, Sakamoto Y, Kiba T, et al. Interspecies hormonal control of host root morphology by parasitic plants. Proc Natl Acad Sci USA. 2017;114:5283–8.
pubmed: 28461500
pmcid: 5441792
doi: 10.1073/pnas.1619078114
Clarke CR, Timko MP, Yoder JI, Axtell MJ, Westwood JH. Molecular dialog between parasitic plants and their hosts. Ann Rev Phytopathol. 2019;57:13.1-13.21.
doi: 10.1146/annurev-phyto-082718-100043
Qin L, Kudla U, Roze EHA, Goverse A, Popeijus H, Nieuwland J, et al. A nematode expansin acting on plants. Nat Brief Commun. 2004;427:30.
Mitreva-Dautova M, Roze E, Overmars H, de Graaff L, Schots A, Helder J, et al. A symbiont-independent endo-1,4-β-Xylanase from the plant-parasitic nematode meloidogyne incognita. MPMI. 2006;19:521–9.
pubmed: 16673939
doi: 10.1094/MPMI-19-0521
Ben-Daniel B-H, Bar-Zvi D, TsrorLahkim L. Pectate lyase affects pathogenicity in natural isolates of Colletotrichum coccodes and in pelA gene-disrupted and gene-overexpressing mutant lines: pectate lyase affects C. coccodes pathogenicity. Molecular Plant Pathology. 2012;13:187–97.
pubmed: 21848609
doi: 10.1111/j.1364-3703.2011.00740.x
Win J, Chaparro-Garcia A, Belhaj K, Saunders DGO, Yoshida K, Dong S, et al. Effector biology of plant-associated organisms: concepts and perspectives. Cold Spring Harb Symp Quant Biol. 2012;77:235–47.
pubmed: 23223409
doi: 10.1101/sqb.2012.77.015933
Białas A, Zess EK, De la Concepcion JC, Franceschetti M, Pennington HG, Yoshida K, et al. Lessons in effector and NLR biology of plant-microbe systems. MPMI. 2018;31:34–45.
pubmed: 29144205
doi: 10.1094/MPMI-08-17-0196-FI
Cook DE, Mesarich CH, Thomma BPHJ. Understanding plant immunity as a surveillance system to detect invasion. Annu Rev Phytopathol. 2015;53:541–63.
pubmed: 26047564
doi: 10.1146/annurev-phyto-080614-120114
Kanyuka K, Rudd JJ. Cell surface immune receptors: the guardians of the plant’s extracellular spaces. Curr Opin Plant Biol. 2019;50:1–8.
pubmed: 30861483
pmcid: 6731392
doi: 10.1016/j.pbi.2019.02.005
Su C, Liu H, Wafula EK, Honaas L, de Pamphilis CW, Timko MP. SHR4z, a novel decoy effector from the haustorium of the parasitic weed Striga gesnerioides, suppresses host plant immunity. New Phytol. 2019;226:891–908.
pubmed: 31788811
pmcid: 7187149
doi: 10.1111/nph.16351
Qiu S, Bradley JM, Zhang P, Chaudhuri R, Blaxter M, Butlin RK, et al. Genome-enabled discovery of candidate virulence loci in Striga hermonthica, a devastating parasite of African cereal crops. New Phytol. 2022;236:622–38.
pubmed: 35699626
pmcid: 9795911
doi: 10.1111/nph.18305
Losner-Goshen D. Pectolytic activity by the haustorium of the parasitic plant Orobanche L. (Orobanchaceae) in host roots. Ann Botany. 1998;81:319–26.
doi: 10.1006/anbo.1997.0563
Honaas LA, Wafula EK, Yang Z, Der JP, Wickett NJ, Altman NS, et al. Functional genomics of a generalist parasitic plant: laser microdissection of host-parasite interface reveals host-specific patterns of parasite gene expression. BMC Plant Biol. 2013;13:9.
pubmed: 23302495
pmcid: 3636017
doi: 10.1186/1471-2229-13-9
Hegenauer V, Slaby P, Körner M, Bruckmüller JB, Burggraf R, Albert I, et al. The tomato receptor CuRe1 senses a cell wall protein to identify Cuscuta as a pathogen. Nat Commun. 2019;11:5299.
doi: 10.1038/s41467-020-19147-4
Brown NA, Antoniw J, Hammond-Kosack KE. The predicted secretome of the plant pathogenic fungus fusarium graminearum: a refined comparative analysis. PLoS ONE. 2012;7:e33731.
pubmed: 22493673
pmcid: 3320895
doi: 10.1371/journal.pone.0033731
Saunders DGO, Win J, Cano LM, Szabo LJ, Kamoun S, Raffaele S. Using hierarchical clustering of secreted protein families to classify and rank candidate effectors of rust fungi. PLoS ONE. 2012;7:e29847.
pubmed: 22238666
pmcid: 3253089
doi: 10.1371/journal.pone.0029847
Guyon K, Balagué C, Roby D, Raffaele S. Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum. BMC Genomics. 2014;15:336.
pubmed: 24886033
pmcid: 4039746
doi: 10.1186/1471-2164-15-336
Thorpe P, Mantelin S, Cock PJ, Blok VC, Eves-van den Akker S, et al. Genomic characterisation of the effector complement of the potato cyst nematode Globodera pallida. BMC Genomics. 2014;15:923.
pubmed: 25342461
pmcid: 4213498
doi: 10.1186/1471-2164-15-923
Gardner M, Dhroso A, Johnson N, Davis EL, Baum TJ, Korkin D, et al. Novel global effector mining from the transcriptome of early life stages of the soybean cyst nematode Heterodera glycines. Sci Rep. 2018;8:2505.
pubmed: 29410430
pmcid: 5802810
doi: 10.1038/s41598-018-20536-5
Catanzariti A-M, Dodds PN, Lawrence GJ, Ayliffe MA, Ellis JG. Haustorially expressed secreted proteins from flax rust are highly enriched for avirulence elicitors. Plant Cell. 2005;18:243–56.
pubmed: 16326930
doi: 10.1105/tpc.105.035980
Cotton JA, Lilley CJ, Jones LM, Kikuchi T, Reid AJ, Thorpe P, et al. The genome and life-stage specific transcriptomes of Globodera pallida elucidate key aspects of plant parasitism by a cyst nematode. Genome Biol. 2014;15:R43.
pubmed: 24580726
pmcid: 4054857
doi: 10.1186/gb-2014-15-3-r43
Vogel A, Schwacke R, Denton AK, Usadel B, Hollmann J, Fischer K, et al. Footprints of parasitism in the genome of the parasitic flowering plant Cuscuta campestris. Nat Commun. 2018;9:2515.
pubmed: 29955043
pmcid: 6023873
doi: 10.1038/s41467-018-04344-z
Estep MC, Gowda BS, Huang K, Timko MP, Bennetzen JL. Genomic characterisation for parasitic weeds of the genus Striga by sample sequence analysis. The Plant Genome. 2012;5:30–41.
doi: 10.3835/plantgenome2011.11.0031
Yang Z, Wafula EK, Honaas LA, Zhang H, Das M, Fernandez-Aparicio M, et al. Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty. Mol Biol Evol. 2015;32:767–90.
pubmed: 25534030
doi: 10.1093/molbev/msu343
Yoshida S, Kee YJ. Large-scale sequencing paves the way for genomic and genetic analyses in parasitic plants. Curr Opin Biotechnol. 2021;70:248–54.
pubmed: 34242992
doi: 10.1016/j.copbio.2021.06.011
Jaouannet M, Magliano M, Arguel MJ, Gourgues M, Evangelisti E, Abad P, et al. The root-knot nematode calreticulin Mi-CRT Is a Key effector in plant defense suppression. MPMI. 2013;26:97–105.
pubmed: 22857385
doi: 10.1094/MPMI-05-12-0130-R
Djamei A, Schipper K, Rabe F, Ghosh A, Vincon V, Kahnt J, et al. Metabolic priming by a secreted fungal effector. Nature. 2011;478:395–8.
pubmed: 21976020
doi: 10.1038/nature10454
Bracher D, Kutchen TM. Strictosidine synthase from Rauvolfia serpentina: analysis of the gene involved in indole alkaloid biosynthesis. Arch Biochem Biophys. 1992;294:717–23.
pubmed: 1567228
doi: 10.1016/0003-9861(92)90746-J
Zhang D, Qi J, Yue J, Huang J, Sun T, Li S, et al. Root parasitic plant orobanche aegyptiaca and shoot parasitic plant Cuscuta australis obtained Brassicaceae-specific strictosidine synthase-like genes by horizontal gene transfer. BMC Plant Biol. 2014;14:19.
pubmed: 24411025
pmcid: 3893544
doi: 10.1186/1471-2229-14-19
Smith JD, Woldemarium MG, Mescher MC, Jander G, De Moraes CM. Glucosinolates from host plants influence growth of the parasitic plant Cuscuta gronovii and its susceptibility to aphid feeding. Plant Physiol. 2016;172:181–97.
pubmed: 27482077
pmcid: 5074629
doi: 10.1104/pp.16.00613
Yoshida S, Kim S, Wafula EK, Tanskanen J, Kim Y-M, Honaas L, et al. Genome sequence of striga asiatica provides insight into the evolution of plant parasitism. Curr Biol. 2019;29:3041-3052.e4.
pubmed: 31522940
doi: 10.1016/j.cub.2019.07.086
Ogawa S, Wakatake T, Spallek T, Ishida J, Sano R, Kurata T, et al. Subtilase activity in intrusive cells mediates haustorium maturation in parasitic plants. Plant Physiol. 2021;185:1381–94.
pubmed: 33793894
doi: 10.1093/plphys/kiaa001
Klepikova AV, Kasianov AS, Gerasimov ES, Logacheva MD, Penin AA. A high resolution map of the Arabidopsis thaliana developmental transcriptome based on RNA-seq profiling. Plant J. 2016;88:1058–70.
pubmed: 27549386
doi: 10.1111/tpj.13312
Sützl L, Foley G, Gillam EMJ, Bodén M, Haltrich D. The GMC superfamily of oxidoreductases revisited: analysis and evolution of fungal GMC oxidoreductases. Biotechnol Biofuels. 2019;12:118.
pubmed: 31168323
pmcid: 6509819
doi: 10.1186/s13068-019-1457-0
Levasseur A, Drula E, Lombard V, Coutinho PM, Henrissat B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. Biotechnol Biofuels. 2013;6:41.
pubmed: 23514094
pmcid: 3620520
doi: 10.1186/1754-6834-6-41
Musser RO, Cipollini DF, Hum-Musser SM, Williams SA, Brown JK, Felton GW. Evidence that the caterpillar salivary enzyme glucose oxidase provides herbivore offense in solanaceous plants. Arch Insect Biochem Physiol. 2005;58:128–37.
pubmed: 15660363
doi: 10.1002/arch.20039
Rao SAK, Carolan JC, Wilkinson TL. Proteomic profiling of cereal aphid saliva reveals both ubiquitous and adaptive secreted proteins. PLoS ONE. 2013;8:e57413.
pubmed: 23460852
pmcid: 3584018
doi: 10.1371/journal.pone.0057413
Kurotani K, Wakatake T, Ichihashi Y, Okayasu K, Sawai Y, Ogawa S, et al. Host-parasite tissue adhesion by a secreted type of β-1,4-glucanase in the parasitic plant Phtheirospermum japonicum. Commun Biol. 2020;3:407.
pubmed: 32733024
pmcid: 7393376
doi: 10.1038/s42003-020-01143-5
Kubicek CP, Starr TL, Glass NL. Plant cell wall-degrading enzymes and their secretion in plant-pathogenic fungi. Annu Rev Phytopathol. 2014;52:427–51.
pubmed: 25001456
doi: 10.1146/annurev-phyto-102313-045831
Bacete L, Mélida H, Miedes E, Molina A. Plant cell wall-mediated immunity: cell wall changes trigger disease resistance responses. Plant J. 2018;93:614–36.
pubmed: 29266460
doi: 10.1111/tpj.13807
Zhao Z, Liu H, Wang C, Xu J-R. Comparative analysis of fungal genomes reveals different plant cell wall degrading capacity in fungi. BMC Genomics. 2013;14:274.
pubmed: 23617724
pmcid: 3652786
doi: 10.1186/1471-2164-14-274
Laohavisit A, Wakatake T, Ishihama N, Mulvey H, Takizawa K, Suzuki T, et al. Quinone perception in plants via leucine-rich-repeat receptor-like kinases. Nature. 2020;587:92–7.
pubmed: 32879491
doi: 10.1038/s41586-020-2655-4
Sun G, Xu Y, Liu H, Sun T, Zhang J, Hettenhausen C, et al. Large-scale gene losses underlie the genome evolution of parasitic plant Cuscuta australis. Nat Commun. 2018;9:2683.
pubmed: 29992948
pmcid: 6041341
doi: 10.1038/s41467-018-04721-8
Conn CE, Bythell-Douglas R, Neumann D, Yoshida S, Whittington B, Westwood JH, et al. Convergent evolution of strigolactone perception enabled host detection in parasitic plants. Science. 2015;349:540–3.
pubmed: 26228149
doi: 10.1126/science.aab1140
Toh S, Holbrook-Smith D, Stogios PJ, Onopriyenko O, Lumba S, Tsuchiya Y, et al. Structure-function analysis identifies highly sensitive strigolactone receptors in Striga. Science. 2015;350:203–7.
pubmed: 26450211
doi: 10.1126/science.aac9476
Luo S, Hu W, Wang Y, Liu B, Yan H, Xiang Y. Genome-wide identification, classification, and expression of phytocyanins in Populus trichocarpa. Planta. 2018;247:1133–48.
pubmed: 29383450
doi: 10.1007/s00425-018-2849-2
Reyt G, Chao Z, Flis P, Salas-González I, Castrillo G, Chao D-Y, et al. Uclacyanin proteins are required for lignified nanodomain formation within casparian strips. Curr Biol. 2020;30:4103-4111.e6.
pubmed: 32857976
pmcid: 7575197
doi: 10.1016/j.cub.2020.07.095
Denancé N, Szurek B, Noël LD. Emerging functions of nodulin-like proteins in non-nodulating plant species. Plant Cell Physiol. 2014;55:469–74.
pubmed: 24470637
doi: 10.1093/pcp/pct198
Bekal S, Niblack TL, Lambert KN. A chorismate mutase from the soybean cyst nematode heterodera glycines shows polymorphisms that correlate with virulence. MPMI. 2003;16:439–46.
pubmed: 12744515
doi: 10.1094/MPMI.2003.16.5.439
Finn RD, Attwood TK, Babbitt PC, Bateman A, Bork P, Bridge AJ, et al. InterPro in 2017—beyond protein family and domain annotations. Nucleic Acids Res. 2017;45:D190–9.
pubmed: 27899635
doi: 10.1093/nar/gkw1107
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.
pubmed: 24451626
pmcid: 3998142
doi: 10.1093/bioinformatics/btu031
Bendtsen DJ, Nielsen H, von Heijne G, Brunak S. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004;340:783–95.
pubmed: 15223320
doi: 10.1016/j.jmb.2004.05.028
Armenteros AJJ, Tsirigos KD, Sønderby CK, Petersen TN, Winther O, Brunak S, et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37:420–3.
doi: 10.1038/s41587-019-0036-z
Krogh A, Larsson B, von Heijne G, Sonnhammer ELL. Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.
pubmed: 11152613
doi: 10.1006/jmbi.2000.4315
R Core Team R: A language and environment for statistical computing. R Foundation for statistical computing, Vienna, Austria. 2022.
Emms DM, Kelly S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019;20:238.
pubmed: 31727128
pmcid: 6857279
doi: 10.1186/s13059-019-1832-y
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol Biol Evol. 2018;35:1547–9.
pubmed: 29722887
pmcid: 5967553
doi: 10.1093/molbev/msy096
Han MV, Thomas GWC, Lugo-Martinez J, Hahn MW. estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol Biol Evol. 2013;30:1987–97.
pubmed: 23709260
doi: 10.1093/molbev/mst100
Cissoko M, Boisnard A, Rodenburg J, Press MC, Scholes JD. New Rice for Africa (NERICA) cultivars exhibit different levels of post-attachment resistance against the parasitic weeds Striga hermonthica and Striga asiatica. New Phytol. 2011;192:952–63.
pubmed: 21883232
doi: 10.1111/j.1469-8137.2011.03846.x
Gurney AL, Slate J, Press MC, Scholes JD. A novel form of resistance in rice to the angiosperm parasite Striga hermonthica. New Phytol. 2006;169:199–208.
pubmed: 16390431
doi: 10.1111/j.1469-8137.2005.01560.x
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
pubmed: 25516281
pmcid: 4302049
doi: 10.1186/s13059-014-0550-8