QTL meta-analysis of root traits in Brassica napus under contrasting phosphorus supply in two growth systems

Quantitative Proteomic Analysis of Brassica Napus Reveals Intersections Between Nutrient Deficiency Responses

Macronutrients such as nitrogen (N), phosphorus (P), potassium (K) and sulphur (S) are critical for plant growth and development. Field-grown canola (Brassica napus L.) is supplemented with fertilizers to maximize plant productivity, while deficiency in these nutrients can cause significant yield loss. A holistic understanding of the interplay between these nutrient deficiency responses in a single study and canola cultivar is thus far lacking, hindering efforts to increase the nutrient use efficiency of this important oil seed crop. To address this, we performed a comparative quantitative proteomic analysis of both shoot and root tissue harvested from soil-grown canola plants experiencing either nitrogen, phosphorus, potassium or sulphur deficiency. Our data provide critically needed insights into the shared and distinct molecular responses to macronutrient deficiencies in canola. Importantly, we find more conserved responses to the four different nutrient deficiencies in canola roots, with more distinct proteome changes in aboveground tissue. Our results establish a foundation for a more comprehensive understanding of the shared and distinct nutrient deficiency response mechanisms of canola plants and pave the way for future breeding efforts.

Identification of quantitative trait loci and candidate genes for pod shatter resistance in Brassica carinata

BACKGROUND

Understanding the genetic control of pod shatter resistance and its association with pod length is crucial for breeding improved pod shatter resistance and reducing pre-harvest yield losses due to extensive shattering in cultivars of Brassica species. In this study, we evaluated a doubled haploid (DH) mapping population derived from an F1 cross between two Brassica carinata parental lines Y-BcDH64 and W-BcDH76 (YWDH), originating from Ethiopia and determined genetic bases of variation in pod length and pod shatter resistance, measured as rupture energy. The YWDH population, its parental lines and 11 controls were grown across three years for genetic analysis.

RESULTS

By using three quantitative trait loci (QTL) analytic approaches, we identified nine genomic regions on B02, B03, B04, B06, B07 and C01 chromosomes for rupture energy that were repeatedly detected across three growing environments. One of the QTL on chromosome B07, flanked with DArTseq markers 100,046,735 and 100,022,658, accounted for up to 27.6% of genetic variance in rupture energy. We observed no relationship between pod length and rupture energy, suggesting that pod length does not contribute to variation in pod shatter resistance. Comparative mapping identified six candidate genes; SHP1 on B6, FUL and MAN on chromosomes B07, IND and NST2 on B08, and MAN7 on C07 that mapped within 0.2 Mb from the QTL for rupture energy.

CONCLUSION

The results suggest that favourable alleles of stable QTL on B06, B07, B08 and C01 for pod shatter resistance can be incorporated into the shatter-prone B. carinata and its related species to improve final seed yield at harvest.

Genome-wide association studies of root system architecture traits in a broad collection of Brassica genotypes

The root systems of Brassica species are complex. Eight root system architecture (RSA) traits, including total root length, total root surface area, root average diameter, number of tips, total primary root length, total lateral root length, total tertiary root length, and basal link length, were phenotyped across 379 accessions representing six Brassica species (B. napus, B. juncea, B. carinata, B. oleracea, B. nigra, and B. rapa) using a semi-hydroponic system and image analysis software. The results suggest that, among the assessed species, B. napus and B. oleracea had the most intricate and largest root systems, while B. nigra exhibited the smallest roots. The two species B. juncea and B. carinata shared comparable root system complexity and had root systems with larger root diameters. In addition, 313 of the Brassica accessions were genotyped using a 19K Brassica single nucleotide polymorphism (SNP) array. After filtering by TASSEL 5.0, 6,213 SNP markers, comprising 5,103 markers on the A-genome (covering 302,504 kb) and 1,110 markers on the C-genome (covering 452,764 kb), were selected for genome-wide association studies (GWAS). Two general linear models were tested to identify the genomic regions and SNPs associated with the RSA traits. GWAS identified 79 significant SNP markers associated with the eight RSA traits investigated. These markers were distributed across the 18 chromosomes of B. napus, except for chromosome C06. Sixty-five markers were located on the A-genome, and 14 on the C-genome. Furthermore, the major marker-trait associations (MTAs)/quantitative trait loci (QTLs) associated with root traits were located on chromosomes A02, A03, and A06. Brassica accessions with distinct RSA traits were identified, which could hold functional, adaptive, evolutionary, environmental, pathological, and breeding significance.

Deep learning-based association analysis of root image data and cucumber yield

The root system is important for the absorption of water and nutrients by plants. Cultivating and selecting a root system architecture (RSA) with good adaptability and ultrahigh productivity have become the primary goals of agricultural improvement. Exploring the correlation between the RSA and crop yield is important for cultivating crop varieties with high-stress resistance and productivity. In this study, 277 cucumber varieties were collected for root system image analysis and yield using germination plates and greenhouse cultivation. Deep learning tools were used to train ResNet50 and U-Net models for image classification and segmentation of seedlings and to perform quality inspection and productivity prediction of cucumber seedling root system images. The results showed that U-Net can automatically extract cucumber root systems with high quality (F1_score ≥ 0.95), and the trained ResNet50 can predict cucumber yield grade through seedling root system image, with the highest F1_score reaching 0.86 using 10-day-old seedlings. The root angle had the strongest correlation with yield, and the shallow- and steep-angle frequencies had significant positive and negative correlations with yield, respectively. RSA and nutrient absorption jointly affected the production capacity of cucumber plants. The germination plate planting method and automated root system segmentation model used in this study are convenient for high-throughput phenotypic (HTP) research on root systems. Moreover, using seedling root system images to predict yield grade provides a new method for rapidly breeding high-yield RSA in crops such as cucumbers.

Genetic control of root morphology in response to nitrogen across rapeseed diversity

Rapeseed (Brassica napus L.) is an oil-containing crop of great economic value but with considerable nitrogen requirement. Breeding root systems that efficiently absorb nitrogen from the soil could be a driver to ensure genetic gains for more sustainable rapeseed production. The aim of this study is to identify genomic regions that regulate root morphology in response to nitrate availability. The natural variability offered by 300 inbred lines was screened at two experimental locations. Seedlings grew hydroponically with low or elevated nitrate levels. Fifteen traits related to biomass production and root morphology were measured. On average across the panel, a low nitrate level increased the root-to-shoot biomass ratio and the lateral root length. A large phenotypic variation was observed, along with important heritability values and genotypic effects, but low genotype-by-nitrogen interactions. Genome-wide association study and bulk segregant analysis were used to identify loci regulating phenotypic traits. The first approach nominated 319 SNPs that were combined into 80 QTLs. Three QTLs identified on the A07 and C07 chromosomes were stable across nitrate levels and/or experimental locations. The second approach involved genotyping two groups of individuals from an experimental F2 population created by crossing two accessions with contrasting lateral root lengths. These individuals were found in the tails of the phenotypic distribution. Co-localized QTLs found in both mapping approaches covered a chromosomal region on the A06 chromosome. The QTL regions contained some genes putatively involved in root organogenesis and represent selection targets for redesigning the root morphology of rapeseed.

Genome-wide association study reveals candidate genes controlling root system architecture under low phosphorus supply at seedling stage in Brassica napus

Optimal root system architecture (RSA) is essential for vigorous growth and yield in crops. Plants have evolved adaptive mechanisms in response to low phosphorus (LP) stress, and one of those is changes in RSA. Here, more than five million single-nucleotide polymorphisms (SNPs) obtained from whole-genome re-sequencing data (WGR) of an association panel of 370 oilseed rape (Brassica napus L.) were used to conduct a genome-wide association study (GWAS) of RSA traits of the panel at LP in "pouch and wick" system. Fifty-two SNPs were forcefully associated with lateral root length (LRL), total root length (TRL), lateral root density (LRD), lateral root number (LRN), mean lateral root length (MLRL), and root dry weight (RDW) at LP. There were significant correlations between phenotypic variation and the number of favorable alleles of the associated loci on chromosomes A06 (chrA06_20030601), C03 (chrC03_3535483), and C07 (chrC07_42348561), respectively. Three candidate genes (BnaA06g29270D, BnaC03g07130D, and BnaC07g43230D) were detected by combining transcriptome, candidate gene association analysis, and haplotype analysis. Cultivar carrying "CCGC" at BnaA06g29270DHap1, "CAAT" at BnaC03g07130DHap1, and "ATC" at BnaC07g43230DHap1 had greater LRL, LRN, and RDW than lines carrying other haplotypes at LP supply. The RSA of a cultivar harboring the three favorable haplotypes was further confirmed by solution culture experiments. These findings define exquisite insights into genetic architectures underlying B. napus RSA at LP and provide valuable gene resources for root breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01411-2.

Integrating genome-wide association studies with selective sweep reveals genetic loci associated with tolerance to low phosphate availability in Brassica napus

Oilseed rape (Brassica napus L.; B. napus) is an important oil crop worldwide. However, the genetic mechanisms of B. napus adaptations to low phosphate (P) stress are largely unknown. In this study, a genome-wide association study (GWAS) identified 68 SNPs significantly associated with seed yield (SY) under low P (LP) availability, and 7 SNPs significantly associated with phosphorus efficiency coefficient (PEC) in two trials. Among these SNPs, two, chrC07__39807169 and chrC09__14194798, were co-detected in two trials, and BnaC07.ARF9 and BnaC09.PHT1;2 were identified as candidate genes of them, respectively, by combining GWAS with quantitative reverse-transcription PCR (qRT-PCR). There were significant differences in the gene expression level of BnaC07.ARF9 and BnaC09.PHT1;2 between P-efficient and -inefficiency varieties at LP. SY_LP had a significant positive correlation with the gene expression level of both BnaC07.ARF9 and BnaC09.PHT1;2. BnaC07.ARF9 and BnaA01.PHR1 could directly bind the promoters of BnaA01.PHR1 and BnaC09.PHT1;2, respectively. Selective sweep analysis was conducted between ancient and derived B. napus, and detected 1280 putative selective signals. Within the selected region, a large number of genes related to P uptake, transport, and utilization were detected, such as purple acid phosphatase (PAP) family genes and phosphate transporter (PHT) family genes. These findings provide novel insights into the molecular targets for breeding P efficiency varieties in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01399-9.

Rapid identification of a major locus qPRL-C06 affecting primary root length in Brassica napus by QTL-seq

BACKGROUND AND AIMS

Brassica napus is one of the most important oilseed crops worldwide. Seed yield of B. napus significantly correlates with the primary root length (PRL). The aims of this study were to identify quantitative trait loci (QTLs) for PRL in B. napus.

METHODS

QTL-seq and conventional QTL mapping were jointly used to detect QTLs associated with PRL in a B. napus double haploid (DH) population derived from a cross between 'Tapidor' and 'Ningyou 7'. The identified major locus was confirmed and resolved by an association panel of B. napus and an advanced backcross population. RNA-seq analysis of two long-PRL lines (Tapidor and TN20) and two short-PRL lines (Ningyou 7 and TN77) was performed to identify differentially expressed genes in the primary root underlying the target QTLs.

KEY RESULTS

A total of 20 QTLs impacting PRL in B. napus grown at a low phosphorus (P) supply were found by QTL-seq. Eight out of ten QTLs affecting PRL at a low P supply discovered by conventional QTL mapping could be detected by QTL-seq. The locus qPRL-C06 identified by QTL-seq was repeatedly detected at both an optimal P supply and a low P supply by conventional QTL mapping. This major constitutive QTL was further confirmed by regional association mapping. qPRL-C06 was delimited to a 0.77 Mb genomic region on chromosome C06 using an advanced backcross population. A total of 36 candidate genes within qPRL-C06 were identified that showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in primary root between the long-PRL and short-PRL lines, including five genes involved in phytohormone biosynthesis and signaling.

CONCLUSIONS

These results both demonstrate the power of the QTL-seq in rapid QTL detection for root traits and will contribute to marker-assisted selective breeding of B. napus cultivars with increased PRL.

Genetic variation in BnGRP1 contributes to low phosphorus tolerance in Brassica napus

Lack of phosphorus (P) is a major environmental factor affecting rapeseed (Brassica napus. L) root growth and development. For breeding purposes, it is crucial to identify the molecular mechanisms of root system architecture (RSA) traits underlying low P tolerance in rapeseed. The natural variations in the glycine-rich protein gene, BnGRP1, were analyzed in the natural population of 400 rapeseed cultivars under low P stress through genome-wide association study (GWAS) and transcriptome analyses. Based on 11 SNP mutations in BnGRP1 sequence, ten types of haplotypes (Hap) were formed. Compared with the other types, the cultivar of BnGRP1Hap1 type in the panel demonstrated the longest root length and heaviest root weight. BnGRP1Hap1 overexpression in rapeseed depicted the ability to enhance its resistance in response to low P tolerance. CRISPR/Cas9-derived BnGRP1Hap4 knockout mutations in rapeseed can lead to sensitivity to low P stress. Furthermore, BnGRP1Hap1 influences the expression of phosphate transporter 1 (PHT1) genes associated with P absorption. Overall, the findings of this study highlight new mechanisms of GRP1 genes in enhancing low P tolerance in rapeseed.

Genome-Wide Association Studies of Root-Related Traits in Brassica napus L. under Low-Potassium Conditions

Roots are essential organs for a plant’s ability to absorb water and obtain mineral nutrients, hence they are critical to its development. Plants use root architectural alterations to improve their chances of absorbing nutrients when their supply is low. Nine root traits of a Brassica napus association panel were explored in hydroponic-system studies under low potassium (K) stress to unravel the genetic basis of root growth in rapeseed. The quantitative trait loci (QTL) and candidate genes for root development were discovered using a multilocus genome-wide association study (ML-GWAS). For the nine traits, a total of 453 significant associated single-nucleotide polymorphism (SNP) loci were discovered, which were then integrated into 206 QTL clusters. There were 45 pleiotropic clusters, and qRTA04-4 and qRTC04-7 were linked to TRL, TSA, and TRV at the same time, contributing 5.25–11.48% of the phenotypic variance explained (PVE) to the root traits. Additionally, 1360 annotated genes were discovered by examining genomic regions within 100 kb upstream and downstream of lead SNPs within the 45 loci. Thirty-five genes were identified as possibly regulating root-system development. As per protein–protein interaction analyses, homologs of three genes (BnaC08g29120D, BnaA07g10150D, and BnaC04g45700D) have been shown to influence root growth in earlier investigations. The QTL clusters and candidate genes identified in this work will help us better understand the genetics of root growth traits and could be employed in marker-assisted breeding for rapeseed adaptable to various conditions with low K levels.

Discovery of Genomic Regions and Candidate Genes Controlling Root Development Using a Recombinant Inbred Line Population in Rapeseed (Brassica napus L.)

Marker-assisted selection enables breeders to quickly select excellent root architectural variations, which play an essential role in plant productivity. Here, ten root-related and shoot biomass traits of a new F₆ recombinant inbred line (RIL) population were investigated under hydroponics and resulted in high heritabilities from 0.61 to 0.83. A high-density linkage map of the RIL population was constructed using a Brassica napus 50k Illumina single nucleotide polymorphism (SNP) array. A total of 86 quantitative trait loci (QTLs) explaining 4.16–14.1% of the phenotypic variances were detected and integrated into eight stable QTL clusters, which were repeatedly detected in different experiments. The codominant markers were developed to be tightly linked with three major QTL clusters, qcA09-2, qcC08-2, and qcC08-3, which controlled both root-related and shoot biomass traits and had phenotypic contributions greater than 10%. Among these, qcA09-2, renamed RT.A09, was further fine-mapped to a 129-kb interval with 19 annotated genes in the B. napus reference genome. By integrating the results of real-time PCR and comparative sequencing, five genes with expression differences and/or amino acid differences were identified as important candidate genes for RT.A09. Our findings laid the foundation for revealing the molecular mechanism of root development and developed valuable markers for root genetic improvement in rapeseed.

Improving phosphate use efficiency in the aquatic crop watercress (Nasturtium officinale)

Watercress is a nutrient-dense leafy green crop, traditionally grown in aquatic outdoor systems and increasingly seen as well-suited for indoor hydroponic systems. However, there is concern that this crop has a detrimental impact on the environment through direct phosphate additions causing environmental pollution. Phosphate-based fertilisers are supplied to enhanced crop yield, but their use may contribute to eutrophication of waterways downstream of traditional watercress farms. One option is to develop a more phosphate use efficient (PUE) crop. This review identifies the key traits for this aquatic crop (the ideotype), for future selection, marker development and breeding. Traits identified as important for PUE are (i) increased root surface area through prolific root branching and adventitious root formation, (ii) aerenchyma formation and root hair growth. Functional genomic traits for improved PUE are (iii) efficacious phosphate remobilisation and scavenging strategies and (iv) the use of alternative metabolic pathways. Key genomic targets for this aquatic crop are identified as: PHT phosphate transporter genes, global transcriptional regulators such as those of the SPX family and genes involved in galactolipid and sulfolipid biosynthesis such as MGD2/3, PECP1, PSR2, PLDζ1/2 and SQD2. Breeding for enhanced PUE in watercress will be accelerated by improved molecular genetic resources such as a full reference genome sequence that is currently in development.

Inoculation with Plant Growth-Promoting Bacteria to Reduce Phosphate Fertilization Requirement and Enhance Technological Quality and Yield of Sugarcane

Phosphorus (P) is a critical nutrient for high sugarcane yields throughout its cultivation cycles, however, a higher amount of P becomes rapidly unavailable to plants due to its adsorption to soil colloids. Some plant growth-promoting bacteria (PGPBs) may be able to enhance P availability to plants and produce phytohormones that contribute to crop development, quality, and yield. Thus, this study aimed to evaluate leaf concentrations of nitrogen (N) and P, yield, and technological quality of sugarcane as a function of different levels of phosphate fertilization associated with inoculation of PGPBs. The experiment was carried out at Ilha Solteira, São Paulo-Brazil. The experimental design was randomized blocks with three replications, consisting of five phosphorus rates (0, 25, 50, 75, and 100% of the recommended P2O5 rate) and eight inoculations, involving three species of PGPBs (Azospirillum brasilense, Bacillus subtilis, and Pseudomonas fluorescens) which were applied combined or in a single application into the planting furrow of RB92579 sugarcane variety. The inoculation of B. subtilis and P. fluorescens provided a higher concentration of leaf P in sugarcane. The P2O5 rates combined with inoculation of bacteria alter technological variables and stalk yield of sugarcane. The excess and lack of phosphate fertilizer is harmful to sugarcane cultivation, regardless of the use of growth-promoting bacteria. We recommend the inoculation with A. brasilense + B. subtilis associated with 45 kg ha-1 of P2O5 aiming at greater stalk yield. This treatment also increases sugar yield, resulting in a savings of 75% of the recommended P2O5 rate, thus being a more efficient and sustainable alternative for reducing sugarcane crop production costs.

Exploiting High-Throughput Indoor Phenotyping to Characterize the Founders of a Structured B. napus Breeding Population

Phenotyping is considered a significant bottleneck impeding fast and efficient crop improvement. Similar to many crops, Brassica napus, an internationally important oilseed crop, suffers from low genetic diversity, and will require exploitation of diverse genetic resources to develop locally adapted, high yielding and stress resistant cultivars. A pilot study was completed to assess the feasibility of using indoor high-throughput phenotyping (HTP), semi-automated image processing, and machine learning to capture the phenotypic diversity of agronomically important traits in a diverse B. napus breeding population, SKBnNAM, introduced here for the first time. The experiment comprised 50 spring-type B. napus lines, grown and phenotyped in six replicates under two treatment conditions (control and drought) over 38 days in a LemnaTec Scanalyzer 3D facility. Growth traits including plant height, width, projected leaf area, and estimated biovolume were extracted and derived through processing of RGB and NIR images. Anthesis was automatically and accurately scored (97% accuracy) and the number of flowers per plant and day was approximated alongside relevant canopy traits (width, angle). Further, supervised machine learning was used to predict the total number of raceme branches from flower attributes with 91% accuracy (linear regression and Huber regression algorithms) and to identify mild drought stress, a complex trait which typically has to be empirically scored (0.85 area under the receiver operating characteristic curve, random forest classifier algorithm). The study demonstrates the potential of HTP, image processing and computer vision for effective characterization of agronomic trait diversity in B. napus, although limitations of the platform did create significant variation that limited the utility of the data. However, the results underscore the value of machine learning for phenotyping studies, particularly for complex traits such as drought stress resistance.

Oilseed Rape Cultivars Show Diversity of Root Morphologies with the Potential for Better Capture of Nitrogen

The worldwide demand for vegetable oils is rising. Oilseed rape (Brassica napus) diversifies cereal dominated crop rotations but requires important nitrogen input. Yet, the root organ is offering an untapped opportunity to improve the nitrogen capture in soil. This study evaluates three culture systems in controlled environment, to observe root morphology and to identify root attributes for superior biomass production and nitrogen use. The phenotypic diversity in a panel of 55 modern winter oilseed rape cultivars was screened in response to two divergent nitrate supplies. Upon in vitro and hydroponic cultures, a large variability for root morphologies was observed. Root biomass and morphological traits positively correlated with shoot biomass or leaf area. The activities of high-affinity nitrate transport systems correlated negatively with the leaf area, while the combined high- and low-affinity systems positively with the total root length. The X-ray computed tomography permitted to visualize the root system in pipes filled with soil. The in vitro root phenotype at germination stage was indicative of lateral root deployment in soil-grown plants. This study highlights great genetic potential in oilseed rape, which could be manipulated to optimize crop root characteristics and nitrogen capture with substantial implications for agricultural production.

Genetic Dissection of Mature Root Characteristics by Genome-Wide Association Studies in Rapeseed (Brassica napus L.)

Roots are complicated quantitative characteristics that play an essential role in absorbing water and nutrients. To uncover the genetic variations for root-related traits in rapeseed, twelve mature root traits of a Brassica napus association panel were investigated in the field within three environments. All traits showed significant phenotypic variation among genotypes, with heritabilities ranging from 55.18% to 79.68%. Genome-wide association studies (GWAS) using 20,131 SNPs discovered 172 marker-trait associations, including 103 significant SNPs (-log10 (p) > 4.30) that explained 5.24-20.31% of the phenotypic variance. With the linkage disequilibrium r² > 0.2, these significant associations were binned into 40 quantitative trait loci (QTL) clusters. Among them, 14 important QTL clusters were discovered in two environments and/or with phenotypic contributions greater than 10%. By analyzing the genomic regions within 100 kb upstream and downstream of the peak SNPs within the 14 loci, 334 annotated genes were found. Among these, 32 genes were potentially associated with root development according to their expression analysis. Furthermore, the protein interaction network using the 334 annotated genes gave nine genes involved in a substantial number of interactions, including a key gene associated with root development, BnaC09g36350D. This research provides the groundwork for deciphering B. napus' genetic variations and improving its root system architecture.

Genomic analysis of ionome-related QTLs in Arabidopsis thaliana

Ionome contributes to maintain cell integrity and acts as cofactors for catalyzing regulatory pathways. Identifying ionome contributing genomic regions provides a practical framework to dissect the genetic architecture of ionomic traits for use in biofortification. Meta-QTL (MQTL) analysis is a robust method to discover stable genomic regions for traits regardless of the genetic background. This study used information of 483 QTLs for ionomic traits identified from 12 populations for MQTL analysis in Arabidopsis thaliana. The selected QTLs were projected onto the newly constructed genetic consensus map and 33 MQTLs distributed on A. thaliana chromosomes were identified. The average confidence interval (CI) of the drafted MQTLs was 1.30 cM, reduced eight folds from a mean CI of 10.88 cM for the original QTLs. Four MQTLs were considered as stable MQTLs over different genetic backgrounds and environments. In parallel to the gene density over the A. thaliana genome, the genomic distribution of MQTLs over the genetic and physical maps indicated the highest density at non- and sub-telomeric chromosomal regions, respectively. Several candidate genes identified in the MQTLs intervals were associated with ion transportation, tolerance, and homeostasis. The genomic context of the identified MQTLs suggested nine chromosomal regions for Zn, Mn, and Fe control. The QTLs for potassium (K) and phosphorus (P) were the most frequently co-located with Zn (78.3%), Mn (76.2%), and Fe (88.2% and 70.6%) QTLs. The current MQTL analysis demonstrates that meta-QTL analysis is cheaper than, and as informative as genome-wide association study (GWAS) in refining the known QTLs.

Genetic Dissection of Root Angle of Brassica napus in Response to Low Phosphorus

Plant root angle determines the vertical and horizontal distribution of roots in the soil layer, which further influences the acquisition of phosphorus (P) in topsoil. Large genetic variability for the lateral root angle (root angle) was observed in a linkage mapping population (BnaTNDH population) and an association panel of Brassica napus whether at a low P (LP) or at an optimal P (OP). At LP, the average root angle of both populations became smaller. Nine quantitative trait loci (QTLs) at LP and three QTLs at OP for the root angle and five QTLs for the relative root angle (RRA) were identified by the linkage mapping analysis in the BnaTNDH population. Genome-wide association studies (GWASs) revealed 11 single-nucleotide polymorphisms (SNPs) significantly associated with the root angle at LP (LPRA). The interval of a QTL for LPRA on A06 (qLPRA-A06c) overlapped with the confidence region of the leading SNP (Bn-A06-p14439400) significantly associated with LPRA. In addition, a QTL cluster on chromosome C01 associated with the root angle and the primary root length (PRL) in the "pouch and wick" high-throughput phenotyping (HTP) system, the root P concentration in the agar system, and the seed yield in the field was identified in the BnaTNDH population at LP. A total of 87 genes on A06 and 192 genes on C01 were identified within the confidence interval, and 14 genes related to auxin asymmetric redistribution and root developmental process were predicted to be candidate genes. The identification and functional analyses of these genes affecting LPRA are of benefit to the cultivar selection with optimal root system architecture (RSA) under P deficiency in Brassica napus.

Characterization of Plant Growth-Promoting Traits and Inoculation Effects on Triticum durum of Actinomycetes Isolates under Salt Stress Conditions

This study aimed to characterize the halotolerant capability, in vitro, of selected actinomycetes strains and to evaluate their competence in promoting halo stress tolerance in durum wheat in a greenhouse experiment. Fourteen isolates were tested for phosphate solubilization, indole acetic acid, hydrocyanic acid, and ammonia production under different salt concentrations (i.e., 0, 0.25, 0.5, 0.75, 1, 1.25, and 1.5 M NaCl). The presence of 1-aminocyclopropane-1-carboxylate deaminase activity was also investigated. Salinity tolerance was evaluated in durum wheat through plant growth and development parameters: shoot and root length, dry and ash-free dry weight, and the total chlorophyll content, as well as proline accumulation. In vitro assays have shown that the strains can solubilize inorganic phosphate and produce indole acetic acid, hydrocyanic acid, and ammonia under different salt concentrations. Most of the strains (86%) had 1-aminocyclopropane-1-carboxylate deaminase activity, with significant amounts of α-ketobutyric acid. In the greenhouse experiment, inoculation with actinomycetes strains improved the morpho-biochemical parameters of durum wheat plants, which also recorded significantly higher content of chlorophylls and proline than those uninoculated, both under normal and stressed conditions. Our results suggest that inoculation of halotolerant actinomycetes can mitigate the negative effects of salt stress and allow normal growth and development of durum wheat plants.

A high-quality Brassica napus genome reveals expansion of transposable elements, subgenome evolution and disease resistance

Rapeseed (Brassica napus L.) is a recent allotetraploid crop, which is well known for its high oil production. Here, we report a high-quality genome assembly of a typical semi-winter rapeseed cultivar, 'Zhongshuang11' (hereafter 'ZS11'), using a combination of single-molecule sequencing and chromosome conformation capture (Hi-C) techniques. Most of the high-confidence sequences (93.1%) were anchored to the individual chromosomes with a total of 19 centromeres identified, matching the exact chromosome count of B. napus. The repeat sequences in the A and C subgenomes in B. napus expanded significantly from 500 000 years ago, especially over the last 100 000 years. These young and recently amplified LTR-RTs showed dispersed chromosomal distribution but significantly preferentially clustered into centromeric regions. We exhaustively annotated the nucleotide-binding leucine-rich repeat (NLR) gene repertoire, yielding a total of 597 NLR genes in B. napus genome and 17.4% of which are paired (head-to-head arrangement). Based on the resequencing data of 991 B. napus accessions, we have identified 18 759 245 single nucleotide polymorphisms (SNPs) and detected a large number of genomic regions under selective sweep among the three major ecotype groups (winter, semi-winter and spring) in B. napus. We found 49 NLR genes and five NLR gene pairs colocated in selective sweep regions with different ecotypes, suggesting a rapid diversification of NLR genes during the domestication of B. napus. The high quality of our B. napus 'ZS11' genome assembly could serve as an important resource for the study of rapeseed genomics and reveal the genetic variations associated with important agronomic traits.

Genome-wide identification of meiotic recombination hot spots detected by SLAF in peanut (Arachis hypogaea L.)

Recombination hot spots (RHP), caused by meiosis, are considered to play crucial roles in improvement and domestication of crop. Cultivated peanut is one of the most important rich-source of oil and protein crops. However, no direct scale of recombination events and RHP have been estimated for peanut. To examine the scale of recombination events and RHP in peanut, a RIL population with 200 lines and a natural population with 49 cultivars were evaluated. The precise integrated map comprises 4837 SLAF markers with genetic length of 2915.46 cM and density of 1.66 markers per cM in whole genome. An average of 30.0 crossover (2.06 cMMb⁻¹) events was detected per RIL plant. The crossover events (CE) showed uneven distribution among B sub-genome (2.32) and A sub-genome (1.85). There were 4.34% and 7.86% of the genome contained large numbers of CE (> 50 cMMb⁻¹) along chromosomes in F₆ and natural population, respectively. High density of CE regions called RHP, showed negative relationship to marker haplotypes conservative region but positive to heatmap of recombination. The genes located within the RHP regions by GO categories showed the responding of environmental stimuli, which suggested that recombination plays a crucial role in peanut adaptation to changing environments

Genetic dissection of the shoot and root ionomes of Brassica napus grown with contrasting phosphate supplies

BACKGROUND AND AIMS

Mineral elements have many essential and beneficial functions in plants. Phosphorus (P) deficiency can result in changes in the ionomes of plant organs. The aims of this study were to characterize the effects of P supply on the ionomes of shoots and roots, and to identify chromosomal quantitative trait loci (QTLs) for shoot and root ionomic traits, as well as those affecting the partitioning of mineral elements between shoot and root in Brassica napus grown with contrasting P supplies.

METHODS

Shoot and root concentrations of 11 mineral elements (B, Ca, Cu, Fe, K, Mg, Mn, Na, P, S and Zn) were investigated by inductively coupled plasma optical emission spectrometry (ICP-OES) in a Brassica napus double haploid population grown at an optimal (OP) and a low phosphorus supply (LP) in an agar system. Shoot, root and plant contents, and the partitioning of mineral elements between shoot and root were calculated.

KEY RESULTS

The tissue concentrations of B, Ca, Cu, K, Mg, Mn, Na, P and Zn were reduced by P starvation, while the concentration of Fe was increased by P starvation in the BnaTNDH population. A total of 133 and 123 QTLs for shoot and root ionomic traits were identified at OP and LP, respectively. A major QTL cluster on chromosome C07 had a significant effect on shoot Mg and S concentrations at LP and was narrowed down to a 2.1 Mb region using an advanced backcross population.

CONCLUSIONS

The tissue concentration and partitioning of each mineral element was affected differently by P starvation. There was a significant difference in mineral element composition between shoots and roots. Identification of the genes underlying these QTLs will enhance our understanding of processes affecting the uptake and partitioning of mineral elements in Brassica napus.

Two-factor ANOVA of SSH and RNA-seq analysis reveal development-associated Pi-starvation genes in oilseed rape

The 5-leaf-stage rape seedlings were more insensitive to Pi starvation than that of the 3-leaf-stage plants, which may be attributed to the higher expression levels of ethylene signaling and sugar-metabolism genes in more mature seedlings. Traditional suppression subtractive hybridization (SSH) and RNA-Seq usually screen out thousands of differentially expressed genes. However, identification of the most important regulators has not been performed to date. Here, we employed two methods, namely, a two-round SSH and two-factor transcriptome analysis derived from the two-factor ANOVA that is commonly used in the statistics, to identify development-associated inorganic phosphate (Pi) starvation-induced genes in Brassica napus. Several of these genes are related to ethylene signaling (such as EIN3, ACO3, ACS8, ERF1A, and ERF2) or sugar metabolism (such as ACC2, GH3, LHCB1.4, XTH4, and SUS2). Although sucrose and ethylene may counteract each other at the biosynthetic level, they may also work synergistically on Pi-starvation-induced gene expression (such as PT1, PT2, RNS1, ACP5, AT4, and IPS1) and root acid phosphatase activation. Furthermore, three new transcription factors that are responsive to Pi starvation were identified: the zinc-finger MYND domain-containing protein 15 (MYND), a Magonashi family protein (MAGO), and a B-box zinc-finger family salt-tolerance protein. This study indicates that the two methods are highly efficient for functional gene screening in non-model organisms.

Genome-wide selection footprints and deleterious variations in young Asian allotetraploid rapeseed

Brassica napus (AACC, 2n = 38) is an important oilseed crop grown worldwide. However, little is known about the population evolution of this species, the genomic difference between its major genetic groups, such as European and Asian rapeseed, and the impacts of historical large-scale introgression events on this young tetraploid. In this study, we reported the de novo assembly of the genome sequences of an Asian rapeseed (B. napus), Ningyou 7, and its four progenitors and compared these genomes with other available genomic data from diverse European and Asian cultivars. Our results showed that Asian rapeseed originally derived from European rapeseed but subsequently significantly diverged, with rapid genome differentiation after hybridization and intensive local selective breeding. The first historical introgression of B. rapa dramatically broadened the allelic pool but decreased the deleterious variations of Asian rapeseed. The second historical introgression of the double-low traits of European rapeseed (canola) has reshaped Asian rapeseed into two groups (double-low and double-high), accompanied by an increase in genetic load in the double-low group. This study demonstrates distinctive genomic footprints and deleterious SNP (single nucleotide polymorphism) variants for local adaptation by recent intra- and interspecies introgression events and provides novel insights for understanding the rapid genome evolution of a young allopolyploid crop.

Temporal genetic patterns of root growth in Brassica napus L. revealed by a low-cost, high-efficiency hydroponic system

Application of a low-cost and high-efficiency hydroponic system in a rapeseed population verified two types of genetic factors ("persistent" and "stage-specific") that control root development. The root system is a vital plant component for nutrient and water acquisition and is targeted to enhance plant productivity. Genetic dissection of the root system generally focuses on a single stage, but roots grow continuously during plant development. To reveal the temporal genetic patterns of root development, we measured nine root-related traits in a rapeseed recombinant inbred line population at six continuous stages during vegetative growth, using a modified hydroponic system with low-cost and high-efficiency features that could synchronize plant growth under controlled conditions. Phenotypic correlation and growth dynamic analysis suggested the existence of two types of genetic factors ("persistent" and "stage-specific") that control root development. Dynamic (unconditional and conditional) quantitative trait loci (QTL) mapping detected 28 stage-specific and 23 persistent QTLs related to root growth. Among them, 13 early stage-specific, 19 persistent and 8 later stage-specific QTLs were detected at 7 DAS (days after sowing), 16 DAS and 5 EL (expanding leaf stage), respectively, providing efficient and adaptable stages for QTL identification. The effective prediction of biomass accumulation using root morphological traits (up to 96.6% or 92.64% at a specific stage or the final stage, respectively) verified that root growth allocation with maximum root uptake area facilitated biomass accumulation. Furthermore, marker-assistant selection, which combined the "persistent" and "stage-specific" QTLs, proved their effectiveness for root improvement with an excellent uptake area. Our results highlight the potential of high-throughput and precise phenotyping to assess the dynamic genetics of root growth and provide new insights into ideotype root system-based biomass breeding.

Shovelomics for phenotyping root architectural traits of rapeseed/canola (Brassica napus L.) and genome-wide association mapping

Root system in plants plays an important role in mining moisture and nutrients from the soil and is positively correlated to yield in many crops including rapeseed/canola (Brassica napus L.). Substantial phenotypic diversity in root architectural traits among the B. napus growth types leads to a scope of root system improvement in breeding populations. In this study, 216 diverse genotypes were phenotyped for five different root architectural traits following shovelomics approach in the field condition during 2015 and 2016. A single nucleotide polymorphism (SNP) marker panel consisting of 30,262 SNPs was used to conduct genome-wide association study to detect marker/trait association. A total of 31 significant marker loci were identified at 0.01 percentile tail P value cutoff for different root traits. Six marker loci for soil-level taproot diameter (R₁Dia), six loci for belowground taproot diameter (R₂Dia), seven loci for number of primary root branches (PRB), eight loci for root angle, and eight loci for root score (RS) were detected in this study. Several markers associated with root diameters R₁Dia and R₂Dia were also associated with PRB and RS. Significant phenotypic correlation between these traits was observed in both environments. Therefore, taproot diameter appears to be a major determinant of the canola root system architecture and can be used as proxy for other root traits. Fifteen candidate genes related to root traits and root development were detected within 100 kbp upstream and downstream of different significant markers. The identified markers associated with different root architectural traits can be considered for marker-assisted selection for root traits in canola in future.

ZmAPRG, an uncharacterized gene, enhances acid phosphatase activity and Pi concentration in maize leaf during phosphate starvation

An uncharacterized gene, ZmAPRG, isolated by map-based cloning, enhances acid phosphatase activity and phosphate concentration in maize leaf during phosphate starvation. Acid phosphatase (APase) plays important roles in the absorption and utilization of phosphate (Pi) during maize growth. The information on genes regulating the acid phosphatase activity (APA) in maize leaves remains obscured. In a previous study, we delimited the quantitative trait locus, QTL-AP9 for APA to a region of about 546 kb. Here, we demonstrate that the GRMZM2G041022 located in the 546 kb region is a novel acid phosphatase-regulating gene (ZmAPRG). Its overexpression significantly increased the APA and Pi concentration in maize and rice leaves. Subcellular localization of ZmAPRG showed that it was anchored on the plasma and nuclear membrane. The transcriptome analysis of maize ZmAPRG overexpressing lines (ZmAPRG OE) revealed 1287 up-regulated and 392 down-regulated genes. Among these, we found APase, protein phosphatase, and phosphate transporter genes, which are known to be implicated in the metabolism and utilization of Pi. We inferred the ZmAPRG functions as an upstream regulation node, directly or indirectly regulating APases, protein phosphatases, and phosphate transporter genes involved in Pi metabolism and utilization in maize. These findings will pave the way for elucidating the mechanism of APase regulation, absorption and utilization of Pi, and would facilitate maize breeding for efficient use of fertilizers.

Acid phosphatase gene GmHAD1 linked to low phosphorus tolerance in soybean, through fine mapping

KEY MESSAGE

Map-based cloning identified GmHAD1, a gene which encodes a HAD-like acid phosphatase, associated with soybean tolerance to low phosphorus stress. Phosphorus (P) deficiency in soils is a major limiting factor for crop growth worldwide. Plants may adapt to low phosphorus (LP) conditions via changes to root morphology, including the number, length, orientation, and branching of the principal root classes. To elucidate the genetic mechanisms for LP tolerance in soybean, quantitative trait loci (QTL) related to root morphology responses to LP were identified via hydroponic experiments. In total, we identified 14 major loci associated with these traits in a RIL population. The log-likelihood scores ranged from 2.81 to 7.43, explaining 4.23-13.98% of phenotypic variance. A major locus on chromosome 08, named qP8-2, was co-localized with an important P efficiency QTL (qPE8), containing phosphatase genes GmACP1 and GmACP2. Another major locus on chromosome 10 named qP10-2 explained 4.80-13.98% of the total phenotypic variance in root morphology. The qP10-2 contains GmHAD1, a gene which encodes an acid phosphatase. In the transgenic soybean hairy roots, GmHAD1 overexpression increased P efficiency by 8.4-16.5% relative to the control. Transgenic Arabidopsis plants had higher biomass than wild-type plants across both short- and long-term P reduction. These results suggest that GmHAD1, an acid phosphatase gene, improved the utilization of organic phosphate by soybean and Arabidopsis plants.

The high-quality genome of Brassica napus cultivar 'ZS11' reveals the introgression history in semi-winter morphotype

Allotetraploid oilseed rape (Brassica napus L.) is an agriculturally important crop. Cultivation and breeding of B. napus by humans has resulted in numerous genetically diverse morphotypes with optimized agronomic traits and ecophysiological adaptation. To further understand the genetic basis of diversification and adaptation, we report a draft genome of an Asian semi-winter oilseed rape cultivar 'ZS11' and its comprehensive genomic comparison with the genomes of the winter-type cultivar 'Darmor-bzh' as well as two progenitors. The integrated BAC-to-BAC and whole-genome shotgun sequencing strategies were effective in the assembly of repetitive regions (especially young long terminal repeats) and resulted in a high-quality genome assembly of B. napus 'ZS11'. Within a short evolutionary period (~6700 years ago), semi-winter-type 'ZS11' and the winter-type 'Darmor-bzh' maintained highly genomic collinearity. Even so, certain genetic differences were also detected in two morphotypes. Relative to 'Darmor-bzh', both two subgenomes of 'ZS11' are closely related to its progenitors, and the 'ZS11' genome harbored several specific segmental homoeologous exchanges (HEs). Furthermore, the semi-winter-type 'ZS11' underwent potential genomic introgressions with B. rapa (A_r ). Some of these genetic differences were associated with key agronomic traits. A key gene of A03.FLC3 regulating vernalization-responsive flowering time in 'ZS11' was first experienced HE, and then underwent genomic introgression event with A_r , which potentially has led to genetic differences in controlling vernalization in the semi-winter types. Our observations improved our understanding of the genetic diversity of different B. napus morphotypes and the cultivation history of semi-winter oilseed rape in Asia.

Genetic variants associated with the root system architecture of oilseed rape (Brassica napus L.) under contrasting phosphate supply

Breeding crops with ideal root system architecture for efficient absorption of phosphorus is an important strategy to reduce the use of phosphate fertilizers. To investigate genetic variants leading to changes in root system architecture, 405 oilseed rape cultivars were genotyped with a 60K Brassica Infinium SNP array in low and high P environments. A total of 285 single-nucleotide polymorphisms were associated with root system architecture traits at varying phosphorus levels. Nine single-nucleotide polymorphisms corroborate a previous linkage analysis of root system architecture quantitative trait loci in the BnaTNDH population. One peak single-nucleotide polymorphism region on A3 was associated with all root system architecture traits and co-localized with a quantitative trait locus for primary root length at low phosphorus. Two more single-nucleotide polymorphism peaks on A5 for root dry weight at low phosphorus were detected in both growth systems and co-localized with a quantitative trait locus for the same trait. The candidate genes identified on A3 form a haplotype ‘BnA3Hap’, that will be important for understanding the phosphorus/root system interaction and for the incorporation into Brassica napus breeding programs.

Incorporating pleiotropic quantitative trait loci in dissection of complex traits: seed yield in rapeseed as an example

A comprehensive linkage atlas for seed yield in rapeseed. Most agronomic traits of interest for crop improvement (including seed yield) are highly complex quantitative traits controlled by numerous genetic loci, which brings challenges for comprehensively capturing associated markers/genes. We propose that multiple trait interactions underlie complex traits such as seed yield, and that considering these component traits and their interactions can dissect individual quantitative trait loci (QTL) effects more effectively and improve yield predictions. Using a segregating rapeseed (Brassica napus) population, we analyzed a large set of trait data generated in 19 independent experiments to investigate correlations between seed yield and other complex traits, and further identified QTL in this population with a SNP-based genetic bin map. A total of 1904 consensus QTL accounting for 22 traits, including 80 QTL directly affecting seed yield, were anchored to the B. napus reference sequence. Through trait association analysis and QTL meta-analysis, we identified a total of 525 indivisible QTL that either directly or indirectly contributed to seed yield, of which 295 QTL were detected across multiple environments. A majority (81.5%) of the 525 QTL were pleiotropic. By considering associations between traits, we identified 25 yield-related QTL previously ignored due to contrasting genetic effects, as well as 31 QTL with minor complementary effects. Implementation of the 525 QTL in genomic prediction models improved seed yield prediction accuracy. Dissecting the genetic and phenotypic interrelationships underlying complex quantitative traits using this method will provide valuable insights for genomics-based crop improvement.

Breeding histories and selection criteria for oilseed rape in Europe and China identified by genome wide pedigree dissection

Selection breeding has played a key role in the improvement of seed yield and quality in oilseed rape (Brassica napus L.). We genotyped Tapidor (European), Ningyou7 (Chinese) and their progenitors with the Brassica 60 K Illumina Infinium SNP array and mapped a total of 29,347 SNP markers onto the reference genome of Darmor-bzh. Identity by descent (IBD) refers to a haplotype segment of a chromosome inherited from a shared common ancestor. IBDs identified on the C subgenome were larger than those on the A subgenome within both the Tapidor and Ningyou7 pedigrees. IBD number and length were greater in the Ningyou7 pedigree than in the Tapidor pedigree. Seventy nine QTLs for flowering time, seed quality and root morphology traits were identified in the IBDs of Tapidor and Ningyou7. Many more candidate genes had been selected within the Ningyou7 pedigree than within the Tapidor pedigree. These results highlight differences in the transfer of favorable gene clusters controlling key traits during selection breeding in Europe and China.

Genetic Dissection of Root Morphological Traits Related to Nitrogen Use Efficiency in Brassica napus L. under Two Contrasting Nitrogen Conditions

As the major determinant for nutrient uptake, root system architecture (RSA) has a massive impact on nitrogen use efficiency (NUE). However, little is known the molecular control of RSA as related to NUE in rapeseed. Here, a rapeseed recombinant inbred line population (BnaZNRIL) was used to investigate root morphology (RM, an important component for RSA) and NUE-related traits under high-nitrogen (HN) and low-nitrogen (LN) conditions by hydroponics. Data analysis suggested that RM-related traits, particularly root size had significantly phenotypic correlations with plant dry biomass and N uptake irrespective of N levels, but no or little correlation with N utilization efficiency (NUtE), providing the potential to identify QTLs with pleiotropy or specificity for RM- and NUE-related traits. A total of 129 QTLs (including 23 stable QTLs, which were repeatedly detected at least two environments or different N levels) were identified and 83 of them were integrated into 22 pleiotropic QTL clusters. Five RM-NUE, ten RM-specific and three NUE-specific QTL clusters with same directions of additive-effect implied two NUE-improving approaches (RM-based and N utilization-based directly) and provided valuable genomic regions for NUE improvement in rapeseed. Importantly, all of four major QTLs and most of stable QTLs (20 out of 23) detected here were related to RM traits under HN and/or LN levels, suggested that regulating RM to improve NUE would be more feasible than regulating N efficiency directly. These results provided the promising genomic regions for marker-assisted selection on RM-based NUE improvement in rapeseed.

Hybrid Performance of an Immortalized F2 Rapeseed Population Is Driven by Additive, Dominance, and Epistatic Effects

Genomics-based prediction of hybrid performance promises to boost selection gain. The main goal of our study was to investigate the relevance of additive, dominance, and epistatic effects for determining hybrid seed yield in a biparental rapeseed population. We re-analyzed 60,000 SNP array and seed yield data points from an immortalized F₂ population comprised of 318 hybrids and 180 parental lines by performing genome-wide QTL mapping and predictions in combination with five-fold cross-validation. Moreover, an additional set of 37 hybrids were genotyped and phenotyped in an independent environment. The decomposition of the phenotypic variance components and the cross-validated results of the QTL mapping and genome-wide predictions revealed that the hybrid performance in rapeseed was driven by a mix of additive, dominance, and epistatic effects. Interestingly, the genome-wide prediction accuracy in the additional 37 hybrids remained high when modeling exclusively additive effects but was severely reduced when dominance or epistatic effects were also included. This loss in accuracy was most likely caused by more pronounced interactions of environments with dominance and epistatic effects than with additive effects. Consequently, the development of robust hybrid prediction models, including dominance and epistatic effects, required much deeper phenotyping in multi-environmental trials.

Seed Quality Traits Can Be Predicted with High Accuracy in Brassica napus Using Genomic Data

Improving seed oil yield and quality are central targets in rapeseed (Brassica napus) breeding. The primary goal of our study was to examine and compare the potential and the limits of marker-assisted selection and genome-wide prediction of six important seed quality traits of B. napus. Our study is based on a bi-parental population comprising 202 doubled haploid lines and a diverse validation set including 117 B. napus inbred lines derived from interspecific crosses between B. rapa and B. carinata. We used phenotypic data for seed oil, protein, erucic acid, linolenic acid, stearic acid, and glucosinolate content. All lines were genotyped with a 60k SNP array. We performed five-fold cross-validations in combination with linkage mapping and four genome-wide prediction approaches in the bi-parental population. Quantitative trait loci (QTL) with large effects were detected for erucic acid, stearic acid, and glucosinolate content, blazing the trail for marker-assisted selection. Despite substantial differences in the complexity of the genetic architecture of the six traits, genome-wide prediction models had only minor impacts on the prediction accuracies. We evaluated the effects of training population size, marker density and phenotyping intensity on the prediction accuracy. The prediction accuracy in the independent and genetically very distinct validation set still amounted to 0.14 for protein content and 0.17 for oil content reflecting the utility of the developed calibration models even in very diverse backgrounds.