Identification of novel and conserved miRNAs involved in pollen development in Brassica campestris ssp. chinensis by high-throughput sequencing and degradome analysis

BMC Genomics

Jiashu Cao¹  Zhiming Ma¹  Ying Liang¹  Meiling Lv¹  Jianxia Jiang¹ 

[1] Laboratory of Cell and Molecular Biology, Institute of Vegetable Science, Zhejiang University, Hangzhou, 310058, China

关键词: Degradome analysis;    Deep sequencing;    High-throughput sequencing;    Pollen development;    miRNAs;    Brassica rapa;    Brassica campestris;   

Others  :  1217850
DOI  :  10.1186/1471-2164-15-146

 received in 2013-09-27, accepted in 2014-02-07,  发布年份 2014

【 摘 要 】

Background

microRNAs (miRNAs) are endogenous, noncoding, small RNAs that have essential regulatory functions in plant growth, development, and stress response processes. However, limited information is available about their functions in sexual reproduction of flowering plants. Pollen development is an important process in the life cycle of a flowering plant and is a major factor that affects the yield and quality of crop seeds.

Results

This study aims to identify miRNAs involved in pollen development. Two independent small RNA libraries were constructed from the flower buds of the male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) of Brassica campestris ssp. chinensis. The libraries were subjected to high-throughput sequencing by using the Illumina Solexa system. Eight novel miRNAs on the other arm of known pre-miRNAs, 54 new conserved miRNAs, and 8 novel miRNA members were identified. Twenty-five pairs of novel miRNA/miRNA* were found. Among all the identified miRNAs, 18 differentially expressed miRNAs with over two-fold change between flower buds of male sterile line (Bcajh97-01A) and male fertile line (Bcajh97-01B) were identified. qRT-PCR analysis revealed that most of the differentially expressed miRNAs were preferentially expressed in flower buds of the male fertile line (Bcajh97-01B). Degradome analysis showed that a total of 15 genes were predicted to be the targets of seven miRNAs.

Conclusions

Our findings provide an overview of potential miRNAs involved in pollen development and interactions between miRNAs and their corresponding targets, which may provide important clues on the function of miRNAs in pollen development.

【 授权许可】

   
2014 Jiang et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

Figure 4.

Figure 5.

【 参考文献 】

• [1]McCormick S: Male gametophyte development. Plant Cell 1993, 5:1265-1275.
• [2]Goldberg RB, Beals TP, Sanders PM: Anther development: basic principles and practical applications. Plant Cell 1993, 5:1217-1229.
• [3]McCormick S: Control of male gametophyte development. Plant Cell 2004, 16(Suppl):S142-S153.
• [4]Ye Q, Zhu W, Li L, Zhang S, Yin Y, Ma H, Wang X: Brassinosteroids control male fertility by regulating the expression of key genes involved in Arabidopsis anther and pollen development. Proc Natl Acad Sci USA 2010, 107:6100-6105.
• [5]Honys D, Twell D: Comparative analysis of the Arabidopsis pollen transcriptome. Plant Physiol 2003, 132:640-652.
• [6]Dobritsa AA, Geanconteri A, Shrestha J, Carlson A, Kooyers N, Coerper D, Urbanczyk-Wochniak E, Bench BJ, Sumner LW, Swanson R, Preuss D: A large-scale genetic screen in Arabidopsis to identify genes involved in pollen exine production. Plant Physiol 2011, 157:947-970.
• [7]Huang L, Cao JS, Zhang AH, Ye YQ, Zhang YC, Liu TT: The polygalacturonase gene BcMF2 from Brassica campestris is associated with intine development. J Exp Bot 2009, 60:301-313.
• [8]Liu L-C, Cao J-S, Yu X-L, Xiang X, Fei Y-J: Expression of an Antisense BcMF3 affects microsporogenesis and pollen tube growth in arabidopsis. Agri Sci China 2006, 5:339-345.
• [9]Liu LC, Xiang X, Cao JS: [BcMF4 gene, encoding a leucine-rich repeat protein, plays a role in male fertility in Chinese cabbage-pak-choi]. Yi Chuan 2006, 28:1428-1434.
• [10]Jiang J, Yu X, Miao Y, Huang L, Yao L, Cao J: Sequence characterization and expression pattern of BcMF21, a novel gene related to pollen development in Brassica campestris ssp. chinensis. Mol Biol Rep 2012, 39:7319-7326.
• [11]Zhang Q, Huang L, Liu T, Yu X, Cao J: Functional analysis of a pollen-expressed polygalacturonase gene BcMF6 in Chinese cabbage (Brassica campestris L. ssp. chinensis Makino). Plant Cell Rep 2008, 27:1207-1215.
• [12]Huang L, Ye Y, Zhang Y, Zhang A, Liu T, Cao J: BcMF9, a novel polygalacturonase gene, is required for both Brassica campestris intine and exine formation. Ann Bot 2009, 104:1339-1351.
• [13]Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for arabidopsis lateral root development. Plant Cell 2005, 17:1376-1386.
• [14]Floyd SK, Bowman JL: Gene regulation: ancient microRNA target sequences in plants. Nature 2004, 428:485-486.
• [15]Wang JW, Park MY, Wang LJ, Koo Y, Chen XY, Weigel D, Poethig RS: miRNA control of vegetative phase change in trees. PLoS Genet 2011, 7:e1002012.
• [16]Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell 2003, 15:2730-2741.
• [17]Chen X: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science 2004, 303:2022-2025.
• [18]Yamasaki H, Abdel-Ghany SE, Cohu CM, Kobayashi Y, Shikanai T, Pilon M: Regulation of copper homeostasis by micro-RNA in Arabidopsis. J Biol Chem 2007, 282:16369-16378.
• [19]Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell 2002, 110:513-520.
• [20]Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci 2007, 12:301-309.
• [21]Katiyar-Agarwal S, Jin H: Role of small RNAs in host-microbe interactions. Annu Rev Phytopathol 2010, 48:225-246.
• [22]Wei LQ, Yan LF, Wang T: Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa. Genome Biol 2011, 12:R53. BioMed Central Full Text
• [23]Chambers C, Shuai B: Profiling microRNA expression in Arabidopsis pollen using microRNA array and real-time PCR. BMC Plant Biol 2009, 9:87. BioMed Central Full Text
• [24]Grant-Downton R, Le Trionnaire G, Schmid R, Rodriguez-Enriquez J, Hafidh S, Mehdi S, Twell D, Dickinson H: MicroRNA and tasiRNA diversity in mature pollen of Arabidopsis thaliana. BMC Genomics 2009, 10:643. BioMed Central Full Text
• [25]Axtell MJ, Bowman JL: Evolution of plant microRNAs and their targets. Trends Plant Sci 2008, 13:343-349.
• [26]Wang XJ, Reyes JL, Chua NH, Gaasterland T: Prediction and identification of Arabidopsis thaliana microRNAs and their mRNA targets. Genome Biol 2004, 5:R65. BioMed Central Full Text
• [27]Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol 2006, 57:19-53.
• [28]Yu X, Wang H, Lu Y, de Ruiter M, Cariaso M, Prins M, van Tunen A, He Y: Identification of conserved and novel microRNAs that are responsive to heat stress in Brassica rapa. J Exp Bot 2012, 63:1025-1038.
• [29]He XF, Fang YY, Feng L, Guo HS: Characterization of conserved and novel microRNAs and their targets, including a TuMV-induced TIR-NBS-LRR class R gene-derived novel miRNA in Brassica. FEBS Lett 2008, 582:2445-2452.
• [30]Zhang J, Zhang S, Han S, Wu T, Li X, Li W, Qi L: Genome-wide identification of microRNAs in larch and stage-specific modulation of 11 conserved microRNAs and their targets during somatic embryogenesis. Planta 2012, 236:647-657.
• [31]Li B, Qin Y, Duan H, Yin W, Xia X: Genome-wide characterization of new and drought stress responsive microRNAs in Populus euphratica. J Exp Bot 2011, 62:3765-3779.
• [32]Heisel SE, Zhang Y, Allen E, Guo L, Reynolds TL, Yang X, Kovalic D, Roberts JK: Characterization of unique small RNA populations from rice grain. PLoS One 2008, 3:e2871.
• [33]Zhao M, Tai H, Sun S, Zhang F, Xu Y, Li WX: Cloning and characterization of maize miRNAs involved in responses to nitrogen deficiency. PLoS One 2012, 7:e29669.
• [34]Zeng QY, Yang CY, Ma QB, Li XP, Dong WW, Nian H: Identification of wild soybean miRNAs and their target genes responsive to aluminum stress. BMC Plant Biol 2012, 12:182. BioMed Central Full Text
• [35]Ma ZR, Coruh C, Axtell MJ: Arabidopsis lyrata Small RNAs: Transient MIRNA and Small Interfering RNA Loci within the Arabidopsis Genus. Plant Cell 2010, 22:1090-1103.
• [36]Cuperus JT, Fahlgren N, Carrington JC: Evolution and functional diversification of MIRNA genes. Plant Cell 2011, 23:431-442.
• [37]Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev 2006, 20:3407-3425.
• [38]Grant-Downton R, Hafidh S, Twell D, Dickinson HG: Small RNA pathways are present and functional in the angiosperm male gametophyte. Mol Plant 2009, 2:500-512.
• [39]Dai X, Zhao PX: psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res 2011, 39:W155-W159.
• [40]Xie F, Zhang B: Target-align: a tool for plant microRNA target identification. Bioinformatics 2010, 26:3002-3003.
• [41]Bonnet E, He Y, Billiau K, Van de Peer Y: TAPIR, a web server for the prediction of plant microRNA targets, including target mimics. Bioinformatics 2010, 26:1566-1568.
• [42]Zhang Z, Yu J, Li D, Liu F, Zhou X, Wang T, Ling Y, Su Z: PMRD: plant microRNA database. Nucleic Acids Res 2010, 38:D806-D813.
• [43]German MA, Pillay M, Jeong DH, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol 2008, 26:941-946.
• [44]Mao W, Li Z, Xia X, Li Y, Yu J: A combined approach of high-throughput sequencing and degradome analysis reveals tissue specific expression of microRNAs and their targets in cucumber. PLoS One 2012, 7:e33040.
• [45]Yang J, Liu X, Xu B, Zhao N, Yang X, Zhang M: Identification of miRNAs and their targets using high-throughput sequencing and degradome analysis in cytoplasmic male-sterile and its maintainer fertile lines of Brassica juncea. BMC Genomics 2013, 14:9. BioMed Central Full Text
• [46]Wu G, Park MY, Conway SR, Wang JW, Weigel D, Poethig RS: The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell 2009, 138:750-759.
• [47]Xing S, Salinas M, Hohmann S, Berndtgen R, Huijser P: miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis. Plant Cell 2010, 22:3935-3950.
• [48]Alonso-Peral MM, Li J, Li Y, Allen RS, Schnippenkoetter W, Ohms S, White RG, Millar AA: The microRNA159-regulated GAMYB-like genes inhibit growth and promote programmed cell death in Arabidopsis. Plant Physiol 2010, 154:757-771.
• [49]Huang L, Cao J, Ye W, Liu T, Jiang L, Ye Y: Transcriptional differences between the male-sterile mutant bcms and wild-type Brassica campestris ssp chinensis reveal genes related to pollen development. Plant Biology 2008, 10:342-355.
• [50]Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res 2003, 31:3406-3415.
• [51]Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant MicroRNAs. Plant Cell 2008, 20:3186-3190.
• [52]Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci 2006, 63:246-254.
• [53]Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ: Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol 2008, 18:758-762.
• [54]Addo-Quaye C, Miller W, Axtell MJ: CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics 2009, 25:130-131.
• [55]Jiang JX, Jiang JJ, Yang YF, Cao JS: Identification of microRNAs potentially involved in male sterility of Brassica campestris ssp chinensis using microRNA array and quantitative RT-PCR assays. Cell Mole Biol Lett 2013, 18:416-432.
• [56]Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25:402-408.