Multi-omics analysis reveals the mechanism of seed coat color formation in Brassica rapa L.

KEYMESSAGE: Multi-omics analysis of the transcriptome, metabolome and genome identified major and minor loci and candidate genes for seed coat color and explored the mechanism of flavonoid metabolites biosynthesis in Brassica rapa. Yellow seed trait is considered an agronomically desirable trait with great potential for improving seed quality of Brassica crops. Mechanisms of the yellow seed trait are complex and not well understood. In this study, we performed an integrated metabolome, transcriptome and genome-wide association study (GWAS) on different B. rapa varieties to explore the mechanisms underlying the seed coat color formation. A total of 2,499 differentially expressed genes and 116 differential metabolites between yellow and black seeds with strong association with the flavonoid biosynthesis pathway was identified. In addition, 330 hub genes involved in the seed coat color formation, and the most significantly differential flavonoids biosynthesis were detected based on weighted gene co-expression network analysis. Metabolite GWAS analysis using the contents of 42 flavonoids in developing seeds of 159 B. rapa lines resulted in the identification of 1,626 quantitative trait nucleotides (QTNs) and 37 chromosomal intervals, including one major locus on chromosome A09. A combination of QTNs detection, transcriptome and functional analyses led to the identification of 241 candidate genes that were associated with different flavonoid metabolites. The flavonoid biosynthesis pathway in B. rapa was assembled based on the identified flavonoid metabolites and candidate genes. Furthermore, BrMYB111 members (BraA09g004490.3C and BraA06g034790.3C) involved in the biosynthesis of taxifolin were functionally analyzed in vitro. Our findings lay a foundation and provide a reference for systematically investigating the mechanism of seed coat color in B. rapa and in the other plants.

Local Housing Market Sentiments and Returns: Evidence from China.

This paper examines the impacts of local housing sentiments on the housing price dynamics of China. With a massive second-hand transaction dataset, we construct monthly local housing sentiment indices for 18 major cities in China from January 2016 to October 2020. We create three sentiment proxies representing the local housing market liquidity and speculative behaviors from the transaction dataset and then use partial least squares (PLS) to extract a recursive look-ahead-bias-free local housing sentiment index for each city considered. The local housing sentiments are shown to have robust predictive powers for future housing returns with a salient short-run underreaction and long-run overreaction pattern. Further analysis shows that local housing sentiment impacts are asymmetric, and housing returns in cities with relatively inelastic housing supply are more sensitive to local housing sentiments. We also document a significant feedback effect between housing returns and market sentiments, indicating the existence of a pricing-sentiment spiral which could potentially enhance the ongoing market fever of Chinese housing markets. The main estimation results are robust to alternative sentiment extraction methods and alternative sentiment proxies, and consistent for the sample period before COVID-19.

Metabolite Profiling and Transcriptome Analysis Provide Insight into Seed Coat Color in Brassica juncea.

The allotetraploid species Brassica juncea (mustard) is grown worldwide as oilseed and vegetable crops; the yellow seed-color trait is particularly important for oilseed crops. Here, to examine the factors affecting seed coat color, we performed a metabolic and transcriptomic analysis of yellow- and dark-seeded B. juncea seeds. In this study, we identified 236 compounds, including 31 phenolic acids, 47 flavonoids, 17 glucosinolates, 38 lipids, 69 other hydroxycinnamic acid compounds, and 34 novel unknown compounds. Of these, 36 compounds (especially epicatechin and its derivatives) accumulated significantly different levels during the development of yellow- and dark-seeded B. juncea. In addition, the transcript levels of BjuDFR, BjuANS,BjuBAN, BjuTT8, and BjuTT19 were closely associated with changes to epicatechin and its derivatives during seed development, implicating this pathway in the seed coat color determinant in B. juncea. Furthermore, we found numerous variations of sequences in the TT8A genes that may be associated with the stability of seed coat color in B. rapa, B. napus, and B. juncea, which might have undergone functional differentiation during polyploidization in the Brassica species. The results provide valuable information for understanding the accumulation of metabolites in the seed coat color of B. juncea and lay a foundation for exploring the underlying mechanism.

Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress.

BACKGROUND: Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea). RESULTS: In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content. CONCLUSIONS: Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

Graphene oxide induced assembly and crumpling of Co3O4 nanoplates.

Cobalt (II, III) oxide (Co3O4) has been widely studied and applied in various fields, however, it suffers from slow mass and electron transfer during applications. Herein, crumpled Co3O4 and Co3O4/reduced graphene oxide (rGO) with tunable 2D-in-3D structures were prepared by combining spray pyrolysis with a graphene oxide (GO) template. The 2D Co3O4 nanoplates were interconnected with each other to form a 3D ball with many wrinkles, resulting in defect enrichment on the abundant boundaries of the nanosheets, which provided more active sites for catalytic reactions. In addition, the unique 2D-in-3D structure allowed fast mass transfer and structural stability. Furthermore, the assembled structure could be understood as being composed of uniformly distributed oxygen-containing functional groups pinning metal cations on the GO surface through electrostatic interaction, and the 2D structure of the GO enabled the in situ converted Co3O4 to grow along the GO surface with excellent dispersion. Taking advantage of the above, the Co3O4/rGO balls demonstrated an excellent oxygen evolution reaction performance, an overpotential of 298 mV at a current density of 10.0 mA cm-2 and a current density of 115.9 mA cm-2 at the overpotential of η = 500 mV.

Genome-Wide Analysis of Phosphorus Transporter Genes in Brassica and Their Roles in Heavy Metal Stress Tolerance.

Phosphorus transporter (PHT) genes encode H2PO4-/H+ co-transporters that absorb and transport inorganic nutrient elements required for plant development and growth and protect plants from heavy metal stress. However, little is known about the roles of PHTs in Brassica compared to Arabidopsis thaliana. In this study, we identified and extensively analyzed 336 PHTs from three diploid (B. rapa, B. oleracea, and B. nigra) and two allotetraploid (B. juncea and B. napus) Brassica species. We categorized the PHTs into five phylogenetic clusters (PHT1-PHT5), including 201 PHT1 homologs, 15 PHT2 homologs, 40 PHT3 homologs, 54 PHT4 homologs, and 26 PHT5 homologs, which are unevenly distributed on the corresponding chromosomes of the five Brassica species. All PHT family genes from Brassica are more closely related to Arabidopsis PHTs in the same vs. other clusters, suggesting they are highly conserved and have similar functions. Duplication and synteny analysis revealed that segmental and tandem duplications led to the expansion of the PHT gene family during the process of polyploidization and that members of this family have undergone purifying selection during evolution based on Ka/Ks values. Finally, we explored the expression profiles of BnaPHT family genes in specific tissues, at various developmental stages, and under heavy metal stress via RNA-seq analysis and qRT-PCR. BnaPHTs that were induced by heavy metal treatment might mediate the response of rapeseed to this important stress. This study represents the first genome-wide analysis of PHT family genes in Brassica species. Our findings improve our understanding of PHT family genes and provide a basis for further studies of BnaPHTs in plant tolerance to heavy metal stress.

Testis-specific calcium-binding protein CBP86-IV (CABYR) binds with phosphoglycerate kinase 2 in vitro and in vivo experiment.

The investigation of the interacting proteins with testis-specific calcium-binding protein CBP86-IV (CABYR) was carried out in human spermatozoa. The total RNA from human spermatozoa was extracted, and the ORF sequence of TSCBP86-IV gene was amplified and cloned into expression vector pET-28a. The positive recombinant clones were transformed into Escherichia coli strain BL21 (DE3) to express fusion protein. Then, co-immunoprecipitation (Co-IP) of TSCBP86-IV was performed in BL21 cell lysate expressing CBP86-IV recombinant protein. The immune complex was captured and identified by mass spectrometry. Reverse Co-IP of potential interacting proteins was performed in human sperm cell lysate. The potential protein interactions were confirmed by yeast two-hybrid system. Thirteen proteins were successfully identified in immune complex from E. coli cell lysate. Phosphoglycerate kinase 2 (PGK2) further showed positive results both in reverse Co-IP and yeast two-hybrid experiments and was confirmed to be interacted with TSCBP86-IV in human sperm cells.

Low-expressed testis-specific calcium-binding protein CBP86-IV (CABYR) is observed in idiopathic asthenozoospermia.

OBJECTIVES: To investigate the expression level of testis-specific calcium-binding protein CBP86-IV in normal and asthenozoospermic human sperm. METHODS: The total RNA was extracted from human sperm, and target cDNA was obtained by reverse transcription-polymerase chain reaction. Then the cDNA was used for quantitative PCR analysis and cloned into the prokaryotic expression vector pET-28a, respectively. The fusion protein was induced and expressed as inclusion body which was used to produce the polyclonal antibody against TSCBP86-IV. The protein expression level of TSCBP86-IV from normal human sperm and idiopathic asthenozoospermic samples was detected by the purified antibody. RESULTS: The experimental results showed that the protein expression of TSCBP86-IV was reduced in idiopathic asthenozoospermia and consistent with the transcriptional changing tendency which was detected by quantitative PCR analysis. CONCLUSIONS: The stable and reliable change of TSCBP86-IV may be taken as a new molecular marker for clinical diagnosis of idiopathic asthenozoospermia.

MgO@polydopamine Nanoparticle-Loaded Photothermal Microneedle Patches Combined with Chitosan Gel Dressings for the Treatment of Infectious Wounds.

As for wound drug delivery, microneedles (MNs) have attracted wide attention. However, while effective at increasing the depth of drug delivery, traditional MNs often have limited drug loads and have difficulty penetrating scabs on wounds. Herein, we develop a drug delivery system combining MgO@polydopamine (MgO@PDA) nanoparticle-loaded photothermal MN patches and chitosan (CS) gel to inhibit the formation of scabs and deliver sufficient drugs into deep tissue. When inserted into the wound, the MN system can keep the wound bed moist and weakly acidic to inhibit the formation of scabs and accelerate wound closure. The released MgO@PDA nanoparticles from both the tips and the backing layer, which immensely increase the drug load, continuously release Mg2+ in the moist, weakly acidic wound bed, promoting tissue migration and the formation of microvessels. MgO@PDA nanoparticles show excellent antibacterial activity under near-infrared irradiation synergized with the CS gel, and the PDA coating can also overcome the adverse effects of oxidative stress. Through in vitro and in vivo experiments, the MN system showed remarkable antibacterial, antioxidant, anti-inflammatory, and pro-angiogenic effects, indicating its potential in the treatment of infectious wounds.

Enhanced Production of Sisomicin in Micromonospora inyoensis by Protoplast Mutagenesis and Fermentation Optimization.

Sisomicin is a broad-spectrum aminoglycoside antibiotic and is the precursor of netilmicin and plazomicin. However, the fermentation level of sisomicin is still low compared with other antibiotics, which restricts the application of sisomicin and its derivatives. In this study, to improve sisomicin production, breeding of high-yielding sisomicin strains was conducted with chemical mutagenesis using Micromonospora inyoensis OG-1 (titer, 1042 U·mL-1) as the starting strain. Protoplast preparation was conducted under optimal conditions, and protoplast mutagenesis was performed with a suitable concentration of diethyl sulfate. Subsequently, a high-yielding and genetically stable strain (H6-32) was obtained by screening, with a sisomicin titer of 1486 U·mL-1 (an increase of 42.6%). Finally, carbon and nitrogen sources were optimized to further improve sisomicin production, and a sisomicin titer of 1780 U·mL-1 was ultimately obtained by controlling the dissolved oxygen level at 30% in a 5-L fermenter, which to the best of our knowledge is the highest reported titer ever achieved by fermentation. Comparative genome analysis showed that a total of 13 genes in the genome of the mutant strain H6-32 were mutated compared to the original strain. This study not only provides a reference for further breeding of high-yielding strains and fermentation optimization, but also enhances our understanding of sisomicin production.

Fabrication of controlled porous and ultrafast dissolution porous microneedles by organic-solvent-free ice templating method.

Porous Microneedles (PMNs) have been widely used in drug delivery and medical diagnosis owing to their abundant interconnected pores. However, the mechanical strength, the use of organic solvent, and drug loading capacity have long been challenging. Herein, a novel strategy of PMNs fabrication based on the Ice Templating Method is proposed that is suitable for insoluble, soluble, and nanosystem drug loading. The preparation process simplifies the traditional microneedle preparation process with a shorter preparation time. It endows the highly tunable porous morphology, enhanced mechanical strength, and rapid dissolution performance. Micro-CT three-dimensional reconstruction was used to better quantify the internal structures of PMNs, and we further established the equivalent pore network model to statistically analyze the internal pore structure parameters of PMNs. In particular, the mechanical strength is mainly negatively correlated with the surface porosity, while the dissolution velocity is mainly positively correlated with the permeability coefficient by the correlation heatmap. The poorly water-soluble Asiatic acid was encapsulated in PMNs in nanostructured lipid carriers, showing prominent hypertrophic scar healing trends. This work offers a quick and easy way of preparation that may be used to expand PMNs function and be introduced in industrial manufacturing development.

Genome-Wide Identification of B-Box Family Genes and Their Potential Roles in Seed Development under Shading Conditions in Rapeseed.

B-box (BBX) proteins, a subfamily of zinc-finger transcription factors, are involved in various environmental signaling pathways. In this study, we conducted a comprehensive analysis of BBX family members in Brassica crops. The 482 BBX proteins were divided into five groups based on gene structure, conserved domains, and phylogenetic analysis. An analysis of nonsynonymous substitutions and (Ka)/synonymous substitutions (Ks) revealed that most BBX genes have undergone purifying selection during evolution. An analysis of transcriptome data from rapeseed (Brassica napus) organs suggested that BnaBBX3d might be involved in the development of floral tissue-specific RNA-seq expression. We identified numerous light-responsive elements in the promoter regions of BnaBBX genes, which were suggestive of participation in light signaling pathways. Transcriptomic analysis under shade treatment revealed 77 BnaBBX genes with significant changes in expression before and after shading treatment. Of these, BnaBBX22e showed distinct expression patterns in yellow- vs. black-seeded materials in response to shading. UPLC-HESI-MS/MS analysis revealed that shading influences the accumulation of 54 metabolites, with light response BnaBBX22f expression correlating with the accumulation of the flavonoid metabolites M46 and M51. Additionally, BnaBBX22e and BnaBBX22f interact with BnaA10.HY5. These results suggest that BnaBBXs might function in light-induced pigment accumulation. Overall, our findings elucidate the characteristics of BBX proteins in six Brassica species and reveal a possible connection between light and seed coat color, laying the foundation for further exploring the roles of BnaBBX genes in seed development.

Flexible microneedles incorporating gold nanorods and tacrolimus for effective synergistic photothermal-chemotherapy of rheumatoid arthritis.

Transdermal drug delivery systems for rheumatoid arthritis (RA) have garnered substantial attention due to their great potential to overcome limitations observed in conventional oral or injection strategies, including limited selectivity and adverse effects on extra-articular tissues. Microneedles (MNs) appear to be highly desirable carriers for transdermal drug delivery of RA. However, microneedles typically are unable to keep up with the flexibility of joints, which decreases the effectiveness of administration. In this study, we developed a flexible microneedles (FMNs) delivery system. And gelatin was employed for the fabrication of flexible backings for microneedles owing to its excellent ductility and biocompatibility. We achieved synergisticphotothermal-chemotherapy of RA by incorporating the chemical drug Tacrolimus (TAC) and the photothermal agent gold nanorods (AuNRs) into dissolving microneedles. Results showed a high mechanical strength of the proposed FMNs. In the animal model of adjuvant-induced arthritis (AA), it is indicated that the prepared FMNs inhibited the expression of related inflammatory cytokines such as IL-1ß and TNF-α while enhancing bone repair and other related factors. Thus, the combination therapy of FMNs-mediated hyperthermia and chemotherapy can serve as a novel and synergistic treatment option for RA.

Methotrexate-loaded hyaluronan-modified liposomes integrated into dissolving microneedles for the treatment of psoriasis.

Methotrexate (MTX) is a first-line drug in treating psoriasis because of its strong anti-proliferation and anti-inflammatory effects. However, systemic administration of MTX will lead to many side effects, such as gastrointestinal irritation, liver and kidney toxicity, etc. Herein, we developed liposome-loaded microneedles (MNs) system to improve transdermal efficiency, which was used to overcome the problems of low transdermal efficiency and poor therapeutic effect of traditional transdermal drug delivery methods. Hyaluronic acid (HA) was modified on the surface of MTX-loaded liposomes. The interaction of HA and CD44 could increase the adhesion of HA-MTX-Lipo to HaCaT cells, thereby promoting the apoptosis or death of HaCaT cells. Results indicated HA-MTX-Lipo MNs could inhibit the development of psoriasis and reduce the degree of skin erythema, scaling, and thickening. The mRNA levels of proinflammatory cytokines such as IL-17A, IL-23, and TNF-α were decreased. The epidermal thickness and proliferative cell-associated antigen Ki67 expression were also reduced. Specifically, the expression of mRNA levels of proinflammatory cytokines was down-regulated. The MNs transdermal delivery of HA-modified-MTX liposomes provided a promising method for treating psoriasis.

Comparative proteomic study between human normal motility sperm and idiopathic asthenozoospermia.

PURPOSE: Idiopathic asthenozoospermia is considered as one of the causes of male infertility and characterized by reduced sperm motility. For a better determination of pathogenic mechanism of asthenozoospermia, the exploration of differentially expressed proteins in normal sperm motility and idiopathic asthenozoospermia was conducted in our study. METHODS: Sperm proteins were extracted and isolated by two-dimensional electrophoresis. All significantly changed protein spots were picked up from 2D gels and identified by tandem mass spectrometry. Sixteen of the thirty-three total differentially expressed protein spots were successfully identified by matrix-assisted laser desorption/ionization time-of-flight tandem mass spectrometry. RESULTS: Sixteen proteins identified belonged to 15 unique protein groups. GRP78, lactoferrin, SPANXB, PGK2, flagellin, DJ-1, XPA binding protein 2, CAB2, GPX4, and GAPDH were the first to be identified as differentially expressed proteins in idiopathic asthenospermia patients. Meanwhile, the analysis of quantitative RT-PCR was carried out to compare the protein levels, and the results indicated that the expression levels of the gene and protein were not entirely consistent. CONCLUSIONS: These experimental results expand the scope of the protein database, generating targets for further investigation of the pathogenic mechanism of idiopathic asthenozoospermia.

[Low-dose testosterone undecanoate capsules combined with tadalafil for late-onset hypogonadism accompanied with ED].

OBJECTIVE: To observe the clinical effect of low-dose testosterone undecanoate capsules combined with tadalafil on late-onset hypogonadism (LOH) accompanied with ED. METHODS: Ninety cases of LOH accompanied with ED who met the inclusion criteria were randomly divided into a control group and a combination therapy group, the former treated with tadalafil and the latter with low-dose testosterone undecanoate capsules combined with tadalafil. The LOH symptoms, IIEF-5 scores, sexual encounter profile (SEP) scores, prostate volumes, and the levels of total testosterone (TT), free testosterone (FT) and prostatic specific antigen (PSA) were recorded and compared between the two groups before and after treatment. RESULTS: The IIEF-5 and SEP scores and the TT and FT levels were 20.6 +/- 3.8, 4.02 +/- 1.08, (15.4 +/- 3.4) nmol/L and (0.391 +/- 0.062) nmol/L, respectively, in the combination therapy group after treatment, significantly higher both than 15.7 +/- 3.9, 1.49 +/- 0.82, (10.1 +/- 1.2) nmol/L and (0.200 +/- 0.045) nmol/L before treatment (P < 0.01) and than 8.6 +/- 3.6, 3.50 +/- 1.21, (10.2 +/- 1.2) nmol/L and (0.210 +/- 0.051) nmol/L in the control group after treatment (P < 0.01). CONCLUSION: Low-dose testosterone undecanoate capsules combined with tadalafil has a definite clinical effect and no obvious adverse reactions in the treatment of LOH accompanied with ED.

High-Speed railways and the spread of Covid-19.

High-speed railways (HSRs) greatly decrease transportation costs and facilitate the movement of goods, services, and passengers across cities. In the context of the Covid-19 pandemic, however, HSRs may contribute to the cross-regional spread of the new coronavirus. This paper evaluates the role of HSRs in spreading Covid-19 from Wuhan to other Chinese cities. We use train frequencies in 1971 and 1990 as instrumental variables. Empirical results from gravity models demonstrate that one more HSR train originating from Wuhan each day before the Wuhan lockdown increases the cumulative number of Covid-19 cases in a city by about 10 percent. The empirical analysis suggests that other transportation modes, including normal-speed trains and airline flights, also contribute to the spread of Covid-19, but their effects are smaller than the effect of HSRs. This paper's findings indicate that transportation infrastructures, especially HSR trains originating from a city where a pandemic broke out, can be important factors promoting the spread of an infectious disease.

Four New Species of Small-Spored Alternaria Isolated from Solanum tuberosum and S. lycopersicum in China.

Small-spored Alternaria species have been frequently isolated from diseased leaves of Solanum plants. To clarify the diversity of small-spored Alternaria species, a total of 118 strains were obtained from leaf samples of S. tuberosum and S. lycopersicum in six provinces of China during 2022-2023. Based on morphological characterization and multi-locus phylogenetic analysis of the internal transcribed spacer of the rDNA region (ITS), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), translation elongation factor 1 alpha (TEF1), RNA polymerase second largest subunit (RPB2), Alternaria major allergen gene (Alt a 1), endopolygalacturonase gene (EndoPG) and an anonymous gene region (OPA10-2), seven species were determined, including four novel species and three known species (A. alternata, A. gossypina and A. arborescens). The novel species were described and illustrated as A. longxiensis sp. nov., A. lijiangensis sp. nov., A. lycopersici sp. nov. and A. solanicola sp. nov.. In addition, the pathogenicity of the seven species was evaluated on potato leaves. The species exhibited various aggressiveness, which could help in disease management.

Identification of the High-Affinity Potassium Transporter Gene Family (HKT) in Brassica U-Triangle Species and Its Potential Roles in Abiotic Stress in Brassica napus L.

Members of the high-affinity potassium transporter (HKT) protein family regulate the uptake and homeostasis of sodium and potassium ions, but little research describes their roles in response to abiotic stresses in rapeseed (Brassica napus L.). In this study, we identified and characterized a total of 36 HKT genes from the species comprising the triangle of U model (U-triangle species): B. rapa, B. nigra, B. oleracea, B. juncea, B. napus, and B. carinata. We analyzed the phylogenetic relationships, gene structures, motif compositions, and chromosomal distributions of the HKT family members of rapeseed. Based on their phylogenetic relationships and assemblage of functional domains, we classified the HKT members into four subgroups, HKT1;1 to HKT1;4. Analysis of the nonsynonymous substitutions (Ka), synonymous substitutions (Ks), and the Ka/Ks ratios of HKT gene pairs suggested that these genes have experienced strong purifying selective pressure after duplication, with their evolutionary relationships supporting the U-triangle theory. Furthermore, the expression profiles of BnaHKT genes varies among potassium, phytohormone and heavy-metal treatment. Their repression provides resistance to heavy-metal stress, possibly by limiting uptake. Our results systematically reveal the characteristics of HKT family proteins and their encoding genes in six Brassica species and lay a foundation for further exploration of the role of HKT family genes in heavy-metal tolerance.

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

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