Analysis of genetic architecture and favorable allele usage of agronomic traits in a large collection of Chinese rice accessions

Genotyping and phenotyping large natural populations provide opportunities for population genomic analysis and genome-wide association studies (GWAS). Several rice populations have been re-sequenced in the past decade; however, many major Chinese rice cultivars were not included in these studies. Here, we report large-scale genomic and phenotypic datasets for a collection mainly comprised of 1,275 rice accessions of widely planted cultivars and parental hybrid rice lines from China. The population was divided into three indica/Xian and three japonica/Geng phylogenetic subgroups that correlate strongly with their geographic or breeding origins. We acquired a total of 146 phenotypic datasets for 29 agronomic traits under multi-environments for different subpopulations. With GWAS, we identified a total of 143 significant association loci, including three newly identified candidate genes or alleles that control heading date or amylose content. Our genotypic analysis of agronomically important genes in the population revealed that many favorable alleles are underused in elite accessions, suggesting they may be used to provide improvements in future breeding efforts. Our study provides useful resources for rice genetics research and breeding.

     

Genome-wide association mapping of salinity tolerance in rice (Oryza sativa)

Salinity tolerance in rice is highly desirable to sustain production in areas rendered saline due to various reasons. It is a complex quantitative trait having different components, which can be dissected effectively by genome-wide association study (GWAS). Here, we implemented GWAS to identify loci controlling salinity tolerance in rice. A custom-designed array based on 6,000 single nucleotide polymorphisms (SNPs) in as many stress-responsive genes, distributed at an average physical interval of <100 kb on 12 rice chromosomes, was used to genotype 220 rice accessions using Infinium high-throughput assay. Genetic association was analysed with 12 different traits recorded on these accessions under field conditions at reproductive stage. We identified 20 SNPs (loci) significantly associated with Na⁺/K⁺ ratio, and 44 SNPs with other traits observed under stress condition. The loci identified for various salinity indices through GWAS explained 5–18% of the phenotypic variance. The region harbouring Saltol, a major quantitative trait loci (QTLs) on chromosome 1 in rice, which is known to control salinity tolerance at seedling stage, was detected as a major association with Na⁺/K⁺ ratio measured at reproductive stage in our study. In addition to Saltol, we also found GWAS peaks representing new QTLs on chromosomes 4, 6 and 7. The current association mapping panel contained mostly indica accessions that can serve as source of novel salt tolerance genes and alleles. The gene-based SNP array used in this study was found cost-effective and efficient in unveiling genomic regions/candidate genes regulating salinity stress tolerance in rice.     

Population structure and genetic basis of the agronomic traits of upland cotton in China revealed by a genome‐wide association study using high‐density SNPs

Gossypium hirsutum L. represents the largest source of textile fibre, and China is one of the largest cotton‐producing and cotton‐consuming countries in the world. To investigate the genetic architecture of the agronomic traits of upland cotton in China, a diverse and nationwide population containing 503 G. hirsutum accessions was collected for a genome‐wide association study (GWAS) on 16 agronomic traits. The accessions were planted in four places from 2012 to 2013 for phenotyping. The CottonSNP63K array and a published high‐density map based on this array were used for genotyping. The 503 G. hirsutum accessions were divided into three subpopulations based on 11 975 quantified polymorphic single‐nucleotide polymorphisms (SNPs). By comparing the genetic structure and phenotypic variation among three genetic subpopulations, seven geographic distributions and four breeding periods, we found that geographic distribution and breeding period were not the determinants of genetic structure. In addition, no obvious phenotypic differentiations were found among the three subpopulations, even though they had different genetic backgrounds. A total of 324 SNPs and 160 candidate quantitative trait loci (QTL) regions were identified as significantly associated with the 16 agronomic traits. A network was established for multieffects in QTLs and interassociations among traits. Thirty‐eight associated regions had pleiotropic effects controlling more than one trait. One candidate gene, Gh_D08G2376, was speculated to control the lint percentage (LP). This GWAS is the first report using high‐resolution SNPs in upland cotton in China to comprehensively investigate agronomic traits, and it provides a fundamental resource for cotton genetic research and breeding.     

A genome-wide association study uncovers novel genomic regions and candidate genes of yield-related traits in upland cotton

Key message

A total of 62 SNPs associated with yield-related traits were identified by a GWAS. Based on significant SNPs, two candidate genes pleiotropically increase lint yield.

Abstract

Improved fibre yield is considered a constant goal of upland cotton (Gossypium hirsutum) breeding worldwide, but the understanding of the genetic basis controlling yield-related traits remains limited. To better decipher the molecular mechanism underlying these traits, we conducted a genome-wide association study to determine candidate loci associated with six yield-related traits in a population of 719 upland cotton germplasm accessions; to accomplish this, we used 10,511 single-nucleotide polymorphisms (SNPs) genotyped by an Illumina CottonSNP63K array. Six traits, including the boll number, boll weight, lint percentage, fruit branch number, seed index and lint index, were assessed in multiple environments; large variation in all phenotypes was detected across accessions. We identified 62 SNP loci that were significantly associated with different traits on chromosomes A07, D03, D05, D09, D10 and D12. A total of 689 candidate genes were screened, and 27 of them contained at least one significant SNP. Furthermore, two genes (Gh_D03G1064 and Gh_D12G2354) that pleiotropically increase lint yield were identified. These identified SNPs and candidate genes provide important insights into the genetic control underlying high yields in G. hirsutum, ultimately facilitating breeding programmes of high-yielding cotton.

     

Genome-wide identification and analysis of rice genes to elucidate morphological agronomic traits

Molecular understanding of morphological agronomic traits is very important to improve grain yield and quality. According to the literature information summarized in Overview of Functionally Characterized Genes in Rice online database, 430 genes related to these traits have been functionally characterized in rice, while the functions of other genes remain to be elucidated. Gene indexed mutants are available for at least half of the genes identified in the rice genome, and are very useful resources to study gene function. To suggest candidate genes for functional studies associated with morphological agronomic traits, we identified genes with tissue/organ-preferred expression patterns through meta-analysis of microarray data, and identified 781 genes for roots, 1,084 for leaves, 1,029 for calluses, 927 for anthers, 241 for embryos, and 343 for endosperms. Additionally, 4,243 genes expressed in all tissue types were allocated to a ubiquitously-expressed gene group (‘housekeeping’ genes). The estimated tissue/organ-preferred and housekeeping genes accounted for 40% of the characterized genes associated with morphological agronomic traits, indicating that identification of tissue/organ-preferred genes is an effective way to provide putative gene function. In this study, we reported the information of gene-indexed mutants for 84% of the identified candidate genes. Our candidate genes and relating indexed mutant resources can potentially be used to improve morphological agronomic traits in rice.     

银川番茄斑萎病毒的分子鉴定

为明确银川番茄(Lycopersicon esculentum)是否遭受了番茄斑萎病毒(TSWV)的危害, 采用国家标准TSWV RT- PCR检测技术对银川番茄上采集的14份疑似感染TSWV病叶样本进行分子鉴定, 对克隆得到的核衣壳蛋白基因N (Nucleocapsid)序列进行多序列比对和系统进化树分析, 随后对PCR阳性样本进行蛋白检测。结果表明, 14份病叶样本中有8份扩增出长度为394 bp的TSWV N基因序列, 且8条序列完全一致; 获得的银川番茄TSWV分离物与云南番茄、中国莴苣(Lactuca sativa)、中国鸢尾(Iris tectorum)和重庆辣椒(Capsicum annuum) TSWV分离物相对近缘, 与山东、黑龙江和北京等地及国外TSWV分离物相对远缘; 利用TSWV的抗体通过Western blot对8个PCR阳性样本进一步检测, 结果也证实8个阳性样本中存在TSWV感染。该研究首次通过分子鉴定及蛋白检测证明银川番茄上存在TSWV感染, 需要加快抗TSWV番茄品种的选育工作。     

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

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

Resequencing 50 accessions of cultivated and wild rice yields markers for identifying agronomically important genes

Rice is a staple crop that has undergone substantial phenotypic and physiological changes during domestication. Here we resequenced the genomes of 40 cultivated accessions selected from the major groups of rice and 10 accessions of their wild progenitors (Oryza rufipogon and Oryza nivara) to >15 × raw data coverage. We investigated genome-wide variation patterns in rice and obtained 6.5 million high-quality single nucleotide polymorphisms (SNPs) after excluding sites with missing data in any accession. Using these population SNP data, we identified thousands of genes with significantly lower diversity in cultivated but not wild rice, which represent candidate regions selected during domestication. Some of these variants are associated with important biological features, whereas others have yet to be functionally characterized. The molecular markers we have identified should be valuable for breeding and for identifying agronomically important genes in rice.     

QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum

Capsaicinoids are unique compounds produced only in peppers (Capsicum spp.). Several studies using classical quantitative trait loci (QTLs) mapping and genomewide association studies (GWAS) have identified QTLs controlling capsaicinoid content in peppers; however, neither the QTLs common to each population nor the candidate genes underlying them have been identified due to the limitations of each approach used. Here, we performed QTL mapping and GWAS for capsaicinoid content in peppers using two recombinant inbred line (RIL) populations and one GWAS population. Whole‐genome resequencing and genotyping by sequencing (GBS) were used to construct high‐density single nucleotide polymorphism (SNP) maps. Five QTL regions on chromosomes 1, 2, 3, 4 and 10 were commonly identified in both RIL populations over multiple locations and years. Furthermore, a total of 109 610 SNPs derived from two GBS libraries were used to analyse the GWAS population consisting of 208 C. annuum‐clade accessions. A total of 69 QTL regions were identified from the GWAS, 10 of which were co‐located with the QTLs identified from the two biparental populations. Within these regions, we were able to identify five candidate genes known to be involved in capsaicinoid biosynthesis. Our results demonstrate that QTL mapping and GBS‐GWAS represent a powerful combined approach for the identification of loci controlling complex traits.     

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

体重性状是肉鸡重要的经济性状。为了寻找可用于京海黄鸡体重性状遗传改良的分子标记及候选基因,本文以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位点,该区域可能是影响京海黄鸡中后期生长体重的重要候选区域。     

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

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

Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton

Key message

Thirty significant associations between 22 SNPs and five plant architecture component traits in Chinese upland cotton were identified via GWAS. Four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits. A candidate gene, Gh_D03G0922, might be responsible for plant height in upland cotton.

Abstract

A compact plant architecture is increasingly required for mechanized harvesting processes in China. Therefore, cotton plant architecture is an important trait, and its components, such as plant height, fruit branch length and fruit branch angle, affect the suitability of a cultivar for mechanized harvesting. To determine the genetic basis of cotton plant architecture, a genome-wide association study (GWAS) was performed using a panel composed of 355 accessions and 93,250 single nucleotide polymorphisms (SNPs) identified using the specific-locus amplified fragment sequencing method. Thirty significant associations between 22 SNPs and five plant architecture component traits were identified via GWAS. Most importantly, four peak SNP loci located on chromosome D03 were simultaneously associated with more plant architecture component traits, and these SNPs were harbored in one linkage disequilibrium block. Furthermore, 21 candidate genes for plant architecture were predicted in a 0.95-Mb region including the four peak SNPs. One of these genes (Gh_D03G0922) was near the significant SNP D03_31584163 (8.40 kb), and its Arabidopsis homologs contain MADS-box domains that might be involved in plant growth and development. qRT-PCR showed that the expression of Gh_D03G0922 was upregulated in the apical buds and young leaves of the short and compact cotton varieties, and virus-induced gene silencing (VIGS) proved that the silenced plants exhibited increased PH. These results indicate that Gh_D03G0922 is likely the candidate gene for PH in cotton. The genetic variations and candidate genes identified in this study lay a foundation for cultivating moderately short and compact varieties in future Chinese cotton-breeding programs.

     

Genome-Wide Association Studies Using Single Nucleotide Polymorphism Markers Developed by Re-Sequencing of the Genomes of Cultivated Tomato

With the aim of understanding relationship between genetic and phenotypic variations in cultivated tomato, single nucleotide polymorphism (SNP) markers covering the whole genome of cultivated tomato were developed and genome-wide association studies (GWAS) were performed. The whole genomes of six tomato lines were sequenced with the ABI-5500xl SOLiD sequencer. Sequence reads covering ∼13.7× of the genome for each line were obtained, and mapped onto tomato reference genomes (SL2.40) to detect ∼1.5 million SNP candidates. Of the identified SNPs, 1.5% were considered to confer gene functions. In the subsequent Illumina GoldenGate assay for 1536 SNPs, 1293 SNPs were successfully genotyped, and 1248 showed polymorphisms among 663 tomato accessions. The whole-genome linkage disequilibrium (LD) analysis detected highly biased LD decays between euchromatic (58 kb) and heterochromatic regions (13.8 Mb). Subsequent GWAS identified SNPs that were significantly associated with agronomical traits, with SNP loci located near genes that were previously reported as candidates for these traits. This study demonstrates that attractive loci can be identified by performing GWAS with a large number of SNPs obtained from re-sequencing analysis.     

SNPs associated with body weight and backfat thickness in two pig breeds identified by a genome-wide association study

Growth and fat deposition are important economic traits due to the influence on production in pigs. In this study, a dataset of 1200 pigs with 345,570 SNPs genotyped by sequencing (GBS) was used to conduct a GWAS with single-marker regression method to identify SNPs associated with body weight and backfat thickness (BFT) and to search for candidate genes in Landrace and Yorkshire pigs. A total of 27 and 13 significant SNPs were associated with body weight and BFT, respectively. In the region of 149.85–149.89 Mb on SSC6, the SNP (SSC6: 149876737) for body weight and the SNP (SSC6: 149876507) for BFT were in the same locus region (a gap of 230 bp). Two SNPs were located in the DOCK7 gene, which is a protein-coding gene that plays an important role in pigmentation. Two SNPs located on SSC8: 54567459 and SSC11: 33043081 were found to overlap weight and BFT; however, no candidate gene was found in these regions. In addition, based on other significant SNPs, two positional candidate genes, NSRP1 and CADPS, were proposed to influence weight. In conclusion, this is the first study report using GBS data to identify the significant SNPs for weight and BFT. A total of four particularly interesting SNPs and one potential candidate genes (DOCK7) were found for these traits in domestic pigs. This study improves our knowledge to better understand the complex genetic architecture of weight and BFT, but further validation studies of these candidate loci and genes are recommended in pigs     

Genome-wide association mapping identifies markers associated with cane yield components and sucrose traits in the Louisiana sugarcane core collection

Sugarcane is an economically important crop for both food and biofuel industries. Marker-assisted breeding in sugarcane is becoming a reality with the recent development and deployment of markers linked with disease resistance genes. Large linkage disequilibrium in sugarcane makes genome-wide association studies (GWAS) a better alternative to biparental mapping to identify markers associated with agronomic traits. GWAS was conducted on a Louisiana core collection to identify marker-trait associations (MTA) for 11 cane yield and sucrose traits using single nucleotide polymorphism (SNP) and insertion-deletion (Indel) markers. Significant (P < .05) MTAs were identified for all traits where the top ranked markers explained up to 15% of the total phenotypic variation. High correlations (0.732 to 0.999) were observed between sucrose traits and 56 markers were found consistent across multiple traits. These markers following validation in more diverse populations could be used in marker-assisted selection of clones in sugarcane breeding program in Louisiana and elsewhere.     

Association between nonsynonymous mutations of starch synthase IIa and starch quality in rice (Oryza sativa)

Starch quality is one of the most important agronomic traits in Asian rice, Oryza sativa. Starch synthase IIa (SsIIa) is a major candidate gene for starch quality variation. Within SsIIa, three nonsynonymous mutations in exon 8 have been shown to affect enzyme activity when expressed in Escherichia coli. To search for the variation in SsIIa that is responsible for starch quality variation in rice, we sequenced the SsIIa exon 8 region and measured starch quality as starch disintegration in alkali for 289 accessions of cultivated rice and 57 accessions of its wild ancestor, Oryza rufipogon. A general linear model and nested clade analysis were used to identify the associations between the three nonsynonymous single nucleotide polymorphisms (SNPs) and starch quality. Among the three nonsynonymous SNPs, we found strong evidence of association at one nucleotide site ('SNP 3'), corresponding to a Leu/Phe replacement at codon 781. A second SNP, corresponding to a Val/Met replacement at codon 737, could potentially show an association with increased sample sizes. Variation in SsIIa enzyme activity is associated with the cohesiveness of rice grains when cooked, and our findings are consistent with selection for more cohesive grains during the domestication of tropical japonica rice.     

小麦芒基因定位及其与农艺性状的相关性分析

芒是位于植物穗上的针状结构,广泛存在于禾本科作物水稻、小麦、高粱和大麦中,不同作物芒的结构存在差异。小麦中,芒对提高穗光合效率和产量、防鸟、抗虫及抗逆有重要作用。前人已经对抑制小麦芒发育的主要基因进行了定位和遗传分析,4个主效基因中仅有B1(Tipped1)基因被克隆。本研究基于人工群体云南3号和偃展1号BC_3F_6群体(YN3/YZ1)和自然群体,分析了芒与其他农艺性状的关系,发现芒对株高和产量性状有显著影响;用小麦660K SNP芯片扫描YN3/YZ1和自然群体,全基因组关联分析(GWAS)显示小麦染色体5AL和6BL存在与芒性状显著相关的基因组区域,分别对应于小麦芒抑制基因B1和B2;长芒和顶芒近等基因系转录组分析发现,在6BL候选区间内有23个差异表达基因。本研究将为进一步克隆B2基因提供参考。     

Association Analysis of the Growth of Black Poplar (<i>Populus Nigra</i> L.) Under Contrasting Nitrogen Levels

The European black poplar (Populus nigra L.) has been used as a germplasm resource for the breeding of new poplar varieties around the world. The identification and screening of its high nitrogen use efficiency genotypes could enable the breeding of new resource-efficient poplar varieties. The accessions were screened using MALDI-TOF MS genotyping technology for ammonium transporter (AMT) and nitrate transporters (NRT) genes against phenotypic data for seedling height and ground diameter traits, in both low and high nitrogen environments. Allele re-sequencing of seven genes related to root development was carried out using the minisequencing method. By cluster analysis, 101 accessions of black poplar were divided into 4 populations, and it was concluded that Central Europe is the origin of the evolution of low-nitrogen and highefficiency populations of European black poplar. Association study between SNP typing and seedling height and ground diameter traits showed that there were significant correlations between four SNP loci and growth traits under the contrasting N levels. We found that SNP3 and SNP4 in the PttAMT1;3 gene were significantly associated with seedling height traits, and that SNP2 and SNP7 in the PttAMT1;2 and PttAMT1;5 genes, respectively, were significantly associated with ground diameter traits. Thus, considerable allelic diversity is present within the candidate genes studied and can be utilized to develop functional markers to select for poplars with improved growth under N stress conditions.     

Genomic consequences of selection and genome-wide association mapping in soybean

Background

Crop improvement always involves selection of specific alleles at genes controlling traits of agronomic importance, likely resulting in detectable signatures of selection within the genome of modern soybean (Glycine max L. Merr.). The identification of these signatures of selection is meaningful from the perspective of evolutionary biology and for uncovering the genetic architecture of agronomic traits.

Results

To this end, two populations of soybean, consisting of 342 landraces and 1062 improved lines, were genotyped with the SoySNP50K Illumina BeadChip containing 52,041 single nucleotide polymorphisms (SNPs), and systematically phenotyped for 9 agronomic traits. A cross-population composite likelihood ratio (XP-CLR) method was used to screen the signals of selective sweeps. A total of 125 candidate selection regions were identified, many of which harbored genes potentially involved in crop improvement. To further investigate whether these candidate regions were in fact enriched for genes affected by selection, genome-wide association studies (GWAS) were conducted on 7 selection traits targeted in soybean breeding (grain yield, plant height, lodging, maturity date, seed coat color, seed protein and oil content) and 2 non-selection traits (pubescence and flower color). Major genomic regions associated with selection traits overlapped with candidate selection regions, whereas no overlap of this kind occurred for the non-selection traits, suggesting that the selection sweeps identified are associated with traits of agronomic importance. Multiple novel loci and refined map locations of known loci related to these traits were also identified.

Conclusions

These findings illustrate that comparative genomic analyses, especially when combined with GWAS, are a promising approach to dissect the genetic architecture of complex traits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1872-y) contains supplementary material, which is available to authorized users.