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

Integrating QTL mapping and GWAS to decipher the genetic mechanisms behind the calcium contents of Brassica napus shoots

Brassica napus is an important oil crop worldwide, and its shoots are rich in vitamin C, calcium, and selenium. Functional oilseed-vegetable-dual-purpose varieties can increase the subsidiary value of B. napus. Consumption of high-calcium B. napus shoots can effectively help provide essential elements to the human body. To investigate the genetic mechanisms underlying the calcium concentrations of B. napus shoots, quantitative trait loci (QTL) mapping, using a population of 189 recombinant inbred lines, and a genome-wide association study, using an association panel of 202 diverse accessions, were performed. A total of 12 QTLs controlling calcium content were identified using the recombinant inbred line population in five environments. Among them, qCaC.22GY-A05-1 was considered the major QTL, with a phenotypic variation of 10.10%. In addition, 228 single nucleotide polymorphisms significantly related to calcium content were identified using the genome-wide association study in six environments, and they were distributed on all of the chromosomes, except A10. Finally, 10 candidate genes involved in regulating calcium absorption and transport in B. napus shoots were identified. However, no overlapping intervals were found through a comprehensive analysis of the two datasets. These results provide valuable information for understanding the genetic control of calcium concentration in B. napus shoots.

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

KEY MESSAGE

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

Dissection of the genetic basis and molecular mechanism of ovule number per ovary in oilseed rape (Brassica napus L.)

Ovule number per ovary (ONPO) determines the maximum potential of seed number per fruit that is a direct component of seed yield in crops. This study aimed to dissect the genetic basis and molecular mechanism of ONPO using a newly developed doubled haploid (DH) population in oilseed rape. In all the four investigated environments, the ONPO of 201 DH lines exhibited normal distribution with a wide variation from 22.6 to 41.8, suggesting quantitative inheritance appropriate for mapping QTL. A skeleton genetic map of 2111 markers within 19 linkage groups was developed, with a total of 1715.71 cM in length and an average of 0.82 cM between markers. Linkage mapping identified ten QTLs that were distributed on eight chromosomes and explained 7.0-15.9% of the phenotypic variance. Among these, four were identical to the reported and two were repeatedly detected with relatively large effects, highlighting their potential for marker-assisted selection. Phytohormone quantification of ovaries (at the ovule initiation stage) from two pools of high and low ONPO lines showed significant differences in the levels of nine sub-types of phytohormones, suggesting their important roles in regulating ovule number. Transcriptomic analysis identified 7689 differentially expressed genes (DEGs) between the two pools, of which nearly half were enriched into functional categories of reported genes regulating ONPO, including protein, RNA, signalling, miscellaneous, development, hormone metabolism, and tetrapyrrole synthesis. Integration of linkage QTL mapping, transcriptome sequencing and BLAST analysis identified 15 homologues of reported ovule number genes and 327 DEGs in the QTL regions, which were considered as direct and potential candidate genes. These findings propose further insights into the genetic basis and molecular mechanisms of ONPO, which will facilitate future gene cloning and genetic improvement for enhancing seed yield in oilseed rape.

Genome-wide screening of the RNase T2 gene family and functional analyses in jujube (Ziziphus jujuba Mill.)

BACKGROUND

Ribonuclease (RNase T2) plays crucial roles in plant evolution and breeding. However, there have been few studies on the RNase T2 gene family in Ziziphus jujuba Mill., one of important dried fruit tree species. Recently, the released sequences of the reference genome of jujube provide a good chance to perform genome-wide identification and characterization of ZjRNase gene family in the jujube.

RESULTS

In this study, we identified four members of RNase T2 in jujube distributed on three chromosomes and unassembled chromosomes. They all contained two conserved sites (CASI and CASII). Analysis of the phylogenetic relationships revealed that the RNase T2 genes in jujube could be divided into two groups: ZjRNase1 and ZjRNase2 belonged to class I, while ZjRNase3 and ZjRNase4 belonged to class II. Only ZjRNase1 and ZjRNase2 expression were shown by the jujube fruit transcriptome analysis. So ZjRNase1 and ZjRNase2 were selected functional verification by overexpression transformation of Arabidopsis. The overexpression of these two genes led to an approximately 50% reduction in seed number, which deserve further attention. Moreover, the leaves of the ZjRNase1 overexpression transgenic lines were curled and twisted. Overexpression of ZjRNase2 resulted in shortened and crisp siliques and the production of trichomes, and no seeds were produced.

CONCLUSION

In summary, these findings will provide new insights into the molecular mechanisms of low number of hybrid seeds in jujube and a reference for the future molecular breeding of jujube.

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.

Genome-wide characterization of ovate family protein gene family associated with number of seeds per silique in Brassica napus

Ovate family proteins (OFPs) were firstly identified in tomato as proteins controlling the pear shape of the fruit. Subsequent studies have successively proved that OFPs are a class of negative regulators of plant development, and are involved in the regulation of complex traits in different plants. However, there has been no report about the functions of OFPs in rapeseed growth to date. Here, we identified the OFPs in rapeseed at the genomic level. As a result, a total of 67 members were obtained. We then analyzed the evolution from Arabidopsis thaliana to Brassica napus, illustrated their phylogenetic and syntenic relationships, and compared the gene structure and conserved domains between different copies. We also analyzed their expression patterns in rapeseed, and found significant differences in the expression of different members and in different tissues. Additionally, we performed a GWAS for the number of seeds per silique (NSPS) in a rapeseed population consisting of 204 natural accessions, and identified a new gene BnOFP13_2 significantly associated with NSPS, which was identified as a novel function of OFPs. Haplotype analysis revealed that the accessions with haplotype 3 had a higher NSPS than other accessions, suggesting that BnOFP13_2 is associated with NSPS. Transcript profiling during the five stages of silique development demonstrated that BnOFP13_2 negatively regulates NSPS. These findings provide evidence for functional diversity of OFP gene family and important implications for oilseed rape breeding.

Genetic dissection of the natural variation of ovule number per ovary in oilseed rape germplasm (Brassica napus L.)

Oilseed rape is one of the world's largest oil and industrial crops, providing humans with various products, such as vegetable oil and biofuel. Ovules are the direct precursors of seeds, and ovule number per ovary (ONPO) largely determines seed number per fruit that affects both yield and fitness of seed crops. The ONPO shows wide variation in oilseed rape, whereas the underlying genes and mechanisms are poorly known. The present study performed the genetic, physiological and transcriptomic analyses of ovule number per ovary using an association panel and the extreme lines. The ONPO of 327 accessions planted in four environments showed a large variation from 19.2 to 43.8, indicating a great potential for the further genetic improvement of ovule number. The genome-wide association study (GWAS) identified a total of 43 significant SNP markers. Further, these SNPs were integrated into 18 association loci, which were distributed on chromosomes A01, A03, A06, A07, A09, C01, C03, C06, C07, and C09, explaining 4.3-11.5% of the phenotypic variance. The ONPO decreased as their appearance order on the inflorescence and was associated with the level of several types of endogenous phytohormones but not related to leaf area and photosynthetic rate. Comparative transcriptomic analysis identified a total of 4,449 DEGs enriched in 30 classes, including DNA, RNA, protein, signaling, transport, development, cell wall, lipid metabolism, and secondary metabolism. Nearly half of DEGs were involved in the known pathways in regulating ovule number, of which 12 were homologous to know ovule number regulating genes, indicating a strong link between the identified DEGs and ovule number. A total of 73 DEGs were located within the genomic regions of association loci, of which six were identified as candidates based on functional annotation. These results provide useful information for the further genetic improvement of ovule and seed number in oilseed rape.

A Preliminary Study for Identifying Quantitative Trait Loci Associated with Seed Production in Radish Using Genotyping-by-Sequencing

The high yield of seeds can reduce the cost of seed production for parental lines, as well as F1 cultivars in radish. The number of seeds per silique and silique length are two important traits among traits determining seed yield, but no study has been conducted on their quantitative trait loci (QTLs) in radish. A high-density linkage map was constructed, based on genotyping-by-sequencing (GBS) of the F2 population, derived from two parental lines, significantly differed by the two traits, which were grown in a controlled environment to minimize the environmental effects. Using the map with 848 SNPs, three significant QTLs were identified, two and one of which were associated with the number of seeds per silique and silique length, respectively. Ortholog analysis was conducted with Arabidopsis thaliana genes, related to the number of seeds per silique, and revealed five radish putative candidate genes. These putative candidate genes appear to be related to ovule, embryo sac, embryo, pollen and seed development, as well as a double fertilization process. The method to pollinate the F2 population, as well as preliminary QTLs and SNPs therein, can be helpful for future QTL studies to improve seed production in radish breeding programs.

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

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.

Genome-wide association study for candidate genes controlling seed yield and its components in rapeseed (Brassica napus subsp. napus)

Genetic improvement of seed yield per plant (SY) is one of the major objectives in Brassica napus breeding programme. SY, being a complex quantitative trait is directly and indirectly influenced by yield-component traits such as siliqua length (SL), number of seeds per siliqua (NSS), and thousand seed weight (TSW). Therefore, concurrent improvement in SL, NSS and TSW can lead to higher SY in B. napus. This study was conducted to identify significant SNPs and putative candidate genes governing SY and its component traits (SL, NSS, TSW). All these traits were evaluated in a diverse set of 200 genotypes representing diversity from wide geographical locations. Of these, a set of 125 genotypes were chosen based on pedigree diversity and multi-location trait variation for genotyping by sequencing (GBS). Best linear unbiased predictors (BLUPs) of all the traits were used for genome-wide association study (GWAS) with 85,126 SNPs obtained from GBS. A total of 16, 18, 27 and 18 SNPs were found to be significantly associated for SL, NSS, TSW and SY respectively. Based on linkage disequilibrium decay analysis, 150 kb genomic region flanking the SNP was used for the identification of underlying candidate genes for each test trait. Important candidate genes involved in phytohormone signaling (WAT1, OSR1, ARR8, CKX1, REM7, REM9, BG1) and seed storage proteins (Cruciferin) were found to have significant influence on seed weight and yield. Genes involved in sexual reproduction and fertilization (PERK7, PERK13, PRK3, GATA15, NFD6) were found to determine the number of seeds per siliqua. Several genes found in this study namely ATS3A, CKX1, SPL2, SPL6, SPL9, WAT1 showed pleiotropic effect with yield component traits. Significant SNPs and putative candidate genes identified for SL, NSS, TSW and SY could be used in marker-assisted breeding for improvement of crop yield in B. napus. Genotypes identified with high SL, NSS, TSW and SY could serve as donors in crop improvement programs in B. napus.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01060-9.

Shotgun proteomics of Brassica rapa seed proteins identifies vicilin as a major seed storage protein in the mature seed

Proteins make up a large percentage of the Brassica seed and are second only to the oil in economic importance with uses for both animal and human nutrition. The most abundant proteins reported in the seeds of Brassica are the seed storage proteins cruciferin and napin, belonging to the 12S globulin and 2S albumin families of proteins, respectively. To gain insight into the Brassica rapa seed proteome and to confirm the presence and relative quantity of proteins encoded by candidate seed storage genes in the mature seed, shotgun proteomics was carried out on protein extracts from seeds of B. rapa inbred line R-o-18. Following liquid chromatography tandem mass spectrometry, a total of 34016 spectra were mapped to 323 proteins, where 233 proteins were identified in 3 out of 4 biological replicates by at least 2 unique peptides. 2S albumin like napin seed storage proteins (SSPs), 11/12S globulin like cruciferin SSPs and 7S globulin like vicilin SSPs were identified in the samples, along with other notable proteins including oil body proteins, namely ten oleosins and two oil body-associated proteins. The identification of vicilin like proteins in the mature B. rapa seed represents the first account of these proteins in the Brassicaceae and analysis indicates high conservation of sequence motifs to other 7S vicilin-like allergenic proteins as well as conservation of major allergenic epitopes in the proteins. This study enriches our existing knowledge on rapeseed seed proteins and provides a robust foundation and rational basis for plant bioengineering of seed storage proteins.

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.