普通菜豆镰孢菌枯萎病抗性种质资源筛选及全基因组关联分析

普通菜豆镰孢菌枯萎病是严重制约菜豆(Phaseolus vulgaris)产量的主要病害之一。采用下胚轴双孔注射法对601份普通菜豆种质资源进行枯萎病抗性鉴定, 共筛选出4份高抗材料。在此基础上, 基于分布在全基因组上的3 765 456个单核苷酸多态性(SNP)标记, 进行全基因组关联分析, 以P<1×10^-5为阈值。结果检测到57个显著关联的SNP位点, 分布于1、2、6、8和11号染色体上; 共获得8个显著关联区域, 其中位于1号染色体上的区域1包含SNP最多(48个), 最显著SNP P值为2.18E-07。在8个显著关联区域中, 共检测到186个基因, 其中157个基因有注释信息, 编码过氧化物酶、抗病蛋白、转录因子和蛋白激酶等。结合KEGG富集分析和序列同源性比对, 鉴定出9个候选基因可能与抗性相关。     

利用传递不平衡分析在多群体中鉴别猪脐疝易感基因

在本实验室前期利用白色杜洛克×二花脸F₂资源家系开展脐疝易感位点全基因组扫描定位的基础上,文章在7号染色体上的SWR1928和10号染色体上的SW830易感标记区域,结合脐疝发病机制在多群体中进行脐疝位置功能候选基因的筛选和易感位点的精细定位。在两个显著关联的微卫星位点区域搜寻到12个位置功能候选基因,采用比较测序法,选取12个候选基因内共计40个SNP位点在白色杜洛克×二花脸资源家系F₂/F₃脐疝群体中进行基因分型,利用Plink v1.07软件对基因型数据进行质量控制和传递不平衡(Transmission disequilibrium test,TDT)分析。结果表明,IL16(Interleukin 16)基因中的g.708C>T位点和CDC73(Cell division cycle 73)基因中的g.10664G>A位点与脐疝的关联性达到显著水平(P<0.05)。对这两个位点在西方商业猪种脐疝患病家系中进行基因分型和TDT验证分析,发现CDC73基因中的g.10664G>A位点仍与猪脐疝呈显著关联(P<0.05)。同时对CDC73基因中与资源家系脐疝呈弱相关的两个SNP位点g.10546A>G和g.10811A>G在西方商业猪种中进行TDT验证分析,发现这两个SNP位点与商业猪种脐疝发生的关联性达到极显著水平(P<0.01)。根据文章的分析结果,结合脐疝发生的生理机制及CDC73基因的生物学功能,推测CDC73基因可能为猪脐疝发生的易感基因。     

利用紧缩线性模型和贝叶斯模型对猪总产仔数和产活仔数性状的全基因组关联研究

全基因组关联分析策略已逐渐成为家畜重要经济性状研究的强有力工具。文章使用猪60K SNP芯片对一个具多胎繁殖性状记录的商业母猪群(n=820)进行分型检测,共计57 814个SNP通过设定质控标准。主成分分析显示群体内不存在显著的群体分层现象,而后分别运用两种统计模型Compressed Mixed Linear Model(GAPIT程序包)、Bayes CPi(GenSel软件)进行第1和第2胎次总产仔数和产活仔数性状的全基因组关联分析。从两种分析方法所得结果中各取最显著的50个SNP位点进行比较:对于第1胎次总产仔数,两种方法分析结果存在31个重合SNP位点,对于第1胎次产活仔数,有20个重合SNP位点;且两种统计分析结果中最显著的SNP位点都在另一方法中得到验证。与第1胎次总产仔数显著关联的SNP位于1、2、3、7、13、16和18号染色体,与第1胎次产活仔数显著关联的SNP位于1、3、4、13和16号染色体上的11个区域内。在1、3、13和16染色体上共有5个区域同时与这两个性状显著关联。与第2胎次总产仔数和产活仔数显著关联的区域主要位于7、10、12、13、14和16号染色体的6个重叠区域内。     

小麦耐低磷相关性状的全基因组关联分析

通过对小麦耐低磷相关性状进行全基因组关联分析(GWAS,genome-wide association study),挖掘与小麦耐低磷性显著相关的单核苷酸多态性标记(SNP,single nucleotide polymorphism)位点及候选基因,为小麦耐低磷性状的遗传基础和分子机制研究提供理论参考。本试验以198份黄淮麦区小麦品种(系)为试验材料,设置低磷和正常磷营养液水培试验,利用小麦35K芯片对分布于小麦全基因组的11896个SNP,采用Q+K关联模型对小麦耐低磷性相关性状进行关联分析。结果表明,小麦耐低磷性状表现出广泛的表型变异,变异系数为15.65%~26.59%,多态性信息含量(PIC,polymorphic information content)为0.095~0.500。群体结构分析表明,试验所用自然群体可分为2个亚群,GWAS共检测到67个与小麦耐低磷相关性状显著关联的SNP位点(P≤0.001),这些位点分布在除3A、3B和3D以外的18条染色体上,单个SNP位点可解释5.826%~9.552%的表型变异。在这些显著位点中有4个SNP位点同时关联到了2个不同的耐低磷性状。对67个SNP位点进行发掘,筛选到7个可能与小麦耐低磷性有关的候选基因。TraesCS6A02G001000和TraesCS6A02G001100在锌指合成中有重要作用;TraesCS6A02G118100可能为低磷胁迫诱导基因;TraesCS5D02G536400、TraesCS1B02G154200和TraesCS5D02G536500与低磷胁迫相关酶类基因家族有关;TraesCS1D02G231200与植物DUF 538结构域蛋白有关,是植物胁迫相关调控蛋白候选基因。     

普通菜豆抗菜豆象性状的全基因组关联分析

普通菜豆(Phaseolus vulgaris)具有丰富的营养价值, 菜豆象(Acanthoscelides obtectus)是危害菜豆的主要害虫, 利用抗虫种质资源防治菜豆象是最安全且经济有效的方法。该研究利用改良的室内人工接虫方法, 对625份普通菜豆种质资源进行2次菜豆象抗性重复鉴定, 筛选出2份抗性稳定且种子受害率均在10%以下的高抗种质。利用种子受害率和蛀孔总数的表型数据, 基于3 767 432个SNP标记进行全基因组关联分析, 鉴定出15个与种子受害率相关的显著关联遗传位点, 8个与蛀孔总数相关的显著关联位点, 解释了4.54%-5.56%的表型变异。在候选位点筛选出包括编码蛋白酶抑制剂、凝集素和过氧化物酶等在内的20个与抗虫防御相关的候选基因。     

大豆需磷关键时期磷高效利用遗传位点挖掘

土壤有效磷缺乏已成为影响大豆产量和品质的重要因素,深入挖掘大豆需磷关键时期磷高效利用遗传位点成为实现其分子遗传改良的重要前提。鉴于此,本研究利用SoySNP6K(5403个SNP标记)通过全基因组关联分析挖掘大豆需磷关键时期磷高效利用10个相关性状遗传位点,结果发现,T1关键时期(四叶期)适磷条件检测到78个关联SNP,以根系干重与植株总干重SNP较多;低磷条件检测到134个关联SNP,以植株总干重检测SNP最多,并在8号、13号、20号染色体分别检测到同时控制地上部干重和总干重、地上部鲜重和干重与总干重SNP簇;T2关键时期(六叶期)适磷条件检测到83个SNP,以株高和地上部干重检测SNP较多,低磷条件检测到53个SNP,以株高和根冠比SNP较多,并在18号染色体检测到同时控制根干重和总干重SNP簇,在11号、16号、18号染色体分别检测到3个一因多效SNP;上述关联SNP中有9个SNP同时在T1与T2时期被检测到,分别与地上部鲜重、干重、根冠比、株高等关联,为大豆磷素高效利用分子标记辅助育种以及候选基因克隆等提供了依据。     

京海黄鸡体重性状全基因组关联分析

体重性状是肉鸡重要的经济性状。为了寻找可用于京海黄鸡体重性状遗传改良的分子标记及候选基因,本文以400只京海黄鸡核心群母鸡为基础,测定了0~14周龄体重,利用简化基因组测序技术(Specific-locus amplified fragment sequencing, SLAF-seq)对京海黄鸡体重性状进行全基因组关联研究(Genome-wide association stndy, GWAS),筛选与京海黄鸡体重性状相关的SNPs位点。结果共检测到100个与京海黄鸡体重相关的SNPs位点,其中15个位点效应达到全基因组显著水平(P<1.87E-06),85个位点效应达到全基因组潜在显著水平(P<3.73E-05)。通过筛选每个显著SNP周围1 Mb区域内的基因,共找到9个可能的候选基因,其中FAM124A(Family with sequence similarity 124A)、QDPR(Quinoid dihydropteridine reductase)、WDR1(WD repeat domain 1)和SLC2A9(Solute carrier family 2 (facilitated glucose transporter), member 9) 4个基因可能是影响体重性状的重要候选基因。同时还发现,4号染色体75.6~80.7 Mb区域集中了大部分与京海黄鸡中后期体重性状显著相关的SNPs位点,该区域可能是影响京海黄鸡中后期生长体重的重要候选区域。     

水稻抽穗期QTL定位及候选基因分析

水稻(Oryza sativa)抽穗期是决定产量和品质的重要性状, 在育种、制种及引种驯化过程中发挥重要作用。将热研2号(O. sativa subsp. japonica cv. ‘Nekken2’)和华占(O. sativa subsp. indica cv. ‘HZ’)杂交获得F₁代, 经连续多代自交得到120个重组自交系(RILs)群体。在常规水肥管理条件下, 对120个RILs株系的抽穗时间进行统计分析。利用已构建好的高密度遗传图谱, 对水稻抽穗期相关性状进行QTL定位分析, 结果共检测到11个QTLs, 分别位于第1、3、4、5、6、8和12号染色体上, 其中1个LOD值高达5.75。通过分析QTLs区间内的候选基因, 筛选出可能影响两亲本抽穗期的相关基因, 并利用实时定量PCR进行基因表达量分析, 发现LOC_Os03g03070、LOC_Os03g50310、LOC_Os03g55389、LOC_Os04g55510、LOC_Os08g07740和LOC_Os08g01670共6个基因在双亲间的表达量差异显著, 其中LOC_Os03g50310在Nekken2中的表达量比HZ高3.6倍。对双亲间候选基因LOC_Os03g50310进行测序分析, 发现该基因在5′UTR、CDS区及3′UTR存在4处差异, 其中CDS区的单核苷酸多态性(SNP)差异引发单个氨基酸改变。研究通过挖掘水稻抽穗期QTL位点为进一步克隆水稻抽穗期相关基因和品种选育提供了新线索。     

冠状动脉粥样硬化性心脏病8号染色体基因扫描

冠状动脉粥样硬化性心脏病(Coronary heart disease,CHD)的全基因组扫描研究在世界各大研究中心展开,关于CHD易感基因位点的报道多集中于1号、3号、9号、11号、16号染色体,对8号染色体的研究报道甚少。在汉族人群中未见关于CHD的8号染色体的基因扫描研究。文章旨在查找汉族人群中CHD易感基因位点,选取8号染色体上间隔10 cM遗传距离的13个微卫星遗传位点,采用DNA混合池的方法对CHD患者组156例和正常对照组1 000例DNA样本进行基因扫描,经卡方检验分析患者组和对照组每个位点的等位基因频率差异。发现在患者组与对照组中D8S264位点(8p23.3-p23.2)及D8S285位点(8q12.1)的等位基因频率差异有统计学意义(P<0.05)。汉族人群中CHD患者8号染色体上8p23.3-p23.2、8q12.1区域可能存在CHD易感基因,需要进行候选基因突变筛查。     

丙酸通路基因多态性与猪肉质及胴体性状的关联分析

为了挖掘新的猪肉品质及胴体性状的候选基因,揭示猪肉质及胴体性状的遗传机制,文章将丙酸代谢通路作为候选通路,将通路内基因与猪肉质及胴体性状进行关联分析。实验采用37头三元杂交商品猪作为研究对象,首次针对丙酸通路中7个基因的36个SNP位点利用SNaPshot方法进行基因分型,分别用最小二乘模型及MB-MDR模型与肉质及胴体性状进行关联分析。结果发现,基因PCCB、MUT、MCEE及ACSS2上的4个SNP位点分别与肌内脂肪含量、背膘厚等性状显著相关(P<0.05),ACSS2与猪脂肪含量显著相关;MCEE及MUT与猪的背膘厚显著相关;PCCB基因与脂重显著相关。通过MB-MDR方法检测到多个SNP位点具有互作效应,并与背膘厚、水分含量、脂肪含量显著相关(P<0.05)。另外,丙酸代谢通路中的基因间的互作效应对猪肉品质有显著影响。     

水稻叶片水势的QTL定位与候选基因分析

为探究叶片水势(LWP)相关基因在水稻(Oryza sativa)抗旱中的作用及其遗传机制,以热研2号(Nekken2)和华占(HZ)为亲本以及构建的120个重组自交系(RILs)群体为实验材料,对水稻分蘖期叶片水势进行检测,并利用前期基于高通量测序构建的分子遗传连锁图谱进行数量性状基因座(QTL)分析。结果表明,共检...     

基于CSSLs群体定位和图位克隆水稻长芒基因GAD1-2

水稻是世界上最早驯化的重要粮食作物之一。水稻芒可以保护水稻种子不被鸟琢食,是水稻重要的驯化性状之一。芒在野生稻中普遍存在,对野生稻的生存和传播至关重要,然而在驯化和人工选择过程中该性状逐渐被淘汰。定位和克隆水稻长芒相关基因是研究水稻芒驯化遗传机制的基础。本研究以籼稻恢复系东南恢810为受体、漳浦野生稻为供体构建的146个染色体片段置换系(chromosome segment substitution lines, CSSLs)为研究材料,调查了146个CSSLs株系和双亲的芒长,结果表明在4个置换系中检测到1个控制水稻芒长主效基因GAD1-2,位于水稻第8号染色体;利用重叠代换作图法,将GAD1-2定位在Ind8-10和RM4936标记之间,遗传距离约为4.75 Mb。选择分离群体中的显性单株,利用开发的标记,最终将GAD1-2 基因定位在两个 Indel 标记之间,两者间的物理距离约为27 kb,该区域内只有两个候选基因Os08g0485500和Os08g0485400。经测序和分析表明,Os08g0485500是GAD1-2的候选基因,GAD1-2在保守的ORF区域存在6个碱基缺失,导致丝氨酸和半胱氨酸这两个氨基酸缺失,从而表现长芒的性状;在Os08g0485500基因位点已克隆了1个控制水稻芒长的GAD1基因,推测GAD1-2与GAD1为等位基因。本研究为进一步理解水稻起源演化和水稻芒长发育基因的遗传机制奠定了基础。     

Comparative genetic analysis of a wheat seed dormancy QTL with rice and <Emphasis Type="Italic">Brachypodium</Emphasis> identifies candidate genes for ABA perception and calcium signaling

Wheat preharvest sprouting (PHS) occurs when seed germinates on the plant before harvest resulting in reduced grain quality. In wheat, PHS susceptibility is correlated with low levels of seed dormancy. A previous mapping of quantitative trait loci (QTL) revealed a major PHS/seed dormancy QTL, QPhs.cnl-2B.1, located on wheat chromosome 2B. A comparative genetic study with the related grass species rice (Oryza sativa L.) and Brachypodium distachyon at the homologous region to the QPhs.cnl-2B.1 interval was used to identify the candidate genes for marker development and subsequent fine mapping. Expressed sequence tags and a comparative mapping were used to design 278 primer pairs, of which 22 produced polymorphic amplicons that mapped to the group 2 chromosomes. Fourteen mapped to chromosome 2B, and ten were located in the QTL interval. A comparative analysis revealed good macrocollinearity between the PHS interval and 3 million base pair (mb) region on rice chromosomes 7 and 3, and a 2.7-mb region on Brachypodium Bd1. The comparative intervals in rice were found to contain three previously identified rice seed dormancy QTL. Further analyses of the interval in rice identified genes that are known to play a role in seed dormancy, including a homologue for the putative Arabidopsis ABA receptor ABAR/GUN5. Additional candidate genes involved in calcium signaling were identified and were placed in a functional protein association network that includes additional proteins critical for ABA signaling and germination. This study provides promising candidate genes for seed dormancy in both wheat and rice as well as excellent molecular markers for further comparative and fine mapping.     

玉米叶形相关性状的Meta-QTL及候选基因分析

叶长、叶宽、叶面积及叶夹角不仅影响玉米(Zea mays)光合效率, 也是株型的重要构成因素。通过对620个叶形QTL进行整合, 构建不同遗传背景下的叶形QTL整合图谱, 利用元分析发掘出22个叶长、22个叶宽、12个叶面积以及17个叶夹角mQTL; 进一步运用生物信息学手段, 确定44个与叶片发育密切相关的候选基因。分析发现, 仅有NAL7-like、YABBY6- like和GRF2等13个基因位于mQTL区间内, 而玉米中已克隆的KNOTTED1、AN3/GIF1、rgd1/lbl1、mwp1及SRL2-like、HYL1-like和CYCB2;4-like等水稻(Oryza sativa)和拟南芥(Arabidopsis thaliana)叶形同源基因位于未被整合的QTL内; 对44个候选基因在叶片长、宽、厚发育过程中基部-末端、中央-边缘、远轴-近轴的调控机理进行归纳分析, 发现玉米中仅有少数几个候选基因被报道, 揭示了叶形发育的部分分子机理。因此, 对玉米叶形相关mQTL/QTL及基因进行全面深入的分析, 不仅有助于增加对其遗传结构的了解, 发掘更多候选基因, 阐明叶形发育和形成的分子机制, 还可为耐密理想株型的分子标记辅助选择提供依据。     

Identifying genetic components controlling fertility in the outcrossing grass species perennial ryegrass (Lolium perenne) by quantitative trait loci analysis and comparative genetics

Mutational load and resource allocation factors and their effects on limiting seed set were investigated in ryegrass by comparative mapping genomics and quantitative trait loci (QTL) analysis in two perennial ryegrass (Lolium perenne) mapping families sharing common genetic markers. Quantitative trait loci for seed-set were identified on chromosome (LG) 7 in both families and on LG4 of the F2/WSC family. On LG7, seed-set and heading date QTLs colocalized in both families and cannot be unequivocally resolved. Comparative genomics suggests that the LG7 region is syntenous to a region of rice LG6 which contains both fertility (S5(n)) and heading date (Hd1, Hd3a) candidate genes. The LG4 region is syntenous to a region of rice LG3 which contains a fertility (S33) candidate gene. QTL maxima for seed-set and heading date on LG4 in the F2/WSC family are separated by c. 8 cm, indicating distinct genetic control. Low seed set is under the control of recessive genes at both LG4 and LG7 locations. The identification of QTLs associated with seed set, a major component of seed yield in perennial ryegrass, indicates that mutational load associated with these genomic regions can be mitigated through marker-assisted selection.     

Dissecting of the FHB resistance QTL on the short arm of wheat chromosome 2D using a comparative genomic approach: from QTL to candidate gene

Colinearity in gene content and order between rice and closely related cereal crops has been a powerful tool for gene identification. Using a comparative genomic approach, we have identified the rice genomic region syntenous to the region of the short arm of wheat chromosome 2D, on which quantitative trait loci (QTLs) for Fusarium head blight (FHB) resistance and for controlling accumulation of the mycotoxin deoxynivalenol (DON) are closely located. Utilizing markers known to reside near the FHB resistance QTL and data from several wheat genetic maps, we have limited the syntenous region to 6.8 Mb of the short arm of rice chromosome 4. From the 6.8-Mb sequence of rice chromosome 4, we found three putative rice genes that could have a role in detoxification of mycotoxins. DNA sequences of these putative rice genes were used in BLAST searches to identify wheat expressed sequence tags (ESTs) exhibiting significant similarity. Combined data from expression analysis and gene mapping of wheat homologues and results of analysis of DON accumulation using doubled haploid populations revealed that a putative gene for multidrug resistance-associated protein (MRP) is a possible candidate for the FHB resistance and/or DON accumulation controlling QTLs on wheat chromosome 2DS and can be used as a molecular marker to eliminate the susceptible allele when the Chinese wheat variety Sumai 3 is used as a resistance source. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genome‐wide association study of callus induction variation to explore the callus formation mechanism of rice

Rice(Oryza sativa) is one of the most widely cultivated food crops, worldwide. Tissue culture is extensively used in rice breeding and functional genome research. The ability to induce callus determines whether a particular rice variety can be subjected to tissue culture and Agrobacterium-mediated transformation. Over the past two decades, many quantitative trait loci(QTLs)related to callus induction traits have been identified;however, individual genes associated with rice callus induction have not been reported. In this study, we characterized three callus-induction traits in a global collection of 510 rice accessions. A genome-wide association study of the rice population in its entirety as well as subpopulations revealed 21 significant loci located in rice callus induction QTLs. We identified three candidate callus induction genes, namely CRL1, Os BMM1, and Os SET1, which Rese are orthologs of Arabidopsis LBD17/LBD29, BBM, and SWN,respectively, which are known to affect callus formation.Furthermore, we predicted that 14 candidate genes might be involved in rice callus induction and showed that RNA interference(RNAi)-mediated disruption of Os IAA10 inhibited callus formation on tissue culture medium.Embryo growth in the Os IAA10 RNAi line was not inhibited by synthetic auxin(2,4-D) treatment, suggesting that Os IAA10 may perceive auxin and activate the expression of downstream genes, such as CRL1, to induce callus formation. The significant loci and candidate genes identified here may provide insight into the mechanism underlying callus formation in rice.     

Genome-wide scan identifies loci associated with classical BSE occurrence

Classical bovine spongiform encephalopathy (BSE) is an acquired prion disease that is invariably fatal in cattle and has been implicated as a significant human health risk. Sequence variations in the coding region of the prion gene (PRNP) have been associated with acquired transmissible spongiform encephalopathy (TSE) susceptibility in mammals; however, this is not the case in cattle. It has been hypothesized that genes, in addition to the prion gene, contribute to genetic susceptibility of acquired TSEs. Accordingly, genetic studies of classical BSE in cattle identified loci other than PRNP that are associated with disease incidence. The objective of this study was to utilize a genome-wide association study to test for genetic loci associated with classical BSE. The samples include 143 BSE affected (case) and 173 unaffected half sib (control) animals collected in the mid 1990s in Southern England. The data analysis identifies loci on two different chromosomes associated with BSE disease occurrence. Most notable is a single nucleotide polymorphism on chromosome 1 at 29.15 Mb that is associated with BSE disease (p = 3.09E-05). Additionally, a locus on chromosome 14, within a cluster of SNPs showed a trend toward significance (p = 5.24E-05). It is worth noting that in a human vCJD study markers on human chromosome 8, a region with shared synteny to the region identified on cattle chromosome 14, were associated with disease. Further, our candidate genes appear to have plausible biological relevance with the known etiology of TSE disease. One of the candidate genes is hypothetical gene LOC521010, similar to FK506 binding protein 2 located on chromosome 1 at 29.32 Mb. This gene encodes a protein that is a member of the immunophilin protein family and is involved in basic cellular processes including protein folding. The chromosomal regions identified in this study and candidate genes within these regions merit further investigation.     

Identification of the duplicated segments in rice chromosomes 1 and 5 by linkage analysis of cDNA markers of known functions

We mapped two loci for ADP-ribosylation factor homologues (ARF1, ARF2) and two loci for cysteine proteinase inhibitors (oryzacystatin-I and -II: OCI, OCII) by linkage analysis of restriction fragment length polymorphism loci in rice (Oryza sativa L.) genomic DNAs using their cDNAs as probes.Oc-1 andArf-2 were found to be closely located to each other on chromosome 1, whileOc-2 andArf-1,both found on chromosome 5, were also located close to each other. The map distances are about 2 cM in both pairs. In each chromosome, theArf locus was located about 27 cM from that of the aldolase gene (Ald-2 in chromosome 1 andAld-1 in chromosome 5). These three genes are in the same order,Ald-Arf-Oc, but in opposite orientations relative to the distal ends of the linkage group. The presence of two sets of three linked genes on chromosomes 1 and 5 strongly suggests a structural similarity of the blocks of the two chromosomes, which probably reflects duplication of the segment. A recent investigation by other workers has shown that these rice blocks correspond to two regions in maize chromosomes 8 and 6, that have previously been shown to share many duplicated nucleotide sequences. It is therefore very likely that the duplication of the region occurred before the divergence of rice and maize during the evolution of the subfamilies of the grasses (Gramineae). In view of a recently discovered possible structural similarity between the small GTP-binding protein superfamily, which includesArf andras proteins, and the cystatin family, the close linkage ofOc andArf loci found in the present study suggests a possible cluster of genes related to the small GTP-binding proteins.