三种绵羊BMP15基因第二外显子多态性及与产羔数的相关性分析

绵羊存在影响多胎性状的不同主效基因,选择影响Romney Hanna绵羊和Cambridge绵羊高繁殖力的骨形态发生蛋白15 (bone morphogenetic protein 15, BMP15)为候选基因,采用PCR-SSCP的方法检测BMP15基因外显子Ⅱ第747位点(T747→C)和755位点(T755→C)在蒙古羊、甘肃高山细毛羊、小尾寒羊三种绵羊母羊中的多态性,同时还研究了上述两处突变对三种绵羊产羔数的影响。表明:(1)一共检测到野生纯合型AA、突变杂合型AB (T747→C)、AC (T755→C)三种不同的基因型,AA为优势基因型,A为优势等位基因;(2)三种基因型在甘肃高山细毛羊中均被检测到,而蒙古羊和小尾寒羊中未检测出AB基因型;(3)突变杂合型蒙古羊(AC)比野生纯合型(AA)的平均产羔数多0.27只(p<0.05)。(4)AC的基因型频率,双羔母羊和多羔母羊均高于单羔母羊。根据以上实验推测,BMP15第755位点发生的T→C突变(AC型)对蒙古羊一胎产双羔影响十分显著,甘肃高山细毛羊中AC基因型的绵羊其产羔数有比AA基因型和AB基因型多的趋势,因此该位点可能是一个影响绵羊高繁殖力潜在的DNA标记。     

四个绵羊品种BMPR-IB基因CDS区864位点多态性及其产羔数相关性分析

为了探索绵羊产羔性状与绵羊BMPR-IB基因多态性位点关系,找寻调控绵羊繁殖力的分子标记,以小尾寒羊(多胎母羊)、湖羊(多胎母羊)、蒙古羊(单,双胎母羊)、和甘肃高山细毛羊(单,双胎母羊)样本为试验材料,运用DNA直接测序及PCR-SSCP技术的方法,对BMPR-IB基因CDS区864位点进行分析。试验羊品种出现3种基因型AA、AB、BB,该基因编码区第846位点发生TC的突变。湖羊母羊以及小尾寒羊母羊的优势基因型均为AA型,而蒙古羊母羊和甘肃高山细毛羊母羊AB型基因频率略高于AA型。4种绵羊的优势等位基因均为A。小尾寒羊、甘肃高山细毛羊、湖羊三种绵羊x2值和G2值均未达到显著水平(p0.05),而蒙古羊的x2值和G2值达到显著水平(p0.05),说明除了蒙古羊其他3种羊均达到Hardy-weinberg平衡状态。四种绵羊的观测值F在BMPR-IB基因中均处于95%置信区间内,且接近于上线。根据各品种间的PIC多态信息含量可知,属于低度多态的是小尾寒羊和湖羊(PIC0.25),而在蒙古羊和甘肃高山细毛羊信息含量处于中度多态(0.25PIC0.5)。显著性检验分析表明:甘肃高山细毛羊与小尾寒羊、湖羊,蒙古羊与小尾寒羊、湖羊中差异呈现显著(p0.05)。BMPR-IB基因编码区第846位点突变在不同繁殖力母羊群体中分布不同,该位点可以作为绵羊潜在分子标记位点。     

河西绒山羊DRB3基因外显子2多态性分析

研究采用PCR-SSCP方法检测1 046只河西绒山羊DRB3基因第2外显子的多态性,并克隆、测序群体内变异产生的各等位基因,分析各等位基因间的遗传关系和进化意义。结果表明,河西绒山羊DRB3基因第2外显子表现了18个等位基因(LG001-LG018)控制的31个基因型,LG005为优势等位基因,LG005*005为优势基因型,18个单倍型序列分析发现44个核苷酸多态位点、29个氨基酸多态位点;与GenBank下载序列比对分析结果表明,16个DRB3的等位基因是本研究首次发现;DRB3基因遗传多态指数中,观察杂合度为0.269 6,多态信息含量为0.869 9,有效等位基因数为1.369 1;经χ2适合性检验,该群体偏离了Hardy-Weinberg平衡状态(P<0.05);18个DRB3基因外显子2的单倍型序列NJ系统发育树呈2支分化趋势。研究认为河西绒山羊DRB3基因第2外显子具有丰富的遗传多态性;河西绒山羊DRB3基因最初由2个等位基因突变分化成两大类等位基因。     

三个绵羊群体MC4R基因多态性及生物信息学分析

旨在探讨绵羊黑素皮质素受体-4(melanocortin-4 receptor,MC4R)的分子机理,采用PCR-SSCP方法对3个绵羊群体(甘肃肉用绵羊新品种群、小尾寒羊和湖羊)的MC4R基因外显子进行多态性检测和生物信息学分析。结果表明,3个绵羊群体均存在3种基因型AA型、AB型和BB型,优势基因型为BB,其中优势等位基因为B;测序结果表明,野生型BB型和突变型AB型相比,AB型个体在该基因编码区第511位点发生G→A突变,第495位发生C→T突变;AA型个体在该基因编码区第511位点发生G→A突变,出现AA的纯合,第495位发生C→T突变,出现CC纯合;3个绵羊群体中小尾寒羊的多态信息含量属于中度多态(0.25PIC0.50),甘肃肉用绵羊新品种群羊和湖羊属于低度多态(PIC0.25);χ2适合性检验表明除湖羊之外,其余2个绵羊品种均处于Hardy-Weinberg平衡状态。生物信息学分析发现MC4R氨基酸序列有明显的疏水性区域,有7个跨膜螺旋区及信号肽,其编码蛋白主要的二级结构元件是α螺旋和无规则卷曲;同源性比对发现绵羊MC4R基因与山羊、牛、野猪、人类及大猩猩的相似度分别为97%、94%、81%、83%及83%,说明MC4R是一个非常保守的蛋白,在绵羊的生长发育中起着重要作用。     

甘肃西门塔尔牛甘州类群MSTN基因外显子多态性分析

采用PCR-SSCP方法检测50头甘肃西门塔尔牛甘州类群MSTN基因3个外显子的遗传多态性,并对群体内各等位基因进行测序,旨为探讨西门塔尔牛MSTN基因与肉质性状关联等方面的研究提供理论依据。结果显示,甘肃西门塔尔牛甘州类群MSTN基因3个外显子中,外显子2受B、C等位基因控制,形成BB、CC和BC 3种基因型,基因型频率分别为0.22(11/50)、0.64(32/50)、0.14(7/50);外显子1和外显子3分别受A等位基因与D等位基因控制,形成AA和DD基因型,基因型频率均为1(50/50)。序列分析表明,甘肃西门塔尔牛甘州类群MSTN基因3个外显子中,外显子2在第41 bp处发生了C→T的突变,但并没有导致氨基酸发生改变,属于同义突变;外显子1和3均无突变。统计结果表明,甘肃西门塔尔牛甘州类群MSTN基因3个外显子中,外显子2为中度多态,其余均无变化。     

八眉猪SLA-DRA基因第4外显子多态性分析

[目的]研究八眉猪SLA-DRA基因第4外显子的多态性,确定其等位基因数以及核苷酸和氨基酸变异位点,探讨其遗传特性。[方法]采用PCR-SSCP和克隆测序的方法对八眉猪SLA-DRA基因第4外显子进行多态性检测。[结果]八眉猪SLA-DRA基因第4外显子共有3种等位基因(A、B、C),形成3种基因型(AA、BB、BC),其中B为优势等位基因,BC为优势基因型。序列对比结果表明,八眉猪SLA-DRA基因第4外显子共有3个核苷酸突变位点(c.4167AG、c.4246AG和c.4282CT),其中4 167 bp的突变导致谷氨酰胺(Q)突变为精氨酸(R)。群体遗传学分析表明,八眉猪SLA-DRA基因第4外显子的多态信息含量为0.3953,杂合度为0.7241,且八眉猪在SLA-DRA基因第4外显子上偏离了Hardy-Weinberg平衡(p0.01)。[结论]八眉猪SLA-DRA基因的第4外显子属于中度多态,并呈现出杂合子优势。     

甘肃河西地区西门塔尔杂交类群NGB基因第3外显子的遗传变异研究

旨在对甘肃河西的临泽、甘州、武威、金昌、高台5个地区283头西门塔尔杂交类群NGB基因第3外显子的遗传多态性及变异特征进行系统分析,采用PCR-SSCP方法检测了283头西门塔尔杂交类群NGB基因第3外显子和部分内含子的多态性,且对群体内各等位基因进行了测序。结果显示,5个地区西门塔尔杂交类群共检测出5个等位基因(A、B、C、D、E),表现为5种基因型(AA、AB、AC、AD、AE)。其中甘州、武威、金昌西门塔尔杂交类群NGB基因均只检测到AA、AB 2种基因型,高台西门塔尔杂交类群检测到AA、AE 2种基因型,临泽西门塔尔杂交类群检测到AA、AB、AC、AD 4种基因型。A等位基因和AA基因型的频率在5个群体中最高,为优势基因和优势基因型。对不同SSCP带型的对应片段进行测序分析,共发现6个核苷酸突变位点(75 bp C→T,78 bp C→G,128 bp G→A,214 bp G→A,232 bp C→T,233 bp G→A),其中第75 bp和第78 bp处的突变位点位于内含子区域,其余4处突变位点均位于外显子区域。第214 bp处的核苷酸突变导致甘氨酸(Gly)突变为丝氨酸(Ser),第232 bp处核苷酸突变导致精氨酸(Arg)突变为色氨酸(Trp),第233 bp处核苷酸突变导致精氨酸(Arg)突变为谷氨酰胺(Gln),经χ2检验结果显示,5个地区的西门塔尔杂交类群在此3个突变位点上都处于Hardy-Weinberg平衡状态(P0.05)。群体遗传学分析结果表明,临泽、甘州、武威、金昌、高台西门塔尔杂交类群的多态信息含量(PIC)分别为0.0582、0.0196、0.0196、0.0161、0.0159,均属于低度多态(PIC0.25)。     

Leptin基因多态性与绵羊季节性发情的相关分析

根据leptin基因在GenBank中的已知序列设计两对引物,采用PCR-SSCP技术在常年发情的湖羊和季节性发情的阿勒泰羊群体中进行单核苷酸多态性(SNPs)检测,对筛查到的SNP位点进行基因型与绵羊季节性发情的关联分析。结果表明,与湖羊相比,阿勒泰羊leptin基因第1内含子上有3个连续碱基TTG的插入和C/T碱基突变;第3外显子3上发生G/T碱基突变,编码氨基酸由缬氨酸变成亮氨酸。Leptin外显子2扩增片段上检测到AA、AB、BB三种基因型,BB基因型在阿勒泰羊群体中属于优势基因型;对两个群体进行基因型频率独立性χ2检验,差异极显著(P0.001),说明BB基因型是影响季节性发情的有利基因型。研究结果提示,绵羊品种中Leptin基因序列的差异性可能是造成绵羊季节性发情的原因之一,可作为常年发情绵羊品种选育的辅助标记。     

采用6个微卫星位点研究6个绵羊品种遗传多样性

利用6个微卫星座位对中国乌珠穆沁羊、小尾寒羊、滩羊、昭通绵羊、呼伦贝尔羊和甘肃高山细毛羊6个绵羊品种,共计280个个体进行遗传多样性分析。计算了6个绵羊品种间的多态信息含量、杂合度和遗传距离,并进行了主成分分析和UPGMA聚类。发现了101个等位基因,平均杂合度0.599-0.691,平均多态信息含量0.609-0.680;甘肃高山细毛羊与其他绵羊品种遗传分歧最大,而呼伦贝尔羊与乌珠穆沁羊间的分歧最小。结果表明6个绵羊品种均具有较高的遗传多样性,品种间的遗传关系与其形成历史、分化及地理分布基本一致。     

红鳍东方鲀(Takifugu rubripes)MC4R基因的多态性分析

采用PCR-SSCP(single strand conformation polymorphism)技术和DNA测序方法分析红鳍东方鲀MC4R(Melanocortin-4receptor)基因编码区多态性。在MC4R基因编码区48 nt和264 nt均发生了碱基的转换突变(G→A),两个突变位点分别位于M1和M2引物扩增产物中。引物M1扩增产物SSCP分析得到两种基因型:AA基因型和AB基因型,并且AA基因型和A等位基因频率明显高于AB基因型和B等位基因。引物M2扩增产物也得到两种基因型:CC基因型和CD基因型,CC基因型和C等位基因频率明显高于CD基因型和D等位基因。遗传变异结果分析表明,两个突变位点均属于低度多态性,而且群体遗传杂合度较低,反映了该群体的遗传一致性较高。     

藏绵羊PPARα基因遗传变异特征及组织表达分析

PPARα(Peroxisome proliferator activated receptorsα)基因参与调控线粒体脂肪酸β氧化,在能量代谢中具有重要作用,是动物高原低氧适应的主要候选基因之一。本研究采用PCR-SSCP和测序方法检测2个品种(藏绵羊,湖羊) PPARα基因P1区核苷酸变异及遗传特征。采用实时荧光定量PCR (quantitative real-time PCR, qRT-PCR)检测藏羊与湖羊心脏、肝脏、脾脏、肺脏、肾脏、骨骼肌6个组织的相对表达量。检测表明:在扩增区域发现2个SNPs (c.515-22AG, c.515-10CT)对应的等位基因A、B。基因频率在2个群体间具有统计学意义(p0.01),藏绵羊的优势基因型为AB,湖羊优势基因型为AA。推测AB基因型个体可能更有利于藏绵羊适应高原低氧环境,可作为藏绵羊低氧适应选育的分子标记。PPARα基因在2个绵羊的6个组织中均有不同程度的表达且具有品种间差异性。PPARα基因在藏绵羊心脏的表达量极显著的低于其余5个组织(p0.01),而2个品种间,藏绵羊各组织的表达量均低于湖羊,其中在心脏、肝脏、肺脏3个组织表达量极显著的低于湖羊(p0.01)。本研究可丰富绵羊基因组学内容,为藏绵羊低氧适应研究提供基础数据。     

DGAT1基因K232A突变位点在3个中国荷斯坦母牛群体中的基因型频率分布及其与产奶量的相关性

二脂酰甘油酰基转移酶1(acylCoA:diacylglycerol acyltransferase,DGAT1)基因是影响奶牛产奶性状的主效基因之一,其第八外显子上的第232位的赖氨酸转化为丙氨酸(K232A)的错义突变是影响产奶量、乳脂和乳蛋白含量的因果突变.本研究利用GenBank中已公布的牛DGAT1基因DNA序列,于K232A位点两端设计引物,采用PCR-RFLP方法鉴定来自新疆和江西的3个中国荷斯坦母牛群体共454个个体DGAT1基因K232A突变位点的基因型,使用SAS9.0软件分析供试群体基因型分布及其与产奶量的相关性.研究结果表明:K等位基因为441 bp的PCR扩增片段,A等位基因为完全切开的202bp和239bp两片段.新疆1、新疆2和南昌金牛3个荷斯坦母牛群体在该位点上K等位基因的频率分别为51.9%、42.4%和27.2%;A等位基因的频率分别为48.1%、57.6%和72.8%.DGAT1基因K232A突变位点与三个母牛群体日平均产奶量的相关性分析发现在3个群体中KK型个体产奶量均低于AA型个体,但不同基因型个体间产奶量差异在统计学上没有显著差异.本研究为进一步探讨DGAT1基因在中国奶牛群体中的遗传效应及分子标记辅助选择培育良种奶牛提供参考.     

合作猪SLA-DQA基因外显子2多态性分析

旨在研究合作猪DQA基因外显子2多态性,确定其等位基因数、核苷酸多态位点、氨基酸多态位点及各个等位基因之间的遗传关系,分析其进化意义。选用PCR-SSCP对439只合作猪SLA-DQA基因外显子2的多态性进行检测;测序群体内因变异而产生的各等位基因序列,并分析序列数据。结果显示,在合作猪SLA-DQA外显子2中发现了7个新等位基因,共18个核苷酸多态位点,10个氨基酸多态位点。合作猪SLA-DQA外显子2具有较丰富的多态性,群体内可能蕴藏着更加丰富的遗传资源;合作猪SLA-DQA外显子2基因最初可能由一个等位基因突变分化成一大类基因;合作猪SLA-DQA外显子2序列与各个猪种的SLA-DQA外显子2序列具有较高的同源性,预示着这些猪种的SLA-DQA外显子2基因最早可能来源于其分歧之前的共同祖先原始序列;新发现的7个SLA-DQA外显子2等位基因,可能由遗传关系较近的两个等位基因突变产生。     

奶牛β-Lg基因第1外显子PCR-SSCP多态性与泌乳性能的相关性

采用PCR-SSCP技术并结合测序对233头奶牛β乳球蛋白(β-Lg)基因5’端部分序列和外显子1全部序列进行了多态性研究,分析了该基因与奶牛泌乳性状的相关性。结果表明:β-Lg基因5’端和外显子1共存在2个等位基因3种基因型,BB型为优势基因型,B为优势等位基因。该群体在这一位点上偏离Hardy-Weinberg平衡状态,多态信息含量(PIC)为0.3548。测序结果显示,与普通牛该基因序列(X14710)相比,B等位基因在2073 bp、2202 bp和2206 bp处发生了G→C、C→T和A→G的碱基突变,其中2202 bp处的C→T突变导致第11位氨基酸由苏氨酸变为异亮氨酸,而A等位基因在3个位点上与X14710相同。最小二乘法分析表明,BB型305 d乳蛋白量显著高于AA型和AB型(P<0.05);AB型305 d乳脂量显著高于AA型(P<0.05),BB型与AB型之间差异不显著(P>0.05);等位基因B为高乳蛋白量和乳脂量的优势基因,可作为奶牛选育的分子遗传标记。     

绵羊H-FABP基因单核苷酸多态性的研究

以小尾寒羊(48只)、宁夏滩羊(121只)、滩寒杂交羊F1(23只)、无角陶塞特(48只)、萨福克(24只)5个绵羊群体为实验材料, 利用PCR-SSCP和DNA测序技术对心脏型脂肪酸结合蛋白(H-FABP)基因(GenBank登录号: AY157617)外显子2和内含子2部分序列进行单核苷酸多态性(SNPs)检测及遗传多态性分析。结果表明: (1)引物2的PCR扩增产物中存在981(G/A)、1014(A/C)、1019(T/C) 和1058 (-/G ) 4个SNP位点, 表现为AA、BB、CC、AB、AC、BC、AD、CD和BD 9种基因型, 其中AA为优势基因型。经χ2检验后, 除滩羊和萨福克羊外, 其他群体的基因频率和基因型频率均处于Hardy-Weinberg平衡状态。群体遗传多态性分析表明: 宁夏滩羊、小尾寒羊和无角陶塞特羊3个群体中的多态信息含量(PIC)均处于0.25和0.50之间, 为中度多态, 萨福克羊和滩寒杂交羊F1为低度多态(PIC<0.25), 表明脂肪酸结合蛋白基因在不同绵羊品种中具有单核苷酸多态性, 可以进一步作为候选基因来分析其与肌内脂肪含量性状的关联性。(2)引物4的PCR扩增产物中检测到1个SNP多态位点为2407(T/C), 表现为HH、Hh和 hh 3种基因型, 基因型频率大小为HH＞Hh＞hh, 经χ2检验后, 在滩羊和无角陶塞特羊中均为达到Hardy-Weinberg平衡状态, 其多态信息含量均为低度多态(PIC<0.25), 而在小尾寒羊、滩寒杂交羊F1和萨福克羊均没有多态出现。     

甘肃不同地区绵羊TLR9基因多态性分析

针对由于急、慢性呼吸道引起的生产性能降低甚至死亡,选择与该病相关的基因TLR9为研究对象,分析群体中该基因的遗传多态性及变异特征,为进一步揭示舍饲绵羊的遗传特性和生产利用提供基础数据。采用PCR-SSCP检测427只表型正常和58只患呼吸道疾病的舍饲绵羊TLR9基因的多态性,测序群体内变异的各等位基因数列,并构建系统发育树以明确舍饲绵羊TLR9基因等位基因之间的遗传关系。结果显示,甘肃地区绵羊在TLR9基因中发现了4个等位基因,共7个核苷酸多态位点,这些位点多是由点突变形成,其中转换4个(占57.14%),颠换3个(占42.86%)。甘肃地区绵羊TLR9基因具有较丰富的多态性;TLR9基因多态性与绵羊呼吸道疾病有一定的相关性。     

猪TLR4基因外显子1新等位基因的分离及遗传变异分析

文章采用PCR-SSCP方法对亚洲野猪、3个引进的商业化品种和10个中国地方猪品种共893个个体TLR4基因外显子1的遗传变异进行了检测,旨在系统分析国内外猪种TLR4基因的多态性,为探讨该基因在免疫和防御系统中发挥的作用提供依据。结果,在猪TLR4基因外显子1中分离到新的等位基因,共检测到3个等位基因,6种基因型。其中杜洛克检测到AA、BB、CC、AB、AC、BC基因型,有杜洛克血统的苏太猪中检测到BB、CC、BC基因型,长白猪、约克夏中检测到CC、BC基因型,野猪及所有10个中国地方猪品种TLR4基因外显子1高度保守,只检测到CC基因型,中国地方猪品种和引进品种TLR4基因外显子1多态性存在极显著的差异。3种基因型中CC型与GenBank中的序列一致,BB和AA基因型分别存在G93C同义突变位点和G194A无义突变位点,这2个变异位点与抗逆性和一般抗病力的关系值得进一步深入研究。     

兔成纤维细胞生长因子5(FGF5)基因SNP及其与产毛量的相关分析

通过PCR产物直接测序的方法, 对荥经长毛兔、天府黑兔以及加利福尼亚兔的FGF5基因外显子1和外显子3进行单核苷酸多态性分析。在外显子1的217位(位点A)检测到由TCT三碱基插入引起的移码突变, 在外显子3的59位(位点B)和3位(位点C)分别发生了错义突变由T→C和同义突变由T→C。通过计算发现各位点不同的基因型和等位基因频率在3个兔品种中存在较大的差异, 位点A、B在长毛兔和肉兔中均有各自的优势基因型和等位基因。各位点基因型与产毛量的最小二乘分析表明, 位点A各基因型的个体在产毛量上差异不显著(P＞0.05), 位点B各基因型个体产毛量的差异极显著(P＜0.01), 位点C各基因型个体产毛量的差异显著(P＜0.05)。初步推断FGF5基因可能是影响长毛兔产毛量潜在的主效基因或者与主效基因连锁, 可作为长毛兔产毛性状连锁分析的候选遗传标记。     

BMP-15基因外显子Ⅰ在三种绵羊中的多态性及与繁殖力相关性分析

为研究绵羊繁殖力与绵羊BMP-15 (bone morphogenetic protein 15)基因的多态位点关系,并寻找调控绵羊繁殖力的分子标记,以甘肃高山细毛羊、蒙古羊、小尾寒羊三种绵羊为研究对象,采用PCR-SSCP技术与DNA碱基测序相结合的方法,检测该基因第一外显子67、92两个位点在上述三个不同品种绵羊中是否存在FecX^H (Q23→Ter)和FecX^I (V31→D)突变,同时根据检测结果与其繁殖力做相关性分析。结果表明,在三种绵羊中既未检测出与Inverdale绵羊相同的FecX^I突变,也未检测出与Romney绵羊相同的FecX^H突变,因此推测BMP-15基因中影响Romney与Inverdale绵羊高繁殖力的突变位点对以上三种绵羊均无显著影响。     

DGAT-1基因K232A突变位点对甘肃地区中国荷斯坦牛泌乳性状的影响

为探索二酰甘油酰基转移酶(DGAT-1)基因多态性与奶牛泌乳性状的相关性,以232头甘肃地区中国荷斯坦牛为实验材料,利用PCR-SSCP技术并结合测序研究DGAT-1基因K232A位点遗传多态性,采用混合动物模型分析DGAT-1基因K232A突变位点对305 d产奶量、305 d乳脂量和305 d乳蛋白量的影响。结果表明:DGAT-1基因K232A位点共存在KK、KA和AA三种基因型,频率分别为0.5086、0.3750和0.1164,等位基因K和A的频率分别为0.6961和0.3039,多态信息含量(PIC)为0.3336,实验群体在这一位点上处于Hardy-Weinberg平衡状态。该位点突变对305 d产奶量的影响达到极显著水平(P<0.01),对305 d乳脂量和305 d乳蛋白量的影响达到显著水平(P<0.05)。最小二乘法分析表明,AA和KA型305 d产奶量极显著高于KK型(P<0.01),KK型乳脂量显著高于AA和KA型(P<0.05),AA型乳蛋白量显著高于KK型(P<0.05);A等位基因是提高产奶量和乳蛋白量的优势基因,而K等位基因是高乳脂量的优势基因。DGAT-1基因K232A位点突变对甘肃地区中国荷斯坦牛泌乳性状有较大的遗传效应,可用于其泌乳性状的分子标记辅助选择。