Search

Genotype × environment interaction (GEI) is one of the key factors affecting breeding value estimation accuracy for agronomic traits in plant breeding. Measures of GEI include fitting prediction models with various kernels to capture the variance resulting from GEI, and characterizing trials into megaenvironment (ME) clusters within which breeding values can be estimated to remove the main GEI effects. However, many of the current approaches require observations of common genotypes across all trials, which is unavailable in most breeding programs. Our study introduces two methods that can be implemented on unbalanced data to categorize trials into clusters, where both need a correlation matrix between trials: one estimated via a factor analytic (FA) model and another estimated via weather variables. The methods were tested using empirical barley (Hordeum vulgare L.) yield data in a commercial breeding program from 102 trials over 5 yr spread across multiple locations in Australia. Leave‐one‐year‐out cross‐validation achieved comparable predictive accuracies using either trials or clusters as the observed variable in GEI FA models (max. 0.45), which was higher than the accuracy achieved using the non‐GEI model (0.37). In the random cross‐validations, accuracies achieved within clusters (0.42–0.64) were mostly comparable with those achieved in the full population (0.62). In the within‐cluster validations, higher predictive accuracies were achieved when the training population was from the same cluster (mean 0.22) than outside of the cluster (mean 0.16). Our proposed methods of characterizing multienvironment trials into clusters provides a novel way to define training populations by reducing the variance resulting from GEI and could be implemented in any plant breeding program. Core Ideas: Genomic prediction implemented in breeding practiceMinimize variance caused by genotype × environment interactionclustering on multienvironment trials [ABSTRACT FROM AUTHOR]

Brassica rapa studies towards metabolic variation have largely been focused on the profiling of the diversity of metabolic compounds in specific crop types or regional varieties, but none aimed to identify genes with regulatory function in metabolite composition. Here we followed a genetical genomics approach to identify regulatory genes for six biosynthetic pathways of health-related phytochemicals, i.e carotenoids, tocopherols, folates, glucosinolates, flavonoids and phenylpropanoids. Leaves from six weeks-old plants of a Brassica rapa doubled haploid population, consisting of 92 genotypes, were profiled for their secondary metabolite composition, using both targeted and LC-MS-based untargeted metabolomics approaches. Furthermore, the same population was profiled for transcript variation using a microarray containing EST sequences mainly derived from three Brassica species: B. napus, B. rapa and B. oleracea. The biochemical pathway analysis was based on the network analyses of both metabolite QTLs (mQTLs) and transcript QTLs (eQTLs). Co-localization of mQTLs and eQTLs lead to the identification of candidate regulatory genes involved in the biosynthesis of carotenoids, tocopherols and glucosinolates. We subsequently focused on the well-characterized glucosinolate pathway and revealed two hotspots of co-localization of eQTLs with mQTLs in linkage groups A03 and A09. Our results indicate that such a large-scale genetical genomics approach combining transcriptomics and metabolomics data can provide new insights into the genetic regulation of metabolite composition of Brassica vegetables. [ABSTRACT FROM AUTHOR]

Genomic selection can increase the rate of genetic gain in crops through accumulation of positive alleles and reduce phenotyping costs by shortening the breeding cycle time. We performed genomic prediction for resistance to wheat rusts in tetraploid wheat accessions using three cross-validation with the objective of predicting: (1) rust resistance when individuals are not tested in all environments/locations, (2) the performance of lines across years, and (3) adult plant resistance (APR) of lines with bivariate models. The rationale for the latter is that seedling assays are faster and could increase prediction accuracy for APR. Predictions were derived from adult plant and seedling responses for leaf rust (Lr), stem rust (Sr) and stripe rust (Yr) in a panel of 391 accessions grown across multiple years and locations and genotyped using 16,483 single nucleotide polymorphisms. Different Bayesian models and genomic best linear unbiased prediction yielded similar accuracies for all traits. Site and year prediction accuracies for Lr and Yr ranged between 0.56–0.71 for Lr and 0.51–0.56 for Yr. While prediction accuracy for Sr was variable across different sites, accuracies for Yr were similar across different years and sites. The changes in accuracies can reflect higher genotype × environment (G × E) interactions due to climate or pathogenic variation. The use of seedling assays in genomic prediction was underscored by significant positive genetic correlations between all stage resistance (ASR) and APR (Lr: 0.45, Sr: 0.65, Yr: 0.50). Incorporating seedling phenotypes in the bivariate genomic approach increased prediction accuracy for all three rust diseases. Our work suggests that the underlying plant-host response to pathogens in the field and greenhouse screens is genetically correlated, but likely highly polygenic and therefore difficult to detect at the individual gene level. Overall, genomic prediction accuracies were in the range suitable for selection in early generations of the breeding cycle. [ABSTRACT FROM AUTHOR]

Control of gene expression is fundamental at every level of cell function. Promoter-proximal pausing and divergent transcription at promoters and enhancers, which are prominent features in animals, have only been studied in a handful of research experiments in plants. PRO-Seq analysis in cassava (Manihot esculenta) identified peaks of transcriptionally engaged RNA polymerase at both the 50 and 30 end of genes, consistent with paused or slowly moving Polymerase. In addition, we identified divergent transcription at intergenic sites. A full genome search for bi-directional transcription using an algorithm for enhancer detection developed in mammals (dREG) identified many intergenic regulatory element (IRE) candidates. These sites showed distinct patterns of methylation and nucleotide conservation based on genomic evolutionary rate profiling (GERP). SNPs within these IRE candidates explained significantly more variation in fitness and root composition than SNPs in chromosomal segments randomly ascertained from the same intergenic distribution, strongly suggesting a functional importance of these sites. Maize GRO-Seq data showed RNA polymerase occupancy at IREs consistent with patterns in cassava. Furthermore, these IREs in maize significantly overlapped with sites previously identified on the basis of open chromatin, histone marks, and methylation, and were enriched for reported eQTL. Our results suggest that bidirectional transcription can identify intergenic genomic regions in plants that play an important role in transcription regulation and whose identification has the potential to aid crop improvement. [ABSTRACT FROM AUTHOR]

Cassava (Manihot esculenta Crantz) is an important security crop that faces severe yield loses due to cassava brown streak disease (CBSD). Motivated by the slow progress of conventional breeding, genetic improvement of cassava is undergoing rapid change due to the implementation of quantitative trait loci mapping, Genome-wide association mapping (GWAS), and genomic selection (GS). In this study, two breeding panels were genotyped for SNP markers using genotyping by sequencing and phenotyped for foliar and CBSD root symptoms at five locations in Uganda. Our GWAS study found two regions associated to CBSD, one on chromosome 4 which co-localizes with a Manihot glaziovii introgression segment and one on chromosome 11, which contains a cluster of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes. We evaluated the potential of GS to improve CBSD resistance by assessing the accuracy of seven prediction models. Predictive accuracy values varied between CBSD foliar severity traits at 3 months after planting (MAP) (0.27-0.32), 6 MAP (0.40-0.42) and root severity (0.31-0.42). For all traits, Random Forest and reproducing kernel Hilbert spaces regression showed the highest predictive accuracies. Our results provide an insight into the genetics of CBSD resistance to guide CBSD marker-assisted breeding and highlight the potential of GS to improve cassava breeding.

Background: Association mapping is a statistical approach combining phenotypic traits and genetic diversity in natural populations with the goal of correlating the variation present at phenotypic and allelic levels. It is essential to separate the true effect of genetic variation from other confounding factors, such as adaptation to different uses and geographical locations. The rapid availability of large datasets makes it necessary to explore statistical methods that can be computationally less intensive and more flexible for data exploration., Methodology/principal Findings: A core collection of 168 Brassica rapa accessions of different morphotypes and origins was explored to find genetic association between markers and metabolites: tocopherols, carotenoids, chlorophylls and folate. A widely used linear model with modifications to account for population structure and kinship was followed for association mapping. In addition, a machine learning algorithm called Random Forest (RF) was used as a comparison. Comparison of results across methods resulted in the selection of a set of significant markers as promising candidates for further work. This set of markers associated to the metabolites can potentially be applied for the selection of genotypes with elevated levels of these metabolites., Conclusions/significance: The incorporation of the kinship correction into the association model did not reduce the number of significantly associated markers. However incorporation of the STRUCTURE correction (Q matrix) in the linear regression model greatly reduced the number of significantly associated markers. Additionally, our results demonstrate that RF is an interesting complementary method with added value in association studies in plants, which is illustrated by the overlap in markers identified using RF and a linear mixed model with correction for kinship and population structure. Several markers that were selected in RF and in the models with correction for kinship, but not for population structure, were also identified as QTLs in two bi-parental DH populations.

Glucosinolates and their breakdown products have been recognized for their effects on plant defense, human health, flavor and taste of cruciferous vegetables. Despite this importance, little is known about the regulation of the biosynthesis and degradation in Brassica rapa. Here, the identification of quantitative trait loci (QTL) for glucosinolate accumulation in B. rapa leaves in two novel segregating double haploid (DH) populations is reported: DH38, derived from a cross between yellow sarson R500 and pak choi variety HK Naibaicai; and DH30, from a cross between yellow sarson R500 and Kairyou Hakata, a Japanese vegetable turnip variety. An integrated map of 1068 cM with 10 linkage groups, assigned to the international agreed nomenclature, is developed based on the two individual DH maps with the common parent using amplified fragment length polymorphism (AFLP) and single sequence repeat (SSR) markers. Eight different glucosinolate compounds were detected in parents and F(1)s of the DH populations and found to segregate quantitatively in the DH populations. QTL analysis identified 16 loci controlling aliphatic glucosinolate accumulation, three loci controlling total indolic glucosinolate concentration and three loci regulating aromatic glucosinolate concentrations. Both comparative genomic analyses based on Arabidopsis-Brassica rapa synteny and mapping of candidate orthologous genes in B. rapa allowed the selection of genes involved in the glucosinolate biosynthesis pathway that may account for the identified QTL.