Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 2. Identification of alleles from unadapted germplasm

Envirotyping within a multi-environment trial allowed identifying genetic determinants of winter oilseed rape yield stability

KEY MESSAGE

A comprehensive environmental characterization allowed identifying stable and interactive QTL for seed yield: QA09 and QC09a were detected across environments; whereas QA07a was specifically detected on the most stressed environments. A main challenge for rapeseed consists in maintaining seed yield while adapting to climate changes and contributing to environmental-friendly cropping systems. Breeding for cultivar adaptation is one of the keys to meet this challenge. Therefore, we propose to identify the genetic determinant of seed yield stability for winter oilseed rape using GWAS coupled with a multi-environmental trial and to interpret them in the light of environmental characteristics. Due to a comprehensive characterization of a multi-environmental trial using 79 indicators, four contrasting envirotypes were defined and used to identify interactive and stable seed yield QTL. A total of four QTLs were detected, among which, QA09 and QC09a, were stable (detected at the multi-environmental trial scale or for different envirotypes and environments); and one, QA07a, was specifically detected into the most stressed envirotype. The analysis of the molecular diversity at QA07a showed a lack of genetic diversity within modern lines compared to older cultivars bred before the selection for low glucosinolate content. The results were discussed in comparison with other studies and methods as well as in the context of breeding programs.

Epistasis and pleiotropy-induced variation for plant breeding

Epistasis refers to nonallelic interaction between genes that cause bias in estimates of genetic parameters for a phenotype with interactions of two or more genes affecting the same trait. Partitioning of epistatic effects allows true estimation of the genetic parameters affecting phenotypes. Multigenic variation plays a central role in the evolution of complex characteristics, among which pleiotropy, where a single gene affects several phenotypic characters, has a large influence. While pleiotropic interactions provide functional specificity, they increase the challenge of gene discovery and functional analysis. Overcoming pleiotropy-based phenotypic trade-offs offers potential for assisting breeding for complex traits. Modelling higher order nonallelic epistatic interaction, pleiotropy and non-pleiotropy-induced variation, and genotype × environment interaction in genomic selection may provide new paths to increase the productivity and stress tolerance for next generation of crop cultivars. Advances in statistical models, software and algorithm developments, and genomic research have facilitated dissecting the nature and extent of pleiotropy and epistasis. We overview emerging approaches to exploit positive (and avoid negative) epistatic and pleiotropic interactions in a plant breeding context, including developing avenues of artificial intelligence, novel exploitation of large-scale genomics and phenomics data, and involvement of genes with minor effects to analyse epistatic interactions and pleiotropic quantitative trait loci, including missing heritability.

Downregulation of the expression of subgenomic chromosome A7 genes promotes plant height in resynthesized allopolyploid Brassica napus

KEY MESSAGE

Homoeolog expression bias and the gene dosage effect induce downregulation of genes on chromosome A7, causing a significant increase in the plant height of resynthesized allopolyploid Brassica napus. Gene expression levels in allopolyploid plants are not equivalent to the simple average of the expression levels in the parents and are associated with several non-additive expression phenomena, including homoeolog expression bias. However, hardly any information is available on the effect of homoeolog expression bias on traits. Here, we studied the effects of gene expression-related characteristics on agronomic traits using six isogenic resynthesized Brassica napus lines across the first ten generations. We found a group of genes located on chromosome A7 whose expression levels were significantly negatively correlated with plant height. They were expressed at significantly lower levels than their homoeologous genes, owing to allopolyploidy rather than inheritance from parents. Homoeolog expression bias resulted in resynthesized allopolyploids with a plant height similar to their female Brassica oleracea parent, but significantly higher than that of the male Brassica rapa parent. Notably, aneuploid lines carrying monosomic and trisomic chromosome A7 had the highest and lowest plant heights, respectively, due to changes in the expression bias of homoeologous genes because of alterations in the gene dosage. These findings suggest that the downregulation of the expression of homoeologous genes on a single chromosome can result in the partial improvement of traits to a significant extent in the nascent allopolyploid B. napus.

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

KEY MESSAGE

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

Discovery of common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits in Brassica napus L

KEY MESSAGE

Common loci and candidate genes for controlling salt-alkali tolerance and yield-related traits were identified in Brassica napus combining QTL mapping with transcriptome under salt and alkaline stresses. The yield of rapeseed (Brassica napus L.) is determined by multiple yield-related traits, which are susceptible to environmental factors. Many yield-related quantitative trait loci (QTLs) have been reported in Brassica napus; however, no studies have been conducted to investigate both salt-alkali tolerance and yield-related traits simultaneously. Here, specific-locus amplified fragment sequencing (SLAF-seq) technologies were utilized to map the QTLs for salt-alkali tolerance and yield-related traits. A total of 65 QTLs were identified, including 30 QTLs for salt-alkali tolerance traits and 35 QTLs for yield-related traits, accounting for 7.61-27.84% of the total phenotypic variations. Among these QTLs, 18 unique QTLs controlling two to four traits were identified by meta-analysis. Six novel and unique QTLs were detected for salt-alkali tolerance traits. By comparing these unique QTLs for salt-alkali tolerance traits with those previously reported QTLs for yield-related traits, seven co-localized chromosomal regions were identified on A09 and A10. Combining QTL mapping with transcriptome of two parents under salt and alkaline stresses, thirteen genes were identified as the candidates controlling both salt-alkali tolerance and yield. These findings provide useful information for future breeding of high-yield cultivars resistant to alkaline and salt stresses.

Genome-Wide Association Studies of Salt Tolerance at the Seed Germination Stage and Yield-Related Traits in Brassica napus L

Salt stress severely affects crop growth and development and reduces the yield of Brassica napus. Exploring natural genetic variations for high salt tolerance in B. napus seedlings is an effective approach to improve productivity under salt stress. Using 10,658 high-quality single nucleotide polymorphic (SNP) markers developed by specific-locus amplified fragment sequencing (SLAF-seq) technology, genome-wide association studies (GWAS) were performed to investigate the genetic basis of salt tolerance and yield-related traits of B. napus. The results revealed that 77 and 497 SNPs were significantly associated with salt tolerance and yield-related traits, of which 40 and 58 SNPs were located in previously reported QTLs/SNPs, respectively. We identified nineteen candidate genes orthologous with Arabidopsis genes known to be associated with salt tolerance and seven potential candidates controlling both salt tolerance and yield. Our study provides a novel genetic resource for the breeding of high-yield cultivars resistant to salt stress.

Transcriptomic analysis of rapeseed (Brassica napus. L.) seed development in Xiangride, Qinghai Plateau, reveals how its special eco-environment results in high yield in high-altitude areas

As one of the most important oil crops, rapeseed (Brassica napus) is cultivated worldwide to produce vegetable oil, animal feed, and biodiesel. As the population grows and the need for renewable energy increases, the breeding and cultivation of high-yield rapeseed varieties have become top priorities. The formation of a high rapeseed yield is so complex because it is influenced not only by genetic mechanisms but also by many environmental conditions, such as climatic conditions and different farming practices. Interestingly, many high-yield areas are located in special eco-environments, for example, in the high-altitude Xiangride area of the Qinghai Plateau. However, the molecular mechanisms underlying the formation of high yields in such a special eco-environment area remain largely unknown. Here, we conducted field yield analysis and transcriptome analysis in the Xiangride area. Compared with the yield and environmental factors in the Xinning area (a low-yielding area), we found that the relatively longer daylight length is the key to high rapeseed yield in the Xiangride area, which leads up to a 52.1% increase in rapeseed yield, especially the increase in thousand seed weight and silique number (SN). Combined with transcriptome H-cluster analysis and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional analyses, we can assume that the grain development of rapeseed in the Xiangride area is ahead of schedule and lasts for a long time, leading to the high-yield results in the Xiangride area, confirmed by the expression analysis by quantitative real-time polymerase chain reaction (qRT-PCR) of yield-related genes. Our results provide valuable information for further exploring the molecular mechanism underlying high yield in special ecological environments and provide a helpful reference for studying seed development characteristics in special-producing regions for Brassica napus.

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

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

Genetic Analysis of Heterosis for Yield Influencing Traits in Brassica juncea Using a Doubled Haploid Population and Its Backcross Progenies

The exploitation of heterosis through hybrid breeding is one of the major breeding objectives for productivity increase in crop plants. This research analyzes the genetic basis of heterosis in Brassica juncea by using a doubled haploid (DH) mapping population derived from F1 between two heterotic inbred parents, one belonging to the Indian and the other belonging to the east European gene pool, and their two corresponding sets of backcross hybrids. An Illumina Infinium Brassica 90K SNP array-based genetic map was used to identify yield influencing quantitative trait loci (QTL) related to plant architecture, flowering, and silique- and seed-related traits using five different data sets from multiple trials, allowing the estimation of additive and dominance effects, as well as digenic epistatic interactions. In total, 695 additive QTL were detected for the 14 traits in the three trials using five data sets, with overdominance observed to be the predominant type of effect in determining the expression of heterotic QTL. The results indicated that the design in the present study was efficient for identifying common QTL across multiple trials and populations, which constitute a valuable resource for marker-assisted selection and further research. In addition, a total of 637 epistatic loci were identified, and it was concluded that epistasis among loci without detectable main effects plays an important role in controlling heterosis in yield of B. juncea.

Integrative analysis of GWAS and transcriptome to reveal novel loci regulation flowering time in semi-winter rapeseed

Flowering is an important turning point from vegetative growth to reproductive growth, and vernalization is an essential condition for the flowering of annual winter plants. To investigate the genetic architecture of flowering time in rapeseed, we used the 60 K Brassica Infinium SNP array to perform a genome-wide analysis of haplotype blocks associated with flowering time in 203 Chinese semi-winter rapeseed inbred lines. Twenty-one haplotype regions carrying one or more candidate genes showed a significant association with flowering time. Interestingly, we detected a SNP (Bn-scaff_22728_1-p285715) located in exon 3 of the BnVIN3-C03 gene that showed a significant association with flowering time on chromosome C03. Based on the SNP alleles A and G, two groups of accessions with early and late flowering time phenotypes were selected, respectively, and PCR amplification and gene expression analysis were combined to reveal the structural variation of the BnVIN3-C03 gene that affected flowering time. Moreover, we found that BnVIN3-C03 inhibited the expression of BnFLC-A02, BnFLC-A03.1, BnFLC-A10 and BnFLC-C03.1, thus modulating the flowering time of Brassica napus. This result provides insight into the genetic improvement of flowering time in B. napus.

Both overlapping and independent loci underlie seed number per pod and seed weight in Brassica napus by comparative quantitative trait loci analysis

Seed number per pod (SNPP) and seed weight (SW) are two components of seed yield in rapeseed (Brassica napus). Here, a natural population of rapeseed was employed for genome-wide association analysis for SNPP and SW across multi-years. A total of 101 and 77 SNPs significantly associated with SNPP and SW with the phenotypic variances (R2) ranging from 1.35 to 29.47% and from 0.78 to 34.58%, respectively. And 43 and 33 homologs of known genes from model plants were located in the 65 and 49 haplotype blocks (HBs) for SNPP and SW, respectively. Notably, we found 5 overlapping loci and 3 sets of loci with collinearity for both SNPP and SW, of which 4 overlapping loci harbored the haplotypes with the same direction of genetic effects on SNPP and SW, indicating high possibility to simultaneously improve SNPP and SW in rapeseed. Our findings revealed both overlapping and independent loci controlling seed number per pod and seed weight in rapeseed.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-021-01232-1.

Mapping of a major QTL controlling plant height using a high-density genetic map and QTL-seq methods based on whole-genome resequencing in Brassica napus

Plant height is a crucial element related to plant architecture that influences the seed yield of oilseed rape (Brassica napus L.). In this study, we isolated a natural B. napus mutant, namely a semi-dwarf mutant (sdw-e), which exhibits a 30% reduction in plant height compared with Zhongshuang 11-HP (ZS11-HP). Quantitative trait locus sequencing (QTL-seq) was conducted using two extreme DNA bulks in F₂ populations in Wuchang-2017 derived from ZS11-HP × sdw-e to identify QTLs associated with plant height. The result suggested that two QTL intervals were located on chromosome A10. The F₂ population consisting of 200 individuals in Yangluo-2018 derived from ZS11-HP × sdw-e was used to construct a high-density linkage map using whole-genome resequencing. The high-density linkage map harbored 4323 bin markers and covered a total distance of 2026.52 cM with an average marker interval of 0.47 cM. The major QTL for plant height named qPHA10 was identified on linkage group A10 by interval mapping and composite interval mapping methods. The major QTL qPHA10 was highly consistent with the QTL-seq results. And then, we integrated the variation sites and expression levels of genes in the major QTL interval to predict the candidate genes. Thus, the identified QTL and candidate genes could be used in marker-assisted selection for B. napus breeding in the future.

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 transition to flowering in winter rapeseed during vernalization

Flowering time is a major determinant of adaptation, fitness and yield in the allopolyploid species rapeseed (Brassica napus). Despite being a close relative to Arabidopsis thaliana, little is known about the timing of floral transition and the genes that govern this process. Winter, semi-winter and spring type plants have important life history characteristics that differ in vernalization requirements for flowering and are important for growing rapeseed in different regions of the world. In this study, we investigated the timing of vernalization-driven floral transition in winter rapeseed and the effect of photoperiod and developmental age on flowering time and vernalization responsiveness. Microscopy and whole transcriptome analyses at the shoot apical meristems of plants grown under controlled conditions showed that floral transition is initiated within few weeks of vernalization. Certain Bna.SOC1 and Bna.SPL5 homeologs were among the induced genes, suggesting that they are regulating the timing of cold-induced floral transition. Moreover, the flowering response of plants with shorter pre-vernalization period correlated with a delayed expression of Bna.SOC1 and Bna.SPL5 genes. In essence, this study presents a detailed analysis of vernalization-driven floral transition and the aspects of juvenility and dormancy and their effect on flowering time in rapeseed.

Genetic dissection of thousand-seed weight and fine mapping of cqSW.A03-2 via linkage and association analysis in rapeseed (Brassica napus L.)

KEY MESSAGE: cqSW.A03-2, one of the six identified quantitative trait loci associated with thousand-seed weight in rapeseed, is mapped to a 61.6-kb region on chromosome A03 and corresponds to the candidate gene BnaA03G37960D. Seed weight is an important factor that determines the seed yield of oilseed rape (Brassica napus L.). To elucidate the genetic mechanism of thousand-seed weight (TSW), quantitative trait locus (QTL) mapping was conducted using a double haploid population derived from the cross between an elite line ZY50 and a pol cytoplasmic male sterility restorer line 7-5. The genetic basis of TSW was dissected into six major QTLs. One major QTL denoted as cqSW.A03-2, which explained 8.46-13.70% of the phenotypic variation, was detected across multiple environments. To uncover the genetic basis of cqSW.A03-2, a set of near-isogenic lines were developed. Based on the test of self-pollinated progenies, cqSW.A03-2 was identified as a single Mendelian factor and the ZY50 allele at cqSW.A03-2 showed a positive effect on TSW. Fine mapping delimited the cqSW.A03-2 locus into a 61.6-kb region, and 18 genes within this region were predicted. Candidate gene association analysis and expression analysis indicated that a histidine kinase gene (BnaA03G37960D) is likely to be the candidate gene for the cqSW.A03-2 locus. Our results may contribute to a better understanding of the molecular mechanism of seed weight regulation and promote the breeding program for yield improvement in rapeseed.

Validation and Characterization of a Seed Number Per Silique Quantitative Trait Locus qSN.A7 in Rapeseed (Brassica napus L.)

Seed number is a key character/trait tightly related to the plant fitness/evolution and crop domestication/improvement. The seed number per silique (SNPS) shows a huge variation from several to more than 30, however the underlying regulatory mechanisms are poorly known, which has hindered its improvement. To answer this question, several representative lines with extreme SNPS were previously subjected to systematic genetic and cytological analyses. The results showed that the natural variation of seed number per silique is mainly controlled by maternal and embryonic genotype, which are co-determined by ovule number per ovary, fertile ovule ratio, ovule fertilization rate, and fertilized ovule development rate. More importantly, we also mapped two repeatable quantitative trait loci (QTLs) for SNPS using the F2:3 population derived from Zhongshuang11 and No. 73290, of which the major QTL qSN.A6 has been fine-mapped. In the current study, the near-isogenic lines (NILs) of qSN.A7 were successfully developed by the successive backcross of F1 with Zhongshuang11. First, the effect of qSN.A7 was validated by evaluating the SNPS of two types of homozygous NILs from BC3F2 population, which showed a significant difference of 2.23 on average. Then, qSN.A7 was successfully fine-mapped from the original 4.237 to 1.389 Mb, using a BC4F2 segregating population of 2,551 individuals. To further clarify the regulatory mechanism of qSN.A7, the two types of homologous NILs were subjected to genetic and cytological analyses. The results showed that the difference in SNPS between the two homologous NILs was determined by the embryonic genotypic effect. Highly accordant with this, no significant difference was observed in ovule number per ovary, ovule fertility, fertilization rate, and pollen fertility between the two homologous NILs. Therefore, the regulatory mechanism of qSN.A7 is completely different from the cloned qSS.C9 and qSN.A6. These results will advance the understanding of SNPS and facilitate gene cloning and molecular breeding in Brassica napus.

Regulation and Subfunctionalization of Flowering Time Genes in the Allotetraploid Oil Crop Brassica napus

Flowering is a vulnerable, but crucial phase in building crop yield. Proper timing of this period is therefore decisive in obtaining optimal yields. However, genetic regulation of flowering integrates many different environmental signals and is therefore extremely complex. This complexity increases in polyploid crops which carry two or more chromosome sets, like wheat, potato or rapeseed. Here, I summarize the current state of knowledge about flowering time gene copies in rapeseed (Brassica napus), an important oil crop with a complex polyploid history and a close relationship to Arabidopsis thaliana. The current data show a high demand for more targeted studies on flowering time genes in crops rather than in models, allowing better breeding designs and a deeper understanding of evolutionary principles. Over evolutionary time, some copies of rapeseed flowering time genes changed or lost their original role, resulting in subfunctionalization of the respective homologs. For useful applications in breeding, such patterns of subfunctionalization need to be identified and better understood.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Fine mapping of a silique length- and seed weight-related gene in Brassica napus

Using microarray analysis combined with map-based cloning, a major locus positively regulating SL and SW was mapped to a 98.47 kb interval on A09 in rapeseed. In rapeseed, seed yield is closely associated with silique-related traits such as silique length (SL) and seed weight (SW). Previously identified quantitative trait loci (QTLs) revealed that SL and SW are complex traits and many QTLs overlap. However, the genetic characterization of the association between SL and SW is poorly understood. In the present study, a BC₃F₃ near isogenic line developed from a short silique plant and the long silique cultivar 'ZS11' was analyzed to identify the locus related to SL. Map-based cloning indicated that a major locus acting as a single Mendelian factor was mapped to a 98.47 kb region on chromosome A09. BLAST analysis and DNA sequencing showed SNP variations and a fragment replacement in the upstream region of the candidate gene BnaA09g55530D may alter gene expression and influence SL. The results showed that this SL locus may also positively affect SW as well as in the 186 rapeseed accessions identified by the associated markers. Therefore, selecting plants with appropriate SL and developing functional markers for the associated gene could play important roles in the molecular breeding of high-yield rapeseed varieties.

The vernalisation regulator FLOWERING LOCUS C is differentially expressed in biennial and annual Brassica napus

Plants in temperate areas evolved vernalisation requirement to avoid pre-winter flowering. In Brassicaceae, a period of extended cold reduces the expression of the flowering inhibitor FLOWERING LOCUS C (FLC) and paves the way for the expression of downstream flowering regulators. As with all polyploid species of the Brassicaceae, the model allotetraploid Brassica napus (rapeseed, canola) is highly duplicated and carries 9 annotated copies of Bna.FLC. To investigate whether these multiple homeologs and paralogs have retained their original function in vernalisation or undergone subfunctionalisation, we compared the expression patterns of all 9 copies between vernalisation-dependent (biennial, winter type) and vernalisation-independent (annual, spring type) accessions, using RT-qPCR with copy-specific primers and RNAseq data from a diversity set. Our results show that only 3 copies – Bna.FLC.A03b, Bna.FLC.A10 and to some extent Bna.FLC.C02 – are differentially expressed between the two growth types, showing that expression of the other 6 copies does not correlate with growth type. One of those 6 copies, Bna.FLC.C03b, was not expressed at all, indicating a pseudogene, while three further copies, Bna.FLC.C03a and Bna.FLC.C09ab, did not respond to cold treatment. Sequence variation at the COOLAIR binding site of Bna.FLC.A10 was found to explain most of the variation in gene expression. However, we also found that Bna.FLC.A10 expression is not fully predictive of growth type.

Different copies of SENSITIVITY TO RED LIGHT REDUCED 1 show strong subfunctionalization in Brassica napus

BACKGROUND

Correct timing of flowering is critical for plants to produce enough viable offspring. In Arabidopsis thaliana (Arabidopsis), flowering time is regulated by an intricate network of molecular signaling pathways. Arabidopsis srr1-1 mutants lacking SENSITIVITY TO RED LIGHT REDUCED 1 (SRR1) expression flower early, particularly under short day (SD) conditions (1). SRR1 ensures that plants do not flower prematurely in such non-inductive conditions by controlling repression of the key florigen FT. Here, we have examined the role of SRR1 in the closely related crop species Brassica napus.

RESULTS

Arabidopsis SRR1 has five homologs in Brassica napus. They can be divided into two groups, where the A02 and C02 copies show high similarity to AtSRR1 on the protein level. The other group, including the A03, A10 and C09 copies all carry a larger deletion in the amino acid sequence. Three of the homologs are expressed at detectable levels: A02, C02 and C09. Notably, the gene copies show a differential expression pattern between spring and winter type accessions of B. napus. When the three expressed gene copies were introduced into the srr1-1 background, only A02 and C02 were able to complement the srr1-1 early flowering phenotype, while C09 could not. Transcriptional analysis of known SRR1 targets in Bna.SRR1-transformed lines showed that CYCLING DOF FACTOR 1 (CDF1) expression is key for flowering time control via SRR1.

CONCLUSIONS

We observed subfunctionalization of the B. napus SRR1 gene copies, with differential expression between early and late flowering accessions of some Bna.SRR1 copies. This suggests involvement of Bna.SRR1 in regulation of seasonal flowering in B. napus. The C09 gene copy was unable to complement srr1-1 plants, but is highly expressed in B. napus, suggesting specialization of a particular function. Furthermore, the C09 protein carries a deletion which may pinpoint a key region of the SRR1 protein potentially important for its molecular function. This is important evidence of functional domain annotation in the highly conserved but unique SRR1 amino acid sequence.

Joint QTL mapping and transcriptome sequencing analysis reveal candidate flowering time genes in Brassica napus L

Background

Optimum flowering time is a key agronomic trait in Brassica napus. To investigate the genetic architecture and genetic regulation of flowering time in this important crop, we conducted quantitative trait loci (QTL) analysis of flowering time in a recombinant inbred line (RIL) population, including lines with extreme differences in flowering time, in six environments, along with RNA-Seq analysis.

Results

We detected 27 QTLs distributed on eight chromosomes among six environments, including one major QTL on chromosome C02 that explained 11–25% of the phenotypic variation and was stably detected in all six environments. RNA-Seq analysis revealed 105 flowering time-related differentially expressed genes (DEGs) that play roles in the circadian clock/photoperiod, autonomous pathway, and hormone and vernalization pathways. We focused on DEGs related to the regulation of flowering time, especially DEGs in QTL regions.

Conclusions

We identified 45 flowering time-related genes in these QTL regions, eight of which are DEGs, including key flowering time genes PSEUDO RESPONSE REGULATOR 7 (PRR7) and FY (located in a major QTL region on C02). These findings provide insights into the genetic architecture of flowering time in B. napus.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5356-8) contains supplementary material, which is available to authorized users.

QTL and Candidate Gene Identification for Silique Length Based on High-Dense Genetic Map in Brassica napus L

Silique length (SL) is an important yield trait and positively correlates with seeds per silique and seed weight. In the present study, two double haploid (DH) populations, established from crosses Zhongshuang11 × R11 (ZR) and R1 × R2 (RR), containing 280 and 95 DH lines, respectively, were used to map quantitative trait loci (QTL) for SL. A high-dense genetic map from ZR population was constructed comprising 14,658 bins on 19 linkage groups, with map length of 2,198.85 cM and an average marker distance of 0.15 cM. Genetic linkage map from RR population was constructed by using 2,046 mapped markers anchored to 19 chromosomes with 2,217-cM map length and an average marker distance of 1.08 cM. Major QTL qSL_ZR_A09 and qSL_RR_A09b on A09 were identified from ZR and RR populations, respectively. Both QTL could be stably detected in four environments. QTL qSL_RR_A09b and qSL_ZR_A09 were located on 68.5-70.8 cM and 91.33-91.94 cM interval with R² values of 14.99-39.07% and 15.00-20.36% in RR and ZR populations, respectively. Based on the physical positions of single nucleotide polymorphism (SNP) markers flanking qSL_ZR_A09 and gene annotation in Arabidopsis, 26 genes were identified with SNP/Indel variation between parents and two genes (BnaA09g41180D and BnaA09g41380D) were selected as the candidate genes. Expression analysis further revealed BnaA09g41180D, encoding homologs of Arabidopsis fasciclin-like arabinogalactan proteins (FLA3), as the most promising candidate gene for qSL_ZR_A09. The QTL identification and candidate gene analysis will provide new insight into the genomic regions controlling SL in Brassica napus as well as candidate genes underlying the QTL.

The effect of INDEHISCENT point mutations on silique shatter resistance in oilseed rape (Brassica napus)

This study elucidates the influence of indehiscent mutations on rapeseed silique shatter resistance. A phenotype with enlarged replum-valve joint area and altered cell dimensions in the dehiscence zone is described. Silique shattering is a major factor reducing the yield stability of oilseed rape (Brassica napus). Attempts to improve shatter resistance often include the use of mutations in target genes identified from Arabidopsis (Arabidopsis thaliana). A variety of phenotyping methods assessing the level of shatter resistance were previously described. However, a comparative and comprehensive evaluation of the methods has not yet been undertaken. We verified the increase of shatter resistance in indehiscent double knock-down mutants obtained by TILLING with a systematic approach comparing three independent phenotyping methods. A positive correlation of silique length and shatter resistance was observed and accounted for in the analyses. Microscopic studies ruled out the influence of different lignification patterns. Instead, we propose a model to explain increased shattering resistance of indehiscent rapeseed mutants by altered cell shapes and sizes within the contact surfaces of replum and valves.

Stable QTLs for Plant Height on Chromosome A09 Identified From Two Mapping Populations in Peanut (Arachis hypogaea L.)

The peanut (Arachis hypogaea L.) is an important grain legume extensively cultivated worldwide, supplying edible oil and protein for human consumption. As in many other crops, plant height is a crucial factor in determining peanut architecture traits and has a unique effect on resistance to lodging and efficiency of mechanized harvesting as well as yield. Currently, the genetic basis underlying plant height remains unclear in peanut, which have hampered marker-assisted selection in breeding. In this study, we conducted a quantitative trait locus (QTL) analysis for peanut plant height by using two recombinant inbred line (RIL) populations including "Yuanza 9102 × Xuzhou 68-4 (YX)" and "Xuhua 13 × Zhonghua 6 (XZ)". In the YX population, 38 QTLs including 10 major QTLs from 9 chromosomes were detected in 4 environments, and 8 consensus QTLs integrated by meta-analysis expressed stably across multiple environments. In the XZ population, 3 major QTLs and seven minor QTLs from 6 chromosomes were detected across 3 environments. Generally, most major QTLs from the two populations were located on pseudomolecule chromosome 9 of Arachis duranesis (A09), indicating there would be key genes on A09 controlling plant height. Further analysis revealed that qPHA09.1a from the XZ population and one consensus QTL, cqPHA09.d from the YX population were co-localized in a reliable 3.4 Mb physical interval on A09, which harbored 161 genes including transcription factors and enzymes related to signaling transduction and cell wall formation. The major and stable QTLs identified in this study may be useful for further gene cloning and identification of molecular markers applicable for breeding.

QTL Alignment for Seed Yield and Yield Related Traits in Brassica napus

Worldwide consumption of oil is increasing with the growing population in need for edible oil and the expansion of industry using biofuels. Then, demand for high yielding varieties of oil crops is always increasing. Brassica napus (rapeseed) is one of the most important oil crop in the world, therefore, increasing rapeseed yield through breeding is inevitable in order to cater for the high demand of vegetable oil and high-quality protein for live stocks. Quantitative trait loci (QTL) analysis is a powerful tool to identify important loci and which is also valuable for molecular marker assisted breeding. Seed-yield (SY) is a complex trait that is controlled by multiple loci and is affected directly by seed weight, seeds per silique and silique number. Some yield-related traits, such as plant height, biomass yield, flowering time, and so on, also affect the SY indirectly. This study reports the assembly of QTLs identified for seed-yield and yield-related traits in rapeseed, in one unique map. A total of 972 QTLs for seed-yield and yield-related were aligned into the physical map of B. napus Darmor-bzh and 92 regions where 198 QTLs overlapped, could be discovered on 16 chromosomes. Also, 147 potential candidate genes were discovered in 65 regions where 131 QTLs overlapped, and might affect nine different traits. At the end, interaction network of candidate genes was studied, and showed nine genes that could highly interact with the other genes, and might have more influence on them. The present results would be helpful to develop molecular markers for yield associated traits and could be used for breeding improvement in B. napus.

Genome-Wide Association Study Reveals Both Overlapping and Independent Genetic Loci to Control Seed Weight and Silique Length in Brassica napus

Seed weight (SW) is one of three determinants of seed yield, which positively correlates with silique length (SL) in Brassica napus (rapeseed). However, the genetic mechanism underlying the relationship between seed weight (SW) and silique length (SL) is largely unknown at present. A natural population comprising 157 inbred lines in rapeseed was genotyped by whole-genome re-sequencing and investigated for SW and SL over four years. The genome-wide association study identified 20 SNPs in significant association with SW on A01, A04, A09, C02, and C06 chromosomes and the phenotypic variation explained by a single locus ranged from 11.85% to 34.58% with an average of 25.43%. Meanwhile, 742 SNPs significantly associated with SL on A02, A03, A04, A07, A08, A09, C01, C03, C04, C06, C07, and C08 chromosomes were also detected and the phenotypic variation explained by a single locus ranged from 4.01 to 48.02% with an average of 33.33%, out of which, more than half of the loci had not been reported in the previous studies. There were 320 overlapping or linked SNPs for both SW and SL on A04, A09, and C06 chromosomes. It indicated that both overlapping and independent genetic loci controlled both SW and SL in B. napus. On the haplotype block on A09 chromosome, the allele variants of a known gene BnaA.ARF18.a controlling both SW and SL were identified in the natural population by developing derived cleaved amplified polymorphic sequence (dCAPS) markers. These findings are valuable for understanding the genetic mechanism of SW and SL and also for rapeseed molecular breeding programs.

Major Co-localized QTL for Plant Height, Branch Initiation Height, Stem Diameter, and Flowering Time in an Alien Introgression Derived Brassica napus DH Population

Plant height (PH), branch initiation height (BIH), and stem diameter (SD) are three stem-related traits that play crucial roles in plant architecture and lodging resistance. Herein, we show one doubled haploid (DH) population obtained from a cross between Y689 (one Capsella bursa-pastoris derived Brassica napus intertribal introgression) and Westar (B. napus cultivar) that these traits were significantly positively correlated with one another and with flowering time (FT). Based on a high-density SNP map, a total of 102 additive quantitative trait loci (QTL) were identified across six environments. Seventy-two consensus QTL and 49 unique QTL were identified using a two-round strategy of QTL meta-analysis. Notably, a total of 19 major QTL, including 11 novel ones, were detected for these traits, which comprised two QTL clusters on chromosomes A02 and A07. Conditional QTL mapping was performed to preliminarily evaluate the genetic basis (pleiotropy or tight linkage) of the co-localized QTL. In addition, QTL by environment interactions (QEI) mapping was performed to verify the additive QTL and estimate the QEI effect. In the genomic regions of all major QTL, orthologs of the genes involved in phytohormone biosynthesis, phytohormone signaling, flower development, and cell differentiation in Arabidopsis were proposed as candidate genes. Of these, BnaA02g02560, an ortholog of Arabidopsis GASA4, was suggested as a candidate gene for PH, SD, and FT; and BnaA02g08490, an ortholog of Arabidopsis GNL, was associated with PH, BIH and FT. These results provide useful information for further genetic studies on stem-related traits and plant growth adaptation.

Dissecting the complex regulation of lodging resistance in Brassica napus

Lodging continues to be a major cause of yield loss in important crop species such as Brassica napus. Understanding the genetic regulation of lodging resistance is therefore of key interest to breeders worldwide. Current strategies aimed at minimising lodging risk involve the incorporation of dwarfing genes or the application of plant growth regulators. However, despite these efforts, lodging continues to be a persistent problem and it is therefore of high interest that novel, complimentary strategies for lodging control are implemented. One approach would be to focus on understanding the genetic properties underlying stem mechanical strength. With this in mind, we screened a training genetic diversity panel of B. napus accession for variation in stem mechanical strength and related traits. Using Associative Transcriptomics, we identified molecular markers for a suite of valuable traits. Using an independent test genetic diversity panel, we show that the methods employed are robust for identification of predictive markers. Furthermore, based on conserved synteny with Arabidopsis thaliana, we are able to provide a biological context to the marker associations detected and provide evidence for a role in pectin methylesterification in contributing to stem mechanical strength in Brassicaceae.

Earliness traits in rapeseed (Brassica napus): SNP loci and candidate genes identified by genome-wide association analysis

Life cycle timing is critical for yield and productivity of Brassica napus (rapeseed) cultivars grown in different environments. To facilitate breeding for earliness traits in rapeseed, SNP loci and underlying candidate genes associated with the timing of initial flowering, maturity and final flowering, as well as flowering period (FP) were investigated in two environments in a diversity panel comprising 300 B. napus inbred lines. Genome-wide association studies (GWAS) using 201,817 SNP markers previously developed from SLAF-seq (specific locus amplified fragment sequencing) revealed a total of 131 SNPs strongly linked (P < 4.96E-07) to the investigated traits. Of these 131 SNPs, 40 fell into confidence intervals or were physically adjacent to previously published flowering time QTL or SNPs. Phenotypic effect analysis detected 35 elite allelic variants for early maturing, and 90 for long FP. Candidate genes present in the same linkage disequilibrium blocks (r²>0.6) or in 100 kb regions around significant trait-associated SNPs were screened, revealing 57 B. napus genes (33 SNPs) orthologous to 39 Arabidopsis thaliana flowering time genes. These results support the practical and scientific value of novel large-scale SNP data generation in uncovering the genetic control of agronomic traits in B. napus, and also provide a theoretical basis for molecular marker-assisted selection of earliness breeding in rapeseed.

A genome-wide survey with different rapeseed ecotypes uncovers footprints of domestication and breeding

Rapeseed (Brassica napus L.) is an important oilseed crop. Despite a short period of domestication and breeding, rapeseed has formed three diverse ecotype groups, namely spring, winter, and semi-winter. However, the genetic changes among the three ecotype groups have remained largely unknown. To detect selective signals, a set of 327 accessions from a worldwide collection were genotyped using a Brassica array, producing 33 186 high-quality single nucleotide polymorphisms (SNPs). Linkage disequilibrium (LD) was unevenly distributed across the genome. A total of 705 (78.2%) weak LD regions were found in the A subgenome, whereas 445 (72.6%) strong LD regions were in the C subgenome. By calculating the nucleotide diversity and population differentiation indices, a total of 198 selective sweeps were identified across ecotype groups, spanning 5.91% (37.9 Mb) of the genome. Within these genome regions, a few known functional genes or loci were found to be in association with environmental adaptability and yield-related traits. In particular, all 12 SNPs detected in significant association with flowering time among accessions were in the selection regions between ecotype groups. These findings provide new insights into the structure of the B. napus genome and uncover the footprints of domestication and breeding.

Development of SSR markers and identification of major quantitative trait loci controlling shelling percentage in cultivated peanut (Arachis hypogaea L.)

A total of 204,439 SSR markers were developed in diploid genomes, and 25 QTLs for shelling percentage were identified in a RIL population across 4 years including five consistent QTLs. Cultivated peanut (Arachis hypogaea L.) is an important grain legume providing edible oil and protein for human nutrition. Genome sequences of its diploid ancestors, Arachis duranensis and A. ipaensis, were reported, but their SSRs have not been well exploited and utilized hitherto. Shelling percentage is an important economic trait and its improvement has been one of the major objectives in peanut breeding programs. In this study, the genome sequences of A. duranensis and A. ipaensis were used to develop SSR markers, and a mapping population (Yuanza 9102 × Xuzhou 68-4) with 195 recombinant inbred lines was used to map QTLs controlling shelling percentage. The numbers of newly developed SSR markers were 84,383 and 120,056 in the A. duranensis and A. ipaensis genomes, respectively. Genotyping of the mapping population was conducted with both newly developed and previously reported markers. QTL analysis using the phenotyping data generated in Wuhan across four consecutive years and genotyping data of 830 mapped loci identified 25 QTLs with 4.46-17.01% of phenotypic variance explained in the four environments. Meta-analysis revealed five consistent QTLs that could be detected in at least two environments. Notably, the consistent QTL cqSPA09 was detected in all four environments and explained 10.47-17.01% of the phenotypic variance. The segregation in the progeny of a residual heterozygous line confirmed that the cpSPA09 locus had additive effect in increasing shelling percentage. These consistent and major QTL regions provide opportunity not only for further gene discovery, but also for the development of functional markers for breeding.

High-density ddRAD linkage and yield-related QTL mapping delimits a chromosomal region responsible for oil content in rapeseed (Brassica napus L.)

Rapeseed (Brassica napus L.) is one of the most important oil crops almost all over the world. Seed-related traits, including oil content (OC), silique length (SL), seeds per silique (SS), and seed weight (SW), are primary targets for oil yield improvement. To dissect the genetic basis of these traits, 192 recombinant inbred lines (RILs) were derived from two parents with distinct oil content and silique length. High-density linkage map with a total length of 1610.4 cM were constructed using 1,329 double-digestion restriction site associated DNA (ddRAD) markers, 107 insertion/deletions (INDELs), and 90 well-distributed simple sequence repeats (SSRs) markers. A total of 37 consensus quantitative trait loci (QTLs) were detected for the four traits, with individual QTL explained 3.1-12.8% of the phenotypic variations. Interestingly, one OC consensus QTL (cqOCA10b) on chromosome A10 was consistently detected in all three environments, and explained 9.8% to 12.8% of the OC variation. The locus was further delimited into an approximately 614 kb genomic region, in which the flanking markers could be further evaluated for marker-assisted selection in rapeseed OC improvement and the candidate genes targeted for map-based cloning and genetic manipulation.

Post-polyploidisation morphotype diversification associates with gene copy number variation

Genetic models for polyploid crop adaptation provide important information relevant for future breeding prospects. A well-suited model is Brassica napus, a recent allopolyploid closely related to Arabidopsis thaliana. Flowering time is a major adaptation trait determining life cycle synchronization with the environment. Here we unravel natural genetic variation in B. napus flowering time regulators and investigate associations with evolutionary diversification into different life cycle morphotypes. Deep sequencing of 35 flowering regulators was performed in 280 diverse B. napus genotypes. High sequencing depth enabled high-quality calling of single-nucleotide polymorphisms (SNPs), insertion-deletions (InDels) and copy number variants (CNVs). By combining these data with genotyping data from the Brassica 60 K Illumina® Infinium SNP array, we performed a genome-wide marker distribution analysis across the 4 ecogeographical morphotypes. Twelve haplotypes, including Bna.FLC.A10, Bna.VIN3.A02 and the Bna.FT promoter on C02_random, were diagnostic for the diversification of winter and spring types. The subspecies split between oilseed/kale (B. napus ssp. napus) and swedes/rutabagas (B. napus ssp. napobrassica) was defined by 13 haplotypes, including genomic rearrangements encompassing copies of Bna.FLC, Bna.PHYA and Bna.GA3ox1. De novo variation in copies of important flowering-time genes in B. napus arose during allopolyploidisation, enabling sub-functionalisation that allowed different morphotypes to appropriately fine-tune their lifecycle.

Co-localization of major quantitative trait loci for pod size and weight to a 3.7 cM interval on chromosome A05 in cultivated peanut (Arachis hypogaea L.)

Background

Cultivated peanut (Arachis hypogaea L.), an important source of edible oil and protein, is widely grown in tropical and subtropical areas of the world. Genetic improvement of yield-related traits is essential for improving yield potential of new peanut varieties. Genomics-assisted breeding (GAB) can accelerate the process of genetic improvement but requires linked markers for the traits of interest. In this context, we developed a recombinant inbred line (RIL) mapping population (Yuanza 9102 × Xuzhou 68-4) with 195 individuals and used to map quantitative trait loci (QTLs) associated with three important pod features, namely pod length, pod width and hundred-pod weight.

Results

QTL analysis using the phenotyping data generated across four environments in two locations and genotyping data on 743 mapped loci identified 15 QTLs for pod length, 11 QTLs for pod width and 16 QTLs for hundred-pod weight. The phenotypic variation explained (PVE) ranged from 3.68 to 27.84%. Thirteen QTLs were consistently detected in at least two environments and three QTLs (qPLA05.7, qPLA09.3 and qHPWA05.6) were detected in all four environments indicating their consistent and stable expression. Three major QTLs, detected in at least three environments, were found to be co-localized to a 3.7 cM interval on chromosome A05, and they were qPLA05.7 for pod length (16.89–27.84% PVE), qPWA05.5 for pod width (13.73–14.12% PVE), and qHPWA05.6 for hundred-pod weight (13.75–26.82% PVE). This 3.7 cM linkage interval corresponds to ~2.47 Mb genomic region of the pseudomolecule A05 of A. duranensis, including 114 annotated genes related to catalytic activity and metabolic process.

Conclusions

This study identified three major consistent and stable QTLs for pod size and weight which were co-localized in a 3.7 cM interval on chromosome A05. These QTL regions not only offer further investigation for gene discovery and development of functional markers but also provide opportunity for deployment of these QTLs in GAB for improving yield in peanut.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3456-x) contains supplementary material, which is available to authorized users.

High Density Linkage Map Construction and QTL Detection for Three Silique-Related Traits in Orychophragmus violaceus Derived Brassica napus Population

Seeds per silique (SS), seed weight (SW), and silique length (SL) are important determinant traits of seed yield potential in rapeseed (Brassica napus L.), and are controlled by naturally occurring quantitative trait loci (QTLs). Mapping QTLs to narrow chromosomal regions provides an effective means of characterizing the genetic basis of these complex traits. Orychophragmus violaceus is a crucifer with long siliques, many SS, and heavy seeds. A novel B. napus introgression line with many SS was previously selected from multiple crosses (B. rapa ssp. chinesis × O. violaceus) × B. napus. In present study, a doubled haploid (DH) population with 167 lines was established from a cross between the introgression line and a line with far fewer SS, in order to detect QTLs for silique-related traits. By screening with a Brassica 60K single nucleotide polymorphism (SNP) array, a high-density linkage map consisting of 1,153 bins and spanning a cumulative length of 2,209.1 cM was constructed, using 12,602 high-quality polymorphic SNPs in the DH population. The average recombination bin densities of the A and C subgenomes were 1.7 and 2.4 cM, respectively. 45 QTLs were identified for the three traits in all, which explained 4.0-34.4% of the total phenotypic variation; 20 of them were integrated into three unique QTLs by meta-analysis. These unique QTLs revealed a significant positive correlation between SS and SL and a significant negative correlation between SW and SS, and were mapped onto the linkage groups A05, C08, and C09. A trait-by-trait meta-analysis revealed eight, four, and seven consensus QTLs for SS, SW, and SL, respectively, and five major QTLs (cqSS.A09b, cqSS.C09, cqSW.A05, cqSW.C09, and cqSL.C09) were identified. Five, three, and four QTLs for SS, SW, and SL, respectively, might be novel QTLs because of the existence of alien genetic loci for these traits in the alien introgression. Thirty-eight candidate genes underlying nine QTLs for silique-related traits were identified.

Genetic mapping of yield traits using RIL population derived from Fuchuan Dahuasheng and ICG6375 of peanut (Arachis hypogaea L.)

The genetic architecture determinants of yield traits in peanut (Arachis hypogaea L.) are poorly understood. In the present study, an effort was made to map quantitative trait loci (QTLs) for yield traits using recombinant inbred lines (RIL). A genetic linkage map was constructed containing 609 loci, covering a total of 1557.48 cM with an average distance of 2.56 cM between adjacent markers. The present map exhibited good collinearity with the physical map of diploid species of Arachis. Ninety-two repeatable QTLs were identified for 11 traits including height of main stem, total branching number, and nine pod- and seed-related traits. Of the 92 QTLs, 15 QTLs were expressed across three environments and 65 QTLs were newly identified. Twelve QTLs for the height of main stem and the pod- and seed-related traits explaining more than 10 % of phenotypic variation showed a great potential for marker-assisted selection in improving these traits. The trait-by-trait meta-analysis revealed 33 consensus QTLs. The consensus QTLs and other QTLs were further integrated into 29 pleiotropic unique QTLs with the confidence interval of 1.86 cM on average. The significant co-localization of QTLs was consistent with the significant phenotypic correlations among these traits. The complexity of the genetic architecture of yield traits was demonstrated. The present QTLs for pod- and seed-related traits could be the most fundamental genetic factors contributing to the yield traits in peanut. The results provide a good foundation for fine mapping, cloning and designing molecular breeding of favorable genes in peanut.

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

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

Genome-Wide Association Study Reveals Candidate Genes for Control of Plant Height, Branch Initiation Height and Branch Number in Rapeseed (Brassica napus L.)

Plant architecture is crucial for rapeseed yield and is determined by plant height (PH), branch initiation height (BIH), branch number (BN) and leaf and inflorescence morphology. In this study, we measured three major factors (PH, BIH, and BN) in a panel of 333 rapeseed accessions across 4 years. A genome-wide association study (GWAS) was performed via Q + K model and the panel was genotyped using the 60 k Brassica Infinium SNP array. We identified seven loci for PH, four for BIH, and five for BN. Subsequently, by determining linkage disequilibrium (LD) decay associated with 38 significant SNPs, we gained 31, 15, and 17 candidate genes for these traits, respectively. We also showed that PH is significantly correlated with BIH, while no other correlation was revealed. Notably, a GA signaling gene (BnRGA) and a flowering gene (BnFT) located on chromosome A02 were identified as the most likely candidate genes associated with PH regulation. Furthermore, a meristem initiation gene (BnLOF2) and a NAC domain transcriptional factor (BnCUC3) that may be associated with BN were identified on the chromosome A07. This study reveals novel insight into the genetic control of plant architecture and may facilitate marker-based breeding for rapeseed.

Development and deployment of a high-density linkage map identified quantitative trait loci for plant height in peanut (Arachis hypogaea L.)

Plant height is one of the most important architecture traits in crop plants. In peanut, the genetic basis of plant height remains ambiguous. In this context, we genotyped a recombinant inbred line (RIL) population with 140 individuals developed from a cross between two peanut varieties varying in plant height, Zhonghua 10 and ICG 12625. Genotyping data was generated for 1,175 SSR and 42 transposon polymorphic markers and a high-density genetic linkage map was constructed with 1,219 mapped loci covering total map length of 2,038.75 cM i.e., accounted for nearly 80% of the peanut genome. Quantitative trait locus (QTL) analysis using genotyping and phenotyping data for three environments identified 8 negative-effect QTLs and 10 positive-effect QTLs for plant height. Among these QTLs, 8 QTLs had a large contribution to plant height that explained ≥10% phenotypic variation. Two major-effect consensus QTLs namely cqPHA4a and cqPHA4b were identified with stable performance across three environments. Further, the allelic recombination of detected QTLs proved the existence of the phenomenon of transgressive segregation for plant height in the RIL population. Therefore, this study not only successfully reported a high-density genetic linkage map of peanut and identified genomic region controlling plant height but also opens opportunities for further gene discovery and molecular breeding for plant height in peanut.

Functional conservation and divergence of GmCHLI genes in polyploid soybean

Polyploidy is prevalent in nature. As the fate of duplicated genes becomes more complicated when the encoded proteins function as oligomers, functional investigations into duplicated oligomer-encoding genes in polyploid genomes will facilitate our understanding of how traits are expressed. In this study, we identified GmCHLI1, a gene encoding the I subunit of magnesium (Mg)-chelatase, which functions in hexamers as responsible for the semi-dominant etiolation phenotype in soybean. Four GmCHLI copies derived from two polyploidy events were identified in the soybean genome. Further investigation with regard to expression patterns indicated that these four copies have diverged into two pairs; mutation in the other copy of the pair that includes GmCHLI1 also resulted in a chlorophyll-deficient phenotype. Protein interaction assays showed that these four GmCHLIs can interact with each other, but stronger interactions were found with mutated subunits. The results indicate that, in polyploidy, deficiency in each copy of duplicated oligomer-encoding genes could result in a mutant phenotype due to hetero-oligomer formation, which is different from the model of allelic dosage or functional redundancy. In addition, we interestingly found an increase in isoflavonoids in the heterozygous etiolated plants, which might be useful for improving soybean seed quality.

Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed

Background

Nitrogen use efficiency is an important breeding trait that can be modified to improve the sustainability of many crop species used in agriculture. Rapeseed is a major oil crop with low nitrogen use efficiency, making its production highly dependent on nitrogen input. This complex trait is suspected to be sensitive to genotype × environment interactions, especially genotype × nitrogen interactions. Therefore, phenotyping diverse rapeseed populations under a dense network of trials is a powerful approach to study nitrogen use efficiency in this crop. The present study aimed to determine the quantitative trait loci (QTL) associated with yield in winter oilseed rape and to assess the stability of these regions under contrasting nitrogen conditions for the purpose of increasing nitrogen use efficiency.

Results

Genome-wide association studies and linkage analyses were performed on two diversity sets and two doubled-haploid populations. These populations were densely genotyped, and yield-related traits were scored in a multi-environment design including seven French locations, six growing seasons (2009 to 2014) and two nitrogen nutrition levels (optimal versus limited). Very few genotype × nitrogen interactions were detected, and a large proportion of the QTL were stable across nitrogen nutrition conditions. In contrast, strong genotype × trial interactions in which most of the QTL were specific to a single trial were found. To obtain further insight into the QTL × environment interactions, genetic analyses of ecovalence were performed to identify the genomic regions contributing to the genotype × nitrogen and genotype × trial interactions. Fifty-one critical genomic regions contributing to the additive genetic control of yield-associated traits were identified, and the structural organization of these regions in the genome was investigated.

Conclusions

Our results demonstrated that the effect of the trial was greater than the effect of nitrogen nutrition levels on seed yield-related traits under our experimental conditions. Nevertheless, critical genomic regions associated with yield that were stable across environments were identified in rapeseed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0432-z) contains supplementary material, which is available to authorized users.

Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis

Seed number per pod (SNPP) is one of the major yield components and breeding targets in rapeseed that shows great variation and is invaluable for genetic improvement. To elucidate the genetic architecture and uncover the mechanism of SNPP, we identified five quantitative trait loci (QTLs) using the BnaZNRIL population, which were integrated with those of previous studies by physical map to demonstrate a complex and relatively complete genetic architecture of SNPP. A major QTL, qSN.A6, was successfully fine-mapped from 1910 to 267 kb using near-isogenic line (NIL). In addition, qSN.A6 exhibited an antagonistic pleiotropy on seed weight (SW), which is caused by a physiological interaction in which SNPP acts "upstream" of SW. Because the negative effect of qSN.A6 on SW cannot fully counteract its positive effect on SNPP, it also enhanced the final yield (17.4%), indicating its great potential for utilization in breeding. The following genetic and cytological experiments further confirmed that the different rate of ovule abortion was responsible for the ~5 seed difference between Zhongshuang11 and NIL-qSN.A6. This systematic approach to dissecting the comprehensive genetic architecture of SNPP and characterizing the underlying mechanism has advanced the understanding of SNPP and will facilitate the development of high-yield cultivars.

Rapid Identification of Candidate Genes for Seed Weight Using the SLAF-Seq Method in Brassica napus

Seed weight is a critical and direct trait for oilseed crop seed yield. Understanding its genetic mechanism is of great importance for yield improvement in Brassica napus breeding. Two hundred and fifty doubled haploid lines derived by microspore culture were developed from a cross between a large-seed line G-42 and a small-seed line 7–9. According to the 1000-seed weight (TSW) data, the individual DNA of the heaviest 46 lines and the lightest 47 lines were respectively selected to establish two bulked DNA pools. A new high-throughput sequencing technology, Specific Locus Amplified Fragment Sequencing (SLAF-seq), was used to identify candidate genes of TSW in association analysis combined with bulked segregant analysis (BSA). A total of 1,933 high quality polymorphic SLAF markers were developed and 4 associated markers of TSW were procured. A hot region of ~0.58 Mb at nucleotides 25,401,885–25,985,931 on ChrA09 containing 91 candidate genes was identified as tightly associated with the TSW trait. From annotation information, four genes (GSBRNA2T00037136001, GSBRNA2T00037157001, GSBRNA2T00037129001 and GSBRNA2T00069389001) might be interesting candidate genes that are highly related to seed weight.

Quantitative trait locus analysis for pod- and kernel-related traits in the cultivated peanut (Arachis hypogaea L.)

BACKGROUND

The cultivated peanut (Arachis hypogaea L.) is an important oil and food crop in the world. Pod- and kernel-related traits are direct factors involved in determining the yield of the peanut. However, the genetic basis underlying pod- and kernel-related traits in the peanut remained largely unknown, which hampered the improvement of peanut through marker-assisted selection. To understand the genetic basis underlying pod- and kernel-related traits in the peanut and provide more useful information for marker-assisted breeding, we conducted quantitative trait locus (QTL) analysis for pod length and width and seed length and width by use of two F2:3 populations derived from cultivar Fuchuan Dahuasheng × ICG 6375 (FI population) and cultivar Xuhua 13 × cultivar Zhonghua 6 (XZ population) in this study.

RESULTS

Two genetic maps containing 347 and 228 polymorphic markers were constructed for FI and XZ populations respectively. In total, 39 QTLs explaining 1.25-26.11% of the phenotypic variations were detected in two populations. For the FI population, 26 QTLs were detected between the two environments, among which twelve were not mapped before. For the XZ population, thirteen QTLs were detected, among which eight were not reported before. One QTL for pod width was repeatedly mapped between the two populations.

CONCLUSION

The QTL analyses for pod length and width and seed length and width were conducted in this study using two mapping populations. Novel QTLs were identified, which included two for pod length, four for pod width, five for seed length and one for seed width in the FI population, and three for pod length, three for pod width and two for seed width in the XZ population. Our results will be helpful for improving pod- and seed-related traits in peanuts through marker-assisted selection.

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

In Brassica napus breeding, traits related to commercial success are of highest importance for plant breeders. However, such traits can only be assessed in an advanced developmental stage. Molecular markers genetically linked to such traits have the potential to accelerate the breeding process of B. napus by marker-assisted selection. Therefore, the objectives of this study were to identify (i) genome regions associated with the examined agronomic and seed quality traits, (ii) the interrelationship of population structure and the detected associations, and (iii) candidate genes for the revealed associations. The diversity set used in this study consisted of 405 B. napus inbred lines which were genotyped using a 6K single nucleotide polymorphism (SNP) array and phenotyped for agronomic and seed quality traits in field trials. In a genome-wide association study, we detected a total of 112 associations between SNPs and the seed quality traits as well as 46 SNP-trait associations for the agronomic traits with a P < 1.28e-05 (Bonferroni correction of α = 0.05) for the inbreds of the spring and winter trial. For the seed quality traits, a single SNP-sulfur concentration in seeds (SUL) association explained up to 67.3% of the phenotypic variance, whereas for the agronomic traits, a single SNP-blossom color (BLC) association explained up to 30.2% of the phenotypic variance. In a basic local alignment search tool (BLAST) search within a distance of 2.5 Mbp around these SNP-trait associations, 62 hits of potential candidate genes with a BLAST-score of ≥100 and a sequence identity of ≥70% to A. thaliana or B. rapa could be found for the agronomic SNP-trait associations and 187 hits of potential candidate genes for the seed quality SNP-trait associations.

Genetic variation and inheritance of phytosterol and oil content in a doubled haploid population derived from the winter oilseed rape Sansibar × Oase cross

Identification of QTL for phytosterol content, oil content, fatty acids content, protein content of defatted meal, and seed weight by multiple interval mapping in a Brassica napus DH population. Phytosterols are minor seed constituents in oilseed rape which have recently drawn wide-interest from the food and nutrition industry due to their health benefit in lowering LDL cholesterol in humans. To understand the genetic basis of phytosterol content and its relationship with other seed quality traits in oilseed rape, QTL mapping was performed in a segregating DH population derived from the cross of two winter oilseed rape varieties, Sansibar and Oase, termed SODH population. Both parental lines are of canola quality which differ in phytosterol and oil content in seed. A genetic map was constructed for SODH population based on a total of 1638 markers organized in 23 linkage groups and covering a map length of 2350 cM with a mean marker interval of 2.0 cM. The SODH population and the parental lines were cultivated at six environments in Europe and were phenotyped for phytosterol content, oil content, fatty acids content, protein content of the defatted meal, and seed weight. Multiple interval mapping identified between one and six QTL for nine phytosterol traits, between two and six QTL for four fatty acids, five QTL for oil content, four QTL for protein content of defatted meal, and three QTL for seed weight. Colocalizations of QTL for different traits were more frequently observed than individual isolated QTL. Major QTL (R (2) ≥ 25%) were all located in the A genome, and the possible candidate genes were investigated by physical localization of the QTL to the reference genome sequence of Brassica rapa.

A genome-wide association study of plant height and primary branch number in rapeseed (Brassica napus)

Crop plant architecture plays a highly important role in its agronomic performance. Plant height (PH) and primary branch number (PB) are two major factors that affect the plant architecture of rapeseed (Brassica napus). Previous studies have shown that these two traits are controlled by multiple quantitative trait loci (QTL); however, QTLs have not been delimited to regions less than 10cM. Genome-wide association study (GWAS) is a highly efficient approach for identifying genetic loci controlling traits at relatively high resolution. In this study, variations in PH and PB of a panel of 472 rapeseed accessions that had previously been analyzed by a 60k SNP array were investigated for three consecutive years and studied by GWAS. Eight QTLs on chromosome A03, A05, A07 and C07 were identified for PH, and five QTLs on A01, A03, A07 and C07 were identified for PB. Although most QTLs have been detected in previous studies based on linkage analyses, the two QTLs of PH on A05 and the QTL of PB on C07 were novel. In the genomic regions close to the GWAS peaks, orthologs of the genes involved in flower development, phytohormone biosynthesis, metabolism and signaling in Arabidopsis were identified.

Genome-Wide Identification of QTL for Seed Yield and Yield-Related Traits and Construction of a High-Density Consensus Map for QTL Comparison in Brassica napus

Seed yield (SY) is the most important trait in rapeseed, is determined by multiple seed yield-related traits (SYRTs) and is also easily subject to environmental influence. Many quantitative trait loci (QTLs) for SY and SYRTs have been reported in Brassica napus; however, no studies have focused on seven agronomic traits simultaneously affecting SY. Genome-wide QTL analysis for SY and seven SYRTs in eight environments was conducted in a doubled haploid population containing 348 lines. Totally, 18 and 208 QTLs for SY and SYRTs were observed, respectively, and then these QTLs were integrated into 144 consensus QTLs using a meta-analysis. Three major QTLs for SY were observed, including cqSY-C6-2 and cqSY-C6-3 that were expressed stably in winter cultivation area for 3 years and cqSY-A2-2 only expressed in spring rapeseed area. Trait-by-trait meta-analysis revealed that the 144 consensus QTLs were integrated into 72 pleiotropic unique QTLs. Among them, all the unique QTLs affected SY, except for uq.A6-1, including uq.A2-3, uq.C1-2, uq.C1-3, uq.C6-1, uq.C6-5, and uq.C6-6 could also affect more than two SYRTs. According to the constructed high-density consensus map and QTL comparison from literatures, 36 QTLs from five populations were co-localized with QTLs identified in this study. In addition, 13 orthologous genes were observed, including five each gene for SY and thousand seed weight, and one gene each for biomass yield, branch height, and plant height. The genomic information of these QTLs will be valuable in hybrid cultivar breeding and in analyzing QTL expression in different environments.

Diverse regulatory factors associate with flowering time and yield responses in winter-type Brassica napus

BACKGROUND

Flowering time, plant height and seed yield are strongly influenced by climatic and day-length adaptation in crop plants. To investigate these traits under highly diverse field conditions in the important oilseed crop Brassica napus, we performed a genome-wide association study using data from diverse agroecological environments spanning three continents.

METHODS

A total of 158 European winter-type B.napus inbred lines were genotyped with 21,623 unique, single-locus single-nucleotide polymorphism (SNP) markers using the Brassica 60 K-SNP Illumina® Infinium consortium array. Phenotypic associations were calculated in the panel over the years 2010-2012 for flowering time, plant height and seed yield in 5 highly diverse locations in Germany, China and Chile, adding up to 11 diverse environments in total.

RESULTS

We identified 101 genome regions associating with the onset of flowering, 69 with plant height, 36 with seed yield and 68 cross-trait regions with potential adaptive value. Within these regions, B.napus orthologs for a number of candidate adaptation genes were detected, including central circadian clock components like CIRCADIAN CLOCK- ASSOCIATED 1 (Bna.CCA1) and the important flowering-time regulators FLOWERING LOCUS T (Bna.FT) and FRUITFUL (Bna.FUL).

DISCUSSION

Gene ontology (GO) enrichment analysis of candidate regions suggested that selection of genes involved in post-transcriptional and epigenetic regulation of flowering time may play a potential role in adaptation of B. napus to highly divergent environments. The classical flowering time regulators Bna.FLC and Bna.CO were not found among the candidate regions, although both show functional variation. Allelic effects were additive for plant height and yield, but not for flowering time. The scarcity of positive minor alleles for yield in this breeding pool points to a lack of diversity for adaptation that could restrict yield gain in the face of environmental change.

CONCLUSIONS

Our study provides a valuable framework to further improve the adaptability and yield stability of this recent allopolyploid crop under changing environments. The results suggest that flowering time regulation within an adapted B. napus breeding pool is driven by a high number of small modulating processes rather than major transcription factors like Bna.CO. In contrast, yield regulation appears highly parallel, therefore yield could be increased by pyramiding positively associated haplotypes.