Agronomic and Seed Quality Traits Dissected by Genome-Wide Association Mapping in Brassica napus

Dissection of quantitative trait nucleotides and candidate genes associated with agronomic and yield-related traits under drought stress in rapeseed varieties: integration of genome-wide association study and transcriptomic analysis

Introduction

An important strategy to combat yield loss challenge is the development of varieties with increased tolerance to drought to maintain production. Improvement of crop yield under drought stress is critical to global food security.

Methods

In this study, we performed multiomics analysis in a collection of 119 diverse rapeseed (Brassica napus L.) varieties to dissect the genetic control of agronomic traits in two watering regimes [well-watered (WW) and drought stress (DS)] for 3 years. In the DS treatment, irrigation continued till the 50% pod development stage, whereas in the WW condition, it was performed throughout the whole growing season.

Results

The results of the genome-wide association study (GWAS) using 52,157 single-nucleotide polymorphisms (SNPs) revealed 1,281 SNPs associated with traits. Six stable SNPs showed sequence variation for flowering time between the two irrigation conditions across years. Three novel SNPs on chromosome C04 for plant weight were located within drought tolerance-related gene ABCG16, and their pleiotropically effects on seed weight per plant and seed yield were characterized. We identified the C02 peak as a novel signal for flowering time, harboring 52.77% of the associated SNPs. The 288-kbps LD decay distance analysis revealed 2,232 candidate genes (CGs) associated with traits. The CGs BIG1-D, CAND1, DRG3, PUP10, and PUP21 were involved in phytohormone signaling and pollen development with significant effects on seed number, seed weight, and grain yield in drought conditions. By integrating GWAS and RNA-seq, 215 promising CGs were associated with developmental process, reproductive processes, cell wall organization, and response to stress. GWAS and differentially expressed genes (DEGs) of leaf and seed in the yield contrasting accessions identified BIG1-D, CAND1, and DRG3 genes for yield variation.

Discussion

The results of our study provide insights into the genetic control of drought tolerance and the improvement of marker-assisted selection (MAS) for breeding high-yield and drought-tolerant varieties.

Multi-omics-driven advances in the understanding of triacylglycerol biosynthesis in oil seeds

Vegetable oils are rich sources of polyunsaturated fatty acids and energy as well as valuable sources of human food, animal feed, and bioenergy. Triacylglycerols, which are comprised of three fatty acids attached to a glycerol backbone, are the main component of vegetable oils. Here, we review the development and application of multiple-level omics in major oilseeds and emphasize the progress in the analysis of the biological roles of key genes underlying seed oil content and quality in major oilseeds. Finally, we discuss future research directions in functional genomics research based on current omics and oil metabolic engineering strategies that aim to enhance seed oil content and quality, and specific fatty acids components according to either human health needs or industrial requirements.

Genetic analysis of the variation for mineral accumulation in the leaves and seeds of natural germplasm of Brassica rapa L. (AA) and the its derived forms extracted from an allotetraploid B.juncea L.(AABB)

Brassica rapa L. (2n = 20; AA) is a vegetable and oilseed crop that is grown all over the world. Its leaves, shoots, and seeds store significant amounts of minerals. We used inductively coupled plasma-optical emission spectroscopy (ICP-OES) to determine the concentrations of eleven minerals in the leaves and seeds of 195 advanced generation inbred lines, of which 92 represented natural (NR) B. rapa and the remaining 103 were derived (DR) from a set of mother genotypes originally extracted from an allotetraploid B. juncea (2n = 36; AABB). The inbred lines differed for the composition of leaf and seed minerals. Leaf concentrations of N, K, Zn, and Se were higher in the DR subpanel as compared to NR subpanel, along with high seed accumulations of K and Se. DArT genotyping and genome wide association mapping led to the identification of SNPs associated with leaf and seed mineral compositions. Chromosomes A03, A05, and A10 harboured the most associated loci. Annotations of the regions adjacent to respective GWAS peaks allowed prediction of genes known for acquisition, transport, and accumulation of minerals and heavy metal detoxification. Transcriptome analysis revealed differential expression patterns of the predicted candidates, with most genes either down-regulated in derived genotypes relative to natural forms or their expression being comparable between the two. General downregulation may be a consequence of extracting B. rapa from allotetraploid B. juncea through genome resection. Some of the identified SNPs may be used as DNA markers for breeding programmes designed to modify the leaf and seed mineral compositions.

Homoeologous non-reciprocal translocation explains a major QTL for seed lignin content in oilseed rape (Brassica napus L.)

A homoeologous non-reciprocal translocation was identified in the major QTL for seed lignin content in the low lignin line SGDH14. The lignin biosynthetic gene PAL4 was deleted. Oilseed rape is a major oil crop and a valuable protein source for animal and human nutrition. Lignin is a non-digestible, major component of the seed coat with negative effect on sensory quality, bioavailability and usage of oilseed rape's protein. Hence, seed lignin reduction is of economic and nutritional importance. In this study, the major QTL for reduced lignin content found on chromosome C05 in the DH population SGDH14 x Express 617 was further examined. SGDH14 had lower seed lignin content than Express 617. Harvested seeds from a F2 population of the same cross were additionally field tested and used for seed quality analysis. The F2 population showed a bimodal distribution for seed lignin content. F2 plants with low lignin content had thinner seed coats compared to high lignin lines. Both groups showed a dark seed colour with a slightly lighter colour in the low lignin group indicating that a low lignin content is not necessarily associated with yellow seed colour. Mapping of genomic long-reads from SGDH14 against the Express 617 genome assembly revealed a homoeologous non-reciprocal translocation (HNRT) in the confidence interval of the major QTL for lignin content. A homologous A05 region is duplicated and replaced the C05 region in SGDH14. As consequence several genes located in the C05 region were lost in SGDH14. Thus, a HNRT was identified in the major QTL region for reduced lignin content in the low lignin line SGDH14. The most promising candidate gene related to lignin biosynthesis on C05, PAL4, was deleted.

Meta-analysis of seed weight QTLome using a consensus and highly dense genetic map in Brassica napus L

KEY MESSAGE

We report here the discovery of high-confidence MQTL regions and of putative candidate genes associated with seed weight in B. napus using a highly dense consensus genetic map and by comparing various large-scale multiomics datasets. Seed weight (SW) is a direct determinant of seed yield in Brassica napus and is controlled by many loci. To unravel the main genomic regions associated with this complex trait, we used 13 available genetic maps to construct a consensus and highly dense map, comprising 40,401 polymorphic markers and 9191 genetic bins, harboring a cumulative length of 3047.8 cM. Then, we performed a meta-analysis using 639 projected SW quantitative trait loci (QTLs) obtained from studies conducted since 1999, enabling the identification of 57 meta-QTLS (MQTLs). The confidence intervals of our MQTLs were 9.8 and 4.3 times lower than the average CIs of the original QTLs for the A and C subgenomes, respectively, resulting in the detection of some key genes and several putative novel candidate genes associated with SW. By comparing the genes identified in MQTL intervals with multiomics datasets and coexpression analyses of common genes, we defined a more reliable and shorter list of putative candidate genes potentially involved in the regulation of seed maturation and SW. As an example, we provide a list of promising genes with high expression levels in seeds and embryos (e.g., BnaA03g04230D, BnaC03g08840D, BnaA10g29580D and BnaA03g27410D) that can be more finely studied through functional genetics experiments or that may be useful for MQTL-assisted breeding for SW. The high-density genetic consensus map and the single nucleotide polymorphism (SNP) physical map generated from the latest B. napus cv. Darmor-bzh v10 assembly will be a valuable resource for further mapping and map-based cloning of other important traits.

BnIR: A multi-omics database with various tools for Brassica napus research and breeding

In the post-genome-wide association study era, multi-omics techniques have shown great power and potential for candidate gene mining and functional genomics research. However, due to the lack of effective data integration and multi-omics analysis platforms, such techniques have not still been applied widely in rapeseed, an important oil crop worldwide. Here, we report a rapeseed multi-omics database (BnIR; http://yanglab.hzau.edu.cn/BnIR), which provides datasets of six omics including genomics, transcriptomics, variomics, epigenetics, phenomics, and metabolomics, as well as numerous "variation-gene expression-phenotype" associations by using multiple statistical methods. In addition, a series of multi-omics search and analysis tools are integrated to facilitate the browsing and application of these datasets. BnIR is the most comprehensive multi-omics database for rapeseed so far, and two case studies demonstrated its power to mine candidate genes associated with specific traits and analyze their potential regulatory mechanisms.

Estimation of genetic diversity using seed storage protein (SSP) profiling in wild and cultivated species of Cicer L

BACKGROUND

The narrow genetic diversity of chickpea is a serious impediment to modern cultivar creation. Seed storage proteins (SSPs) are stable and have minimal or no degradation when subjected to isolation and SDS-PAGE.

METHODS AND RESULTS

We have characterized SSPs of 436 chickpea genotypes, belonging to nine annual Cicer species, originated from 47 countries by SDS-PAGE and determined the extent of genetic diversity in chickpea through clustering. Based on scoring, a total of 44 bands (10 to 170 kDa) were identified, which were all polymorphic. The least appeared protein bands were 11, 160 and 170 kDa where band of 11 and 160 kDa was present exclusively in wild type. Five bands were present in < 10% of genotypes. Bands appeared in 200-300 genotypes were suggested less polymorphic, on contrary bands present in 10-150 genotypes were suggested more polymorphic. Polymorphism of protein bands in context to their potential functions reported in literature were explored and suggested that the glubulins were most and glutelins were least abundant, whereas albumins with their known role in stress tolerance can be used as marker in chickpea breeding. Cluster analysis produced 14 clusters, interestingly three clusters contained only Pakistani genotypes and thus Pakistani genotypes appeared as a separate entity from the rest of the genotypes.

CONCLUSION

Our results indicate that SDS-PAGE of SSPs is a powerful technique in determining the genetic diversity plus it is easily adaptable, due to its cost effectiveness in comparison to other genomics tools.

Transcriptome and Regional Association Analyses Reveal the Effects of Oleosin Genes on the Accumulation of Oil Content in Brassica napus

Rapeseed stores lipids in the form of oil bodies. Oil bodies in the seeds of higher plants are surrounded by oleosins. Adjusting oleosin protein levels can prevent the fusion of oil bodies and maintain oil body size during seed development. However, oil contents are affected by many factors, and studies on the complex molecular regulatory mechanisms underlying the variations in seed oil contents of B. napus are limited. In this study, a total of 53 BnOLEO (B. napus oleosin) genes were identified in the genome of B. napus through a genome-wide analysis. The promoter sequences of oleosin genes consisted of various light-, hormone-, and stress-related cis-acting elements, along with transcription factor (TF) binding sites, for 25 TF families in 53 BnOLEO genes. The differentially expressed oleosin genes between two high- and two low-oil-content accessions were explored. BnOLEO3-C09, BnOLEO4-A02, BnOLEO4-A09, BnOLEO2-C04, BnOLEO1-C01, and BnOLEO7-A03 showed higher expressions in the high-oil-content accessions than in low-oil-content accessions, at 25, 35, and 45 days after pollination (DAP) in two different environments. A regional association analysis of 50 re-sequenced rapeseed accessions was used to further analyze these six BnOLEO genes, and it revealed that the nucleotide variations in the BnOLEO1-C01 and BnOLEO7-A03 gene regions were related to the phenotypic variations in seed oil content. Moreover, a co-expression network analysis revealed that the BnOLEO genes were directly linked to lipid/fatty acid metabolism, TF, lipid transport, and carbohydrate genes, thus forming a molecular network involved in seed oil accumulation. These favorable haplotypes can be utilized in molecular marker-assisted selection in order to further improve seed oil contents in rapeseed.

Genome-wide association, RNA-seq and iTRAQ analyses identify candidate genes controlling radicle length of wheat

The radicle, present in the embryo of a seed, is the first root to emerge at germination, and its rapid growth is essential for establishment and survival of the seedling. However, there are few studies on the critical mechanisms underlying radicle and then radicle length in wheat seedlings, despite its importance as a food crop throughout the world. In the present study, 196 wheat accessions from the Huanghuai Wheat Region were screened to measure radicle length under 4 hydroponic culture environments over 3 years. Different expression genes and proteins (DEGs/DEPs) between accessions with extremely long [Yunong 949 (WRL1), Zhongyu 9,302 (WRL2)] and short roots [Yunong 201 (WRS1), Beijing 841 (WRS2)] were identified in 12 sets of root tissue samples by RNA-seq and iTRAQ (Isobaric tags for relative and absolute quantification). Phenotypic results showed that the elongation zone was significantly longer in root accessions with long roots compared to the short-rooted accessions. A genome-wide association study (GWAS) identified four stable chromosomal regions significantly associated with radicle length, among which 1A, 4A, and 7A chromosomes regions explained 7.17% to12.93% of the phenotypic variation. The omics studies identified the expression patterns of 24 DEGs/DEPs changed at both the transcriptional and protein levels. These DEGs/DEPs were mainly involved in carbon fixation in photosynthetic organisms, photosynthesis and phenylpropanoid biosynthesis pathways. TraesCS1A02G104100 and TraesCS2B02G519100 were involved in the biosynthesis of tricin-lignins in cell walls and may affect the extension of cell walls in the radicle elongation zone. A combination of GWAS and RNA-seq analyses revealed 19 DEGs with expression changes in the four accessions, among which, TraesCS1A02G422700 (a cysteine-rich receptor-like protein kinase 6, CRK6) also showed upregulation in the comparison group by RNA-seq, iTRAQ, and qRT-PCR. BSMV-mediated gene silencing also showed that TaCRK6 improves root development in wheat. Our data suggest that TaCRK6 is a candidate gene regulating radicle length in wheat.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Genomic Regions Associated With Seed Meal Quality Traits in Brassica napus Germplasm

The defatted Brassica napus (rapeseed) meal can be high-protein feed for livestock as the protein value of rapeseed meal is higher than that of the majority of other vegetable proteins. Extensive work has already been carried out on developing canola rapeseed where the focus was on reducing erucic acid and glucosinolate content, with less consideration to other antinutritional factors such as tannin, phytate, sinapine, crude fiber, etc. The presence of these antinutrients limits the use and marketing of rapeseed meals and a significant amount of it goes unused and ends up as waste. We investigated the genetic architecture of crude protein, methionine, tryptophan, total phenols, β-carotene, glucosinolates (GLSs), phytate, tannins, sinapine, and crude fiber content of defatted seed meal samples by conducting a genome-wide association study (GWAS), using a diversity panel comprising 96 B. napus genotypes. Genotyping by sequencing was used to identify 77,889 SNPs, spread over 19 chromosomes. Genetic diversity and phenotypic variations were generally high for the studied traits. A total of eleven genotypes were identified which showed high-quality protein, high antioxidants, and lower amount of antinutrients. A significant negative correlation between protein and limiting amino acids and a significant positive correlation between GLS and phytic acid were observed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits and quantile-quantile (QQ) plots were generated to allow the best-fit algorithm. Annotation of genomic regions around associated SNPs helped to predict various trait-related candidates such as ASP2 and EMB1027 (amino acid biosynthesis); HEMA2, GLU1, and PGM (tryptophan biosynthesis); MS3, CYSD1, and MTO1 (methionine biosynthesis); LYC (β-carotene biosynthesis); HDR and ISPF (MEP pathway); COS1 (riboflavin synthesis); UGT (phenolics biosynthesis); NAC073 (cellulose and hemicellulose biosynthesis); CYT1 (cellulose biosynthesis); BGLU45 and BGLU46 (lignin biosynthesis); SOT12 and UGT88A1 (flavonoid pathway); and CYP79A2, DIN2, and GSTT2 (GLS metabolism), etc. The functional validation of these candidate genes could confirm key seed meal quality genes for germplasm enhancement programs directed at improving protein quality and reducing the antinutritional components in B. napus.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

SNP- and Haplotype-Based GWAS of Flowering-Related Traits in Brassica napus

Traits related to flowering time are the most promising agronomic traits that directly impact the seed yield and oil quality of rapeseed (Brassica napus L.). Developing early flowering and maturity rapeseed varieties is an important breeding objective in B. napus. Many studies have reported on days to flowering, but few have reported on budding, bolting, and the interval between bolting and DTF. Therefore, elucidating the genetic architecture of QTLs and genes regulating flowering time, we presented an integrated investigation on SNP and haplotype-based genome-wide association study of 373 diverse B. napus germplasm, which were genotyped by the 60K SNP array and were phenotyped in the four environments. The results showed that a total of 15 and 37 QTLs were detected from SNP and haplotype-based GWAS, respectively. Among them, seven QTL clusters were identified by haplotype-based GWAS. Moreover, three and eight environmentally stable QTLs were detected by SNP-GWAS and haplotype-based GWAS, respectively. By integrating the above two approaches and by co-localizing the four traits, ten (10) genomic regions were under selection on chromosomes A03, A07, A08, A10, C06, C07, and C08. Interestingly, the genomic regions FT.A07.1, FT.A08, FT.C06, and FT.C07 were identified as novel. In these ten regions, a total of 197 genes controlling FT were detected, of which 14 highly expressed DEGs were orthologous to 13 Arabidopsis thaliana genes after integration with transcriptome results. In a nutshell, the above results uncovered the genetic architecture of important agronomic traits related to flowering time and provided a basis for multiple molecular marker-trait associations in B. napus.

Plant monounsaturated fatty acids: Diversity, biosynthesis, functions and uses

Monounsaturated fatty acids are straight-chain aliphatic monocarboxylic acids comprising a unique carbon‑carbon double bond, also termed unsaturation. More than 50 distinct molecular structures have been described in the plant kingdom, and more remain to be discovered. The evolution of land plants has apparently resulted in the convergent evolution of non-homologous enzymes catalyzing the dehydrogenation of saturated acyl chain substrates in a chemo-, regio- and stereoselective manner. Contrasted enzymatic characteristics and different subcellular localizations of these desaturases account for the diversity of existing fatty acid structures. Interestingly, the location and geometrical configuration of the unsaturation confer specific characteristics to these molecules found in a variety of membrane, storage, and surface lipids. An ongoing research effort aimed at exploring the links existing between fatty acid structures and their biological functions has already unraveled the importance of several monounsaturated fatty acids in various physiological and developmental contexts. What is more, the monounsaturated acyl chains found in the oils of seeds and fruits are widely and increasingly used in the food and chemical industries due to the physicochemical properties inherent in their structures. Breeders and plant biotechnologists therefore develop new crops with high monounsaturated contents for various agro-industrial purposes.

The Bnapus50K array: a quick and versatile genotyping tool for Brassica napus genomic breeding and research

Rapeseed is a globally cultivated commercial crop, primarily grown for its oil. High-density single nucleotide polymorphism (SNP) arrays are widely used as a standard genotyping tool for rapeseed research, including for gene mapping, genome-wide association studies, germplasm resource analysis, and cluster analysis. Although considerable rapeseed genome sequencing data have been released, DNA arrays are still an attractive choice for providing additional genetic data in an era of high-throughput whole-genome sequencing. Here, we integrated re-sequencing DNA array data (32,216, 304 SNPs) from 505 inbred rapeseed lines, allowing us to develop a sensitive and efficient genotyping DNA array, Bnapus50K, with a more consistent genetic and physical distribution of probes. A total of 42,090 high-quality probes were filtered and synthesized, with an average distance between adjacent SNPs of 8 kb. To improve the practical application potential of this array in rapeseed breeding, we also added 1,618 functional probes related to important agronomic traits such as oil content, disease resistance, male sterility, and flowering time. The additional probes also included those specifically for detecting genetically modified material. These probes show a good detection efficiency and are therefore useful for gene mapping, along with crop variety improvement and identification. The novel Bnapus50K DNA array developed in this study could prove to be a quick and versatile genotyping tool for B. napus genomic breeding and research.

Identification and Development of KASP Markers for Novel Mutant BnFAD2 Alleles Associated With Elevated Oleic Acid in Brassica napus

The fatty acid desaturase FAD2 genes are the main contributors to oleic acid content, and different FAD2 alleles can result in different oleic acid contents in rapeseed oil. Hence, identification of allelic variation in FAD2 is an extremely desirable breeding goal. By performing QTL mapping using 190 F_2:3 lines genotyped by genome-wide single nucleotide polymorphism (SNP) markers assayed by the Brassica 60 K Infinium BeadChip Array, four quantitative trait loci (QTL) for C18:1 content were mapped on chromosomes A01, A05, A09 and C05 over 3 years in a population segregating for oleic acid content. Two BnFAD2 genes on A05 and C05 were anchored within the QTL intervals, explaining 45-52 and 15-44% of the observed variation for C18:1 content. Sequence polymorphisms between the corresponding coding regions of the parental lines found two single-nucleotide polymorphisms (SNPs) in BnFAD2.A05 and BnFAD2.C05, respectively, which led to the amino acid changes (C421T and G1073E) in the corresponding proteins. The mutation sites of Bnfad2.A05 and Bnfad2.C05 alleles were located within the second H-box and near the third H-box motif of the protein, respectively, and were found to be novel mutant alleles. Lines resulting from the combination of these two alleles contained up to 88% oleic acid in their seed oil, compared with 63% in wild-type controls. Two competitive allele-specific PCR (KASP) markers based on these two mutation sites were successfully developed and validated in segregating F₂ populations. These markers will facilitate breeding for ultra-high seed oleic acid content in oilseed rape.

Mapping of seed quality traits in the C genome of Brassica napus by using a population carrying genome content of B. oleracea and their effect on other traits

Increasing seed oil and protein contents and reducing the content of seed glucosinolates (GSLs) in Brassica oilseed crops are important objectives in breeding. By using an oilseed rape (B. napus L.) doubled-haploid (DH) population carrying genome content introgressed from Chinese kale (B. oleracea L.), we mapped quantitative trait loci (QTL) for these seed quality traits and investigated their effect on other traits including seed yield. A stable QTL for seed oil content was identified on chromosome C5 at 40-42 Mb position and a QTL for seed GSL content was identified on C9 at 7-8 Mb position. The C5 and C9 QTL alleles for high oil and GSL contents were derived from Chinese kale, demonstrating that high-oil QTL allele can be found in the parental species of oilseed rape. The low-GSL QTL allele of C9 exerted a significant positive effect on seed protein content, demonstrating that selection for this QTL allele contributed to higher protein content in canola seed. These two QTL were not affected by field environmental conditions and did not exert a significant effect on days to flowering and seed yield. Thus, the genomic regions and the molecular markers identified in this study should be useful in molecular breeding of the seed quality traits in oilseed rape.

Genetics of days to flowering, maturity and plant height in natural and derived forms of Brassica rapa L

Genome wide association studies enabled prediction of many candidate genes for flowering, maturity and plant height under differing day-length conditions. Some genes were envisaged only from derived B. rapa. Flowering and plant height are the key life history traits. These are crucial for adaptation and productivity. Current investigations aimed to examine genotypic differences governing days to flowering, maturity and plant height under contrasting day-length conditions; and identify genomic regions governing the observed phenotypic variations. An association panel comprising 195 inbred lines, representing natural (NR) and derived (DR) forms of Brassica rapa (AA; 2n = 20), was evaluated at two sowing dates and two locations, representing different day-length regimes. Derived B. rapa is a unique pre-breeding material extracted from B. juncea (AABB; 2n = 36). Population structure analysis, using DArT genotypes established derived B. rapa as a genetic resource distinct from natural B. rapa. Genome wide association studies facilitated detection of many trait associated SNPs. Chromosomes A03, A05 and A09 harboured majority of these. Functional annotation of the associated SNPs and surrounding genome space(s) helped to predict 43 candidate genes. Many of these were predicted under specific day-length conditions. Important among these were the genes encoding floral meristem identity (SPL3, SPL15, AP3, BAM2), photoperiodic responses (COL2, AGL18, SPT, NF-YC4), gibberellic acid biosynthesis (GA1) and regulation of flowering (EBS). Some of the predicted genes were detected for DR subpanel alone. Genes controlling hormones, auxins and gibberellins appeared important for the regulation of plant height. Many of the significant SNPs were located on chromosomes harbouring previously reported QTLs and candidate genes. The identified loci may be used for marker-assisted selection after due validation.

Enhancing the Nutritional Quality of Major Food Crops Through Conventional and Genomics-Assisted Breeding

Nutritional stress is making over two billion world population malnourished. Either our commercially cultivated varieties of cereals, pulses, and oilseed crops are deficient in essential nutrients or the soils in which these crops grow are becoming devoid of minerals. Unfortunately, our major food crops are poor sources of micronutrients required for normal human growth. To overcome the problem of nutritional deficiency, greater emphasis should be laid on the identification of genes/quantitative trait loci (QTLs) pertaining to essential nutrients and their successful deployment in elite breeding lines through marker-assisted breeding. The manuscript deals with information on identified QTLs for protein content, vitamins, macronutrients, micro-nutrients, minerals, oil content, and essential amino acids in major food crops. These QTLs can be utilized in the development of nutrient-rich crop varieties. Genome editing technologies that can rapidly modify genomes in a precise way and will directly enrich the nutritional status of elite varieties could hold a bright future to address the challenge of malnutrition.

Association Mapping of Seed Quality Traits Under Varying Conditions of Nitrogen Application in Brassica juncea L. Czern & Coss

Indian mustard (Brassica juncea) is a major source of vegetable oil in the Indian subcontinent. The seed cake left after the oil extraction is used as livestock feed. We examined the genetic architecture of oil, protein, and glucosinolates by conducting a genome-wide association study (GWAS), using an association panel comprising 92 diverse genotypes. We conducted trait phenotyping over 2 years at two levels of nitrogen (N) application. Genotyping by sequencing was used to identify 66,835 loci, covering 18 chromosomes. Genetic diversity and phenotypic variations were high for the studied traits. Trait performances were stable when averaged over years and N levels. However, individual performances differed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits. Population structure, principal components (PCs) analysis, and discriminant analysis of principal components (DAPCs) were included as covariates to overcome the bias due to the population stratification. We identified 16, 23, and 27 loci associated with oil, protein, and glucosinolates, respectively. We also established LD patterns and haplotype structures for the candidate genes. The average block sizes were larger on A-genome chromosomes as compared to the B- genome chromosomes. Genetic associations differed over N levels. However, meta-analysis of GWAS datasets not only improved the power to recognize associations but also helped to identify common SNPs for oil and protein contents. Annotation of the genomic region around the identified SNPs led to the prediction of 21 orthologs of the functional candidate genes related to the biosynthesis of oil, protein, and glucosinolates. Notable among these are: LACS5 (A09), FAD6 (B05), ASN1 (A06), GTR2 (A06), CYP81G1 (B06), and MYB44 (B06). The identified loci will be very useful for marker-aided breeding for seed quality modifications in B. juncea.

Towards model-driven characterization and manipulation of plant lipid metabolism

Plant lipids have versatile applications and provide essential fatty acids in human diet. Therefore, there has been a growing interest to better characterize the genetic basis, regulatory networks, and metabolic pathways that shape lipid quantity and composition. Addressing these issues is challenging due to context-specificity of lipid metabolism integrating environmental, developmental, and tissue-specific cues. Here we systematically review the known metabolic pathways and regulatory interactions that modulate the levels of storage lipids in oilseeds. We argue that the current understanding of lipid metabolism provides the basis for its study in the context of genome-wide plant metabolic networks with the help of approaches from constraint-based modeling and metabolic flux analysis. The focus is on providing a comprehensive summary of the state-of-the-art of modeling plant lipid metabolic pathways, which we then contrast with the existing modeling efforts in yeast and microalgae. We then point out the gaps in knowledge of lipid metabolism, and enumerate the recent advances of using genome-wide association and quantitative trait loci mapping studies to unravel the genetic regulations of lipid metabolism. Finally, we offer a perspective on how advances in the constraint-based modeling framework can propel further characterization of plant lipid metabolism and its rational manipulation.

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.

Dissection of genetic architecture for glucosinolate accumulations in leaves and seeds of Brassica napus by genome-wide association study

Glucosinolates (GSLs), whose degradation products have been shown to be increasingly important for human health and plant defence, compose important secondary metabolites found in the order Brassicales. It is highly desired to enhance pest and disease resistance by increasing the leaf GSL content while keeping the content low in seeds of Brassica napus, one of the most important oil crops worldwide. Little is known about the regulation of GSL accumulation in the leaves. We quantified the levels of 9 different GSLs and 15 related traits in the leaves of 366 accessions and found that the seed and leaf GSL content were highly correlated (r = 0.79). A total of 78 loci were associated with GSL traits, and five common and eleven tissue-specific associated loci were related to total leaf and seed GSL content. Thirty-six candidate genes were inferred to be involved in GSL biosynthesis. The candidate gene BnaA03g40190D (BnaA3.MYB28) was validated by DNA polymorphisms and gene expression analysis. This gene was responsible for high leaf/low seed GSL content and could explain 30.62% of the total leaf GSL variation in the low seed GSL panel and was not fixed during double-low rapeseed breeding. Our results provide new insights into the genetic basis of GSL variation in leaves and seeds and may facilitate the metabolic engineering of GSLs and the breeding of high leaf/low seed GSL content in B. napus.

Strong temporal dynamics of QTL action on plant growth progression revealed through high-throughput phenotyping in canola

A major challenge of plant biology is to unravel the genetic basis of complex traits. We took advantage of recent technical advances in high-throughput phenotyping in conjunction with genome-wide association studies to elucidate genotype-phenotype relationships at high temporal resolution. A diverse Brassica napus population from a commercial breeding programme was analysed by automated non-invasive phenotyping. Time-resolved data for early growth-related traits, including estimated biovolume, projected leaf area, early plant height and colour uniformity, were established and complemented by fresh and dry weight biomass. Genome-wide SNP array data provided the framework for genome-wide association analyses. Using time point data and relative growth rates, multiple robust main effect marker-trait associations for biomass and related traits were detected. Candidate genes involved in meristem development, cell wall modification and transcriptional regulation were detected. Our results demonstrate that early plant growth is a highly complex trait governed by several medium and many small effect loci, most of which act only during short phases. These observations highlight the importance of taking the temporal patterns of QTL/allele actions into account and emphasize the need for detailed time-resolved analyses to effectively unravel the complex and stage-specific contributions of genes affecting growth processes that operate at different developmental phases.

Integrating GWAS, linkage mapping and gene expression analyses reveals the genetic control of growth period traits in rapeseed (Brassica napus L.)

BACKGROUND

Brassica napus is one of the most important oilseed crops, and also an important biofuel plant due to its low air pollution and renewability. Growth period are important traits that affect yield and are crucial for its adaptation to different environments in B. napus.

RESULTS

To elucidate the genetic basis of growth period traits, genome-wide association analysis (GWAS) and linkage mapping were employed to detect the quantitative trait loci (QTL) for days to initial flowering (DIF), days to final flowering (DFF), flowering period (FP), maturity time (MT), and whole growth period (GP). A total of 146 SNPs were identified by association mapping, and 83 QTLs were identified by linkage mapping using the RIL population. Among these QTLs, 19 were pleiotropic SNPs related to multiple traits, and six (q18DFF.A03-2, q18MT.A03-2, q17DFF.A05-1, q18FP.C04, q17DIF.C05 and q17GP.C09) were consistently detected using both mapping methods. Additionally, we performed RNA sequencing to analyze the differential expression of gene (DEG) transcripts between early- and late-flowering lines selected from the RIL population, and the DEGs were integrated with association mapping and linkage analysis to confirm their roles in the growth period. Consequently, 12 candidate genes associated with growth period traits were identified in B. napus. Among these genes, seven have polymorphic sites in the coding sequence and the upstream 2-kb sequence based on the resequencing data. The haplotype BnaSOC1.A05-Haplb and BnaLNK2.C06-Hapla showed more favorable phenotypic traits.

CONCLUSIONS

The candidate genes identified in this study will contribute to our genetic understanding of growth period traits and can be used as targets for target mutations or marker-assisted breeding for rapeseed adapted to different environments.

A novel quantitative trait locus on chromosome A9 controlling oleic acid content in Brassica napus

One of the most important goals in the breeding of oilseed crops, including Brassica napus, is to improve the quality of edible vegetable oil, which is mainly determined by the seed fatty acid composition, particularly the C18:1 content. Previous studies have indicated that the C18:1 content is a polygenic trait, and no stable quantitative trait loci (QTLs) except for FAD2 have been reported. By performing a GWAS using 375 low erucic acid B. napus accessions genotyped with the Brassica 60K SNP array and constructing a high-density SNP-based genetic map of a 150 DH population, we identified a novel QTL on the A9 chromosome. The novel locus could explain 11.25%, 5.72% and 6.29% of phenotypic variation during three consecutive seasons and increased the C18:1 content by approximately 3%-5%. By fine mapping and gene expression analysis, we found three potential candidate genes and verified the fatty acids in a homologous gene mutant of Arabidopsis. A metal ion-binding protein was found to be the most likely candidate gene in the region. Thus, the C18:1 content can be further increased to about 80% with this novel locus together with FAD2 mutant allele without compromise of agronomic performance. A closely linked marker, BnA129, for this novel QTL (OLEA9) was developed so that we can effectively identify materials with high C18:1 content at an early growth stage by marker-assisted selection. Our results may also provide new insight for understanding the complex genetic mechanism of fatty acid metabolism.

Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus

BACKGROUND

Increasing seed oil content is one of the most important targets for rapeseed (Brassica napus) breeding. However, genetic mechanisms of mature seed oil content in Brassica napus (B. napus) remain little known. To identify oil content-related genes, a genome-wide association study (GWAS) was performed using 588 accessions.

RESULTS

High-throughput genome resequencing resulted in 385,692 high-quality single nucleotide polymorphism (SNPs) with a minor allele frequency (MAF) > 0.05. We identified 17 loci that were significantly associated with seed oil content, among which 12 SNPs were distributed on the A3 (11 loci) and A1 (one loci) chromosomes, and five novel significant SNPs on the C5 (one loci) and C7 (four loci) chromosomes, respectively. Subsequently, we characterized differentially expressed genes (DEGs) between the seeds and silique pericarps on main florescences and primary branches of extremely high- and low-oil content accessions (HO and LO). A total of 64 lipid metabolism-related DEGs were identified, 14 of which are involved in triacylglycerols (TAGs) biosynthesis and assembly. Additionally, we analyzed differences in transcription levels of key genes involved in de novo fatty acid biosynthesis in the plastid, TAGs assembly and lipid droplet packaging in the endoplasmic reticulum (ER) between high- and low-oil content B. napus accessions.

CONCLUSIONS

The combination of GWAS and transcriptome analyses revealed seven candidate genes located within the confidence intervals of significant SNPs. Current findings provide valuable information for facilitating marker-based breeding for higher seed oil content in B. napus.

Genetic dissection of grain iron and zinc concentrations in lentil (Lens culinaris Medik.)

Iron (Fe) and zinc (Zn) deficiencies are wide spread in South Asia and Africa. Biofortification of food crops is a viable means of addressing micronutrient deficiencies. Lentil is an important pulse crop that provides affordable source of proteins, minerals, fibre and carbohydrates for micronutrient deficient countries. An association mapping (AM) panel of 96 diverse lentil genotypes from India and Mediterranean region was evaluated for three seasons and genotyped using 80 polymorphic simple-sequence repeat (SSR) markers for identification of the markers associated with grain Fe and Zn concentrations. A Bayesian model based clustering identified five subpopulations, adequately explaining the genetic structure of the AM panel. The linkage disequilibrium (LD) analysis using mixed linear model (MLM) identified two SSR markers, GLLC106 and GLLC108, associated with grain Fe concentration explaining 17% and 6% phenotypic variation, respectively and three SSR markers (PBALC 364, PBALC 92 and GLLC592) associated with grain Zn concentration, explaining 6%, 8% and 13% phenotypic variation, respectively. The identified SSRs exhibited consistent performance across three seasons and have potential for utilization in lentil molecular breeding programme.

Identification of SNP loci and candidate genes related to four important fatty acid composition in Brassica napus using genome wide association study

Rapeseed oil (canola, Brassica napus L.) is an important healthy vegetable oil throughout the world, the nutritional and economical value of which largely depends on its seed fatty acid composition. In this study, based on 201,187 SNP markers developed from the SLAF-seq (specific locus amplified fragment sequencing), a genome wide association study of four important fatty acid content traits (erucic acid, oleic acid, linoleic acid and linolenic acid) in a panel of 300 inbred lines of rapeseed in two environments (JXAU and JXRIS) was carried out. A total of 148 SNP loci significantly associated with these traits were detected by MLM model analysis respectively, and 30 SNP loci on A08 and C03 chromosomes were detected in three traits of erucic acid, oleic acid and linoleic acid contents simultaneously. Furthermore, 108 highly favorable alleles for increasing oleic acid and linoleic acid content, also for decreasing erucic acid content simultaneously were observed. By a basic local alignment search tool (BLAST) search with in a distance of 100 Kb around these significantly SNP-trait associations, we identified 20 orthologs of the functional candidate genes related to fatty acid biosynthesis, including the known vital fatty acid biosynthesis genes of BnaA.FAE1 and BnaC. FAE1 on the A08 and C03 chromosomes, and other potential candidate genes involving in the fatty acid biosynthesis pathway, such as the orthologs genes of FAD2, LACS09, KCS17, CER4, TT16 and ACBP5. This study lays a basis for uncovering the genetic variations and the improvement of fatty acid composition in B. napus.

Effect of breeding on nitrogen use efficiency-associated traits in oilseed rape

Oilseed rape is one of the most important dicotyledonous field crops in the world, where it plays a key role in productive cereal crop rotations. However, its production requires high nitrogen fertilization and its nitrogen footprint exceeds that of most other globally important crops. Hence, increased nitrogen use efficiency (NUE) in this crop is of high priority for sustainable agriculture. We report a comprehensive study of macrophysiological characteristics associated with breeding progress, conducted under contrasting nitrogen fertilization levels in a large panel of elite oilseed rape varieties representing breeding progress over the past 20 years. The results indicate that increased plant biomass at flowering, along with increases in primary yield components, have increased NUE in modern varieties. Nitrogen uptake efficiency has improved through breeding, particularly at high nitrogen. Despite low heritability, the number of seeds per silique is associated positively with increased yield in modern varieties. Seed weight remains unaffected by breeding progress; however, recent selection for high seed oil content and for high seed yields appears to have promoted a negative correlation (r= -0.39 at high and r= -0.49 at low nitrogen) between seed weight and seed oil concentration. Overall, our results reveal valuable breeding targets to improve NUE in oilseed rape.

Systems biology of seeds: decoding the secret of biochemical seed factories for nutritional security

Seeds serve as biochemical factories of nutrition, processing, bio-energy and storage related important bio-molecules and act as a delivery system to transmit the genetic information to the next generation. The research pertaining towards delineating the complex system of regulation of genes and pathways related to seed biology and nutrient partitioning is still under infancy. To understand these, it is important to know the genes and pathway(s) involved in the homeostasis of bio-molecules. In recent past with the advent and advancement of modern tools of genomics and genetic engineering, multi-layered 'omics' approaches and high-throughput platforms are being used to discern the genes and proteins involved in various metabolic, and signaling pathways and their regulations for understanding the molecular genetics of biosynthesis and homeostasis of bio-molecules. This can be possible by exploring systems biology approaches via the integration of omics data for understanding the intricacy of seed development and nutrient partitioning. These information can be exploited for the improvement of biologically important chemicals for large-scale production of nutrients and nutraceuticals through pathway engineering and biotechnology. This review article thus describes different omics tools and other branches that are merged to build the most attractive area of research towards establishing the seeds as biochemical factories for human health and nutrition.

Association Between SSR Molecular Markers and Chemical and Sensory Traits of Cacao Samples Using Multiple Regression Analysis

The sensory properties of cacao beans are linked to the chemical composition of the seeds, and both characteristics are the partial results of its allelic composition. Therefore, it is useful to search for molecular markers associated with these traits. We perform multiple regression analysis to associate previously generated data of alleles generated with 12 SSR (of cultivated cacao trees) with data obtained from chemical and sensory characterization (of beans) of plants grown in the southern region from Mexico. When the association was significant, the mathematical models for predictive purposes were proposed. All phenotypic traits evaluated showed equations with setting values R² > 0.5. All chemical characters tested have a significant association with at least two alleles (P < 0.05). In addition, the fat content was associated with six molecular markers (mTcCIR03₂₀₉, mTcCIR12₁₈₈, mTcCIR19₂₈₆, mTcCIR07₁₅₀, mTcCIR19₃₁₀). The most common allele was mTcCIR12₁₈₈, which was associated with the contents of eicosanoic acid, moisture, fat and total polyphenols content. The mTcCIR28₃₆₂ allele is associated with sensory characters bitterness, musty odor, and roasted odor. These alleles could be useful as molecular markers of chemical and sensory characteristics of cacao samples.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Genetic changes in a novel breeding population of Brassica napus synthesized from hundreds of crosses between B. rapa and B. carinata

Introgression of genomic variation between and within related crop species is a significant evolutionary approach for population differentiation, genome reorganization and trait improvement. Using the Illumina Infinium Brassica 60K SNP array, we investigated genomic changes in a panel of advanced generation new-type Brassica napus breeding lines developed from hundreds of interspecific crosses between 122 Brassica rapa and 74 Brassica carinata accessions, and compared them with representative accessions of their three parental species. The new-type B. napus population presented rich genetic diversity and abundant novel genomic alterations, consisting of introgressions from B. rapa and B. carinata, novel allelic combinations, reconstructed linkage disequilibrium patterns and haplotype blocks, and frequent deletions and duplications (nonrandomly distributed), particularly in the C subgenome. After a much shorter, but very intensive, selection history compared to traditional B. napus, a total of 15 genomic regions with strong selective sweeps and 112 genomic regions with putative signals of selective sweeps were identified. Some of these regions were associated with important agronomic traits that were selected for during the breeding process, while others were potentially associated with restoration of genome stability and fertility after interspecific hybridization. Our results demonstrate how a novel method for population-based crop genetic improvement can lead to rapid adaptation, restoration of genome stability and positive responses to artificial selection.

Integrated QTL and eQTL Mapping Provides Insights and Candidate Genes for Fatty Acid Composition, Flowering Time, and Growth Traits in a F2 Population of a Novel Synthetic Allopolyploid Brassica napus

Brassica napus (B. napus, AACC), is an economically important allotetraploid crop species that resulted from hybridization between two diploid species, Brassica rapa (AA) and Brassica olereacea (CC). We have created one new synthetic B. napus genotype Da-Ae (AACC) and one introgression line Da-Ol-1 (AACC), which were used to generate an F₂ mapping population. Plants in this F₂ mapping population varied in fatty acid content, flowering time, and growth-related traits. Using quantitative trait locus (QTL) mapping, we aimed to determine if Da-Ae and Da-Ol-1 provided novel genetic variation beyond what has already been found in B. napus. Making use of the genotyping information generated from RNA-seq data of these two lines and their F₂ mapping population of 166 plants, we constructed a genetic map consisting of 2,021 single nucleotide polymorphism markers that spans 2,929 cM across 19 linkage groups. Besides the known major QTL identified, our high resolution genetic map facilitated the identification of several new QTL contributing to the different fatty acid levels, flowering time, and growth-related trait values. These new QTL probably represent novel genetic variation that existed in our new synthetic B. napus strain. By conducting genome-wide expression variation analysis in our F₂ mapping population, genetic regions that potentially regulate many genes across the genome were revealed. A FLOWERING LOCUS C gene homolog, which was identified as a candidate regulating flowering time and multiple growth-related traits, was found underlying one of these regions. Integrated QTL and expression QTL analyses also helped us identified candidate causative genes associated with various biological traits through expression level change and/or possible protein function modification.

Genome-wide regression models considering general and specific combining ability predict hybrid performance in oilseed rape with similar accuracy regardless of trait architecture

KEY MESSAGE

Genomic prediction using the Brassica 60 k genotyping array is efficient in oilseed rape hybrids. Prediction accuracy is more dependent on trait complexity than on the prediction model. In oilseed rape breeding programs, performance prediction of parental combinations is of fundamental importance. Due to the phenomenon of heterosis, per se performance is not a reliable indicator for F₁-hybrid performance, and selection of well-paired parents requires the testing of large quantities of hybrid combinations in extensive field trials. However, the number of potential hybrids, in general, dramatically exceeds breeding capacity and budget. Integration of genomic selection (GS) could substantially increase the number of potential combinations that can be evaluated. GS models can be used to predict the performance of untested individuals based only on their genotypic profiles, using marker effects previously predicted in a training population. This allows for a preselection of promising genotypes, enabling a more efficient allocation of resources. In this study, we evaluated the usefulness of the Illumina Brassica 60 k SNP array for genomic prediction and compared three alternative approaches based on a homoscedastic ridge regression BLUP and three Bayesian prediction models that considered general and specific combining ability (GCA and SCA, respectively). A total of 448 hybrids were produced in a commercial breeding program from unbalanced crosses between 220 paternal doubled haploid lines and five male-sterile testers. Predictive ability was evaluated for seven agronomic traits. We demonstrate that the Brassica 60 k genotyping array is an adequate and highly valuable platform to implement genomic prediction of hybrid performance in oilseed rape. Furthermore, we present first insights into the application of established statistical models for prediction of important agronomical traits with contrasting patterns of polygenic control.

Integration of omics approaches to understand oil/protein content during seed development in oilseed crops

Oilseed crops, especially soybean (Glycine max) and canola/rapeseed (Brassica napus), produce seeds that are rich in both proteins and oils and that are major sources of energy and nutrition worldwide. Most of the nutritional content in the seed is accumulated in the embryo during the seed filling stages of seed development. Understanding the metabolic pathways that are active during seed filling and how they are regulated are essential prerequisites to crop improvement. In this review, we summarize various omics studies of soybean and canola/rapeseed during seed filling, with emphasis on oil and protein traits, to gain a systems-level understanding of seed development. Currently, most (80-85%) of the soybean and rapeseed reference genomes have been sequenced (950 and 850 megabases, respectively). Parallel to these efforts, extensive omics datasets from different seed filling stages have become available. Transcriptome and proteome studies have detected preponderance of starch metabolism and glycolysis enzymes to be the possible cause of higher oil in B. napus compared to other crops. Small RNAome studies performed during the seed filling stages have revealed miRNA-mediated regulation of transcription factors, with the suggestion that this interaction could be responsible for transitioning the seeds from embryogenesis to maturation. In addition, progress made in dissecting the regulation of de novo fatty acid synthesis and protein storage pathways is described. Advances in high-throughput omics and comprehensive tissue-specific analyses make this an exciting time to attempt knowledge-driven investigation of complex regulatory pathways.

Investigation of the Genetic Diversity and Quantitative Trait Loci Accounting for Important Agronomic and Seed Quality Traits in Brassica carinata

Brassica carinata (BBCC) is an allotetraploid in Brassicas with unique alleles for agronomic traits and has huge potential as source for biodiesel production. To investigate the genome-wide molecular diversity, population structure and linkage disequilibrium (LD) pattern in this species, we genotyped a panel of 81 accessions of B. carinata with genotyping by sequencing approach DArTseq, generating a total of 54,510 polymorphic markers. Two subpopulations were exhibited in the B. carinata accessions. The average distance of LD decay (r2 = 0.1) in B subgenome (0.25 Mb) was shorter than that of C subgenome (0.40 Mb). Genome-wide association analysis (GWAS) identified a total of seven markers significantly associated with five seed quality traits in two experiments. To further identify the quantitative trait loci (QTL) for important agronomic and seed quality traits, we phenotyped a doubled haploid (DH) mapping population derived from the "YW" cross between two parents (Y-BcDH64 and W-BcDH76) representing from the two subpopulations. The YW DH population and its parents were grown in three contrasting environments; spring (Hezheng and Xining, China), semi-winter (Wuhan, China), and spring (Wagga Wagga, Australia) across 5 years for QTL mapping. Genetic bases of phenotypic variation in seed yield and its seven related traits, and six seed quality traits were determined. A total of 282 consensus QTL accounting for these traits were identified including nine major QTL for flowering time, oleic acid, linolenic acid, pod number of main inflorescence, and seed weight. Of these, 109 and 134 QTL were specific to spring and semi-winter environment, respectively, while 39 consensus QTL were identified in both contrasting environments. Two QTL identified for linolenic acid (B3) and erucic acid (C7) were validated in the diverse lines used for GWAS. A total of 25 QTL accounting for flowering time, erucic acid, and oleic acid were aligned to the homologous QTL or candidate gene regions in the C genome of B. napus. These results would not only provide insights for genetic improvement of this species, but will also identify useful genetic variation hidden in the Cc subgenome of B. carinata to improve canola cultivars.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

A Major Locus for Manganese Tolerance Maps on Chromosome A09 in a Doubled Haploid Population of Brassica napus L

Soil acidity poses a major threat to productivity of several crops; mainly due to the prevalence of toxic levels of Al³⁺ and Mn²⁺. Crop productivity could be harnessed on acid soils via the development of plant varieties tolerant to phytotoxic levels of these cations. In this study, we investigated the extent of natural variation for Mn²⁺ tolerance among ten parental lines of the Australian and International canola mapping populations. Response to Mn²⁺ toxicity was measured on the bases of cotyledon chlorosis, shoot biomass, and leaf area in nutrient solution under control (9 μM of MnCl₂⋅4H₂O) and Mn treatment (125 μM of MnCl₂⋅4H₂O). Among parental lines, we selected Darmor-bzh and Yudal that showed significant and contrasting variation in Mn²⁺ tolerance to understand genetic control and identify the quantitative trait loci (QTL) underlying Mn²⁺ tolerance. We evaluated parental lines and their doubled haploid (DH) progenies (196 lines) derived from an F₁ cross, Darmor-bzh/Yudal for Mn²⁺ tolerance. Mn²⁺-tolerant genotypes had significantly higher shoot biomass and leaf area compared to Mn²⁺-sensitive genotypes. A genetic linkage map based on 7,805 DArTseq markers corresponding to 2,094 unique loci was constructed and further utilized for QTL identification. A major locus, BnMn²⁺.A09 was further mapped with a SNP marker, Bn-A09-p29012402 (LOD score of 34.6) accounting for most of the variation in Mn²⁺ tolerance on chromosome A09. This is the first report on the genomic localization of a Mn²⁺ tolerance locus in B. napus. Additionally, an ortholog of A. thaliana encoding for cation efflux facilitator transporter was located within 3,991 bp from significant SNP marker associated with BnMn²⁺.A09. A suite of genome sequence based markers (DArTseq and Illumina Infinium SNPs) flanking the BnMn²⁺.A09 locus would provide an invaluable tool for various molecular breeding applications to improve canola production and profitability on Mn²⁺ toxic soils.

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.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.

Identification of stable QTLs for seed oil content by combined linkage and association mapping in Brassica napus

Seed oil content is an important agricultural trait in rapeseed breeding. Although numerous quantitative trait locus (QTL) have been identified, most of them cannot be applied in practical breeding mainly due to environmental instability or large confidence intervals. The purpose of this study was to identify and validate high quality and more stable QTLs by combining linkage mapping and genome-wide association study (GWAS). For linkage mapping, we constructed two F₂ populations from crosses of high-oil content (∼50%) lines 6F313 and 61616 with a low-oil content (∼40%) line 51070. Two high density linkage maps spanned 1987cM (1659 bins) and 1856cM (1746 bins), respectively. For GWAS, we developed more than 34,000 high-quality SNP markers based on 227 accessions. Finally, 40 QTLs and 29 associations were established by linkage and association mapping in different environments. After merging the results, 32 consensus QTLs were obtained and 7 of them were identified by both mapping methods. Seven overlapping QTLs covered an average confidence interval of 183kb and explained the phenotypic variation of 10.23 to 24.45%. We further developed allele-specific PCR primers to identify each of the seven QTLs. These stable QTLs should be useful in gene cloning and practical breeding application.