Genome-wide Association In Oilseed Brassicas Research Articles

Prospects and challenges for genome-wide association and genomic selection in oilseed Brassica speciesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

were submitted forpublication following the conference ‘‘Exploiting Genome-wide Association in Oilseed Brassicas: a model for geneticimprovement of major OECD crops for sustainable futurefarming’’, held at The University of Western Australia, 9–12 November 2009, which was sponsored by the OECDCo-operative Research Programme on Biological ResourceManagement for Sustainable Agricultural Systems

Generating and maintaining diversity at the elite level in crop breedingThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

Most breeding programs develop elite genotypes that are well adapted to the normal range of environmental conditions in the target production region. These elite lines have similar essential alleles for desirable end use characteristics, agronomics, disease resistance, and adaptation in the target region. The genetic makeup of these elite lines is complex. Intermating among the elite lines will often produce new variability through recombination with minimal risk of introducing new undesirable features, and is the source of most new cultivars. Eventually, this variation will be exhausted and new alleles must be introduced into the elite breeding population. Introducing desirable alleles from exotic germplasm may "pollute" the elite gene pool with undesirable alleles. Exotic germplasm may also disrupt essential allele combinations for adaptation, quality, and agronomic performance. New desirable alleles from exotic germplasm can be introgressed into an elite population in a systematic way through limited backcrossing with a minimal disturbance to the finely tuned elite background. Combining recurrent selection within elite germplasm with a systematic introgression from exotic germplasm in the recurrent introgressive population enrichment (RIPE) system has created an open-ended, continually improving, and sustainable elite population breeding system, which is simple, effective, and a regular source of new cultivars.

Future tools for association mapping in crop plantsThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Association mapping currently relies on the identification of genetic markers. Several technologies have been adopted for genetic marker analysis, with single nucleotide polymorphisms (SNPs) being the most popular where a reasonable quantity of genome sequence data are available. We describe several tools we have developed for the discovery, annotation, and visualization of molecular markers for association mapping. These include autoSNPdb for SNP discovery from assembled sequence data; TAGdb for the identification of gene specific paired read Illumina GAII data; CMap3D for the comparison of mapped genetic and physical markers; and BAC and Gene Annotator for the online annotation of genes and genomic sequences.

Whole-genome association mapping in elite inbred crop varietiesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

We have previously shown that linkage disequilibrium (LD) in the elite cultivated barley ( Hordeum vulgare ) gene pool extends, on average, for <1–5 cM. Based on this information, we have developed a platform for whole genome association studies that comprises a collection of elite lines that we have characterized at 3060 genome-wide single nucleotide polymorphism (SNP) marker loci. Interrogating this data set shows that significant population substructure is present within the elite gene pool and that diversity and LD vary considerably across each of the seven barley chromosomes. However, we also show that a subpopulation comprised of only the two-rowed spring germplasm is less structured and well suited to whole genome association studies without the need for extensive statistical intervention to account for structure. At the current marker density, the two-rowed spring population is suited for fine mapping simple traits that are located outside of the genetic centromeres with a resolution that is sufficient for candidate gene identification by exploiting conservation of synteny with fully sequenced model genomes and the emerging barley physical map.

Association mapping for phenological, morphological, and quality traits in canola quality winter rapeseed (Brassica napusL.)This article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

QTL mapping by association analysis has recently gained interest in plant breeding research as an alternative to QTL mapping in segregating populations from biparental crosses. In a first experiment on whole-genome association analysis in rapeseed, 684 mapped AFLP markers were tested for association with 14 traits in a set of 84 canola quality winter rapeseed cultivars. For association analysis a general linear model was used. By testing significance of marker-trait associations against a false discovery rate of 0.2, between 1 and 34 associated markers were found for 10 of the 14 traits. Taking into account linkage disequilibrium between the significant markers, these markers represent between 1 and 22 putative QTL for the respective traits. The minimum phenotypic variance explained by the QTL for the different traits ranged from 15% to 53%. A subset of 27 markers were significantly associated with two or more traits. These markers were predominantly shared between traits that were significantly correlated at the phenotypic level. The results show clearly that in rapeseed, QTL mapping by association analysis is a viable alternative to QTL mapping in segregating populations.

Genome-wide association analyses of common wheat (Triticum aestivum L.) germplasm identifies multiple loci for aluminium resistanceThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”

Aluminium (Al3+) toxicity restricts productivity and profitability of wheat (Triticum aestivum L.) crops grown on acid soils worldwide. Continued gains will be obtained by identifying superior alleles and novel Al3+ resistance loci that can be incorporated into breeding programs. We used association mapping to identify genomic regions associated with Al3+ resistance using 1055 accessions of common wheat from different geographic regions of the world and 178 polymorphic diversity arrays technology (DArT) markers. Bayesian analyses based on genetic distance matrices classified these accessions into 12 subgroups. Genome-wide association analyses detected markers that were significantly associated with Al3+ resistance on chromosomes 1A, 1B, 2A, 2B, 2D, 3A, 3B, 4A, 4B, 4D, 5B, 6A, 6B, 7A, and 7B. Some of these genomic regions correspond to previously identified loci for Al3+ resistance, whereas others appear to be novel. Among the markers targeting TaALMT1 (the major Al3+-resistance gene located on chromosome 4D), those that detected alleles in the promoter explained most of the phenotypic variance for Al3+ resistance, which is consistent with this region controlling the level of TaALMT1 expression. These results demonstrate that genome-wide association mapping cannot only confirm known Al3+-resistance loci, such as those on chromosomes 4D and 4B, but they also highlight the utility of this technique in identifying novel resistance loci.

Genome-wide association and genomic selection in animal breedingThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Results from genome-wide association studies in livestock, and humans, has lead to the conclusion that the effect of individual quantitative trait loci (QTL) on complex traits, such as yield, are likely to be small; therefore, a large number of QTL are necessary to explain genetic variation in these traits. Given this genetic architecture, gains from marker-assisted selection (MAS) programs using only a small number of DNA markers to trace a limited number of QTL is likely to be small. This has lead to the development of alternative technology for using the available dense single nucleotide polymorphism (SNP) information, called genomic selection. Genomic selection uses a genome-wide panel of dense markers so that all QTL are likely to be in linkage disequilibrium with at least one SNP. The genomic breeding values are predicted to be the sum of the effect of these SNPs across the entire genome. In dairy cattle breeding, the accuracy of genomic estimated breeding values (GEBV) that can be achieved and the fact that these are available early in life have lead to rapid adoption of the technology. Here, we discuss the design of experiments necessary to achieve accurate prediction of GEBV in future generations in terms of the number of markers necessary and the size of the reference population where marker effects are estimated. We also present a simple method for implementing genomic selection using a genomic relationship matrix. Future challenges discussed include using whole genome sequence data to improve the accuracy of genomic selection and management of inbreeding through genomic relationships.

Association mapping of seed oil content in Brassica napus and comparison with quantitative trait loci identified from linkage mappingThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Association mapping has been used increasingly in natural populations with rich genetic diversity to detect DNA-based markers that are associated with important agronomic traits. Brassica napus is an important oil crop with limited genetic diversity. "New-type" B. napus that is introgressed with subgenomic components from related species has been developed to broaden the genetic basis of "traditional" B. napus. In this study, new-type B. napus lines and a collection of traditional B. napus varieties from different countries were used as two different populations to evaluate seed oil content and to determine the efficacy of association mapping by comparison with previous study of linkage mapping. Relatively rich genetic diversity, but a higher level of linkage disequilibrium was observed in the new-type B. napus as compared with the traditional B. napus. Similarly, a larger variation in oil content and a greater number of associated markers were detected in the population of new-type B. napus. Meanwhile, more than half of the genetic loci, to which the associated markers corresponded, were located within the quantitative trait loci intervals identified previously in linkage mapping experiments, which demonstrated the power of association mapping in B. napus.

Statistical genetic issues for genome-wide association studiesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Genotyping technology now allows the rapid and affordable generation of million-SNP profiles for humans, leading to considerable activity in association mapping. Similar activity is anticipated for many plant species, including Brassica. These plant association mapping activities will require the same care in quality control and quality assurance as for humans. The subsequent analyses may draw upon the same body of theory that is described here in the language of quantitative genetics.

Analysis of yield and oil from a series of canola breeding trials. Part II. Exploring variety by environment interaction using factor analysisThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Exploring and exploiting variety by environment (V × E) interaction is one of the major challenges facing plant breeders. In paper I of this series, we presented an approach to modelling V × E interaction in the analysis of complex multi-environment trials using factor analytic models. In this paper, we develop a range of statistical tools which explore V × E interaction in this context. These tools include graphical displays such as heat-maps of genetic correlation matrices as well as so-called E-scaled uniplots that are a more informative alternative to the classical biplot for large plant breeding multi-environment trials. We also present a new approach to prediction for multi-environment trials that include pedigree information. This approach allows meaningful selection indices to be formed either for potential new varieties or potential parents.

Exploring and exploiting epigenetic variation in cropsThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

This review addresses the mechanisms by which epigenetic variation modulates plant gene regulation and phenotype. In particular we explore the scope for harnessing such processes within the context of crop genetic improvement. We focus on the role of DNA methylation as an epigenetic mark that contributes to epiallelic diversity and modulation of gene regulation. We outline the prevalence and distribution of epigenetic marks in relation to eukaryote developmental processes, and in particular identify where this may be relevant to crop traits both in terms of specific developmental stages and in relation to physiological responses to environmental change. Recent whole genome surveys have identified specific characteristics of the distribution of DNA methylation within plant genomes. Together with greater understanding of the mode of action of different maintenance and de novo methyltransferases, this provides an opportunity to modulate DNA methylation status at specific loci as an intervention strategy in crop genetic improvement. We discuss alternative approaches that may be suitable for harnessing such induced epiallelic variation. Most of the discussion is associated with Brassica crops, which demonstrate considerable morphological plasticity, segmental chromosomal duplication, and polyploidy.

Simultaneously accounting for population structure, genotype by environment interaction, and spatial variation in marker–trait associations in sugarcaneThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Few association mapping studies have simultaneously accounted for population structure, genotype by environment interaction (GEI), and spatial variation. In this sugarcane association mapping study we tested models accounting for these factors and identified the impact that each model component had on the list of markers declared as being significantly associated with traits. About 480 genotypes were evaluated for cane yield and sugar content at three sites and scored with DArT markers. A mixed model was applied in analysis of the data to simultaneously account for the impacts of population structure, GEI, and spatial variation within a trial. Two forms of the DArT marker data were used in the analysis: the standard discrete data (0, 1) and a continuous DArT score, which is related to the marker dosage. A large number of markers were significantly associated with cane yield and sugar content. However, failure to account for population structure, GEI, and (or) spatial variation produced both type I and type II errors, which on the one hand substantially inflated the number of significant markers identified (especially true for failing to account for GEI) and on the other hand resulted in failure to detect markers that could be associated with cane yield or sugar content (especially when failing to account for population structure). We concluded that association mapping based on trials from one site or analysis that failed to account for GEI would produce many trial-specific associated markers that would have low value in breeding programs.

Development of a high density integrated reference genetic linkage map for the multinational Brassica rapa Genome Sequencing ProjectThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

We constructed a high-density Brassica rapa integrated linkage map by combining a reference genetic map of 78 doubled haploid lines derived from Chiifu-401-42× Kenshin (CKDH) and a new map of 190 F2 lines derived from Chiifu-401-42× rapid cycling B. rapa (CRF2). The integrated map contains 1017 markers and covers 1262.0 cM of the B. rapa genome, with an average interlocus distance of 1.24 cM. High similarity of marker order and position was observed among the linkage groups of the maps with few short-distance inversions. In total, 155 simple sequence repeat (SSR) markers, anchored to 102 new bacterial artificial chromosomes (BACs) and 146 intron polymorphic (IP) markers were mapped in the integrated map, which would be helpful to align the sequenced BACs in the ongoing multinational Brassica rapa Genome Sequencing Project (BrGSP). Further, comparison of the B. rapa consensus map with the 10 B. juncea A-genome linkage groups by using 98 common IP markers showed high-degree colinearity between the A-genome linkage groups, except for few markers showing inversion or translocation. Suggesting that chromosomes are highly conserved between these Brassica species, although they evolved independently after divergence. The sequence information coming out of BrGSP would be useful for B. juncea breeding. and the identified Arabidopsis chromosomal blocks and known quantitative trait loci (QTL) information of B. juncea could be applied to improve other Brassica crops including B. rapa.

Analysis of yield and oil from a series of canola breeding trials. Part I. Fitting factor analytic mixed models with pedigree informationThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

In this paper multiplicative mixed models have been used for the analysis of multi-environment trial (MET) data for canola oil and grain yield. Information on pedigrees has been included to allow for the modelling of additive and nonadditive genetic effects. The MET data set included a total of 19 trials (synonymous with sites or environments), which were sown across southern Australia in 2007 and 2008. Each trial was designed as a p-rep design using DiGGeR with the default prespecified spatial model. Lines in their first year of testing were unreplicated, whereas there were two or three replications of advanced lines or varieties. Pedigree information on a total of 578 entries was available, and there were 69 entries that had unknown pedigrees. The degree of inbreeding varied from 0 (55 entries) to nearly fully inbred (337 entries). Subsamples of 2 g harvested grain were taken from each plot for determination of seed oil percentage by near infrared reflectance spectroscopy. The MET analysis for both yield and oil modelled genetic effects in different trials using factor analytic models and the residual plot effects for each trial were modelled using spatial techniques. Models in which pedigree information was included provided significantly better fits to both yield and oil data.

Towards unambiguous transcript mapping in the allotetraploid Brassica napusThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

The architecture of the Brassica napus genome is marked by its evolutionary origins. The genome of B. napus was formed from the hybridization of two closely related diploid Brassica species, both of which evolved from an hexaploid ancestor. The extensive whole genome duplication events in its near and distant past result in the allotetraploid genome of B. napus maintaining multiple copies of most genes, which predicts a highly complex and redundant transcriptome that can confound any expression analyses. A stringent assembly of 142,399 B. napus expressed sequence tags allowed the development of a well-differentiated set of reference transcripts, which were used as a foundation to assess the efficacy of available tools for identifying and distinguishing transcripts in B. napus; including microarray hybridization and 3' anchored sequence tag capture. Microarray platforms cannot distinguish transcripts derived from the two progenitors or close homologues, although observed differential expression appeared to be biased towards unique transcripts. The use of 3' capture enhanced the ability to unambiguously identify homologues within the B. napus transcriptome but was limited by tag length. The ability to comprehensively catalogue gene expression in polyploid species could be transformed by the application of cost-efficient next generation sequencing technologies that will capture millions of long sequence tags.

Design of aBrassica rapacore collection for association mapping studiesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

A Brassica rapa collection of 239 accessions, based on two core collections representing different morphotypes from different geographical origins, is presented and its use for association mapping is illustrated for flowering time. We analyzed phenotypic variation of leaf and seed pod traits, plant architecture, and flowering time using data collected from three field experiments and evaluated the genetic diversity with a set of SSR markers. The Wageningen University and Research Centre (WUR) and the Vavilov Research Institute of Plant Industry (VIR) core collections had similar representations of most morphotypes, as illustrated by the phenotypic and genetic variation within these groups. The analysis of population structure revealed five subgroups in the collection, whereas previous studies of the WUR core collection indicated four subgroups; the fifth group identified consisted mainly of oil accessions from the VIR core collection, winter oils from Pakistan, and a number of other types. A very small group of summer oils is described, that is not related to other oil accessions. A candidate gene approach was chosen for association mapping of flowering time with a BrFLC1 biallelic CAPS marker and a BrFLC2 multiallelic SSR marker. The two markers were significantly associated with flowering time, but their effects were confined to certain morphotypes and (or) alleles. Based on these results, we discuss the optimal design for an association mapping population and the need to fix the heterogeneous accessions to facilitate phenotyping and genotyping.

SNP discovery by amplicon sequencing and multiplex SNP genotyping in the allopolyploid species Brassica napusThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Oilseed rape (Brassica napus) is an allotetraploid species consisting of two genomes, derived from B. rapa (A genome) and B. oleracea (C genome). The presence of these two genomes makes single nucleotide polymorphism (SNP) marker identification and SNP analysis more challenging than in diploid species, as for a given locus usually two versions of a DNA sequence (based on the two ancestral genomes) have to be analyzed simultaneously during SNP identification and analysis. One hundred amplicons derived from expressed sequence tag (ESTs) were analyzed to identify SNPs in a panel of oilseed rape varieties and within two sister species representing the ancestral genomes. A total of 604 SNPs were identified, averaging one SNP in every 42 bp. It was possible to clearly discriminate SNPs that are polymorphic between different plant varieties from SNPs differentiating the two ancestral genomes. To validate the identified SNPs for their use in genetic analysis, we have developed Illumina GoldenGate assays for some of the identified SNPs. Through the analysis of a number of oilseed rape varieties and mapping populations with GoldenGate assays, we were able to identify a number of different segregation patterns in allotetraploid oilseed rape. The majority of the identified SNP markers can be readily used for genetic mapping, showing that amplicon sequencing and Illumina GoldenGate assays can be used to reliably identify SNP markers in tetraploid oilseed rape and to convert them into successful SNP assays that can be used for genetic analysis.

Regional association analysis delineates a sequenced chromosome region influencing antinutritive seed meal compounds in oilseed rapeThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

This study describes the use of regional association analyses to delineate a sequenced region of a Brassica napus chromosome with a significant effect on antinutritive seed meal compounds in oilseed rape. A major quantitative trait locus (QTL) influencing seed colour, fibre content, and phenolic compounds was mapped to the same position on B. napus chromosome A9 in biparental mapping populations from two different yellow-seeded × black-seeded B. napus crosses. Sequences of markers spanning the QTL region identified synteny to a sequence contig from the corresponding chromosome A9 in Brassica rapa. Remapping of sequence-derived markers originating from the B. rapa sequence contig confirmed their position within the QTL. One of these markers also mapped to a seed colour and fibre QTL on the same chromosome in a black-seeded × black-seeded B. napus cross. Consequently, regional association analysis was performed in a genetically diverse panel of dark-seeded, winter-type oilseed rape accessions. For this we used closely spaced simple sequence repeat (SSR) markers spanning the sequence contig covering the QTL region. Correction for population structure was performed using a set of genome-wide SSR markers. The identification of QTL-derived markers with significant associations to seed colour, fibre content, and phenolic compounds in the association panel enabled the identification of positional and functional candidate genes for B. napus seed meal quality within a small segment of the B. rapa genome sequence.

Exploiting genome-wide association in oilseed Brassica speciesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.

Exploiting genome-wide association in oilseed <i>Brassica</i> speciesThis article is one of a selection of papers from the conference “Exploiting Genome-wide Association in Oilseed Brassicas: a model for genetic improvement of major OECD crops for sustainable farming”.