玉米秸秆消化率全基因组关联分析

48h体外干物质消化率（48h in vitro dry matter digestibility, 48h IVDMD）是衡量青贮玉米品质的重要指标。为了初步探究玉米秸秆消化率的分子遗传机理,以341份玉米自交系为材料,于2018年在沈阳和通辽种植,收获后测定秸秆48h IVDMD。利用全基因组重测序获得的6 276 612个高质量SNPs进行全基因组关联分析,共检测到153个与玉米秸秆消化率显著相关的SNPs位点（P<1.0×10^-6）,4个SNPs显著水平在P<1.0×10^-8以上;共找到38个秸秆消化率的候选基因,主要涉及细胞生长发育、防御反应和信号转导等生物学功能。     

瓜类褪绿黄化病毒及其系统进化分析

【目的】检测新疆南疆7个县/市瓜类褪绿黄化病毒(Cucurbit chlorotic yellows virus, CCYV),进行系统发育分析,为新疆南疆CCYV预警和防控提供参考。【方法】于2019年,从中国新疆南疆7个县/市采集带有黄化特征的51份甜瓜叶片样品,RT-PCR进行CCYV的检测,进行特异性片段的克隆、测序和系统发育分析。【结果】中国新疆巴楚县、阿克苏市、莎车县、伽师县、疏勒县、疏附县均未检测到CCYV,而洛浦县的13份样本中,8份检出为CCYV阳性,检出率为61.53%,洛浦县CCYV特异性片段克隆测序获得了685 bp的核酸序列,与GenBank中的CCYV分离物外壳蛋白(Coat protein, CP)的一致性达到100%。所得序列与中国新疆(吐鲁番)、中国其他省、苏丹、日本、塞浦路斯、黎巴嫩的CP基因聚在I组(Group),沙特阿拉伯的分离物聚在Ⅱ组,伊朗的分离物聚在Ⅲ组。【结论】CCYV在中国新疆南疆洛浦县甜瓜产区发生,与中国新疆吐鲁番、中国其他省及周围国家的CCYV分离物亲缘关系很近,且群体遗传变异很小。     

猪乙型脑炎病毒GZ0409-31分离株E基因遗传进化分析

本研究采用RT-PCR方法对分离自贵州省发病猪睾丸组织中的流行性乙型脑炎病毒GZ0409-31株的E基因进行克隆和序列测定,利用生物信息学软件对E基因核苷酸序列和推导的氨基酸序列进行分子遗传进化关系分析。结果表明,乙型脑炎病毒GZ0409-31毒株与基因Ⅲ型的乙脑病毒毒株属同一分支,遗传进化关系较近。     

5种蔷薇科树种多酚氧化酶比较基因组学分析

蔷薇科植物中有多种具有重要价值的水果,如苹果、梨、桃、草莓和黑树莓等。这些水果普遍容易发生由多酚氧化酶(PPO)介导的酶促褐变而导致重大经济损失,从全基因组角度分析比较PPO基因家族,有助于加深对PPO基因家族和功能的认识。采用比较基因组学的方法,对这5种植物的PPO基因家族进行了基因鉴定、染色体定位、编码蛋白的亚细胞定位、内含子统计分析、基因系统进化、基因倍增与丢失等特征分析。从这5种植物中,共计鉴定出了42个PPO基因,其中6个基因被认定是假基因;绝大多数基因编码的蛋白定位于叶绿体,2个定位于线粒体,3个是分泌型蛋白;PPO基因在染色体上有串联重复和散布两种形式;9个基因具有内含子,聚类结果显示可以将它们分为6个类型,每个基因类型在进化过程中都发生过基因倍增和丢失;同型基因内含子的位置和大小具有相似性。这些结果揭示蔷薇科植物PPO基因内含子是伴随着基因倍增产生的,基因倍增也是推动PPO基因多样化的重要动力,PPO基因倍增和丢失差异导致PPO基因数量在不同物种之间产生差异。     

香蕉与B基因组相关的SCAR标记研究

香蕉(Musa L.)是由2个二倍体野生种Musa acuminata Colla(AA基因型)和Musa balbisiana Colla(BB基因型)种内或种间杂交进化而来,其B基因组中带有重要的优良基因。利用与香蕉B基因组相关的gypsy-IRAP分子标记,成功开发了一对SCAR引物,适用于鉴定尖叶蕉(AAw)、长梗蕉(BB)、香牙蕉(AAA)、贡蕉(AAcv)、大蕉、粉蕉(ABB)、粉大蕉(ABB)、龙牙蕉(AAB)以及四倍体香蕉(AAAB)等是否含有B基因组。     

大花君子兰叶绿体基因组及其特征

采用Illumina MiSeq测序平台对大花君子兰（Clivia miniata）叶片总DNA进行测序，通过组装获得了其叶绿体基因组（cpDNA）全长序列（158 114 bp）。对其cpDNA注释得到135个基因，包含87个蛋白编码基因、40个tRNA基因和8个rRNA基因。采用生物信息学方法对获得的cpDNA进行简单序列重复（SSR）分析和密码子偏好性分析。结果显示：①大花君子兰cpDNA中共有61个SSR位点，其中单核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸和六核苷酸重复数分别为38、9、2、8、3和1个，多数SSR分布在基因间隔区；②大花君子兰cpDNA密码子偏爱以A或U（T）结尾，亮氨酸使用频率最高，半胱氨酸使用频率最低。基于24种植物的cpDNA全长和23种植物的叶绿体ycf2基因序列进行系统发育分析，结果显示大花君子兰与石蒜科植物在同一分支，显示最近的亲缘关系，支持大花君子兰属于石蒜科。基于叶绿体ycf2的系统发育分析结果与基于cpDNA全长的系统发育分析研究结果大部分相同，支持ycf2基因可以代替cpDNA全长用于植物系统发育分析。     

芸薹属及其他异源多倍体植物的基因表达特征研究进展

远缘杂交及异源多倍化导致许多重要作物的起源与进化，而芸薹属栽培异源四倍种是研究作物异源多倍化的模式系统之一。异源多倍体是如何调控及协调来自不同二倍体祖先的不同基因组的遗传行为及基因表达，是过去二十年间的研究热点和重要的生物学问题。利用不断发展的分子生物学技术，一方面揭示出芸薹属及其他多倍体物种基因组表现出动态性质，即在形成初期及长期进化过程中持续发生遗传及表观遗传的变化；另一方面发现异源多倍化过程中伴随着大量的基因表达模式改变，包括非加性表达、超亲表达、表达水平显性、部分同源偏向表达、基因剂量平衡效应等现象。上述基因组结构、表观遗传改变以及基因表达模式的调控，使新产生的多倍体得以成功进化为新物种。     

全基因组关联分析阈值的快速算法

逐个SNP的全基因组关联分析中的统计量不仅依赖于遗传效应,而且依赖于SNP,因此统计量不能直接应用于推断无遗传效应的零假设中。检测不同样本的QTN,常用的统计量除了卡方统计量,还有t统计量、F统计量和标准正态统计量。首先给出了各个统计量之间的关系,接下来针对冗余参数背景下的全基因组关联分析,提出了检验统计量阈值的快速计算方法。再次利用获得的高通量SNP的统计概率构建卡方统计量,进而估计临界值。最后模拟不同样本的阈值并与文献上提出的阈值算法比较。大量模拟实验证明,提出的方法快速而且有效。     

具有III型分泌系统的Pseudomonas moraviensis FP1761植物益生功能及其全基因组分析

 本研究从节瓜根际分离到一株具有防病促生功能的荧光类假单胞菌FP1761。该菌株对部分植物病原真菌和细菌具有拮抗能力,能够解钾、解有机磷和无机磷,可产生氨、蛋白酶、嗜铁素、吲哚乙酸。生物测定表明菌株FP1761可显著促进小麦生长。生理生化、平均核苷酸相似度、16S rDNA和多基因分析将FP1761鉴定为摩拉维亚假单胞菌(Pseudomonas moraviensis)。菌株FP1761基因组草图全长6.12 Mb,(G+C)含量为59.9%,共编码5467个基因序列。将该菌株与种内3个代表性菌株进行泛基因组和核心基因组分析,共产生 4 357个共有基因,菌株FP1761特有基因327个。利用antiSMASH对菌株次生代谢基因簇进行预测,发现其含有8个潜在的次生代谢产物基因簇。其中两个基因簇与嗜铁素pyoverdine合成相关,未见聚酮类合成基因。基因组分析发现,该菌株具有与病原性假单胞菌相似的III型分泌系统,但丢失了效应蛋白调控因子hrpS和转运相关的hrpH和hrpK1基因。对全基因组扫描,菌株FP1761仅保留了病原性假单胞菌的保守效应蛋白AvrE和HopAA1-1。FP1761是目前已发现的唯一具有III型分泌系统的摩拉维亚假单胞菌。本研究表明摩拉维亚假单胞菌FP1761具有潜在的植物防病促生功能,但其III型分泌系统与植物益生互作机制有待进一步解析。     

矮牡丹与芍药属其他5个种叶绿体基因组特征的比较

【目的】中国特有的野生牡丹一直被国内外视为珍贵的种质资源。野生矮牡丹被认为是现代栽培牡丹品种重要的祖先种之一,开展矮牡丹在内的芍药属植物叶绿体基因组(cp DNA)特征分析对阐明牡丹系统进化、培育和改良栽培品种具有重要的理论与实践价值。【方法】在矮牡丹叶绿体高通量测序的基础上,从NCBI数据库下载凤丹牡丹、大花黄牡丹、滇牡丹、川赤芍和草芍药的cp DNA数据,利用Geneious 8. 0、EMBOSS 6. 4. 0等软件,对芍药属6个种的cp DNA进行对比分析。【结果】矮牡丹cp DNA序列共152 628 bp,共有112个基因,使用25 988个密码子,编码蛋白78个;有19个基因(包括4个rRNA、7个tRNA、8个蛋白编码基因)在IR(反向重复)区重复。共搜索到143个SSR位点,单核苷酸重复基序位点最多,为116个(占81. 12%),没有六核苷酸重复基序。尽管芍药属叶绿体基因组比较保守,但不同种间IR和LSC(大单拷贝区)的边界位置仍有一定变化,凤丹牡丹LSC/IR的rpl2基因有718 bp延伸至LSC区域,而其他种的rpl2基因均完整地位于IR区。【结论】从基因组大小和基因内容来看,芍药属cp DNA高度保守;草芍药与滇牡丹cp DNA最大,为152 698 bp;凤丹牡丹cp DNA最小,为152 153bp。矮牡丹cp DNA蛋白编码基因密码子偏好使用A/T碱基,SSRs位点的碱基组成也偏好使用A/T碱基,143个SSRs位点中,A/T组成的位点有134个;矮牡丹cp DNA SSRs分布具有不均匀性,14个SSR位点位于IR区段,103个位于LSC区段,26个位于SSC区段。IR边界分析显示,芍药属LSC/IRb的边界变化是IR区扩张与收缩的主要原因。研究结果为芍药属植物的系统进化与栽培起源等研究提供支持,对芍药属植物分子标记开发及优良品种选育具有参考价值。