Heterosis is a complex phenomenon in which hybrids show better phenotypic characteristics than their parents do. Chinese cabbage (Brassica rapaL. spp.pekinensis) is a popular leafy crop species, hybrids of which are widely used in commercial production; however, the molecular basis of heterosis for biomass of Chinese cabbage is poorly understood. We characterized heterosis in a Chinese cabbage F₁hybrid cultivar and its parental lines from the seedling stage to the heading stage; marked heterosis of leaf weight and biomass yield were observed. Small RNA sequencing revealed 63 and 50 differentially expressed microRNAs (DEMs) at the seedling and early-heading stages, respectively. The expression levels of the majority of miRNA clusters in the F₁hybrid were lower than the mid-parent values (MPVs). Using degradome sequencing, we identified 1,819 miRNA target genes. Gene ontology (GO) analyses demonstrated that the target genes of the MPV-DEMs and low parental expression level dominance (ELD) miRNAs were significantly enriched in leaf morphogenesis, leaf development, and leaf shaping. Transcriptome analysis revealed that the expression levels of photosynthesis and chlorophyll synthesis-related MPV-DEGs (differentially expressed genes) were significantly different in the F₁hybrid compared to the parental lines, resulting in increased photosynthesis capacity and chlorophyll content in the former. Furthermore, expression of genes known to regulate leaf development was also observed at the seedling stage. Arabidopsis plants overexpressingBrGRF4.2and bra-miR396 presented increased and decreased leaf sizes, respectively. These results provide new insight into the regulation of target genes and miRNA expression patterns in leaf size and heterosis for biomass ofB. rapa.

The genes in polyphyllins pathway mixed with other steroid biosynthetic genes form an extremely complex biosynthetic network in Paris polyphylla with a giant genome. The lack of genomic data and tissue specificity causes the study of the biosynthetic pathway notably difficult. Here, we report an effective method for the prediction of key genes of polyphyllin biosynthesis. Full-length transcriptome from eight different organs via hybrid sequencing of next generation sequencingand third generation sequencing platforms annotated two 2,3-oxidosqualene cyclases (OSCs), 216 cytochrome P450s (CYPs), and 199 UDP glycosyltransferases (UGTs). Combining metabolic differences, gene-weighted co-expression network analysis, and phylogenetic trees, the candidate ranges of OSC, CYP, and UGT genes were further narrowed down to 2, 15, and 24, respectively. Beside the three previously characterized CYPs, we identified the OSC involved in the synthesis of cycloartenol and the UGT (PpUGT73CR1) at the C-3 position of diosgenin and pennogenin in P. polyphylla. This study provides an idea for the investigation of gene cluster deficiency biosynthesis pathways in medicinal plants.

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 (BrJMJ18^Par) 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 BrJMJ18^Par can modulate the expression of BrFLC3, one of the five BrFLC orthologs. Furthermore, ChIP-seq and transcriptome data reveal that BrJMJ18^Par 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.

Brassica downy mildew, a severe disease caused by Hyaloperonospora brassicae, can cause enormous economic losses in Chinese cabbage (Brassica rapa L. ssp. pekinensis) production. Although some research has been reported recently concerning the underlying resistance to this disease, no studies have identified or characterized long noncoding RNAs involved in this defense response. In this study, using high-throughput RNA sequencing, we analyzed the disease-responding mRNAs and long noncoding RNAs in two resistant lines (T12–19 and 12–85) and one susceptible line (91–112). Clustering and Gene Ontology analysis of differentially expressed genes (DEGs) showed that more DEGs were involved in the defense response in the two resistant lines than in the susceptible line. Different expression patterns and proposed functions of differentially expressed long noncoding RNAs among T12–19, 12–85, and 91–112 indicated that each has a distinct disease response mechanism. There were significantly more cis- and trans-functional long noncoding RNAs in the resistant lines than in the susceptible line, and the genes regulated by these RNAs mostly participated in the disease defense response. Furthermore, we identified a candidate resistance-related long noncoding RNA, MSTRG.19915, which is a long noncoding natural antisense transcript of a MAPK gene, BrMAPK15. Via an agroinfiltration-mediated transient overexpression system and virus-induced gene silencing technology, BrMAPK15 was indicated to have a greater ability to defend against pathogens. MSTRG.19915-silenced seedlings showed enhanced resistance to downy mildew, probably because of the upregulated expression of BrMAPK15. This research identified and characterized long noncoding RNAs involved in resistance to downy mildew, laying a foundation for future in-depth studies of disease resistance mechanisms in Chinese cabbage.

Increasing plant resistance to Verticillium wilt (VW), which causes massive losses of Brassica rapa crops, is a challenge worldwide. However, few causal genes for VW resistance have been identified by forward genetic approaches, resulting in limited application in breeding. We combine a genome-wide association study in a natural population and quantitative trait locus mapping in an F2 population and identify that the MYB transcription factor BrMYB108 regulates plant resistance to VW. A 179 bp insertion in the BrMYB108 promoter alters its expression pattern during Verticillium longisporum (VL) infection. High BrMYB108 expression leads to high VL resistance, which is confirmed by disease resistance tests using BrMYB108 overexpression and loss-of-function mutants. Furthermore, we verify that BrMYB108 confers VL resistance by regulating reactive oxygen species (ROS) generation through binding to the promoters of respiratory burst oxidase genes (Rboh). A loss-of-function mutant of AtRbohF in Arabidopsis shows significant susceptibility to VL. Thus, BrMYB108 and its target ROS genes could be used as targets for genetic engineering for VL resistance of B. rapa.