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.

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.

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.

Decoding early stress signaling waves in living plants using nanosensor multiplexing.

Increased exposure to environmental stresses due to climate change have adversely affected plant growth and productivity. Upon stress, plants activate a signaling cascade, involving multiple molecules like H2O2, and plant hormones such as salicylic acid (SA) leading to resistance or stress adaptation. However, the temporal ordering and composition of the resulting cascade remains largely unknown. In this study we developed a nanosensor for SA and multiplexed it with H2O2 nanosensor for simultaneous monitoring of stress-induced H2O2 and SA signals when Brassica rapa subsp. Chinensis (Pak choi) plants were subjected to distinct stress treatments, namely light, heat, pathogen stress and mechanical wounding. Nanosensors reported distinct dynamics and temporal wave characteristics of H2O2 and SA generation for each stress. Based on these temporal insights, we have formulated a biochemical kinetic model that suggests the early H2O2 waveform encodes information specific to each stress type. These results demonstrate that sensor multiplexing can reveal stress signaling mechanisms in plants, aiding in developing climate-resilient crops and pre-symptomatic stress diagnoses.

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.

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.

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.

Effect of CAX1a TILLING mutations on photosynthesis performance in salt-stressed Brassica rapa plants.

Salinity is an important environmental factor that reduces plant productivity in many world regions. It affects negatively photosynthesis causing a growth reduction. Likewise, calcium (Ca2+) is crucial in plant stress response. Therefore, the modification of Ca2+ cation exchangers (CAX) transporters could be a potential strategy to increase plant tolerance to salinity. Using Targeting Induced Local Lesions in Genomes (TILLING), researchers generated three mutants of Brassica rapa CAX1a transporter: BraA.cax1a-7, BraA.cax1a-4, and BraA.cax1a-12. The aim of this study was to test the effect of those mutations on salt tolerance focusing on the response to the photosynthesis process. Thus, the three BraA.cax1a mutants and the parental line (R-o-18) were grown under salinity conditions, and parameters related to biomass, photosynthesis performance, glucose-6-phosphate dehydrogenase (G6PDH, EC 1.1.1.49), and soluble carbohydrates were measured. BraA.cax1a-4 provided higher biomass and a better photosynthetic performance manifested by higher water use efficiency (WUE), Fv/Fm, electron fluxes, and Rubisco (EC 4.1.1.39) values. In addition, BraA.cax1a-4 presented increased osmotic protection through myo-inositol accumulation. On the other hand, BraA.cax1a-7 produced some negative effects on photosynthesis performance and lower G6PDH and Rubisco accumulations. Therefore, this study points out BraA.cax1a-4 as a useful mutation to improve photosynthetic performance in plants grown under saline conditions.

Haploid and Doubled Haploid Plant Production in Brassica rapa L. subsp. pekinensis Via Microspore Culture.

The production of haploid and doubled haploid plants is a biotechnological tool that shortens the breeding process of new cultivars in many species. Doubled haploid plants are homozygous at every locus and they can be utilized as parents to produce F1 hybrids. In this chapter, we describe a protocol for the production of doubled haploid plants in Brassica rapa L. subsp. pekinensis using androgenesis induced by isolated microspore cultures.

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.

CELL CYCLE SEITCH 52 regulates tillering by interacting with LATERAL SUPPRESSOR in non-heading Chinese cabbage.

With the discovery of essential genes regulating tillering, such as MONOCULM 1 (MOC1) in rice and LATERAL SUPPRESSOR (LAS in Arabidopsis, LS in tomato), research on tillering mechanisms has made great progress; however, the study of tillering in non-heading Chinese cabbage (NHCC) is rare. Here, we report that BcLAS, as a member of the GRAS family, plays an important role in the tillering of NHCC during its vegetative growth. BcLAS was almost not expressed in other examed parts except leaf axils throughout life. When the expression of BcLAS was silenced utilizing virus-induced gene silencing (VIGS) technology, we found that the tiller number of 'Maertou' decreased sharply. In 'Suzhouqing', overexpression of BcLAS significantly promoted tillering. BcCCS52, the orthologue to CELL CYCLE SEITCH 52 (CCS52), interacts with BcLAS. Downregulation of the expression of BcCCS52 promoted tillering of 'Suzhouqing'; therefore, we conclude that BcCCS52 plays a negative role in tillering regulation. Our findings reveal the tillering regulation mechanism of NHCCs at the vegetative stage and report an orthologue of CCS52 regulating tillering in NHCC.

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.

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.

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.

Effects of long light exposure and drought stress on plant growth and glucosinolate production in pak choi (Brassica rapa subsp. chinensis).

Glucosinolates (GLs), found in Brassicaceae family, are precursor metabolites with anti-cancer properties. Increased GLs have been studied under various environmental growth conditions. Pak choi (Brassica rapa subsp. chinensis) is a GL-rich vegetable. We hypothesize that long exposure to light and drought will increase the biomass of, and GL production in, pak choi. The experiment was conducted for 6 weeks. Long light exposure (20 h/day) increased, whilst drought exposure (12 h/week) decreased the plant growth. The plants exposed to a combination of drought and long light conditions showed similar growth pattern as control plants. GL production increased at week 6 in plants exposed to long light, while drought exposure had no impact on GL production, with the exception of glucoraphanin. Significant positive correlations were observed between plant growth and GL yield with accumulated light exposure time. Our findings suggest that long exposure to light can be used to increase both the biomass and GL production in pak choi.

Hydrogen peroxide is involved in strigolactone induced low temperature stress tolerance in rape seedlings (Brassica rapa L.).

Strigolactone (SL) is a plant hormone that can improve plant stress resistance by regulating physiological processes and gene expression. GR24 is a synthetic strigolactone, which can also be used as a plant growth regulator. In this paper, the effects of exogenous GR24 on the growth and development of rape (Brassica rapa L.) under low temperature (4 °C) were studied. The results showed that low temperature (4 °C) inhibited the growth of rape seedlings, and exogenous GR24 significantly alleviated the effect of low temperature stress on rape seedlings. Compared with 4 °C treatment, GR24 + 4 °C treatment can increase the cell viability, soluble protein and proline content, enhance antioxidant enzyme activity, inhibit the production of reactive oxygen species (ROS), improve photosynthesis, and reduce the relative conductivity of rape seedlings. Further research shows that H2O2 plays a central role in improving the cold resistance of rape seedlings by GR24. qRT-PCR results indicated that GR24 significantly increased the expression of genes. Mainly includes antioxidant enzyme genes, nicotinamide adenine dinucleotide phosphate (NADPH) oxidase genes, mitogen-activated protein kinase (MAPK) genes and cold-regulated genes. These results indicate that GR24 improves the cold tolerance of plants by regulating the expression of related genes. RNA-seq analysis revealed that there were 152 differentially expressed genes (DGEs) in T (4 °C)_vs_ST (GR24 + 4 °C), including 100 up-regulated genes and 52 down-regulated genes. These DEGs play an important role in carbon metabolism pathway, oxidative phosphorylation pathway, antioxidant activity and photosynthesis pathways. We selected 11 differentially expressed genes for qRT-PCR verification, and the verification results were consistent with RNA-seq results.

Mechanism of self/nonself-discrimination in Brassica self-incompatibility.

Self-incompatibility (SI) is a breeding system that promotes cross-fertilization. In Brassica, pollen rejection is induced by a haplotype-specific interaction between pistil determinant SRK (S receptor kinase) and pollen determinant SP11 (S-locus Protein 11, also named SCR) from the S-locus. Although the structure of the B. rapa S9-SRK ectodomain (eSRK) and S9-SP11 complex has been determined, it remains unclear how SRK discriminates self- and nonself-SP11. Here, we uncover the detailed mechanism of self/nonself-discrimination in Brassica SI by determining the S8-eSRK-S8-SP11 crystal structure and performing molecular dynamics (MD) simulations. Comprehensive binding analysis of eSRK and SP11 structures reveals that the binding free energies are most stable for cognate eSRK-SP11 combinations. Residue-based contribution analysis suggests that the modes of eSRK-SP11 interactions differ between intra- and inter-subgroup (a group of phylogenetically neighboring haplotypes) combinations. Our data establish a model of self/nonself-discrimination in Brassica SI.

Natural antisense transcripts of MIR398 genes suppress microR398 processing and attenuate plant thermotolerance.

MicroRNAs (miRNAs) and natural antisense transcripts (NATs) control many biological processes and have been broadly applied for genetic manipulation of eukaryotic gene expression. Still unclear, however, are whether and how NATs regulate miRNA production. Here, we report that the cis-NATs of MIR398 genes repress the processing of their pri-miRNAs. Through genome-wide analysis of RNA sequencing data, we identify cis-NATs of MIRNA genes in Arabidopsis and Brassica. In Arabidopsis, MIR398b and MIR398c are coexpressed in vascular tissues with their antisense genes NAT398b and NAT398c, respectively. Knock down of NAT398b and NAT398c promotes miR398 processing, resulting in stronger plant thermotolerance owing to silencing of miR398-targeted genes; in contrast, their overexpression activates NAT398b and NAT398c, causing poorer thermotolerance due to the upregulation of miR398-targeted genes. Unexpectedly, overexpression of MIR398b and MIR398c activates NAT398b and NAT398c. Taken together, these results suggest that NAT398b/c repress miR398 biogenesis and attenuate plant thermotolerance via a regulatory loop.