2022 年 3 期目录

  目的  木棉是我国南方地区重要的观赏植物,研究不同时期花蕾和花朵的花瓣基因表达差异,揭示花色变异的遗传调控机制,为建立花色定向育种技术提供科学依据。  方法  以木棉不同发育时期深红色花和黄色花的花瓣为研究对象,利用Illumina HiSeqTM 4000开展转录组测序;分别采用DESeq2和Kyoto Encyclopedia of Genes and Genomes（KEGG）数据库进行差异表达基因鉴定和通路富集分析。  结果  测序共获得75 190个单基因,4 772个单基因能被公共数据库注释;基因功能富集分析显示,基因主要富集在基本功能预测、信号转导机制和转录后修饰、蛋白折叠、分子伴侣等通路。共获得不同时期差异表达基因10 397个,显著富集在29个生物学通路,主要包括光合作用、代谢和植物激素信号转导通路等。其中,参与苯丙烷生物合成通路的差异表达基因有72个,参与类胡萝卜素、黄酮类生物合成通路和苯丙氨酸代谢通路的差异表达基因分别为25个,这些通路和基因均与花青素的生物合成相关;另有4个差异表达基因显著富集在甜菜碱的生物合成通路中。qRT-PCR数据验证了转录组数据的可靠性。  结论  （1）光合作用、代谢、植物激素信号转导、黄酮类生物合成、苯丙烷生物合成、苯丙氨酸代谢及能量代谢等通路相关基因可能参与花瓣的发育过程。（2）黄酮类生物合成通路相关基因在深红色花发育的花蕾中期与花朵期显著高表达,可能是花色呈现红色的主要原因。（3）类胡萝卜素生物合成通路的关键合成酶基因的部分家族成员在黄色花发育的花蕾期和花蕾中期显著高表达,可能是导致花色呈现黄色的主要原因。      

大白菜花发育不同时期的转录组研究

[目的]研究大白菜花发育不同时期的基因表达变化,为进一步完善成花调控分子机理提供依据。[方法]利用Illumina HiSeq4000测序平台对大白菜花芽分化1级和5级的茎尖生长点进行转录组测序,比较花发育过程基因表达的差异。[结果]发现2 859个差异表达基因,且差异表达基因显著富集在苯丙素生物合成、苯丙氨酸代谢、黄酮类化合物生物合成、脂肪酸延长和植物昼夜节律这5条代谢通路上;研究表达差异倍数在10倍以上的基因功能,筛选出14个可能影响大白菜花发育的新基因。[结论]大白菜花发育过程苯丙素生物合成和苯丙氨酸代谢活动有所减弱,而黄酮类化合物生物合成、脂肪酸延长和植物昼夜节律活动有所增强;筛选出的14个基因可能在大白菜花发育过程有重要作用,可进一步研究。     

藏医临床常用绿绒蒿药材的鉴别研究简

鉴别藏医临床常用的3种“欧贝”类绿绒蒿药材.采用原植物和药材性状比较,光镜(OM)法观察根、花葶、叶、花的内部组织结构显微特征,扫描电镜(SEM)法观察花粉粒形态以及薄层层析法(TLC)对3种绿绒蒿进行鉴别.结果表明:1)原植物和性状特征可根据花颜色不同或植株大小区分;2)显微组织结构特征可根据根的木质化程度和叶的主脉维管束差异进行鉴别;3)粉末特征以花粉囊内壁细胞的特征最为明显;4)花粉粒表面特征:全缘叶绿绒蒿和红花绿绒蒿的花粉粒为球形或近球形,无萌发孔,全缘叶绿绒蒿花粉粒上的刺状突起较红花绿绒蒿钝,红花绿绒蒿花粉粒表面比全缘叶绿绒蒿平滑;川西绿绒蒿花粉粒类型为3孔沟,表面具刺状雕纹.得出3种绿绒蒿的原植物、性状、内部显微组织结构、粉末特征、花粉粒表面特征等区别明显,可为藏医临床常用的3种“欧贝”类绿绒蒿药材的鉴别提供参考依据.     

藏医临床常用绿绒蒿药材的鉴别研究

鉴别藏医临床常用的3种"欧贝"类绿绒蒿药材.采用原植物和药材性状比较,光镜(OM)法观察根、花葶、叶、花的内部组织结构显微特征,扫描电镜(SEM)法观察花粉粒形态以及薄层层析法(TLC)对3种绿绒蒿进行鉴别.结果表明:1)原植物和性状特征可根据花颜色不同或植株大小区分;2)显微组织结构特征可根据根的木质化程度和叶的主脉维管束差异进行鉴别;3)粉末特征以花粉囊内壁细胞的特征最为明显;4)花粉粒表面特征:全缘叶绿绒蒿和红花绿绒蒿的花粉粒为球形或近球形,无萌发孔,全缘叶绿绒蒿花粉粒上的刺状突起较红花绿绒蒿钝,红花绿绒蒿花粉粒表面比全缘叶绿绒蒿平滑;川西绿绒蒿花粉粒类型为3孔沟,表面具刺状雕纹.得出3种绿绒蒿的原植物、性状、内部显微组织结构、粉末特征、花粉粒表面特征等区别明显,可为藏医临床常用的3种"欧贝"类绿绒蒿药材的鉴别提供参考依据.     

桃嫩枝扦插不定根发育的分子机制

利用高通量测序技术对桃砧木‘GF677’生根发育的4个时期进行转录组测序和数据分析，为阐明桃不定根发育的分子机制奠定理论基础。以扦插后0(对照)、7(激活期)、14(愈伤形成期)、21(不定根形成期)、28 d(不定根伸长期)插穗基部的韧皮部为材料进行RNA-sep测序并研究其生根机理。结果表明，生根过程中不同时期共鉴定出25 656个差异表达基因，其中上调13 166个，下调12 490个。GO功能注释表明，差异表达基因主要富集在代谢过程、细胞过程、细胞、细胞组分、细胞器黏合物和催化活性等代谢途径中。KEGG富集结果显示，差异表达基因主要响应糖酵解、半胱氨酸和蛋氨酸、精氨酸和脯氨酸、核糖体和氨基酸生物合成。应用实时荧光定量qRT-PCR对主要通路的15个DEGs表达模式进行验证，其中13个基因表达水平的变化与转录组基因丰度的变化基本一致，表明转录组数据具有较高的可靠性。总之，糖酵解通路关键基因显著上调，参与激活期的能量补给；蛋氨酸代谢通路诱导基因的表达，参与韧皮部维管组织的分化与形成；精氨酸代谢激活NOS基因的表达，有利于NO的传递，实现根源信号向地上部的运输，进而调控桃的不定根发育...     

花魔芋响应软腐病菌侵染初期的转录组分析

为探讨花魔芋抵御软腐病害的分子机制,分别以染病初期和健康的花魔芋球茎为材料进行转录组测序。结果表明,健康和染病初期花魔芋两组出现3222个差异表达基因,其中2660个基因上调表达,562个基因下调表达,差异表达基因主要涉及细胞组分、生物学过程及分子功能等生理生化过程。表达量上调和下调最大的前10个基因包含编码MiAMP1抗菌蛋白、热激蛋白、胰蛋白酶与蛋白酶抑制剂、胰凝乳蛋白酶抑制剂等与抗病相关的基因。GO功能分类和KEGG通路分析表明,类黄酮和苯丙烷类物质在花魔芋抵御软腐病菌侵染初期发挥重要作用。值得关注的是,硒化合物代谢、维生素B6代谢、类黄酮生物合成、N-聚糖生物合成及链霉素生物合成通路等抗病相关的差异表达基因在花魔芋染病初期全部或大部分受诱导表达。利用高通量测序获得大量花魔芋转录组信息,有助于挖掘与花魔芋抗软腐病相关的关键基因,也为魔芋抗病分子育种提供理论参考。     

基于转录组测序的铁皮石斛黄酮代谢途径及相关基因解析

【目的】分析铁皮石斛Dendrobium officinale黄酮类化合物的生物合成途径及相关基因,为铁皮石斛黄酮类化合物代谢调控、药用价值的开发研究提供参考。【方法】利用Illumina HiSeq 4000测序平台对铁皮石斛2个生长阶段茎、叶进行高通量转录组测序,对组装获得的unigenes进行功能注释和黄酮类化合物的生物合成相关基因解析。【结果】铁皮石斛黄酮类化合物代谢相关Unigenes 48个,涉及14个酶。5个CHS相关Unigenes具有CHSlike保守结构域,活性位点氨基酸残基(Cys-His-Asn)高度保守,丙二酰辅酶A结合位点和产物结合位点的部分氨基酸残基发生变异。CHI(Unigene0013781)属于类型I查尔酮异构酶,异构化6′-羟基查尔酮为5-羟基黄烷酮。表达分析表明,2个生长时期的茎和叶,CHS(Unigene0008250)、 CHI(Unigene0013781)和F3H的表达量都高于CHS(Unigene0012884)和C3′H。【结论】通过对转录组数据分析,共获得铁皮石斛黄酮类化合物代谢相关Unigenes48个,涉及14个酶;5个CHS相关Unigenes中,2个编码查尔酮合酶,3个编码联苄合酶;CHI(Unigene0013781)属于类型I查尔酮异构酶。     

不同发育时期核桃内果皮的差异代谢物分析

为探明核桃内果皮发育过程中代谢物的变化和木质化发育规律,以不同发育期的‘赞美’核桃内果皮为材料进行了代谢物测定,同时结合转录组信息,研究不同发育时期核桃内果皮中代谢物质的变化及相关基因的转录与表达。结果表明:不同发育期核桃内果皮中代谢物具有显著差异。通过高效液相色谱-质谱联用(UPLC-MS)对样品进行检测分析共鉴定出7 837个代谢物,对其进行筛选后,共得到2 633个差异代谢物,包括核苷酸及其衍生物、脂质、酚酸、氨基酸及其衍生物、生物碱、糖苷、有机酸、木脂素、黄酮类等,被注释到35个代谢通路中,其中玉米素生物合成、组氨酸代谢、核黄素代谢、硫胺素代谢、苯丙氨酸和酪氨酸及色氨酸生物合成通路最显著;硫代葡萄糖苷生物合成、嘌呤代谢、卟啉和叶绿素代谢和苯丙烷类生物合成通路富集到最多差异代谢物。花后45 d时,大部分的差异代谢物的含量较花后20 d时降低。转录组与代谢组联合分析发现,二者具有22个相同的代谢通路,包括玉米素生物合成通路等;激素代谢相关通路和苯丙烷生物合成相关通路中差异代谢物变化与酶编码基因的表达情况一致。核桃内果皮代谢物分析为揭示内果皮发育规律,提高其利用价值提供了理论依据。     

茶树花不同发育时期的转录组分析

为了揭示茶树花不同发育时期生长发育和生化成分代谢的分子机制,以‘黄棪’茶树品种为研究对象,通过高通量测序技术对茶树花两个发育时期(花苞期、开放期)的转录组进行比较分析.结果显示:转录组测序共获得135 039 271个过滤序列,GC含量为44.30%～45.35%;茶树花两个发育时期共获得18 352个差异表达基因,其中,上调的基因有8 019个,下调的基因有10 333个.对差异表达基因进行基因本体论(GO)及京都基因和基因组百科全书(KEGG)分析,发现其主要富集在与茶树花生长发育、生化成分形成有关的植物激素信号转导、光合作用—天线蛋白、植物的昼夜节律、糖胺聚糖降解等代谢途径上.本研究结果表明,茶树花的发育受到植物激素的调控和光周期的诱导,与氨基酸、可溶性总糖有关的代谢途径随着花的发育发生了变化,该结果为进一步研究茶树花发育过程中的动态变化提供了参考.     

藤茶高通量转录组分析及黄酮类化合物合成相关基因挖掘

[目的]分析藤茶高通量转录组序列,从中挖掘出黄酮类化合物合成相关基因,为进一步揭示藤茶黄酮类化合物生物合成调控机制提供理论参考.[方法]分别采集藤茶的幼叶和成熟叶,提取其总RNA构建cDNA文库,采用Il-lumina HiSeqTM 4000高通量测序平台对藤茶叶片进行转录组测序,经过滤处理后运用Trinity组装,将获得的Unigene与Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG 7个数据库进行比对注释,并预测Unigenes的编码区序列(CDS);基于KEGG信号通路富集分析,发掘藤茶黄酮类化合物合成相关基因.[结果]藤茶叶片转录组测序获得82126236条原始测序序列(Raw reads),过滤处理后得到80156972条高质量序列(Clean reads),进一步组装拼接得到92472条Unige-nes,平均长度为1208 bp,N50长度为1780 bp,其中,至少在1个数据库注释的Unigenes有84217条,占Unigenes总数的91.07%,有8944条Unigenes在7个数据库均被注释,占Unigenes总数的9.67%.在GO数据库成功注释的41116条Unige-nes可分为生物学过程、细胞组分和分子功能三大类,共56个小类;在KOG数据库注释的14553条Unigenes可分成25类,其中,一般功能预测注释成功的Unigenes最多(1946条);其次是翻译后修饰、蛋白质翻转、分子伴侣(1776条),参与次生代谢物质的生物合成、转运和降解的Unigenes较少,仅有319条;KEGG信号通路富集分析发现,共有15262条Unigenes注释到128条KEGG信号通路,以注释为代谢的Unigenes最多,为8694条,其中筛选获得有98个黄酮类化合物合成相关基因,分别编码苯丙烷代谢通路的3种关键酶和类黄酮代谢通路的14种关键酶.藤茶叶片转录组Unigenes与Swiss-Prot和Nr数据库比对,获得52582条CDS序列,ESTScan 3.0.3预测获得35535条CDS序列.[结论]藤茶在细胞过程、代谢过程、单有机体过程、细胞和细胞部分、结合和催化活性能力分布的基因较丰富,在一般功能、翻译、翻译后修饰、蛋白质翻转及分子伴侣的基因表达量较高,具有较强的碳水化合物代谢能力.多种关键酶基因参与藤茶黄酮类化合物的生物合成,推测其生物合成途径存在多条分支,调控机制也较复杂.     

桑葚幼果的落果与正常果的果柄转录组分析

  目的  果实的脱落是一种常见的生理现象，会在一定程度上限制果实的产量。探究桑葚幼果在脱落过程中的代谢及生化反应，可为桑葚的脱落机制研究提供参考。  方法  通过高通量测序对白桑Morus alba幼果在落果和正常果的果柄离区进行转录组测序，筛选差异表达基因，并对其进行功能注释和代谢通路分析，采用实时荧光定量PCR技术对转录组数据进行验证。  结果  经华大基因转录组测序分析共计获得262.92 Mb的原始序列；桑葚幼果果柄在2种状态下共筛选到10 481个差异表达基因，其中，有5 239个差异表达基因上调表达，5 242个差异表达基因下调表达。经基因本体数据库(GO)功能富集分析，共发现37个显著性条目，大多数差异基因具有催化活性、膜的组成成分、氧化还原酶活性、膜的内在成分等功能；京都基因与基金组百科全书(KEGG)富集分析发现:大多数差异表达基因集中在柠檬酸循环、植物激素信号转导途径、黄酮类生物合成等代谢通路上。  结论  桑葚在脱落过程中受植物激素的调控，氨基酸、黄酮类次生代谢物和碳水化合物等物质也能调控果实脱落。图9表2参35     

大白菜倍性变异引起花发育变化的microRNA调控机制研究

为探索microRNAs(miRNAs)在大白菜花发育中的调控机制,以大白菜DH系FT为试材,通过游离小孢子培养创制出同一基因组的二倍体和四倍体纯系,并稳定遗传。与二倍体相比,四倍体植株在形态上表现出生长发育缓慢、开花时间延迟、花器官整体变大、花瓣颜色更深等特征。利用sRNA-seq技术,分别构建了上述二倍体和四倍体的small RNA(sRNA)文库。通过与植物数据库miRBase比对,分别获得92个和95个成熟的bra-miRNAs,新预测的分别有70个和74个novel miRNAs。通过差异表达分析,筛选获得了34个差异表达的miRNAs,并预测了其287个靶基因。根据基因功能注释分析,筛选出9个与花发育相关的miRNAs,其靶基因涉及类黄酮合成、激素调控、细胞生长、开花时间调控等多方面。通过对287个靶基因进行GO功能富集和KEGG代谢通路分析,其中179个靶基因GO功能富集到810个条目上。最显著差异的条目有:核酸结合、细胞内组分、细胞器、细胞内的细胞器及DNA结合。KEGG代谢通路分析表明,有126个靶基因参与到39条代谢通路上。富集基因显著且数量最多的通路是代谢途径,其次是氰基氨基酸代谢、木质素生物合成以及淀粉和蔗糖代谢。通过对涉及影响四倍体植株花发育变化的靶基因功能进行分析,讨论其9个miRNAs在花发育过程中可能参与的调控作用,推测出有6个miRNAs(bra-miR156a-5p、bra-miR172b-5p、brami R172c-3p、bra-miR172c-5p、bra-miR172d-5p和novel_87)在花期通过互相调控多个靶基因导致开花时间的延迟;有2个miRNAs(bra-miR395d-3p、novel_122)参与调控生长素和脱落酸的变化,可能导致花瓣变大及生长发育变缓等表型特征;brami R395d-3p调控的靶基因还涉及类胡萝卜素生物合成,同时另一个miRNA (bra-miR9552a-3p)可能参与类黄酮生物合成的调控,它们的调控可能是使花瓣颜色变深的主要原因。本研究为揭示由大白菜倍性变异引起的花发育表型变化的分子机制奠定了基础。     

刺瓣绿绒蒿一新变型

[目的]介绍刺瓣绿绒蒿的一个新变型。[方法]在青海省花石峡采集绿绒蒿属藏药时发现了刺瓣绿绒蒿(Meconopsis horridulavar.spinulifera.)的一个白花新变型,与原变种(总状绿绒蒿Meconopsis horridulavar.racemosa)进行了比较。[结果]在青海省玉树州玛多县花石峡生境为海拔4 200 m的小河岸边阳坡草地上发现了刺瓣绿绒蒿的白花新变型,与原变种的区别在于花冠白色;花瓣两面中下部疏生细刺;花柱具4棱,棱呈膜质翅状,宽约1.5 mm;花丝窄线形。[结论]发现的刺瓣绿绒蒿新变型与原变型的主要区别是植物花冠呈白色。     

盐胁迫下园林植物彩叶树响应菌根共生的比较转录组分析

土壤盐分是典型的非生物胁迫因素之一,严重影响着植物的生长发育。为了解盐胁迫下鸡爪槭(Acer palmatum)对接种丛枝菌根真菌的应答分子机制,对接种丛枝菌根处理(AM)、盐胁迫处理(SS)、盐胁迫接种丛枝菌根处理(AS)及对照处理(CK)进行转录组测序分析。结果表明,基于转录组测序技术(RNA-Seq)测序共获得4 672个新基因;将SS与AS处理进行比较时,鉴定出455个差异表达基因(DEGs),其中286个基因上调表达,169个基因下调表达。对获得的DEGs进行功能注释及富集分析,GO(gene ontology)分析结果表明,盐胁迫下接种丛枝菌根真菌涉及蛋白质生物合成相关过程、ATP生物合成过程、蛋白质谷胱甘肽化、细胞分化调控、氮同化相关过程、呼吸电子传递链以及类胡萝卜素代谢等生物过程。KEGG(kyoto encyclopedia of genes and genomes)富集分析结果表明,主要涉及苯丙烷、黄酮类、二芳基庚烷类及姜酚的次生代谢通路。综合GO、KEGG富集分析结果可知,这些基因主要参与植物细胞内部环境的改善、氮代谢相关过程和宿主光保护机制。实时荧光定量PCR(...     

西方蜜蜂工蜂不同虫态发育的转录组学分析

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（DEGs）进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因（51.86%上调,48.14%下调）,在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因（57.51%上调,42.49%下调）,在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因（52.95%上调,47.05%下调）。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。     

多穗柯MYB转录因子的鉴定及适应光因子调控黄酮类合成成员的筛选

为了研究多穗柯(Lithocarpus polystachyus Rehd)MYB转录因子在光胁迫时调控多穗柯黄酮类物质合成的作用，以不同光处理条件的多穗柯为试验材料，通过转录组学测序和生物信息学分析方法，对多穗柯MYB转录因子和多穗柯黄酮类物质合成关键基因进行了分析。结果标明：在多穗柯的转录组学测序数据中，鉴定得到60个LpMYB转录因子，按照DNA结合结构域的特点，将这些转录因子分为4个亚族(1R-MYB,R2R3-MYB,3R-MYB,4R-MYB),其中R2R3-MYB数量最多，占总LpMYB转录因子的43.33%。多穗柯大部分MYB转录因子能够适应红光胁迫；长时间的光照、适当的光照度提高了LpMYB转录因子的表达。鉴定得到10个参与调控黄酮类物质生物合成的S4、S5、S6和S7亚族LpMYB转录因子。多穗柯黄酮类物质合成相关基因与LpMYB转录因子的相关性分析结果表明，LpMYB转录因子在调控黄酮化合物的合成中发挥着重要的作用，LpMYB转录因子能够正调控黄酮类物质合成基因PAL、C4H、CHS3、ANS的表达，从而提高多穗柯的品质。     

基于RNA-seq技术对不同品种猪背最长肌差异表达基因的筛选与注释

【目的】筛选皖南花猪和大约克猪背最长肌组织差异表达基因。【方法】采集皖南花猪和大约克猪背最长肌(各3头),测定其肉质性状指标,并提取其RNA,采用高通量转录组测序(RNA-seq)技术进行测序,对所获序列进行GO和Pathway显著性富集分析。【结果】测序后皖南花猪3个样本得到的总序列数分别为65 584 126,64 327 738和59 244 502,其中有效序列达到94.4%,94.3%和94.2%;大约克猪3个样本得到的总序列数分别为57 969 134,68 254 168和72 298 142,其中有效序列达到93.8%,93.7%和93.8%。测序饱和度良好,证明测序数据真实可靠。差异表达分析结果显示,共筛选到347个差异表达基因,其中包含94个上调基因,253个下调基因。GO功能显著性富集和Pathway显著性富集分析发现,差异表达基因富集在与糖酵解代谢和肌肉发育相关的生物学过程中以及与脂肪酸代谢和生长性状、胴体性状相关的信号通路上。【结论】获得了猪背最长肌组织基因表达谱,筛选出了一些与肉质性状和生长性状有关的候选基因。