陆地棉叶绿体基因组密码子使用偏性的分析

 【目的】通过分析陆地棉叶绿体基因组密码子使用偏性,探讨影响密码子偏性形成的主要因素。【方法】以陆地棉叶绿体基因组为研究对象,利用Mobyle及其它软件,分析叶绿体50个基因密码子的使用模式。【结果】陆地棉叶绿体基因组密码子GC1与GC2之间的相关系数达到极显著水平,GC3含量最低,与GC12的相关系数为0.14,双尾检验未达到显著水平;单个基因ENC比值多分布在-0.05—0.05;对应性分析,第一轴上显示了10.27%的差异,与GC3、ENC的相关系数分别为0.223和-0.147,均未达到显著水平。【结论】陆地棉叶绿体密码子第三位偏好使用嘧啶,而最优密码子多以A或T结尾;密码子使用模式可能是由选择和突变及其它因素共同作用形成的。     

黄丹木姜子叶绿体基因组特征分析

【目的】分析黄丹木姜子（Litsea elongata）叶绿体基因组特征,为木姜子属物种鉴定、遗传多样性分析和资源保护提供理论参考。【方法】基于Illumina HiSeq 2000高通量测序平台对黄丹木姜子叶绿体基因组进行测序,利用GeSeq在线工具对叶绿体基因组进行注释,并利用REPuter、MISA、CodonW和IQ-TREE等生物信息学软件对其基因组结构、基因数目、序列重复、密码子使用偏性和系统发育进行分析。【结果】黄丹木姜子叶绿体基因组全长为154028 bp,具有典型的四分结构,编码126个基因,其中蛋白编码基因82个,rRNA基因 8个,tRNA基因 36个。叶绿体基因组的注释基因中,有9个基因含1个内含子,有3个基因含有2个内含子,其余基因均不含内含子;44个基因编码蛋白参与光合作用信号途径,21个基因编码蛋白构成了核糖体大小亚基。黄丹木姜子叶绿体基因组含有32对长序列重复和90个SSR位点,其中,正向重复和回文重复最多,均为12对,反向重复和互补重复分别为6和2对;95.56%的SSR位点位于单拷贝区［大单拷贝区（LSC）和小单拷贝区（SSC）］,仅4.44%的SSR位点位于反向重复区（IR）。黄丹木姜子叶绿体蛋白编码基因GC含量为39.14%,GC3s为27.95%,平均有效密码子数（ENC）为49.04,说明其密码子偏性弱;相对同义密码子使用度（RSCU）大于1.00的密码子31个,其中13个以A结尾,16个以U（T）结尾。系统发育进化树分析结果显示,木姜子属的14个物种聚为两组,其中黄丹木姜子和10种木姜子属植物聚在一个组,与日本木姜子的亲缘关系最近。【结论】黄丹木姜子叶绿体基因组结构保守,偏好A或U（T）结尾的密码子,鉴定的SSR位点可用于物种鉴定和群体遗传学研究。     

香蕉基因组密码子使用偏好性分析

[目的]分析香蕉基因组的密码子组成及使用偏好性,探讨影响密码子偏好性形成的主要因素,为提高香蕉外源基因的表达水平及转基因抗病育种提供参考.[方法]以香蕉基因组的36242条高置信蛋白编码基因CDS序列为研究对象,运用CodonW 1.4.4统计分析香蕉基因组的密码子组成及使用参数,确定最优密码子,并分析密码子使用参数间相关性.[结果]从香蕉基因组数据中筛选获得36242个高置信蛋白编码基因CDS序列,平均长度为1035 bp,GC含量为3.0%~75.8%,其中低于20.0%的仅13个序列,全基因组中GC总含量为50.4%;同义密码子第3位出现G或C的频率为52.9%,比出现A或T的频率高.香蕉基因组的有效密码子数(ENC)介于20.0~61.0,平均为50.7;共有17个最优密码子,其中有15个密码子的第3位是G或C;基因编码区的长度和ENC存在正相关,随着基因编码区长度的增加,对以G或C结尾的密码子使用偏好性逐渐降低,且编码区长度为400~600 bp的基因具有较高的基因表达水平.[结论]香蕉基因组中多数基因的密码子使用偏好性较弱,但少部分基因具有强偏好性,偏好使用以G或C结尾的密码子,且偏好性受核苷酸组成、基因突变及自然选择等因素的影响.     

无刺龙舌兰叶绿体基因组特征及密码子偏好性分析

【目的】分析无刺龙舌兰叶绿体基因组特征及密码子偏好性,为无刺龙舌兰叶绿体相关基因的表达、修饰和物种进化研究提供参考。【方法】对无刺龙舌兰的叶绿体基因组进行测序、组装和注释,分析密码子偏好性及其影响因素,并通过建立高、低基因表达库,筛选出最优密码子。基于20个已发表的龙舌兰科植物叶绿体基因组数据构建系统发育进化树。【结果】无刺龙舌兰叶绿体基因组总长157579 bp,大单拷贝区(LSC)、小单拷贝区(SSC)和2个反向重复区(IRa和IRb)的长度分别为85940、18279和26680 bp,GC含量为37.8%,包括135个基因(85个蛋白编码基因、38个tRNA基因、8个rRNA基因及4个未知功能的基因),从中筛选出51个长度大于300 bp的基因编码区(CDS)序列,其有效密码子数(ENC)均大于41.0。GC1、GC2、GC3和GC3s含量分别为46.75%、39.61%、29.19%和26.06%,说明密码子第3位多以A/T结尾。GCall与GC1、GC2和GC3均呈极显著相关(P<0.01,下同),但GC3与GC1和GC2均无显著相关性(P>0.05,下同),表明密码子第1、2位的碱基组成相似,但与第3位的相似度不高。选择和突变是导致叶绿体基因组密码子偏好性的主要因素。筛选出14个多以A/U结尾的最优密码子。无刺龙舌兰与克雷塔罗丝兰和西地格丝兰为姊妹关系,自荐值为100%。【结论】无刺龙舌兰叶绿体基因组为保守的四分体结构,叶绿体基因组密码子偏好性较弱,主要受选择和突变等多因素影响。基于植物叶绿体基因组构建系统发育进化树在物种的分类鉴定及确定各物种间系统发育关系的研究中是一种准确、可靠的方法。     

灯盏花叶绿体基因组密码子偏好性分析

【目的】了解灯盏花[Erigeron breviscapus (Vant.) Hand-Mazz]叶绿体基因组密码子的使用偏好性。【方法】以灯盏花叶绿体基因组全序列为研究对象,利用CodonW 1.4.2软件和CUSP程序等对筛选出的47条CDS (coding DNA sequence)进行分析,并进行中性绘图、ENC-plot绘图和PR2-plot绘图分析。【结果】叶绿体基因组密码子中不同位置的GC碱基含量比值为GC_1GC_2GC_3,说明密码子末位碱基偏性以A/T (U)结尾;有效密码子数的均值为47.18,密码子使用偏性较弱;3种绘图分析结果发现:影响灯盏花叶绿体基因组密码子偏好性的主要因素是自然选择,也受到其他因素(如突变)的影响;最优密码子分析中确定了UUU、UUA和UUG等18个最优密码子,且大多数以A和U结尾,有且仅有1个以G结尾。【结论】灯盏花叶绿体基因组密码子第3位偏好性以A或U碱基结尾,密码子使用模式受选择和其他因素共同影响,其中选择的作用较大。该结果可为后续利用基因工程手段对外源基因密码子改造、提高所转基因在灯盏花叶绿体中的表达提供参考。     

H9N2禽流感病毒全基因组密码子使用偏好性及影响因素分析

【目的】探究H9N2禽流感病毒(AIV)全基因组的密码子使用偏好性及影响因素。【方法】选取2010—2018年国内H9N2 AIV流行毒株的全基因组为研究对象,分析其碱基组成特性、最优密码子、密码子使用偏好性的影响因素以及病毒对宿主密码子使用模式的适应性。【结果】H9N2 AIV的全基因组中AU含量高于GC。大部分最优密码子以A或U结尾,有效密码子数(ENC)平均值为52.86,提示存在密码子使用偏好性但偏好性较低。密码子使用偏好性主要受到突变压力和自然选择的共同作用,其中自然选择(所占比例为61.79%～76.15%)作用大于突变压力(所占比例为23.85%～38.21%)。H9N2 AIV对人Homo sapiens的密码子适应指数平均值为0.739～0.741,提示H9N2AIV禽流感病毒可能已适应人类的密码子使用模式。【结论】本研究为H9N2 AIV的基因进化分析、已有疫苗的密码子优化和新型疫苗(密码子去优化疫苗)研制提供了理论依据。     

竹节参转录组使用密码子偏好性分析

竹节参是我国珍稀濒危中药材,研究其基因密码子使用模式,可为利用基因工程技术实现人参皂苷异源生物合成及竹节参分子育种改良提供理论依据。以竹节参转录组测序结果为数据来源,筛选编码蛋白基因序列(codingsequence,简称CDS)碱基数不小于300 bp的11 199条完整开放阅读框序列作为研究对象,利用Codon和SPSS软件分别统计竹节参基因密码子GC含量、密码子第3位的(C+G)含量(GC3)和密码子第1、第2位(G+C)含量的平均值(GC12)、同义密码子的相对使用度(RSCU)、有效密码子数(ENC)等密码子偏好性指标,通过中性绘图(GC12 vs.GC3)、PR2绘图和ENC-GC3s绘图分析影响竹节参密码子使用模式的因素。结果表明,竹节参基因的平均GC、GC12和GC3s含量分别为44. 67%、46. 97%和39. 80%,其密码子使用模式受到突变和选择等多重因素的影响,确定了31个竹节参最优密码子,除了UUG外,其余最优密码子均以A或T结尾。竹节参密码子使用模式与大肠杆菌和酿酒酵母相比差异较大,选取毕赤酵母作为竹节参基因的异源表达宿主更为合适。     

灵芝属真菌线粒体基因组特征及进化

【目的】分析灵芝属（Ganoderma）真菌的线粒体基因组特征及进化,为灵芝属物种分类、分子进化和系统发育分析提供理论依据。【方法】基于灵芝属真菌15个线粒体基因组序列,利用MEGA X、MISA、mVISTA、MAFFT、DnaSP、PAML X和IQ-TREE等生物信息学软件对基因组特征、序列多态性、简单重复序列（SSR）、基因进化和系统发育进行分析。【结果】灵芝属真菌线粒体基因组全长为50603~124588 bp,GC含量为25.4%~27.3%,含有15个保守的蛋白编码基因（PCG）、2个rRNA基因和25~29个tRNA基因。SSR主要由AT构成,单核苷酸重复类型比例最高,其次为三核苷酸重复和四核苷酸重复。种间线粒体基因组序列差异较大,非编码区的变异水平高于编码区,nad6、nad3和cob基因编码序列的变异度较高,内含子长度与线粒体基因组大小呈显著正相关。15个保守的线粒体蛋白编码基因主要受纯化选择影响,其中cob、cox1和nad2基因含有正选择位点。基因编码偏好A/T含量高的密码子,27个高频密码子中,13个以A结尾,14个以T结尾。系统发育分析结果显示,灵芝属真菌主要分为2个聚类组,其中紫芝、狭长孢灵芝和G.wiiroense聚为一组;喜热灵芝、白肉灵芝和铁杉灵芝聚成一支,与树舌灵芝、梅氏灵芝、四川灵芝和亮盖灵芝构成姊妹类群,共同构成另一组。【结论】灵芝属真菌的线粒体基因组在进化过程中发生明显的遗传变异,基因组长度主要与内含子插入和删除有关,蛋白编码基因密码子使用偏性强。     

瘟病毒属3种病毒多聚蛋白酶基因同义密码子偏爱性分析

【目的】分析黄病毒科瘟病毒属的猪瘟病毒(CSFV)、牛病毒性腹泻病毒(BVDV)和羊边界病毒(BDV)3种病毒多聚蛋白酶基因密码子的使用特点,初步探讨3种病毒多聚蛋白酶基因的进化及对不同宿主的适应策略。【方法】运用EMBOSS软件的CHIPS和CUSP程序及CodonW在线分子生物学软件,分析CSFV、BVDV和BDV的有效密码子数(ENc),碱基GC含量(Gcall),密码子第1、2、3位GC含量(GC1、GC2、GC3),相对同义密码子使用度(RSCU)及同义密码子第3位的GC含量(GC3s),并对3种病毒多聚蛋白酶基因的ENc进行绘图分析。【结果】CSFV、BVDV、BDV 3种病毒多聚蛋白酶基因的ENc值分别为51.205,51.323和52.580,表明3种病毒多聚蛋白酶基因的密码子偏好程度均较低;上述3种病毒多聚蛋白酶基因的GC含量分别为46.84%,45.99%和45.71%,均不超过50%;除BDV多聚蛋白酶基因同义密码子含量表现为GC1>GC3>GC2外,CSFV和BVDV均表现为GC3>GC1>GC2,且除CSFV的GC3为52.02%外,所有病毒的GC1、GC2、GC3均小于50%,可见3种病毒间的碱基组成(GC含量)及ENc值差异较小。RSCU分析表明,与CSFV偏爱使用G/C结尾的密码子不同,BVDV和BDV更偏爱使用以A结尾的密码子,这可能与病毒对宿主的适应策略有关。ENc绘图分析(Nc-plot)说明,碱基组成对密码子使用模式影响较小。【结论】CSFV、BVDV、BDV 3种病毒多聚蛋白酶基因密码子的偏好性较低,这可能是病毒进化过程中对宿主的一种适应策略。     

芍药花色调控基因的密码子使用模式及其影响因素分析

【目的】芍药花色的优劣影响其观赏价值和商业价值,研究芍药花色调控基因的密码子使用偏好性和密码子使用模式的影响因素,为芍药花色调控基因在m RNA翻译、转基因设计、新基因表达与功能预测以及分子生物进化研究提供参考。【方法】根据前期芍药花色嵌合体品种‘金辉’转录组测序筛选的6 345个芍药花色调控基因,并根据CDS序列特征和大于300 bp原则进行过滤后最终获得的2 234个基因序列作为研究对象,利用Mobyle软件计算GC含量、第1与2位密码子的平均GC含量(GC12)、第3位密码子的GC含量(GC3s)、有效密码子数ENC、密码子适应指数CAI、相对同义密码子使用度RSCU等密码子偏性指标,其次进行中性绘图(GC12 vs.GC3)、ENC-GC3s绘图以及PR2(Parity Rule 2)绘图分析,并运用多元统计分析方法探讨突变压力和选择作用对密码子使用模式的影响程度,最后以5%CAI值作为高、低表达样本组,计算这两个样本组的同义密码子相对使用度,利用卡方检验Chi-square test分析两组之间的显著性差异来确定最优密码子。【结果】芍药花色相关基因的密码子GC3s含量为46.37%,大部分基因GC含量主要分布在30%—55%;中性绘图分析表明GC3s与GC12呈极显著的正相关(R2=0.202,P0.01);ENC-GC3s绘图表明大部分基因分布在标准曲线周围,也有一部分基因分布在标准曲线下方较远的位置,同时大部分基因(ENCexp-ENCobs)/ENCexp比值集中分布在0.0—0.4;PR2绘图分析显示密码子第三位T的使用频率高于A,C使用频率高于G,表明嘌呤(A和G)与嘧啶(T和C)的使用频率并不均衡;对应性分析COA(Correspondence Analysis)表明,第一轴上显示了38.09%的差异,其他3个轴分别为18.42%、15.09%、14.59%,表明芍药花色调控基因的密码子使用模式评价以第一轴(Axis 1)为主;突变压力和选择作用分析发现,第一主轴与GC3s、CAI的相关系数均达到极显著正相关(R2=0.736,P0.01;R2=0.286,P0.01);利用△RSCU和卡方显著性检验的方法,确定了21个为芍药花色调控相关基因的最优密码子,其中18个以G或C结尾,仅CGU、GGU等2个密码子以U结尾。【结论】芍药花色调控基因的最优密码子多数以G/C结尾,并且密码子使用模式主要受到碱基差异(R2=0.736)和基因表达水平(R2=0.286)共同作用的影响,其中碱基差异占主导因素。本研究了解了芍药花色调控基因的密码子使用模式情况,为通过密码子改造开展芍药花色遗传改良以及分子进化研究提供了一定的理论依据。     

濒危植物峨眉凤仙花叶绿体基因组分析

【目的】对峨眉凤仙花（Impatiens omeiana Hook.f.）叶绿体基因组的结构特征及系统发育进行研究，为其资源保护及开发利用提供理论依据。【方法】基于峨眉凤仙花叶绿体基因组序列，利用生物信息学软件，对叶绿体基因组进行组装、注释、基因特征、序列重复和系统发育分析。【结果】峨眉凤仙花完整的叶绿体基因组长度为152 527bp，共有130个基因，包括86个蛋白质编码基因、8个rRNA基因和36个tRNA基因，GC含量为37%，且具有保守的四分体结构，包括大单拷贝区、小单拷贝区各1个和2个相同的反向重复区域，其长度分别为83 150、17 903、25 737 bp，其中13个基因有1个内含子，2个基因有2个内含子。特征分析表明：峨眉凤仙花叶绿体全基因组中共检测到76个SSR序列，且多以A/T单核苷酸序列为主，其长度为10～91 bp；检测到50 842个密码子，其中以亮氨酸（Leu）最多，色氨酸（Tyr）最少；密码子偏好性分析显示33个RSCU≥1的密码子多数以A/U结尾。通过邻接法（NJ）构建系统发育树发现，峨眉凤仙花与贵州凤仙花亲缘关系最近，均属于棒凤仙花亚属植物。【结论】...     

五指毛桃叶绿体基因组结构与序列特征分析

【目的】阐明药食两用植物五指毛桃的叶绿体基因组结构及其系统进化关系,为其资源利用和产品开发等研究提供基因组信息。【方法】采用高通量测序技术对五指毛桃叶绿体基因组测序,并借助生物信息工具和软件进行序列拼接、注释、比对和系统进化分析。【结果】五指毛桃叶绿体基因组为环状双链四分体分子,长度160 340 bp, GC含量为35.9%,包含130个基因。该基因组含有26 679个密码子,偏好使用以A/T结尾的密码子;共有93个简单重复序列,其中以A/T形成的单、二核苷酸重复为优势重复基序。五指毛桃与其他榕属植物相比,叶绿体基因组序列同源性较高,碱基变异位点数量较少,主要分布在非编码区域如rpoB-trnC、trnT-trnL、rpl32-trnL等,与薜荔具有最高的序列相似度。【结论】五指毛桃叶绿体基因组具有典型的高等植物叶绿体基因组结构和基因组成,存在密码子偏好性,包含类型丰富的简单重复序列,且与同属植物薜荔的亲缘关系最近。     

13种植物actin基因的密码子使用特性分析

[目的]分析13种植物actin基因的密码子组成、密码子偏性及聚类关系,了解其密码子使用模式及影响密码子使用的因素,为深入研究分子进化及物种进化提供参考.[方法]运用CodonW 1.4.4软件分析13种植物的肌动蛋白基因(actin)密码子组成及使用参数,并对影响密码子偏性的因素进行研究.应用MEGA 4.1对13条基因的CDS序列进行聚类分析,采用SPSS 20.0进行密码子偏性的聚类分析.[结果]双子叶植物中actin基因的GC含量为45.0%~51.2%、GC3s含量为36.3%~53.8%,单子叶植物中GC含量为53.0%~58.3%、GC3s含量为60.9%~75.8%;actin基因在单子叶植物中偏爱G/C结尾的密码子,双子叶植物中偏爱A/T结尾的密码子.GC和GC3s与有效密码子数(ENC)呈极显著负相关(P<0.01),相关系数均为-0.906.ENC绘图分析结果表明,actin基因密码子偏性同时受突变和选择压力影响,单子叶植物受选择压力影响的程度大于双子叶植物.基于actin基因密码子偏性的聚类将单子叶植物高粱、玉米聚为一类,水稻、大麦、竹聚为一类,8种双子叶植物聚为一类.[结论]actin基因密码子偏性与碱基组成密切相关,其密码子偏性在单、双子叶植物间存在差异,依据密码子偏性的聚类能在一定程度上反映物种间的亲缘关系.     

峨眉含笑叶绿体基因组密码子的使用模式分析

为分析峨眉含笑叶绿体基因组密码子的使用偏性,揭示影响其密码子使用偏性的主要因素,首先从NCBI数据库中下载峨眉含笑叶绿体基因组完整序列,用于CDS序列筛选,然后利用CodonW 1.4.2和EMBOSS在线软件对峨眉含笑密码子的组成成分进行分析,从而得到GC含量、ENC、RSCU等参数,最后做中性、ENC和PR2绘图分析,研究导致密码子偏倚形成的主要原因,以期取得最优密码子。结果表明：峨眉含笑叶绿体基因组CDS序列的密码子在第三位上GC含量较低,平均值为31.26%。ENC为36.36～57.89,表明密码子的偏性较弱,ENC与GC₂和GC₃均呈极显著相关。有32个密码子的RSCU > 1,其中以G/C结尾的有3个,A/U结尾的有29个。中性绘图、ENC−plot和PR2−plot分析显示：由于选择压力的存在,使峨眉含笑叶绿体基因组密码子在一定程度上表现出偏好性使用。对应性分析显示：编码核糖体蛋白的基因位点分布比较零散,说明该类基因密码子的使用模式相差较大,在峨眉含笑叶绿体基因组中的表达率较高。该研究明确了峨眉含笑叶绿体基因组密码子使用偏性主要受选择压力的影响,并最终从峨眉含笑叶绿体基因组中筛选出UUU、CUU和UCU等14个最优密码子,且均以碱基A和U结尾。     

猕猴桃查尔酮合成酶基因CHS密码子偏好性分析

为了解猕猴桃查尔酮合成酶基因(Chalcone Synthase,CHS)的密码子使用特征,采用CodonW、SPSS、MEGA软件和EMBOSS在线程序等分析猕猴桃CHS密码子偏好性,并构建单双子叶植物CHS系统发育树。研究结果表明,猕猴桃CHS有效密码子数(ENc)、总G和C含量(GC)和密码子第3位上G或C含量(GC3s)分别为57.21、0.533和0.607;同义密码子相对使用度(RSCU)大于1.0的密码子数为30,但不存在偏好性极强的密码子;聚类分析结果发现,基于密码子使用偏好性分类可将单双子叶进行准确归类,而CHS序列相似性聚类结果并不理想;密码子使用频率分析表明,猕猴桃CHS与模式生物拟南芥和大肠杆菌基因组间密码子使用偏好性差异相对较小。研究认为,猕猴桃CHS密码子偏好性较弱,但倾向使用G和C,且以G或C结尾的密码子;单双子叶植物CHS间存在严格的密码子使用法则;此外,拟南芥和大肠杆菌可能是猕猴桃CHS遗传转化和异源表达较为理想的受体系统。     

楸树叶绿体基因组密码子偏性分析

【目的】分析楸树叶绿体基因组密码子使用偏性,为楸树乃至梓属植物叶绿体基因组学、物种进化特征及遗传多样性分析等研究提供理论依据。【方法】利用Codon W 1.4.2及EMBOSS expiorer的cusp和chips软件对51个楸树叶绿体基因的编码区（CDS）序列进行密码子相关参数、中性绘图、ENC-plot和PR2-plot等分析,并通过构建高、低表达偏性基因库,以ENC值和RSCU值为参考标准,最终筛选最优密码子。【结果】51个基因密码子的平均GC_all含量为38.48%,密码子不同位置GC含量为:GC₁（47.03%）>GC₂（39.93%）>GC₃（28.09%）,说明第3位的GC并非均匀分布;ENC为34.93~56.50,平均为46.87。GC₂与GC₁呈极显著正相关（P<0.01）,而GC₃与GC₁和GC₂无显著正相关（P>0.05）;ENC与GC_all、GC₁、GC₃和N均呈显著正相关（P<0.05,下同）;密码子数（N）与GC3和ENC呈显著正相关。RSCU>1的密码子共有30个,其中,以A或U结尾的密码子有29个。中性绘图分析结果显示,GC₃与GC₁₂相关系数为0.0188,二者无显著相关性,回归系数为0.143。ENC-plot分析结果显示,大部分基因偏离了ENC预期标准曲线,少数分布在预期曲线附近;有7个密码子分布在GC₃含量-0.05~0.15区间,其余区间共有44个密码子,而这44个基因离ENC标准曲线较远。PR2-plot分析结果显示,落在中线上的基因较少,平面图右下方的基因分布较多,表明在各碱基使用频率方面存在偏性,即T>A,G>C。各碱基使用频率存在偏性（T>A,G>C）,说明基因组密码子的偏性更多受选择效应的影响。最终筛选出7个为最优密码子UUU、CUU、CCU、AAA、GAU、CGA和GGU。【结论】楸树叶绿体基因组密码子偏性较弱,偏好使用A或U结尾的密码子,主要受选择压力的影响,受突变等其他因素影响较小。     

金银花大毛花叶绿体基因组密码子的偏好性分析

通过分析金银花大毛花叶绿体基因组密码子使用偏好性,探讨影响其密码子偏性形成的主要因素,为金银花系统发育及叶绿体基因组密码子进化研究提供参考。利用生物信息学手段对金银花大毛花50个叶绿体蛋白编码基因的密码子使用偏好性进行分析。金银花大毛花叶绿体基因组GC含量为39.18%,有效密码子数(ENc)为48.81,密码子的使用偏好性较弱;GC与GC1、GC2、GC3、GC3s、GC和CAI之间的相关系数均达到显著或极显著水平;同义密码子的相对使用频率(RSCU)大于1的密码子有30个,其中,28个以A/T结尾;综合ENc-plot分析、PR2-plot分析、中性绘图分析及对应性分析表明,在金银花大毛花叶绿体基因组密码子使用偏好性的形成过程中,选择压力和突变及其他因素共同发挥着重要作用。基于相对同义密码子使用分析,确定了30个高频密码子和5个最优密码子,最优密码子均以A或T结尾。本研究明确了金银花大毛花叶绿体基因组密码子使用偏好性受选择压力和突变等因素的共同影响,并筛选出5个最优密码子。     

猫尾草叶绿体基因组结构及其解析

【目的】阐明岭南药食两用植物猫尾草的叶绿体基因组结构特征，为猫尾草资源的保护、利用和开发提供理论依据。【方法】采用高通量测序技术开展猫尾草叶绿体基因组测序，并通过生物信息方法进行拼接、注释、解析以及系统进化分析。【结果】猫尾草叶绿体基因组是149 774 bp的环状双链四段式分子，包含128个基因，GC含量为38.2%。猫尾草叶绿体基因组含有26 015个密码子，偏好以A或T结尾；存在110个简单重复序列，以A或T单核苷酸重复居多。序列比对和进化分析显示猫尾草与同属狸尾豆的亲缘关系最近。【结论】首次报道猫尾草叶绿体基因组的全序列，并明确其结构特点，为猫尾草的栽培育种、遗传多样性和资源利用等奠定基础。     

瓜螺线粒体全基因组密码子偏好性分析

【目的】探究瓜螺（Melo melo）线粒体全基因组密码子使用偏好性,明确影响密码子偏好性的因素,为提高基因表达效率及了解瓜螺线粒体基因组的特性和进化方向提供科学依据。【方法】以瓜螺线粒体全基因组序列为研究对象,从中筛选出11条长度大于300 bp且以ATG为起始密码子的非重复基因序列,利用CodonW 1.4.2、Excel 2017及SPSS 22.0等软件进行瓜螺线粒体全基因组密码子偏好性分析。【结果】不同基因序列密码子位置所对应的G+C含量不同,其中,GC₁（密码子第1位所代表的G+C含量）、GC₂（密码子第2位所代表的G+C含量）、GC₃（密码子第3位所代表的G+C含量）的变化范围分别为35.20%~51.90%、16.40%~40.60%和13.10%~34.90%,对应平均值为41.15%、32.29%和23.71%,GC₁和GC₂明显高于GC₃;密码子适应指数（CAI）在0.110~0.180,平均为0.137;密码子偏好性指数（CBI）在-0.272~-0.090,平均为-0.175;最优密码子使用频率（FOP）在0.216~0.339,平均为0.277;有效密码子数（ENC）在36.84~53.75,平均为45.59;总平均亲水性（GRAVY）在0.480~1.394,平均为0.968。有30个同义密码子的相对使用度（RSCU）大于1.00,且主要以A/U结尾。中性绘图分析发现,GC₁和GC₂平均值（GC₁₂）与GC₃的相关系数为-0.438,相关性不显著（P>0.05,下同）;ENC-plot绘图分析结果显示,各散点（基因）均位于标准曲线附近,其中,NAD4、COX2和NAD3基因位于标准曲线上方,COX1、ATP6、NAD1、NAD6、COB、NAD5、COX3及NAD2等8个基因位于标准曲线下方;在对应性分析中,第一、第二、第三、第四向量轴的贡献率分别23.12%、15.19%、13.17%和9.63%,前4轴累计差异为61.12%。综合高频密码子与高表达密码子,最终筛选出4个最优密码子（UUU、GUU、ACU和CAU）,且均以U结尾。【结论】瓜螺线粒体全基因组密码子偏好以A/U结尾的形式,筛选出的4个最优密码子（UUU、GUU、ACU和CAU）均以U结尾。除了自然选择在瓜螺线粒体全基因组密码子偏好性形成中占主导地位外,突变压力和碱基组成也影响其密码子偏好。因此,可通过瓜螺目的基因密码子优化,有效提高外源基因表达且优化性状决定基因,从而加速瓜螺优良亲本培育的进程。     

基于RNA-seq数据的栽培种花生SSR位点鉴定和标记开发

【目的】鉴定花生RNA-seq数据中的SSR位点,明确转录组中SSR位点的分布和结构特点,开发与花生基因相关联的SSR标记,为花生重要功能基因的挖掘、等位变异研究和分子标记辅助育种奠定基础。【方法】根据栽培种花生全生育期中22种不同类型的组织RNA-Seq数据,使用MISA软件分析SSR位点分布及特征,采用Primer3设计基因关联的SSR引物,并利用电子PCR软件对引物的质量进行检测,随机合成38对引物,进行多态性检测。【结果】从52 280条转录本中共鉴定19 143个SSR位点,分布于14 084条转录本,发生频率为26.94%。重复单元类型为单核苷酸—五核苷酸,以单核苷酸和三核苷酸为重复单元的SSR位点数最多,分别占位点总数的39.24%和38.40%。各重复单元优势基序类型分别为A/T、AG/CT、AAG/CTT、AAAG/CTTT和AACAC/GTGTT,占所在重复单元中的比例分别为97.62%、72.01%、30.96%、24.59%和16.67%。重复单元的重复次数为5—47次,单个SSR位点的长度的分布范围为10—47 bp,基序长度主要集中在10—14 bp;复合SSR位点的长度范围为21—249 bp,以31—40 bp为主。鉴定的SSR位点中共有13 477个SSR位点可以进行引物设计,其中5 020条转录本序列对应到特定的基因,共包含5 859个可进行引物设计的SSR位点,这些SSR位点在A基因组和B基因组共20条染色体上不均匀分布,其中B03染色体上SSR位点最多,为484个。对特定基因SSR引物进行电子PCR检测,在A.duranensis、A.ipaensis、A.hypogaea基因组中有效扩增位点分别为4 468、4 929和10 188个,有效引物数分别为3 968(67.74%)、4 232(72.25%)和5 174(88.33%)对,在A. hypogaea基因组中,SSR引物扩增位点主要以2个位点为主,其中,有1 477对引物单位点扩增。根据SSR引物扩增位点在栽培种花生基因组中的位置信息绘制了SSR位点的物理图谱。在38对SSR引物中,共有35对(92.1%)SSR引物可以扩增出清晰的条带,其中,有11对(28.9%)SSR引物在2个品种间扩增出差异条带。【结论】鉴定了13 477个可进行标记开发的SSR位点,开发、检测了5 859个基因相关SSR标记,在栽培种花生基因组中具有较高的扩增效率,并构建了基因相关SSR位点的物理图谱。