Background The R2R3MYB proteins comprise one of the largest families of transcription factors in plants. Although genome-wide analysis of this family has been carried out in some species, little is known about R2R3MYB genes in cucumber (Cucumis sativus L.). Principal Findings This study has identified 55 R2R3MYB genes in the latest cucumber genome and the CsR2R3MYB family contained the smallest number of identified genes compared to other species that have been studied due to the absence of recent gene duplication events. These results were also supported by genome distribution and gene duplication analysis. Phylogenetic analysis showed that they could be classified into 11 subgroups. The evolutionary relationships and the intron - exon organizations that showed similarities with Arabidopsis, Vitis and Glycine R2R3MYB proteins were also analyzed and suggested strong gene conservation but also the expansions of particular functional genes during the evolution of the plant species. In addition, we found that 8 out of 55 (~14.54%) cucumber R2R3MYB genes underwent alternative splicing events, producing a variety of transcripts from a single gene, which illustrated the extremely high complexity of transcriptome regulation. Tissue-specific expression profiles showed that 50 cucumber R2R3MYB genes were expressed in at least one of the tissues and the other 5 genes showed very low expression in all tissues tested, which suggested that cucumber R2R3MYB genes took part in many cellular processes. The transcript abundance level analysis during abiotic conditions (NaCl, ABA and low temperature treatments) identified a group of R2R3MYB genes that responded to one or more treatments. Conclusions This study has produced a comparative genomics analysis of the cucumber R2R3MYB gene family and has provided the first steps towards the selection of CsR2R3MYB genes for cloning and functional dissection that can be used in further studies to uncover their roles in cucumber growth and development.
References
[1]
Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, et al. (2010) MYB transcription factors in Arabidopsis. Trends Plant Sci 15: 573–581.
[2]
Rosinski JA, Atchley WR (1998) Molecular evolution of the Myb family of transcription factors: evidence for polyphyletic origin. J Mol Evol 46: 74–83.
[3]
Jin H, Martin C (1999) Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol 41: 577–585.
[4]
Zhang LC, Zhao GY, Jia JZ, Liu x, Kong XY (2012) Molecular characterization of 60 isolated wheat MYB genes and analysis of their expression during abiotic stress. J Exp Bot 63: 203–214.
[5]
Stracke R, Werber M, Weisshaar B (2001) The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol 4: 447–456.
[6]
Ogata K, Morikawa S, Nakamura H, Sekikawa A, Inoue T, et al. (1994) Solution structure of a specific DNA complex of the Myb DNA-binding domain with cooperative recognition helices. Cell 79: 639–648.
[7]
Matus JT, Aquea F, Arce-Johnson P (2008) Analysis of the grape MYB R2R3 subfamily reveals expanded wine quality-related clades and conserved gene structure organization across Vitis and Arabidopsis genomes. BMC Plant Biol 8: 83.
[8]
Dias AP, Braun EL, McMullen MD, Grotewold E (2003) Recently Duplicated Maize R2R3 Myb Genes Provide Evidence for Distinct Mechanisms of Evolutionary Divergence after Duplication. Plant Physiol 131: 610–620.
[9]
Chen YH, Yang XY, He K, Liu MH, Li JG, et al. (2006) The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol 60: 107–124.
[10]
Wilkins O, Nahal H, Foong J, Provart NJ, Campbell MM (2009) Expansion and diversification of the Populus R2R3-MYB family of transcription factors. Plant Physiol 149: 981–993.
[11]
Nesi N, Jond C, Debeaujon I, Caboche M, Lepiniec L (2001) The Arabidopsis TT2 gene encodes an R2R3 MYB domain protein that acts as a key determinant for proanthocyanidin accumulation in developing seed. Plant Cell 13: 2099–2114.
[12]
Baudry A, Heim MA, Dubreucq B, Caboche M, Weisshaar B, et al. (2004) TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant J 39: 366–80.
[13]
Gonzalez A, Zhao M, Leavitt JM, Lloyd AM (2008) Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J 53: 814–827.
[14]
Gonzalez A, Mendenhall J, Huo Y, Lloyd A (2009) TTG1 complex MYBs, MYB5 and TT2, control outer seed coat differentiation. Dev. Biol. 325: 412–421.
[15]
Li SF, Milliken ON, Pham H, Seyit R, Napoli R, et al. (2009) The Arabidopsis MYB5 transcription factor regulates mucilage synthesis, seed coat development, and trichome morphogenesis. Plant Cell 21: 72–89.
[16]
Oppenheimer DG, Herman PL, Sivakumaran S, Esch J, Marks MD (1991) A myb gene required for leaf trichome differentiation in Arabidopsis is expressed in stipules. Cell 67: 483–493.
[17]
Glover BJ, Perez-Rodriguez M, Martin C (1998) Development of several epidermal cell types can be specified by the same MYB-related plant transcription factor. Development 125: 3497–3508.
[18]
Payne CT, Clement J, Arnold D, Lioyd A (1999) Heterologous MYB genes distinct from GL1 enhance trichome production when overexpressed in Nicotiana tabacum. Development 126: 671–682.
[19]
Lee MM (2001) Schiefelbein J (2001) Developmentally distinct MYB genes encode functionally equivalent proteins in Arabidopsis. Development 128: 1539–1546.
[20]
Ito M, Araki S, Matsunaga S, Itoh T, Nishihama R, et al. (2001) G2/Mphase–specific transcription during the plant cell cycle is mediated by c-Myb–Like transcription factors. Plant Cell 13: 1891–1905.
[21]
Araki S, Ito M, Soyano T, Nishihama R, Machida Y (2004) Mitotic cyclins stimulate the activity of c-Myb-like factors for transactivation of G2/M phase-specific genes in tobacco. J Biol Chem 279: 32979–88.
[22]
Brownfield L. Hafidh S, Borg M, Sidorova A, Mori T, et al. (2009) A plant germline-specific integrator of sperm specification and cell cycle progression. PLoS Genet 5: e1000430.
[23]
Millar AA, Gubler F (2005) The Arabidopsis GAMYB-like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell 17: 705–721.
[24]
Müller D, Schmitz G, Theres K (2006) Blind homologous R2R3 Myb genes control the pattern of lateral meristem initiation in Arabidopsis. Plant Cell 18: 586–597.
[25]
Keller T, Abbott J, Moritz T, Doerner P (2006) Arabidopsis REGULATOR OF AXILLARY MERISTEMS1 controls a leaf axil stem cell niche and modulates vegetative development. Plant Cell 18: 598–611.
[26]
Urao T, Yamaguchi-Shinozaki K, Urao S, Shinozaki K (1993) An Arabidopsis myb homolog is induced by dehydration stress and its gene product binds to the conserved MYB recognition sequence. Plant Cell 5: 1529–1539.
[27]
Daniel X, Lacomme C, Morel JB, Roby D (1999) A novel myb oncogene homologue in Arabidopsis thaliana related to hypersensitive cell death. Plant J 20: 57–66.
[28]
Zhu J, Verslues PE, Zheng X, Lee BH, Zhan X, et al. (2005) HOS10 encodes an R2R3-type MYB transcription factor essential for cold acclimation in plants. Proc Natl Acad Sci USA 102: 9966–71.
[29]
Abe H, Urao T, Ito T, Seki M, Shinozaki K, et al. (2003) Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 15: 63–78.
[30]
Mengiste T, Chen X, Salmeron J, Dietrich R (2003) The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell 15: 2551–65.
[31]
Abe H, Urao T, Ito T, Seki M, Shinozaki K, et al. (2003) Arabidopsis AtMYC2 (bHLH) and AtMYB2 (MYB) function as transcriptional activators in abscisic acid signaling. Plant Cell 15: 63–78.
[32]
Abe H, Yamaguchi-Shinozaki K, Urao T, Iwasaki T, Hosokawa D, et al. (1997) Role of Arabidopsis MYC and MYB homologs in drought- and abscisic acid-regulated gene expression. Plant Cell 9: 1859–1868.
[33]
Gocal GF, Sheldon CC, Gubler F, Moritz T, Bagnall DJ, et al. (2001) GAMYB-like genes, flowering, and gibberellin signaling in Arabidopsis. Plant Physiol 127: 1682–1693.
[34]
Jin H, Cominelli E, Bailey P, Parr A, Mehrtens F, et al. (2000) Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. EMBO J 19: 6150–6161.
[35]
Hartmann U, Sagasser M, Mehrtens F, Stracke R, Weisshaar B (2005) Differential combinatorial interactions of cis-acting elements recognized by R2R3-MYB, BZIP, and BHLH factors control light-responsive and tissue-specific activation of phenylpropanoid biosynthesis genes. Plant Mol Biol 57: 155–171.
[36]
Urao T, Noji M, Yamaguchi-Shinozaki K, Shinozaki K (1996) A transcriptional activation domain of AtMYB2, a drought-inducible Arabidopsis Myb-related protein. Plant J 10: 1145–1148.
[37]
Jung C, Seo JS, Han SW, Koo YJ, Kim CH, et al. (2008) Overexpression of AtMYB44 enhances stomatal closure to confer abiotic stress tolerance in transgenic Arabidopsis. Plant Physiol 146: 623–635.
[38]
Liu SQ, Xu L, Jia ZQ, Xu Y, Yang Q, et al. (2004) Genetic association of ETHYLENE-INSENSITIVE3-like sequence with the sex-determining M locus in cucumber (Cucumis sativus L.). Theor Appl Genet 117: 927–933.
[39]
Shi QH, Ding F, Wang XF, Wei M (2007) Exogenous nitric oxide protect cucumber roots against oxidative stress induced by salt stress. Plant Physiol Biochem 45: 542–550.
[40]
K？obus G, Janicka-Russak M (2004) Modulation by cytosolic components of proton pump activities in plasma membrane and tonoplast from Cucumis sativus roots during salt stress. Physiol Plant 12: 84–92.
[41]
Janoudi AK, Widders IE, Flore JA (1993) Water Deficits and Environmental Factors Affect Photosynthesis in Leaves of Cucumber (Cucumis sativus). J Am Soc Hortic Sci 118: 366–370.
[42]
Lee SH, Singh AP, Chung GC (2004a) Rapid accumulation of hydrogen peroxide in cucumber roots due to exposure to low temperature appears to mediate decreases in water transport. J Exp Bot 55: 1733–1741.
[43]
Lee SH, Singh AP, Chung GC, Ahn SJ, Noh EK, et al. (2004b) Exposure of roots of cucumber (Cucumis sativus) to low temperature severely reduces root pressure, hydraulic conductivity and active transport of nutrients. Physiol Plant 120: 413–420.
[44]
Huang SW, Li RQ, Zhang ZH, Li L, Gu XF, et al. (2009) The genome of the cucumber, Cucumis sativus L. Nat Genet. 475: 1–7.
[45]
Du H, Yang SS, Feng BR, Liu L, Tang YX, et al. (2012) Genome-wide analysis of the MYB transcription factor superfamily in soybean. BMC Plant Biol 12: 106.
[46]
Cheng H, Song S, Xiao L, Soo HM, Cheng Z, et al. (2009) Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLoS Genet 5: e1000440.
[47]
Mandaokar A, Browse J (2009) MYB108 acts together with MYB24 to regulate jasmonate-mediated stamen maturation in Arabidopsis. Plant Physiol. 149: 851–862.
[48]
Baumann K, Perez-Rodriguez M, Bradley D, Venail J, Bailey P, et al. (2007) Control of cell and petal morphogenesis by R2R3 MYB transcription factors. Development 134: 1691–1701.
[49]
Jakoby MJ, Falkenhan D, Mader MT, Brininstool G, Wischnitzki E, et al. (2008) Transcriptional profiling of mature Arabidopsis trichomes reveals that NOECK encodes the MIXTA like transcriptional regulator MYB106. Plant Physiol 148: 1583–1602.
[50]
Li XX, Duan XP, Jiang HX, Sun YJ, Zhang DB, et al. (2006) Genome-Wide Analysis of Basic/Helix-Loop-Helix Transcription Factor Family in Rice and Arabidopsis. Plant Physiol 141: 1167–1184.
[51]
Sharp PA (1981) Speculations on RNA splicing. Cell 23: 643–646.
[52]
Jiang CZ, Gu X, Peterson T (2004) Identification of conserved gene structures and carboxy-terminal motifs in the Myb gene family of Arabidopsis and Oryza sativa L. ssp. Indica.. Genome Biol 5: R46.
[53]
Kent WJ, Baertsch R, Hinrichs A, Miller W, Haussler D (2003) Evolution’s cauldron: duplication, deletion, and rearrangement in the mouse and human genomes. Proc Natl Acad Sci USA 100: 11484–11489.
[54]
Cannon SB, Mitra A, Baumgarten A, Young ND, May G (2004) The roles of segmental and tandem gene duplication in the evolution of large gene families in Arabidopsis thaliana. BMC Plant Biol 4: 10.
[55]
Mehan MR, Freimer NB, Ophoff RA (2004) A genome-wide survey of segmental duplications that mediate common human genetic variation of chromosomal architecture. Hum Genomics 1: 335–344.
[56]
Schauser L, Wieloch W, Stougaard J (2005) Evolution of NIN-like proteins in Arabidopsis, rice, and Lotus japonicus. J Mol Evol 60: 229–237.
[57]
Ling J, Jiang WJ, Zhang Y, Yu HJ, Mao ZC, et al. (2011) Genome-wide analysis of WRKY gene family in Cucumis sativus.. BMC Genomics 12: 471.
[58]
Hu LF, Liu SQ (2012) Genome-wide analysis of the MADS-box gene family in cucumber. Genome 55: 245–256.
[59]
Liu SQ, Liu XH, Jiang LW (2011) Genome-wide identification, phylogeny and expression analysis of the lipoxygenase gene family in cucumber. Genet Mol Res 10: 2613–2636.
[60]
Hu LF, Liu SQ (2011) Genome-wide identification and phylogenetic analysis of the ERF gene family in cucumbers. Genet Mol Biol 34: 624–633.
[61]
Holub EB (2001) The arms race is ancient history in Arabidopsis, the wildflower. Nat Rev Genet 2: 516–527.
[62]
Li JG, Li XJ, Guo L, Qu LJ, Gu HY, et al. (2006) A subgroup of MYB transcription factor genes undergoes highly conserved alternative splicing in Arabidopsis and rice. J Exp Bot 57: 1263–1273.
[63]
Wang L, Wessler SR (1998) Inefficient reinitiation is responsible for upstream open reading frame-mediated translational repression of the maize R gene. Plant Cell 10: 1733–1746.
[64]
Bailey-Serres J (1999) Selective translation of cytoplasmic mRNAs in plants. Trends Plant Sci 4: 142–148.
[65]
Locatelli F, Magnani E, Vighi C, Lanzanova C, Coraggio I (2002) Inhibitory effect of myb7 uORF on downstream gene expression in homologous (rice) and heterologous (tobacco) systems. Plant Mol Biol. 48: 309–318.
[66]
Wiese A, Elzinga N, Wobbes B, Smeekens S (2004) A conserved upstream open reading frame mediates sucrose-induced repression of translation. Plant Cell 16: 1717–1729.
[67]
Lee TH, Kim SJ, Kang SW, Lee KK, Rhee SG, et al. (2000) Molecular cloning and characterization of the mouse peroxiredoxin V gene. Biochem Biophys Res Commun 270: 356–362.
[68]
Kozak M (2002) Pushing the limits of the scanning mechanism for initiation of translation. Gene 299: 1–34.
[69]
Larsen LK, Amri EZ, Mandrup S, Pacot C, Kristiansen K (2002) Genomic organization of the mouse peroxisome proliferatoractivated receptor beta/delta gene: alternative promoter usage and splicing yield transcripts exhibiting differential translational efficiency. Biochem J 366: 767–775.
[70]
Taylor JS, Raes J (2004) Duplication and divergence: The evolution of new genes and old ideas. Annu Rev Genet 38: 615–643.
[71]
Nozawa M, Kawahara Y, Masatoshi Nei M (2007) Genomic drift and copy number variation of sensory receptor genes in humans. Proc Natl Acad Sci USA 104: 20421–20426.
[72]
Serna L, Martin C (2006) Trichomes: different regulatory networks lead to convergent structures. Trends Plant Sci 11: 274–280.
[73]
Yang CX, Li HX, Zhang JH, Luo ZD, Ye ZB, et al. (2011) A regulatory gene induces trichome formation and embryo lethality in tomato. Proc Natl Acad Sci USA 108: 11836–41.
[74]
Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C (2000) Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell 12: 2383–2394.
[75]
Gonzalez A, Zhao M, Leavitt JM, Lloyd AM (2007) Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/Myb transcriptional complex in Arabidopsis seedlings. Plant J 53: 814–27.
[76]
Bowers JE, Chapman BA, Rong J, Paterson AH (2003) Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature 422: 433–438.
[77]
Grubb CD, Abel S (2006) Glucosinolate metabolism and its control. Trends Plant Sci 11: 89–100.
[78]
Halkier BA, Gershenzon J (2006) Biology and biochemistry of glucosinolates. Annu Rev Plant Biol 57: 303–33.
[79]
Lorkovic ZJ, Wieczorek Kirk DA, Lambermon MH, Filipowicz W (2000) Pre-mRNA splicing in higher plants. Trends Plant Sci 5: 160–167.
[80]
Kazan K (2003) Alternative splicing and proteome diversity in plants: the tip of the iceberg has just emerged. Trends Plant Sci 8: 468–471.
[81]
Letunic I, Doerks T, Bork P (2009) SMART 6: recent updates and new developments. Nucleic Acids Res 37: D229–232.
[82]
Crooks GE, Hon G, Chandonia JM, Brenner SE (2004) WebLogo: a sequence logo generator. Genome Research 14: 1188–1190.
[83]
Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24: 1596–1599.
[84]
Guo AY, Zhu QH, Chen X, Luo JC (2007) GSDS: a gene structure display server. Yi Chuan 29: 1023–1026.
[85]
He HS, Dong Q, Shao YH, Jiang HY, Zhu SW, et al.. (2012) Genome-wide survey and characterization of the WRKY gene family in Populus trichocarpa. Plant Cell Rep DOI: 10.1007/s00299–012–1241–0.
[86]
Li Q, Wang XF, Ma LY, Wei M, Shi QH, et al. (2012) Molecular characterization of a cucumber nitrate reductase (CsNR) gene under NO3- stress. Mol Biol Rep 39: 4283–4290.
[87]
Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using Real-Time Quantitative PCR and the 2 ？ΔΔCt method. Methods 25: 402–408.
