A comparative analysis of small RNAs between two Upland cotton backcross inbred lines with different fiber length: Expression and distribution
|
摘要
The cotton fiber is the most important raw material for the textile industry and an ideal model system for studying cell elongation. However, the genetic variation of fiber elongation in relation to miRNA is poorly understood. A high-throughput comparative RNA-seq of two lines differing in fiber length(FL) from a backcross inbred line(BIL)population of G. hirsutum × G. barbadense revealed differentially expressed(DE) miRNAs and their targets in rapidly elongating fibers. A real-time quantitative PCR analysis was further performed to validate the results. A total of 463(including 47 DE) miRNAs were identified,and seven DE miRNAs were co-localized with seven FL quantitative trait loci(QTL)identified in the G. hirsutum × G. barbadense population. Of 82(including 21 DE) targets identified, nine(including one DE) were also co-localized with the seven FL QTL. The relationship between the allopolyploid and its diploid ancestral species with respect to miRNAs and their targets was also characterized. These results will facilitate the understanding of the molecular genetic mechanism of fiber elongation with regards to miRNAs in cotton.
The cotton fiber is the most important raw material for the textile industry and an ideal model system for studying cell elongation. However, the genetic variation of fiber elongation in relation to miRNA is poorly understood. A high-throughput comparative RNA-seq of two lines differing in fiber length(FL) from a backcross inbred line(BIL)population of G. hirsutum × G. barbadense revealed differentially expressed(DE) miRNAs and their targets in rapidly elongating fibers. A real-time quantitative PCR analysis was further performed to validate the results. A total of 463(including 47 DE) miRNAs were identified,and seven DE miRNAs were co-localized with seven FL quantitative trait loci(QTL)identified in the G. hirsutum × G. barbadense population. Of 82(including 21 DE) targets identified, nine(including one DE) were also co-localized with the seven FL QTL. The relationship between the allopolyploid and its diploid ancestral species with respect to miRNAs and their targets was also characterized. These results will facilitate the understanding of the molecular genetic mechanism of fiber elongation with regards to miRNAs in cotton.
The cotton fiber is the most important raw material for the textile industry and an ideal model system for studying cell elongation. However, the genetic variation of fiber elongation in relation to miRNA is poorly understood. A high-throughput comparative RNA-seq of two lines differing in fiber length(FL) from a backcross inbred line(BIL)population of G. hirsutum × G. barbadense revealed differentially expressed(DE) miRNAs and their targets in rapidly elongating fibers. A real-time quantitative PCR analysis was further performed to validate the results. A total of 463(including 47 DE) miRNAs were identified,and seven DE miRNAs were co-localized with seven FL quantitative trait loci(QTL)identified in the G. hirsutum × G. barbadense population. Of 82(including 21 DE) targets identified, nine(including one DE) were also co-localized with the seven FL QTL. The relationship between the allopolyploid and its diploid ancestral species with respect to miRNAs and their targets was also characterized. These results will facilitate the understanding of the molecular genetic mechanism of fiber elongation with regards to miRNAs in cotton.
引文
[1]A.S.Basra,C.P.Malik,Development of the cotton fiber,Int.Rev.Cytol.89(1984)65-113.
[2]H.J.Kim,B.A.Triplett,Cotton fiber growth in planta and in vitro.Models for plant cell elongation and cell wall biogenesis,Plant Physiol.127(2001)1361-1366.
[3]J.J.Lee,A.W.Woodward,Z.J.Chen,Gene expression changes and early events in cotton fibre development,Ann.Bot.100(2007)1391-1401.
[4]Y.L.Ruan,D.J.Llewellyn,R.T.Furbank,The control of single-celled cotton fiber elongation by developmentally reversible gating of plasmodesmata and coordinated expression of sucrose and K+transporters and expansin,Plant Cell 13(2001)47-60.
[5]Y.M.Qin,Y.X.Zhu,How cotton fibers elongate:a tale of linear cell-growth mode,Curr.Opin.Plant Biol.14(2011)106-111.
[6]Y.Q.Huang,J.Wang,L.D.Zhang,K.J.Zuo,A cotton annexin protein AnxGb6 regulates fiber elongation through its interaction with actin 1,PLoS One 8(2013),e66160.
[7]M.Zhang,X.L.Zheng,S.Q.Song,Q.W.Zeng,L.Hou,D.M.Li,J.Zhao,Y.Wei,X.B.Li,M.Luo,Y.H.Xiao,X.Y.Luo,J.F.Zhang,C.B.Xiang,Y.Pei,Spatiotemporal manipulation of auxin biosynthesis in cotton ovule epidermal cells enhances fiber yield and quality,Nat.Biotechnol.29(2011)453-458.
[8]S.J.Ji,Y.C.Lu,J.X.Feng,G.Wei,J.Li,Y.H.Shi,Q.Fu,D.Liu,J.C.Luo,Y.X.Zhu,Isolation and analyses of genes preferentially expressed during early cotton fiber development by subtractive PCR and cDNA array,Nucleic Acids Res.31(2003)2534-2543.
[9]P.B.Kwak,Q.Q.Wang,X.S.Chen,C.X.Qiu,Z.M.Yang,Enrichment of a set of microRNAs during the cotton fiber development,BMC Genomics 10(2009)457.
[10]M.X.Pang,A.W.Woodward,V.Agarwal,X.Y.Guan,M.Ha,V.Ramachandran,X.M.Chen,B.A.Triplett,D.M.Stelly,Z.J.Chen,Genome-wide analysis reveals rapid and dynamic changes in miRNA and siRNA sequence and expression during ovule and fiber development in allotetraploid cotton(Gossypium hirsutum L.),Genome Biol.10(2009)R122.
[11]Z.M.Wang,W.Xue,C.J.Dong,L.G.Jin,S.M.Bian,C.Wang,X.Y.Wu,J.Y.Liu,A comparative miRNAome analysis reveals seven fiber initiation-related and 36 novel miRNAs in developing cotton ovules,Mol.Plant 5(2012)889-900.
[12]H.Zhang,Q.Wan,W.X.Ye,Y.D.Lv,H.T.Wu,T.Z.Zhang,Genome-wide analysis of small rNA and novel microRNAdiscovery during fiber and seed initial development in Gossypium hirsutum.L.PLoS One 8(2013)e69743.
[13]Q.Li,X.Jin,Y.X.Zhu,Identification and analyses of miRNAgenes in allotetraploid Gossypium hirsutum fiber cells based on the sequenced diploid G.raimondii genome,J.Genet.Genomics 39(2012)351-360.
[14]W.Xue,Z.M.Wang,M.J.Du,Y.D.Liu,J.Y.Liu,Genome-wide analysis of small RNAs reveals eight fiber elongation-related and 257 novel microRNAs in elongating cotton fiber cells,BMC Genomics 14(2013)629.
[15]X.Y.Guan,M.X.Pang,G.Nah,X.L.Shi,W.X.Ye,D.M.Stelly,Z.J.Chen,miR828 and miR858 regulate homoeologous MYB2gene functions in Arabidopsis trichome and cotton fibre development,Nat.Commun.5(2014)3050.
[16]N.Liu,L.Tu,L.Wang,H.Hu,J.Xu,X.Zhang,MicroRNA 157-targeted SPL genes regulate floral organ size and ovule production in cotton,BMC Plant Biol.17(2017)7.
[17]J.W.Yu,K.Zhang,S.Y.Li,S.X.Yu,H.H.Zhai,M.Wu,X.L.Li,S.L.Fan,M.Z.Song,D.G.Yang,Y.H.Li,J.F.Zhang,Mapping quantitative trait loci for lint yield and fiber quality across environments in a Gossypium hirsutum×Gossypium barbadense backcross inbred line population,Theor.Appl.Genet.126(2013)275-287.
[18]X.M.Chen,W.H.Gao,J.F.Zhang,X.L.Zhang,Z.X.Lin,Linkage mapping and expression analysis of miRNAs and their target genes during fiber development in cotton,BMC Genomics 14(2013)706.
[19]M.X.Pang,S.F.Hanson,J.F.Zhang,Cloning and expression studies of novel small RNAs in tetraploid cotton,Plant Mol.Biol.Report.30(2012)710-718.
[20]M.X.Pang,C.Z.Xing,N.Adams,L.Rodriguez-Uribe,S.E.Hughs,S.F.Hanson,J.F.Zhang,Comparative expression of miRNA genes and miRNA-based AFLP marker analysis in cultivated tetraploid cottons,J.Plant Physiol.168(2011)824-830.
[21]R.Q.Li,Y.R.Li,K.Kristiansen,J.Wang,SOAP:short oligonucleotide alignment program,Bioinformatics 24(2008)713-714.
[22]B.C.Meyers,M.J.Axtell,B.Bartel,D.P.Bartel,D.Baulcombe,J.L.Bowman,X.F.Cao,J.C.Carrington,X.M.Chen,P.J.Green,S.Griffiths-Jones,S.E.Jacobsen,A.C.Mallory,R.A.Martienssen,R.S.Poethig,Y.J.Qi,H.Vaucheret,O.Voinnet,Y.Watanabe,D.Weigel,J.K.Zhu,Criteria for annotation of plant microRNAs,Plant Cell 20(2008)3186-3190.
[23]X.Y.Yang,L.C.Wang,D.J.Yuan,K.Lindsey,X.L.Zhang,Small RNA and degradome sequencing reveal complex miRNAregulation during cotton somatic embryogenesis,J.Exp.Bot.64(2013)1521-1536.
[24]J.I.Said,J.A.Knapka,M.Z.Song,J.F.Zhang,Cotton QTLdb:a cotton QTL database for QTL analysis,visualization,and comparison between Gossypium hirsutum and G.hirsutum×G.barbadense populations,Mol.Gen.Genomics.290(2015)1615-1625.
[25]J.I.Said,M.Z.Song,H.T.Wang,Z.X.Lin,X.L.Zhang,D.D.Fang,J.F.Zhang,A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G.hirsutum×G.barbadense populations,Mol.Gen.Genomics.290(2015)1003-1025.
[26]E.Allen,Z.X.Xie,A.M.Gustafson,J.C.Carrington,microRNA-directed phasing during trans-acting siRNA biogenesis in plants,Cell 121(2005)207-221.
[27]R.Schwab,J.F.Palatnik,M.Riester,C.Schommer,M.Schmid,D.Weigel,Specific effects of microRNAs on the plant transcriptome,Dev.Cell 8(2005)517-527.
[28]A.Conesa,S.Gotz,J.M.Garcia-Gomez,J.Terol,M.Talon,M.Robles,Blast2GO:a universal tool for annotation,visualization and analysis in functional genomics research,Bioinformatics21(2005)3674-3676.
[29]M.A.German,M.Pillay,D.H.Jeong,A.Hetawal,S.J.Luo,P.Janardhanan,V.Kannan,L.A.Rymarquis,K.Nobuta,R.German,E.De Paoli,C.Lu,G.Schroth,B.C.Meyers,P.J.Green,Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends,Nat.Biotechnol.26(2008)941-946.
[30]C.F.Chen,D.A.Ridzon,A.J.Broomer,Z.H.Zhou,D.H.Lee,J.T.Nguyen,M.Barbisin,N.L.Xu,V.R.Mahuvakar,M.R.Andersen,K.Q.Lao,K.J.Livak,K.J.Guegler,Real-time quantification of microRNAs by stem-loop RT-PCR,Nucleic Acids Res.33(2005)e179.
[31]X.H.Li,M.Wu,G.Y.Liu,W.F.Pei,H.H.Zhai,J.W.Yu,J.F.Zhang,S.X.Yu,Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton,BMC Genomics 18(2017)427.
[32]R.D.Morin,G.Aksay,E.Dolgosheina,H.A.Ebhardt,V.Magrini,E.R.Mardis,S.C.Sahinalp,P.J.Unrau,Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa,Genome Res.18(2008)571-584.
[33][33]C.Z.Zhao,H.Xia,T.P.Frazier,Y.Y.Yao,Y.P.Bi,A.Q.Li,M.J.Li,C.S.Li,B.H.Zhang,X.J.Wang,Deep sequencing identifies novel and conserved microRNAs in peanuts(Arachis hypogaea L.),BMC Plant Biol.10(2010)3.
[34]Q.X.Song,Y.F.Liu,X.Y.Hu,W.K.Zhang,B.A.Ma,S.Y.Chen,J.S.Zhang,Identification of miRNAs and their target genes in developing soybean seeds by deep sequencing,BMC Plant Biol.11(2011)5.
[35]G.Szittya,S.Moxon,D.M.Santos,R.Jing,M.P.S.Fevereiro,V.Moulton,T.Dalmay,High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNAfamilies,BMC Genomics 9(2008)593.
[36]Z.J.Yin,Y.Li,X.L.Han,F.F.Shen,Genome-wide profiling of miRNAs and other small non-coding RNAs in the Verticillium dahliae-inoculated cotton roots,PLoS One 7(2012),e35765.
[37]F.G.Li,G.Y.Fan,C.R.Lu,G.H.Xiao,C.S.Zou,R.J.Kohel,Z.Y.Ma,H.H.Shang,X.F.Ma,J.Y.Wu,X.M.Liang,G.Huang,R.G.Percy,K.Liu,W.H.Yang,W.B.Chen,X.M.Du,C.C.Shi,Y.L.Yuan,W.W.Ye,X.Liu,X.Y.Zhang,W.Q.Liu,H.L.Wei,S.J.Wei,G.D.Huang,X.L.Zhang,S.J.Zhu,H.Zhang,F.M.Sun,X.F.Wang,J.Liang,J.H.Wang,Q.He,L.H.Huang,J.Wang,J.J.Cui,G.L.Song,K.B.Wang,X.Xu,J.Z.Yu,Y.X.Zhu,S.X.Yu,Genome sequence of cultivated Upland cotton(Gossypium hirsutum TM-1)provides insights into genome evolution,Nat.Biotechnol.33(2015)524-530.
[38]T.Z.Zhang,Y.Hu,W.K.Jiang,L.Fang,X.Y.Guan,J.D.Chen,J.B.Zhang,C.A.Saski,B.E.Scheffler,D.M.Stelly,A.M.HulseKemp,Q.Wan,B.L.Liu,C.X.Liu,S.Wang,M.Q.Pan,Y.K.Wang,D.W.Wang,W.X.Ye,L.J.Chang,W.P.Zhang,Q.X.Song,R.C.Kirkbride,X.Y.Chen,E.Dennis,D.J.Llewellyn,D.G.Peterson,P.Thaxton,D.C.Jones,Q.Wang,X.Y.Xu,H.Zhang,H.T.Wu,L.Zhou,G.F.Mei,S.Q.Chen,Y.Tian,D.Xiang,X.H.Li,J.Ding,Q.Y.Zuo,L.N.Tao,Y.C.Liu,J.Li,Y.Lin,Y.Y.Hui,Z.S.Cao,C.P.Cai,X.F.Zhu,Z.Jiang,B.L.Zhou,W.Z.Guo,R.Q.Li,Z.J.Chen,Sequencing of allotetraploid cotton(Gossypium hirsutum L.acc.TM-1)provides a resource for fiber improvement,Nat.Biotechnol.33(2015)531-537.
[39]K.B.Wang,Z.W.Wang,F.G.Li,W.W.Ye,J.Y.Wang,G.L.Song,Z.Yue,L.Cong,H.H.Shang,S.L.Zhu,C.S.Zou,Q.Li,Y.L.Yuan,C.R.Lu,H.L.Wei,C.Y.Gou,Z.Q.Zheng,Y.Yin,X.Y.Zhang,K.Liu,B.Wang,C.Song,N.Shi,R.J.Kohel,R.G.Percy,J.Z.Yu,Y.X.Zhu,J.Wang,S.X.Yu,The draft genome of a diploid cotton Gossypium raimondii,Nat.Genet.44(2012)1098-1103.
[40]F.G.Li,G.Y.Fan,K.B.Wang,F.M.Sun,Y.L.Yuan,G.L.Song,Q.Li,Z.Y.Ma,C.R.Lu,C.S.Zou,W.B.Chen,X.M.Liang,H.H.Shang,W.Q.Liu,C.C.Shi,G.H.Xiao,C.Y.Gou,W.W.Ye,X.Xu,X.Y.Zhang,H.L.Wei,Z.F.Li,G.Y.Zhang,J.Y.Wang,K.Liu,R.J.Kohel,R.G.Percy,J.Z.Yu,Y.X.Zhu,J.Wang,S.X.Yu,Genome sequence of the cultivated cotton Gossypium arboreum,Nat.Genet.46(2014)567-572.
[41]S.F.Altschul,W.Gish,W.Miller,E.W.Myers,D.J.Lipman,Basic local alignment search tool,J.Mol.Biol.215(1990)403-410.
[42]C.Addo-Quaye,T.W.Eshoo,D.P.Bartel,M.J.Axtell,Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome,Curr.Biol.18(2008)758-762.
[43]S.Griffiths-Jones,H.K.Saini,S.van Dongen,A.J.Enright,miRBase:tools for microRNA genomics,Nucleic Acids Res.36(2008)D154-D158.
[44]O.Voinnet,Origin,biogenesis,and activity of plant microRNAs,Cell 136(2009)669-687.
[45]N.Liu,L.L.Tu,W.X.Tang,W.H.Gao,K.Lindsey,X.L.Zhang,Small RNA and degradome profiling reveals a role for miRNAs and their targets in the developing fibers of Gossypium barbadense,Plant J.80(2014)331-344.
[46]S.J.Szalma,B.M.Hostert,J.R.LeDeaux,C.W.Stuber,J.B.Holland,QTL mapping with near-isogenic lines in maize,Theor.Appl.Genet.114(2007)1211-1228.
[47]R.Kooke,E.Wijnker,J.J.Keurentjes,Backcross populations and near isogenic lines,Methods Mol.Biol.871(2012)3-16.
[48]L.Gong,A.Kakrana,S.Arikit,B.C.Meyers,J.F.Wendel,Composition and expression of conserved microRNA genes in diploid cotton(Gossypium)species,Genome Biol.Evol.5(2013)2449-2459.
[49]X.Y.Guan,Q.X.Song,Z.J.Chen,Polyploidy and small RNAregulation of cotton fiber development,Trends Plant Sci.19(2014)516-528.
[50]G.Baum,S.Lev-Yadun,Y.Fridmann,T.Arazi,H.Katsnelson,M.Zik,H.Fromm,Calmodulin binding to glutamate decarboxylase is required for regulation of glutamate and GABA metabolism and normal development in plants,EMBO J.15(1996)2988-2996.
[51]M.L.Preuss,D.B.Delmer,The cotton kinesin-like calmodulinbinding protein associates with cortical microtubules in cotton fibers,Plant Physiol.132(2003)154-160.
[52]C.M.Shan,X.X.Shangguan,B.Zhao,X.F.Zhang,L.M.Chao,C.Q.Yang,L.J.Wang,H.Y.Zhu,Y.D.Zeng,W.Z.Guo,B.L.Zhou,G.J.Hu,X.Y.Guan,Z.J.Chen,J.F.Wendel,T.Z.Zhang,X.Y.Chen,Control of cotton fibre elongation by a homeodomain transcription factor GhHOX3,Nat.Commun.5(2014)5519.
[53]A.L.Silverstone,H.S.Jung,A.Dill,H.Kawaide,Y.Kamiya,T.P.Sun,Repressing a repressor:gibberellin-induced rapid reduction of the RGA protein in Arabidopsis,Plant Cell 13(2001)1555-1565.
[54]G.Hagen,T.Guilfoyle,Auxin-responsive gene expression:genes,promoters and regulatory factors,Plant Mol.Biol.49(2002)373-385.
[55]E.W.Wang,B.M.Kessler,A.Borodovsky,B.F.Cravatt,M.Bogyo,H.L.Ploegh,R.Glas,Integration of the ubiquitin-proteasome pathway with a cytosolic oligopeptidase activity,Proc.Natl.Acad.Sci.U.S.A.97(2000)9990-9995.
[2]H.J.Kim,B.A.Triplett,Cotton fiber growth in planta and in vitro.Models for plant cell elongation and cell wall biogenesis,Plant Physiol.127(2001)1361-1366.
[3]J.J.Lee,A.W.Woodward,Z.J.Chen,Gene expression changes and early events in cotton fibre development,Ann.Bot.100(2007)1391-1401.
[4]Y.L.Ruan,D.J.Llewellyn,R.T.Furbank,The control of single-celled cotton fiber elongation by developmentally reversible gating of plasmodesmata and coordinated expression of sucrose and K+transporters and expansin,Plant Cell 13(2001)47-60.
[5]Y.M.Qin,Y.X.Zhu,How cotton fibers elongate:a tale of linear cell-growth mode,Curr.Opin.Plant Biol.14(2011)106-111.
[6]Y.Q.Huang,J.Wang,L.D.Zhang,K.J.Zuo,A cotton annexin protein AnxGb6 regulates fiber elongation through its interaction with actin 1,PLoS One 8(2013),e66160.
[7]M.Zhang,X.L.Zheng,S.Q.Song,Q.W.Zeng,L.Hou,D.M.Li,J.Zhao,Y.Wei,X.B.Li,M.Luo,Y.H.Xiao,X.Y.Luo,J.F.Zhang,C.B.Xiang,Y.Pei,Spatiotemporal manipulation of auxin biosynthesis in cotton ovule epidermal cells enhances fiber yield and quality,Nat.Biotechnol.29(2011)453-458.
[8]S.J.Ji,Y.C.Lu,J.X.Feng,G.Wei,J.Li,Y.H.Shi,Q.Fu,D.Liu,J.C.Luo,Y.X.Zhu,Isolation and analyses of genes preferentially expressed during early cotton fiber development by subtractive PCR and cDNA array,Nucleic Acids Res.31(2003)2534-2543.
[9]P.B.Kwak,Q.Q.Wang,X.S.Chen,C.X.Qiu,Z.M.Yang,Enrichment of a set of microRNAs during the cotton fiber development,BMC Genomics 10(2009)457.
[10]M.X.Pang,A.W.Woodward,V.Agarwal,X.Y.Guan,M.Ha,V.Ramachandran,X.M.Chen,B.A.Triplett,D.M.Stelly,Z.J.Chen,Genome-wide analysis reveals rapid and dynamic changes in miRNA and siRNA sequence and expression during ovule and fiber development in allotetraploid cotton(Gossypium hirsutum L.),Genome Biol.10(2009)R122.
[11]Z.M.Wang,W.Xue,C.J.Dong,L.G.Jin,S.M.Bian,C.Wang,X.Y.Wu,J.Y.Liu,A comparative miRNAome analysis reveals seven fiber initiation-related and 36 novel miRNAs in developing cotton ovules,Mol.Plant 5(2012)889-900.
[12]H.Zhang,Q.Wan,W.X.Ye,Y.D.Lv,H.T.Wu,T.Z.Zhang,Genome-wide analysis of small rNA and novel microRNAdiscovery during fiber and seed initial development in Gossypium hirsutum.L.PLoS One 8(2013)e69743.
[13]Q.Li,X.Jin,Y.X.Zhu,Identification and analyses of miRNAgenes in allotetraploid Gossypium hirsutum fiber cells based on the sequenced diploid G.raimondii genome,J.Genet.Genomics 39(2012)351-360.
[14]W.Xue,Z.M.Wang,M.J.Du,Y.D.Liu,J.Y.Liu,Genome-wide analysis of small RNAs reveals eight fiber elongation-related and 257 novel microRNAs in elongating cotton fiber cells,BMC Genomics 14(2013)629.
[15]X.Y.Guan,M.X.Pang,G.Nah,X.L.Shi,W.X.Ye,D.M.Stelly,Z.J.Chen,miR828 and miR858 regulate homoeologous MYB2gene functions in Arabidopsis trichome and cotton fibre development,Nat.Commun.5(2014)3050.
[16]N.Liu,L.Tu,L.Wang,H.Hu,J.Xu,X.Zhang,MicroRNA 157-targeted SPL genes regulate floral organ size and ovule production in cotton,BMC Plant Biol.17(2017)7.
[17]J.W.Yu,K.Zhang,S.Y.Li,S.X.Yu,H.H.Zhai,M.Wu,X.L.Li,S.L.Fan,M.Z.Song,D.G.Yang,Y.H.Li,J.F.Zhang,Mapping quantitative trait loci for lint yield and fiber quality across environments in a Gossypium hirsutum×Gossypium barbadense backcross inbred line population,Theor.Appl.Genet.126(2013)275-287.
[18]X.M.Chen,W.H.Gao,J.F.Zhang,X.L.Zhang,Z.X.Lin,Linkage mapping and expression analysis of miRNAs and their target genes during fiber development in cotton,BMC Genomics 14(2013)706.
[19]M.X.Pang,S.F.Hanson,J.F.Zhang,Cloning and expression studies of novel small RNAs in tetraploid cotton,Plant Mol.Biol.Report.30(2012)710-718.
[20]M.X.Pang,C.Z.Xing,N.Adams,L.Rodriguez-Uribe,S.E.Hughs,S.F.Hanson,J.F.Zhang,Comparative expression of miRNA genes and miRNA-based AFLP marker analysis in cultivated tetraploid cottons,J.Plant Physiol.168(2011)824-830.
[21]R.Q.Li,Y.R.Li,K.Kristiansen,J.Wang,SOAP:short oligonucleotide alignment program,Bioinformatics 24(2008)713-714.
[22]B.C.Meyers,M.J.Axtell,B.Bartel,D.P.Bartel,D.Baulcombe,J.L.Bowman,X.F.Cao,J.C.Carrington,X.M.Chen,P.J.Green,S.Griffiths-Jones,S.E.Jacobsen,A.C.Mallory,R.A.Martienssen,R.S.Poethig,Y.J.Qi,H.Vaucheret,O.Voinnet,Y.Watanabe,D.Weigel,J.K.Zhu,Criteria for annotation of plant microRNAs,Plant Cell 20(2008)3186-3190.
[23]X.Y.Yang,L.C.Wang,D.J.Yuan,K.Lindsey,X.L.Zhang,Small RNA and degradome sequencing reveal complex miRNAregulation during cotton somatic embryogenesis,J.Exp.Bot.64(2013)1521-1536.
[24]J.I.Said,J.A.Knapka,M.Z.Song,J.F.Zhang,Cotton QTLdb:a cotton QTL database for QTL analysis,visualization,and comparison between Gossypium hirsutum and G.hirsutum×G.barbadense populations,Mol.Gen.Genomics.290(2015)1615-1625.
[25]J.I.Said,M.Z.Song,H.T.Wang,Z.X.Lin,X.L.Zhang,D.D.Fang,J.F.Zhang,A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G.hirsutum×G.barbadense populations,Mol.Gen.Genomics.290(2015)1003-1025.
[26]E.Allen,Z.X.Xie,A.M.Gustafson,J.C.Carrington,microRNA-directed phasing during trans-acting siRNA biogenesis in plants,Cell 121(2005)207-221.
[27]R.Schwab,J.F.Palatnik,M.Riester,C.Schommer,M.Schmid,D.Weigel,Specific effects of microRNAs on the plant transcriptome,Dev.Cell 8(2005)517-527.
[28]A.Conesa,S.Gotz,J.M.Garcia-Gomez,J.Terol,M.Talon,M.Robles,Blast2GO:a universal tool for annotation,visualization and analysis in functional genomics research,Bioinformatics21(2005)3674-3676.
[29]M.A.German,M.Pillay,D.H.Jeong,A.Hetawal,S.J.Luo,P.Janardhanan,V.Kannan,L.A.Rymarquis,K.Nobuta,R.German,E.De Paoli,C.Lu,G.Schroth,B.C.Meyers,P.J.Green,Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends,Nat.Biotechnol.26(2008)941-946.
[30]C.F.Chen,D.A.Ridzon,A.J.Broomer,Z.H.Zhou,D.H.Lee,J.T.Nguyen,M.Barbisin,N.L.Xu,V.R.Mahuvakar,M.R.Andersen,K.Q.Lao,K.J.Livak,K.J.Guegler,Real-time quantification of microRNAs by stem-loop RT-PCR,Nucleic Acids Res.33(2005)e179.
[31]X.H.Li,M.Wu,G.Y.Liu,W.F.Pei,H.H.Zhai,J.W.Yu,J.F.Zhang,S.X.Yu,Identification of candidate genes for fiber length quantitative trait loci through RNA-Seq and linkage and physical mapping in cotton,BMC Genomics 18(2017)427.
[32]R.D.Morin,G.Aksay,E.Dolgosheina,H.A.Ebhardt,V.Magrini,E.R.Mardis,S.C.Sahinalp,P.J.Unrau,Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa,Genome Res.18(2008)571-584.
[33][33]C.Z.Zhao,H.Xia,T.P.Frazier,Y.Y.Yao,Y.P.Bi,A.Q.Li,M.J.Li,C.S.Li,B.H.Zhang,X.J.Wang,Deep sequencing identifies novel and conserved microRNAs in peanuts(Arachis hypogaea L.),BMC Plant Biol.10(2010)3.
[34]Q.X.Song,Y.F.Liu,X.Y.Hu,W.K.Zhang,B.A.Ma,S.Y.Chen,J.S.Zhang,Identification of miRNAs and their target genes in developing soybean seeds by deep sequencing,BMC Plant Biol.11(2011)5.
[35]G.Szittya,S.Moxon,D.M.Santos,R.Jing,M.P.S.Fevereiro,V.Moulton,T.Dalmay,High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNAfamilies,BMC Genomics 9(2008)593.
[36]Z.J.Yin,Y.Li,X.L.Han,F.F.Shen,Genome-wide profiling of miRNAs and other small non-coding RNAs in the Verticillium dahliae-inoculated cotton roots,PLoS One 7(2012),e35765.
[37]F.G.Li,G.Y.Fan,C.R.Lu,G.H.Xiao,C.S.Zou,R.J.Kohel,Z.Y.Ma,H.H.Shang,X.F.Ma,J.Y.Wu,X.M.Liang,G.Huang,R.G.Percy,K.Liu,W.H.Yang,W.B.Chen,X.M.Du,C.C.Shi,Y.L.Yuan,W.W.Ye,X.Liu,X.Y.Zhang,W.Q.Liu,H.L.Wei,S.J.Wei,G.D.Huang,X.L.Zhang,S.J.Zhu,H.Zhang,F.M.Sun,X.F.Wang,J.Liang,J.H.Wang,Q.He,L.H.Huang,J.Wang,J.J.Cui,G.L.Song,K.B.Wang,X.Xu,J.Z.Yu,Y.X.Zhu,S.X.Yu,Genome sequence of cultivated Upland cotton(Gossypium hirsutum TM-1)provides insights into genome evolution,Nat.Biotechnol.33(2015)524-530.
[38]T.Z.Zhang,Y.Hu,W.K.Jiang,L.Fang,X.Y.Guan,J.D.Chen,J.B.Zhang,C.A.Saski,B.E.Scheffler,D.M.Stelly,A.M.HulseKemp,Q.Wan,B.L.Liu,C.X.Liu,S.Wang,M.Q.Pan,Y.K.Wang,D.W.Wang,W.X.Ye,L.J.Chang,W.P.Zhang,Q.X.Song,R.C.Kirkbride,X.Y.Chen,E.Dennis,D.J.Llewellyn,D.G.Peterson,P.Thaxton,D.C.Jones,Q.Wang,X.Y.Xu,H.Zhang,H.T.Wu,L.Zhou,G.F.Mei,S.Q.Chen,Y.Tian,D.Xiang,X.H.Li,J.Ding,Q.Y.Zuo,L.N.Tao,Y.C.Liu,J.Li,Y.Lin,Y.Y.Hui,Z.S.Cao,C.P.Cai,X.F.Zhu,Z.Jiang,B.L.Zhou,W.Z.Guo,R.Q.Li,Z.J.Chen,Sequencing of allotetraploid cotton(Gossypium hirsutum L.acc.TM-1)provides a resource for fiber improvement,Nat.Biotechnol.33(2015)531-537.
[39]K.B.Wang,Z.W.Wang,F.G.Li,W.W.Ye,J.Y.Wang,G.L.Song,Z.Yue,L.Cong,H.H.Shang,S.L.Zhu,C.S.Zou,Q.Li,Y.L.Yuan,C.R.Lu,H.L.Wei,C.Y.Gou,Z.Q.Zheng,Y.Yin,X.Y.Zhang,K.Liu,B.Wang,C.Song,N.Shi,R.J.Kohel,R.G.Percy,J.Z.Yu,Y.X.Zhu,J.Wang,S.X.Yu,The draft genome of a diploid cotton Gossypium raimondii,Nat.Genet.44(2012)1098-1103.
[40]F.G.Li,G.Y.Fan,K.B.Wang,F.M.Sun,Y.L.Yuan,G.L.Song,Q.Li,Z.Y.Ma,C.R.Lu,C.S.Zou,W.B.Chen,X.M.Liang,H.H.Shang,W.Q.Liu,C.C.Shi,G.H.Xiao,C.Y.Gou,W.W.Ye,X.Xu,X.Y.Zhang,H.L.Wei,Z.F.Li,G.Y.Zhang,J.Y.Wang,K.Liu,R.J.Kohel,R.G.Percy,J.Z.Yu,Y.X.Zhu,J.Wang,S.X.Yu,Genome sequence of the cultivated cotton Gossypium arboreum,Nat.Genet.46(2014)567-572.
[41]S.F.Altschul,W.Gish,W.Miller,E.W.Myers,D.J.Lipman,Basic local alignment search tool,J.Mol.Biol.215(1990)403-410.
[42]C.Addo-Quaye,T.W.Eshoo,D.P.Bartel,M.J.Axtell,Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome,Curr.Biol.18(2008)758-762.
[43]S.Griffiths-Jones,H.K.Saini,S.van Dongen,A.J.Enright,miRBase:tools for microRNA genomics,Nucleic Acids Res.36(2008)D154-D158.
[44]O.Voinnet,Origin,biogenesis,and activity of plant microRNAs,Cell 136(2009)669-687.
[45]N.Liu,L.L.Tu,W.X.Tang,W.H.Gao,K.Lindsey,X.L.Zhang,Small RNA and degradome profiling reveals a role for miRNAs and their targets in the developing fibers of Gossypium barbadense,Plant J.80(2014)331-344.
[46]S.J.Szalma,B.M.Hostert,J.R.LeDeaux,C.W.Stuber,J.B.Holland,QTL mapping with near-isogenic lines in maize,Theor.Appl.Genet.114(2007)1211-1228.
[47]R.Kooke,E.Wijnker,J.J.Keurentjes,Backcross populations and near isogenic lines,Methods Mol.Biol.871(2012)3-16.
[48]L.Gong,A.Kakrana,S.Arikit,B.C.Meyers,J.F.Wendel,Composition and expression of conserved microRNA genes in diploid cotton(Gossypium)species,Genome Biol.Evol.5(2013)2449-2459.
[49]X.Y.Guan,Q.X.Song,Z.J.Chen,Polyploidy and small RNAregulation of cotton fiber development,Trends Plant Sci.19(2014)516-528.
[50]G.Baum,S.Lev-Yadun,Y.Fridmann,T.Arazi,H.Katsnelson,M.Zik,H.Fromm,Calmodulin binding to glutamate decarboxylase is required for regulation of glutamate and GABA metabolism and normal development in plants,EMBO J.15(1996)2988-2996.
[51]M.L.Preuss,D.B.Delmer,The cotton kinesin-like calmodulinbinding protein associates with cortical microtubules in cotton fibers,Plant Physiol.132(2003)154-160.
[52]C.M.Shan,X.X.Shangguan,B.Zhao,X.F.Zhang,L.M.Chao,C.Q.Yang,L.J.Wang,H.Y.Zhu,Y.D.Zeng,W.Z.Guo,B.L.Zhou,G.J.Hu,X.Y.Guan,Z.J.Chen,J.F.Wendel,T.Z.Zhang,X.Y.Chen,Control of cotton fibre elongation by a homeodomain transcription factor GhHOX3,Nat.Commun.5(2014)5519.
[53]A.L.Silverstone,H.S.Jung,A.Dill,H.Kawaide,Y.Kamiya,T.P.Sun,Repressing a repressor:gibberellin-induced rapid reduction of the RGA protein in Arabidopsis,Plant Cell 13(2001)1555-1565.
[54]G.Hagen,T.Guilfoyle,Auxin-responsive gene expression:genes,promoters and regulatory factors,Plant Mol.Biol.49(2002)373-385.
[55]E.W.Wang,B.M.Kessler,A.Borodovsky,B.F.Cravatt,M.Bogyo,H.L.Ploegh,R.Glas,Integration of the ubiquitin-proteasome pathway with a cytosolic oligopeptidase activity,Proc.Natl.Acad.Sci.U.S.A.97(2000)9990-9995.