甘蓝型油菜每角果粒数全基因组关联分析

为挖掘甘蓝型油菜每角果粒数显著关联单核苷酸多态性（SNP）位点及相关候选基因。本研究以300份甘蓝型油菜自交系为试验材料,对甘蓝油菜每角果粒数进行一年两地表型考察,并结合该群体前期开发的201 817个SNPs标记,采用一般线性模型（GLM）和混合线性模型（MLM）进行全基因组关联分析（GWAS）,此外,对性状显著关联SNP位点两侧100 kb区域内相关候选基因进行功能预测。300份甘蓝型油菜每角果粒数在两地均表现出广泛的表型变异,筛选出2份每角果粒数较多的油菜种质资源。基于GLM模型检测到39个与油菜每角果粒数显著关联SNPs,采用MLM分析发现,两地共检测到的3个每角果粒数显著关联SNPs位点均在GLM检测到。8个位点附近找到CIK,ERF022和EDE1等19个拟南芥已报道角果籽粒发育相关的同源基因。研究结果有助于解析甘蓝型油菜每角果粒数的遗传基础,为研究每角果粒数的调控机制、指导每角果粒数的遗传改良奠定基础。     

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

通过对小麦耐低磷相关性状进行全基因组关联分析(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结构域蛋白有关,是植物胁迫相关调控蛋白候选基因。     

甘蓝型油菜遗传图谱的构建及单株产量构成因素的QTL分析

采用常规品系04-1139与高产多角果品系05-1054构建的F2代群体为作图群体, 运用SSR(Simple sequence repeat)和SRAP(Sequence-related amplified polymorphism)构建分子标记遗传图谱并对甘蓝型油菜单株产量构成因素进行QTL分析。遗传图谱包含200个分子标记, 分布于19个连锁群上, 总长度1 700.23 cM, 标记间的平均距离8.50 cM。采用复合区间作图法(Composite interval mapping, CIM)对单株产量构成因素(单株有效角果数、每果粒数和千粒重)进行QTL分析, 共检测到12个QTL: 其中单株有效角果数4个QTL, 分别解释表型变异为35.64%、12.96%、28.71%和34.02%; 每果粒数获得5个QTL, 分别解释表型变异为8.41%、7.87%、24.37%、8.57%和14.31%; 千粒重获得３个QTL, 分别解释表型变异为2.33%、1.81%和1.86%。结果表明: 同一性状的等位基因增效作用可以同时来自高值亲本和低值亲本; 文章中与主效QTL连锁的标记可用于油菜产量性状的分子标记辅助选择和聚合育种。     

小麦根部耐盐性状全基因组关联分析

为了解小麦耐盐相关性状的遗传机理,挖掘与小麦耐盐性显著相关的SNP位点及候选基因,本研究利用浓度200 mmol/L的NaCl溶液和正常营养液对全国300份小麦品种(系)进行耐盐性试验,并利用小麦90 K芯片对分布于小麦全基因组的16650个SNP,采用Q+K关联混合模型对小麦最长根长、根干重、根鲜重、根平均直径、根尖数、根表面积、根体积和总根长等8个根部耐盐性相关性状进行全基因组关联分析(GWAS,genome-wide association study)。研究结果表明,小麦根部性状表现出广泛的表型变异,变异系数为24.3%~50.0%,多态性信息含量(PIC,polymorphic information content)为0.170~0.562,全基因组LD衰减距离为6 Mb;群体结构分析表明,试验所用300份小麦品种(系)可分为3个亚群,亚群1包含143个(47.67%)试验材料,主要来自河南、陕西和四川;亚群2包含74个(24.67%)试验材料,主要来自北京;亚群3包含83个(27.67%)试验材料,主要来自河南。GWAS共检测到77个与小麦耐盐相关性状显著关联的SNP位点(P≤0.001),这些位点分布在小麦除6D外的20条染色体上,单个SNP位点可解释3.70%~19.45%的表型变异,其中位于1A、3A、4A、7A、3D和5D染色体上的RAC875_c13169_459等6个位点同时关联到2个或2个以上性状,贡献率为3.78%~19.45%;对77个SNP位点进行发掘,筛选到17个可能与小麦耐盐性有关的候选基因。TraesCS5B01G031800(阳离子反转运蛋白)在Na+等阳离子转运中起重要作用,TraesCS5A01G329000(防御素)可以在阻断Na+等阳离子进入过程中起作用,TraesCS2A01G079000(重复富脯氨酸细胞壁蛋白)在细胞壁的形成中起重要作用,这些候选基因可作为耐盐性重要基因。     

玉米穗部性状及其一般配合力的关联分析

穗部性状是影响玉米产量的重要性状,一般配合力是评价玉米自交系利用价值的重要指标。为解析穗部性状及其一般配合力的遗传基础,本研究对248份玉米自交系组成的自然群体和以其中100份自交系为母本按照NCⅡ遗传交配设计与4个测验种(Mo17、昌7-2、E28和郑58)组配的400份F1杂交组合的穗部性状进行研究,并利用分布于全基因组的83057个SNP标记进行穗部性状及其一般配合力的关联分析。结果表明,穗长、穗粗2个穗部性状基因型间、环境间差异达极显著水平,其广义遗传率分别为81.22%和87.70%。母本间、父本间及不同杂交组合间穗长、穗粗差异均达极显著水平,在基因型方差中特殊配合力贡献率较大。利用2年2点4个环境下的数据分别进行关联分析,检测到34个性状SNP关联,利用BLUP值检测到7个性状SNP关联。这些性状SNP关联可解释的表型变异为0.01%～19.42%,其中有5个性状SNP关联的表型贡献率大于10%,未检测到穗部性状本身与一般配合力性状的相同SNP位点。基于该群体的LD衰减距离在显著关联SNP位点上下游各120 kb范围内进行候选基因搜索,共发现158个候选基因,推测可能的候选基因涉及泛素代谢相关基因(GRMZM2G360374、GRMZM2G049568、GRMZM2G178120),β半乳糖苷酶(GRMZM2G178106),丝氨酸苏氨酸蛋白激酶(GRMZM2G127050),赖氨酸和组氨酸特异性转运体(GRMZM2G116004)。研究结果为解析玉米穗长和穗粗及其一般配合力的遗传基础和分子辅助选择育种提供了参考。     

复杂疾病的遗传易感基因区域的精细定位

以单核苷酸多态性(Single-nucleotide polymorphism, SNP)为遗传标记, 采用全基因组关联研究(Genome-wide association studies, GWAS)的策略, 已经在660多种疾病(或性状)中发现了3800多个遗传易感基因区域。但是, 其中最显著关联的遗传变异或致病性的遗传变异位点及其生物学功能并不完全清楚。这些位点的鉴定有助于阐明复杂疾病的生物学机制, 以及发现新的疾病标记物。后GWAS时代的主要任务之一就是通过精细定位研究找到复杂疾病易感基因区域内最显著关联的易感位点或致病性的易感位点并阐明其生物学功能。针对常见变异, 可通过推断或重测序增加SNP密度, 寻找最显著关联的SNP位点, 并通过功能元件分析、表达数量性状位点(Expression quantitative trait locus, eQTL)分析和单体型分析等方法寻找功能性的SNP位点和易感基因。针对罕见变异, 则可采用重测序、罕见单体型分析、家系分析和负荷检验等方法进行精细定位。文章对这些策略和所面临的问题进行了综述。     

中国人群肝功能检测指标全基因组关联分析研究

肝功能检测(liver function test,LFTs)指标是受遗传和环境影响的复杂性状,具有个体差异性。为系统性研究中国人群全基因组范围内单核苷酸多态性(single nucleotide polymorphism,SNP)与肝功能指标之间的联系,本研究利用英国生物银行(UK Biobank)中1653名中国人的基因分型数据和表型数据为研究对象,利用PLINK软件进行全关联分析研究(genome-wide association study,GWAS),发现229个SNP与中国人群血液中的总胆红素(total bilirubin,TB)相关,27个SNP与中国人群血液中碱性磷酸酶(alkaline phosphatase,ALP)相关,36个SNP与中国人群血液中的γ-谷氨酰转肽酶(γ-glutamyl transpeptidase,GGT)相关,1个SNP与中国人群血液中的门冬氨酸氨基转移酶(aspartate transaminase,AST)相关,最显著的位点中有11个位点是新的LFTs关联位点。通过功能基因组分析,发现这些位点的临床意义(如吉尔伯特综合征),确定了候选基因(UGT1A,ABO,GGT1),为从遗传角度理解中国人群LFTs的个体差异性和肝功能指标临床精准检测提供了前期研究基础。     

羊重要性状全基因组关联分析研究进展

全基因组关联分析(genome-wide association study, GWAS)是一种复杂性状功能基因鉴定的分析策略,已成为挖掘畜禽重要经济性状候选基因的重要手段。随着绵羊和山羊基因组完成和公布,以及不同密度的SNP (single nucleotide polymorphism)芯片的推出并进行商业化推广,不仅大大丰富了羊标记辅助选择可利用的分子标记,而且还为开展重要性状的分子机理的探索提供了重要技术支撑。本文主要针对羊角、羊毛、羊奶、生长发育、肉质、繁殖和疾病等重要性状的GWAS研究所用的群体、主要研究方法和研究结果进行了综述,并对GWAS方法研究现状进行了归纳,以期为进一步利用GWAS进行羊的各种性状的遗传基础研究提供参考。     

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

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

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

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

西藏野生油菜种群间角果趋同生长格局研究

西藏地处我国西南边陲,素有"世界屋脊"和"地球第三极"之称。由于地质史独特,地形地貌复杂,气候带全,土壤类型繁多,野生植被多样,凡此种种,西藏油菜生境具有全球最典型的立体生境特色,其生态环境千差万别,堪称全球之最,独特而复杂的油菜生境,产生了丰富多样的野生白菜型油菜(Brassica Campestris)和野生芥菜型油菜(Brassica juncea)遗传资源。本文根据植物种群的构件理论,采用大样本取样方法,对西藏野生白菜型油菜和野生芥菜型油菜分枝进行了角果种群统计的生长分析。结果表明,西藏野生白菜型油菜和野生芥菜型油菜一级分枝的角果数、角果生物量、角果生物量比例均呈现Peel-Reed模型变化规律,角果宽度与角果长度之间的生长关系均呈现Yield Density模型变化规律,显示出角果种群形成与生长过程的趋同生长格局与物质分配策略。     

Major changes in the cell wall during silique development in Arabidopsis thaliana

Fruit development is a highly complex process, which involves major changes in plant metabolism leading to cell growth and differentiation. Changes in cell wall composition and structure play a major role in modulating cell growth. We investigated the changes in cell wall composition and the activities of associated enzymes during the dry fruit development of the model plant Arabidopsis thaliana. Silique development is characterized by several specific phases leading to fruit dehiscence and seed dispersal. We showed that early phases of silique growth were characterized by specific changes in non-cellulosic sugar content (rhamnose, arabinose, xylose, galactose and galacturonic acid). Xyloglucan oligosaccharide mass profiling further showed a strong increase in O-acetylated xyloglucans over the course of silique development, which could suggest a decreased capacity of xyloglucans to be associated with each other or to cellulose. The degree of methylesterification, mediated by the activity of pectin methylesterases (PMEs), decreased over the course of silique growth and dehiscence. The major changes in cell wall composition revealed by our analysis suggest that it could be major determinants in modulating cell wall rheology leading to growth or growth arrest.     

Induction of vivipary in Arabidopsis by silique culture: implications for seed dormancy and germination

Culture of excised fruits (siliques) of different ages of Arabidopsis thaliana in a solidified mineral salt medium supplemented with vitamins, myo-inositol, and 3% sucrose induces vivipary. Whereas early stage and immature embryos complete their full development before germinating viviparously in seeds enclosed in the silique, mature green embryos enclosed in green ovules germinate without further growth in culture. Vivipary is not observed in cultured siliques enclosing brown ovules with yellowish mature embryos inside. Suggestive of a role for abscisic acid in preventing vivipary on the mother plant, addition of the hormone to the culture medium is found to inhibit vivipary in cultured siliques. Although dried green ovules enclosing mature embryos require a cold treatment for germination, undried ovules of the same age do not germinate even after a cold treatment. This indicates that mature embryos enclosed in green ovules that germinate viviparously are cold resistant and have not become dormant at the time of culture of siliques. The circumvention by silique culture of a cold treatment and light exposure normally required for germination of isolated seeds of A. thaliana provides new possibilities to study the molecular biology of vivipary and seed germination in this model plant.     

Precise editing of CLAVATA genes in Brassica napus L. regulates multilocular silique development

Multilocular silique is a desirable agricultural trait with great potential for the development of high‐yield varieties of Brassica. To date, no spontaneous or induced multilocular mutants have been reported in Brassica napus, which likely reflects its allotetraploid nature and the extremely low probability of the simultaneous random mutagenesis of multiple gene copies with functional redundancy. Here, we present evidence for the efficient knockout of rapeseed homologues of CLAVATA3 (CLV3) for a secreted peptide and its related receptors CLV1 and CLV2 in the CLV signalling pathway using the CRISPR/Cas9 system and achieved stable transmission of the mutations across three generations. Each BnCLV gene has two copies located in two subgenomes. The multilocular phenotype can be recovered only in knockout mutations of both copies of each BnCLV gene, illustrating that the simultaneous alteration of multiple gene copies by CRISPR/Cas9 mutagenesis has great potential in generating agronomically important mutations in rapeseed. The mutagenesis efficiency varied widely from 0% to 48.65% in T₀ with different single‐guide RNAs (sgRNAs), indicating that the appropriate selection of the sgRNA is important for effectively generating indels in rapeseed. The double mutation of BnCLV3 produced more leaves and multilocular siliques with a significantly higher number of seeds per silique and a higher seed weight than the wild‐type and single mutant plants, potentially contributing to increased seed production. We also assessed the efficiency of the horizontal transfer of Cas9/gRNA cassettes by pollination. Our findings reveal the potential for plant breeding strategies to improve yield traits in currently cultivated rapeseed varieties.     

Comprehensive expression profiling of the pectin methylesterase gene family during silique development in Arabidopsis thaliana

Pectin methylesterases (PME, EC. 3.1.1.11) are enzymes that demethylesterify plant cell wall pectins in muro. In Arabidopsis thaliana, putative PME proteins are thought to be encoded by a 66-member gene family. This study used real-time RT-PCR to gain an overview of the expression of the entire family at eight silique developmental stages, in flower buds and in vegetative tissue in the Arabidopsis. Only 15% of the PMEs were not expressed at any of the developmental stages studied. Among expressed PMEs, expression data could be clustered into five distinct groups: 19 PMEs highly or uniquely expressed in floral buds, 4 PMEs uniquely expressed at mid-silique developmental stages, 16 PMEs highly or uniquely expressed in silique at late developmental stages, 16 PMEs mostly ubiquitously expressed, and 1 PME with a specific expression pattern, i.e. not expressed during early silique development. Comparison of expression and phylogenetic profiles showed that, within phylogenetic group 2, all but one PME belong to the floral bud expression group. Similar results were shown for a subset of one of the phylogenetic group, which differed from others by containing most of the PMEs that do not possess any PRO part next to their catalytic part. Expression data were confirmed by two promoter:GUS transgenic plant analysis revealing a PME expressed in pollen and one in young seeds. Our results highlight the high diversity of PME expression profiles. They are discussed with regard to the role of PMEs in fruit development and cell growth.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Water-soluble chlorophyll protein (WSCP) of Arabidopsis is expressed in the gynoecium and developing silique

Water-soluble chlorophyll protein (WSCP) has been found in many Brassicaceae, most often in leaves. In many cases, its expression is stress-induced, therefore, it is thought to be involved in some stress response. In this work, recombinant WSCP from Arabidopsis thaliana (AtWSCP) is found to form chlorophyll-protein complexes in vitro that share many properties with recombinant or native WSCP from Brassica oleracea, BoWSCP, including an unusual heat resistance up to 100°C in aqueous solution. A polyclonal antibody raised against the recombinant apoprotein is used to identify plant tissues expressing AtWSCP. The only plant organs containing significant amounts of AtWSCP are the gynoecium in open flowers and the septum of developing siliques, specifically the transmission tract. In fully grown but still green siliques, the protein has almost disappeared. Possible implications for AtWSCP functions are discussed.     

In vitro seed maturation in Brassica rapa L.: Relationship of silique atmosphere to storage reserve deposition

The Brassica rapa L. silique is a self-contained environment that maintains hypoxia around the developing seeds, and in which carbon dioxide accumulates to very high concentrations (>30,000 ppm). How the silique microenvironment modulates the composition and amount of storage reserves in the seeds is of interest because of the important agricultural role played by canola (B. rapa and Brassica napus) as an oilseed. Because of the small volume and dynamic nature of this microenvironment in Brassica, a standardized system was needed to study the environmental role played in storage reserve deposition. For this purpose we have developed a silique culture system that permits maturation of seed in vitro. Siliques excised from plants just 11 days after pollination complete the ripening of their seeds after 20 days of culture in light (200 μmol/m²/s) on MS medium containing 30 g/l sucrose, 0.25 mg/l BAP, and 0.025 mg/l NAA. Cytochemical localization and biochemical analyses revealed that storage reserves were affected by the in vitro maturation system. Although following a comparable ripening timeline to that occurring on the plant, and producing fully germinable seeds, in vitro maturation resulted in a 40% reduction in seed weight and the mature seeds contained decreased lipid, but increased protein, starch and soluble carbohydrates. To study the internal atmosphere surrounding the seeds, we developed a method to capture silique gases in helium with subsequent quantification of O₂ and CO₂ in the sample by gas chromatography. Analysis of the internal silique atmosphere showed that in vitro siliques provided seeds with a less oxygenated environment than they experience attached to the plant. Carbon dioxide concentrations remained high later into the maturation sequence in vitro than on the plant. When sampling gases from siliques attached to plants, we found multiple samples from the same plant resulted in higher variance than when only a single silique was sampled, suggesting that connection to the plant directly influences internal silique gases. Lower O₂ in the in vitro siliques was correlated with depressed lipid content in their mature seeds, supporting the conclusion that oxygen availability limits lipid accumulation. Previous studies showed how environmental factors influence Brassica embryos grown in tissue culture. These systems fail to preserve the component of metabolic regulation that is enforced by the silique wall tissues. Our in vitro maturation system provides a useful tool for specialized investigations since both the gaseous and hormonal environments can be readily manipulated.     

Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.)

Silique length (SL) and seed weight (SW) are two important yield-related traits controlled by quantitative trait loci (QTL) in oilseed rape (Brassica napus L.). The genetic bases underlying these two traits are largely unknown at present. In this study, we conducted QTL analyses for SL and SW using 186 recombinant inbred lines (RILs) derived from a cross between S1, an EMS mutant with extremely long siliques and large seeds, and S2, an inbred line with regular silique length and seed size. RILs were grown in Wuhan in the 2008/09 (SS09) and 2009/10 (SS10) growing seasons, and mean SL and SW for each line were investigated. Ten non-redundant QTL were identified for SL. Of these, a major QTL, cqSLA9, consistently explained as much as 53.4% of SL variation across environments. The others are minor QTL and individually explained less than 10% of the SL variation. Nine non-redundant QTL were identified for SW. Of which, one major QTL, cqSWA9, explained as much as 28.2% of the total SW variation in the SS09 and SS10 environments. In addition, three additive by additive interactions with small effects were detected for SL, and no interactions were detected for SW. Interestingly, the two major QTL, cqSLA9 for SL and cqSWA9 for SW colocalized in the same chromosomal region and were integrated into a unique QTL, uqA9. The S1 allele at this locus increases both SL and SW, suggesting that uqA9 has pleiotropic effects on both SL and SW. The existence and effect of uqA9 was confirmed in genetically different RILs derived from the cross between S1 and No2127, a resynthesized DH line having regular silique length and seed size. Individuals in one residual heterozygous line for cqSLA9 showed significant difference in silique length. The results in this study revealed that silique length in the S1 mutant is mainly controlled by the cqSLA9 locus, which will be suitable for fine mapping and marker-assisted selection in rapeseed breeding for high yield.