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.

Too hot to handle: temperature-induced plasticity influences pollinator behaviour and plant fitness.

Increased temperature can induce plastic changes in many plant traits. However, little is known about how these changes affect plant interactions with insect pollinators and herbivores, and what the consequences for plant fitness and selection are. We grew fast-cycling Brassica rapa plants at two temperatures (ambient and increased temperature) and phenotyped them (floral traits, scent, colour and glucosinolates). We then exposed plants to both pollinators (Bombus terrestris) and pollinating herbivores (Pieris rapae). We measured flower visitation, oviposition of P. rapae, herbivore development and seed output. Plants in the hot environment produced more but smaller flowers, with lower UV reflectance and emitted a different volatile blend with overall lower volatile emission. Moreover, these plants received fewer first-choice visits by bumblebees and butterflies, and fewer flower visits by butterflies. Seed production was lower in hot environment plants, both because of a reduction in flower fertility due to temperature and because of the reduced visitation of pollinators. The selection on plant traits changed in strength and direction between temperatures. Our study highlights an important mechanism by which global warming can change plant-pollinator interactions and negatively impact plant fitness, as well as potentially alter plant evolution through changes in phenotypic selection.

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.

Optimizing oilseed rape growth: Exploring the effect of foliar biostimulants on the interplay among metabolism, phenology, and yield.

The current agricultural system is in search of new strategies to achieve a more sustainable production while keeping or even increasing crop yield and quality. In this scenario, the application of biostimulants constitutes a potent solution. In the current study, the impact of a blue-green microalgal extract (MB) and a pig tissue hydrolysate (PTH) on rapeseed plants' development was characterized. Obtained results revealed a positive effect on yield parameters of plants treated with MB and, especially, PTH; this was associated to an improvement on the photosynthetic performance. Moreover, this study remarked the effects of biostimulants on plant phenology through their pivotal role in modulating developmental processes. More specifically, proteomic, metabolomic, and hormone content analyses revealed distinct alterations associated with the acceleration of phenology induced by biostimulant application. Additionally, some antioxidant enzymes and stress-related compounds were up-regulated upon MB and PTH treatments, indicating enhanced plant defense mechanisms in response to accelerated phenological transitions. Such findings highlight the intricate interplay between biostimulants and plant physiology, wherein biostimulants orchestrate rapid developmental changes, ultimately influencing growth dynamics. Altogether, the current study reveals that the application of both MB and PTH biostimulants promoted rapeseed plant phenology and productivity associated with an improvement in the photosynthetic machinery while boosting other physiological and molecular mechanisms.

Identification and Functional Exploration of BraGASA Genes Reveal Their Potential Roles in Drought Stress Tolerance and Sexual Reproduction in Brassica rapa L. ssp. pekinensis.

Gibberellic acid-stimulated Arabidopsis sequences (GASAs) are a subset of the gibberellin (GA)-regulated gene family and play crucial roles in various physiological processes. However, the GASA genes in Brassica rapa have not yet been documented. In this study, we identified and characterized 16 GASA genes in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Analysis of the conserved motifs revealed significant conservation within the activation segment of BraGASA genes. This gene family contains numerous promoter elements associated with abiotic stress tolerance, including those for abscisic acid (ABA) and methyl jasmonate (MeJA). Expression profiling revealed the presence of these genes in various tissues, including roots, stems, leaves, flowers, siliques, and callus tissues. When plants were exposed to drought stress, the expression of BraGASA3 decreased notably in drought-sensitive genotypes compared to their wild-type counterparts, highlighting the potentially crucial role of BraGASA3 in drought stress. Additionally, BraGASAs exhibited various functions in sexual reproduction dynamics. The findings contribute to the understanding of the function of BraGASAs and provide valuable insights for further exploration of the GASA gene function of the BraGASA gene in Chinese cabbage.

Integrated miRNA and mRNA Transcriptome Analysis Reveals Regulatory Mechanisms in the Response of Winter Brassica rapa to Drought Stress.

Drought is a major abiotic stress factor that reduces agricultural productivity. Understanding the molecular regulatory network of drought response in winter rape is of great significance for molecular Brassica rapa. In order to comprehensively analyze the network expression of DEGs and DEMIs in winter rape under drought stress, in this study we used Longyou 7 as the experimental material to identify DEGs and DEMIs related to drought stress by transcriptome and miRNA sequencing. A total of 14-15 key differential mRNA genes related to drought stress and biological stress were screened out under different treatments in the three groups. and 32 differential miRNAs were identified through targeted regulatory relationships, and the mRNA expression of 20 target genes was negatively regulated by the targeting regulatory relationship. It is mainly enriched in starch and sucrose metabolism, carbon metabolism and other pathways. Among them, gra-MIR8731-p3_2ss13GA18GA regulated the expression of multiple mRNAs in the three treatments. miRNA is mainly involved in the drought resistance of Chinese cabbage winter rape by regulating the expression of target genes, such as starch and sucrose metabolism, amino acid biosynthesis, and carbon metabolism. These miRNAs and their target genes play an indispensable role in winter rapeseed drought stress tolerance regulation.

A premature stop codon in BrFLC2 transcript results in early flowering in oilseed-type Brassica rapa plants.

KEY MESSAGE: Nonsense-mediated mRNA decay (NMD)-mediated degradation of BrFLC2 transcripts is the main cause of rapid flowering of oilseed-type B. rapa 'LP08' plants. Many Brassica species require vernalization (long-term winter-like cooling) for transition to the reproductive stage. In the past several decades, scientific efforts have been made to discern the molecular mechanisms underlying vernalization in many species. Thus, to identify the key regulators required for vernalization in Brassica rapa L., we constructed a linkage map composed of 7833 single nucleotide polymorphism markers using the late-flowering Chinese cabbage (B. rapa L. ssp. pekinensis) inbred line 'Chiifu' and the early-flowering yellow sarson (B. rapa L. ssp. trilocularis) line 'LP08' and identified a single major QTL on the upper-arm of the chromosome A02. In addition, we compared the transcriptomes of the lines 'Chiifu' and 'LP08' at five vernalization time points, including both non-vernalized and post-vernalization conditions. We observed that BrFLC2 was significantly downregulated in the early flowering 'LP08' and had two deletion sites (one at 4th exon and the other at 3' downstream region) around the BrFLC2 genomic region compared with the BrFLC2 genomic region in 'Chiifu'. Large deletion at 3' downstream region did not significantly affect transcription of both sense BrFLC2 transcript and antisense transcript, BrFLC2as along vernalization time course. However, the other deletion at 4th exon of BrFLC2 resulted in the generation of premature stop codon in BrFLC2 transcript in LP08 line. Cycloheximide treatment of LP08 line showed the de-repressed level of BrFLC2 in LP08, suggesting that low transcript level of BrFLC2 in LP08 might be caused by nonsense-mediated mRNA decay removing the nonsense transcript of BrFLC2. Collectively, this study provides a better understanding of the molecular mechanisms underlying floral transition in B. rapa.

Evidence that miR168a contributes to salinity tolerance of Brassica rapa L. via mediating melatonin biosynthesis.

Melatonin is a master regulator of diverse biological processes, including plant's abiotic stress responses and tolerance. Despite the extensive information on the role of melatonin in response to abiotic stress, how plants regulate endogenous melatonin content under stressful conditions remains largely unknown. In this study, we computationally mined Expressed Sequence Tag (EST) libraries of salinity-exposed Chinese cabbage (Brassica rapa) to identify the most reliable differentially expressed miRNA and its target gene(s). In light of these analyses, we found that miR168a potentially targets a key melatonin biosynthesis gene, namely O-METHYLTRANSFERASE 1 (OMT1). Accordingly, molecular and physiochemical evaluations were performed in a separate salinity experiment using contrasting B. rapa genotypes. Then, the association between B. rapa salinity tolerance and changes in measured molecular and physiochemical characteristics was determined. Results indicated that the expression profiles of miR168a and OMT1 significantly differed between B. rapa genotypes. Moreover, the expression profiles of miR168a and OMT1 significantly correlated with more melatonin content, robust antioxidant activities, and better ion homeostasis during salinity stress. Our results suggest that miR168a plausibly mediates melatonin biosynthesis, mainly through the OMT1 gene, under salinity conditions and thereby contributes to the salinity tolerance of B. rapa. To our knowledge, this is the first report on the role of miR168a and OMT1 in B. rapa salinity response.

Direct evidence for a new mode of plant defense against insects via a novel polygalacturonase-inhibiting protein expression strategy.

Plant cell wall-associated polygalacturonase-inhibiting proteins (PGIPs) are widely distributed in the plant kingdom. They play a crucial role in plant defense against phytopathogens by inhibiting microbial polygalacturonases (PGs). PGs hydrolyze the cell wall polysaccharide pectin and are among the first enzymes to be secreted during plant infection. Recent studies demonstrated that herbivorous insects express their own PG multi-gene families, raising the question whether PGIPs also inhibit insect PGs and protect plants from herbivores. Preliminary evidence suggested that PGIPs may negatively influence larval growth of the leaf beetle Phaedon cochleariae (Coleoptera: Chrysomelidae) and identified BrPGIP3 from Chinese cabbage (Brassica rapa ssp. pekinensis) as a candidate. PGIPs are predominantly studied in planta because their heterologous expression in microbial systems is problematic and instability and aggregation of recombinant PGIPs has complicated in vitro inhibition assays. To minimize aggregate formation, we heterologously expressed BrPGIP3 fused to a glycosylphosphatidylinositol (GPI) membrane anchor, immobilizing it on the extracellular surface of insect cells. We demonstrated that BrPGIP3_GPI inhibited several P. cochleariae PGs in vitro, providing the first direct evidence of an interaction between a plant PGIP and an animal PG. Thus, plant PGIPs not only confer resistance against phytopathogens, but may also aid in defense against herbivorous beetles.

Cold acclimation can specifically inhibit chlorophyll biosynthesis in young leaves of Pakchoi.

BACKGROUND: Leaf color is an important trait in breeding of leafy vegetables. Y-05, a pakchoi (Brassica rapa ssp. chinensis) cultivar, displays yellow inner (YIN) and green outer leaves (GOU) after cold acclimation. However, the mechanism of this special phenotype remains elusive. RESULTS: We assumed that the yellow leaf phenotype of Y-05 maybe caused by low chlorophyll content. Pigments measurements and transmission electron microscopy (TEM) analysis showed that the yellow phenotype is closely related with decreased chlorophyll content and undeveloped thylakoids in chloroplast. Transcriptomes and metabolomes sequencing were next performed on YIN and GOU. The transcriptomes data showed that 4887 differentially expressed genes (DEGs) between the YIN and GOU leaves were mostly enriched in the chloroplast- and chlorophyll-related categories, indicating that the chlorophyll biosynthesis is mainly affected during cold acclimation. Together with metabolomes data, the inhibition of chlorophyll biosynthesis is contributed by blocked 5-aminolevulinic acid (ALA) synthesis in yellow inner leaves, which is further verified by complementary and inhibitory experiments of ALA. Furthermore, we found that the blocked ALA is closely associated with increased BrFLU expression, which is indirectly altered by cold acclimation. In BrFLU-silenced pakchoi Y-05, cold-acclimated leaves still showed green phenotype and higher chlorophyll content compared with control, meaning silencing of BrFLU can rescue the leaf yellowing induced by cold acclimation. CONCLUSIONS: Our findings suggested that cold acclimation can indirectly promote the expression of BrFLU in inner leaves of Y-05 to block ALA synthesis, resulting in decreased chlorophyll content and leaf yellowing. This study revealed the underlying mechanisms of leaves color change in cold-acclimated Y-05.

Brassica rapa L. activates macrophages via Toll-like receptors.

Macrophages can initiate innate immune responses against microbes and cancer. The aim of this study was to elucidate the effects of Brassica rapa L. on macrophages. The production of interleukin (IL)-6, tumor necrosis factor (TNF)-α, and interferon-γ induced by the insoluble fraction of B. rapa L. was decreased in macrophage-depleted spleen cells compared with controls. The insoluble fraction of B. rapa L. induced expression of H-2Kb, I-Ab, CD40, and CD86, production of cytokines and nitric oxide, and phagocytic activity in RAW264 cells. After treatment with the insoluble fraction, IL-6 and TNF-α production was significantly decreased by anti-Toll-like receptor (TLR)2 mAb or polymyxin B compared with the control. Furthermore, insoluble fraction-mediated cytokine production was significantly lower in peritoneal macrophages from TLR2-/- and TLR4-/- mice compared with wild-type mice. These results suggest that B. rapa L. is a potentially effective immunomodulator for activating macrophages to prevent infections.

Comprehensive Analysis of Arabinogalactan Protein-Encoding Genes Reveals the Involvement of Three BrFLA Genes in Pollen Germination in Brassica rapa.

Arabinogalactan proteins (AGPs) are a superfamily of hydroxyproline-rich glycoproteins that are massively glycosylated, widely implicated in plant growth and development. No comprehensive analysis of the AGP gene family has been performed in Chinese cabbage (Brassica rapa ssp. chinensis). Here, we identified a total of 293 putative AGP-encoding genes in B. rapa, including 25 classical AGPs, three lysine-rich AGPs, 30 AG-peptides, 36 fasciclin-like AGPs (FLAs), 59 phytocyanin-like AGPs, 33 xylogen-like AGPs, 102 other chimeric AGPs, two non-classical AGPs and three AGP/extensin hybrids. Their protein structures, phylogenetic relationships, chromosomal location and gene duplication status were comprehensively analyzed. Based on RNA sequencing data, we found that 73 AGP genes were differentially expressed in the floral buds of the sterile and fertile plants at least at one developmental stage in B. rapa, suggesting a potential role of AGPs in male reproductive development. We further characterized BrFLA2, BrFLA28 and BrFLA32, three FLA members especially expressed in anthers, pollen grains and pollen tubes. BrFLA2, BrFLA28 and BrFLA32 are indispensable for the proper timing of pollen germination under high relative humidity. Our study greatly extends the repertoire of AGPs in B. rapa and reveals a role for three members of the FLA subfamily in pollen germination.

Nitrogen Signaling Genes and SOC1 Determine the Flowering Time in a Reciprocal Negative Feedback Loop in Chinese Cabbage (Brassica rapa L.) Based on CRISPR/Cas9-Mediated Mutagenesis of Multiple BrSOC1 Homologs.

Precise flowering timing is critical for the plant life cycle. Here, we examined the molecular mechanisms and regulatory network associated with flowering in Chinese cabbage (Brassica rapa L.) by comparative transcriptome profiling of two Chinese cabbage inbred lines, "4004" (early bolting) and "50" (late bolting). RNA-Seq and quantitative reverse transcription PCR (qPCR) analyses showed that two positive nitric oxide (NO) signaling regulator genes, nitrite reductase (BrNIR) and nitrate reductase (BrNIA), were up-regulated in line "50" with or without vernalization. In agreement with the transcription analysis, the shoots in line "50" had substantially higher nitrogen levels than those in "4004". Upon vernalization, the flowering repressor gene Circadian 1 (BrCIR1) was significantly up-regulated in line "50", whereas the flowering enhancer genes named SUPPRESSOR OF OVEREXPRESSION OF CONSTANCE 1 homologs (BrSOC1s) were substantially up-regulated in line "4004". CRISPR/Cas9-mediated mutagenesis in Chinese cabbage demonstrated that the BrSOC1-1/1-2/1-3 genes were involved in late flowering, and their expression was mutually exclusive with that of the nitrogen signaling genes. Thus, we identified two flowering mechanisms in Chinese cabbage: a reciprocal negative feedback loop between nitrogen signaling genes (BrNIA1 and BrNIR1) and BrSOC1s to control flowering time and positive feedback control of the expression of BrSOC1s.

BrPP5.2 Overexpression Confers Heat Shock Tolerance in Transgenic Brassica rapa through Inherent Chaperone Activity, Induced Glucosinolate Biosynthesis, and Differential Regulation of Abiotic Stress Response Genes.

Plant phosphoprotein phosphatases are ubiquitous and multifarious enzymes that respond to developmental requirements and stress signals through reversible dephosphorylation of target proteins. In this study, we investigated the hitherto unknown functions of Brassica rapa protein phosphatase 5.2 (BrPP5.2) by transgenic overexpression of B. rapa lines. The overexpression of BrPP5.2 in transgenic lines conferred heat shock tolerance in 65-89% of the young transgenic seedlings exposed to 46 °C for 25 min. The examination of purified recombinant BrPP5.2 at different molar ratios efficiently prevented the thermal aggregation of malate dehydrogenase at 42 °C, thus suggesting that BrPP5.2 has inherent chaperone activities. The transcriptomic dynamics of transgenic lines, as determined using RNA-seq, revealed that 997 and 1206 (FDR < 0.05, logFC ≥ 2) genes were up- and down-regulated, as compared to non-transgenic controls. Statistical enrichment analyses revealed abiotic stress response genes, including heat stress response (HSR), showed reduced expression in transgenic lines under optimal growth conditions. However, most of the HSR DEGs were upregulated under high temperature stress (37 °C/1 h) conditions. In addition, the glucosinolate biosynthesis gene expression and total glucosinolate content increased in the transgenic lines. These findings provide a new avenue related to BrPP5.2 downstream genes and their crucial metabolic and heat stress responses in plants.

Ethylene negatively mediates self-incompatibility response in Brassica rapa.

Self-incompatibility (SI) is a genetic mechanism most flowering plants adopted to reject self-pollen thus avoid inbreeding. In the Brassicaceae, self-pollen recognition triggers downstream signaling pathways to reject self-pollen. However, the downstream signaling pathways are not very clear. Here we show that ethylene negatively mediates self-incompatibility response of Chinese cabbage (Brassica rapa L. ssp. Pekinensis) via PCD in papilla cells. We found that ethylene signaling genes were upregulated after cross-pollination. Treating stigmas with ethylene, or suppressing the expression of a negative regulator of ethylene signaling, CONSTITUTIVE TRIPLE RESPONSE 1 (CTR1), caused PCD in papilla cells and broke down the self-incompatibility. On the other hand, treating stigmas with ethylene inhibitors, or suppressing the expression of ethylene-responsive factors (ERFs), inhibited PCD in papilla cells and the compatible pollination. Our study identified an additional signaling pathway mediating self-incompatibility responses in the Brassicaceae and also developed a new method in overcoming self-incompatibility to improve the efficiency of inbred line propagation in agriculture practice.

Ectopic expression of nucleolar DEAD-Box RNA helicase OsTOGR1 confers improved heat stress tolerance in transgenic Chinese cabbage.

KEY MESSAGE: The DEAD-Box RNA helicase OsTOGR1 positively regulates heat stress tolerance in Chinese cabbage. Non-heading Chinese cabbage (Brassica rapa L. ssp. chinensis) is primarily cultivated vegetable crop in Asian countries. Heat stress is one of the major threats for its growth and yield. Numerous regulatory genes in various crops have shown to contribute thermotolerance. Among them, Thermotolerant growth required 1 (TOGR1) is an important DEAD-box RNA helicase. To examine whether its role is conserved in other crops, we constructed pCAMBIA1300-pHSP:OsTOGR1 expression vector driven by the rice small heat shock protein promoter (pHSP17.9) and successfully produced transgenic non-heading Chinese cabbage plants expressing OsTOGR1 gene via Agrobacterium-mediated vacuum infiltration transformation. In total, we generated three independent transgenic cabbage lines expressing TOGR1 gene. Expression and integration of TOGR1 was confirmed by PCR, RT-PCR and qPCR in T1 and T2 generations. The relative leaf electrical conductivity of transgenic seedlings was reduced subjected to high temperature (38 °C) compared to heat shock treatment (46 °C). In addition, hypocotyl length of transgenic seedlings increased compared to wild-type plants under high temperature and heat shock treatment. Furthermore, the transgenic plants exhibited higher chlorophyll content than wild-type plants under high temperature and heat shock treatment. The transgenic seeds displayed better germination under heat shock treatment. Tested heat stress-responsive genes were also up-regulated in the transgenic plants subjected to high temperature or heat shock treatment. To the best of our knowledge, this is the first report on describing the role of DAED-Box RNA helicases in improving heat stress tolerance of transgenic plants.

Biosynthesis of cabbage phytoalexins from indole glucosinolate.

Brassica crop species are prolific producers of indole-sulfur phytoalexins that are thought to have an important role in plant disease resistance. These molecules are conspicuously absent in the model plant Arabidopsis thaliana, and little is known about the enzymatic steps that assemble the key precursor brassinin. Here, we report the minimum set of biosynthetic genes required to generate cruciferous phytoalexins starting from the well-studied glucosinolate pathway. In vitro biochemical characterization revealed an additional role for the previously described carbon-sulfur lyase SUR1 in processing cysteine-isothiocyanate conjugates, as well as the S-methyltransferase DTCMT that methylates the resulting dithiocarbamate, together completing a pathway to brassinin. Additionally, the ß-glucosidase BABG that is present in Brassica rapa but absent in Arabidopsis was shown to act as a myrosinase and may be a determinant of plants that synthesize phytoalexins from indole glucosinolate. Transient expression of the entire pathway in Nicotiana benthamiana yields brassinin, demonstrating that the biosynthesis of indole-sulfur phytoalexins can be engineered into noncruciferous plants. The identification of these biosynthetic enzymes and the heterologous reconstitution of the indole-sulfur phytoalexin pathway sheds light on an important pathway in an edible plant and opens the door to using metabolic engineering to systematically quantify the impact of cruciferous phytoalexins on plant disease resistance and human health.

[Genomic Markers Associated with Cold-Hardiness in Brassica rapa L.]

Brassica rapa L. is a valuable and widespread species, but its cultivation in risk farming areas requires high-quality cold-hardy varieties to be developed. Mechanisms of the cold stress response in plants involve expression of numerous genes, including ribosomal ones, and are related to plant chromosome variability. FISH- and PCR-based methods were used to study intraspecific chromosome variability in the number and localization of 45S and 5S rDNA clusters and also to examine a set of molecular markers associated with cold-hardiness in winter B. rapa cultivars from high-risk farming areas. Several SSR (Na10-CO3 and BrgMS5339-1) and SCAR (BoCCA1-F/BoCCA1-1R1 and BoCCA1-F/BoCCA1-2R1) markers were identified as suitable for diagnosing cold-resistant and cold-susceptible genotypes in B. rapa. Compared with fodder cultivars, oilseed and leaf cultivars were shown to have more molecular markers associated with cold-hardiness and a higher level of polymorphism for the chromosomal distribution of 45S and 5S rDNAs, including chromosome heteromorphism. Thus, the least cold-resistant genotypes were found to display the lowest level of chromosome variability in the distribution of the 45S and 5S rDNA clusters and vice versa. The findings could be useful for the development of new cold-tolerant B. rapa varieties.

Systematic identification of long non-coding RNAs during pollen development and fertilization in Brassica rapa.

The importance of long non-coding RNAs (lncRNAs) in plant development has been established, but a systematic analysis of lncRNAs expressed during pollen development and fertilization has been elusive. We performed a time series of RNA-seq experiments at five developmental stages during pollen development and three different time points after pollination in Brassica rapa and identified 12 051 putative lncRNAs. A comprehensive view of dynamic lncRNA expression networks underpinning pollen development and fertilization was provided. B. rapa lncRNAs share many common characteristics of lncRNAs: relatively short length, low expression but specific in narrow time windows, and low evolutionary conservation. Gene modules and key lncRNAs regulating reproductive development such as exine formation were uncovered. Forty-seven cis-acting lncRNAs and 451 trans-acting lncRNAs were revealed to be highly coexpressed with their target protein-coding genes. Of particular importance are the discoveries of 14 lncRNAs that were highly coexpressed with 10 function-known pollen-associated coding genes. Fifteen lncRNAs were predicted as endogenous target mimics for 13 miRNAs, and two lncRNAs were proved to be functional target mimics for miR160 after experimental verification and shown to function in pollen development. Our study provides the systematic identification of lncRNAs during pollen development and fertilization in B. rapa and forms the foundation for future genetic, genomic, and evolutionary studies.

Maternal components of RNA-directed DNA methylation are required for seed development in Brassica rapa.

Small RNAs trigger repressive DNA methylation at thousands of transposable elements in a process called RNA-directed DNA methylation (RdDM). The molecular mechanism of RdDM is well characterized in Arabidopsis, yet the biological function remains unclear, as loss of RdDM in Arabidopsis causes no overt defects, even after generations of inbreeding. It is known that 24 nucleotide Pol IV-dependent siRNAs, the hallmark of RdDM, are abundant in flowers and developing seeds, indicating that RdDM might be important during reproduction. Here we show that, unlike Arabidopsis, mutations in the Pol IV-dependent small RNA pathway cause severe and specific reproductive defects in Brassica rapa. High rates of abortion occur when seeds have RdDM mutant mothers, but not when they have mutant fathers. Although abortion occurs after fertilization, RdDM function is required in maternal somatic tissue, not in the female gametophyte or the developing zygote, suggesting that siRNAs from the maternal soma might function in filial tissues. We propose that recently outbreeding species such as B. rapa are key to understanding the role of RdDM during plant reproduction.