Unraveling the complex trait of crop yield with quantitative trait loci mapping in Brassica napus

Multi-trait multi-environment genomic prediction of preliminary yield trial in pulse crop

Phenotypic selection of complex traits such as seed yield and protein in the preliminary yield trial (PYT) is often constrained by limited seed availability, resulting in trials with few environments and minimal to no replications. Multi-trait multi-environment enabled genomic prediction (MTME-GP) offers a valuable alternative to predict missing phenotypes of selection candidates for multiple traits and diverse environments. In this study, we assessed the efficiency of MTME-GP for improving seed protein and seed yield in field pea, the top two breeding targets but highly antagonistic traits in pulse crop. We utilized a set of 300 selection candidates in the PYT that virtually represented all possible families of the North Dakota State University field pea breeding program. Selection candidates were evaluated in three diverse, contrasting environments, as indicated by a range of heritability. Using whole- and split-environment cross validation schemes, MTME-GP had higher predictive ability than a standard additive G-BLUP model. Integrating a range of overlapping genotypes in between environments showed improvement on the predictive ability of the MTME-GP model but tends to plateau at 50%-80% training set size. Regardless of the cross-validation scheme, accuracy was among the lowest in stressed environments, presumably due to low heritability for seed protein and yield. This study provided insights into the potential of MTME-GP in a public pulse crop breeding program. The MTME-GP framework can be further improved with more testing environments and integration of additional orthogonal information in the early stages of the breeding pipeline.

Exploring silique number in Brassica napus L.: Genetic and molecular advances for improving yield

Silique number is a crucial yield-related trait for the genetic enhancement of rapeseed (Brassica napus L.). The intricate molecular process governing the regulation of silique number involves various factors. Despite advancements in understanding the mechanisms regulating silique number in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the molecular processes involved in controlling silique number in rapeseed remain largely unexplored. In this review, we identify candidate genes and review the roles of genes and environmental factors in regulating rapeseed silique number. We use genetic regulatory networks for silique number in Arabidopsis and grain number in rice to uncover possible regulatory pathways and molecular mechanisms involved in regulating genes associated with rapeseed silique number. A better understanding of the genetic network regulating silique number in rapeseed will provide a theoretical basis for the genetic improvement of this trait and genetic resources for the molecular breeding of high-yielding rapeseed.

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.

Molecular mapping and candidate gene identification of two major quantitative trait loci associated with silique length in oilseed rape (Brassica napus L.)

Rapeseed is a significant global source of plant oil. Silique size, particularly silique length (SL), impacts rapeseed yield. SL is a typical quantitative trait controlled by multiple genes. In our previous study, we constructed a DH population of 178 families known as the 158A-SGDH population. In this study, through SL QTL mapping, we identified twenty-six QTL for SL across five replicates in two environments. A QTL meta-analysis revealed eight consensus QTL, including two major QTL: cqSL.A02-1 (11.32-16.44% of PVE for SL), and cqSL.C06-1 (10.90-11.95% of PVE for SL). Based on biparental resequencing data and microcollinearity analysis of target regions in Brassica napus and Arabidopsis, we identified 11 candidate genes at cqSL.A02-1 and 6 candidate genes at cqSL.C06-1, which are potentially associated with silique development. Furthermore, transcriptome analysis of silique valves from both parents on the 14th, 21st, and 28th days after pollination (DAP) combined with gene function annotation revealed three significantly differentially expressed genes at cqSL.A02-1, BnaA02G0058500ZS, BnaA02G0060100ZS, and BnaA02G0060900ZS. Only the gene BnaC06G0283800ZS showed significant differences in parental transcription at cqSL.C06-1. Two tightly linked insertion-deletion markers for the cqSL.A02-1 and cqSL.C06-1 loci were developed. Using these two QTL, we generated four combinations: A02SGDH284C06158A, A02SGDH284C06SGDH284, A02158AC06158A, and A02158AC06SGDH284. Subsequent analysis identified an ideal QTL combination, A02158AC06SGDH284, which exhibited the longest SL of this type, reaching 6.06 ± 0.10 cm, significantly surpassing the other three combinations. The results will provide the basis for the cloning of SL-related genes of rapeseed, along with the development of functional markers of target genes and the breeding of rapeseed varieties.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01464-x.

Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

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.

BnPLP1 Positively Regulates Flowering Time, Plant Height, and Main Inflorescence Length in Brassica napus

Protein prenylation mediated by the Arabidopsis thaliana PLURIPETALA (AtPLP) gene plays a crucial role in plant growth, development, and environmental response by adding a 15-carbon farnesyl group or one to two 20-carbon geranylgeranyl groups onto one to two cysteine residues at the C-terminus of the target protein. However, the homologous genes and their functions of AtPLP in rapeseed are unclear. In this study, bioinformatics analysis and gene cloning demonstrated the existence of two homologous genes of AtPLP in the Brassica napus L. genome, namely, BnPLP1 and BnPLP2. Evolutionary analysis revealed that BnPLP1 originated from the B. rapa L. genome, while BnPLP2 originated from the B. oleracea L. genome. Genetic transformation analysis revealed that the overexpression of BnPLP1 in Arabidopsis plants exhibited earlier flowering initiation, a prolonged flowering period, increased plant height, and longer main inflorescence length compared to the wild type. Contrarily, the downregulation of BnPLP1 expression in B. napus plants led to delayed flowering initiation, shortened flowering period, decreased plant height, and reduced main inflorescence length compared to the wild type. These findings indicate that the BnPLP1 gene positively regulates flowering time, plant height, and main inflorescence length. This provides a new gene for the genetic improvement of flowering time and plant architecture in rapeseed.

Transcriptomic Profiling of Shoot Apical Meristem Aberrations in the Multi-Main-Stem Mutant (ms) of Brassica napus L

Rapeseed (Brassica napus L.) is a globally important oilseed crop with various uses, including the consumption of its succulent stems as a seasonal vegetable, but its uniaxial branching habit limits the stem yield. Therefore, developing a multi-stem rapeseed variety has become increasingly crucial. In this study, a natural mutant of the wild type (ZY511, Zhongyou511) with stable inheritance of the multi-stem trait (ms) was obtained, and it showed abnormal shoot apical meristem (SAM) development and an increased main stem number compared to the WT. Histological and scanning electron microscopy analyses revealed multiple SAMs in the ms mutant, whereas only a single SAM was found in the WT. Transcriptome analyses showed significant alterations in the expression of genes involved in cytokinin (CK) biosynthesis and metabolism pathways in the ms mutant. These findings provide insight into the mechanism of multi-main-stem formation in Brassica napus L. and lay a theoretical foundation for breeding multi-main-stem rapeseed vegetable varieties.

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.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r ² ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.

Genetic analysis and QTL mapping for silique density in rapeseed (Brassica napus L.)

KEY MESSAGE

Genetic models, QTLs and candidate gene for silique density on main inflorescence of rapeseed were identified. Silique density is one of the critical factors to determine seed yield and plant architecture in rapeseed (Brassica napus L.); however, the genetic control of this trait is largely unknown. In this study, the genetic model for silique density on main inflorescence (SDMI) of rapeseed was estimated according to the phenotypic data of P₁ (an inbreed line with high SDMI), P₂ (an inbreed line with low SDMI), F₁, F₂, BC₁P₁ and BC₁P₂ populations, revealing that SDMI is probably controlled by multi-minor genes with or without major gene. The QTLs for SDMI and its component characters including silique number on main inflorescence (SNMI) and main inflorescence length (MIL) were consequently mapped from a DH population derived from P₁ and P₂ by using a genetic linkage map constructed by restriction site-associated DNA sequencing (RAD seq) technology. A total of eight, 14 and three QTLs were identified for SDMI, SNMI and MIL under three environments, respectively, with an overlap among SDMI and SNMI in 55.7-75.4 cm on linkage group C06 which corresponding to 11.6-27.3 Mb on chromosome C06. Genomic resequencing was further conducted between a high- and a low-SDMI pool constructed from the DH population, and QTL-seq analysis identified a 0.15 Mb interval (25.98-26.13 Mb) from the C06-QTL region aforementioned. Transcriptome sequencing and qRT-PCR identified one possible candidate gene (BnARGOS) from the 0.15 Mb interval. This study will provide novel insights into the genetic basis of SD in rapeseed.

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.

Molecular mapping and characterization of QTLs for grain quality traits in a RIL population of US rice under high nighttime temperature stress

Elevated nighttime temperatures resulting from climate change significantly impact the rice crop worldwide. The rice (Oryza sativa L.) plant is highly sensitive to high nighttime temperature (HNT) during grain-filling (reproductive stage). HNT stress negatively affects grain quality traits and has a major impact on the value of the harvested rice crop. In addition, along with grain dimensions determining rice grain market classes, the grain appearance and quality traits determine the rice grain market value. During the last few years, there has been a major concern for rice growers and the rice industry over the prevalence of rice grains opacity and the reduction of grain dimensions affected by HNT stress. Hence, the improvement of heat-stress tolerance to maintain grain quality of the rice crop under HNT stress will bolster future rice value in the market. In this study, 185 F₁₂-recombinant inbred lines (RILs) derived from two US rice cultivars, Cypress (HNT-tolerant) and LaGrue (HNT-sensitive) were screened for the grain quality traits grain length (GL), grain width (GW), and percent chalkiness (%chalk) under control and HNT stress conditions and evaluated to identify the genomic regions associated with the grain quality traits. In total, there were 15 QTLs identified; 6 QTLs represented under control condition explaining 3.33% to 8.27% of the phenotypic variation, with additive effects ranging from - 0.99 to 0.0267 on six chromosomes and 9 QTLs represented under HNT stress elucidating 6.39 to 51.53% of the phenotypic variation, with additive effects ranging from - 8.8 to 0.028 on nine chromosomes for GL, GW, and % chalk. These 15 QTLs were further characterized and scanned for natural genetic variation in a japonica diversity panel (JDP) to identify candidate genes for GL, GW, and %chalk. We found 6160 high impact single nucleotide polymorphisms (SNPs) characterized as such depending on their type, region, functional class, position, and proximity to the gene and/or gene features, and 149 differentially expressed genes (DEGs) in the 51 Mbp genomic region comprising of the 15 QTLs. Out of which, 11 potential candidate genes showed high impact SNP associations. Therefore, the analysis of the mapped QTLs and their genetic dissection in the US grown Japonica rice genotypes at genomic and transcriptomic levels provide deep insights into genetic variation beneficial to rice breeders and geneticists for understanding the mechanisms related to grain quality under heat stress in rice.

Phenotyping of Silique Morphology in Oilseed Rape Using Skeletonization with Hierarchical Segmentation

Silique morphology is an important trait that determines the yield output of oilseed rape (Brassica napus L.). Segmenting siliques and quantifying traits are challenging because of the complicated structure of an oilseed rape plant at the reproductive stage. This study aims to develop an accurate method in which a skeletonization algorithm was combined with the hierarchical segmentation (SHS) algorithm to separate siliques from the whole plant using 3-dimensional (3D) point clouds. We combined the L1-median skeleton with the random sample consensus for iteratively extracting skeleton points and optimized the skeleton based on information such as distance, angle, and direction from neighborhood points. Density-based spatial clustering of applications with noise and weighted unidirectional graph were used to achieve hierarchical segmentation of siliques. Using the SHS, we quantified the silique number (SN), silique length (SL), and silique volume (SV) automatically based on the geometric rules. The proposed method was tested with the oilseed rape plants at the mature stage grown in a greenhouse and field. We found that our method showed good performance in silique segmentation and phenotypic extraction with R ² values of 0.922 and 0.934 for SN and total SL, respectively. Additionally, SN, total SL, and total SV had the statistical significance of correlations with the yield of a plant, with R values of 0.935, 0.916, and 0.897, respectively. Overall, the SHS algorithm is accurate, efficient, and robust for the segmentation of siliques and extraction of silique morphological parameters, which is promising for high-throughput silique phenotyping in oilseed rape breeding.

Linkage and association mapping of ovule number per ovary (ON) in oilseed rape (Brassica napus L.)

Ovule number (ON) produced during flower development determines the maximum number of seeds per silique and thereby affects crop productivity; however, the genetic basis of ON remains poorly understood in oilseed rape (Brassica napus). In this study, we genetically dissected the ON variations in a double haploid (DH) population and in natural population (NP) by linkage mapping and genome-wide association analysis. Phenotypic analysis showed that ON displayed normal distribution in both populations with the broad-sense heritability of 0.861 (DH population) and 0.930 (natural population). Linkage mapping identified 5 QTLs related to ON, including qON-A03, qON-A07, qON-A07-2, qON-A10, and qON-C06. Genome-wide association studies (GWAS) revealed 214, 48, and 40 significant single-nucleotide polymorphisms (SNPs) by individually using the single-locus model GLM and the multiple-locus model MrMLM and FASTMrMLM. The phenotypic variation explained (PVE) by these QTLs and SNPs ranged from 2.00-17.40% to 5.03-7.33%, respectively. Integration of the results from both strategies identified four consensus genomic regions associated with ON from the chromosomes A03, A07, and A10. Our results preliminarily resolved the genetic basis of ON and provides useful molecular markers for plant yield improvement in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01355-7.

Integrated multi-locus genome-wide association studies and transcriptome analysis for seed yield and yield-related traits in Brassica napus

Rapeseed (Brassica napus L.), the third largest oil crop, is an important source of vegetable oil and biofuel for the world. Although the breeding and yield has been improved, rapeseed still has the lowest yield compared with other major crops. Thus, increasing rapeseed yield is essential for the high demand of vegetable oil and high-quality protein for live stocks. Silique number per plant (SN), seed per pod (SP), and 1000-seed weight (SW) are the three important factors for seed yield in rapeseed. Some yield-related traits, including plant height (PH), flowering time (FT), primary branch number (BN) and silique number per inflorescence (SI) also affect the yield per plant (YP). Using six multi-locus genome-wide association study (ML-GWAS) approaches, a total of 908 yield-related quantitative trait nucleotides (QTNs) were identified in a panel consisting of 403 rapeseed core accessions based on whole-genome sequencing. Integration of ML-GWAS with transcriptome analysis, 79 candidate genes, including BnaA09g39790D (RNA helicase), BnaA09g39950D (Lipase) and BnaC09g25980D (SWEET7), were further identified and twelve genes were validated by qRT-PCRs to affect the SW or SP in rapeseed. The distribution of superior alleles from nineteen stable QTNs in 20 elite rapeseed accessions suggested that the high-yielding accessions contained more superior alleles. These results would contribute to a further understanding of the genetic basis of yield-related traits and could be used for crop improvement in B. napus.

QTL analysis and candidate gene prediction for seed density per silique by QTL-seq and RNA-seq in spring Brassica napus L

Seed density per silique (SD) is an important agricultural trait and plays an important role in the yield performance of Brassica napus L. (B. napus). In this study, a genetic linkage map was constructed using a double haploid (DH) population with 213 lines derived from a cross between a low SD line No. 935 and a high SD line No. 3641, and a total of 1,098,259 SNP (single-nucleotide polymorphisms) markers and 2,102 bins were mapped to 19 linkage groups. Twenty-eight QTLs for SD were detected on chromosomes A02, A04, A05, A09, C02, C03, C06, and C09 of B. napus, of which eight QTLs were on chromosome A09 and explained 5.89%-13.24% of the phenotypic variation. Furthermore, a consistent QTL for SD on chromosome A09, cqSD-A9a, was identified in four environments by QTL meta-analysis, explaining 10.68% of the phenotypic variation. In addition, four pairs of epistatic interactions were detected in the DH population via QTL epistasis analysis, indicating that SD is controlled not only by additive effects but also by epistatic effects that play an important role in spring B. napus., but with little environmental effect. Moreover, 18 closely linked SSR markers for cqSD-A9a were developed, as a result, it was mapped to a 1.86Mb (7.80-9.66 Mb) region on chromosome A09. A total of 13 differentially expressed genes (DEGs) were screened in the candidate interval by RNA-seq analysis, which were differentially expressed in buds, leaves and siliques both between and siliques both between two parents and two pools of extremely high-SD and low-SD lines in the DH population. Three of 13 DEGs were possible candidate genes that might control SD: BnaA09g14070D, which encodes a callose synthase that plays an important role in development and stress responses; BnaA09g14800D, a plant synaptic protein that encodes a membrane component; and BnaA09g18250D, which is responsible for DNA binding, transcriptional regulation, and sequence-specific DNA binding and is involved in the response to growth hormone stimulation. Overall, these results lay a foundation for fine mapping and gene cloning for SD in B. napus.

Genetic structure and molecular mechanism underlying the stalk lodging traits in maize (Zea mays L.)

Stalk lodging seriously affects yield and quality of crops, and it can be caused by several factors, such as environments, developmental stages, and internal chemical components of plant stalks. Breeding of stalk lodging-resistant varieties is thus an important task for maize breeders. To better understand the genetic basis underlying stalk lodging resistance, several methods such as quantitative trait locus (QTL) mapping and genome-wide association study (GWAS) have been used to mine potential gene resources. Based on different types of genetic populations and mapping methods, many significant loci associated with stalk lodging resistance have been identified so far. However, few work has been performed to compare and integrate these reported genetic loci. In this study, we first collected hundreds of QTLs and quantitative trait nucleotides (QTNs) related to stalk lodging traits in maize. Then we mapped and integrated the QTLs and QTNs in maize genome to identify overlapped hotspot regions. Based on the genomic confidence intervals harboring these overlapped hotspot regions, we predicted candidate genes related to stalk lodging traits. Meanwhile, we mapped reported genes to these hotspot regions. Finally, we constructed molecular regulatory networks underlying stalk lodging resistance in maize. Collectively, this study provides not only useful genetic loci for deeply exploring molecular mechanisms of stalk lodging resistance traits, but also potential candidate genes and targeted strategies for improving stalk lodging resistance to increase crop yields in future.

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.

Systematic trait dissection in oilseed rape provides a comprehensive view, further insight, and exact roadmap for yield determination

Background

Yield is the most important and complex trait that is influenced by numerous relevant traits with very complicated interrelations. While there are a large number of studies on the phenotypic relationship and genetic basis of yield traits, systematic studies with further dissection focusing on yield are limited. Therefore, there is still lack of a comprehensive and in-depth understanding of the determination of yield.

Results

In this study, yield was systematically dissected at the phenotypic, genetic to molecular levels in oilseed rape (Brassica napus L.). The analysis of correlation, network, and principal component for 21 traits in BnaZN-RIL population showed that yield was determined by a complex trait network with key contributors. The analysis of the constructed high-density single nucleotide polymorphism (SNP) linkage map revealed the concentrated distribution of distorted and heterozygous markers, likely due to selection on genes controlling the growth period and yield heterosis. A total of 134 consensus quantitative trait loci (QTL) were identified for 21 traits, of which all were incorporated into an interconnecting QTL network with dozens of hub-QTL. Four representative hub-QTL were further dissected to the target or candidate genes that governed the causal relationships between the relevant traits.

Conclusions

The highly consistent results at the phenotypic, genetic, and molecular dissecting demonstrated that yield was determined by a multilayer composite network that involved numerous traits and genes showing complex up/down-stream and positive/negative regulation. This provides a systematic view, further insight, and exact roadmap for yield determination, which represents a significant advance toward the understanding and dissection of complex traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13068-022-02134-w.

Salt tolerance QTLs of an endemic rice landrace, Horkuch at seedling and reproductive stages

Salinity has a significant negative impact on production of rice. To cope with the increased soil salinity due to climate change, we need to develop salt tolerant rice varieties that can maintain their high yield. Rice landraces indigenous to coastal Bangladesh can be a great resource to study the genetic basis of salt adaptation. In this study, we implemented a QTL analysis framework with a reciprocal mapping population developed from a salt tolerant landrace Horkuch and a high yielding rice variety IR29. Our aim was to detect genetic loci that contributes to the salt adaptive responses of the two different developmental stages of rice which are very sensitive to salinity stress. We identified 14 QTLs for 9 traits and found that most are unique to specific developmental stages. In addition, we detected a significant effect of the cytoplasmic genome on the QTL model for some traits such as leaf total potassium and filled grain weight. This underscores the importance of considering cytoplasm-nuclear interaction for breeding programs. Finally, we identified QTLs co-localization for multiple traits that highlights the possible constraint of multiple QTL selection for breeding programs due to different contributions of a donor allele for different traits.

Genome-wide association study identifies candidate genes and favorable haplotypes for seed yield in Brassica napus

Oilseed rape (Brassica napus L.) is one of the most essential oil crops. Genetic improvement of seed yield (SY) is a major aim of B. napus breeding. Several studies have been reported on the genetic mechanisms of SY of B. napus. Here, a genome-wide association study (GWAS) of SY was conducted using a panel of 403 natural accessions of B. napus, with more than five million high-quality single-nucleotide polymorphisms (SNPs). A total of 1773 significant SNPs were detected associated with SY, and 783 significant SNPs were co-located with previously reported QTLs. The lead SNPs chrA01__8920351 and chrA02__4555979 were jointly detected in Trial 2_2 and Trial 2_mean value, and in Trial 1_2 and Trial 1_mean value, respectively. Subsequently, two candidate genes of BnaA01g17200D and BnaA02g08680D were identified through combining transcriptome, candidate gene association analysis, and haplotype analysis. BnaA09g10430D detected through lead SNP chrA09__5160639 was associated with SY of B. napus. Our results provide valuable information for studying the genetic control of seed yield in B. napus and valuable genes, haplotypes, and cultivars resources for the breeding of high seed yield B. napus cultivars.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01332-6.

Genome-wide analysis of the JAZ subfamily of transcription factors and functional verification of BnC08.JAZ1-1 in Brassica napus

BACKGROUND

JAZ subfamily plays crucial roles in growth and development, stress, and hormone responses in various plant species. Despite its importance, the structural and functional analyses of the JAZ subfamily in Brassica napus are still limited.

RESULTS

Comparing to the existence of 12 JAZ genes (AtJAZ1-AtJAZ12) in Arabidopsis, there are 28, 31, and 56 JAZ orthologues in the reference genome of B. rapa, B. oleracea, and B. napus, respectively, in accordance with the proven triplication events during the evolution of Brassicaceae. The phylogenetic analysis showed that 127 JAZ proteins from A. thaliana, B. rapa, B. oleracea, and B. napus could fall into five groups. The structure analysis of all 127 JAZs showed that these proteins have the common motifs of TIFY and Jas, indicating their conservation in Brassicaceae species. In addition, the cis-element analysis showed that the main motif types are related to phytohormones, biotic and abiotic stresses. The qRT-PCR of the representative 11 JAZ genes in B. napus demonstrated that different groups of BnJAZ individuals have distinct patterns of expression under normal conditions or treatments with distinctive abiotic stresses and phytohormones. Especially, the expression of BnJAZ52 (BnC08.JAZ1-1) was significantly repressed by abscisic acid (ABA), gibberellin (GA), indoleacetic acid (IAA), polyethylene glycol (PEG), and NaCl treatments, while induced by methyl jasmonate (MeJA), cold and waterlogging. Expression pattern analysis showed that BnC08.JAZ1-1 was mainly expressed in the vascular bundle and young flower including petal, pistil, stamen, and developing ovule, but not in the stem, leaf, and mature silique and seed. Subcellular localization showed that the protein was localized in the nucleus, in line with its orthologues in Arabidopsis. Overexpression of BnC08.JAZ1-1 in Arabidopsis resulted in enhanced seed weight, likely through regulating the expression of the downstream response genes involved in the ubiquitin-proteasome pathway and phospholipid metabolism pathway.

CONCLUSIONS

The systematic identification, phylogenetic, syntenic, and expression analyses of BnJAZs subfamily improve our understanding of their roles in responses to stress and phytohormone in B. napus. In addition, the preliminary functional validation of BnC08.JAZ1-1 in Arabidopsis demonstrated that this subfamily might also play a role in regulating seed weight.

Study on exogenous application of thidiazuron on seed size of Brassica napus L

Thidiazuron (TDZ) is a novel and efficient cytokinin commonly used in tissue culture, and numerous studies have demonstrated that TDZ can increase berry size. However, no study to date has explored the effect of TDZ on seed size of Brassica napus and the mechanism. To shed light on the effect of TDZ on the seed size of B. napus, four different concentrations of TDZ were applied to B. napus. Results indicated that TDZ treatment could increase the seed diameter and silique length of B. napus to varying degrees and 100 and 200 μmol/L TDZ treatments were the most effective with a 3.6 and 4.6% increase in seed diameter, respectively. In addition, the yield of B. napus was also substantially increased under TDZ treatment. On the other hand, confocal micrographs of embryos and cotyledon cells suggested that embryos and their cotyledon epidermal cells treated with 200 μmol/L TDZ were obviously larger in size than the control. Furthermore, TDZ promoted the upregulation of some key maternal tissue growth-related genes, including two G-protein signaling genes (AGG3 and RGA1) and two transcriptional regulators (ANT and GS2). The expression analysis of genes related to the auxin metabolic pathways, G-protein signaling, endosperm growth and transcriptional regulators confirmed that treatment with TDZ negatively regulated the key genes ABI5, AGB1, AP2, ARF2, and ARF18 during bud development stage and florescence. The results strongly suggested that TDZ might regulate the transcriptional levels of key genes involved in auxin metabolic pathways, G-protein signaling, endosperm growth and transcriptional regulators, which resulted in bigger cotyledon epidermal cells and seed size in B. napus. This study explored the mechanism of TDZ treatment on the seed size of B. napus and provided an important reference for improving rapeseed yield.

Genome-wide analysis of the serine carboxypeptidase-like (SCPL) proteins in Brassica napus L

The serine carboxypeptidase-like protein (SCPL) family plays a key part in plant growth, development and stress responses. However, the serine carboxypeptidase-like (SCPL) proteins in Brassica napus L. (B. napus) have not been reported yet. Here, we identified a total of 117 putative SCPL genes in B. napus, which were unevenly distributed on all 19 chromosomes and were divided into three groups (carboxypeptidase Ⅰ to Ⅲ) according to their phylogenetic relationships. Synteny and duplication analysis revealed that the SCPL gene family of B. napus was amplified during allopolyploidization, in which the whole genome triplication and dispersed duplication played critical roles. After the separation of Brassica and Arabidopsis lineages, orthologous gene analysis showed that many SCPL genes were lost during the evolutionary process in B. rapa, B. oleracea and B. napus. Subsequently, the analyses of the gene structure, conserved motifs, cis-element and expression patterns showed that the members in the same group were highly conserved. Furthermore, candidate gene based association study suggested the role of BnSCPL52 in controlling seed number per silique, seed weight and silique length and a CAPS marker was developed to distinguish different haplotypes. Our results provide an overview of rapeseed SCPL genes that enable us for further functional research and benefit the marker-assisted breeding in Brassica napus.

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.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

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.

QTL Identification for Stem Fiber, Strength and Rot Resistance in a DH Population from an Alien Introgression of Brassica napus

Stem fiber, stem strength and stem-rot resistance are important agronomic traits in Brassica napus. To understand the molecular mechanism that controls the stem-related traits, we investigated the stem lignin (ADL), cellulose (Cel), hemicellulose (Hem) content, S/G monolignol ratio (SG), stem breaking force (BF), breaking strength (F) and Sclerotinia sclerotiorum resistance (SSR). Each trait was significantly positively or negatively correlated with more than three of the other six traits. QTL mapping for ADL, Cel, Hem, SG, BF, F and SSR were performed using a doubled haploid population derived from an intertribal B. napus introgression line 'Y689' crossed with B. napus cv. 'Westar'. A total of 67 additive QTL were identified and integrated into 55 consensus QTL by meta-analysis. Among the 55 consensus QTL, 23 (41.8%) QTL were co-located and were integrated into 11 unique QTL. The QTL by environment (Q × E) interactions were analyzed and 22 combined QTL were identified. In addition, candidate genes within the QTL intervals were proposed based on the known function of Arabidopsis orthologs. These results provided valuable information for improving lodging resistance, S. sclerotiorum resistance and mechanized harvesting of B. napus.

BnaC7.ROT3, the causal gene of cqSL-C7, mediates silique length by affecting cell elongation in Brassica napus

Siliques are a major carbohydrate source of energy for later seed development in rapeseed (Brassica napus). Thus, silique length has received great attention from breeders. We previously detected a novel quantitative trait locus cqSL-C7 that controls silique length in B. napus. Here, we further validated the cqSL-C7 locus and isolated its causal gene (BnaC7.ROT3) by map-based cloning. In 'Zhongshuang11' (parent line with long siliques), BnaC7.ROT3 encodes the potential cytochrome P450 monooxygenase CYP90C1, whereas in 'G120' (parent line with short siliques), a single nucleotide deletion in the fifth exon of BnaC7.ROT3 results in a loss-of-function truncated protein. Sub-cellular localization and expression pattern analysis revealed that BnaC7.ROT3 is a membrane-localized protein mainly expressed in leaves, flowers and siliques. Cytological observations showed that the cells in silique walls of BnaC7.ROT3-transformed positive plants were longer than those of transgene-negative plants in the background of 'G120', suggesting that BnaC7.ROT3 affects cell elongation. Haplotype analysis demonstrated that most alleles of BnaC7.ROT3 are favorable in B. napus germplasms, and its homologs may also be involved in silique length regulation. Our findings provide novel insights into the regulatory mechanisms of natural silique length variations and valuable genetic resources for the improvement of silique length in rapeseed.

Epigenetics and its role in effecting agronomical traits

Climate-resilient crops with improved adaptation to the changing climate are urgently needed to feed the growing population. Hence, developing high-yielding crop varieties with better agronomic traits is one of the most critical issues in agricultural research. These are vital to enhancing yield as well as resistance to harsh conditions, both of which help farmers over time. The majority of agronomic traits are quantitative and are subject to intricate genetic control, thereby obstructing crop improvement. Plant epibreeding is the utilisation of epigenetic variation for crop development, and has a wide range of applications in the field of crop improvement. Epigenetics refers to changes in gene expression that are heritable and induced by methylation of DNA, post-translational modifications of histones or RNA interference rather than an alteration in the underlying sequence of DNA. The epigenetic modifications influence gene expression by changing the state of chromatin, which underpins plant growth and dictates phenotypic responsiveness for extrinsic and intrinsic inputs. Epigenetic modifications, in addition to DNA sequence variation, improve breeding by giving useful markers. Also, it takes epigenome diversity into account to predict plant performance and increase crop production. In this review, emphasis has been given for summarising the role of epigenetic changes in epibreeding for crop improvement.

Adaptation Strategies to Improve the Resistance of Oilseed Crops to Heat Stress Under a Changing Climate: An Overview

Temperature is one of the decisive environmental factors that is projected to increase by 1. 5°C over the next two decades due to climate change that may affect various agronomic characteristics, such as biomass production, phenology and physiology, and yield-contributing traits in oilseed crops. Oilseed crops such as soybean, sunflower, canola, peanut, cottonseed, coconut, palm oil, sesame, safflower, olive etc., are widely grown. Specific importance is the vulnerability of oil synthesis in these crops against the rise in climatic temperature, threatening the stability of yield and quality. The natural defense system in these crops cannot withstand the harmful impacts of heat stress, thus causing a considerable loss in seed and oil yield. Therefore, a proper understanding of underlying mechanisms of genotype-environment interactions that could affect oil synthesis pathways is a prime requirement in developing stable cultivars. Heat stress tolerance is a complex quantitative trait controlled by many genes and is challenging to study and characterize. However, heat tolerance studies to date have pointed to several sophisticated mechanisms to deal with the stress of high temperatures, including hormonal signaling pathways for sensing heat stimuli and acquiring tolerance to heat stress, maintaining membrane integrity, production of heat shock proteins (HSPs), removal of reactive oxygen species (ROS), assembly of antioxidants, accumulation of compatible solutes, modified gene expression to enable changes, intelligent agricultural technologies, and several other agronomic techniques for thriving and surviving. Manipulation of multiple genes responsible for thermo-tolerance and exploring their high expressions greatly impacts their potential application using CRISPR/Cas genome editing and OMICS technology. This review highlights the latest outcomes on the response and tolerance to heat stress at the cellular, organelle, and whole plant levels describing numerous approaches applied to enhance thermos-tolerance in oilseed crops. We are attempting to critically analyze the scattered existing approaches to temperature tolerance used in oilseeds as a whole, work toward extending studies into the field, and provide researchers and related parties with useful information to streamline their breeding programs so that they can seek new avenues and develop guidelines that will greatly enhance ongoing efforts to establish heat stress tolerance in oilseeds.

QTL Analysis of Five Silique-Related Traits in Brassica napus L. Across Multiple Environments

As an important physiological and reproductive organ, the silique is a determining factor of seed yield and a breeding target trait in rapeseed (Brassica napus L.). Genetic studies of silique-related traits are helpful for rapeseed marker-assisted high-yield breeding. In this study, a recombinant inbred population containing 189 lines was used to perform a quantitative trait loci (QTLs) analysis for five silique-related traits in seven different environments. As a result, 120 consensus QTLs related to five silique-related traits were identified, including 23 for silique length, 25 for silique breadth, 29 for silique thickness, 22 for seed number per silique and 21 for silique volume, which covered all the chromosomes, except C5. Among them, 13 consensus QTLs, one, five, two, four and one for silique length, silique breadth, silique thickness, seed number per silique and silique volume, respectively, were repeatedly detected in multiple environments and explained 4.38-13.0% of the phenotypic variation. On the basis of the functional annotations of Arabidopsis homologous genes and previously reported silique-related genes, 12 potential candidate genes underlying these 13 QTLs were screened and found to be stable in multiple environments by analyzing the re-sequencing results of the two parental lines. These findings provide new insights into the gene networks affecting silique-related traits at the QTL level in rapeseed.

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

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

Uncovering Trait Associations Resulting in Maximal Seed Yield in Winter and Spring Oilseed Rape

Seed yield is a complex trait for many crop species including oilseed rape (OSR) (Brassica napus), the second most important oilseed crop worldwide. Studies have focused on the contribution of distinct factors in seed yield such as environmental cues, agronomical practices, growth conditions, or specific phenotypic traits at the whole plant level, such as number of pods in a plant. However, how female reproductive traits contribute to whole plant level traits, and hence to seed yield, has been largely ignored. Here, we describe the combined contribution of 33 phenotypic traits within a B. napus diversity set population and their trade-offs at the whole plant and organ level, along with their interaction with plant level traits. Our results revealed that both Winter OSR (WOSR) and Spring OSR (SOSR); the two more economically important OSR groups in terms of oil production; share a common dominant reproductive strategy for seed yield. In this strategy, the main inflorescence is the principal source of seed yield, producing a good number of ovules, a large number of long pods with a concomitantly high number of seeds per pod. Moreover, we observed that WOSR opted for additional reproductive strategies than SOSR, presenting more plasticity to maximise seed yield. Overall, we conclude that OSR adopts a key strategy to ensure maximal seed yield and propose an ideal ideotype highlighting crucial phenotypic traits that could be potential targets for breeding.

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.

Fine mapping and candidate gene analysis of a major locus controlling ovule abortion and seed number per silique in Brassica napus L

A major QTL controlling ovule abortion and SN was fine-mapped to a 80.1-kb region on A8 in rapeseed, and BnaA08g07940D and BnaA08g07950D are the most likely candidate genes. The seed number per silique (SN), an important yield determining trait of rapeseed, is the final consequence of a complex developmental process including ovule initiation and the subsequent ovule/seed development. To explore the genetic mechanism regulating the natural variation of SN and its related components, quantitative trait locus (QTL) mapping was conducted using a doubled haploid (DH) population derived from the cross between C4-146 and C4-58B, which showed significant differences in SN and aborted ovule number (AON), but no obvious differences in ovule number (ON). QTL analysis identified 19 consensus QTLs for six SN-related traits across three environments. A novel QTL on chromosome A8, un.A8, which associates with multiple traits, except for ON, was stably detected across the three environments. This QTL explained more than 50% of the SN, AON and percentage of aborted ovules (PAO) variations as well as a moderate contribution on silique length (SL) and thousand seed weight (TSW). The C4-146 allele at the locus increases SN and SL but decreases AON, PAO and TSW. Further fine mapping narrowed down this locus into an 80.1-kb interval flanked by markers BM1668 and BM1672, and six predicted genes were annotated in the delimited region. Expression analyses and DNA sequencing showed that two homologs of Arabidopsis photosystem I subunit F (BnaA08g07940D) and zinc transporter 10 precursor (BnaA08g07950D) were the most promising candidate genes underlying this locus. These results provide a solid basis for cloning un.A8 to reduce the ovule abortion and increase SN in the yield improvement of rapeseed.

Genome-wide association analysis for yield-related traits at the R6 stage in a Chinese soybean mini core collection

BACKGROUND

Soybean (Glycine max (L.) Merr.) is an economically important crop for vegetable oil and protein production, and yield is a critical trait for grain/vegetable uses of soybean. However, our knowledge of the genes controlling the vegetable soybean yield remains limited.

OBJECTIVE

To better understand the genetic basis of the vegetable soybean yield.

METHODS

The 100-pod fresh weight (PFW), 100-seed fresh weight (SFW), kernel percent (KP) and moisture content of fresh seeds (MCFS) at the R6 stage are four yield-related traits for vegetable soybean. We investigated a soybean mini core collection composed of 224 germplasm accessions for four yield-related traits in two consecutive years. Based on 1514 single nucleotide polymorphisms (SNPs), genome-wide association studies (GWAS) were conducted using a mixed linear model (MLM).

RESULTS

Extensive phenotypic variation existed in the soybean mini core collection and significant positive correlations were shown among most of traits. A total of 16 SNP markers for PFW, SFW, KP and MCFS were detected in all environments via GWAS. Nine SNP markers were repeatedly identified in two environments. Among these markers, eight were located in or near regions where yield-related QTLs have been reported in previous studies, and one was a novel genetic locus identified in this study. In addition, we conducted candidate gene analysis to the large-effect SNP markers, a total of twelve genes were proposed as potential candidate genes of soybean yield at the R6 stage.

CONCLUSION

These results will be beneficial for understanding the genetic basis of soybean yield at the R6 stage and facilitating the pyramiding of favourable alleles for future high-yield breeding by marker-assisted selection in vegetable soybean.

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.

Construction of a high density genetic map for hexaploid kiwifruit (Actinidia chinensis var. deliciosa) using genotyping by sequencing

Commercially grown kiwifruit (genus Actinidia) are generally of two sub-species which have a base haploid genome of 29 chromosomes. The yellow-fleshed Actinidia chinensis var. chinensis, is either diploid (2n = 2x = 58) or tetraploid (2n = 4x = 116) and the green-fleshed cultivar A. chinensis var. deliciosa “Hayward,” is hexaploid (2n = 6x = 174). Advances in breeding green kiwifruit could be greatly sped up by the use of molecular resources for more efficient and faster selection, for example using marker-assisted selection (MAS). The key genetic marker that has been implemented for MAS in hexaploid kiwifruit is for gender testing. The limited marker-trait association has been reported for other polyploid kiwifruit for fruit and production traits. We have constructed a high-density linkage map for hexaploid green kiwifruit using genotyping-by-sequence (GBS). The linkage map obtained consists of 3686 and 3940 markers organized in 183 and 176 linkage groups for the female and male parents, respectively. Both parental linkage maps are co-linear with the A. chinensis “Red5” reference genome of kiwifruit. The linkage map was then used for quantitative trait locus (QTL) mapping, and successfully identified QTLs for king flower number, fruit number and weight, dry matter accumulation, and storage firmness. These are the first QTLs to be reported and discovered for complex traits in hexaploid kiwifruit.

A tetrasomic inheritance model and likelihood-based method for mapping quantitative trait loci in autotetraploid species

Dissecting the genetic architecture of quantitative traits in autotetraploid species is a methodologically challenging task, but a pivotally important goal for breeding globally important food crops, including potato and blueberry, and ornamental species such as rose. Mapping quantitative trait loci (QTLs) is now a routine practice in diploid species but is far less advanced in autotetraploids, largely due to a lack of analytical methods that account for the complexities of tetrasomic inheritance. We present a novel likelihood-based method for QTL mapping in outbred segregating populations of autotetraploid species. The method accounts properly for sophisticated features of gene segregation and recombination in an autotetraploid meiosis. It may model and analyse molecular marker data with or without allele dosage information, such as that from microarray or sequencing experiments. The method developed outperforms existing bivalent-based methods, which may fail to model and analyse the full spectrum of experimental data, in the statistical power of QTL detection, and accuracy of QTL location, as demonstrated by an intensive simulation study and analysis of data sets collected from a segregating population of potato (Solanum tuberosum). The study enables QTL mapping analysis to be conducted in autotetraploid species under a rigorous tetrasomic inheritance model.

Targeted mutagenesis of EOD3 gene in Brassica napus L. regulates seed production

Seed size and number are central to the evolutionary fitness of plants and are also crucial for seed production of crops. However, the molecular mechanisms of seed production control are poorly understood in Brassica crops. Here, we report the gene cloning, expression analysis, and functional characterization of the EOD3/CYP78A6 gene in rapeseed. BnaEOD3 has four copies located in two subgenomes, which exhibited a steady higher expression during seed development with differential expression among copies. The targeted mutations of BnaEOD3 gene were efficiently generated by stable transformation of the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeat) vector. These mutations were stably transmitted to T1 and T2 generations and a large collection of homozygous mutants with combined loss-of-function alleles across four BnaEOD3 copies were created for phenotyping. All mutant T1 lines had shorter siliques, smaller seeds, and an increased number of seeds per silique, in which the quadrable mutants showed the most significant changes in these traits. Consequently, the seed weight per plant in the quadrable mutants increased by 13.9% on average compared with that of wild type, indicating that these BnaEOD3 copies have redundant functions in seed development in rapeseed. The phenotypes of the different allelic combinations of BnaEOD3 copies also revealed gene functional differentiation among the two subgenomes. Cytological observations indicated that the BnaEOD3 could act maternally to promote cotyledon cell expansion and proliferation to regulate seed growth in rapeseed. Collectively, our findings reveal the quantitative involvement of the different BnaEOD3 copies function in seed development, but also provided valuable resources for rapeseed breeding programs.

Genome wide association analyses to understand genetic basis of flowering and plant height under three levels of nitrogen application in Brassica juncea (L.) Czern & Coss

Timely transition to flowering, maturity and plant height are important for agronomic adaptation and productivity of Indian mustard (B. juncea), which is a major edible oilseed crop of low input ecologies in Indian subcontinent. Breeding manipulation for these traits is difficult because of the involvement of multiple interacting genetic and environmental factors. Here, we report a genetic analysis of these traits using a population comprising 92 diverse genotypes of mustard. These genotypes were evaluated under deficient (N75), normal (N100) or excess (N125) conditions of nitrogen (N) application. Lower N availability induced early flowering and maturity in most genotypes, while high N conditions delayed both. A genotyping-by-sequencing approach helped to identify 406,888 SNP markers and undertake genome wide association studies (GWAS). 282 significant marker-trait associations (MTA's) were identified. We detected strong interactions between GWAS loci and nitrogen levels. Though some trait associated SNPs were detected repeatedly across fertility gradients, majority were identified under deficient or normal levels of N applications. Annotation of the genomic region (s) within ± 50 kb of the peak SNPs facilitated prediction of 30 candidate genes belonging to light perception, circadian, floral meristem identity, flowering regulation, gibberellic acid pathways and plant development. These included over one copy each of AGL24, AP1, FVE, FRI, GID1A and GNC. FLC and CO were predicted on chromosomes A02 and B08 respectively. CDF1, CO, FLC, AGL24, GNC and FAF2 appeared to influence the variation for plant height. Our findings may help in improving phenotypic plasticity of mustard across fertility gradients through marker-assisted breeding strategies.

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.

Identification of loci and candidate gene GmSPX-RING1 responsible for phosphorus efficiency in soybean via genome-wide association analysis

BACKGROUND

Phosphorus (P) is an essential element in maintaining high biomass and yield in crops. Soybean [Glycine max (L.) Merr.] requires a large amount of P during growth and development. Improvement of P efficiency and identification of P efficiency genes are important strategies for increasing soybean yield.

RESULTS

Genome-wide association analysis (GWAS) with NJAU 355 K SoySNP array was performed to identify single nucleotide polymorphisms (SNPs) significantly associated with three shoot P efficiency-related traits of a natural population of 211 cultivated soybeans and relative values of these traits under normal P (+P) condition and P deficiency (-P) condition. A total of 155 SNPs were identified significantly associated with P efficiency-related traits. SNPs that were significantly associated with shoot dry weight formed a SNP cluster on chromosome 11, while SNPs that were significantly associated with shoot P concentration formed a SNP cluster on chromosome 10. Thirteen haplotypes were identified based on 12 SNPs, and Hap9 was considered as the optimal haplotype. Four SNPs (AX-93636685, AX-93636692, AX-93932863, and AX-93932874) located on chromosome 10 were identified to be significantly associated with shoot P concentration under +P condition in two hydroponic experiments. Among these four SNPs, two of them (AX-93636685 and AX-93932874) were also significantly associated with the relative values of shoot P concentration under two P conditions. One SNP AX-93932874 was detected within 5'-untranslated region of Glyma.10 g018800, which contained SPX and RING domains and was named as GmSPX-RING1. Furthermore, the function research of GmSPX-RING1 was carried out in soybean hairy root transformation. Compared with their respective controls, P concentration in GmSPX-RING1 overexpressing transgenic hairy roots was significantly reduced by 32.75% under +P condition; In contrast, P concentration in RNA interference of GmSPX-RING1 transgenic hairy roots was increased by 38.90 and 14.51% under +P and -P conditions, respectively.

CONCLUSIONS

This study shows that the candidate gene GmSPX-RING1 affects soybean phosphorus efficiency by negatively regulating soybean phosphorus concentration in soybean hairy roots. The SNPs and candidate genes identified should be potential for improvement of P efficiency in future soybean breeding programs.

Association analysis between agronomic traits and AFLP markers in a wide germplasm of proso millet (Panicum miliaceum L.) under normal and salinity stress conditions

BACKGROUND

Proso millet is a highly nutritious cereal considered an essential component of processed foods. It is also recognized with high water-use efficiency as well as short growing seasons. This research was primarily aimed at investigating the genetic diversity among genotypes based on evaluating those important traits proposed in previous researches under both normal and salinity- stress conditions. Use of Amplified fragment length polymorphism (AFLP) molecular markers as well as evaluating the association between markers and the investigated traits under both conditions was also another purpose of this research.

RESULTS

According to the phenotypic correlation coefficients, the seed yield had the highest correlation with the forage and biological yields under both conditions. By disintegrating those traits investigated under normal and salinity-stress conditions into principal component analysis, it was found that the first four principal components justified more than 59.94 and 62.48% of the whole variance, respectively. The dendrogram obtained by cluster analysis displayed three groups of genotypes under both normal and salinity- stress conditions. Then, association analyses were conducted on 143 proso millet genotypes and 15 agronomic traits as well as 514 polymorphic AFLP markers (out of 866 created bands) generated by 11 primer combinations (out of the initial 20 primer combinations) EcoRI/MseI. The results obtained by mixed linear model (MLM) indicated that under normal conditions, the M14/E10-45 and M14/E10-60 markers had strong associations with seed yield. A similar trend was also observed for M14/E10-45 and M14/E11-44 markers in relation to forage yield. On the other hand, M14/E10-14, M14/E10-64 markers (for seed yield) and M14/E10-64 marker (for forage yield), had significant and stable association in all environments under salinity-stress conditions. Moreover, a number of markers showed considerable associations and stability under both normal and salinity stress conditions.

CONCLUSIONS

According to the analysis of phenotypic data, the wide germplasm of Iranian proso millet has significant variation in terms of measured traits. It can be concluded that markers showing strong associations with traits under salinity-stress conditions are suitable candidates to be used in future marker-assisted selection (MAS) studies to improve salinity-resistance genotypes of Panicum miliaceum in arid and semiarid areas.

Development and validation of high-throughput and low-cost STARP assays for genes underpinning economically important traits in wheat

We developed and validated 56 gene-specific semi-thermal asymmetric reverse PCR (STARP) markers for 46 genes of important wheat quality, biotic and abiotic stress resistance, grain yield, and adaptation-related traits for marker-assisted selection in wheat breeding. Development of high-throughput, low-cost, gene-specific molecular markers is important for marker-assisted selection in wheat breeding. In this study, we developed 56 gene-specific semi-thermal asymmetric reverse PCR (STARP) markers for wheat quality, tolerance to biotic and abiotic stresses, grain yield, and adaptation-related traits. The STARP assays were validated by (1) comparison of the assays with corresponding diagnostic STS/CAPS markers on 40 diverse wheat cultivars and (2) characterization of allelic effects based on the phenotypic and genotypic data of three segregating populations and 305 diverse wheat accessions from China and 13 other countries. The STARP assays showed the advantages of high-throughput, accuracy, flexibility, simple assay design, low operational costs, and platform compatibility. The state-of-the-art assays of this study provide a robust and reliable molecular marker toolkit for wheat breeding programs.

Meta-QTL analysis and identification of candidate genes for drought tolerance in bread wheat (Triticum aestivum L.)

Meta-QTL (MQTL) analysis for drought tolerance was undertaken in bread wheat to identify consensus and robust MQTLs using 340 known QTLs from 11 earlier studies; 13 MQTLs located on 6 chromosomes (1D, 3B, 5A, 6D, 7A and 7D) were identified, with maximum of 4 MQTLs on chromosome 5A. Mean confidence intervals for MQTLs were much narrower (mean, 6.01 cM; range 2.07-19.46 cM), relative to those in original QTLs (mean, 13.6 cM; range, 1.0-119.1 cM). Two MQTLs, namely MQTL4 and MQTL12, were major MQTLs with potential for use in marker-assisting breeding. As many as 228 candidate genes (CGs) were also identified using 6 of the 13 MQTLs. In-silico expression analysis of these 228 CGs allowed identification of 14 important CGs, with + 3 to - 8 fold change in expression under drought (relative to normal conditions) in a tolerant cv. named TAM107. These CGs encoded proteins belonging to the following families: NAD-dependent epimerase/dehydratase, protein kinase, NAD(P)-binding domain protein, heat shock protein 70 (Hsp70), glycosyltransferase 2-like, etc. Important MQTLs and CGs identified in the present study should prove useful for future molecular breeding and for the study of molecular basis of drought tolerance in cereals in general and wheat in particular.

Rapid Chlorophyll a Fluorescence Light Response Curves Mechanistically Inform Photosynthesis Modeling

Crop improvement is crucial to ensuring global food security under climate change, and hence there is a pressing need for phenotypic observations that are both high throughput and improve mechanistic understanding of plant responses to environmental cues and limitations. In this study, chlorophyll a fluorescence light response curves and gas-exchange observations are combined to test the photosynthetic response to moderate drought in four genotypes of Brassica rapa The quantum yield of PSII (ϕ PSII ) is here analyzed as an exponential decline under changing light intensity and soil moisture. Both the maximum ϕ PSII and the rate of ϕ PSII decline across a large range of light intensities (0-1,000 μmol photons m-2 s-1; β PSII ) are negatively affected by drought. We introduce an alternative photosynthesis model (β PSII model) incorporating parameters from rapid fluorescence response curves. Specifically, the model uses β PSII as an input for estimating the photosynthetic electron transport rate, which agrees well with two existing photosynthesis models (Farquhar-von Caemmerer-Berry and Yin). The β PSII model represents a major improvement in photosynthesis modeling through the integration of high-throughput fluorescence phenotyping data, resulting in gained parameters of high mechanistic value.

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.