Prospecting sugarcane resistance to Sugarcane yellow leaf virus by genome-wide association

Key message

Using GWAS approaches, we detected independent resistant markers in sugarcane towards a vectored virus disease. Based on comparative genomics, several candidate genes potentially involved in virus/aphid/plant interactions were pinpointed.

Abstract

Yellow leaf of sugarcane is an emerging viral disease whose causal agent is a Polerovirus, the Sugarcane yellow leaf virus (SCYLV) transmitted by aphids. To identify quantitative trait loci controlling resistance to yellow leaf which are of direct relevance for breeding, we undertook a genome-wide association study (GWAS) on a sugarcane cultivar panel (n = 189) representative of current breeding germplasm. This panel was fingerprinted with 3,949 polymorphic markers (DArT and AFLP). The panel was phenotyped for SCYLV infection in leaves and stalks in two trials for two crop cycles, under natural disease pressure prevalent in Guadeloupe. Mixed linear models including co-factors representing population structure fixed effects and pairwise kinship random effects provided an efficient control of the risk of inflated type-I error at a genome-wide level. Six independent markers were significantly detected in association with SCYLV resistance phenotype. These markers explained individually between 9 and 14 % of the disease variation of the cultivar panel. Their frequency in the panel was relatively low (8–20 %). Among them, two markers were detected repeatedly across the GWAS exercises based on the different disease resistance parameters. These two markers could be blasted on Sorghum bicolor genome and candidate genes potentially involved in plant–aphid or plant–virus interactions were localized in the vicinity of sorghum homologs of sugarcane markers. Our results illustrate the potential of GWAS approaches to prospect among sugarcane germplasm for accessions likely bearing resistance alleles of significant effect useful in breeding programs.     

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

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

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

为明确银川番茄(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番茄品种的选育工作。     

甘蔗杂交品种核心种质重要农艺性状评价及亲缘关系分析

以107份甘蔗杂交品种核心种质为研究材料,ROC22为对照品种,通过新植和宿根2年的田间试验,对蔗茎公顷产量、甘蔗蔗糖分、公顷含糖量等性状进行最小显著差数法(LSD)两两比较和Duncan多重比较评价,筛选优异材料,利用SSR分子标记建立优异材料与我国骨干亲本的亲缘关系。结果表明36份材料在蔗茎公顷产量、甘蔗蔗糖分、公顷含糖量等指标上优于对照品种,可作为亲本进行杂交,且蔗茎公顷产量对公顷含糖量的影响远大于甘蔗蔗糖分对公顷含糖量的影响;20对SSR引物扩增得到292个标记,其中283个为多态性标记;优异材料与我国骨干亲本相似性系数范围在0.384~0.590之间,平均为0.437,存在较大的遗传差异;UPGMA聚类可将所有材料划分为5个类群,蔗茎公顷产量和公顷含糖量表现优异的材料在各类群都有分布,甘蔗蔗糖分表现优异的材料主要集中在Ⅱ、Ⅲ、Ⅳ类群,Ⅲ类群更集中;本研究为甘蔗种质创新、遗传育种提供了优异亲本材料,为杂交组合的配制提供重要指导。     

Detection and validation of EST-SSR markers associated with sugar-related traits in sugarcane using linkage and association mapping

Sugar-related traits are of great importance in sugarcane breeding. In the present study, quantitative trait loci (QTL) mapping validated with association mapping was used to identify expressed sequence tag-simple sequence repeats (EST-SSRs) associated with sugar-related traits. For linkage mapping, 524 EST-SSRs, 241 Amplified Fragment Length Polymorphisms, and 10 genomic SSR markers were mapped using 283 F₁ progenies derived from an interspecific cross. Six regions were identified using Multiple QTL Mapping, and 14 unlinked markers using single marker analysis. Association analysis was performed on a set of 200 accessions, based on the mixed linear model. Validation of the EST-SSR markers using association mapping within the target QTL genomic regions identified two EST-SSR markers showing a putative relationship with uridine diphosphate (UDP) glycosyltransferase, and beta-amylase, which are associated with pol and sugar yield. These functional markers can be used for marker-assisted selection of sugarcane.     

A mixed model QTL analysis for sugarcane multiple-harvest-location trial data

Sugarcane-breeding programs take at least 12 years to develop new commercial cultivars. Molecular markers offer a possibility to study the genetic architecture of quantitative traits in sugarcane, and they may be used in marker-assisted selection to speed up artificial selection. Although the performance of sugarcane progenies in breeding programs are commonly evaluated across a range of locations and harvest years, many of the QTL detection methods ignore two- and three-way interactions between QTL, harvest, and location. In this work, a strategy for QTL detection in multi-harvest-location trial data, based on interval mapping and mixed models, is proposed and applied to map QTL effects on a segregating progeny from a biparental cross of pre-commercial Brazilian cultivars, evaluated at two locations and three consecutive harvest years for cane yield (tonnes per hectare), sugar yield (tonnes per hectare), fiber percent, and sucrose content. In the mixed model, we have included appropriate (co)variance structures for modeling heterogeneity and correlation of genetic effects and non-genetic residual effects. Forty-six QTLs were found: 13 QTLs for cane yield, 14 for sugar yield, 11 for fiber percent, and 8 for sucrose content. In addition, QTL by harvest, QTL by location, and QTL by harvest by location interaction effects were significant for all evaluated traits (30 QTLs showed some interaction, and 16 none). Our results contribute to a better understanding of the genetic architecture of complex traits related to biomass production and sucrose content in sugarcane.     

Conversion of AFLP markers to high-throughput markers in a complex polyploid,sugarcane

The application of DNA markers linked to traits of commercial value in sugarcane may increase the efficiency of sugarcane breeding. The majority of markers generated for quantitative trait locus mapping in sugarcane have been single sequence repeats or AFLPs (amplified fragment length polymorphisms). Since AFLP markers are not adapted for large-scale implementation in plant breeding, our objective was to assess the feasibility of converting AFLP markers to fast, cheap and reliable PCR-based assays in a complex polyploid, sugarcane. Three AFLP markers were selected on the basis of an association to resistance to the fungal pathogen Ustilago scitaminea, the causal agent of smut in sugarcane. We developed an approach which enabled the identification of polymorphisms in these AFLP markers. Towards this goal, we employed GenomeWalking and 454 sequencing to isolate sequences adjacent to the linked AFLP markers and identify SNP (single nucleotide polymorphisms) haplotypes present in the homo(eo)logous chromosomes of sugarcane. One AFLP marker was converted to a cleavage amplified polymorphic sequence marker, another to a SCAR (sequence characteristered amplified region) marker and the final AFLP marker to a SNP PCR-based assay. However, validation of each of the markers in 240 genotypes resulted in 99, 90 and 60% correspondence with the original AFLP marker. These experiments indicate that even in a complex polyploid such as sugarcane, polymorphisms identified by AFLP can be converted to high-throughput marker systems, but due to the complexity this would only be carried out for high-value markers. In some cases, the polymorphisms identified are not transferable to more sequence-specific PCR applications.     

Marker-trait associations in Virginia Tech winter barley identified using genome-wide mapping

Genome-wide association studies (GWAS) provide an opportunity to examine the genetic architecture of quantitatively inherited traits in breeding populations. The objectives of this study were to use GWAS to identify chromosome regions governing traits of importance in six-rowed winter barley (Hordeum vulgare L.) germplasm and to identify single-nucleotide polymorphisms (SNPs) markers that can be implemented in a marker-assisted breeding program. Advanced hulled and hulless lines (329 total) were screened using 3,072 SNPs as a part of the US. Barley Coordinated Agricultural Project (CAP). Phenotypic data collected over 4 years for agronomic and food quality traits and resistance to leaf rust (caused by Puccinia hordei G. Otth), powdery mildew [caused by Blumeria graminis (DC.) E.O. Speer f. sp. hordei Em. Marchal], net blotch (caused by Pyrenophora teres), and spot blotch [caused by Cochliobolus sativus (Ito and Kuribayashi) Drechsler ex Dastur] were analyzed with SNP genotypic data in a GWAS to determine marker-trait associations. Significant SNPs associated with previously described quantitative trait loci (QTL) or genes were identified for heading date on chromosome 3H, test weight on 2H, yield on 7H, grain protein on 5H, polyphenol oxidase activity on 2H and resistance to leaf rust on 2H and 3H, powdery mildew on 1H, 2H and 4H, net blotch on 5H, and spot blotch on 7H. Novel QTL also were identified for agronomic, quality, and disease resistance traits. These SNP-trait associations provide the opportunity to directly select for QTL contributing to multiple traits in breeding programs.     

Sugarcane biotechnology: The challenges and opportunities

Summary Commercial sugarcane, belonging to the genus Saccharum (Poaceae), is an important industrial crop accounting for nearly 70% of sugar produced worldwide. Compared to other major crops, efforts to improve sugarcane are limited and relatively recent, with the first introduction of interspecific hybrids about 80 yr ago. Progress in traditional breeding of sugareane, a highly polyploid and frequently aneuploid plant, is impeded by its narrow gene pool, complex genome, poor fertility, and the long breeding/selection cycle. These constraints, however, make sugarcane a good candidate for molecular breeding. In the past decade considerable progress has been made in understanding and manipulating the sugarcane genome using various biotechnological and cell biological approaches. Notable among them are the creation of transgenic plants with improved agronomic or other important traits, advances in genomics and molecular markers, and progress in understanding the molecular aspects of sucrose transport and accumulation. More recently, substantial effort has been directed towards developing sugarcane as a biofactory for high-value products. While these achievements are commendable, a greater understanding of the sugarcane genome, and cell and whole plant physiology, will accelerate the implementation of commercially significant biotechnology outcomes. We anticipate that the rapid advancements in molecular biology and emerging biotechnology innovations would play a significant role in the future sugarcane crop improvement programs and offer many new opportunities to develop it as a new-generation industrial crop.     

Genomic Selection and Association Mapping in Rice (Oryza sativa): Effect of Trait Genetic Architecture,Training Population Composition,Marker Number and Statistical Model on Accuracy of Rice Genomic Selection in Elite,Tropical Rice Breeding Lines

Genomic Selection (GS) is a new breeding method in which genome-wide markers are used to predict the breeding value of individuals in a breeding population. GS has been shown to improve breeding efficiency in dairy cattle and several crop plant species, and here we evaluate for the first time its efficacy for breeding inbred lines of rice. We performed a genome-wide association study (GWAS) in conjunction with five-fold GS cross-validation on a population of 363 elite breeding lines from the International Rice Research Institute''s (IRRI) irrigated rice breeding program and herein report the GS results. The population was genotyped with 73,147 markers using genotyping-by-sequencing. The training population, statistical method used to build the GS model, number of markers, and trait were varied to determine their effect on prediction accuracy. For all three traits, genomic prediction models outperformed prediction based on pedigree records alone. Prediction accuracies ranged from 0.31 and 0.34 for grain yield and plant height to 0.63 for flowering time. Analyses using subsets of the full marker set suggest that using one marker every 0.2 cM is sufficient for genomic selection in this collection of rice breeding materials. RR-BLUP was the best performing statistical method for grain yield where no large effect QTL were detected by GWAS, while for flowering time, where a single very large effect QTL was detected, the non-GS multiple linear regression method outperformed GS models. For plant height, in which four mid-sized QTL were identified by GWAS, random forest produced the most consistently accurate GS models. Our results suggest that GS, informed by GWAS interpretations of genetic architecture and population structure, could become an effective tool for increasing the efficiency of rice breeding as the costs of genotyping continue to decline.     

Multiple locus genome-wide association studies for important economic traits of oil palm

Palm oil has a balanced fatty acid composition and has no trans fat. As a result, its use in food has increased as food-labeling laws have changed to specify trans fat content. Increasing oil production is the main goal in oil palm breeding. Genetic mapping and genomic studies in palm trees are necessary to understand the genetic architecture of economic traits of importance for palm oil production. To help achieve this, we sampled 422 oil palms from MPOB (Malaysian Palm Oil Board)­Angola germplasm collection and measured 13 economic traits from these palms. Multi-locus genome-wide association studies (GWAS) were conducted using least absolute shrinkage and selection operator (LASSO) and genome-wide efficient mixed model analysis. We identified 19 quantitative trait loci (QTLs) for 8 traits. Of these, four Angola-specific QTLs associated with bunch components were detected on chromosomes 4, 8, and 11. These QTLs are potentially useful for introgression of desirable genes from the Angola palms to advanced breeding populations for improvement of bunch and oil yield traits. The majority of the QTLs were detected by LASSO-A, in which the p values of individual markers were calculated based on bootstrapped standard errors. Many of the detected QTLs are nearby known QTLs detected from linkage studies reported by other research groups. We also conducted genomic selection (GS) for the 13 traits and concluded that GS can be an effective tool for oil palm breeding. This is the first GWAS and GS study conducted on oil palm germplasm from Angola, and the results can be very useful in oil palm genetic studies and breeding.     

Markers associated with stalk number and suckering in sugarcane colocate with tillering and rhizomatousness QTLs in sorghum.

Two important factors influencing sugar yield, the primary focus of sugarcane plant breeding programs, are stalk number and suckering. Molecular markers linked to both of these traits are sought to assist in the identification of high sugar yield, high stalk number, low-suckering sugarcane clones. In this preliminary mapping study, 108 progeny from a biparental cross involving two elite Australian sugarcane clones were evaluated at two sites for two years for both stalk number and suckering. A total of 258 DNA markers, including both restriction fragment length polymorphisms (RFLPs) and radio-labelled amplified fragments (RAFs), were scored and evaluated using single-factor analysis. Sixteen (7 RFLPs and 9 RAFs) and 14 (6 RFLPs and 8 RAFs) markers were identified that were significantly associated (P < 0.01) with stalk number and suckering, respectively, across both years and sites. The seven and six RFLP markers associated with stalk number and suckering, respectively, were generated by eight different RFLP probes, of which seven had been mapped in sorghum and (or) sugarcane. Of significant interest was the observation that all seven RFLP probes could be shown to be located within or near QTLs associated with tillering and rhizomatousness in sorghum. This observation highlights the usefulness of comparative mapping between sorghum and sugarcane and suggests that the identification of useful markers for stalk number and suckering in sugarcane would be facilitated by focussing on sorghum QTLs associated with related traits.     

Genome-wide association study reveals the genetic basis of growth trait in yellow catfish with sexual size dimorphism

Sexual size dimorphism has been widely observed in a large number of animals including fish species. Genome-wide association study (GWAS) is a powerful tool to dissect the genetic basis of complex traits, whereas the sex-differences in the genomics of animal complex traits have been ignored in the GWAS analysis. Yellow catfish (Pelteobagrus fulvidraco) is an important aquaculture fish in China with significant sexual size dimorphism. In this study, GWAS was conducted to identify candidate SNPs and genes related to body length (BL) and body weight (BW) in 125 female yellow catfish from a breeding population. In total, one BL-related SNP and three BW-related SNPs were identified to be significantly associated with the traits. Besides, one of these SNPs (Chr15:19195072) was shared in both the BW and BL traits in female yellow catfish, which was further validated in 185 male individuals and located on the exon of stat5b gene. Transgenic yellow catfish and zebrafish that expressed yellow catfish stat5b showed increased growth rate and reduction of sexual size dimorphism. These results not only reveal the genetic basis of growth trait and sexual size dimorphism in fish species, but also provide useful information for the marker-assisted breeding in yellow catfish.     

Genome-wide association study for oat (Avena sativa L.) beta-glucan concentration using germplasm of worldwide origin

Detection of quantitative trait loci (QTL) controlling complex traits followed by selection has become a common approach for selection in crop plants. The QTL are most often identified by linkage mapping using experimental F₂, backcross, advanced inbred, or doubled haploid families. An alternative approach for QTL detection are genome-wide association studies (GWAS) that use pre-existing lines such as those found in breeding programs. We explored the implementation of GWAS in oat (Avena sativa L.) to identify QTL affecting β-glucan concentration, a soluble dietary fiber with several human health benefits when consumed as a whole grain. A total of 431 lines of worldwide origin were tested over 2?years and genotyped using Diversity Array Technology (DArT) markers. A mixed model approach was used where both population structure fixed effects and pair-wise kinship random effects were included. Various mixed models that differed with respect to population structure and kinship were tested for their ability to control for false positives. As expected, given the level of population structure previously described in oat, population structure did not play a large role in controlling for false positives. Three independent markers were significantly associated with β-glucan concentration. Significant marker sequences were compared with rice and one of the three showed sequence homology to genes localized on rice chromosome seven adjacent to the CslF gene family, known to have β-glucan synthase function. Results indicate that GWAS in oat can be a successful option for QTL detection, more so with future development of higher-density markers.     

Genome-wide prediction models that incorporate de novo GWAS are a powerful new tool for tropical rice improvement

To address the multiple challenges to food security posed by global climate change, population growth and rising incomes, plant breeders are developing new crop varieties that can enhance both agricultural productivity and environmental sustainability. Current breeding practices, however, are unable to keep pace with demand. Genomic selection (GS) is a new technique that helps accelerate the rate of genetic gain in breeding by using whole-genome data to predict the breeding value of offspring. Here, we describe a new GS model that combines RR-BLUP with markers fit as fixed effects selected from the results of a genome-wide-association study (GWAS) on the RR-BLUP training data. We term this model GS + de novo GWAS. In a breeding population of tropical rice, GS + de novo GWAS outperformed six other models for a variety of traits and in multiple environments. On the basis of these results, we propose an extended, two-part breeding design that can be used to efficiently integrate novel variation into elite breeding populations, thus expanding genetic diversity and enhancing the potential for sustainable productivity gains.     

Identification of loci controlling forage yield and nutritive value in diploid alfalfa using GBS-GWAS

Key message

We attempted to identify genomic regions controlling forage yield and nutritive value in alfalfa. Several candidate genes and associated genetic markers were identified that could potentially be useful for alfalfa breeding to more efficiently develop improved cultivars.

Abstract

Alfalfa is one of the most widely cultivated forage legumes worldwide and improving alfalfa forage yield and nutritive value is a major global breeding goal. Genotyping-by-sequencing (GBS) provides cost-effective molecular marker genotyping for genome-wide association studies (GWAS). Using more than 15,000 genome-wide single nucleotide polymorphisms (SNP) identified from GBS, we conducted a GWAS to investigate forage yield and nutritive value-related traits. We have detected a number of associations for all the traits evaluated and a number of associations detected were located on the Medicago truncatula genome. The SNP in a coding region of a cell wall biosynthesis gene was associated with several cell wall-related traits, and we suggest that it may be the causative polymorphism. Two other SNPs residing in meristematic development and early growth genes were found to associate with the total biomass yield. None of the SNPs associated with regrowth after harvest or with spring regrowth were mapped to the M. truncatula genome, possibly reflecting the fact that M. truncatula is an annual species related to alfalfa that typically has limited ability to regrow. The alleles we identify with the major impact on forage yield and nutritive value can be rapidly incorporated into our breeding program.

     

Characterisation of genome regions incorporated from an important wild relative into Australian sugarcane

Mandalay is an important Saccharum spontaneum clone used historically in Australian sugarcane breeding programs, and has given rise to many valuable cultivars. In order to better understand the genetic contribution of Mandalay to Australian varieties and elite parental material, a combined pedigree and quantitative trait loci (QTL) mapping approach was undertaken. A genetic map containing 400 single-dose markers was constructed for the Australian sugarcane clone MQ77340, one parent of an Australian sugarcane population (Q117 × MQ77340), using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers. This cultivar was selected because it is a direct descendent of Mandalay; its grandparents are Korpi, a S. officinarum clone, and Mandalay. The 400 markers were scattered onto 101 linkage groups (LGs) with an estimated map length of 3582 cM. The ancestral origin of all of the markers was determined with approximately 25% of the markers shown to originate from Mandalay, and a similar percentage from Korpi. Of the 101 LGs, 65 contained markers originating from Mandalay and/or Korpi. QTL analysis was undertaken using the map and 3 years of field data for three sugar-related traits (pol, brix, and CCS) and using single year field data for fibre, stalk weight and cane and sugar yield. Markers from both Mandalay and Korpi were found to be associated with both positive and negative effects on all of the traits analysed.