A systematic regulatory network related to bulbil formation in Lilium lancifolium based on metabolome and transcriptome analyses.
Lilium lancifolium
Bulbil formation
Hormone
Transcription factors
Transcriptome
Journal
BMC plant biology
ISSN: 1471-2229
Titre abrégé: BMC Plant Biol
Pays: England
ID NLM: 100967807
Informations de publication
Date de publication:
16 Oct 2024
16 Oct 2024
Historique:
received:
25
06
2024
accepted:
30
09
2024
medline:
16
10
2024
pubmed:
16
10
2024
entrez:
15
10
2024
Statut:
epublish
Résumé
Lilium lancifolium is a special wild triploid species native to China and can produce abundant bulbils on its stem under natural conditions, which is very valuable to study bulbil organogenesis in plants. Although similar to the lateral and tillering principles, the molecular mechanism underlying bulbil formation has remained incompletely understood. The metabolome and transcriptome of L. lancifolium bulbils across four development stages were analyzed. The pairwise comparison of metabolomes across the four stages identified 17 differential hormones, predominantly auxin (IAA), cytokinin (CK), and jasmonic acid (JA). Short Time-series Expression Miner (STEM) trend analysis of differential genes revealed four significant trends across these stages. The KEGG enrichment analysis of the four clusters highlighted pathways, such as plant hormone signal transduction, which were speculated to play a crucial role in development stages. these pathways were speculated to play a crucial role in development stages. To explore the key differential expressed genes and transcription factors associated with bulbil occurrence, two periods were focused on: Ll_UN and Ll_DN, which represented the stages with and without bulbils, respectively. Through correlation analysis and qRT-PCR analysis, 11 candidate differentially expressed genes and 27 candidate transcription factors were selected. By spraying exogenous hormones to validate these candidates, LlbHLH128, LlTIFY10A, LlbHLH93, and LlMYB108, were identified as the key genes for L. lancifolium bulbils. A regulatory network of L. lancifolium bulbil development was predicted. LlTIFY10A and LlbHLH93 might be involved in the JA and auxin signal transduction pathways, which jointly formed a regulatory network to affect the occurrence of L. lancifolium bulbil. This study not only provided more information about the differentially expressed genes and metabolites through transcriptome and metabolomics analyses, but also provided a clearer understanding of the effect of hormones on bulbil formation in lily.
Sections du résumé
BACKGROUND
BACKGROUND
Lilium lancifolium is a special wild triploid species native to China and can produce abundant bulbils on its stem under natural conditions, which is very valuable to study bulbil organogenesis in plants. Although similar to the lateral and tillering principles, the molecular mechanism underlying bulbil formation has remained incompletely understood.
RESULTS
RESULTS
The metabolome and transcriptome of L. lancifolium bulbils across four development stages were analyzed. The pairwise comparison of metabolomes across the four stages identified 17 differential hormones, predominantly auxin (IAA), cytokinin (CK), and jasmonic acid (JA). Short Time-series Expression Miner (STEM) trend analysis of differential genes revealed four significant trends across these stages. The KEGG enrichment analysis of the four clusters highlighted pathways, such as plant hormone signal transduction, which were speculated to play a crucial role in development stages. these pathways were speculated to play a crucial role in development stages. To explore the key differential expressed genes and transcription factors associated with bulbil occurrence, two periods were focused on: Ll_UN and Ll_DN, which represented the stages with and without bulbils, respectively. Through correlation analysis and qRT-PCR analysis, 11 candidate differentially expressed genes and 27 candidate transcription factors were selected. By spraying exogenous hormones to validate these candidates, LlbHLH128, LlTIFY10A, LlbHLH93, and LlMYB108, were identified as the key genes for L. lancifolium bulbils.
CONCLUSION
CONCLUSIONS
A regulatory network of L. lancifolium bulbil development was predicted. LlTIFY10A and LlbHLH93 might be involved in the JA and auxin signal transduction pathways, which jointly formed a regulatory network to affect the occurrence of L. lancifolium bulbil. This study not only provided more information about the differentially expressed genes and metabolites through transcriptome and metabolomics analyses, but also provided a clearer understanding of the effect of hormones on bulbil formation in lily.
Identifiants
pubmed: 39407139
doi: 10.1186/s12870-024-05654-9
pii: 10.1186/s12870-024-05654-9
doi:
Substances chimiques
Plant Growth Regulators
0
Types de publication
Journal Article
Langues
eng
Sous-ensembles de citation
IM
Pagination
969Subventions
Organisme : Liaoning Lily Resource Repository, Liaoning Department of Agriculture and Rural Affairs; Lily Germplasm Resource Nursery
ID : 000000
Organisme : Liaoning Lily Resource Repository, Liaoning Department of Agriculture and Rural Affairs; Lily Germplasm Resource Nursery
ID : 000000
Informations de copyright
© 2024. The Author(s).
Références
Jiang FX, Huang YX, Wang L, Zhou P, Sun DJ, Li X, Chen QB. Exploring the mystery of plant bulblet. Mol Plant Breed. 2017;15(1):346–52. https://doi.org/10.13271/j.mpb.015.000346 .
doi: 10.13271/j.mpb.015.000346
Yang PP, Xu LF, Xu H, Tang YC, He GR, Cao YW, Feng YY, Yuan SX, Ming J. Histological and transcriptomic analysis during bulbil formation in Lilium lancifolium. Front Plant Sci. 2017;8:1508. https://doi.org/10.3389/fpls.2017.01508 .
doi: 10.3389/fpls.2017.01508
Fan JP, Bing W, Yan FX, Liu JS, Chen MX. Study on morphological characteristics and physiological changes of bulbils of Lilium lancifolium. J Northeast Agricultural Univ. 2019;50(2):18–27. https://doi.org/10.19720/j.cnki.issn.1005-9369.2019.2.0003 .
doi: 10.19720/j.cnki.issn.1005-9369.2019.2.0003
Bennett T, Leyser O. Something on the side: axillary meristems and plant development. Plant Mol Biol. 2006;60(6):843–54. https://doi.org/10.1007/s11103-005-2763-4 .
doi: 10.1007/s11103-005-2763-4
Wang Q, Hasson A, Rossmann S, Theres K. Divide et impera: boundaries shape the plant body and initiate new meristems. New Phytol. 2016;209(2):485–98. https://doi.org/10.1111/nph.13641 .
doi: 10.1111/nph.13641
Luo R, Tian C, Xia LS, Chen HL, Du SS, Luo X. Cloning and expression analysis of PteAUX1 gene related to bulblet development in Pinellia ternata. Molecular Plant Breeding: 2022; 1–19. https://kns.cnki.net/kcms/detail/46.1068.S.20210909.1656.026.html
He GR, Yang PP, Tang YC, Cao YW, Qi XY, Xu LF, Ming J. Mechanism of exogenous cytokinins inducing bulbil formation in Lilium lancifolium in vitro. Plant Cell Rep. 2020;39(1):861–72. https://doi.org/10.1007/s00299-020-02535-x .
doi: 10.1007/s00299-020-02535-x
He GR, Yang PP, Cao YW, Tang YC, Wang L, Song M, Wang J, Xu LF, Ming J. Cytokinin type-B response regulators promote bulbil initiation in Lilium lancifolium. Int J Mol Sci. 2021;22(7):3320. https://doi.org/10.3390/ijms22073320 .
doi: 10.3390/ijms22073320
He GR, Cao YW, Wang J, Song M, Bi MM, Tang YH, Xu LF, Ming J, Yang PP. WUSCHEL-related homeobox genes cooperate with cytokinin to promote bulbil formation in Lilium lancifolium. Plant Physiol. 2022;190(1):387–402. https://doi.org/10.1093/plphys/kiac259 .
doi: 10.1093/plphys/kiac259
Yang ML, Jiao YL. Regulation of axillary meristem initiation by transcription factors and plant hormones. Front Plant Sci. 2016;18(7):183. https://doi.org/10.3389/fpls.2016.00183 .
doi: 10.3389/fpls.2016.00183
Wang Y, Li J. Molecular basis of plant architecture. Annu Rev Plant Biol. 2008;59:253–79. https://doi.org/10.1146/annurev.arplant.59.032607.092902 .
doi: 10.1146/annurev.arplant.59.032607.092902
Gomez-Roldan V, Fermas S, Brewer PB, Puech-Pagès V, Dun EA, Pillot JP, Letisse F, Matusova R, Danoun S, Portais JC, Bouwmeester H, Bécard G, Beveridge CA, Rameau C, Rochange SF. Strigolactone inhibition of shoot branching. Nature. 2008;455:189–94. https://doi.org/10.1038/nature07271 .
doi: 10.1038/nature07271
Janssen BJ, Drummond RS, Snowden KC. Regulation of axillary shoot development. Curr Opin Plant Biol. 2014;17:28–35. https://doi.org/10.1016/j.pbi.2013.11.004 .
doi: 10.1016/j.pbi.2013.11.004
Yamaguchi S. Gibberellin metabolism and its regulation. Annu Rev Plant Biol. 2008;59:225–51. https://doi.org/10.1146/annurev.arplant.59.032607.092804 .
doi: 10.1146/annurev.arplant.59.032607.092804
Hong SY, Sun B, Straub D, Blaakmeer A, Mineri L, Koch J, Brinch-Pedersen H, Holme IB, Burow M, Lyngs Jørgensen HJ, Albà MM, Wenkel S. Heterologous microProtein expression identifies LITTLE NINJA, a dominant regulator of jasmonic acid signaling. Proc. Natl. Acad. Sci. 2020; 117, 26197–26205. https://doi.org/10.1073/pnas.2005198117
Waldie T, Leyser O. Cytokinin targets auxin transport to promote shoot branching. Plant Physiol. 2018;177(2):803–18. https://doi.org/10.1104/pp.17.01691 .
doi: 10.1104/pp.17.01691
Wang J, Yan DW, Yuan TT, Gao X, Lu YT. A gain-of-function mutation in IAA8 alters Arabidopsis floral organ development by change of jasmonic acid level. Plant Mol Biol. 2013;82:71–83. https://doi.org/10.1007/s11103-013-0039-y .
doi: 10.1007/s11103-013-0039-y
Jiang FX, Wei PW, Wu S, Jiang XY, Shi JT, Chen QB. Transcriptional analysis of the pearl buds on the leaves of Ornithogalum arabicum L. Mol Plant Breed. 2017;15(2):519–31. https://doi.org/10.13271/j.mpb.015.000519 .
doi: 10.13271/j.mpb.015.000519
Jiang FX, Huang YX, Zhou P, Sun XL, Song X, Chen QB. Cloning and functional analysis of QtKN1 gene of Ornithogalum arabicum L. Mol Plant Breed. 2018;16(6):1777–83. https://doi.org/10.13271/j.mpb.016.001777 .
doi: 10.13271/j.mpb.016.001777
Jiang FX, Wei KY, Sun XL, Song X, Wen HY, Chen QB. Cloning and characteristic analysis of QtWOX8 gene of Ornithogalum arabicum L. Mol Plant Breed. 2018;16(14):4600–6. https://doi.org/10.13271/j.mpb.016.004600 .
doi: 10.13271/j.mpb.016.004600
Jiang FX, Huang YX, Zhou P, Zhao J, Wen HY, Zhang L. Cloning and functional analysis of OtDCAF8 gene of Ornithogalum arabicum L. Mol Plant Breed. 2019;17(21):7010–6. https://doi.org/10.13271/j.mpb.017.007010 .
doi: 10.13271/j.mpb.017.007010
Rosin FM, Hart JK, Van Onckelen H, Hannapel DJ. Suppression of a vegetative MADS box gene of potato activates axillary meristem development. Plant Physiol. 2003;131(4):1613–22. https://doi.org/10.1104/pp.102.012500 .
doi: 10.1104/pp.102.012500
Jasinski S, Piazza P, Craft J, Hay A, Woolley L, Rieu I, Phillips A, Hedden P, Tsiantis M. KNOX action in Arabidopsis is mediated by coordinate regulation of cytokinin and gibberellin activities. Curr Biol. 2005;15(17):1560–5. https://doi.org/10.1016/j.cub.2005.07.023 .
doi: 10.1016/j.cub.2005.07.023
Abraham-Juárez MJ, Martínez-Hernández A, Leyva-González MA, Herrera-Estrella L, Simpson J. Class I KNOX genes are associated with organogenesis during bulbil formation in Agave tequilana. J Exp Bot. 2010;61(14):4055–67. https://doi.org/10.1093/jxb/erq215 .
doi: 10.1093/jxb/erq215
Delgado Sandoval S, del C, Abraham-Juárez MJ, Simpson J. Agave tequilana MADS genes show novel expression patterns in meristems, developing bulbils and floral organs. Sex Plant Reprod. 2012;25(1):11–26. https://doi.org/10.1007/s00497-011-0176-x .
doi: 10.1007/s00497-011-0176-x
Liu D, Luo R, Chen HL, Xia LS, Gao HX. Transcriptome and expression profile analysis revealed that PEBP gene family members were involved in the regulation of bulblet development of Pinellia cordata. J Guizhou Univ (Natural Sci Edition). 2018;35(2):48–53. https://doi.org/10.15958/j.cnki.gdxbzrb.2018.02.10 .
doi: 10.15958/j.cnki.gdxbzrb.2018.02.10
Yang PP, Xu H, Xu LF, Tang YC, Ming J. Cloning and expression analysis of LlAGO1 in Lilium lancifolium. Acta Horticulturae Sinica. 2018;45:784–94. https://doi.org/10.16420/jissn0513-353x2017-0601 .
doi: 10.16420/jissn0513-353x2017-0601
Bohmert K, Camus I, Bellini C, Bouchez D, Caboche M, Benning C. AGO1 defines a novel locus of Arabidopsis controlling leaf development. EMBO J. 1998;17(1):170–80. https://doi.org/10.1093/emboj/17.1.170 .
doi: 10.1093/emboj/17.1.170
Yang F, Wang Q, Schmitz G, Müller D, Theres K. The bHLH protein ROX acts in concert with RAX1 and LAS to modulate axillary meristem formation in Arabidopsis. Plant J. 2012;71(1):61–70. https://doi.org/10.1111/j.1365-313X.2012.04970.x .
doi: 10.1111/j.1365-313X.2012.04970.x
Dong LL, Yang TY, Gao D, Wang T, Deng XY. PhMYB37 promotes shoot branching in Petunia. Genes. 2022;13(11):2064. https://doi.org/10.3390/genes13112064 .
doi: 10.3390/genes13112064
Schmitz G, Tillmann E, Carriero F, Fiore C, Cellini F, Theres K. The tomato blind gene encodes a MYB transcription factor that controls the formation of lateral meristems. Proc Natl Acad Sci U S A. 2002;99(2):1064–9. https://doi.org/10.1073/pnas.022516199 .
doi: 10.1073/pnas.022516199
Song JJ, Chen YJ, Li X, Ma QQ, Liu QL, Pan YJ, Jiang BB. Cloning and functional verification of CmRAX2 gene associated with chrysanthemum lateral branches development. Genes. 2022;27(5):779. https://doi.org/10.3390/genes13050779 .
doi: 10.3390/genes13050779
Müller D, Schmitz G, Theres K. Blind homologous R2R3 myb genes control the pattern of lateral meristem initiation in Arabidopsis. Plant Cell. 2006;18(3):586–97. https://doi.org/10.1105/tpc.105.038745 .
doi: 10.1105/tpc.105.038745
Guo DS, Zhang JZ, Wang XL, Han X, Wei BY, Wang JQ, Li BX, Yu H, Huang QP, Gu HY, Qu LJ, Qin GJ, The. WRKY transcription factor WRKY71/EXB1 controls shoot branching by transcriptionally regulating RAX genes in Arabidopsis. Plant Cell. 2015,27(11):3112-27. https://doi.org/10.1105/tpc.15.00829
Li G, Serek M, Gehl C. Physiological changes besides the enhancement of pigmentation in Petunia hybrida caused by overexpression of PhAN2, an R2R3-MYB transcription factor. Plant Cell Rep. 2023;42(3):609–27. https://doi.org/10.1007/s00299-023-02983-1 .
doi: 10.1007/s00299-023-02983-1
Yang H, Xue Q, Zhang ZZ, Du JY, Yu DY, Huang F. GmMYB181, a soybean R2R3-MYB protein, increases branch number in transgenic Arabidopsis. Front Plant Sci. 2018;9:1027. https://doi.org/10.3389/fpls.2018.01027 .
doi: 10.3389/fpls.2018.01027
Jia TQ, Zhang KD, Li F, Huang YF, Fan MM, Huang T. The AtMYB2 inhibits the formation of axillary meristem in Arabidopsis by repressing RAX1 gene under environmental stresses. Plant Cell Rep. 2020;39(12):1755–65. https://doi.org/10.1007/s00299-020-02602-3 .
doi: 10.1007/s00299-020-02602-3
Lee DK, Geisler M, Springer PS. LATERAL ORGAN FUSION1 and LATERAL ORGAN FUSION2 function in lateral organ separation and axillary meristem formation in Arabidopsis. Development. 2009;136(14):2423–32. https://doi.org/10.1242/dev.031971 .
doi: 10.1242/dev.031971
Zhou HP, He JX, Zhang YY, Zhao HY, Sun X, Chen X, Liu XR, Zheng Y, Lin HH. RHA2b-mediated MYB30 degradation facilitates MYB75-regulated, sucrose-induced anthocyanin biosynthesis in Arabidopsis seedlings. Plant Commun. 2024;5(3):100744. https://doi.org/10.1016/j.xplc.2023.100744 .
doi: 10.1016/j.xplc.2023.100744
Heim MA, Jakoby M, Werber M, Martin C, Weisshaar B, Bailey PC. The basic helix-loop-helix transcription factor family in plants: a genome-wide study of protein structure and functional diversity. Mol Biol Evol. 2003;20:735–47. https://doi.org/10.1093/molbev/msg088 .
doi: 10.1093/molbev/msg088
Kurbidaeva A, Ezhova T, Novokreshchenova M. Arabidopsis thaliana ICE2 gene: phylogeny, structural evolution and functional diversification from ICE1. Plant Sci. 2014;229:10–22. https://doi.org/10.1016/j.plantsci.2014.08.011 .
doi: 10.1016/j.plantsci.2014.08.011
Lin YJ, Li MJ, Hsing HC, Chen TK, Yang TT, Ko SS. activatorivator 1, encoding a bHLH, mediates axillary bud development and spike initiation in PhalaenAphroditerodite. Int J Mol Sci. 2019;30(21):5406. https://doi.org/10.3390/ijms20215406 .
doi: 10.3390/ijms20215406
Holalu SV, Reddy SK, Blackman BK, Finlayson SA. Phytochrome interacting factors 4 and 5 regulate axillary branching via bud abscisic acid and stem auxin signalling. Plant Cell Environ. 2020;43(9):2224–38. https://doi.org/10.1111/pce.13824 .
doi: 10.1111/pce.13824
Basile A, Fambrini M, Tani C, Shukla V, Licausi F, Pugliesi C. The Ha-ROXL gene is required for initiation of axillary and floral meristems in sunflower. Genesis. 2019;57(9):e23307. https://doi.org/10.1002/dvg.23307 .
doi: 10.1002/dvg.23307
Yao H, Skirpan A, Wardell B, Matthes MS, Best NB, McCubbin T, Durbak A, Smith T, Malcomber S, McSteen P. The barren stalk2 gene is required for axillary meristem development in maize. Mol Plant. 2019;12(3):374–89. https://doi.org/10.1016/j.molp.2018.12.024 .
doi: 10.1016/j.molp.2018.12.024
Mach J. Rice axillary meristem formation requires directional movement of LAX PANICLE1 protein. Plant Cell. 2009;21(4):1027. https://doi.org/10.1105/tpc.109.210410 .
doi: 10.1105/tpc.109.210410
Lv YP, Zhang XF, Hu YJ, Liu S, Yin YB, Wang XX. BOS1 is a basic helix-loop-helix transcription factor involved in regulating panicle development in rice. Front Plant Sci. 2023;14:1162828. https://doi.org/10.3389/fpls.2023.1162828 .
doi: 10.3389/fpls.2023.1162828
Tabuchi H, Zhang Y, Hattori S, Omae M, Shimizu-Sato S, Oikawa T, Qian Q, Nishimura M, Kitano H, Xie H, Fang X, Yoshida H, Kyozuka J, Chen F, Sato Y. LAX PANICLE2 of rice encodes a novel nuclear protein and regulates the formation of axillary meristems. Plant Cell. 2011;23(9):3276–87. https://doi.org/10.1105/tpc.111.088765 .
doi: 10.1105/tpc.111.088765
He KR, Du JC, Han X, Li HQ, Kui MY, Zhang JP, Huang ZC, Fu QT, Jiang YJ, Hu YR. PHOSPHATE STARVATION RESPONSE1 (PHR1) interacts with JASMONATE ZIM-DOMAIN (JAZ) and MYC2 to modulate phosphate deficiency-induced jasmonate signaling in Arabidopsis. Plant Cell. 2023;35(6):2132–56. https://doi.org/10.1093/plcell/koad057 .
doi: 10.1093/plcell/koad057
Hu YR, Jiang LQ, Wang F, Yu DQ. Jasmonate regulates the INDUCER OF CBF EXPRESSION-C-REPEAT BINDING FACTOR/DRE BINDING FACTOR1 cascade and freezing tolerance in Arabidopsis. Plant Cell. 2013;25(8):2907–24. https://doi.org/10.1105/tpc.113.112631 .
doi: 10.1105/tpc.113.112631
Grunewald W, Vanholme B, Pauwels L, Plovie E, Inzé D, Gheysen G, Goossens A. Expression of the Arabidopsis jasmonate signalling repressor JAZ1/TIFY10A is stimulated by auxin. EMBO Rep. 2009;10(8):923–8. https://doi.org/10.1038/embor.2009.103 .
doi: 10.1038/embor.2009.103
Floková K, Tarkowská D, Miersch O, Strnad M, Wasternack C, Novák O. UHPLC-MS/MS based target profiling of stress-induced phytohormones. Phytochemistry. 2014;105:147–57. https://doi.org/10.1016/j.phytochem.2014.05.015 .
doi: 10.1016/j.phytochem.2014.05.015
Li Y, Zhou CX, Yan XJ, Zhang JR, Xu JL. Simultaneous analysis of ten phytohormones in Sargassum Horneri by high-performance liquid chromatography with electrospray ionization tandem mass spectrometry. J Sep Sci. 2016;39(10):1804–13. https://doi.org/10.1002/jssc.201501239 .
doi: 10.1002/jssc.201501239
Pan XQ, Welti R, Wang XM. Quantitative analysis of major plant hormones in crude plant extracts by high-performance liquid chromatography-mass spectrometry. Nat Protoc. 2010;5(6):986–92. https://doi.org/10.1038/nprot.2010.37 .
doi: 10.1038/nprot.2010.37
Cui KY, Lin YY, Zhou XY, Li SC, Liu H, Zeng F, Zhu F, Ouyang GF, Zeng ZX. Comparison of sample pretreatment methods for the determination of multiple phytohormones in plant samples by liquid chromatography-electrospray ionization-tandem mass spectrometry. Microchem J. 2015;121:25–31. https://doi.org/10.1016/j.microc.2015.02.004 .
doi: 10.1016/j.microc.2015.02.004
Šimura J, Antoniadi I, Široká J, Tarkowská D, Strnad M, Ljung K, Novák O. Plant hormonomics: multiple phytohormone profiling by targeted metabolomics. Plant Physiol. 2018;177(2):476–89. https://doi.org/10.1104/pp.18.00293 .
doi: 10.1104/pp.18.00293
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads EMBnet. Journal. 2011;10–2. https://doi.org/10.14806/EJ.17.1.200 . ISSN 2226–6089.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644–52. https://doi.org/10.1038/nbt.1883 .
doi: 10.1038/nbt.1883
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60. https://doi.org/10.1038/nmeth.3176 .
doi: 10.1038/nmeth.3176
Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417–9. https://doi.org/10.1038/nmeth.4197 .
doi: 10.1038/nmeth.4197
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. https://doi.org/10.1038/nmeth.1226 .
doi: 10.1038/nmeth.1226
Robinson MD, McCarthy DJ, Smyth GK. Bioinformatics. 2010;26(1):139–40. https://doi.org/10.1093/bioinformatics/btp616 . edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.
Chen MM, Zhang RJ, Cha Q, Yang LY, Li Xin, Yin LQ, Zhang YC. Induction and development of lily somatic embryos and screening of real-time quantitative PCR internal reference genes in different tissues. Mol Plant Breed. 2018;16(15):4982–90. https://doi.org/10.13271/j.mpb.016.004982 .
doi: 10.13271/j.mpb.016.004982