Comparative transcriptome analysis reveals transcriptional regulation of anthocyanin biosynthesis in purple radish (Raphanus sativus L.).

Radish exhibits significant variation in color, particularly in sprouts, leaves, petals, fleshy roots, and other tissues, displaying a range of hues such as green, white, red, purple, and black. Although extensive research has been conducted on the color variation of radish, the underlying mechanism behind the variation in radish flower color remains unclear. To date, there is a lack of comprehensive research investigating the variation mechanism of radish sprouts, leaves, fleshy roots, and flower organs. This study aims to address this gap by utilizing transcriptome sequencing to acquire transcriptome data for white and purple radish flowers. Additionally, the published transcriptome data of sprouts, leaves, and fleshy roots were incorporated to conduct a systematic analysis of the regulatory mechanisms underlying anthocyanin biosynthesis in these four radish tissues. The comparative transcriptome analysis revealed differential expression of the anthocyanin biosynthetic pathway genes DFR, UGT78D2, TT12 and CPC in the four radish tissues. Additionally, the WGCNA results identified RsDFR.9c and RsUGT78D2.2c as hub genes responsible for regulating anthocyanin biosynthesis. By integrating the findings from the comparative transcriptome analysis, WGCNA, and anthocyanin biosynthetic pathway-related gene expression patterns, it is hypothesized that genes RsDFR.9c and RsUGT78D2.2c may serve as pivotal regulators of anthocyanins in the four radish tissues. Furthermore, the tissue-specific expression of the four copies of RsPAP1 is deemed crucial in governing anthocyanin synthesis and accumulation. Our results provide new insights into the molecular mechanism of anthocyanin biosynthesis and accumulation in different tissues of radish.

Metabolic and transcriptome analysis of dark red taproot in radish (Raphanus sativus L.).

The red color in radish taproots is an important quality index and is mainly affected by anthocyanins. However, the metabolite components and gene expression underlying dark red taproot color formation in radish remain elusive. In this study, the metabolites and gene expression patterns affecting anthocyanin biosynthesis were monitored in the dark red taproots. Comparative analysis of anthocyanin metabolites between dark red taproots and white taproots indicated that pelargonin and pelargonidin 3-O-beta-D-glucoside were the most promising dark red pigments responsible for the coloration of the taproots. Transcriptomic analysis of gene expression between dark red taproots and white taproots revealed that most of genes involved in the anthocyanin biosynthesis pathway were up-regulated in dark red taproots. In particular, RsCHS and RsDFR were the two most up-regulated genes in the dark red taproots. Moreover, the higher coexpression of two R2R3-Myb transcription factors, RsMYB1 and RsMYB2, may contribute to dark red color formation. Our work documents metabolomic and transcriptomic changes related to the dark red color formation in taproots radish and provides valuable data for anthocyanin-rich radish breeding.

Genetic and Molecular Characterization of a Self-Compatible Brassica rapa Line Possessing a New Class II S Haplotype.

Most flowering plants have evolved a self-incompatibility (SI) system to maintain genetic diversity by preventing self-pollination. The Brassica species possesses sporophytic self-incompatibility (SSI), which is controlled by the pollen- and stigma-determinant factors SP11/SCR and SRK. However, the mysterious molecular mechanism of SI remains largely unknown. Here, a new class II S haplotype, named BrS-325, was identified in a pak choi line '325', which was responsible for the completely self-compatible phenotype. To obtain the entire S locus sequences, a complete pak choi genome was gained through Nanopore sequencing and de novo assembly, which provided a good reference genome for breeding and molecular research in B. rapa. S locus comparative analysis showed that the closest relatives to BrS-325 was BrS-60, and high sequence polymorphism existed in the S locus. Meanwhile, two duplicated SRKs (BrSRK-325a and BrSRK-325b) were distributed in the BrS-325 locus with opposite transcription directions. BrSRK-325b and BrSCR-325 were expressed normally at the transcriptional level. The multiple sequence alignment of SCRs and SRKs in class II S haplotypes showed that a number of amino acid variations were present in the contact regions (CR II and CR III) of BrSCR-325 and the hypervariable regions (HV I and HV II) of BrSRK-325s, which may influence the binding and interaction between the ligand and the receptor. Thus, these results suggested that amino acid variations in contact sites may lead to the SI destruction of a new class II S haplotype BrS-325 in B. rapa. The complete SC phenotype of '325' showed the potential for practical breeding application value in B. rapa.

Characterization of a Common S Haplotype BnS-6 in the Self-Incompatibility of Brassica napus.

Self-incompatibility (SI) is a pollen-stigma recognition system controlled by a single and highly polymorphic genetic locus known as the S-locus. The S-locus exists in all Brassica napus (B. napus, AACC), but natural B. napus accessions are self-compatible. About 100 and 50 S haplotypes exist in Brassica rapa (AA) and Brassica oleracea (CC), respectively. However, S haplotypes have not been detected in B. napus populations. In this study, we detected the S haplotype distribution in B. napus and ascertained the function of a common S haplotype BnS-6 through genetic transformation. BnS-1/BnS-6 and BnS-7/BnS-6 were the main S haplotypes in 523 B. napus cultivars and inbred lines. The expression of SRK in different S haplotypes was normal (the expression of SCR in the A subgenome affected the SI phenotype) while the expression of BnSCR-6 in the C subgenome had no correlation with the SI phenotype in B. napus. The BnSCR-6 protein in BnSCR-6 overexpressed lines was functional, but the self-compatibility of overexpressed lines did not change. The low expression of BnSCR-6 could be a reason for the inactivation of BnS-6 in the SI response of B. napus. This study lays a foundation for research on the self-compatibility mechanism and the SI-related breeding in B. napus.

A novel negative-stranded RNA virus of the order Bunyavirales identified in Brassica campestris L. ssp. chinensis.

Here, we report the full-length genome sequence of a novel cogu-like virus identified in Brassica campestris L. ssp. Chinensis (B. campestris), an economically important vegetable in China. This virus, tentatively named "Brassica campestris chinensis coguvirus 1" (BCCoV1), has a bipartite genome that consists of two RNA molecules (RNA1 and RNA2). The negative-stranded (ns) RNA1 is 6757 nt in length, encoding the putative RNA-dependent RNA polymerase (RdRp), and the ambisense RNA2 is 3061 nt long, encoding the putative movement protein (MP) and nucleocapsid protein (NP). A homology search of the RdRp, MP, and NP showed that they are closely related to five other recently discovered negative-stranded RNA (nsRNA) viruses infecting plants, belonging to the new genus Coguvirus. Phylogenetic analysis of the 252-kDa RdRp confirmed the classification of this virus, showing that BCCoV1 possibly belongs to the genus Coguvirus, family Phenuiviridae, order Bunyavirales. The present study improves our understanding of the viral diversity in B. campestris and the evolution of nsRNA viruses.

Complete nucleotide sequence of a novel partitivirus from Brassica campestris L. ssp. chinensis.

In the present work, we report the discovery and complete genome sequence of a novel partitivirus identified from Brassica campestris L. ssp. chinensis, which we have named "Brassica campestris chinensis cryptic virus 1" (BCCV1). Next-generation sequencing (NGS) combined with adapter-ligation-mediated amplification allowed assembly of the full-length genome sequence of BCCV1. The genome of BCCV1 contains two dsRNA segments, dsRNA1 (1595 bp) and dsRNA2 (1591 bp), which encode a conserved RNA-dependent RNA polymerase (RdRp) and a putative capsid protein (CP), respectively. Homology searches and phylogenetic analysis of the 479-aa RdRp and 438-aa CP showed that BCCV1 is a new member of the genus Deltapartitivirus, family Partitiviridae. This is the first report of the identification of a member of the family Partitiviridae in Brassica campestris L. ssp. chinensis.

Comparative transcript profiling and cytological observation of the newly bred recessive genic male sterility non-heading Chinese cabbage (Brassica rapa ssp. chinensis) line WS24-3A.

BACKGROUND: WS24-3A is a newly bred non-heading Chinese cabbage genic male-sterile line, in which sterility is controlled by a recessive gene, designated as Bra2ms. WS24-3A has been used for hybrid breeding. OBJECTIVE: To reveal the underlying molecular mechanisms responsible for the sterility of WS24-3A. METHODS: Cytological observation of the process of sterile/fertile anther development was performed to determine the tissue and stage in which sterility occurs. Phenotyping and transcriptomic analyses were performed to identify differentially expressed genes (DEGs) between sterile and fertile flower buds at different stages. RESULTS: Cytological analysis revealed no tetrads at stage 7 or at later stages of anther development, and the degradation of callose was delayed. Abnormal meiocytes were surrounded by sustaining callose that degenerated gradually in WS24-3A. Comparative transcript profiling identified 3282 DEGs during three anther developmental stages, namely, pre-meiotic anther, meiotic anther, and anthers with single-celled pollen stage. The difference in DEG percentage between up-regulated and down-regulated at meiotic anther stage was obviously larger than at the other two stages; further, most DEGs are important for male meiosis, callose synthesis and dissolution, and tapetum development. Ten DEGs were found to be involved in anther and pollen development, which were analyzed by quantitative PCR. CONCLUSION: Bra2ms affected gene expression in meiocytes and associated with callose synthesis, degradation and tapetum development. Our results provide clues to elucidate the molecular mechanism of genic male sterility in non-heading Chinese cabbage.

Genome-wide analysis of long non-coding RNAs unveils the regulatory roles in the heat tolerance of Chinese cabbage (Brassica rapa ssp.chinensis).

Long non-coding RNAs (lncRNAs) mediate important epigenetic regulation in various biological processes related to the stress response in plants. However, the systematic analysis of the lncRNAs expressed in Brassica rapa under heat stress has been elusive. In this study, we performed a genome-wide analysis of the lncRNA expression profiles in non-heading Chinese cabbage leaves using strand-specific RNA-sequencing. A total of 4594 putative lncRNAs were identified with a comprehensive landscape of dynamic lncRNA expression networks under heat stress. Co-expression networks of the interactions among the differentially expressed lncRNAs, mRNAs and microRNAs revealed that several phytohormones were associated with heat tolerance, including salicylic acid (SA) and brassinosteroid (BR) pathways. Of particular importance is the discovery of 25 lncRNAs that were highly co-expressed with 10 heat responsive genes. Thirty-nine lncRNAs were predicted as endogenous target mimics (eTMs) for 35 miRNAs, and five of them were validated to be involved in the heat tolerance of Chinese cabbage. Heat responsive lncRNA (TCONS_00048391) is an eTM for bra-miR164a, that could be a sponge for miRNA binding and may be a competing endogenous RNA (ceRNA) for the target gene NAC1 (Bra030820), affecting the expression of bra-miR164a in Chinese cabbage. Thus, these findings provide new insights into the functions of lncRNAs in heat tolerance and highlight a set of candidate lncRNAs for further studies in non-heading Chinese cabbage.

Time-Course Transcriptome Analysis of Compatible and Incompatible Pollen-Stigma Interactions in Brassica napus L.

Brassica species exhibit both compatible and incompatible pollen-stigma interactions, however, the underlying molecular mechanisms remain largely unknown. Here, RNA-seq technology was applied in a comprehensive time-course experiment (2, 5, 10, 20, and 30 min) to explore gene expression during compatible/incompatible pollen-stigma interactions in stigma. Moderate changes of gene expression were observed both in compatible pollination (PC) and incompatible pollination (PI) within 10 min, whereas drastic changes showed up by 30 min, especially in PI. Stage specific DEGs [Differentially Expressed Gene(s)] were identified, and signaling pathways such as stress response, defense response, cell wall modification and others were found to be over-represented. In addition, enriched genes in all samples were analyzed as well, 293 most highly expressed genes were identified and annotated. Gene Ontology and metabolic pathway analysis revealed 10 most highly expressed genes and 37 activated metabolic pathways. According to the data, downstream components were activated in signaling pathways of both compatible and incompatible responses, and incompatible response had more complicated signal transduction networks. This study provides more detailed molecular information at different time points after compatible and incompatible pollination, deepening our knowledge about pollen-stigma interactions.

Helitron-like transposons contributed to the mating system transition from out-crossing to self-fertilizing in polyploid Brassica napus L.

The mating system transition in polyploid Brassica napus (AACC) from out-crossing to selfing is a typical trait to differentiate it from their diploid progenitors. Elucidating the mechanism of mating system transition has profound consequences for understanding the speciation and evolution in B. napus. Functional complementation experiment has shown that the insertion of 3.6 kb into the promoter of self-incompatibility male determining gene, BnSP11-1 leads to its loss of function in B. napus. The inserted fragment was found to be a non-autonomous Helitron transposon. Further analysis showed that the inserted 3.6 kb non-autonomous Helitron transposon was widely distributed in B. napus accessions which contain the S haplotype BnS-1. Through promoter deletion analysis, an enhancer and a putative cis-regulatory element (TTCTA) that were required for spatio-temporal specific expression of BnSP11-1 were identified, and both might be disrupted by the insertion of Helitron transposon. We suggested that the insertion of Helitron transposons in the promoter of BnSP11-1 gene had altered the mating system and might facilitated the speciation of B. napus. Our findings have profound consequences for understanding the self-compatibility in B. napus as well as for the trait variations during evolutionary process of plant polyploidization.

Tribenuron-Methyl Induces Male Sterility through Anther-Specific Inhibition of Acetolactate Synthase Leading to Autophagic Cell Death.

Tribenuron-methyl (TM) is a powerful sulfonylurea herbicide that inhibits branched-chain amino acid (BCAA) biosynthesis by targeting the catalytic subunit (CSR1) of acetolactate synthase (ALS). Selective induction of male sterility by foliar spraying of TM at low doses has been widely used for hybrid seed production in rapeseed (Brassica napus); however, the underlying mechanism remains unknown. Here, we report greater TM accumulation and subsequent stronger ALS inhibition and BCAA starvation in anthers than in leaves and stems after TM application. Constitutive or anther-specific expression of csr1-1D (a CSR1 mutant) eliminated anther-selective ALS inhibition and reversed the TM-induced male sterile phenotype in both rapeseed and Arabidopsis. The results of TM daub-stem experiments, combined with the observations of little TM accumulation in anthers and reversion of TM-induced male sterility by targeted expression of the TM metabolism gene Bel in either the mesophyll or phloem, suggested that foliar-sprayed TM was polar-transported to anthers mainly through the mesophyll and phloem. Microscopy and immunoblotting revealed that autophagy, a bulk degradation process induced during cell death, was elevated in TM-induced male sterile anthers and by anther-specific knockdown of ALS. Moreover, TM-induced pollen abortion was significantly inhibited by the autophagy inhibitor 3-MA. These data suggested that TM was polar-transported to anthers, resulting in BCAA starvation via anther-specific ALS inhibition and, ultimately, autophagic cell death in anthers.

Gene expression and genetic analysis reveal diverse causes of recessive self-compatibility in Brassica napus L.

BACKGROUND: Brassica napus (AACC) is self-compatible, although its ancestor species Brassica rapa (AA) and Brassica oleracea (CC) are self-incompatible. Most B.napus accessions have dominant self-compatibility (SC) resulting from an insertion of 3.6 kb in the promoter region of BnSCR-1 on the A genome, while recessive SC in B.napus has rarely been observed. Expression and cloning of SRK and SCR genes and genetic analysis were carried out to dissect bases of recessive SC in B.napus. RESULTS: Eleven accessions were screened to identify stable recessive SC and had the S genotype BnS-7 on the A genome and BnS-6 on the C genome similarly to BrS-29 and BoS-15, respectively. In eight SC accessions, BnSCR-7 and BnSCR-6 were nearly undetectable and harbored no structural mutations in the promoters, while SRK genes were expressed at normal levels and contained intact CDS, with the exception of BnSRK-7 in line C32. SRK and SCR genes were expressed normally but their CDSs had no mutations in three SC accessions. In self-incompatible S-1300 and 11 F1 hybrids, SRK genes and BnSCR-1300 transcripts were present at high levels, while expression of the BnSCR-7 and BnSCR-6 were absent. Plants of S genotype S1300S1300 were completely SI, while SI phenotypes of SBnS-7SBnS-7 and S1300SBnS-7 plants were segregated in BC1 and F2 populations. CONCLUSIONS: The recessive SC in eight accessions is caused by the loss of function of BnSCR-7 and BnSCR-6 in pollen. Translational repression contributes to the recessive SC in three accessions, whose SRK and SCR genes were expressed normally and had identical CDSs to BrS-29 or BoS-15. SI in 11 F1 hybrids relies on the expression of BnSCR-1300 rather than SRK genes. Other factor(s) independent of the S locus are involved in recessive SC. Therefore, diverse causes underlie recessive SC in B. napus, yielding insight into these complex mechanisms.

Characterization of Sucrose transporter alleles and their association with seed yield-related traits in Brassica napus L.

BACKGROUND: Sucrose is the primary photosynthesis product and the principal translocating form within higher plants. Sucrose transporters (SUC/SUT) play a critical role in phloem loading and unloading. Photoassimilate transport is a major limiting factor for seed yield. Our previous research demonstrated that SUT co-localizes with yield-related quantitative trait loci. This paper reports the isolation of BnA7.SUT1 alleles and their promoters and their association with yield-related traits. RESULTS: Two novel BnA7.SUT1 genes were isolated from B. napus lines 'Eagle' and 'S-1300' and designated as BnA7.SUT1.a and BnA7.SUT1.b, respectively. The BnA7.SUT1 protein exhibited typical SUT features and showed high amino acid homology with related species. Promoters of BnA7.SUT1.a and BnA7.SUT1.b were also isolated and classified as pBnA7.SUT1.a and pBnA7.SUT1.b, respectively. Four dominant sequence-characterized amplified region markers were developed to distinguish BnA7.SUT1.a and BnA7.SUT1.b. The two genes were estimated as alleles with two segregating populations (F2 and BC1) obtained by crossing '3715'×'3769'. BnA7.SUT1 was mapped to the A7 linkage group of the TN doubled haploid population. In silico analysis of 55 segmental BnA7.SUT1 alleles resulted three BnA7.SUT1 clusters: pBnA7.SUT1.a- BnA7.SUT1.a (type I), pBnA7.SUT1.b- BnA7.SUT1.a (type II), and pBnA7.SUT1.b- BnA7.SUT1.b (type III). Association analysis with a diverse panel of 55 rapeseed lines identified single nucleotide polymorphisms (SNPs) in promoter and coding domain sequences of BnA7.SUT1 that were significantly associated with one of three yield-related traits: number of effective first branches (EFB), siliques per plant (SP), and seed weight (n = 1000) (TSW) across all four environments examined. SNPs at other BnA7.SUT1 sites were also significantly associated with at least one of six yield-related traits: EFB, SP, number of seeds per silique, seed yield per plant, block yield, and TSW. Expression levels varied over various tissue/organs at the seed-filling stage, and BnA7.SUT1 expression positively correlated with EFB and TSW. CONCLUSIONS: Sequence, mapping, association, and expression analyses collectively showed significant diversity between the two BnA7.SUT1 alleles, which control some of the phenotypic variation for branch number and seed weight in B. napus consistent with expression levels. The associations between allelic variation and yield-related traits may facilitate selection of better genotypes in breeding.