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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Targeted mutagenesis of BnTTG1 homologues generated yellow-seeded rapeseed with increased oil content and seed germination under abiotic stress.

Yellow seed is one desirable trait with great potential to improve seed oil quality and yield. The present study surveys the redundant role of BnTTG1 genes in the proanthocyanidins (PA) biosynthesis, oil content and abiotic stress resistance. Stable yellow seed mutants were generated after mutating BnTTG1 by CRISPR/Cas9 genome editing system. Yellow seed phenotype could be obtained only when both functional homologues of BnTTG1 were simultaneously knocked out. Homozygous mutants of BnTTG1 homologues showed decreased thickness and PA accumulation in seed coat. Transcriptome and qRT-PCR analysis indicated that BnTTG1 mutation inhibited the expression of genes involved in phenylpropanoid and flavonoid biosynthetic pathways. Increased seed oil content and alteration of fatty acid (FA) composition were observed in homozygous mutants of BnTTG1 with enriched expression of genes involved in FA biosynthesis pathway. In addition, target mutation of BnTTG1 accelerated seed germination rate under salt and cold stresses. Enhanced seed germination capacity in BnTTG1 mutants was correlated with the change of expression level of ABA responsive genes. Overall, this study elucidated the redundant role of BnTTG1 in regulating seed coat color and established an efficient approach for generating yellow-seeded oilseed rape genetic resources with increase oil content, modified FA composition and resistance to multiple abiotic stresses.

Genome-wide characterization of SDR gene family and its potential role in seed dormancy of Brassica napus L.

Rapeseed (Brassica napus L.) with short or no dormancy period are easy to germinate before harvest (pre-harvest sprouting, PHS). PHS has seriously decreased seed weight and oil content in B. napus. Short-chain dehydrogenase/ reductase (SDR) genes have been found to related to seed dormancy by promoting ABA biosynthesis in rice and Arabidopsis. In order to clarify whether SDR genes are the key factor of seed dormancy in B. napus, homology sequence blast, protein physicochemical properties, conserved motif, gene structure, cis-acting element, gene expression and variation analysis were conducted in present study. Results shown that 142 BnaSDR genes, unevenly distributed on 19 chromosomes, have been identified in B. napus genome. Among them, four BnaSDR gene clusters present in chromosome A04、A05、C03、C04 were also identified. These 142 BnaSDR genes were divided into four subfamilies on phylogenetic tree. Members of the same subgroup have similar protein characters, conserved motifs, gene structure, cis-acting elements and tissue expression profiles. Specially, the expression levels of genes in subgroup A, B and C were gradually decreased, but increased in subgroup D with the development of seeds. Among seven higher expressed genes in group D, six BnaSDR genes were significantly higher expressed in weak dormancy line than that in nondormancy line. And the significant effects of BnaC01T0313900ZS and BnaC03T0300500ZS variation on seed dormancy were also demonstrated in present study. These findings provide a key information for investigating the function of BnaSDRs on seed dormancy in B. napus.

Simulation-based establishment of base pools for a hybrid breeding program in winter rapeseed.

KEY MESSAGE: Simulation planned pre-breeding can increase the efficiency of starting a hybrid breeding program. Starting a hybrid breeding program commonly comprises a grouping of the initial germplasm in two pools and subsequent selection on general combining ability. Investigations on pre-breeding steps before starting the selection on general combining ability are not available. Our goals were (1) to use computer simulations on the basis of DNA markers and testcross data to plan crosses that separate genetically two initial germplasm pools of rapeseed, (2) to carry out the planned crosses, and (3) to verify experimentally the pool separation as well as the increase in testcross performance. We designed a crossing program consisting of four cycles of recombination. In each cycle, the experimentally generated material was used to plan the subsequent crossing cycle with computer simulations. After finishing the crossing program, the initially overlapping pools were clearly separated in principal coordinate plots. Doubled haploid lines derived from the material of crossing cycles 1 and 2 showed an increase in relative testcross performance for yield of about 5% per cycle. We conclude that simulation-designed pre-breeding crossing schemes, that were carried out before the general combining ability-based selection of a newly started hybrid breeding program, can save time and resources, and in addition conserve more of the initial genetic variation than a direct start of a hybrid breeding program with general combining ability-based selection.

Rapid genomic evolution in Brassica rapa with bumblebee selection in experimental evolution.

BACKGROUND: Insect pollinators shape rapid phenotypic evolution of traits related to floral attractiveness and plant reproductive success. However, the underlying genomic changes remain largely unknown despite their importance in predicting adaptive responses to natural or to artificial selection. Based on a nine-generation experimental evolution study with fast cycling Brassica rapa plants adapting to bumblebees, we investigate the genomic evolution associated with the previously observed parallel phenotypic evolution. In this current evolve and resequencing (E&R) study, we conduct a genomic scan of the allele frequency changes along the genome in bumblebee-pollinated and hand-pollinated plants and perform a genomic principal component analysis (PCA). RESULTS: We highlight rapid genomic evolution associated with the observed phenotypic evolution mediated by bumblebees. Controlling for genetic drift, we observe significant changes in allelic frequencies at multiple loci. However, this pattern differs according to the replicate of bumblebee-pollinated plants, suggesting putative non-parallel genomic evolution. Finally, our study underlines an increase in genomic variance implying the putative involvement of multiple loci in short-term pollinator adaptation. CONCLUSIONS: Overall, our study enhances our understanding of the complex interactions between pollinator and plants, providing a stepping stone towards unravelling the genetic basis of plant genomic adaptation to biotic factors in the environment.