Genome-wide association mapping and Identification of candidate genes for fatty acid composition in Brassica napus L. using SNP markers

A Golden2-like transcription factor regulates Brassica napus seed vigor after artificial aging

KEY MESSAGE

A novel GLK transcription factor, BnaA03.G2-like, regulates seed vigor after aging by affecting FAs composition in B. napus. Seeds, a most crucial materials for crop production, constitute the basis of agriculture. Seed aging usually occurs during storage, especially for seeds with high oil content. However, the genes and mechanisms underlying seed aging in Brassica napus, a widely grown oilseed crop, remain largely unexplored. In this study, nine SNPs associated with seed vigor after aging were identified through a genome-wide association study, and transgenic seeds confirmed that BnaA03.G2-like negatively regulated seed vigor after aging. Haplotype and gene-based association analyses revealed that S_A3.16900957 and S_A3.16901008 were two core variants located within BnaA03.G2-like. Transcriptome analysis and gas chromatography results suggested that BnaA03.G2-like regulated fatty acid biosynthesis, and the linoleic acid content significantly increased in the seeds of BnaA03.G2-like overexpression lines but decreased in the seeds of BnaA03.G2-like RNA interference lines. Further correlation analysis and malondialdehyde determination revealed that BnaA03.G2-like regulated seed vigor by affecting the linoleic acid content in mature seeds and the degree of membrane lipid peroxidation during seed aging. This study initially reported the function of BnaA03.G2-like and its favorable alleles, providing genetic resources for seed antiaging breeding in B. napus.

Increasing oil content in Brassica oilseed species

Brassica oilseed species are the third most important in the world, providing approximately 15 % of the total vegetable oils. Three species (Brassica rapa, B. juncea, B. napus) dominate with B. napus being the most common in Canada, China and Europe. Originally, B. napus was a crop producing seed with high erucic acid content, which still persists today, to some extent, and is used for industrial purposes. In contrast, cultivars which produce seed used for food and feed are low erucic acid cultivars which also have reduced glucosinolate content. Because of the limit to agricultural land, recent efforts have been made to increase productivity of oil crops, including Brassica oilseed species. In this article, we have detailed research in this regard. We have covered modern genetic, genomic and metabolic control analysis approaches to identifying potential targets for the manipulation of seed oil content. Details of work on the use of quantitative trait loci, genome-wide association and comparative functional genomics to highlight factors influencing seed oil accumulation are given and functional proteins which can affect this process are discussed. In summary, a wide variety of inputs are proving useful for the improvement of Brassica oilseed species, as major sources of global vegetable oil.

Identification of transcription factor BnHDG4-A08 as a novel candidate associated with the accumulation of oleic, linoleic, linolenic, and erucic acid in Brassica napus

KEY MESSAGE

We screened 47 significantly associated haplotype blocks for oleic, linoleic, linolenic, and erucic acid, with 17 blocks influencing multiple traits. A novel candidate of transcription factor BnHDG4 A08 influencing oleic, linoleic, linolenic, and erucic acid was identified, by a joint strategy of haplotype-based genome-wide association study, genomic resequencing, gene cloning, and co-expression network Fatty acid (FA) composition determines the quality and economic value of rapeseed oil (Brassica napus). However, the molecular network of FAs is unclear. In the current study, multi-strategies of haplotype-based genome-wide association study (GWAS), genomic resequencing, gene cloning, and co-expression network were joint to reveal novel genetic factors influencing FA accumulation in rapeseed. We identified 47 significantly associated haplotype blocks for oleic, linoleic, linolenic, and erucic acid, with 17 blocks influencing multiple traits, using a haplotype-based GWAS with phenotype data from 203 Chinese semi-winter accessions. A total of 61 rapeseed orthologs involved in acyl-lipid metabolism, carbohydrate metabolism, or photosynthesis were identified in these 17 blocks. Among these genes, BnHDG4-A08, encoding a class IV homeodomain leucine-zipper transcription factor, exhibited two single-nucleotide polymorphisms (SNPs) in the exon and intron, with significant associations with oleic, linoleic, linolenic, and erucic acid. Gene cloning further validated two SNPs in the exon of BnHDG4-A08 in a population with 75 accessions, leading to two amino acid changes (T372A and P366L) and significant variation of oleic, linoleic, linolenic, and erucic acid. A competitive allele-specific PCR (KASP) marker based on the SNPs was successfully developed and validated. Moreover, 98 genes exhibiting direct interconnections and high weight values with BnHDG4-A08 were identified through co-expression network analysis using transcriptome data from 13 accessions. Our study identified a novel FA candidate of transcription factor BnHDG4-A08 influencing oleic, linoleic, linolenic, and erucic acid. This gene provides a potential promising gene resource for the novel mechanistic understanding of transcription factors regulating FA accumulation.

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.

Genomic and transcriptome analyses reveal potential contributors to erucic acid biosynthesis in seeds of rapeseed (Brassica napus)

KEY MESSAGE

Through comprehensive genomic and transcriptomic analyses, we identified a set of 23 genes that act up- or downstream of erucic acid content (EAC) production in rapeseed seeds. We selected example genes to showcase the distribution of single nucleotide polymorphisms, haplotypes associated with EAC phenotypes, and the creation of molecular markers differentiating low EAC and high EAC genotypes. Erucic acid content (EAC) is a crucial trait in rapeseed, with low LEAC oil recognized for its health benefits and high EA oil holding industrial value. Despite its significance, the genomic consequences of intensive LEAC-cultivar selection and the genetic basis underlying EA regulation remain largely unexplored. To address this knowledge gap, we conducted selective signal analyses, genome-wide association studies (GWAS), and transcriptome analyses. Our investigation unveiled the genetic footprints resulting from LEAC selection in germplasm populations, drawing attention to specific loci that contribute to enriching diversity. By integrating GWAS and transcriptome analyses, we identified a set of 23 genes that play a significant role in determining EAC in seeds or are downstream consequences of EA-level alterations. These genes have emerged as promising candidates for elucidating the potential mechanisms governing EAC in rapeseed. To exemplify the findings, we selected specific genes to demonstrate the distribution of single nucleotide polymorphisms and haplotypes associated with different EAC phenotypes. Additionally, we showcased to develop molecular markers distinguishing between LEAC and high EAC genotypes.

GWAS determined genetic loci associated with callus induction in oil palm tissue culture

KEY MESSAGE

GWAS identified six loci at 25 kb downstream of WAK2, a crucial gene for cell wall and callus formation, enabling development of a SNP marker for enhanced callus induction potential. Efficient callus induction is vital for successful oil palm tissue culture, yet identifying genomic loci and markers for early detection of genotypes with high potential of callus induction remains unclear. In this study, immature male inflorescences from 198 oil palm accessions (dura, tenera and pisifera) were used as explants for tissue culture. Callus induction rates were collected at one-, two- and three-months after inoculation (C1, C2 and C3) as phenotypes. Resequencing generated 11,475,258 high quality single nucleotide polymorphisms (SNPs) as genotypes. GWAS was then performed, and correlation analysis revealed a positive association of C1 with both C2 (R = 0.81) and C3 (R = 0.50), indicating that C1 could be used as the major phenotype for callus induction rate. Therefore, only significant SNPs (P ≤ 0.05) in C1 were identified to develop markers for screening individuals with high potential of callus induction. Among 21 significant SNPs in C1, LD block analysis revealed six SNPs on chromosome 12 (Chr12) potentially linked to callus formation. Subsequently, 13 SNP markers were identified from these loci and electrophoresis results showed that marker C-12 at locus Chr12_12704856 can be used effectively to distinguish the GG allele, which showed the highest probability (69%) of callus induction. Furthermore, a rapid SNP variant detection method without electrophoresis was established via qPCR-based melting curve analysis. Our findings facilitated marker-assisted selection for specific palms with high potential of callus induction using immature male inflorescence as explant, aiding ortet palm selection in oil palm tissue culture.

Hub Genes in Stable QTLs Orchestrate the Accumulation of Cottonseed Oil in Upland Cotton via Catalyzing Key Steps of Lipid-Related Pathways

Upland cotton is the fifth-largest oil crop in the world, with an average supply of nearly 20% of vegetable oil production. Cottonseed oil is also an ideal alternative raw material to be efficiently converted into biodiesel. However, the improvement in kernel oil content (KOC) of cottonseed has not received sufficient attention from researchers for a long time, due to the fact that the main product of cotton planting is fiber. Previous studies have tagged QTLs and identified individual candidate genes that regulate KOC of cottonseed. The regulatory mechanism of oil metabolism and accumulation of cottonseed are still elusive. In the current study, two high-density genetic maps (HDGMs), which were constructed based on a recombinant inbred line (RIL) population consisting of 231 individuals, were used to identify KOC QTLs. A total of forty-three stable QTLs were detected via these two HDGM strategies. Bioinformatic analysis of all the genes harbored in the marker intervals of the stable QTLs revealed that a total of fifty-one genes were involved in the pathways related to lipid biosynthesis. Functional analysis via coexpression network and RNA-seq revealed that the hub genes in the co-expression network that also catalyze the key steps of fatty acid synthesis, lipid metabolism and oil body formation pathways (ACX4, LACS4, KCR1, and SQD1) could jointly orchestrate oil accumulation in cottonseed. This study will strengthen our understanding of oil metabolism and accumulation in cottonseed and contribute to KOC improvement in cottonseed in the future, enhancing the security and stability of worldwide food supply.

Genome-wide association analysis of tan spot disease resistance in durum wheat accessions from Tunisia

Background: Tunisia harbors a rich collection of unexploited durum wheat landraces (Triticum durum ssp. durum) that have been gradually replaced by elite cultivars since the 1970s. These landraces represent an important potential source for broadening the genetic background of elite durum wheat cultivars and for the introgression of novel genes for key traits, including disease resistance, into these cultivars. Methods: In this study, single nucleotide polymorphism (SNP) markers were used to investigate the genetic diversity and population structure of a core collection of 235 durum wheat accessions consisting mainly of landraces. The high phenotypic and genetic diversity of the fungal pathogen Pyrenophora tritici-repentis (cause of tan spot disease of wheat) in Tunisia allowed the assessment of the accessions for tan spot resistance at the adult plant stage under field conditions over three cropping seasons. A genome-wide association study (GWAS) was performed using a 90k SNP array. Results: Bayesian population structure analysis with 9191 polymorphic SNP markers classified the accessions into two groups, where groups 1 and 2 included 49.79% and 31.49% of the accessions, respectively, while the remaining 18.72% were admixtures. Principal coordinate analysis, the unweighted pair group method with arithmetic mean and the neighbor-joining method clustered the accessions into three to five groups. Analysis of molecular variance indicated that 76% of the genetic variation was among individuals and 23% was between individuals. Genome-wide association analyses identified 26 SNPs associated with tan spot resistance and explained between 8.1% to 20.2% of the phenotypic variation. The SNPs were located on chromosomes 1B (1 SNP), 2B (4 SNPs), 3A (2 SNPs), 3B (2 SNPs), 4A (2 SNPs), 4B (1 SNP), 5A (2 SNPs), 5B (4 SNPs), 6A (5 SNPs), 6B (2 SNPs), and 7B (1 SNP). Four markers, one on each of chromosomes 1B, and 5A, and two on 5B, coincided with previously reported SNPs for tan spot resistance, while the remaining SNPs were either novel markers or closely related to previously reported SNPs. Eight durum wheat accessions were identified as possible novel sources of tan spot resistance that could be introgressed into elite cultivars. Conclusion: The results highlighted the significance of chromosomes 2B, 5B, and 6A as genomic regions associated with tan spot resistance.

Factors Affecting the Quality of Canola Grains and Their Implications for Grain-Based Foods

Canola, Brassica napus L., is a major oilseed crop that has various uses in the food, feed, and industrial sectors. It is one of the most widely produced and consumed oilseeds in the world because of its high oil content and favorable fatty acid composition. Canola grains and their derived products, such as canola oil, meal, flour, and bakery products, have a high potential for food applications as they offer various nutritional and functional benefits. However, they are affected by various factors during the production cycle, post-harvest processing, and storage. These factors may compromise their quality and quantity by affecting their chemical composition, physical properties, functional characteristics, and sensory attributes. Therefore, it is important to optimize the production and processing methods of canola grains and their derived products to ensure their safety, stability, and suitability for different food applications. This literature review provides a comprehensive overview of how these factors affect the quality of canola grains and their derived products. The review also suggests future research needs and challenges for enhancing canola quality and its utilization in food.

High-Throughput Association Mapping in Brassica napus L.: Methods and Applications

Oil seed rape (Braasica napus L.) is ranked second among oil seed crops cultivated globally for edible oil for human, and seed cake for animal consumption. Recent genetic and genomics advancements highlighted the diversity that exists within B. napus, which is largely discovered using the most promising genetic markers called single nucleotide polymorphism (SNP). Their calling rate is also enhanced to ~100 folds after the continuous advancements in the next generation sequencing (NGS) technologies. As the high throughput of NGS resulted in multi-Giga bases data, the detailed quality control (QC) prior to downstream analyses is a pre-requisite. It mainly involved the removal of false positives, missing proportions, filtering of low-quality SNPs, and adjustments of minor-allele frequency and heterozygosity. After marker-trait association, for conformation of target SNPs, validations of SNPs can be performed using various methods, especially allele-specific PCR assay-based methods have been utilized for SNP genotyping of genes targeting agronomic traits and somaclonal variations occurred during transgenic studies. In the present study, the authors mainly argue on the genotypic progress, and pipelines/methods that are being used for detection, calling, filtering, and validation of SNPs. Also, insight is provided into the application of SNPs in linkage and association mapping, including QTL mapping and genome-wide association studies targeting mainly developmental traits related to the root system and plant architecture, flowering time, silique, and oil quality. Briefly, the present study provides the recent information and recommendations on the SNP genotyping methods and its applications, which can be useful for marker-assisted breeding in B. napus and other crops.

Genome-Wide Association Study of Glucosinolate Metabolites (mGWAS) in Brassica napus L

Glucosinolates (GSLs) are secondary plant metabolites that are enriched in rapeseed and related Brassica species, and they play important roles in defense due to their anti-nutritive and toxic properties. Here, we conducted a genome-wide association study of six glucosinolate metabolites (mGWAS) in rapeseed, including three aliphatic glucosinolates (m145 gluconapin, m150 glucobrassicanapin and m151 progoitrin), one aromatic glucosinolate (m157 gluconasturtiin) and two indole glucosinolates (m165 indolylmethyl glucosinolate and m172 4-hydroxyglucobrassicin), respectively. We identified 113 candidate intervals significantly associated with these six glucosinolate metabolites. In the genomic regions linked to the mGWAS peaks, 187 candidate genes involved in glucosinolate biosynthesis (e.g., BnaMAM1, BnaGGP1, BnaSUR1 and BnaMYB51) and novel genes (e.g., BnaMYB44, BnaERF025, BnaE2FC, BnaNAC102 and BnaDREB1D) were predicted based on the mGWAS, combined with analysis of differentially expressed genes. Our results provide insight into the genetic basis of glucosinolate biosynthesis in rapeseed and should facilitate marker-based breeding for improved seed quality in Brassica species.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Clone and Function Verification of the OPR gene in Brassica napus Related to Linoleic Acid Synthesis

BACKGROUND

Fatty acid composition and content affect rapeseed oil quality. Fatty acid synthesis-related genes in rapeseed have been studied globally by researchers. Nevertheless, rapeseed oil is mainly composed of seven different fatty acids (FA), and each fatty acid was regulated by different genes. Furthermore, different FA affect each other, which needs continuous and in-depth research to obtain more clear results in Brassica napus.

RESULTS

In this paper, broad-scale miRNA expression profiles were constructed and 21 differentially expressed miRNAs were detected. GO enrichment analysis showed that most up-regulated proteins were involved in transcription factor activity and catalytic activity. KEGG pathway enrichment analysis indicated that 20 pathways involving 36 target genes were enriched, of which the bna00592 pathway may be involved in fatty acid metabolism. The results were verified using a quantitative real-time PCR (RT-qPCR) analysis, we found that the target gene of bna-miR156b > c > g was the OPR (12-oxo-phytodienoic acid reductase). Four copies of OPR gene were found, and the over-expression vectors (pCAMBIA1300-35 s-OPR and pCAMBIA1300-RNAi-OPR) were constructed to verify their functions. In T1 and T2 generation, the content of linoleic acid (LA) increased significantly in OE but deceased in OPRi.

CONCLUSIONS

This is the first study to provide four copies of the OPR gene that regulates LA metabolism, can be used for the molecular mechanism of LA and optimizing fatty acid profiles in oilseed for breeding programs.

Genome-wide association studies dissect the G × E interaction for agronomic traits in a worldwide collection of safflowers (Carthamus tinctorius L.)

Genome-wide association studies were conducted using a globally diverse safflower (Carthamus tinctorius L.) Genebank collection for grain yield (YP), days to flowering (DF), plant height (PH), 500 seed weight (SW), seed oil content (OL), and crude protein content (PR) in four environments (sites) that differed in water availability. Phenotypic variation was observed for all traits. YP exhibited low overall genetic correlations (rGoverall) across sites, while SW and OL had high rGoverall and high pairwise genetic correlations (rGij) across all pairwise sites. In total, 92 marker-trait associations (MTAs) were identified using three methods, single locus genome-wide association studies (GWAS) using a mixed linear model (MLM), the Bayesian multi-locus method (BayesR), and meta-GWAS. MTAs with large effects across all sites were detected for OL, SW, and PR, and MTAs specific for the different water stress sites were identified for all traits. Five MTAs were associated with multiple traits; 4 of 5 MTAs were variously associated with the three traits of SW, OL, and PR. This study provided insights into the phenotypic variability and genetic architecture of important safflower agronomic traits under different environments.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01295-8.

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.

CRISPR/Cas9-Targeted Mutagenesis of BnaFAE1 Genes Confers Low-Erucic Acid in Brassica napus

Rapeseed (Brassica napus) is an important oilseed crop widely planted in the world, providing substantial edible oil and other nutrients for mankind. The composition of fatty acids affects the edible and processing quality of vegetable oils, among which erucic acid (EA) is potentially to cause health problems. Therefore, low erucic acid (LEA) has always been a breeding trait of B. napus. Fatty acid elongase 1 (FAE1) plays a decisive role in the synthesis of EA. There are two functional homologous copies of FAE1 on the A08 and C03 chromosomes in B. napus. In this study, we used CRISPR/Cas9 technology to create targeted mutations on these two homologous copies of BnaFAE1 in three B. napus germplasms with high EA (>30%) and high oil (>50%). Our results show that the EA content was significantly reduced by more than 10 percentage points in the mutant of BnaC03.FAE1 (c03), while the double mutation of BnaA08.FAE1 and BnaC03.FAE1 (a08c03) resulted in nearly zero EA in three BnaFAE1-edited germplasms, and the oleic acid content was increased in different degrees. In addition, knockout of BnaA08.FAE1 or/and BnaC03.FAE1 mildly decreased seed oil content, but had no significant effect on other agronomic traits. In general, we successfully created low EA germplasms of B. napus, which provides a feasible way for future low EA breeding.

Genome-Wide Association Reveals Trait Loci for Seed Glucosinolate Accumulation in Indian Mustard (Brassica juncea L.)

Glucosinolates (GSLs) are sulphur- and nitrogen-containing secondary metabolites implicated in the fitness of Brassicaceae and appreciated for their pungency and health-conferring properties. In Indian mustard (Brassica juncea L.), GSL content and composition are seed-quality-determining traits affecting its economic value. Depending on the end use, i.e., condiment or oil, different GSL levels constitute breeding targets. The genetic control of GSL accumulation in Indian mustard, however, is poorly understood, and current knowledge of GSL biosynthesis and regulation is largely based on Arabidopsis thaliana. A genome-wide association study was carried out to dissect the genetic architecture of total GSL content and the content of two major GSLs, sinigrin and gluconapin, in a diverse panel of 158 Indian mustard lines, which broadly grouped into a South Asia cluster and outside-South-Asia cluster. Using 14,125 single-nucleotide polymorphisms (SNPs) as genotyping input, seven distinct significant associations were discovered for total GSL content, eight associations for sinigrin content and 19 for gluconapin. Close homologues of known GSL structural and regulatory genes were identified as candidate genes in proximity to peak SNPs. Our results provide a comprehensive map of the genetic control of GLS biosynthesis in Indian mustard, including priority targets for further investigation and molecular marker development.

Comprehensive genomics and expression analysis of eceriferum (CER) genes in sunflower (Helianthus annuus)

Sunflower occupies the fourth position among oilseed crops the around the world. Eceriferum (CER) is an important gene family that plays critical role in very-long-chain fatty acids elongation and biosynthesis of epicuticular waxes under both biotic and abiotic stress conditions. The aim of present study was to investigate the effect of sunflower CER genes during drought stress condition. Thus, comparative analysis was undertaken for sunflower CER genes with Arabidopsis genome to determine phylogenetic relationship, chromosomal mapping, gene structures, gene ontology and conserved motifs. Furthermore, we subjected the sunflower cultivars under drought stress and used qRT-PCR analysis to explore the expression pattern of CER genes during drought conditions. We identified thirty-seven unevenly distributed CER genes in the sunflower genome. The phylogenetic analysis revealed that CER genes were grouped into seven clades in Arabidopsis, Helianthus annuus, and Gossypium hirsutum. Expression analysis showed that genes CER10 and CER60 were upregulated in sunflower during drought conditions, indicating that these genes are activated during drought stress. The results obtained will serve to characterize the CER gene family in sunflower and exploit the role of these genes in wax biosynthesis under limited water conditions.

Key message

Cuticular waxes protect the plants from drought stress, so we observed the expression of wax bio synthesis genes in recently sequences genome of Helianthus annuus. We observed that expression of wax biosynthesis genes CER10 and CER60 was upregulated when the plants were subjected to drought stress.

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.

Breeding Canola (Brassica napus L.) for Protein in Feed and Food

Interest in canola (Brassica napus L.). In response to this interest, scientists have been tasked with altering and optimizing the protein production chain to ensure canola proteins are safe for consumption and economical to produce. Specifically, the role of plant breeders in developing suitable varieties with the necessary protein profiles is crucial to this interdisciplinary endeavour. In this article, we aim to provide an overarching review of the canola protein chain from the perspective of a plant breeder, spanning from the genetic regulation of seed storage proteins in the crop to advancements of novel breeding technologies and their application in improving protein quality in canola. A review on the current uses of canola meal in animal husbandry is presented to underscore potential limitations for the consumption of canola meal in mammals. General discussions on the allergenic potential of canola proteins and the regulation of novel food products are provided to highlight some of the challenges that will be encountered on the road to commercialization and general acceptance of canola protein as a dietary protein source.

Advances and Challenges for QTL Analysis and GWAS in the Plant-Breeding of High-Yielding: A Focus on Rapeseed

Yield is one of the most important agronomic traits for the breeding of rapeseed (Brassica napus L), but its genetic dissection for the formation of high yield remains enigmatic, given the rapid population growth. In the present review, we review the discovery of major loci underlying important agronomic traits and the recent advancement in the selection of complex traits. Further, we discuss the benchmark summary of high-throughput techniques for the high-resolution genetic breeding of rapeseed. Biparental linkage analysis and association mapping have become powerful strategies to comprehend the genetic architecture of complex agronomic traits in crops. The generation of improved crop varieties, especially rapeseed, is greatly urged to enhance yield productivity. In this sense, the whole-genome sequencing of rapeseed has become achievable to clone and identify quantitative trait loci (QTLs). Moreover, the generation of high-throughput sequencing and genotyping techniques has significantly enhanced the precision of QTL mapping and genome-wide association study (GWAS) methodologies. Furthermore, this study demonstrates the first attempt to identify novel QTLs of yield-related traits, specifically focusing on ovule number per pod (ON). We also highlight the recent breakthrough concerning single-locus-GWAS (SL-GWAS) and multi-locus GWAS (ML-GWAS), which aim to enhance the potential and robust control of GWAS for improved complex traits.

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.

Genotyping and lipid profiling of 601 cultivated sunflower lines reveals novel genetic determinants of oil fatty acid content

BACKGROUND

Sunflower is an important oilseed crop domesticated in North America approximately 4000 years ago. During the last century, oil content in sunflower was under strong selection. Further improvement of oil properties achieved by modulating its fatty acid composition is one of the main directions in modern oilseed crop breeding.

RESULTS

We searched for the genetic basis of fatty acid content variation by genotyping 601 inbred sunflower lines and assessing their lipid and fatty acid composition. Our genome-wide association analysis based on the genotypes for 15,483 SNPs and the concentrations of 23 fatty acids, including minor fatty acids, revealed significant genetic associations for eleven of them. Identified genomic regions included the loci involved in rare fatty acids variation on chromosomes 3 and 14, explaining up to 34.5% of the total variation of docosanoic acid (22:0) in sunflower oil.

CONCLUSIONS

This is the first large scale implementation of high-throughput lipidomic profiling to sunflower germplasm characterization. This study contributes to the genetic characterization of Russian sunflower collections, which made a substantial contribution to the development of sunflower as the oilseed crop worldwide, and provides new insights into the genetic control of oil composition that can be implemented in future studies.

Toward the identification of molecular markers associated with phytochemical traits in the Iranian sumac (Rhus coriaria L.) population

Sumac (Rhus coriaria L.) is one of the important forest species dispersed in the northwest of Iran. It is one of the important spice in Iran and the Middle East because of active components containing organic acids, phenolic acids, flavonoids, anthocyanins, tannins and terpenoids. This study aimed to investigate population structure and linkage disequilibrium (LD) extent within Rhus coriaria L. genotypes using ISSR markers and identify molecular markers associated with phytochemical traits using association analysis. In the molecular part of the experiment, the genetic diversity of 75 sumac genotypes from five different areas of northwest Iran was assessed by 18 ISSR primers. In the phenotypic assessment, the fruits of the sumac genotypes were analyzed using HPLC-LC/MSMS for determining phytochemical components including maleic acid, ellagic acid, maleic acid hexoside, gallic acid, coumaric acid, quercetin, caftaric acid, and linoleic acid. The phenotypic data analysis revealed the great phenotypic diversity among and within Iranian sumac populations for the studied phytochemical traits. The studied sumac genotypes were divided into two subpopulations based on molecular marker-based structure analysis. A significant level of LD was observed in 11.64% of the ISSR marker pairs (p < .05). The mixed linear model procedure showed that 12 loci had a significant association with investigated traits. The ISSR loci identified in this study can be potentially used in marker-assisted selection in sumac breeding programs.

Identification of resistance loci against new pathotypes of Plasmodiophora brassicae in Brassica napus based on genome-wide association mapping

Genetic resistance is a successful strategy for management of clubroot (Plasmodiophora brassicae) of brassica crops, but resistance can break down quickly. Identification of novel sources of resistance is especially important when new pathotypes arise. In the current study, the reaction of 177 accessions of Brassica napus to four new, virulent pathotypes of P. brassicae was assessed. Each accession was genotyped using genotyping by sequencing to identify and map novel sources of clubroot resistance using mixed linear model (MLM) analysis. The majority of accessions were highly susceptible (70–100 DSI), but a few accessions exhibited strong resistance (0–20 DSI) to pathotypes 5X (21 accessions), 3A (8), 2B (7), and 3D (15), based on the Canadian Clubroot Differential system. In total, 301,753 SNPs were mapped to 19 chromosomes. Population structure analysis indicated that the 177 accessions belong to seven major populations. SNPs were associated with resistance to each pathotype using MLM. In total, 13 important SNP loci were identified, with 9 SNPs mapped to the A-genome and 4 to the C-genome. The SNPs were associated with resistance to pathotypes 5X (2 SNPs), 3A (4), 2B (5) and 3D (6). A Blast search of 1.6 Mb upstream and downstream from each SNP identified 13 disease-resistance genes or domains. The distance between a SNP locus and the nearest resistance gene ranged from 0.04 to 0.74 Mb. The resistant lines and SNP markers identified in this study can be used to breed for resistance to the most prevalent new pathotypes of P. brassicae in Canada.

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.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

The pineapple MADS-box gene family and the evolution of early monocot flower

Unlike the flower of the model monocot rice, which has diverged greatly from the ancestral monocot flower, the pineapple (Ananas comosus) flower is more typical of monocot flowers. Here, we identified 43 pineapple genes containing MADS-box domains, including 11 type I and 32 type II genes. RNA-seq expression data generated from five pineapple floral organs (sepals, petals, stamens, pistils, and ovules) and quantitative real-time PCR revealed tissue-specific expression patterns for some genes. We found that AcAGL6 and AcFUL1 were mainly expressed in sepals and petals, suggesting their involvement in the regulation of these floral organs. A pineapple ‘ABCDE’ model was proposed based on the phylogenetic analysis and expression patterns of MADS-box genes. Unlike rice and orchid with frequent species-specific gene duplication and subsequent expression divergence, the composition and expression of the ABCDE genes were conserved in pineapple. We also found that AcSEP1/3, AcAG, AcAGL11a/b/c, and AcFUL1 were highly expressed at different stages of fruit development and have similar expression profiles, implicating these genes’ role in fruit development and ripening processes. We propose that the pineapple flower can be used as a model for studying the ancestral form of monocot flowers to investigate their development and evolutionary history.

Genome-wide association study reveals a patatin-like lipase relating to the reduction of seed oil content in Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) is an important oil crop world-widely cultivated, and seed oil content (SOC) is one of the most important traits for rapeseed. To increase SOC, many efforts for promoting the function of genes on lipid biosynthesis pathway have been previously made. However, seed oil formation is a dynamic balance between lipid synthesis and breakdown. It is, therefore, also reasonable to weaken or eliminate the function of genes involved in lipid degradation for a higher final SOC.

RESULTS

We applied a genome-wide association study (GWAS) on SOC in a collection of 290 core germplasm accessions. A total of 2,705,480 high-quality SNPs were used in the GWAS, and we identified BnaC07g30920D, a patatin-like lipase (PTL) gene, that was associated with SOC. In particular, six single-nucleotide-polymorphisms (SNPs) in the promoter region of BnaC07g30920D were associated with the significant reduction of SOC, leading to a 4.7-6.2% reduction of SOCs. We performed in silico analysis to show a total of 40 PTLs, which were divided into four clades, evenly distributed on the A and C subgenomes of Brassica napus. RNA-seq analysis unveiled that BnPTLs were preferentially expressed in reproductive tissues especially maturing seeds.

CONCLUSIONS

We identified BnaC07g30920D, a BnPTL gene, that was associated with SOC using GWAS and performed in silico analysis of 40 PTLs in Brassica napus. The results enrich our knowledge about the SOC formation in rapeseed and facilitate the future study in functional characterization of BnPTL genes.

New Transcriptome-Based SNP Markers for Noug (Guizotia abyssinica) and Their Conversion to KASP Markers for Population Genetics Analyses

The development and use of genomic resources are essential for understanding the population genetics of crops for their efficient conservation and enhancement. Noug (Guizotia abyssinica) is an economically important oilseed crop in Ethiopia and India. The present study sought to develop new DNA markers for this crop. Transcriptome sequencing was conducted on two genotypes and 628 transcript sequences containing 959 single nucleotide polymorphisms (SNPs) were developed. A competitive allele-specific PCR (KASP) assay was developed for the SNPs and used for genotyping of 24 accessions. A total of 554 loci were successfully genotyped across the accessions, and 202 polymorphic loci were used for population genetics analyses. Polymorphism information content (PIC) of the loci varied from 0.01 to 0.37 with a mean of 0.24, and about 49% of the loci showed significant deviation from the Hardy-Weinberg equilibrium. The mean expected heterozygosity was 0.27 suggesting moderately high genetic variation within accessions. Low but significant differentiation existed among accessions (FST = 0.045, p < 0.0001). Landrace populations from isolated areas may have useful mutations and should be conserved and used in breeding this crop. The genomic resources developed in this study were shown to be useful for population genetics research and can also be used in, e.g., association genetics.

The Use of Genetic and Gene Technologies in Shaping Modern Rapeseed Cultivars (Brassica napus L.)

Since their domestication, Brassica oilseed species have undergone progressive transformation allied with the development of breeding and molecular technologies. The canola (Brassica napus) crop has rapidly expanded globally in the last 30 years with intensive innovations in canola varieties, providing for a wider range of markets apart from the food industry. The breeding efforts of B. napus, the main source of canola oil and canola meal, have been mainly focused on improving seed yield, oil quality, and meal quality along with disease resistance, abiotic stress tolerance, and herbicide resistance. The revolution in genetics and gene technologies, including genetic mapping, molecular markers, genomic tools, and gene technology, especially gene editing tools, has allowed an understanding of the complex genetic makeup and gene functions in the major bioprocesses of the Brassicales, especially Brassica oil crops. Here, we provide an overview on the contributions of these technologies in improving the major traits of B. napus and discuss their potential use to accomplish new improvement targets.

Breeding and biotechnological interventions for trait improvement: status and prospects

Present review describes the molecular tools and strategies deployed in the trait discovery and improvement of major crops. The prospects and challenges associated with these approaches are discussed. Crop improvement relies on modulating the genes and genomic regions underlying key traits, either directly or indirectly. Direct approaches include overexpression, RNA interference, genome editing, etc., while breeding majorly constitutes the indirect approach. With the advent of latest tools and technologies, these strategies could hasten the improvement of crop species. Next-generation sequencing, high-throughput genotyping, precision editing, use of space technology for accelerated growth, etc. had provided a new dimension to crop improvement programmes that work towards delivering better varieties to cope up with the challenges. Also, studies have widened from understanding the response of plants to single stress to combined stress, which provides insights into the molecular mechanisms regulating tolerance to more than one stress at a given point of time. Altogether, next-generation genetics and genomics had made tremendous progress in delivering improved varieties; however, the scope still exists to expand its horizon to other species that remain underutilized. In this context, the present review systematically analyses the different genomics approaches that are deployed for trait discovery and improvement in major species that could serve as a roadmap for executing similar strategies in other crop species. The application, pros, and cons, and scope for improvement of each approach have been discussed with examples, and altogether, the review provides comprehensive coverage on the advances in genomics to meet the ever-growing demands for agricultural produce.

Deciphering the transcriptional regulatory networks that control size, color, and oil content in Brassica rapa seeds

BACKGROUND

Brassica rapa is an important oilseed and vegetable crop species and is the A subgenome donor of two important oilseed Brassica crops, Brassica napus and Brassica juncea. Although seed size (SZ), seed color (SC), and oil content (OC) substantially affect seed yield and quality, the mechanisms regulating these traits in Brassica crops remain unclear.

RESULTS

We collected seeds from a pair of B. rapa accessions with significantly different SZ, SC, and OC at seven seed developmental stages (every 7 days from 7 to 49 days after pollination), and identified 28,954 differentially expressed genes (DEGs) from seven pairwise comparisons between accessions at each developmental stage. K-means clustering identified a group of cell cycle-related genes closely connected to variation in SZ of B. rapa. A weighted correlation analysis using the WGCNA package in R revealed two important co-expression modules comprising genes whose expression was positively correlated with SZ increase and negatively correlated with seed yellowness, respectively. Upregulated expression of cell cycle-related genes in one module was important for the G₂/M cell cycle transition, and the transcription factor Bra.A05TSO1 seemed to positively stimulate the expression of two CYCB1;2 genes to promote seed development. In the second module, a conserved complex regulated by the transcription factor TT8 appear to determine SC through downregulation of TT8 and its target genes TT3, TT18, and ANR. In the third module, WRI1 and FUS3 were conserved to increase the seed OC, and Bra.A03GRF5 was revealed as a key transcription factor on lipid biosynthesis. Further, upregulation of genes involved in triacylglycerol biosynthesis and storage in the seed oil body may increase OC. We further validated the accuracy of the transcriptome data by quantitative real-time PCR of 15 DEGs. Finally, we used our results to construct detailed models to clarify the regulatory mechanisms underlying variations in SZ, SC, and OC in B. rapa.

CONCLUSIONS

This study provides insight into the regulatory mechanisms underlying the variations of SZ, SC, and OC in plants based on transcriptome comparison. The findings hold great promise for improving seed yield, quality and OC through genetic engineering of critical genes in future molecular breeding.

High-Density Genetic Linkage Mapping of Lepidium Based on Genotyping-by-Sequencing SNPs and Segregating Contig Tag Haplotypes

Lepidium campestre has been targeted for domestication as future oilseed and catch crop. Three hundred eighty plants comprising genotypes of L. campestre, Lepidium heterophyllum, and their interspecific F2 mapping population were genotyped using genotyping by sequencing (GBS), and the generated polymorphic markers were used for the construction of high-density genetic linkage map. TASSEL-GBS, a reference genome-based pipeline, was used for this analysis using a draft L. campestre whole genome sequence. The analysis resulted in 120,438 biallelic single-nucleotide polymorphisms (SNPs) with minor allele frequency (MAF) above 0.01. The construction of genetic linkage map was conducted using MSTMap based on phased SNPs segregating in 1:2:1 ratio for the F2 individuals, followed by genetic mapping of segregating contig tag haplotypes as dominant markers against the linkage map. The final linkage map consisted of eight linkage groups (LGs) containing 2,330 SNP markers and spanned 881 Kosambi cM. Contigs (10,302) were genetically mapped to the eight LGs, which were assembled into pseudomolecules that covered a total of ∼120.6 Mbp. The final size of the pseudomolecules ranged from 9.4 Mbp (LG-4) to 20.4 Mpb (LG-7). The following major correspondence between the eight Lepidium LGs (LG-1 to LG-8) and the five Arabidopsis thaliana (At) chromosomes (Atx-1-Atx-5) was revealed through comparative genomics analysis: LG-1&2_Atx-1, LG-3_Atx-2&3, LG-4_Atx-2, LG-5_Atx-2&Atx-3, LG-6_Atx-4&5, LG-7_Atx-4, and LG-8_Atx-5. This analysis revealed that at least 66% of the sequences of the LGs showed high collinearity with At chromosomes. The sequence identity between the corresponding regions of the LGs and At chromosomes ranged from 80.6% (LG-6) to 86.4% (LG-8) with overall mean of 82.9%. The map positions on Lepidium LGs of the homologs of 24 genes that regulate various traits in A. thaliana were also identified. The eight LGs revealed in this study confirm the previously reported (1) haploid chromosome number of eight in L. campestre and L. heterophyllum and (2) chromosomal fusion, translocation, and inversion events during the evolution of n = 8 karyotype in ancestral species shared by Lepidium and Arabidopsis to n = 5 karyotype in A. thaliana. This study generated highly useful genomic tools and resources for Lepidium that can be used to accelerate its domestication.

GWAS and co-expression network combination uncovers multigenes with close linkage effects on the oleic acid content accumulation in Brassica napus

Background

Strong artificial and natural selection causes the formation of highly conserved haplotypes that harbor agronomically important genes. GWAS combination with haplotype analysis has evolved as an effective method to dissect the genetic architecture of complex traits in crop species.

Results

We used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with oleic acid (C18:1) in rapeseed. Six haplotype regions were identified as significantly associated with oleic acid (C18:1) that mapped to chromosomes A02, A07, A08, C01, C02, and C03. Additionally, whole-genome sequencing of 50 rapeseed accessions revealed three genes (BnmtACP2-A02, BnABCI13-A02 and BnECI1-A02) in the A02 chromosome haplotype region and two genes (BnFAD8-C02 and BnSDP1-C02) in the C02 chromosome haplotype region that were closely linked to oleic acid content phenotypic variation. Moreover, the co-expression network analysis uncovered candidate genes from these two different haplotype regions with potential regulatory interrelationships with oleic acid content accumulation.

Conclusions

Our results suggest that several candidate genes are closely linked, which provides us with an opportunity to develop functional haplotype markers for the improvement of the oleic acid content in rapeseed.

ddRAD sequencing-based identification of inter-genepool SNPs and association analysis in Brassica juncea

BACKGROUND

Narrow genetic base, complex allo-tetraploid genome and presence of repetitive elements have led the discovery of single nucleotide polymorphisms (SNPs) in Brassica juncea (AABB; 2n = 4x = 36) at a slower pace. Double digest RAD (ddRAD) - a genome complexity reduction technique followed by NGS was used to generate a total of 23 million paired-end reads from three genotypes each of Indian (Pusa Tarak, RSPR-01 and Urvashi) and Exotic (Donskaja IV, Zem 1 and EC287711) genepools.

RESULTS

Sequence data analysis led to the identification of 10,399 SNPs in six genotypes at a read depth of 10x coverage among the genotypes of two genepools. A total of 44 hyper-variable regions (nucleotide variation hotspots) were also found in the genome, of which 93% were found to be a part of coding genes/regions. The functionality of the identified SNPs was estimated by genotyping a subset of SNPs on MassARRAY® platform among a diverse set of B. juncea genotypes. SNP genotyping-based genetic diversity and population studies placed the genotypes into two distinct clusters based mostly on the place of origin. The genotypes were also characterized for six morphological traits, analysis of which revealed a significant difference in the mean values between Indian and Exotic genepools for six traits. The association analysis for six traits identified a total of 45 significant marker-trait associations on 11 chromosomes of A- and B- group of progenitor genomes.

CONCLUSIONS

Despite narrow diversity, the ddRAD sequencing was able to identify large number of nucleotide polymorphisms between the two genepools. Association analysis led to the identification of common SNPs/genomic regions associated between flowering and maturity traits, thereby underscoring the possible role of common chromosomal regions-harboring genes controlling flowering and maturity in Brassica juncea.

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 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.

Comparative transcriptomics analysis uncovers alternative splicing events and molecular markers in cabbage (Brassica oleracea L.)

Alternative splicing (AS) events were identified and verified in cabbage by comparative transcriptome analysis. The corresponding markers were developed and the germplasm resources were identified. Alternative splicing (AS) is a central regulatory mechanism that greatly contributes to plant gene expression and transcriptome diversity. A large body of evidence has shown that AS complexity is relevant for plant development, evolution, complexity, and adaptation. Both insertion/deletion (InDel) and single nucleotide polymorphism (SNP) are typically co-dominant inheritance markers and have abundant polymorphisms. These have been widely used for marker-assisted selection, genetic mapping, and germplasm identification in plants. However, little is known about the molecular mechanisms underlying AS events and the development of markers including SNP and InDel from the cabbage transcriptome. In this study, three cabbage transcriptome datasets were collected and aligned to the cabbage reference genome to analyze AS events and marker development. 31,524 AS events were identified from three cabbage genotypes, accounting for 20.8% of the total cabbage genes. Alternative 3' splice site donor (A3SS) was the most frequent type of the four main AS events in cabbage. 70,475 InDels and 706,269 SNPs were identified with average frequencies of 1 InDel/6.9 kb and 1 SNP/0.7 kb, respectively. 71,942 potential SSRs were identified in 53,129 assembled unigenes with a density of 1 SSR/6.8 kb. The ratio of SNPs with synonymous/non-synonymous mutations was 1:0.65. 142 InDels and 36 SNPs were randomly selected and validated via Sanger sequencing and polymorphism was found among 66.2% of the InDels and 78.6% of the SNPs. Furthermore, 35 informative InDel markers were successfully used for genetic diversity analysis on 36 cabbage accessions. These results facilitate understanding of the molecular regulation mechanism underlying AS events in cabbage. They also provide molecular marker resource data for genetic mapping construction and germplasm identification, and facilitate the genetic improvement of cabbage via breeding.

CropSNPdb: a database of SNP array data for Brassica crops and hexaploid bread wheat

Advances in sequencing technology have led to a rapid rise in the genomic data available for plants, driving new insights into the evolution, domestication and improvement of crops. Single nucleotide polymorphisms (SNPs) are a major component of crop genomic diversity, and are invaluable as genetic markers in research and breeding programs. High-throughput SNP arrays, or 'SNP chips', can generate reproducible sets of informative SNP markers and have been broadly adopted. Although there are many public repositories for sequencing data, which are routinely uploaded, there are no formal repositories for crop SNP array data. To make SNP array data more easily accessible, we have developed CropSNPdb (http://snpdb.appliedbioinformatics.com.au), a database for SNP array data produced by the Illumina Infinium™ hexaploid bread wheat (Triticum aestivum) 90K and Brassica 60K arrays. We currently host SNPs from datasets covering 526 Brassica lines and 309 bread wheat lines, and provide search, download and upload utilities for users. CropSNPdb provides a useful repository for these data, which can be applied for a range of genomics and molecular crop-breeding activities.

Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement

Brassica napus (2n = 4x = 38, AACC) is an important allopolyploid crop derived from interspecific crosses between Brassica rapa (2n = 2x = 20, AA) and Brassica oleracea (2n = 2x = 18, CC). However, no truly wild B. napus populations are known; its origin and improvement processes remain unclear. Here, we resequence 588 B. napus accessions. We uncover that the A subgenome may evolve from the ancestor of European turnip and the C subgenome may evolve from the common ancestor of kohlrabi, cauliflower, broccoli, and Chinese kale. Additionally, winter oilseed may be the original form of B. napus. Subgenome-specific selection of defense-response genes has contributed to environmental adaptation after formation of the species, whereas asymmetrical subgenomic selection has led to ecotype change. By integrating genome-wide association studies, selection signals, and transcriptome analyses, we identify genes associated with improved stress tolerance, oil content, seed quality, and ecotype improvement. They are candidates for further functional characterization and genetic improvement of B. napus.

Brassica napus is a globally important oil crop, but the origin of the allotetraploid genome and its improvement process are largely unknown. Here, the authors take a population genetic approach to resolve its origin and evolutionary history, and identify candidate genes related to important agricultural traits.

Evolutionary aspects of direct or indirect selection for seed size and seed metabolites in Brassica juncea and diploid progenitor species

Seed size and seed metabolites have been the targets of direct or indirect selection during domestication and subsequent crop breeding. Understanding these traits and associated genetics can prove very useful for plant translational research. Large germplasm assemblage (235) of Brassica juncea and its progenitor species (B. rapa and B. nigra) was evaluated to establish seed trait variations for seed size and seed metabolites. Seeds were smallest in B. nigra and largest in B. juncea. Australian B. juncea and Indian B. rapa var brown sarson types averaged more seed oil content. Seed size and oil content were generally higher in modern cultivars in comparison to the land races. Allelic diversity for known associated genes for seed-size and oil-content (AP2, ARF2, TTG2, GRF2, GL2, CYP78A5, CYP78A6, MINI3, IKU2, IKU1, BRI1, DGAT, GPDH, LPAAT, GPAT and DA1) was studied so as to infer the effect of domestication on seed traits. Three genes (IKU1, IKU2, AP2) in B. rapa, two (TTG2 and GL2) in B. nigra and two (IKU1 and GRF2) in natural B. juncea were identified as targets of selection on the basis of Fst outlier and/or sequence diversity tests. We report parallel divergence for seed traits between B. juncea and B. rapa. Directional selection appeared stronger for seed-size as compared to correlated seed metabolites. Positive selection on seed-size is likely to have played a significant role in structuring regional variation in the germplasm.

Identification and Functional Analysis of Two New Mutant BnFAD2 Alleles That Confer Elevated Oleic Acid Content in Rapeseed

Rapeseed (Brassica napus L.) is a vital oil crop worldwide. High oleic acid content is a desirable quality trait for rapeseed oil, which makes it more beneficial to human health. However, many germplasm resources with high oleic acid content in rapeseed have not been evaluated with regard to their genotypes, making it difficult to select the best strains with this trait for the breeding of high oleic acid rapeseed variety. This work was to explore the gene-regulation mechanism of this trait using a new super-high oleic acid content (∼85%) line N1379T as genetic material. In this study, the sequences of four homologous fatty acid desaturase (BnFAD2) genes were compared between super-high (∼85%, N1379T) and normal (∼63%) oleic acid content lines. Results showed that there were two single-nucleotide polymorphisms (SNPs) in BnFAD2-1 and BnFAD2-2, respectively, which led to the amino acid changes (E106K and G303E) in the corresponding proteins. Functional analysis of both genes in yeast confirmed that these SNPs were loss-of-function mutations, thus limiting the conversion of oleic acid to linoleic acid and resulting in the considerable accumulation of oleic acid. Moreover, two specific cleaved amplified polymorphic sequences (CAPS) markers for the two SNPs were developed to identify genotypes of each line in the F₂ and BC₁ populations. Furthermore, these two mutant loci of BnFAD2-1 and BnFAD2-2 genes were positively associated with elevated oleic acid levels and had a similar effect with regard to the increase of oleic acid content. Taken together, these two novel SNPs in two different BnFAD2 genes jointly regulated the high oleic acid trait in this special germplasm. The study provided insight into the genetic regulation involved in oleic acid accumulation and highlighted the use of new alleles of BnFAD2-1 and BnFAD2-2 in breeding high oleic acid rapeseed varieties.

Dissection of the genetic architecture of three seed-quality traits and consequences for breeding in Brassica napus

Genome-wide association studies (GWASs) combining high-throughput genome resequencing and phenotyping can accelerate the dissection of genetic architecture and identification of genes for plant complex traits. In this study, we developed a rapeseed genomic variation map consisting of 4 542 011 SNPs and 628 666 INDELs. GWAS was performed for three seed-quality traits, including erucic acid content (EAC), glucosinolate content (GSC) and seed oil content (SOC) using 3.82 million polymorphisms in an association panel. Six, 49 and 17 loci were detected to be associated with EAC, GSC and SOC in multiple environments, respectively. The mean total contribution of these loci in each environment was 94.1% for EAC and 87.9% for GSC, notably higher than that for SOC (40.1%). A high correlation was observed between phenotypic variance and number of favourable alleles for associated loci, which will contribute to breeding improvement by pyramiding these loci. Furthermore, candidate genes were detected underlying associated loci, based on functional polymorphisms in gene regions where sequence variation was found to correlate with phenotypic variation. Our approach was validated by detection of well-characterized FAE1 genes at each of two major loci for EAC on chromosomes A8 and C3, along with MYB28 genes at each of three major loci for GSC on chromosomes A9, C2 and C9. Four novel candidate genes were detected by correlation between GSC and SOC and observed sequence variation, respectively. This study provides insights into the genetic architecture of three seed-quality traits, which would be useful for genetic improvement of B. napus.

Unconditional and conditional QTL analyses of seed fatty acid composition in Brassica napus L

BACKGROUND

The fatty acid composition of B. napus' seeds determines the oil's nutritional and industrial values, and affects seed germination. Many studies have reported correlations among C16:0, C18:0, C18:1, C18:2 and C18:3 based on phenotypic data; however, the genetic basis of the fatty acid composition in B. napus is still not well understood.

RESULTS

In this study, unconditional and conditional quantitative trail locus (QTL) mapping analyses were conducted using a recombinant inbred line in six environments. In total, 21 consensus QTLs each for C16:0, C18:0 and C18:2, 16 for C18:1 and 22 for C18:3 were detected by unconditional mapping. The QTLs with overlapping confidence intervals were integrated into 71 pleiotropically unique QTLs by meta-analysis. Two major QTLs, uuqA5-6 and uuqA5-7, simultaneously affected the fatty acids, except C18:0, in most of environments, with the homologous genes fatty acid desaturase 2 (FAD2) and glycerol-3-phosphate sn-2-acyltransferase 5 (GPAT5) occurring in the confidence interval of uuqA5-6, while phosphatidic acid phosphohydrolase 1 (PAH1) was assigned to uuqA5-7. Moreover, 49, 30, 48, 60 and 45 consensus QTLs were detected for C16:0, C18:0, C18:1, C18:2 and C18:3, respectively, by the conditional mapping analysis. In total, 128 unique QTLs were subsequently integrated from the 232 conditional consensus QTLs. A comparative analysis revealed that 63 unique QTLs could be identified by both mapping methodologies, and 65 additional unique QTLs were only identified in conditional mapping.

CONCLUSIONS

Thus, conditional QTL mapping for fatty acids may uncover numerous additional QTLs that were inhibited by the effects of other traits. These findings provide useful information for better understanding the genetic relationships among fatty acids at the QTL level.

Genome-wide association study reveals candidate genes influencing lipids and diterpenes contents in Coffea arabica L

Lipids, including the diterpenes cafestol and kahweol, are key compounds that contribute to the quality of coffee beverages. We determined total lipid content and cafestol and kahweol concentrations in green beans and genotyped 107 Coffea arabica accessions, including wild genotypes from the historical FAO collection from Ethiopia. A genome-wide association study was performed to identify genomic regions associated with lipid, cafestol and kahweol contents and cafestol/kahweol ratio. Using the diploid Coffea canephora genome as a reference, we identified 6,696 SNPs. Population structure analyses suggested the presence of two to three groups (K = 2 and K = 3) corresponding to the east and west sides of the Great Rift Valley and an additional group formed by wild accessions collected in western forests. We identified 5 SNPs associated with lipid content, 4 with cafestol, 3 with kahweol and 9 with cafestol/kahweol ratio. Most of these SNPs are located inside or near candidate genes related to metabolic pathways of these chemical compounds in coffee beans. In addition, three trait-associated SNPs showed evidence of directional selection among cultivated and wild coffee accessions. Our results also confirm a great allelic richness in wild accessions from Ethiopia, especially in accessions originating from forests in the west side of the Great Rift Valley.

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.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.