A photorespiratory bypass increases plant growth and seed yield in biofuel crop Camelina sativa
Biotechnology for Biofuels - 2015
Tóm tắt
Camelina sativa is an oilseed crop with great potential for biofuel production on marginal land. The seed oil from camelina has been converted to jet fuel and improved fuel efficiency in commercial and military test flights. Hydrogenation-derived renewable diesel from camelina is environmentally superior to that from canola due to lower agricultural inputs, and the seed meal is FDA approved for animal consumption. However, relatively low yield makes its farming less profitable. Our study is aimed at increasing camelina seed yield by reducing carbon loss from photorespiration via a photorespiratory bypass. Genes encoding three enzymes of the Escherichia coli glycolate catabolic pathway were introduced: glycolate dehydrogenase (GDH), glyoxylate carboxyligase (GCL) and tartronic semialdehyde reductase (TSR). These enzymes compete for the photorespiratory substrate, glycolate, convert it to glycerate within the chloroplasts, and reduce photorespiration. As a by-product of the reaction, CO2 is released in the chloroplast, which increases photosynthesis. Camelina plants were transformed with either partial bypass (GDH), or full bypass (GDH, GCL and TSR) genes. Transgenic plants were evaluated for physiological and metabolic traits. Expressing the photorespiratory bypass genes in camelina reduced photorespiration and increased photosynthesis in both partial and full bypass expressing lines. Expression of partial bypass increased seed yield by 50–57 %, while expression of full bypass increased seed yield by 57–73 %, with no loss in seed quality. The transgenic plants also showed increased vegetative biomass and faster development; they flowered, set seed and reached seed maturity about 1 week earlier than WT. At the transcriptional level, transgenic plants showed differential expression in categories such as respiration, amino acid biosynthesis and fatty acid metabolism. The increased growth of the bypass transgenics compared to WT was only observed in ambient or low CO2 conditions, but not in elevated CO2 conditions. The photorespiratory bypass is an effective approach to increase photosynthetic productivity in camelina. By reducing photorespiratory losses and increasing photosynthetic CO2 fixation rates, transgenic plants show dramatic increases in seed yield. Because photorespiration causes losses in productivity of most C3 plants, the bypass approach may have significant impact on increasing agricultural productivity for C3 crops.
Từ khóa
Tài liệu tham khảo
Budin JT, Breene WM, Putnam DH. Some compositional properties of camelina (Camelina-Sativa L. Crantz) seeds and oils. J Am Oil Chem Soc. 1995;72(3):309–15.
Zubr J. Qualitative variation of Camelina sativa seed from different locations. Ind Crop Prod. 2003;17(3):161–9. doi:10.1016/S0926-6690(02)00091-2.
Vollmann J, Moritz T, Kargl C, Baumgartner S, Wagentristl H. Agronomic evaluation of camelina genotypes selected for seed quality characteristics. Ind Crop Prod. 2007;26(3):270–7. doi:10.1016/j.indcrop.2007.03.017.
Moser BR. Biodiesel from alternative oilseed feedstocks: camelina and field pennycress. Biofuels. 2012;3(2):193–209. doi:10.4155/BFS.12.6.
Bernardo A, Howard-Hildige R, O’Connell A, Nichol R, Ryan J, Rice B, et al. Camelina oil as a fuel for diesel transport engines. Ind Crop Prod. 2003;17(3):191–7. doi:10.1016/S0926-6690(02)00098-5.
Frohlich A, Rice B. Evaluation of Camelina sativa oil as a feedstock for biodiesel production. Ind Crop Prod. 2005;21(1):25–31. doi:10.1016/j.indcrop.2003.12.004.
Shonnard DR, Williams L, Kalnes TN. Camelina-derived jet fuel and diesel: sustainable advanced biofuels. Environ Prog Sustain. 2010;29(3):382–92. doi:10.1002/Ep.10461.
Natelson RH, Wang W, William RL, Zering KD. Technoeconomic analysis of jet fuel production from hydrolysis, decarboxylation, and reforming of camelina oil. Biomass Bioenergy. 2015;75:23–34.
Helsel ZR. Camelina for biofuel production. Farm Energy. 2013. http://www.extension.org/pages/69879/camelina-for-biofuel-production#.VipEkivLJjZ.
Foulke T, Geiger, M, Hess, B. Research summary: economics of dryland camelina biofuels production farm energy. 2014. http://www.extension.org/pages/70514/research-summary:-economics-of-dryland-camelina-biofuels-production.
Bota J, Flexas J, Keys AJ, Loveland J, Parry MAJ, Medrano H. CO(2)/O(2) specificity factor of ribulose-1,5-bisphosphate carboxylase/oxygenase in grapevines (Vitis vinifera L.): first in vitro determination and comparison to in vivo estimations. Vitis. 2002;41(4):163–8.
Jordan DB, Ogren WL. Species variation in kinetic-properties of ribulose 1,5-bisphosphate carboxylase oxygenase. Arch Biochem Biophys. 1983;227(2):425–33. doi:10.1016/0003-9861(83)90472-1.
Sharkey TD, Berry JA, Sage RF. Regulation of photosynthetic electron-transport in Phaseolus vulgaris L., as determined by room-temperature chlorophyll a fluorescence. Planta. 1988;176(3):415–24. doi:10.1007/BF00395423.
Anderson LE. Chloroplast and cytoplasmic enzymes. II. Pea leaf triose phosphate isomerases. Biochim Biophys Acta. 1971;235(1):237–44.
Peterhansel C, Horst I, Niessen M, Blume C, Kebeish R, Kurkcuoglu S, et al. Photorespiration. Arabidopsis Book Am Soc Plant Biol. 2010;8:e0130. doi:10.1199/tab.0130.
Kozaki A, Takeba G. Photorespiration protects C3 plants from photooxidation. Nature. 1996;384(6609):557–60. doi:10.1038/384557a0.
Fahnenstich H, Flugge UI, Maurino VG. Arabidopsis thaliana overexpressing glycolate oxidase in chloroplasts: H(2)O(2)-induced changes in primary metabolic pathways. Plant Signal Behav. 2008;3(12):1122–5.
Xu H, Zhang J, Zeng J, Jiang L, Liu E, Peng C, et al. Inducible antisense suppression of glycolate oxidase reveals its strong regulation over photosynthesis in rice. J Exp Bot. 2009;60(6):1799–809. doi:10.1093/jxb/erp056.
Bloom AJ, Burger M, Rubio-Asensio JS, Cousins AB. Carbon dioxide enrichment inhibits nitrate assimilation in wheat and Arabidopsis. Science. 2010;328(5980):899–903. doi:10.1126/science.1186440.
Rachmilevitch S, Cousins AB, Bloom AJ. Nitrate assimilation in plant shoots depends on photorespiration. Proc Natl Acad Sci USA. 2004;101(31):11506–10. doi:10.1073/pnas.0404388101.
Kebeish R, Niessen M, Thiruveedhi K, Bari R, Hirsch HJ, Rosenkranz R, et al. Chloroplastic photorespiratory bypass increases photosynthesis and biomass production in Arabidopsis thaliana. Nat Biotechnol. 2007;25(5):593–9.
Maier A, Fahnenstich H, vonCaemmerer S, Engqvist MK, Weber AP, Flugge UI, et al. Transgenic introduction of a glycolate oxidative cycle into A. thaliana chloroplasts leads to growth improvement. Front Plant Sci. 2012;3:38. doi:10.3389/fpls.2012.00038.
Nolke G, Houdelet M, Kreuzaler F, Peterhansel C, Schillberg S. The expression of a recombinant glycolate dehydrogenase polyprotein in potato (Solanum tuberosum) plastids strongly enhances photosynthesis and tuber yield. Plant Biotechnol J. 2014;. doi:10.1111/pbi.12178.
Xin CP, Tholen D, Devloo V, Zhu XG. The benefits of photorespiratory bypasses: How can they work? Plant Physiol. 2015;167(2):574–85. doi:10.1104/pp.114.248013.
Kearney PC, Tolbert NE. Appearance of glycolate and related products of photosynthesis outside of chloroplasts. Arch Biochem Biophys. 1962;98:164–71.
Kisaki T, Tolbert NE. Glycolate and glyoxylate metabolism by isolated peroxisomes or chloroplasts. Plant Physiol. 1969;44(2):242–50.
Lawyer AL, Cornwell KL, Gee SL, Bassham JA. Glyoxylate and glutamate effects on photosynthetic carbon metabolism in isolated chloroplasts and mesophyll cells of spinach. Plant Physiol. 1983;72(2):420–5.
Mulligan RM, Wilson B, Tolbert NE. Effects of glyoxylate on photosynthesis by intact chloroplasts. Plant Physiol. 1983;72(2):415–9.
Oliver DJ. The effect of glyoxylate on photosynthesis and photorespiration by isolated soybean mesophyll cells. Plant Physiol. 1980;65(5):888–92.
Oliver DJ, Zelitch I. Increasing photosynthesis by inhibiting photorespiration with glyoxylate. Science. 1977;196(4297):1450–1.
Lord JM. Glycolate oxidoreductase in Escherichia coli. Biochim Biophys Acta. 1972;267(2):227–37. doi:10.1016/0005-2728(72)90111-9.
Gotto AM, Kornberg HL. Metabolism of C2 compounds in micro-organisms. 7. Preparation and properties of crystalline tartronic semialdehyde reductase. Biochem J. 1961;81(2):273.
Yu C, Liang Z, Huang AHC. Glyoxylate transamination in intact leaf peroxisomes. Plant Physiol. 1984;75(1):7–12. doi:10.1104/Pp.75.1.7.
Berry JA, Osmond CB, Lorimer GH. Fixation of O(2) during photorespiration: kinetic and steady-state studies of the photorespiratory carbon oxidation cycle with intact leaves and isolated chloroplasts of C(3) plants. Plant Physiol. 1978;62(6):954–67.
Novitskaya L, Trevanion SJ, Driscoll S, Foyer CH, Noctor G. How does photorespiration modulate leaf amino acid contents? A dual approach through modelling and metabolite analysis. Plant Cell Environ. 2002;25(7):821–35. doi:10.1046/j.1365-3040.2002.00866.x.
Farquhar GD, von Caemmerer S, Berry JA. A biochemical model of photosynthetic CO2 assimilation in leaves of C3 species. Planta. 1980;149(1):78–90. doi:10.1007/BF00386231.
Sharkey TD, Bernacchi CJ, Farquhar GD, Singsaas EL. Fitting photosynthetic carbon dioxide response curves for C(3) leaves. Plant Cell Environ. 2007;30(9):1035–40. doi:10.1111/j.1365-3040.2007.01710.x.
Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, et al. TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003;34(2):374–8.
Laisk A, Loreto F. Determining photosynthetic parameters from leaf CO2 exchange and chlorophyll fluorescence (ribulose-1,5-bisphosphate carboxylase/oxygenase specificity factor, dark respiration in the light, excitation distribution between photosystems, alternative electron transport rate, and mesophyll diffusion resistance. Plant Physiol. 1996;110(3):903–12.
Walker BJ, Ort DR. Improved method for measuring the apparent CO photocompensation point resolves the impact of multiple internal conductances to CO to net gas exchange. Plant Cell Environ. 2015;. doi:10.1111/pce.12562.
Busch FA. Current methods for estimating the rate of photorespiration in leaves. Plant Biol. 2013;15(4):648–55. doi:10.1111/j.1438-8677.2012.00694.x.
Szecowka M, Heise R, Tohge T, Nunes-Nesi A, Vosloh D, Huege J, et al. Metabolic fluxes in an illuminated Arabidopsis rosette. Plant Cell. 2013;25(2):694–714. doi:10.1105/tpc.112.106989.
Ma F, Jazmin LJ, Young JD, Allen DK. Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation. Proc Natl Acad Sci USA. 2014;111(47):16967–72. doi:10.1073/pnas.1319485111.
Hashimoto M, Negi J, Young J, Israelsson M, Schroeder JI, Iba K. Arabidopsis HT1 kinase controls stomatal movements in response to CO2. Nat Cell Biol. 2006;8(4):391–7. doi:10.1038/ncb1387.
Blume C, Behrens C, Eubel H, Braun HP, Peterhansel C. A possible role for the chloroplast pyruvate dehydrogenase complex in plant glycolate and glyoxylate metabolism. Phytochemistry. 2013;95:168–76. doi:10.1016/j.phytochem.2013.07.009.
Reid EE, Thompson P, Lyttle CR, Dennis DT. Pyruvate dehydrogenase complex from higher plant mitochondria and proplastids. Plant Physiol. 1977;59(5):842–8.
Davey PA, Hunt S, Hymus GJ, DeLucia EH, Drake BG, Karnosky DF, et al. Respiratory oxygen uptake is not decreased by an instantaneous elevation of [CO2], but is increased with long-term growth in the field at elevated [CO2]. Plant Physiol. 2004;134(1):520–7. doi:10.1104/pp.103.030569.
Leakey AD, Xu F, Gillespie KM, McGrath JM, Ainsworth EA, Ort DR. Genomic basis for stimulated respiration by plants growing under elevated carbon dioxide. Proc Natl Acad Sci USA. 2009;106(9):3597–602. doi:10.1073/pnas.0810955106.
Bouma TJ, Devisser R, Vanleeuwen PH, Dekock MJ, Lambers H. The respiratory energy-requirements involved in nocturnal carbohydrate export from starch-storing mature source leaves and their contribution to leaf dark respiration. J Exp Bot. 1995;46(290):1185–94. doi:10.1093/jxb/46.9.1185.
Morcuende R, Perez P, MartinezCarrasco R. Short-term feedback inhibition of photosynthesis in wheat leaves supplied with sucrose and glycerol at two temperatures. Photosynthetica. 1997;33(2):179–88. doi:10.1023/A:1022199927264.
Franck N, Vaast P, Genard M, Dauzat J. Soluble sugars mediate sink feedback down-regulation of leaf photosynthesis in field-grown Coffea arabica. Tree Physiol. 2006;26(4):517–25.
Quereix A, Dewar RC, Gaudillere JP, Dayau S, Valancogne C. Sink feedback regulation of photosynthesis in vines: measurements and a model. J Exp Bot. 2001;52(365):2313–22.
Paul MJ, Foyer CH. Sink regulation of photosynthesis. J Exp Bot. 2001;52(360):1383–400.
Stitt M. Metabolic regulation of photosynthesis. In: Baker NR, editor. Photosynthesis and the environment. vol 5. Netherlands: Springer; 1996. p. 151–90. doi:10.1007/0-306-48135-9_6.
Pritchard SG, Rogers HH, Prior SA, Peterson CM. Elevated CO2 and plant structure: a review. Glob Change Biol. 1999;5(7):807–37. doi:10.1046/j.1365-2486.1999.00268.x.
Watling JR, Press MC, Quick WP. Elevated CO(2) induces biochemical and ultrastructural changes in leaves of the C(4) cereal sorghum. Plant Physiol. 2000;123(3):1143–52.
Tosens T, Niinemets U, Vislap V, Eichelmann H, Castro Diez P. Developmental changes in mesophyll diffusion conductance and photosynthetic capacity under different light and water availabilities in Populus tremula: how structure constrains function. Plant Cell Environ. 2012;35(5):839–56. doi:10.1111/j.1365-3040.2011.02457.x.
Li Y, Ren BB, Ding L, Shen QR, Peng SB, Guo SW. Does chloroplast size influence photosynthetic nitrogen use efficiency? PloS ONE. 2013;8(4):e62036. doi:10.1371/journal.pone.0062036.
Evans JR, Voncaemmerer S, Setchell BA, Hudson GS. The relationship between CO2 transfer conductance and leaf anatomy in transgenic tobacco with a reduced content of Rubisco. Aust J Plant Physiol. 1994;21(4):475–95.
Muller O, Oguchi R, Hirose T, Werger MJ, Hikosaka K. The leaf anatomy of a broad-leaved evergreen allows an increase in leaf nitrogen content in winter. Physiol Plant. 2009;136(3):299–309. doi:10.1111/j.1399-3054.2009.01224.x.
Oguchi R, Hikosaka K, Hirose T. Does the photosynthetic light-acclimation need change in leaf anatomy? Plant Cell Environ. 2003;26(4):505–12. doi:10.1046/j.1365-3040.2003.00981.x.
Chen A, Lichstein JW, Osnas JL, Pacala SW. Species-independent down-regulation of leaf photosynthesis and respiration in response to shading: evidence from six temperate tree species. PLoS One. 2014;9(4):e91798. doi:10.1371/journal.pone.0091798.
Long SP, Ainsworth EA, Leakey AD, Nosberger J, Ort DR. Food for thought: lower-than-expected crop yield stimulation with rising CO2 concentrations. Science. 2006;312(5782):1918–21. doi:10.1126/science.1114722.
Long SP, Marshall-Colon A, Zhu XG. Meeting the global food demand of the future by engineering crop photosynthesis and yield potential. Cell. 2015;161(1):56–66. doi:10.1016/j.cell.2015.03.019.
Cassman KG, Dobermann A, Walters DT. Agroecosystems, nitrogen-use efficiency, and nitrogen management. Ambio. 2002;31(2):132–40.
Roberts TL, Ross WJ, Norman RJ, Slaton NA, Wilson CE. Predicting nitrogen fertilizer needs for rice in arkansas using alkaline hydrolyzable- nitrogen. Soil Sci Soc Am J. 2011;75(3):1161–71. doi:10.2136/sssaj2010.0145.
Crous KY, Walters MB, Ellsworth DS. Elevated CO2 concentration affects leaf photosynthesis-nitrogen relationships in Pinus taeda over nine years in FACE. Tree Physiol. 2008;28(4):607–14.
Walker AP, Beckerman AP, Gu LH, Kattge J, Cernusak LA, Domingues TF, et al. The relationship of leaf photosynthetic traits—V-cmax and J(max)—to leaf nitrogen, leaf phosphorus, and specific leaf area: a meta-analysis and modeling study. Ecol Evol. 2014;4(16):3218–35. doi:10.1002/ece3.1173.
Liang C, Liu X, Yiu SM, Lim BL. De novo assembly and characterization of Camelina sativa transcriptome by paired-end sequencing. BMC Genom. 2013;14(1):146. doi:10.1186/1471-2164-14-146.
Ow DW, Jacobs JD, Howell SH. Functional regions of the cauliflower mosaic virus 35S RNA promoter determined by use of the firefly luciferase gene as a reporter of promoter activity. Proc Natl Acad Sci USA. 1987;84(14):4870–4.
Malik K, Wu K, Li XQ, Martin-Heller T, Hu M, Foster E, et al. A constitutive gene expression system derived from the tCUP cryptic promoter elements. Theor Appl Genet. 2002;105(4):505–14. doi:10.1007/s00122-002-0926-0.
An YQC, Meagher RB. Strong Expression and Conserved Regulation of ACT2 in Arabidopsis thaliana and Physcomitrella patens. Plant Mol Biol Rep. 2010;28(3):481–90. doi:10.1007/s11105-009-0171-7.
Lee DW, Kim JK, Lee S, Choi S, Kim S, Hwang I. Arabidopsis nuclear-encoded plastid transit peptides contain multiple sequence subgroups with distinctive chloroplast-targeting sequence motifs. Plant Cell. 2008;20(6):1603–22.
Nelson BK, Cai X, Nebenfuhr A. A multicolored set of in vivo organelle markers for co-localization studies in Arabidopsis and other plants. Plant J. 2007;51(6):1126–36. doi:10.1111/j.1365-313X.2007.03212.x.
Dalal J, Yalamanchili R, La Hovary C, Ji M, Rodriguez-Welsh M, Aslett D, et al. A novel gateway-compatible binary vector series (PC-GW) for flexible cloning of multiple genes for genetic transformation of plants. Plasmid. 2015;81:55–62. doi:10.1016/j.plasmid.2015.06.003.
Earley KW, Haag JR, Pontes O, Opper K, Juehne T, Song K, et al. Gateway-compatible vectors for plant functional genomics and proteomics. Plant J. 2006;45(4):616–29. doi:10.1111/j.1365-313X.2005.02617.x.
Cao J, Duan XL, Mcelroy D, Wu R. Regeneration of herbicide resistant transgenic rice plants following microprojectile-mediated transformation of suspension-culture cells. Plant Cell Rep. 1992;11(11):586–91.
Christou P, Ford TL, Kofron M. Production of transgenic rice (Oryza sativa L.) plants from agronomically important indica and japonica varieties via electric-discharge particle-acceleration of exogenous DNA into immature zygotic embryos. Nat Biotechnol. 1991;9(10):957–62. doi:10.1038/Nbt1091-957.
Lu C, Kang J. Generation of transgenic plants of a potential oilseed crop Camelina sativa by Agrobacterium-mediated transformation. Plant Cell Rep. 2008;27(2):273–8. doi:10.1007/s00299-007-0454-0.
Li L, Zhou Y, Cheng X, Sun J, Marita JM, Ralph J, et al. Combinatorial modification of multiple lignin traits in trees through multigene cotransformation. Proc Natl Acad Sci USA. 2003;100(8):4939–44. doi:10.1073/pnas.0831166100.
Dalal J, Land E, Vasani N, He L, Smith C, Rodriguez-Welsh M et al. Methods for RNA profiling of gravi-responding plant tissues. In: Blancaflor EB, editor. Plant gravitropism: methods and protocols. vol 1309. New York: Springer; 2015. p. 97–117. doi:10.1007/978-1-4939-2697-8_9.
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.
Joly D, Carpentier R. Rapid isolation of intact chloroplasts from spinach leaves. Methods Mol Biol. 2011;684:321–5. doi:10.1007/978-1-60761-925-3_24.
Li Y, Schellhorn HE. Rapid kinetic microassay for catalase activity. J Biomol Tech JBT. 2007;18(4):185–7.
Liu Y, Koh CM, Sun L, Ji L. Tartronate semialdehyde reductase defines a novel rate-limiting step in assimilation and bioconversion of glycerol in Ustilago maydis. PLoS One. 2011;6(1):e16438. doi:10.1371/journal.pone.0016438.
Wuyts N, Palauqui JC, Conejero G, Verdeil JL, Granier C, Massonnet C. High-contrast three-dimensional imaging of the Arabidopsis leaf enables the analysis of cell dimensions in the epidermis and mesophyll. Plant Methods. 2010;6:17. doi:10.1186/1746-4811-6-17.
Porra RJ, Thompson WA, Kriedemann PE. Determination of accurate extinction coefficients and simultaneous-equations for assaying chlorophyll-a and chlorophyll-B extracted with 4 different solvents—verification of the concentration of chlorophyll standards by atomic-absorption spectroscopy. Biochim Biophys Acta. 1989;975(3):384–94. doi:10.1016/S0005-2728(89)80347-0.
Gromova M, Roby C. Toward Arabidopsis thaliana hydrophilic metabolome: assessment of extraction methods and quantitative 1H NMR. Physiol Plant. 2010;140(2):111–27. doi:10.1111/j.1399-3054.2010.01387.x.
Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Physiol. 1959;37(8):911–7.
Winter H, Robinson DG, Heldt HW. subcellular volumes and metabolite concentrations in barley leaves. Planta. 1993;191(2):180–90.
Li Y, Beisson F, Pollard M, Ohlrogge J. Oil content of Arabidopsis seeds: the influence of seed anatomy, light and plant-to-plant variation. Phytochemistry. 2006;67(9):904–15. doi:10.1016/j.phytochem.2006.02.015.
Christie WW, Brechany EY, Jie MSFLK, Bakare O. Mass-spectral characterization of picolinyl and methyl-ester derivatives of isomeric thia fatty-acids. Biol Mass Spectrom. 1991;20(10):629–35. doi:10.1002/bms.1200201007.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8(8):1494–512. doi:10.1038/nprot.2013.084.
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005;21(18):3674–6. doi:10.1093/bioinformatics/bti610.
Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25(14):1754–60. doi:10.1093/bioinformatics/btp324.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40. doi:10.1093/bioinformatics/btp616.