Biotic interactions promote local adaptation to soil in plants.

Although different ecological factors shape adaptative evolution in natural habitats, we know little about how their interactions impact local adaptation. Here we used eight generations of experimental evolution with outcrossing Brassica rapa plants as a model system, in eight treatment groups that varied in soil type, herbivory (with/without aphids), and pollination mode (hand- or bumblebee-pollination), to study how biotic interactions affect local adaptation to soil. First, we show that several plant traits evolved in response to biotic interactions in a soil-specific way. Second, using a reciprocal transplant experiment, we demonstrate that significant local adaptation to soil-type evolved in the "number of open flowers", a trait used as a fitness proxy, but only in plants that evolved with herbivory and bee pollination. Whole genome re-sequencing of experimental lines revealed that biotic interactions caused a 10-fold increase in the number of SNPs across the genome with significant allele frequency change, and that alleles with opposite allele frequency change in different soil types (antagonistic pleiotropy) were most common in plants with an evolutionary history of herbivory and bee pollination. Our results demonstrate that the interaction with mutualists and antagonists can facilitate local adaptation to soil type through antagonistic pleiotropy.

Temperature-dependent jumonji demethylase modulates flowering time by targeting H3K36me2/3 in Brassica rapa.

Global warming has a severe impact on the flowering time and yield of crops. Histone modifications have been well-documented for their roles in enabling plant plasticity in ambient temperature. However, the factor modulating histone modifications and their involvement in habitat adaptation have remained elusive. In this study, through genome-wide pattern analysis and quantitative-trait-locus (QTL) mapping, we reveal that BrJMJ18 is a candidate gene for a QTL regulating thermotolerance in thermotolerant B. rapa subsp. chinensis var. parachinensis (or Caixin, abbreviated to Par). BrJMJ18 encodes an H3K36me2/3 Jumonji demethylase that remodels H3K36 methylation across the genome. We demonstrate that the BrJMJ18 allele from Par (BrJMJ18Par) influences flowering time and plant growth in a temperature-dependent manner via characterizing overexpression and CRISPR/Cas9 mutant plants. We further show that overexpression of BrJMJ18Par can modulate the expression of BrFLC3, one of the five BrFLC orthologs. Furthermore, ChIP-seq and transcriptome data reveal that BrJMJ18Par can regulate chlorophyll biosynthesis under high temperatures. We also demonstrate that three amino acid mutations may account for function differences in BrJMJ18 between subspecies. Based on these findings, we propose a working model in which an H3K36me2/3 demethylase, while not affecting agronomic traits under normal conditions, can enhance resilience under heat stress in Brassica rapa.

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa.

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

Brassica rapa CURLY LEAF is a major H3K27 methyltransferase regulating flowering time.

MAIN CONCLUSION: In Brassica rapa, the epigenetic modifier BraA.CLF orchestrates flowering by modulating H3K27me3 levels at the floral integrator genes FT, SOC1, and SEP3, thereby influencing their expression. CURLY LEAF (CLF) is the catalytic subunit of the plant Polycomb Repressive Complex 2 that mediates the trimethylation of histone H3 lysine 27 (H3K27me3), an epigenetic modification that leads to gene silencing. While the function of CURLY LEAF (CLF) has been extensively studied in Arabidopsis thaliana, its role in Brassica crops is barely known. In this study, we focused on the Brassica rapa homolog of CLF and found that the loss-of-function mutant braA.clf-1 exhibits an accelerated flowering together with pleiotropic phenotypic alterations compared to wild-type plants. In addition, we carried out transcriptomic and H3K27me3 genome-wide analyses to identify the genes regulated by BraA.CLF. Interestingly, we observed that several floral regulatory genes, including the B. rapa homologs of FT, SOC1 and SEP3, show reduced H3K27me3 levels and increased transcript levels compared to wild-type plants, suggesting that they are direct targets of BraA.CLF and key players in regulating flowering time in this crop. In addition, the results obtained will enhance our understanding of the epigenetic mechanisms regulating key developmental traits and will aid to increase crop yield by engineering new Brassica varieties with different flowering time requirements.

Combined Transcriptomics and Metabolomics Analysis Reveals Profenofos-Induced Invisible Injury in Pakchoi (Brassica rapa L.) through Inhibition of Carotenoid Accumulation.

Profenofos insecticide poses risks to nontarget organisms including mammals and hydrobionts, and its effects on crops are not known. This study examined the invisible toxicity of profenofos on pakchoi (Brassica rapa L.), using transcriptome and metabolome analyses. Profenofos inhibited the photosynthetic efficiency and light energy absorption by leaves and severely damaged the chloroplasts, causing the accumulation of reactive oxygen species (ROS). Metabolomic analysis confirmed that profenofos promoted the conversion of ß-carotene into abscisic acid (ABA), as evidenced by the upregulation of the carotenoid biosynthesis pathway genes: zeaxanthin epoxidase (ZEP), 9-cis-epoxycarotenoid dioxygenase (NCED3), and xanthoxin dehydrogenase (XanDH). The inhibitory effects on carotenoid accumulation, photosynthesis, and increased ABA and ROS contents of the leaves led to invisible injury and stunted growth of the pakchoi plants. The findings of this study revealed the toxicological risk of profenofos to nontarget crops and provide guidance for the safe use of insecticides.

Assembly and comparative analysis of the complete mitochondrial genome of Brassica rapa var. Purpuraria.

BACKGROUND: Purple flowering stalk (Brassica rapa var. purpuraria) is a widely cultivated plant with high nutritional and medicinal value and exhibiting strong adaptability during growing. Mitochondrial (mt) play important role in plant cells for energy production, developing with an independent genetic system. Therefore, it is meaningful to assemble and annotate the functions for the mt genome of plants independently. Though there have been several reports referring the mt genome of in Brassica species, the genome of mt in B. rapa var. purpuraria and its functional gene variations when compared to its closely related species has not yet been addressed. RESULTS: The mt genome of B. rapa var. purpuraria was assembled through the Illumina and Nanopore sequencing platforms, which revealed a length of 219,775 bp with a typical circular structure. The base composition of the whole B. rapa var. purpuraria mt genome revealed A (27.45%), T (27.31%), C (22.91%), and G (22.32%). 59 functional genes, composing of 33 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes, were annotated. The sequence repeats, codon usage, RNA editing, nucleotide diversity and gene transfer between the cp genome and mt genome were examined in the B. rapa var. purpuraria mt genome. Phylogenetic analysis show that B. rapa var. Purpuraria was closely related to B. rapa subsp. Oleifera and B. juncea. Ka/Ks analysis reflected that most of the PCGs in the B. rapa var. Purpuraria were negatively selected, illustrating that those mt genes were conserved during evolution. CONCLUSIONS: The results of our findings provide valuable information on the B.rapa var. Purpuraria genome, which might facilitate molecular breeding, genetic variation and evolutionary researches for Brassica species in the future.

Genome-Wide Identification, Expression Analysis and Functional Study of DELLA Genes in Chinese Cabbage (Brassica rapa L. ssp. pekinensis).

BACKGROUND: DELLA protein is a crucial factor which played pivotal roles in regulating numerous intriguing biological processes in plant development and abiotic stress responses. However, little is known about the function and information of DELLA protein in Chinese cabbage. METHODS: Using 5 DELLA gene sequences in Arabidopsis Thaliana as probes, 5 DELLA genes in Chinese cabbage were identified by Blast search in Chinese cabbage database (Brassica database (BRAD)). The National Center for Biotechnology Information (NCBI), ExPaSy, SWISS-MODEL, DNAMAN, MEGA 11, PlantCARE were used to identify and analyze the DELLA gene family of Chinese cabbage. Gene expression was analyzed by quantitative real-time polymerase chain reaction (qRT-PCR). The function of BraA10gRGL3 was verified by overexpression and phenotypic analysis of BraA10gRGL3 and yeast hybrid. RESULTS: In this study, 5 BraDELLAs homologous to Arabidopsis thaliana were identified and cloned based on the Brassica database, namely, BraA02gRGL1, BraA05gRGL2, BraA10gRGL3, BraA06gRGA and BraA09gRGA. All BraDELLAs contain the DELLA, TVHYNP, and GRAS conserved domains. Cis-element analysis revealed that the promoter regions of these 5 DELLA genes all contain light-responsive elements, TCT motif, I-box, G-box, and box 4, which are associated with GA signaling. Transcriptome analysis results proved that the expression of BraA02gRGL1, BraA05gRGL2, and BraA10gRGL3 in Y2 at different growth stages were lower than them in Y7, which is consistent with the phenotype that Y7 exhibited stronger stress tolerance than Y2. It is worth emphasizing that even through the overexpression of BraA10gRGL3-Y7 in Arabidopsis resulted in smaller leaf size and lower fresh weight compared to the wild type (WT) Arabidopsis: Columbia, a stronger response to abiotic stresses was observed in BraA10gRGL3-Y7. It indicated that BraA10gRGL3-Y7 can improve the stress resistance of plants by inhibiting their growth. Moreover, the yeast two-hybrid experiment confirmed that BraA10gRGL3-Y7 can interacted with BraA05gGID1a-Y7, BraA04gGID1b1, BraA09gGID1b3-Y2, and BraA06gGID1c, whereas BraA10gRGL3-Y2 cannot interact with any BraGID1. CONCLUSIONS: Collectively, BraDELLAs play important role in plant development and response to abiotic stress. The differences in amino acid sequences between BraA10gRGL3-Y2 and BraA10gRGL3-Y7 may result in variations in their protein binding sites, thus affecting their interaction with the BraGID1 family proteins. This systematic analysis lays the foundation for further study of the functional characteristics of DELLA genes of Chinese cabbage.

Variations of Major Glucosinolates in Diverse Chinese Cabbage (Brassica rapa ssp. pekinensis) Germplasm as Analyzed by UPLC-ESI-MS/MS.

In this study, the variability of major glucosinolates in the leaf lamina of 134 Chinese cabbage accessions was investigated using Acquity ultra-performance liquid chromatography (UPLC-ESI-MS/MS). A total of twenty glucosinolates were profiled, of which glucobrassicanapin and gluconapin were identified as the predominant glucosinolates within the germplasm. These two glucosinolates had mean concentration levels above 1000.00 µmol/kg DW. Based on the principal component analysis, accessions IT186728, IT120044, IT221789, IT100417, IT278620, IT221754, and IT344740 were separated from the rest in the score plot. These accessions exhibited a higher content of total glucosinolates. Based on the VIP values, 13 compounds were identified as the most influential and responsible for variation in the germplasm. Sinigrin (r = 0.73), gluconapin (r = 0.78), glucobrassicanapin (r = 0.70), epiprogoitrin (r = 0.73), progoitrin (r = 0.74), and gluconasturtiin (r = 0.67) all exhibited a strong positive correlation with total glucosinolate at p < 0.001. This indicates that each of these compounds had a significant influence on the overall glucosinolate content of the various accessions. This study contributes valuable insights into the metabolic diversity of glucosinolates in Chinese cabbage, providing potential for breeding varieties tailored to consumer preferences and nutritional demands.

Identification and fine mapping of Brmmd1 gene controlling recessive genic male sterility in Brassica rapa L.

Recessive genic male sterility (RGMS) provides an effective approach for the commercial exploitation of heterosis, especially in Brassica crops. Although some artificial RGMS mutants have been reported in B. rapa, no causal genes derived from these natural mutants have been identified so far. In this study, a spontaneous RGMS mutant Bcajh97-01A derived from the 'Aijiaohuang' line traced back to the 1980 s was identified. Genetic analysis revealed that the RGMS trait was controlled by a single locus in the Bcajh97-01A/B system. Bulk segregant analysis (BSA) in combination with linkage analysis was employed to delimit the causal gene to an approximate 129 kb interval on chromosome A02. The integrated information of transcriptional levels and the predicted genes in the target region indicated that the Brmmd1 (BraA02g017420) encoding a PHD-containing nuclear protein was the most likely candidate gene. A 374 bp miniature inverted-repeat transposable element (MITE) was inserted into the first exon to prematurely stop the Brmmd1 gene translation, thus blocking the normal expression of this gene at the tetrad stage in the Bcajh97-01A. Additionally, a co-segregating structure variation (SV) marker was developed to rapidly screen the RGMS progenies from Bcajh97-01A/B system. Our findings reveal that BraA02g017420 is the causal gene responsible for the RGMS trait. This study lays a foundation for marker-assisted selection and further molecular mechanism exploration of pollen development in B. rapa.

Functional Differentiation of the Duplicated Gene BrrCIPK9 in Turnip (Brassica rapa var. rapa).

Gene duplication is a key biological process in the evolutionary history of plants and an important driving force for the diversification of genomic and genetic systems. Interactions between the calcium sensor calcineurin B-like protein (CBL) and its target, CBL-interacting protein kinase (CIPK), play important roles in the plant's response to various environmental stresses. As a food crop with important economic and research value, turnip (Brassica rapa var. rapa) has been well adapted to the environment of the Tibetan Plateau and become a traditional crop in the region. The BrrCIPK9 gene in turnip has not been characterized. In this study, two duplicated genes, BrrCIPK9.1 and BrrCIPK9.2, were screened from the turnip genome. Based on the phylogenetic analysis, BrrCIPK9.1 and BrrCIPK9.2 were found located in different sub-branches on the phylogenetic tree. Real-time fluorescence quantitative PCR analyses revealed their differential expression levels between the leaves and roots and in response to various stress treatments. The differences in their interactions with BrrCBLs were also revealed by yeast two-hybrid analyses. The results indicate that BrrCIPK9.1 and BrrCIPK9.2 have undergone Asparagine-alanine-phenylalanine (NAF) site divergence during turnip evolution, which has resulted in functional differences between them. Furthermore, BrrCIPK9.1 responded to high-pH (pH 8.5) stress, while BrrCIPK9.2 retained its ancestral function (low K+), thus providing further evidence of their functional divergence. These functional divergence genes facilitate turnip's good adaptation to the extreme environment of the Tibetan Plateau. In summary, the results of this study reveal the characteristics of the duplicated BrrCIPK9 genes and provide a basis for further functional studies of BrrCBLs-BrrCIPKs in turnip.

Deciphering the possible role of RNA-helicase genes mechanism in response to abiotic stresses in rapeseed (Brassica napus L.).

BACKGROUND: Plants mediate several defense mechanisms to withstand abiotic stresses. Several gene families respond to stress as well as multiple transcription factors to minimize abiotic stresses without minimizing their effects on performance potential. RNA helicase (RH) is one of the foremost critical gene families that can play an influential role in tolerating abiotic stresses in plants. However, little knowledge is present about this protein family in rapeseed (canola). Here, we performed a comprehensive survey analysis of the RH protein family in rapeseed (Brassica napus L.). RESULTS: A total of 133 BnRHs genes have been discovered in this study. By phylogenetic analysis, RHs genes were divided into one main group and a subgroup. Examination of the chromosomal position of the identified genes showed that most of the genes (27%) were located on chromosome 3. All 133 identified sequences contained the main DEXDC domain, the HELICC domain, and a number of sub-domains. The results of biological process studies showed that about 17% of the proteins acted as RHs, 22% as ATP binding, and 14% as mRNA binding. Each part of the conserved motifs, communication network, and three-dimensional structure of the proteins were examined separately. The results showed that the RWC in leaf tissue decreased with higher levels of drought stress and in both root and leaf tissues sodium concentration was increased upon increased levels of salt stress treatments. The proline content were found to be increased in leaf and root with the increased level of stress treatment. Finally, the expression patterns of eight selected RHs genes that have been exposed to drought, salinity, cold, heat and cadmium stresses were investigated by qPCR. The results showed the effect of genes under stress. Examination of gene expression in the Hayola #4815 cultivar showed that all primers except primer #79 had less expression in both leaves and roots than the control level. CONCLUSIONS: New finding from the study have been presented new insights for better understanding the function and possible mechanism of RH in response to abiotic stress in rapeseed.

The evolution and functional divergence of FT-related genes in controlling flowering time in Brassica rapa ssp. rapa.

KEY MESSAGE: The BrrFT paralogues exhibit distinct expression patterns and play different roles in regulating flowering time, and BrrFT4 competes with BrrFT1 and BrrFT2 to interact with BrrFD proteins. Flowering time is an important agricultural trait for Brassica crops, and early bolting strongly affects the yield and quality of Brassica rapa ssp. rapa. Flowering Locus T paralogues play an important role in regulating flowering time. In this study, we identified FT-related genes in turnip by phylogenetic classification, and four BrrFT homoeologs that shared with high identities with BraFT genes were isolated. The different gene structures, promoter binding sites, and expression patterns observed indicated that these genes may play different roles in flowering time regulation. Further genetic and biochemical experiments showed that as for FT-like paralogues, BrrFT2 acted as the key floral inducer, and BrrFT1 seems to act as a mild 'florigen' protein. However, BrrFT4 acts as a floral repressor and antagonistically regulates flowering time by competing with BrrFT1 and BrrFT2 to bind BrrFD proteins. BrrFT3 may have experienced loss of function via base shift mutation. Our results revealed the potential roles of FT-related genes in flowering time regulation in turnip.

Investigating the effect of drought stress and methanol spraying on the influential genes in the Calvin cycle and photorespiration of rapeseed (Brassica napus).

Due to global warming and changes in precipitation patterns, many regions are prone to permanent drought. Rapeseed (Brassica napus ) is one of the main sources of edible oils worldwide, and its production and yield are affected by drought. In this study, gene expression alterations under drought stress are investigated with bioinformatics studies to examine evolutionary relations of conserved motifs structure and interactions among Calvin cycle and photorespiration pathways key genes in drought-tolerant (SLM046) and drought-sensitive (Hayola308) genotypes of rapeseed. Investigating the conservation and evolutionary relationships revealed high conservation in motifs of FBPase, PRK, GlyK and NADP-ME enzymes. The analysis of protein interactions showed the correlation between FTRC, FBPase1, PRKX1, GlyKX2 and NADP-ME4 genes. Furthermore, in rapeseed, for the GlyKX2 and NADP-ME4 genes, four microRNAs of the miR172 family and four members of the miR167 family were identified as post-transcriptional regulators, respectively. The expression of ferredoxin thioredoxin reductase, fructose-1,6-bisphosphatase genes, phosphoribulokinase, glycerate kinase and malic enzyme 4 genes in the two rapeseed genotypes were evaluated by real-time qPCR method under 72h of drought stress and methanol foliar application. As a result, the highest expression levels of FTRC, PRKX1, GlyKX2, NADP-ME4 and FBPase1 were observed in methanol foliar application on the SLM046 genotype at 24h. In contrast, in methanol foliar application on the Hayola308 genotype, the highest expression levels of FTRC, PRKX1, GlyKX2, NADP-ME4 and FBPase1 were observed 8h after the treatment. Our study illustrated that methanol foliar application enhanced plant tolerance under drought stress.

Investigation of CaMV-host co-evolution through synonymous codon pattern.

Cauliflower mosaic virus (CaMV) has a double-stranded DNA genome and is globally distributed. The phylogeny tree of 121 CaMV isolates was categorized into two primary groups, with Iranian isolates showing the greatest genetic variations. Nucleotide A demonstrated the highest percentage (36.95%) in the CaMV genome and the dinucleotide odds ratio analysis revealed that TC dinucleotide (1.34 ≥ 1.23) and CG dinucleotide (0.63 ≤ 0.78) are overrepresented and underrepresented, respectively. Relative synonymous codon usage (RSCU) analysis confirmed codon usage bias in CaMV and its hosts. Brassica oleracea and Brassica rapa, among the susceptible hosts of CaMV, showed a codon adaptation index (CAI) value above 0.8. Additionally, relative codon deoptimization index (RCDI) results exhibited the highest degree of deoptimization in Raphanus sativus. These findings suggest that the genes of CaMV underwent codon adaptation with its hosts. Among the CaMV open reading frames (ORFs), genes that produce reverse transcriptase and virus coat proteins showed the highest CAI value of 0.83. These genes are crucial for the creation of new virion particles. The results confirm that CaMV co-evolved with its host to ensure the optimal expression of its genes in the hosts, allowing for easy infection and effective spread. To detect the force behind codon usage bias, an effective number of codons (ENC)-plot and neutrality plot were conducted. The results indicated that natural selection is the primary factor influencing CaMV codon usage bias.

Identification of Clubroot (Plasmodiophora brassicae) Resistance Loci in Chinese Cabbage (Brassica rapa ssp. pekinensis) with Recessive Character.

The soil-borne pathogen Plasmodiophora brassicae is the causal agent of clubroot, a major disease in Chinese cabbage (Brassica rapa ssp. pekinensis). The host's resistance genes often confer immunity to only specific pathotypes and may be rapidly overcome. Identification of novel clubroot resistance (CR) from germplasm sources is necessary. In this study, Bap246 was tested by being crossed with different highly susceptible B. rapa materials and showed recessive resistance to clubroot. An F2 population derived from Bap246 × Bac1344 was used to locate the resistance Quantitative Trait Loci (QTL) by Bulk Segregant Analysis Sequencing (BSA-Seq) and QTL mapping methods. Two QTL on chromosomes A01 (4.67-6.06 Mb) and A08 (10.42-11.43 Mb) were found and named Cr4Ba1.1 and Cr4Ba8.1, respectively. Fifteen and eleven SNP/InDel markers were used to narrow the target regions in the larger F2 population to 4.67-5.17 Mb (A01) and 10.70-10.84 Mb (A08), with 85 and 19 candidate genes, respectively. The phenotypic variation explained (PVE) of the two QTL were 30.97% and 8.65%, respectively. Combined with gene annotation, mutation site analysis, and real-time quantitative polymerase chain reaction (qRT-PCR) analysis, one candidate gene in A08 was identified, namely Bra020861. And an insertion and deletion (InDel) marker (co-segregated) named Crr1-196 was developed based on the gene sequence. Bra013275, Bra013299, Bra013336, Bra013339, Bra013341, and Bra013357 in A01 were the candidate genes that may confer clubroot resistance in Chinese cabbage. The resistance resource and the developed marker will be helpful in Brassica breeding programs.

Mutations in BrABCG26, encoding an ATP-binding cassette transporter, are responsible for male sterility in Chinese cabbage (Brassica rapa L. ssp. pekinensis).

KEY MESSAGE: The gene BrABCG26 responsible for male sterility of Chinese cabbage was confirmed by two allelic mutants. Male-sterile lines are an important way of heterosis utilization in Chinese cabbage. In this study, two allelic male-sterile mutants msm3-1 and msm3-2 were obtained from a Chinese cabbage double haploid (DH) line 'FT' by using EMS-mutagenesis. Compared to the wild-type 'FT,' the stamens of mutants were completely degenerated and had no pollen, and other characters had no obvious differences. Cytological observation revealed that the failure of vacuolation of the mononuclear microspore, accompanied by abnormal tapetal degradation, resulted in anther abortion in mutants. Genetic analysis showed that a recessive gene controlled the mutant trait. MutMap combined with kompetitive allele specific PCR genotyping analyses showed that BraA01g038270.3C, encoding a transporter ABCG26 that played a vital role in pollen wall formation, was the candidate gene for msm3-1, named BrABCG26. Compared with wild-type 'FT,' the mutations existed on the second exon (C to T) and the sixth exon (C to T) of BrABCG26 gene in mutants msm3-1 and msm3-2, leading to the loss-of-function truncated protein, which verified the BrABCG26 function in stamen development. Subcellular localization and expression pattern analysis indicated that BrABCG26 was localized in the nucleus and was expressed in all organs, with the highest expression in flower buds. Compared to the wild-type 'FT,' the expressions of BrABCG26 were significantly reduced in flower buds and anthers of mutants. Promoter activity analysis showed that a strong GUS signal was detected in flower buds. These results indicated that BrABCG26 is responsible for the male sterility of msm3 mutants in Chinese cabbage.

Morpho-Physiochemical Indices and Transcriptome Analysis Reveal the Role of Glucosinolate and Erucic Acid in Response to Drought Stress during Seed Germination of Rapeseed.

The global expansion of rapeseed seed quality has been focused on maintaining glucosinolate (GSL) and erucic acid (EA) contents. However, the influence of seed GSL and EA contents on the germination process under drought stress remains poorly understood. Herein, 114 rapeseed accessions were divided into four groups based on GSL and EA contents to investigate their performance during seed imbibition under drought stress. Our results revealed significant variations in seed germination-related traits, particularly with higher GSL and EA, which exhibited higher germination % (G%) and lower mean germination time (MGT) under drought stress conditions. Moreover, osmoregulation, enzymatic system and hormonal regulation were improved in high GSL and high EA (HGHE) versus low GSL and low EA (LGLE) seeds, indicating the essential protective role of GSL and EA during the germination process in response to drought stress. The transcriptional regulation mechanism for coordinating GSL-EA-related pathways in response to drought stress during seed imbibition was found to involve the differential expression of sugar metabolism-, antioxidant-, and hormone-related genes with higher enrichment in HGHE compared to LGLE seeds. GO enrichment analysis showed higher variations in transcription regulator activity and DNA-binding transcription factors, as well as ATP and microtubule motor activity in GSL-EA-related pathways. Furthermore, KEGG analysis identified cellular processes, environmental information processing, and metabolism categories, with varied gene participation between GSL, EA and GSL-EA-related pathways. For further clarification, QY7 (LGLE) seeds were primed with different concentrations of GSL and EA under drought stress conditions. The results showed that 200 µmol/L of GSL and 400 µmol/L of EA significantly improved G%, MGT, and seedling fresh weight, besides regulating stress and fatty acid responsive genes during the seed germination process under drought stress conditions. Conclusively, exogenous application of GSL and EA is considered a promising method for enhancing the drought tolerance of LGLE seeds. Furthermore, the current investigation could provide a theoretical basis of GSL and EA roles and their underlying mechanisms in stress tolerance during the germination process.

Homoeologous exchanges contribute to branch angle variations in rapeseed: Insights from transcriptome, QTL-seq and gene functional analysis.

Branch angle (BA) is a critical morphological trait that significantly influences planting density, light interception and ultimately yield in plants. Despite its importance, the regulatory mechanism governing BA in rapeseed remains poorly understood. In this study, we generated 109 transcriptome data sets for 37 rapeseed accessions with divergent BA phenotypes. Relative to adaxial branch segments, abaxial segments accumulated higher levels of auxin and exhibited lower expression of six TCP1 homologues and one GA20ox3. A co-expression network analysis identified two modules highly correlated with BA. The modules contained homologues to known BA control genes, such as FUL, YUCCA6, TCP1 and SGR3. Notably, a homoeologous exchange (HE), occurring at the telomeres of A09, was prevalent in large BA accessions, while an A02-C02 HE was common in small BA accessions. In their corresponding regions, these HEs explained the formation of hub gene hotspots in the two modules. QTL-seq analysis confirmed that the presence of a large A07-C06 HE (~8.1 Mb) was also associated with a small BA phenotype, and BnaA07.WRKY40.b within it was predicted as candidate gene. Overexpressing BnaA07.WRKY40.b in rapeseed increased BA by up to 20°, while RNAi- and CRISPR-mediated mutants (BnaA07.WRKY40.b and BnaC06.WRKY40.b) exhibited decreased BA by up to 11.4°. BnaA07.WRKY40.b was exclusively localized to the nucleus and exhibited strong expression correlations with many genes related to gravitropism and plant architecture. Taken together, our study highlights the influence of HEs on rapeseed plant architecture and confirms the role of WRKY40 homologues as novel regulators of BA.

Two aquaporins, PIP1;1 and PIP2;1, mediate the uptake of neonicotinoid pesticides in plants.

Neonicotinoids (NEOs), a large class of organic compounds, are a type of commonly used pesticide for crop protection. Their uptake and accumulation in plants are prerequisites for their intra- and intercellular movements, transformation, and function. Understanding the molecular mechanisms that underpin NEO uptake by plants is crucial for effective application, which remains elusive. Here, we demonstrate that NEOs enter plant cells primarily through the transmembrane symplastic pathway and accumulate mainly in the cytosol. Two plasma membrane intrinsic proteins discovered in Brassica rapa, BraPIP1;1 and BraPIP2;1, were found to encode aquaporins (AQPs) that are highly permeable to NEOs in different plant species and facilitate NEO subcellular diffusion and accumulation. Their conserved transport function was further demonstrated in Xenopus laevis oocyte and yeast assays. BraPIP1;1 and BraPIP2;1 gene knockouts and interaction assays suggested that their proteins can form functional heterotetramers. Assessment of the potential of mean force indicated a negative correlation between NEO uptake and the energy barrier of BraPIP1;1 channels. This study shows that AQPs transport organic compounds with greater osmolarity than previously thought, providing new insight into the molecular mechanisms of organic compound uptake and facilitating innovations in systemic pesticides.

Identification of Yellow Seed Color Genes Using Bulked Segregant RNA Sequencing in Brassica juncea L.

Yellow seed breeding is an effective method to improve oil yield and quality in rapeseed (Brassica napus L.). However, naturally occurring yellow-seeded genotypes have not been identified in B. napus. Mustard (Brassica juncea L.) has some natural, yellow-seeded germplasms, yet the molecular mechanism underlying this trait remains unclear. In this study, a BC9 population derived from the cross of yellow seed mustard "Wuqi" and brown seed mustard "Wugong" was used to analyze the candidate genes controlling the yellow seed color of B. juncea. Subsequently, yellow-seeded (BY) and brown-seeded (BB) bulks were constructed in the BC9 population and subjected to bulked segregant RNA sequencing (BSR-Seq). A total of 511 differentially expressed genes (DEGs) were identified between the brown and yellow seed bulks. Enrichment analysis revealed that these DEGs were involved in the phenylpropanoid biosynthetic process and flavonoid biosynthetic process, including key genes such as 4CL, C4H, LDOX/TT18, PAL1, PAL2, PAL4, TT10, TT12, TT4, TT8, BAN, DFR/TT3, F3H/TT6, TT19, and CHI/TT5. In addition, 111,540 credible single-nucleotide polymorphisms (SNPs) and 86,319 INDELs were obtained and used for quantitative trait locus (QTL) identification. Subsequently, two significant QTLs on chromosome A09, namely, qSCA09-3 and qSCA09-7, were identified by G' analysis, and five DEGs (BjuA09PAL2, BjuA09TT5, BjuA09TT6, BjuA09TT4, BjuA09TT3) involved in the flavonoid pathway were identified as hub genes based on the protein-to-protein network. Among these five genes, only BjuA09PAL2 and BjuA09F3H had SNPs between BY and BB bulks. Interestingly, the majority of SNPs in BjuA09PAL2 were consistent with the SNPs identified between the high-quality assembled B. juncea reference genome "T84-66" (brown-seed) and "AU213" (yellow-seed). Therefore, BjuA09PAL2, which encodes phenylalanine lyase, was considered as the candidate gene associated with yellow seed color of B. juncea. The identification of a novel gene associated with the yellow seed coloration of B. juncea through this study may play a significant role in enhancing yellow seed breeding in rapeseed.