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

Genome-wide association study reveals genetic loci for seed density per silique in rapeseed (Brassica napus L.)

KEY MESSAGE

Two stable QTLs controlling seed density per silique were detected on chromosomes A09 and C05 in rapeseed via GWAS, and ARF18 was the only causal gene of QTL qSDPS-A09. Seed density per silique (SDPS) is a key agronomic trait that directly or indirectly affects seed yield in rapeseed (Brassica napus L.). Exploring the genetic control of SDPS is beneficial for increasing rapeseed production. In this study, we evaluated the SDPS phenotypes of 413 rapeseed cultivars (lines) across five natural environments and genotyped them by resequencing. A GWAS analysis was performed using 5,277,554 high-quality variants with the MLM_PCA + K and FarmCPU models. A total of 51 loci were identified to be significantly (p < - log10(1.88 × 10-6)) associated with SDPS, of which 5 were detected in all environments (except for SNP-2095656) by both GWAS models. Among the five loci, three were located on chromosome A09, whereas the other two loci were located on chromosome C05. The three loci on chromosome A09 and the two loci on chromosome C05 were physically close to each other. Therefore, only the two common candidate QTLs were integrated and named QTL qSDPS-A09 (320 kb) and qSDPS-C05 (331.48 kb), respectively. Sixty-seven and forty-eight candidate genes were initially identified on A09 and C05 and then narrowed down to 17 and 13 candidate genes, respectively, via LD block analyses. Gene-based association studies, haplotype analyses and expression analyses confirmed that three homologs of Arabidopsis auxin-response factor 18 (BnaA09G0559300ZS) was the most likely candidate genes underlying the QTL qSDPS-A09. ARF18Hap4 was identified as a favorable haplotype for high SDPS. These findings will aid in elucidating the genetic and molecular mechanisms of SDPS and promoting genetic modifications in rapeseed breeding.

Effect of Photoperiod on Ascorbic Acid Metabolism Regulation and Accumulation in Rapeseed (Brassica napus L.) Seedlings

Ascorbic acid (AsA) is an important antioxidant for human health. The concept of "oil-vegetable-duel-purpose" can significantly enhance the economic benefits of the rapeseed industry. Rapeseed, when utilized as a vegetable, serves as a valuable food source of AsA. In this study, we integrated transcriptome and metabolome analyses, along with substrate feeding, to identify the L-galactose pathway as the primary source for AsA production, which is primarily regulated by light. Through seven different photoperiod treatments from 12 h/12 h (light/dark) to 24 h/0 h, we found that AsA content increased with longer photoperiods, as well as chlorophyll, carotenoids, and soluble sugars. However, an excessively long photoperiod led to photooxidative stress, which negatively affected biomass accumulation in rapeseed seedlings and subsequently impacted the total accumulation of AsA. Furthermore, different enzymes respond differently to different photoperiods. Analysis of the correlation between the expression levels of AsA biosynthesis-related genes and AsA content highlighted a dynamic balancing mechanism of AsA metabolism in response to different photoperiods. The study revealed that the 16 h/8 h photoperiod is optimal for long-term AsA accumulation in rapeseed seedlings. However, extending the photoperiod before harvest can enhance AsA content without compromising yield. These findings offer novel insights into an effective strategy for the biofortification of AsA in rapeseed.

Structural variation reshapes population gene expression and trait variation in 2,105 Brassica napus accessions

Although individual genomic structural variants (SVs) are known to influence gene expression and trait variation, the extent and scale of SV impact across a species remain unknown. In the present study, we constructed a reference library of 334,461 SVs from genome assemblies of 16 representative morphotypes of neopolyploid Brassica napus accessions and detected 258,865 SVs in 2,105 resequenced genomes. Coupling with 5 tissue population transcriptomes, we uncovered 285,976 SV-expression quantitative trait loci (eQTLs) that associate with altered expression of 73,580 genes. We developed a pipeline for the high-throughput joint analyses of SV-genome-wide association studies (SV-GWASs) and transcriptome-wide association studies of phenomic data, eQTLs and eQTL-GWAS colocalization, and identified 726 SV-gene expression-trait variation associations, some of which were verified by transgenics. The pervasive SV impact on how SV reshapes trait variation was demonstrated with the glucosinolate biosynthesis and transport pathway. The study highlighting the impact of genome-wide and species-scale SVs provides a powerful methodological strategy and valuable resources for studying evolution, gene discovery and breeding.

Allelic Variation of BnaFTA2 and BnaFTC6 Is Associated With Flowering Time and Seasonal Crop Type in Rapeseed (Brassica napus L.)

Different ecological types of rapeseed (Brassica napus L.), including winter, spring, and semi-winter cultivars, exhibit varying flowering times and cannot be planted in the same cultivation areas. FLOWERING LOCUS T (FT) plays a key role in regulating flowering. In allotetraploid B. napus six copies of FT (BnaFT) have been reported. However, there is uncertainty about how the translated products of each paralog, as well as cis-allelic variations at each locus, contribute functionally to flowering time and define specific crop types. In this study, we confirm that BnaFT exhibit distinct expression patterns in different crop types of rapeseed. Using the CRISPR/Cas9 gene editing system, we provide functional evidence that the mutants between Bnaft paralogues affects the regulation of flowering time. Furthermore, we identify a new haplotype of BnaFT.A2 that is associated with early flowering time, although this appears necessary but not sufficient to confer a spring type phenotype. Three haplotypes of BnaFT.C6 were further identified and associated with both flowering time and crop types. We speculate that variations in both BnaFT.A2 and BnaFT.C6 may have undergone diversifying selection during the divergence of seasonal crop types in rapeseed.

A complex receptor locus confers responsiveness to necrosis and ethylene-inducing like peptides in Brassica napus

Brassica crops are susceptible to diseases which can be mitigated by breeding for resistance. MAMPs (microbe-associated molecular patterns) are conserved molecules of pathogens that elicit host defences known as pattern-triggered immunity (PTI). Necrosis and Ethylene-inducing peptide 1-like proteins (NLPs) are MAMPs found in a wide range of phytopathogens. We studied the response to BcNEP2, a representative NLP from Botrytis cinerea, and showed that it contributes to disease resistance in Brassica napus. To map regions conferring NLP response, we used the production of reactive oxygen species (ROS) induced during PTI across a population of diverse B. napus accessions for associative transcriptomics (AT), and bulk segregant analysis (BSA) on DNA pools created from a cross of NLP-responsive and non-responsive lines. In silico mapping with AT identified two peaks for NLP responsiveness on chromosomes A04 and C05 whereas the BSA identified one peak on A04. BSA delimited the region for NLP-responsiveness to 3 Mbp, containing ~245 genes on the Darmor-bzh reference genome and four co-segregating KASP markers were identified. The same pipeline with the ZS11 genome confirmed the highest-associated region on chromosome A04. Comparative BLAST analysis revealed unannotated clusters of receptor-like protein (RLP) homologues on ZS11 chromosome A04. However, no specific RLP homologue conferring NLP response could be identified. Our results also suggest that BR-SIGNALLING KINASE1 may be involved with modulating the NLP response. Overall, we demonstrate that responsiveness to NLP contributes to disease resistance in B. napus and define the associated genomic location. These results can have practical application in crop improvement.

The story of a decade: Genomics, functional genomics, and molecular breeding in Brassica napus

Rapeseed (Brassica napus L.) is one of the major global sources of edible vegetable oil and is also used as a feed and pioneer crop and for sightseeing and industrial purposes. Improvements in genome sequencing and molecular marker technology have fueled a boom in functional genomic studies of major agronomic characters such as yield, quality, flowering time, and stress resistance. Moreover, introgression and pyramiding of key functional genes have greatly accelerated the genetic improvement of important traits. Here we summarize recent progress in rapeseed genomics and genetics, and we discuss effective molecular breeding strategies by exploring these findings in rapeseed. These insights will extend our understanding of the molecular mechanisms and regulatory networks underlying agronomic traits and facilitate the breeding process, ultimately contributing to more sustainable agriculture throughout the world.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Comprehensive analysis of annexin gene family and its expression in response to branching architecture and salt stress in crape myrtle

BACKGROUND

Annexin (ANN) is calcium (Ca2+)-dependent and phospholipid binding protein family, which is involved in plant growth and development and response to various stresses. However, little known about ANN genes were identified from crape myrtle, an ornamental horticultural plant widely cultivated in the world.

RESULTS

Here, 9 LiANN genes were identified from Lagerstroemia indica, and their characterizations and functions were investigated in L. indica for the first time. The LiANN genes were divided into 2 subfamilies. The gene structure, chromosomal location, and collinearity relationship were also explored. In addition, the GO annotation analysis of these LiANNs indicated that they are enriched in molecular functions, cellular components, and biological processes. Moreover, transcription factors (TFs) prediction analysis revealed that bHLH, MYB, NAC, and other TFs can interact with the LiANN promoters. Interestingly, the LiANN2/4/6-9 were demonstrated to play critical roles in the branching architecture of crape myrtle. Furthermore, the LiANN2/6/8/9 were differentially expressed under salt treatment, and a series of TFs regulating LiANN2/6/8/9 expression were predicted to play essential roles in salt resistance.

CONCLUSIONS

These results shed light on profile and function of the LiANN gene family, and lay a foundation for further studies of the LiANN genes.

Lipid Metabolism and Improvement in Oilseed Crops: Recent Advances in Multi-Omics Studies

Oilseed crops are rich in plant lipids that not only provide essential fatty acids for the human diet but also play important roles as major sources of biofuels and indispensable raw materials for the chemical industry. The regulation of lipid metabolism genes is a major factor affecting oil production. In this review, we systematically summarize the metabolic pathways related to lipid production and storage in plants and highlight key research advances in characterizing the genes and regulatory factors influencing lipid anabolic metabolism. In addition, we integrate the latest results from multi-omics studies on lipid metabolism to provide a reference to better understand the molecular mechanisms underlying oil anabolism in oilseed crops.

The high-quality sequencing of the Brassica rapa 'XiangQingCai' genome and exploration of genome evolution and genes related to volatile aroma

'Vanilla' (XQC, brassica variety chinensis) is an important vegetable crop in the Brassica family, named for its strong volatile fragrance. In this study, we report the high-quality chromosome-level genome sequence of XQC. The assembled genome length was determined as 466.11 Mb, with an N50 scaffold of 46.20 Mb. A total of 59.50% repetitive sequences were detected in the XQC genome, including 47 570 genes. Among all examined Brassicaceae species, XQC had the closest relationship with B. rapa QGC ('QingGengCai') and B. rapa Pakchoi. Two whole-genome duplication (WGD) events and one recent whole-genome triplication (WGT) event occurred in the XQC genome in addition to an ancient WGT event. The recent WGT was observed to occur during 21.59-24.40 Mya (after evolution rate corrections). Our findings indicate that XQC experienced gene losses and chromosome rearrangements during the genome evolution of XQC. The results of the integrated genomic and transcriptomic analyses revealed critical genes involved in the terpenoid biosynthesis pathway and terpene synthase (TPS) family genes. In summary, we determined a chromosome-level genome of B. rapa XQC and identified the key candidate genes involved in volatile fragrance synthesis. This work can act as a basis for the comparative and functional genomic analysis and molecular breeding of B. rapa in the future.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

Comparative Transcriptome Analysis Points to the Biological Processes of Hybrid Incompatibility between Brassica napus and B. oleracea

Improving Brassica napus via introgression of the genome components from its parental species, B. oleracea and B. rapa, is an important breeding strategy. Interspecific hybridization between B. napus and B. rapa is compatible with high rate of survival ovules, while the hybridization between B. napus and B. oleracea is incompatible with the high occurrence of embryo abortion. To understand the diverse embryo fate in the two interspecific hybridizations, here, the siliques of B. napus pollinated with B. oleracea (AE) and B. rapa (NE) were employed for transcriptome sequencing at 8 and 16 days after pollination. Compared to NE and the parental line of B. napus, more specific differentially expressed genes (DEGs) (1274 and 1698) were obtained in AE and the parental line of B. napus at 8 and 16 days after pollination (DAP). These numbers were 51 and 5.8 times higher than the number of specific DEGs in NE and parental line of B. napus at 8 and 16 DAP, respectively, suggesting more complex transcriptional changes in AE. Most of DEGs in the terms of cell growth and cell wall formation exhibited down-regulated expression patterns (96(down)/131(all) in AE8, 174(down)/235(all) in AE16), while most of DEGs in the processes of photosynthesis, photorespiration, peroxisome, oxidative stress, and systemic acquired resistance exhibited up-regulated expression patterns (222(up)/304(all) in AE8, 214(up)/287(all) in AE16). This is in accordance with a high level of reactive oxygen species (ROS) in the siliques of B. napus pollinated with B. oleracea. Our data suggest that the disorder of plant hormone metabolism, retardation of cell morphogenesis, and the accumulation of ROS may be associated with hybrid incompatibility between B. napus and B. oleracea.

Genome-wide investigation and expression profiling of LOR gene family in rapeseed under salinity and ABA stress

The Brassica napus (B. napus) LOR (Lurp-One-Related) gene family is a little-known gene family characterized by a conserved LOR domain in the proteins. Limited research in Arabidopsis showed that LOR family members played important roles in Hyaloperonospora parasitica (Hpa) defense. Nevertheless, there is a paucity of research investigating the role of the LOR gene family towards their responses to abiotic stresses and hormone treatments. This study encompassed a comprehensive survey of 56 LOR genes in B. napus, which is a prominent oilseed crop that holds substantial economic significance in China, Europe, and North America. Additionally, the study evaluated the expression profiles of these genes in response to salinity and ABA stress. Phylogenetic analysis showed that 56 BnLORs could be divided into 3 subgroups (8 clades) with uneven distribution on 19 chromosomes. 37 out of 56 BnLOR members have experienced segmental duplication and 5 of them have undergone tandem repeats events with strong evidence of purifying selection. Cis-regulatory elements (CREs) analysis indicated that BnLORs involved in process such as light response, hormone response, low temperature response, heat stress response, and dehydration response. The expression pattern of BnLOR family members revealed tissue specificity. RNA-Seq and qRT-PCR were used to validate BnLOR gene expression under temperature, salinity and ABA stress, revealing that most BnLORs showed inducibility. This study enhanced our comprehension of the B. napus LOR gene family and could provide valuable information for identifying and selecting genes for stress resistant breeding.

Fine mapping and causal gene identification of a novel QTL for early flowering by QTL-seq, Target-seq and RNA-seq in spring oilseed rape

A novel quantitative trait locus for early flowering in spring oilseed rape, BnaC08cqDTF, was mapped to an 86-kb region on chromosome C08, and its causal gene, CRY2, was uncovered. Days to flowering is a very important agronomic and adaptive trait of Brassica napus oilseed rape (AACC, 2n = 38). We previously identified BnaC08cqDTF as a novel candidate quantitative trait locus (QTL) for early flowering in spring oilseed rape. Here, we present fine mapping of the locus and a study of its causal gene. Initial mapping was performed by QTL sequencing of DNA pools of BC₃F₂ plants with extreme flowering times derived from crosses between the spring-type cv. No. 4512 (early flowering) and cv. No. 5246 (late flowering), along with fine mapping by target sequencing of the BC₃F₂ and BC₄F₂ populations. Fine mapping narrowed down BnaC08cqDTF to an 86-kb region on chromosome C08. The region harbored fifteen genes. After comparative analyses of the DNA sequences for mutation between A and C syntenic regions and detected by RNA-seq and qRT-PCR between the two parents, we found that BnaC08G0010400ZS harbors an A/G nonsynonymous mutation in exon 3. This single nucleotide polymorphism (SNP) haplotype was also correlated with early flowering in a 256 accession panel. BnaC08G0010400ZS is a homolog of the AT1G04400 gene (CRY2) in Arabidopsis. The analyses of transgenic Arabidopsis verified that BnaC08G0010400ZS is responsible for early flowering. Our results contribute to a better understanding of the genetic control mechanism of early flowering in spring Brassica napus and will promote the breeding for early mature varieties.

Structural variations and environmental specificities of flowering time-related genes in Brassica napus

We found that the flowering time order of accessions in a genetic population considerably varied across environments, and homolog copies of essential flowering time genes played different roles in different locations. Flowering time plays a critical role in determining the life cycle length, yield, and quality of a crop. However, the allelic polymorphism of flowering time-related genes (FTRGs) in Brassica napus, an important oil crop, remains unclear. Here, we provide high-resolution graphics of FTRGs in B. napus on a pangenome-wide scale based on single nucleotide polymorphism (SNP) and structural variation (SV) analyses. A total of 1337 FTRGs in B. napus were identified by aligning their coding sequences with Arabidopsis orthologs. Overall, 46.07% of FTRGs were core genes and 53.93% were variable genes. Moreover, 1.94%, 0.74%, and 4.49% FTRGs had significant presence-frequency differences (PFDs) between the spring and semi-winter, spring and winter, and winter and semi-winter ecotypes, respectively. SNPs and SVs across 1626 accessions of 39 FTRGs underlying numerous published qualitative trait loci were analyzed. Additionally, to identify FTRGs specific to an eco-condition, genome-wide association studies (GWASs) based on SNP, presence/absence variation (PAV), and SV were performed after growing and observing the flowering time order (FTO) of plants in a collection of 292 accessions at three locations in two successive years. It was discovered that the FTO of plants in a genetic population changed a lot across various environments, and homolog copies of some key FTRGs played different roles in different locations. This study revealed the molecular basis of the genotype-by-environment (G × E) effect on flowering and recommended a pool of candidate genes specific to locations for breeding selection.

Fine mapping of genes controlling pigment accumulation in oilseed rape (Brassica napus L.)

Purple/red appearance is one of the common phenotypic variations in leaves, stems, and siliques of oilseed rape (Brassica napus L.) but very rare in flowers. In this study, the causal genes for the purple/red traits in stems and flowers in two accessions of oilseed rape (DH_PR and DH_GC001, respectively) derived from the wide hybridization were fine mapped, and candidate genes were determined by methods combined with bulked segregant analysis (BSA) and RNA-seq analysis. Both traits of purple stem and red flowers were mapped to the locus as AtPAP2 homologous genes (BnaPAP2.C6a and BnaPAP2.A7b, respectively) belonging to the R2R3-MYB family. Sequence comparisons of full-length allelic genes revealed several InDels and SNPs in intron 1 as well as exons, and completely different promoter region of BnaPAP2.C6a and a 211 bp insertion was identified in the promoter region of BnaPAP2.A7b of DH_GC001. Our results not only contribute to a better understanding of anthocyanin inheritance in B. napus, but also provide a useful toolbox for future breeding of cultivars with purple/red traits through the combination of different functional alleles and homologs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01365-5.

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.

QTL mapping for seed density per silique in Brassica napus

Seed density per silique (SDPS) and valid silique length (VSL) are two important yield-influencing traits in rapeseed. SDPS has a direct or indirect effect on rapeseed yield through its effect on seed per silique. In this study, a quantitative trait locus (QTL) for SDPS was detected on chromosome A09 using the QTL-seq approach and confirmed via linkage analysis in the mapping population obtained from 4263 × 3001 cross. Furthermore, one major QTL for SDPS (qSD.A9-1) was mapped to a 401.8 kb genomic interval between SSR markers Nys9A190 and Nys9A531. In the same genomic region, a QTL (qSL.A9) linked to VSL was also detected. The phenotypic variation of qSD.A9-1 and qSL.A9 was 53.1% and 47.6%, respectively. Results of the additive and dominant effects demonstrated that the expression of genes controlling SDPS and VSL were derived from a different parent in this population. Subsequently, we identified 56 genes that included 45 specific genes with exonic (splicing) variants. Further analysis identified specific genes containing mutations that may be related to seed density as well as silique length. These genes could be used for further studies to understand the details of these traits of rapeseed.

GROP: A genomic information repository for oilplants

Biomass energy is an essential component of the agriculture economy and represents an important and particularly significant renewable energy source in the fight against fossil fuel depletion and global warming. The recognition that many plants naturally synthesize hydrocarbons makes these oil plants indispensable resources for biomass energy, and the advancement of next-generation sequencing technology in recent years has now made available mountains of data on plants that synthesize oil. We have utilized a combination of bioinformatic protocols to acquire key information from this massive amount of genomic data and to assemble it into an oil plant genomic information repository, built through website technology, including Django, Bootstrap, and echarts, to create the Genomic Information Repository for Oil Plants (GROP) portal (http://grop.site/) for genomics research on oil plants. The current version of GROP integrates the coding sequences, protein sequences, genome structure, functional annotation information, and other information from 18 species, 22 genome assemblies, and 46 transcriptomes. GROP also provides BLAST, genome browser, functional enrichment, and search tools. The integration of the massive amounts of oil plant genomic data with key bioinformatics tools in a database with a user-friendly interface allows GROP to serve as a central information repository to facilitate studies on oil plants by researchers worldwide.

Transposon insertions within alleles of BnaFT.A2 are associated with seasonal crop type in rapeseed

We identified two new transposon insertions within the promoter of BnaFT.A2 in addition to an existing 288 bp MITE within the second intron. Each insertion event corresponds to a distinct BnaFT.A2 haplotype and is closely associated with established crop seasonal ecotypes. Florigen, encoded by FLOWERING LOCUS T (FT), plays key roles not only as a flowering hormone, but also a universal growth factor affecting several aspects of plant architecture. In rapeseed, BnaFT.A2 has been revealed as one of the major loci associated with flowering time and different ecotypes. However, it is unclear how allelic variations of BnaFT.A2 affect its function in flowering time regulation and beyond. In this study, we confirmed an existing 288 bp miniature inverted-repeat transposable element (MITE) insertion within the second intron and identified two new insertions within the promoter of BnaFT.A2-a 3971 bp CACTA and a 1079 bp Helitron. Each insertion event corresponds to a distinct BnaFT.A2 haplotype and is closely associated with established crop seasonal ecotypes. These alleles have similar tissue-specific expression patterns but discrete transcriptional patterns tightly associated with rapeseed flowering time and ecotype. RNAi lines and mutants of BnaFT.A2 flowered significantly later than controls. Differentially expressed genes (DEGs), identified in transcriptomic profiling of seedling leaves from two loss-of-function mutants (Bnaft.a2-L1 and Bnaft.a2-L2) compared with controls, indicated significant enrichment for hormone metabolic genes and roles related to plant cell wall synthesis and photosynthesis. Plants with loss-of-function BnaFT.A2 had smaller leaves and lower net photosynthetic rate compared to controls. These findings not only further clarify the genetic basis of flowering time variation and ecotype formation in B. napus, but also provide an additional toolbox for genetic improvement of seasonal adaptation and production.

Oil plant genomes: current state of the science

Vegetable oils are an indispensable nutritional component of the human diet as well as important raw materials for a variety of industrial applications such as pharmaceuticals, cosmetics, oleochemicals, and biofuels. Oil plant genomes are highly diverse, and their genetic variation leads to a diversity in oil biosynthesis and accumulation along with agronomic traits. This review discusses plant oil biosynthetic pathways, current state of genome assembly, polyploidy and asymmetric evolution of genomes of oil plants and their wild relatives, and research progress of pan-genomics in oil plants. The availability of complete high-resolution genomes and pan-genomes has enabled the identification of structural variations in the genomes that are associated with the diversity of agronomic and environment fitness traits. These and future genomes also provide powerful tools to understand crop evolution and to harvest the rich natural variations to improve oil crops for enhanced productivity, oil quality, and adaptability to changing environments.

A large-scale population based organelle pan-genomes construction and phylogeny analysis reveal the genetic diversity and the evolutionary origins of chloroplast and mitochondrion in Brassica napus L

BACKGROUND

Allotetraploid oilseed rape (Brassica napus L.) is an important worldwide oil-producing crop. The origin of rapeseed is still undetermined due to the lack of wild resources. Despite certain genetic architecture and phylogenetic studies have been done focus on large group of Brassica nuclear genomes, the organelle genomes information under global pattern is largely unknown, which provide unique material for phylogenetic studies of B. napus. Here, based on de novo assemblies of 1,579 B. napus accessions collected globally, we constructed the chloroplast and mitochondrial pan-genomes of B. napus, and investigated the genetic diversity, phylogenetic relationships of B. napus, B. rapa and B. oleracea.

RESULTS

Based on mitotype-specific markers and mitotype-variant ORFs, four main cytoplasmic haplotypes were identified in our groups corresponding the nap, pol, ole, and cam mitotypes, among which the structure of chloroplast genomes was more conserved without any rearrangement than mitochondrial genomes. A total of 2,092 variants were detected in chloroplast genomes, whereas only 326 in mitochondrial genomes, indicating that chloroplast genomes exhibited a higher level of single-base polymorphism than mitochondrial genomes. Based on whole-genome variants diversity analysis, eleven genetic difference regions among different cytoplasmic haplotypes were identified on chloroplast genomes. The phylogenetic tree incorporating accessions of the B. rapa, B. oleracea, natural and synthetic populations of B. napus revealed multiple origins of B. napus cytoplasm. The cam-type and pol-type were both derived from B. rapa, while the ole-type was originated from B. oleracea. Notably, the nap-type cytoplasm was identified in both the B. rapa population and the synthetic B. napus, suggesting that B. rapa might be the maternal ancestor of nap-type B. napus.

CONCLUSIONS

The phylogenetic results provide novel insights into the organelle genomic evolution of Brassica species. The natural rapeseeds contained at least four cytoplastic haplotypes, of which the predominant nap-type might be originated from B. rapa. Besides, the organelle pan-genomes and the overall variation data offered useful resources for analysis of cytoplasmic inheritance related agronomical important traits of rapeseed, which can substantially facilitate the cultivation and improvement of rapeseed varieties.

Enhancing canola breeding by editing a glucosinolate transporter gene lacking natural variation

A low seed glucosinolate resource was developed in polyploid B. napus using a method that identifies the functions of genes with rare or no genetic variation.

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.

Discrimination of Transgenic Canola (Brassica napus L.) and their Hybrids with B. rapa using Vis-NIR Spectroscopy and Machine Learning Methods

In recent years, the rapid development of genetically modified (GM) technology has raised concerns about the safety of GM crops and foods for human health and the ecological environment. Gene flow from GM crops to other crops, especially in the Brassicaceae family, might pose a threat to the environment due to their weediness. Hence, finding reliable, quick, and low-cost methods to detect and monitor the presence of GM crops and crop products is important. In this study, we used visible near-infrared (Vis-NIR) spectroscopy for the effective discrimination of GM and non-GM Brassica napus, B. rapa, and F1 hybrids (B. rapa X GM B. napus). Initially, Vis-NIR spectra were collected from the plants, and the spectra were preprocessed. A combination of different preprocessing methods (four methods) and various modeling approaches (eight methods) was used for effective discrimination. Among the different combinations, the Savitzky-Golay and Support Vector Machine combination was found to be an optimal model in the discrimination of GM, non-GM, and hybrid plants with the highest accuracy rate (100%). The use of a Convolutional Neural Network with Normalization resulted in 98.9%. The same higher accuracy was found in the use of Gradient Boosted Trees and Fast Large Margin approaches. Later, phenolic acid concentration among the different plants was assessed using GC-MS analysis. Partial least squares regression analysis of Vis-NIR spectra and biochemical characteristics showed significant correlations in their respective changes. The results showed that handheld Vis-NIR spectroscopy combined with chemometric analyses could be used for the effective discrimination of GM and non-GM B. napus, B. rapa, and F1 hybrids. Biochemical composition analysis can also be combined with the Vis-NIR spectra for efficient discrimination.

QTL mapping for low temperature germination in rapeseed

Rapeseed, a major oil crop in the world, is easily affected by low-temperature stress. A low temperature delays seed germination and increases seedling mortality, adversely affecting rapeseed growth and production. In the present study, a tolerant cultivar (Huyou21) was crossed with a susceptible genotype (3429) to develop a mapping population consisting of 574 F₂ progenies and elucidate the genetic mechanisms of seed germination under low temperatures. Two quantitative trait loci (QTL) for low-temperature germination (LTG) were detected, one on chromosome A09 (named qLTGA9-1) and the other on chromosome C01 (named qLTGC1-1), using the QTL-seq approach and confirmed via linkage analysis in the mapping population. Further, qLTGA9-1 was mapped to a 341.86 kb interval between the SSR markers Nys9A212 and Nys9A215. In this region, 69 genes including six specific genes with moderate or high effect function variants were identified based on the Ningyou7 genome sequence. Meanwhile, qLTGC1-1 was mapped onto a 1.31 Mb interval between SSR markers Nys1C96 and Nys1C117. In this region, 133 genes including five specific genes with moderate effect function variants were identified. These specific genes within the two QTL could be used for further studies on cold tolerance and as targets in rapeseed breeding programs.

A genome variation map provides insights into the genetics of walnut adaptation and agronomic traits

BACKGROUND

Common walnut (Juglans regia L.) is one of the top four most consumed nuts in the world due to its health benefits and pleasant taste. Despite its economic importance, the evolutionary history and genetic control of its adaptation and agronomic traits remain largely unexplored.

RESULTS

We report a comprehensive walnut genomic variation map based on whole-genome resequencing of 815 walnut accessions. Evolutionary analyses suggest that Chinese J. regia diverged from J. sigillata with extensive hybridizations after the split of the two species. In contrast to annual crops, the genetic diversity and heterozygous deleterious mutations of Chinese common walnut trees have continued to increase during the improvement process. Selective sweep analyses identify 902 genes uniquely selected in the improved common walnut compared to its progenitor population. Five major-effect loci are identified to be involved in walnut adaptations to temperature, precipitation, and altitude. Genome-wide association studies reveal 27 genomic loci responsible for 18 important agronomic traits, among which JrFAD2 and JrANR are the potentially major-effect causative genes controlling linoleic acid content and color of the endopleura of the nut, respectively.

CONCLUSIONS

The largest genomic resource for walnuts to date has been generated and explored in this study, unveiling their evolutionary history and cracking the genetic code for agronomic traits and environmental adaptation of this economically crucial crop tree.

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.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

A Modified Meiotic Recombination in Brassica napus Largely Improves Its Breeding Efficiency

Simple Summary

The selection of varieties more resilient to disease and climate change requires generating new genetic diversity for breeding. The main mechanism for reshuffling genetic information is through the recombination of chromosomes during meiosis. We showed in oilseed rape (Brassica napus, AACC, 2n = 4x = 38), which is a natural hybrid formed from a cross between turnip (B. rapa, AA, 2n = 2x = 20) and cabbage (B. oleracea, CC, 2n = 2x = 18), that there is significantly more crossovers occurring along the entire A chromosomes in allotriploid AAC (crossbetween B. napus and B. rapa) than in diploid AA or allotetraploid AACC hybrids. We demonstrated that these allotriploid AAC hybrids are highly efficient to introduce new variability within oilseed rape varieties, notably by enabling the introduction of small genomic regions carrying genes controlling agronomically interesting traits.

Abstract

Meiotic recombination is the main tool used by breeders to generate biodiversity, allowing genetic reshuffling at each generation. It enables the accumulation of favorable alleles while purging deleterious mutations. However, this mechanism is highly regulated with the formation of one to rarely more than three crossovers, which are not randomly distributed. In this study, we showed that it is possible to modify these controls in oilseed rape (Brassica napus, AACC, 2n = 4x = 38) and that it is linked to AAC allotriploidy and not to polyploidy per se. To that purpose, we compared the frequency and the distribution of crossovers along A chromosomes from hybrids carrying exactly the same A nucleotide sequence, but presenting three different ploidy levels: AA, AAC and AACC. Genetic maps established with 202 SNPs anchored on reference genomes revealed that the crossover rate is 3.6-fold higher in the AAC allotriploid hybrids compared to AA and AACC hybrids. Using a higher SNP density, we demonstrated that smaller and numerous introgressions of B. rapa were present in AAC hybrids compared to AACC allotetraploid hybrids, with 7.6 Mb vs. 16.9 Mb on average and 21 B. rapa regions per plant vs. nine regions, respectively. Therefore, this boost of recombination is highly efficient to reduce the size of QTL carried in cold regions of the oilseed rape genome, as exemplified here for a QTL conferring blackleg resistance.

Genome structural evolution in Brassica crops

The cultivated Brassica species include numerous vegetable and oil crops of global importance. Three genomes (designated A, B and C) share mesohexapolyploid ancestry and occur both singly and in each pairwise combination to define the Brassica species. With organizational errors (such as misplaced genome segments) corrected, we showed that the fundamental structure of each of the genomes is the same, irrespective of the species in which it occurs. This enabled us to clarify genome evolutionary pathways, including updating the Ancestral Crucifer Karyotype (ACK) block organization and providing support for the Brassica mesohexaploidy having occurred via a two-step process. We then constructed genus-wide pan-genomes, drawing from genes present in any species in which the respective genome occurs, which enabled us to provide a global gene nomenclature system for the cultivated Brassica species and develop a methodology to cost-effectively elucidate the genomic impacts of alien introgressions. Our advances not only underpin knowledge-based approaches to the more efficient breeding of Brassica crops but also provide an exemplar for the study of other polyploids.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

A 598-bp InDel Variation in the Promoter Region of Bna.SOC1.A05 Is Predominantly Present in Winter Type Rapeseeds

During rapeseed domestication and breeding, genetic diversity allowed to adapt it to different eco-geographical regions and to shape its useful traits. Structural variations (SVs), including presence/absence variations (PAVs), are thought to play a major role in the genetic diversity and phenotypic plasticity of rapeseed. In this study, we detected a 598-bp PAV within the promoter region of an Arabidopsis ortholog of a major flowering time gene and a downstream target of FLC, SOC1, which is one of the first genes that are upregulated in rapeseed during vernalization. Further analysis showed that the insertion is present predominantly in winter types while absent in spring types. The 589-bp sequence is present only in the A sub-genome indicating that it originated from Brassica rapa. Since the genomic region around Bna.SOC1.A05 showed a strong reduction in nucleotide diversity, the insertion might represent a larger selected sweep for rapeseed adaptation. Cis-element analysis showed that the insertion contains an ACGTG box, which is the strongest binding motif for the HY5 transcription factor in Arabidopsis. In addition, expression analyses showed that mRNA levels of Bna.SOC1.A05 were lower in accessions carrying the insertion compared to the ones that had no insertion.

Asymmetric epigenome maps of subgenomes reveal imbalanced transcription and distinct evolutionary trends in Brassica napus

The complexity of the epigenome landscape and transcriptional regulation is significantly increased during plant polyploidization, which drives genome evolution and contributes to the increased adaptability to diverse environments. However, a comprehensive epigenome map of Brassica napus is still unavailable. In this study, we performed integrative analysis of five histone modifications, RNA polymerase II occupancy, DNA methylation, and transcriptomes in two B. napus lines (2063A and B409), and established global maps of regulatory elements, chromatin states, and their dynamics for the whole genome (including the An and Cn subgenomes) in four tissue types (young leaf, flower bud, silique, and root) of these two lines. Approximately 65.8% of the genome was annotated with different epigenomic signals. Compared with the Cn subgenome, the An subgenome possesses a higher level of active epigenetic marks and lower level of repressive epigenetic marks. Genes from subgenome-unique regions contribute to the major differences between the An and Cn subgenomes. Asymmetric histone modifications between homeologous gene pairs reflect their biased expression patterns. We identified a novel bivalent chromatin state (with H3K4me1 and H3K27me3) in B. napus that is associated with tissue-specific gene expression. Furthermore, we observed that different types of duplicated genes have discrepant patterns of histone modification and DNA methylation levels. Collectively, our findings provide a valuable epigenetic resource for allopolyploid plants.

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

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

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Genome-Wide Identification and Characterization of the IGT Gene Family in Allotetraploid Rapeseed (Brassica napus L.)

IGT family genes function critically to regulate lateral organ orientation in plants. However, little information is available about this family of genes in Brassica napus. In this study, 27 BnIGT genes were identified on 16 chromosomes and divided into seven clades, namely LAZY1∼LAZY6 and TAC1 (Tiller Angle Control 1), based on their phylogenetic relationships. Duplication analysis revealed that 91.1% of the gene pairs were derived from whole-genome duplication. Most BnIGT genes had a similar structural pattern with one or two very short exons followed by a long and a shorter exon. Common and specific motifs were identified among the seven clades, and motif 1, containing the family-specific GφL(A/T)IGT sequence, was observed in all clades except LAZY5. Three types of cis-elements pertinent to transcription factor binding, light responses, and hormone signaling were detected in the BnIGT promoters. Intriguingly, more than half of the BnIGT genes exhibited no or very low expression in various tissues, and the LAZY1 and TAC1 clade members showed distinct tissue expression preferences. Coexpression analysis revealed that the LAZY1 members had strong associations with cell wall biosynthesis genes. This analysis provides a deeper understanding of the BnIGT gene family and will facilitate further deduction of their role in regulating plant architecture in B. napus.

Identification of QTLs associated with potassium use efficiency and underlying candidate genes by whole-genome resequencing of two parental lines in Brassica napus

Breeding crops that acquire and/or utilize potassium (K) more effectively could reduce the use of K fertilizers. Sixteen traits affecting K use efficiency (KUE) at the seedling stage were investigated in a B. napus double haploid population grown at an optimal K supply (OK) and a low K supply (LK) in a hydroponic culture system. In total, 50 and 62 QTLs associated with these traits were identified at OK and LK, respectively. A total of 25 orthologues of 23 Arabidopsis genes regulating K transport were identified in the confidence intervals of nine QTLs impacting shoot dry weight at LK, and 22 of these showed variations in coding sequences and/or exhibited significant differences in mRNA abundances in roots at LK between the two parental lines. This study provided insights to the genetic basis of KUE in B. napus, which will accelerate the breeding of K-efficient rapeseed cultivars by marker-assisted selection.

Fine Mapping of a Locus Underlying the Ectopic Blade-Like Outgrowths on Leaf and Screening Its Candidate Genes in Rapeseed (Brassica napus L.)

Leaf is an important organ for higher plants, and the shape of it is one of the crucial traits of crops. In this study, we investigated a unique aberrant leaf morphology trait in a mutational rapeseed material, which displayed ectopic blade-like outgrowths on the adaxial side of leaf. The abnormal line 132000B-3 was crossed with the normal line 827-3. Based on the F_{2 : 3} family, we constructed two DNA pools (normal pool and abnormal pool) by the bulked segregant analysis (BSA) method and performed whole genome re-sequencing (WGR), obtaining the single-nucleotide polymorphism (SNP) and insertion/deletion (InDel) data. The SNP-index method was used to calculate the Δ(SNP/InDel-index), and then an association region was identified on chromosome A10 with a length of 5.5 Mbp, harboring 1048 genes totally. Subsequently, the fine mapping was conducted by using the penta-primer amplification refractory mutation system (PARMS), and the associated region was narrowed down to a 35.1-kbp segment, containing only seven genes. These seven genes were then analyzed according to their annotations and finally, BnA10g0422620 and BnA10g0422610, orthologs of LATE MERISTEM IDENTITY1 (LMI1) gene from Arabidopsis and REDUCED COMPLEXITY (RCO) gene from its relative Cardamine hirsuta, respectively, were identified as the candidate genes responding to this blade-like outgrowth trait in rapeseed. This study provides a novel perspective into the leaf formation in Brassica plants.

BnaGVD: A genomic variation database of rapeseed (Brassica napus)

Rapeseed (Brassica napus L.) is a typical polyploid crop and one of the most important oilseed crops worldwide. With the rapid progress on high-throughput sequencing technologies and the reduction of sequencing cost, large-scale genomic data of a specific crop have become available. However, raw sequence data are mostly deposited in the sequence read archive of the National Center of Biotechnology Information (NCBI) and the European Nucleotide Archive (ENA), which is freely accessible to all researchers. Extensive tools for practical purposes should be developed to efficiently utilize these large raw data. Here, we report a web-based rapeseed genomic variation database (BnaGVD, http://rapeseed.biocloud.net/home) from which genomic variations, such as single nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) across a world-wide collection of rapeseed accessions, can be referred. The current release of the BnaGVD contains 34,591,899 high-quality SNPs and 12,281,923 high-quality InDels and provides search tools to retrieve genomic variations and gene annotations across 1,007 accessions of worldwide rapeseed germplasm. We implement a variety of built-in tools (e.g., BnaGWAS, BnaPCA, and BnaStructure) to help users perform in-depth analyses. We recommend this web resource for accelerating studies on the functional genomics and screening of molecular markers for rapeseed breeding.

High-throughput phenotyping accelerates the dissection of the dynamic genetic architecture of plant growth and yield improvement in rapeseed

Rapeseed is the second most important oil crop species and is widely cultivated worldwide. However, overcoming the 'phenotyping bottleneck' has remained a significant challenge. A clear goal of high-throughput phenotyping is to bridge the gap between genomics and phenomics. In addition, it is important to explore the dynamic genetic architecture underlying rapeseed plant growth and its contribution to final yield. In this work, a high-throughput phenotyping facility was used to dynamically screen a rapeseed intervarietal substitution line population during two growing seasons. We developed an automatic image analysis pipeline to quantify 43 dynamic traits across multiple developmental stages, with 12 time points. The time-resolved i-traits could be extracted to reflect shoot growth and predict the final yield of rapeseed. Broad phenotypic variation and high heritability were observed for these i-traits across all developmental stages. A total of 337 and 599 QTLs were identified, with 33.5% and 36.1% consistent QTLs for each trait across all 12 time points in the two growing seasons, respectively. Moreover, the QTLs responsible for yield indicators colocalized with those of final yield, potentially providing a new mechanism of yield regulation. Our results indicate that high-throughput phenotyping can provide novel insights into the dynamic genetic architecture of rapeseed growth and final yield, which would be useful for future genetic improvements in rapeseed.

Transposon insertions within alleles of BnaFLC.A10 and BnaFLC.A2 are associated with seasonal crop type in rapeseed

In Brassicaceae, the requirement for vernalization is conferred by high expression of FLOWERING LOCUS C (FLC). The expression of FLC is known to be repressed by prolonged exposure to cold. Rapeseed (Brassica napus L.) cultivars can be classified into spring, winter, and semi-winter crop types, depending on their respective vernalization requirements. In addition to two known distinct transposon insertion events, here we identified a 4.422 kb hAT and a 5.625 kb long interspersed nuclear element transposon insertion within BnaFLC.A10, and a 810 bp miniature inverted-repeat transposable element (MITE) in BnaFLC.A2. Quantitative PCR demonstrated that these insertions lead to distinct gene expression patterns and contribute differentially to the vernalization response. Transgenic and haplotype analysis indicated that the known 621 bp MITE in the promoter region of BnaFLC.A10 is a transcriptional enhancer that appears to be the main determinant of rapeseed vernalization, and has contributed to the adaptation of rapeseed in winter cultivation environments. In the absence of this transposon insertion, the functional allele of BnaFLC.A2 is a major determinant of vernalization demand. Thus, the combination of BnaFLC.A10 carrying the 621 bp MITE insertion and a functional BnaFLC.A2 appears necessary to establish the winter rapeseed crop phenotype.

High-density genetic variation maps reveal the correlation between asymmetric interspecific introgressions and improvement of agronomic traits in Upland and Pima cotton varieties developed in Xinjiang, China

The two new world tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, are cultivated worldwide and are characterised by a high yield and superior fibre quality, respectively. Historical genetic introgression has been reported between them; however, the existence of introgression and its genetic effects on agronomic traits remain unclear with regard to independent breeding of G. hirsutum (Upland cotton) and G. barbadense (Pima cotton) elite cultivars. We collected 159 G. hirsutum and 70 G. barbadense cultivars developed in Xinjiang, China, along with 30 semi-wild accessions of G. hirsutum, to perform interspecific introgression tests, intraspecific selection analyses and genome-wide association studies (GWAS) with fibre quality and yield component traits in multiple environments. In total, we identified seven interspecific introgression events and 52 selective sweep loci in G. hirsutum, as well as 17 interspecific introgression events and 19 selective sweep loci in G. barbadense. Correlation tests between agronomic traits and introgressions showed that introgression loci were mutually beneficial for the improvement of fibre quality and yield traits in both species. In addition, the phenotypic effects of four interspecific introgression events could be detected by intraspecific GWAS, with Gb_INT13 significantly improving fibre yield in G. barbadense. The present study describes the landscape of genetic introgression and selection between the two species, and highlights the genetic effects of introgression among populations, which can be used for future improvement of fibre yield and quality in G. barbadense and G. hirsutum, respectively.

Genome-wide identification and characterization of SnRK family genes in Brassica napus

BACKGROUND

Sucrose non-fermenting 1 related protein kinases (SnRK) play crucial roles in responding to biotic and abiotic stresses through activating protein phosphorylation pathways. However, little information of SnRK genes was available in Brassica napus, one of important oil crops. Recently, the released sequences of the reference genome of B.napus provide a good chance to perform genome-wide identification and characterization of BnSnRK gene family in the rapeseed.

RESULTS

Totally 114 SnRK genes distributed on 19 chromosomes were identified in the genome of B.napus and classified into three subfamilies on the basis of phylogenetic analysis and the domain types. According to gene structure and motif composition analysis, the BnSnRK sequences showed obvious divergence among three subfamilies. Gene duplication and synteny between the genomes of the rapeseed and Arabidopsis were also analyzed to provide insights into the evolutionary characteristics of BnSnRK family genes. Cis-element analysis revealed that BnSnRKs may response to diverse environmental stresses. Moreover, the expression patterns of BnSnRKs in various tissues and under diverse abiotic stresses were distinct difference. Besides, Single Nucleotide Polymorphisms (SNP) distribution analysis suggests the function disparity of BnSnRK family genes in different genotypes of the rapeseed.

CONCLUSION

We examined genomic structures, evolution features, expression patterns and SNP distribution of 114 BnSnRKs. The results provide valuable information for functional characterization of BnSnRK genes in future studies.

Fine-mapping and transcriptome analysis of a candidate gene controlling plant height in Brassica napus L

BACKGROUND

Brassica napus provides approximately 13-16% of global vegetable oil for human consumption and biodiesel production. Plant height (PH) is a key trait that affects plant architecture, seed yield and harvest index. However, the genetic mechanism of PH in B. napus is poorly understood.

RESULTS

A dwarf mutant df59 was isolated from a large-scale screening of an ethyl methanesulphonate-mutagenized rapeseed variety Ningyou 18. A genetic analysis showed that the dwarfism phenotype was controlled by one semi-dominant gene, which was mapped on C9 chromosome by quantitative trait loci sequencing analysis and designated as BnaDwf.C9. To fine-map BnaDwf.C9, two F2 populations were constructed from crosses between conventional rapeseed cultivars (Zhongshuang 11 and Holly) and df59. BnaDwf.C9 was fine-mapped to the region between single-nucleotide polymorphism (SNP) markers M14 and M4, corresponding to a 120.87-kb interval of the B. napus 'Darmor-bzh' genome. Within this interval, seven, eight and nine annotated or predicted genes were identified in "Darmor-bzh", "Ningyou 7" and "Zhongshuang 11" reference genomes, respectively. In addition, a comparative transcriptome analysis was performed using stem tips from Ningyou 18 and df59 at the stem elongation stage. In total, 3995 differentially expressed genes (DEGs) were identified. Among them, 118 DEGs were clustered in plant hormone-related signal transduction pathways, including 81 DEGs were enriched in auxin signal transduction. Combining the results of fine-mapping and transcriptome analyses, BnaC09g20450D was considered a candidate gene for BnaDwf.C9, which contains a SNP that co-segregated in 4746 individuals. Finally, a PCR-based marker was developed based on the SNP in BnaC09g20450D.

CONCLUSIONS

The combination of quantitative trait loci sequencing, fine-mapping and genome-wide transcriptomic analysis revealed one candidate gene located within the confidence interval of 120.87-kb region. This study provides a new genetic resource for semi-dwarf breeding and new insights into understanding the genetic architecture of PH in B. napus.

Fine mapping and candidate gene analysis of a seed glucosinolate content QTL, qGSL-C2, in rapeseed (Brassica napus L.)

QTL mapping and candidate gene analysis indicate that allelic variations in BnaC2.MYB28 resulted from homeologous exchange and determine difference in seed glucosinolate content. A low seed glucosinolate content has long been an important breeding objective in rapeseed improvement. However, the molecular mechanisms underlying seed GSL content variations remain to be elucidated in allotetraploid Brassica napus. Here, we developed a double haploid population from a cross between two B. napus accessions that possess relatively low, but significantly different seed GSL contents and identified a major QTL, qGSL-C2, on chromosome C02 that explains 30.88-72.87% of the phenotypic variation observed in five environments. Using near-isogenic lines, we further delimited qGSL-C2 to a physical region of 49 kb on the B. rapa chromosome A02 which is highly homologous to the target C02 interval. Among five candidate genes, BnaC2.MYB28, a homologue of the Arabidopsis MYB28 encoding a putative R2R3-MYB-type transcription factor functioning in aliphatic methionine-derived GSL synthesis, was most likely to be the target gene underlying the QTL. Sequence analysis revealed multiple insertion/deletion and SNP variations in the genomic region between the alleles of the NILs. Furthermore, the allelic variations in BnaC2.MYB28 in the natural B. napus population were significantly associated with seed GSL content. Remarkably, the phylogenetic analysis and sequence comparison suggested that while the BnaC2.MYB28 allele from the parental line G120 was inherited from B. oleracea BolC2.MYB28, its counterpart from the other parent, 9172, most likely evolved from B. rapa BraA2.MYB28 via possible homeologous exchange. Our study promotes greater understanding of the molecular regulation of seed GSL content and provides useful molecular markers for seed GSL improvement in B. napus.

Eight high-quality genomes reveal pan-genome architecture and ecotype differentiation of Brassica napus

Rapeseed (Brassica napus) is the second most important oilseed crop in the world but the genetic diversity underlying its massive phenotypic variations remains largely unexplored. Here, we report the sequencing, de novo assembly and annotation of eight B. napus accessions. Using pan-genome comparative analysis, millions of small variations and 77.2-149.6 megabase presence and absence variations (PAVs) were identified. More than 9.4% of the genes contained large-effect mutations or structural variations. PAV-based genome-wide association study (PAV-GWAS) directly identified causal structural variations for silique length, seed weight and flowering time in a nested association mapping population with ZS11 (reference line) as the donor, which were not detected by single-nucleotide polymorphisms-based GWAS (SNP-GWAS), demonstrating that PAV-GWAS was complementary to SNP-GWAS in identifying associations to traits. Further analysis showed that PAVs in three FLOWERING LOCUS C genes were closely related to flowering time and ecotype differentiation. This study provides resources to support a better understanding of the genome architecture and acceleration of the genetic improvement of B. napus.