生物信息学分析方法I: 全基因组关联分析概述

全基因组关联分析(GWAS)是动植物复杂性状相关基因定位的常用手段。高通量基因分型技术的应用极大地推动了GWAS的发展。在植物中, 利用GWAS不仅能够以较高的分辨率在全基因组水平鉴定出各种自然群体特定性状相关的基因或区间, 而且可揭示表型变异的遗传架构全景图。目前, 人们利用GWAS分析方法已在拟南芥(Arabidopsis thaliana)、水稻(Oryza sativa)、小麦(Triticum aestivum)、玉米(Zea mays)和大豆(Glycine max)等模式植物和重要农作物品系中发掘出与各种性状显著相关的数量性状座位(QTL)及其候选基因位点, 阐明了这些性状的遗传基础, 并为揭示这些性状背后的分子机理提供候选基因, 也为作物高产优质品种的选育提供了理论依据。该文对GWAS的方法、影响因素及数据分析流程进行了详细描述, 以期为相关研究提供参考。     

360度群体遗传变异扫描——大豆泛基因组研究

大豆(Glycine max)是重要的油料和蛋白作物, 其丰富的遗传变异为生物学性状挖掘和育种改良提供了重要的资源基础。然而, 单个基因组信息无法全面揭示种质资源的遗传变异, 泛基因组研究为解决这一不足提供了新方案。近日, 中国科学院遗传与发育生物学研究所田志喜和梁承志研究团队从2 898份大豆种质中选取26份代表性材料, 并整合已有的3个基因组, 构建了包含野生和栽培大豆的泛基因组和图基因组(graph-based genome), 鉴定了整个群体的绝大多数结构变异数据集, 确定了大豆种质的核心、非必需和个体特异的基因集。利用这些数据系统地揭示了生育期位点E3的等位基因变异和基因融合事件、种皮颜色基因I的单体型和演化关系以及结构变异对铁离子转运基因表达和地区适应性选择的影响。该研究为作物基因组学研究提供了一个新的模式, 同时将加速推动大豆遗传变异的鉴定、性状解析和种质创新。     

是金子无论在何处都发光: 玉米和水稻驯化中的趋同选择

野生植物的驯化为人类定居与文明起源奠定了重要基础。在世界范围内不同地区生活的古人类分别对当地不同的野生植物进行了驯化, 而经过驯化的作物常常表现出相似的驯化综合性状。在基因组层面上对趋同选择规律的解析, 可为作物育种提供重要信息与遗传资源。近日, 中国农业大学杨小红/李建生和华中农业大学严建兵领衔的团队从单基因和全基因组2个层次系统解析了玉米(Zea mays)和水稻(Oryza sativa)趋同选择的遗传基础, 发现玉米KRN2与水稻OsKRN2受到了趋同选择, 并通过相似的途径调控玉米与水稻的粒数与产量。他们还发现玉米与水稻在全基因组范围内存在大量趋同选择同源基因对(gene pair), 这些基因在淀粉代谢、糖及辅酶合成等途径特异富集。该研究不仅克隆了在玉米与水稻中均具有重要育种价值的趋同选择同源基因对KRN2/OsKRN2, 而且在全基因组水平上揭示了玉米与水稻趋同选择的规律, 为进一步解析驯化综合性状形成的分子机理及其在育种中的应用奠定了重要理论基础。     

鱼和熊掌的选择: 反向重复序列变异介导的玉米环境适应与产量平衡

作物育种的目标是找到产量和抗性的最佳平衡点, 其中涉及“鱼和熊掌”二者兼得的选择策略。哪些逆境负调控位点影响产量性状, 以及如何调控等是突破育种瓶颈的重要科学问题。近百年来, 高产玉米(Zea mays)育种使玉米单产不断提高, 同时现代玉米品种对干旱的敏感性也呈现出增强趋势, 故而存在高产稳产的潜在风险。可对于这一现象背后确切的遗传机制却知之甚少, 从而限制了既高产又高抗玉米新品种的培育。玉米的非生物胁迫抗性与产量性状均为多基因控制的复杂数量性状, 涉及全基因组范围内大量基因的表达与调控。玉米基因组内存在大量的小RNA (sRNA), 其对基因表达起精细调控作用, 但人们对sRNA调控作物环境胁迫应答与产量性状机制的理解仍然有限。近日, 华中农业大学代明球课题组与李林和李峰两个课题组合作, 基于对338份玉米关联群体在不同环境下的sRNA表达组分析, 鉴定到大量干旱应答的sRNA, 以及调控这些sRNA表达的遗传位点(eQTL); 并克隆了8号染色体上1个干旱特异性eQTL热点DRESH8。生物信息学分析显示, DRESH8是1个由转座子组成的长度约为21.4 kb的反向重复序列(TE-IR)。DRESH8通过产生小干扰RNA (siRNA)介导抗旱基因的转录后沉默, 并间接抑制产量负调控因子的表达, 在负调控干旱应答的同时正调控产量性状。进一步研究发现, DRESH8在玉米驯化和改良过程中受到了人工选择。据此, 他们认为DRESH8可能是玉米平衡抗旱性和产量的关键遗传位点。该研究在全基因组水平上揭示了作物调控产量和环境胁迫抗性平衡的关键遗传机制, 同时也鉴定到大量IR位点, 为未来“高抗、高产”玉米设计育种提供了有价值的操控靶点。     

水稻OsJMJ718基因可选择性多聚腺苷酸化序列的 克隆及生殖发育期表达模式

可选择性多聚腺苷酸化是真核生物重要的基因调控机制之一, 通过形成不同长度的3'端非翻译区影响信使RNA的稳定性、定位和翻译效率, 从而增加转录本的复杂度。已有研究表明, 拟南芥(Arabidopsis thaliana)中参与DNA去甲基化调控IBM1基因的可选择性多聚腺苷酸化加工受染色质调节因子EDM2调控, 从而影响拟南芥基因组数千基因的CHG甲基化水平, 但该类调控机制是否在其它物种中同样存在仍然未知。以水稻(Oryza sativa)基因组中IBM1同源基因OsJMJ718为研究对象, 利用生物信息学分析和3'RACE实验, 发现IBM1同源基因也存在可选择性多聚腺苷酸化修饰, OsJMJ718基因可能存在9个可选择性多聚腺苷酸化序列。序列比对分析表明, NCBI网站现存日本晴OsJMJ718基因组3'末端序列与9522和明辉63等其它生态型基因组序列组成可能不同。荧光实时定量PCR分析表明, OsJMJ718的9个转录本在水稻生殖发育阶段呈现不同的动态表达模式, 其中TVX5转录本表达量最高。研究获得的OsJMJ718基因可选择性多聚腺苷酸化序列信息及相关的表达模式分析为进一步揭示水稻OsJMJ718基因的可选择性多聚腺苷酸化分子机制和生物学功能奠定了基础。     

全基因组关联分析实现水稻粒型自然变异的分子解析

全基因组关联分析(GWAS)近年来被广泛应用于解析生物自然变异的遗传基础。但限于其遗传定位精度, 在水稻(Oryza sativa)遗传学研究中, 该方法尚无法取代传统的图位克隆法在克隆复杂性状调控基因中的作用。近期, 中国科学家在应用GWAS等大数据来克隆控制水稻粒长和粒重等复杂性状的QTL方面取得了新突破。     

野生大豆耐逆分子调控机制研究进展

野生大豆(Glycine soja)起源于中国, 是栽培大豆(G. max)的近缘祖先, 逆境适应能力强, 是研究耐逆分子机制和挖掘耐逆关键调控基因的优良材料。该文综述了野生大豆耐逆基因组、转录组和蛋白质组等组学研究进展, 总结了近年来类受体蛋白激酶、转录因子、离子通道和氧化还原在野生大豆耐逆应答中的调控作用及机制, 为耐逆作物新品种培育提供了新思路。     

敲降VPS28基因对中国荷斯坦奶牛乳脂合成的调控

本课题组前期通过GWAS研究,发现VPS28基因在荷斯坦奶牛乳腺组织中特异性高表达,且其5′-UTR的突变位点-58C>T与乳脂性状关联,但其对乳脂性状的调控机理尚未明确。本研究为了明确VPS28基因及其突变位点-58C>T对乳脂的调控机理,首先利用启动子活性分析检测突变位点-58C>T对VPS28基因的影响,发现该突变位点显著降低VPS28基因启动子活性;然后利用RNA干扰技术敲降奶牛原代乳腺上皮细胞中VPS28基因表达量,检测VPS28通路和乳脂合成相关基因mRNA表达量以及细胞中脂肪滴形态,分析结果发现敲降VPS28基因可降低泛素化-溶酶体和泛素化-蛋白酶体通路基因和乳脂合成相关基因的表达量,并提高细胞中甘油三酯的合成,预示VPS28基因可能通过泛素化-溶酶体和泛素化-蛋白酶体途径调控乳脂生成。本研究结果在转录组水平揭示VPS28基因对乳脂合成的调控机制,为奶牛乳脂性状的分子育种研究提供参考依据。     

ZmFBL41 Chang7-2: 玉米抗纹枯病的关键利器

由真菌Rhizoctonia solani引起的纹枯病严重危害玉米(Zea mays)和水稻(Oryza sativa)等作物的安全生产。R. solani的宿主范围广且抗源少, 加之相关的抗性机制研究有限, 导致纹枯病的危害长期得不到有效控制。近期, 中国科学家通过对318份玉米自交系进行全基因组关联分析, 筛选到1个与纹枯病抗性相关的、编码F-box结构域蛋白的候选基因ZmFBL41 (GRMZM2G109140)。ZmFBL41蛋白是SCF (SKP1-Cullin-F-box) E3泛素连接酶复合体的一员, 能介导复合体对肉桂醇脱氢酶ZmCAD的降解, 从而降低木质素的积累, 使玉米易感纹枯病。玉米抗病自交系Chang7-2中, 蛋白ZmFBL41 ^Chang7-2因2个关键氨基酸的变异, 不能结合并降解底物ZmCAD, 使木质素含量增加, 从而提高玉米对纹枯病的抗性。该研究率先揭示了SCF复合体可通过降解肉桂醇脱氢酶来调控植物免疫反应的新型分子机制, 为提高玉米及其它作物对纹枯病的抗性提供了重要理论依据和基因资源。     

ZmFBL41 Chang7-2: 玉米抗纹枯病的关键利器

由真菌Rhizoctonia solani引起的纹枯病严重危害玉米(Zea mays)和水稻(Oryza sativa)等作物的安全生产。R. solani的宿主范围广且抗源少, 加之相关的抗性机制研究有限, 导致纹枯病的危害长期得不到有效控制。近期, 中国科学家通过对318份玉米自交系进行全基因组关联分析, 筛选到1个与纹枯病抗性相关的、编码F-box结构域蛋白的候选基因ZmFBL41 (GRMZM2G109140)。ZmFBL41蛋白是SCF (SKP1-Cullin-F-box) E3泛素连接酶复合体的一员, 能介导复合体对肉桂醇脱氢酶ZmCAD的降解, 从而降低木质素的积累, 使玉米易感纹枯病。玉米抗病自交系Chang7-2中, 蛋白ZmFBL41 ^Chang7-2因2个关键氨基酸的变异, 不能结合并降解底物ZmCAD, 使木质素含量增加, 从而提高玉米对纹枯病的抗性。该研究率先揭示了SCF复合体可通过降解肉桂醇脱氢酶来调控植物免疫反应的新型分子机制, 为提高玉米及其它作物对纹枯病的抗性提供了重要理论依据和基因资源。     

Development of soybean aphid genomic SSR markers using next generation sequencing

Simple sequence repeats (SSRs) or microsatellites are very useful molecular markers, owing to their locus-specific codominant and multiallelic nature, high abundance in the genome, and high rates of transferability across species. The soybean aphid (Aphis glycines Matsumura) has become the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America, since it was first found in the Midwest of the United States in 2000. Biotypes of the soybean aphid capable of colonizing newly developed aphid-resistant soybean cultivars have been recently discovered. Genetic resources, including molecular markers, to study soybean aphids are severely lacking. Recently developed next generation sequencing platforms offer opportunities for high-throughput and inexpensive genome sequencing and rapid marker development. The objectives of this study were (i) to develop and characterize genomic SSR markers from soybean aphid genomic sequences generated by next generation sequencing technology and (ii) to evaluate the utility of the SSRs for genetic diversity or relationship analyses. In total 128 SSR primer pairs were designed from sequences generated by Illumina GAII from a reduced representation library of A. glycines. Nearly 94% (120) of the primer pairs amplified SSR alleles of expected size and 24 SSR loci were polymorphic among three aphid samples from three populations. The polymorphic SSRs were successfully used to differentiate among 24 soybean aphids from Ohio and South Dakota. Sequencing of PCR products of two SSR markers from four aphid samples revealed that the allelic polymorphism was due to variation in the SSR repeats among the aphids. These markers should be particularly useful for genetic differentiation among aphids collected from soybean fields at different localities and regions. These SSR markers provide the soybean aphid research community with the first set of PCR-based codominant markers developed from the genomic sequences of A. glycines.     

Nuclear DNA Content Diversity in Chinese Soybean Introductions

Intraspecific nuclear DNA content has been documented in variousangiosperm species. The purpose of this study was to determinethe nuclear DNA content variation in soybean (Glycine max(L.)Merr.). Several studies have suggested that DNA content variationexists in soybean. Ninety soybean lines representing diversegeographic locations in China were analysed by flow cytometry.Nuclei were isolated and stained with either the fluorochromeDAPI or PI. After analysis, it was determined that PI stainednuclei more accurately assessed the total DNA content in soybean.A 12% variation in nuclear DNA content was observed among the90 lines. The amount of nuclear DNA in the lines was withinthe DNA range of United States cultivars previously examined.Nuclear DNA content variation in soybean is much less than thevariation reported in maize. These results could be due to thelack of polymorphism in large chromosomal elements containinglarge blocks of repetitive DNA. In addition, unlike maize, theamounts of DNA variation did not decrease as a result of moreintensive breeding in United States cultivars. Intraspecificnuclear DNA content variability is very different between thesetwo major agronomic species. Glycine max(L.) Merr.; germplasm; genome size     

大豆耐盐机理及相关基因分子标记

大豆耐盐涉及多种生理代谢途径.耐盐大豆能够通过Cl-排除、控制Na 的吸收和转运、合成渗透调节物质、改变细胞膜膜脂组分及相关酶类的活性等多种形式来适应盐胁迫;野生大豆群体具有盐腺,从形态结构上适应盐逆境;大豆-根瘤菌共生体在盐胁迫下通过互作来提高整体的耐盐性.分子生物学技术应用于大豆耐盐研究,已获得了一些与耐盐相关基因连锁的分子标记.广泛搜集筛选大豆栽培种和野生种资源,利用分子生物学技术和基因工程提高大豆耐盐性,将成为未来大豆耐盐研究的主要内容.     

大豆耐盐机理及相关基因分子标记

大豆耐盐涉及多种生理代谢途径。耐盐大豆能够通过Cl-排除、控制Na+的吸收和转运、合成渗透调节物质、改变细胞膜膜脂组分及相关酶类的活性等多种形式来适应盐胁迫;野生大豆群体具有盐腺,从形态结构上适应盐逆境;大豆-根瘤菌共生体在盐胁迫下通过互作来提高整体的耐盐性。分子生物学技术应用于大豆耐盐研究,已获得了一些与耐盐相关基因连锁的分子标记。广泛搜集筛选大豆栽培种和野生种资源,利用分子生物学技术和基因工程提高大豆耐盐性,将成为未来大豆耐盐研究的主要内容。     

Allelic Combinations of Soybean Maturity Loci E1, E2, E3 and E4 Result in Diversity of Maturity and Adaptation to Different Latitudes

Soybean cultivars are extremely diverse in time to flowering and maturation as a result of various photoperiod sensitivities. The underlying molecular genetic mechanism is not fully clear, however, four maturity loci E1, E2, E3 and E4 have been molecularly identified. In this report, cultivars were selected with various photoperiod sensitivities from different ecological zones, which covered almost all maturity groups (MG) from MG 000 to MG VIII and MG X adapted from latitude N 18° to N 53°. They were planted in the field under natural daylength condition (ND) in Beijing, China or in pots under different photoperiod treatments. Maturity-related traits were then investigated. The four E maturity loci were genotyped at the molecular level. Our results suggested that these four E genes have different impacts on maturity and their allelic variations and combinations determine the diversification of soybean maturity and adaptation to different latitudes. The genetic mechanisms underlying photoperiod sensitivity and adaptation in wild soybean seemed unique from those in cultivated soybean. The allelic combinations and functional molecular markers for the four E loci will significantly assist molecular breeding towards high productivity.     

Genetics of Soybean-Heterodera glycines Interactions

The soybean cyst nematode, Heterodera glycines, is one of the most economically important pathogens of soybean. Effective management of the nematode is often dependent on the planting of resistant soybean cultivars. During the past 40 years, more than 60 soybean genotypes and plant introductions (PI) have been reported as resistant to H. glycines. About 130 modern soybean cultivars registered in the United States are resistant to certain races of H. glycines. Several resistance genes have been identified and genetically mapped; however, resistance levels in many soybean cultivars are not durable. Some older cultivars are no longer resistant to certain H. glycines populations in many production areas, especially if a soybean monoculture has been practiced. Past soybean registration reports show that all resistant cultivars developed in public institutions from the mid-1960s to the present have been derived from five PIs. This narrow genetic background is fragile. To further complicate the issue, soybean-H. glycines genetic interactions are complex and poorly understood. Studies to identify soybean resistance genes sometimes have overlapped, and the same genes may have been reported several times and designated by different names. Nevertheless, many potential resistance genes in existing germplasm resources have not yet been characterized. Clearly, it is necessary to identify new resistance genes, develop more precise selection methods, and integrate these resistance genes into new cultivars. Rational deployment of resistant cultivars is critical to future sustained soybean production.     

A PCR-based diagnostic tool for distinguishing grape skin color mutants

Berry skin color mutants are phenotypically different from their original cultivars, but they show identical molecular profile if analysed by using microsatellite markers. This work gives an easy, inexpensive and quick diagnostic tool to discriminate these somatic variants. We distinguished some grape (Vitis vinifera L.) skin color mutants from white to red or pink and from black to grey, pink or white and we investigated their molecular bases by single-strand conformational polymorphism (SSCP), single base primer extension and coding sequence analysis of anthocyanin biosynthetic enzyme genes and by polymerase chain reaction (PCR) analysis of VvmybA1 regulatory gene. Analyses of structural genes did not reveal polymorphisms between wild type and mutant cultivars but only among different varieties, whereas the study of VvmybA1 regulatory gene has given important outcomes for color mutants characterisation. The discrimination between white wild type and its derived colored mutant and between black wild type and white mutant has been obtained through a simple test of amplification for presence/absence. The discrimination between black wild type and less colored mutant has occurred through a quantitative result on agarose gel confirmed by real-time PCR analysis: the amount of functional allele in less colored somatic variants genome was about one-fourth of the correspondent quantity in original black cultivars genome.     

Glyphosate affects photosynthesis in first and second generation of glyphosate-resistant soybeans

The crop area planted to conventional soybeans has decreased annually while that planted to glyphosate-resistant (RR) soybean has drastically increased mainly due to the wide adoption of glyphosate in current weed management systems. With the extensive use of glyphosate, many farmers have noted visual plant injury in RR soybean varieties after glyphosate application. A new generation designated as “second generation—RR2” has been recently developed and these RR2 cultivars already are commercially available for farmers and promoted as higher yielding relative to the previous RR cultivars. However, little information is currently available about the performance of RR2 soybean beyond commercial and farmer testimonial data. Thus, an evaluation of different glyphosate rates applied in different growth stages of the first and second generation of RR soybeans, revealed a significant decrease in photosynthesis. In general, increased glyphosate rate and late applications (V6) pronounced decrease photosynthetic parameters and consequently decreased in leaf area and shoot biomass production. In contrast, low rate and early applications were less damage for the RR soybean plants, suggesting that with early applications (V2), plants probably have more time to recover from glyphosate or its metabolites effects regarding late applications.