Differences in the transcriptional immune response to Albugo candida between white rust resistant and susceptible cultivars in Brassica rapa L.

Albugo candida causing white rust disease decreases the yield of Brassica rapa vegetables greatly. Resistant and susceptible cultivars in B. rapa vegetables have different immune responses against A. candida inoculation, however, the mechanism of how host plants respond to A. candida is still unknown. Using RNA-sequencing, we identified differentially expressed genes (DEGs) between A. candida inoculated [48 and 72 h after inoculation (HAI)] and non-inoculated samples in resistant and susceptible cultivars of komatsuna (B. rapa var. perviridis). Functional DEGs differed between the resistant and susceptible cultivars in A. candida inoculated samples. Salicylic acid (SA) responsive genes tended to be changed in their expression levels by A. candida inoculation in both resistant and susceptible cultivars, but different genes were identified in the two cultivars. SA-dependent systemic acquired resistance (SAR) involving genes were upregulated following A. candida inoculation in the resistant cultivar. Particular genes categorized as SAR that changed expression levels overlapped between A. candida and Fusarium oxysporum f. sp. conglutinans inoculated samples in resistant cultivar, suggesting a role for SAR in defense response to both pathogens particularly in the effector-triggered immunity downstream pathway. These findings will be useful for understanding white rust resistance mechanisms in B. rapa.

Whole-Genome DNA Methylation Analysis in Brassica rapa subsp. perviridis in Response to Albugo candida Infection.

DNA methylation is an epigenetic mark associated with several mechanisms in plants including immunity mechanisms. However, little is known about the regulatory role of DNA methylation in the resistance response of Brassica species against fungal diseases. White rust, caused by the fungus Albugo candida, is one of the most widespread and destructive diseases of all the cultivated Brassica species, particularly Brassica rapa L. and Brassica juncea (L.) Czern and Coss. Here, we investigate whole-genome DNA methylation modifications of B. rapa subsp. perviridis in response to white rust. As a result, 233 and 275 differentially methylated regions (DMRs) in the susceptible cultivar "Misugi" and the resistant cultivar "Nanane" were identified, respectively. In both cultivars, more than half of the DMRs were associated with genes (DMR-genes). Gene expression analysis showed that 13 of these genes were also differentially expressed between control and infected samples. Gene ontology enrichment analysis of DMR genes revealed their involvement in various biological processes including defense mechanisms. DMRs were unevenly distributed around genes in susceptible and resistant cultivars. In "Misugi," DMRs tended to be located within genes, while in "Nanane," DMRs tended to be located up and downstream of the genes. However, CG DMRs were predominantly located within genes in both cultivars. Transposable elements also showed association with all three sequence contexts of DMRs but predominantly with CHG and CHH DMRs in both cultivars. Our findings indicate the occurrence of DNA methylation modifications in B. rapa in response to white rust infection and suggest a potential regulatory role of DNA methylation modification in defense mechanisms which could be exploited to improve disease resistance.

Development of a New DNA Marker for Fusarium Yellows Resistance in Brassica rapa Vegetables.

In vegetables of Brassica rapa L., Fusarium oxysporum f. sp. rapae (For) or F. oxysporum f. sp. conglutinans (Foc) cause Fusarium yellows. A resistance gene against Foc (FocBr1) has been identified, and deletion of this gene results in susceptibility (focbr1-1). In contrast, a resistance gene against For has not been identified. Inoculation tests showed that lines resistant to Foc were also resistant to For, and lines susceptible to Foc were susceptible to For. However, prediction of disease resistance by a dominant DNA marker on FocBr1 (Bra012688m) was not associated with disease resistance of For in some komatsuna lines using an inoculation test. QTL-seq using four F2 populations derived from For susceptible and resistant lines showed one causative locus on chromosome A03, which covers FocBr1. Comparison of the amino acid sequence of FocBr1 between susceptible and resistant alleles (FocBr1 and FocBo1) showed that six amino acid differences were specific to susceptible lines. The presence and absence of FocBr1 is consistent with For resistance in F2 populations. These results indicate that FocBr1 is essential for For resistance, and changed amino acid sequences result in susceptibility to For. This susceptible allele is termed focbr1-2, and a new DNA marker (focbr1-2m) for detection of the focbr1-2 allele was developed.

Characterization of Histone H3 Lysine 4 and 36 Tri-methylation in Brassica rapa L.

Covalent modifications of histone proteins act as epigenetic regulators of gene expression. We report the distribution of two active histone marks (H3K4me3 and H3K36me3) in 14-day leaves in two lines of Brassica rapa L. by chromatin immunoprecipitation sequencing. Both lines were enriched with H3K4me3 and H3K36me3 marks at the transcription start site, and the transcription level of a gene was associated with the level of H3K4me3 and H3K36me3. H3K4me3- and H3K36me3-marked genes showed low tissue-specific gene expression, and genes with both H3K4me3 and H3K36me3 had a high level of expression and were constitutively expressed. Bivalent active and repressive histone modifications such as H3K4me3 and H3K27me3 marks or antagonistic coexistence of H3K36me3 and H3K27me3 marks were observed in some genes. Expression may be susceptible to changes by abiotic and biotic stresses in genes having both H3K4me3 and H3K27me3 marks. We showed that the presence of H3K36me3 marks was associated with different gene expression levels or tissue specificity between paralogous paired genes, suggesting that H3K36me3 might be involved in subfunctionalization of the subgenomes.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa.

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

The transcriptional response to salicylic acid plays a role in Fusarium yellows resistance in Brassica rapa L.

KEY MESSAGE: Fusarium yellows resistant and susceptible lines in Brassica rapa showed different salicylic acid responses; the resistant line showed a similar response to previous reports, but the susceptible line differed. Fusarium yellows caused by Fusarium oxysporum f. sp. conglutinans (Foc) is an important disease. Previous studies showed that genes related to salicylic acid (SA) response were more highly induced following Foc infection in Brassica rapa Fusarium yellows resistant lines than susceptible lines. However, SA-induced genes have not been identified at the whole genome level and it was unclear whether they were up-regulated by Foc inoculation. Transcriptome analysis with and without SA treatment in the B. rapa Fusarium yellows susceptible line 'Misugi' and the resistant line 'Nanane' was performed to obtain insights into the relationship between SA sensitivity/response and Fusarium yellows resistance. 'Nanane's up-regulated genes were related to SA response and down-regulated genes were related to jasmonic acid (JA) or ethylene (ET) response, but differentially expressed genes in 'Misugi' were not. This result suggests that Fusarium yellows resistant and susceptible lines have a different SA response and that an antagonistic transcription between SA and JA/ET responses was found only in a Fusarium yellows resistant line. SA-responsive genes were induced by Foc inoculation in Fusarium yellows resistant (RJKB-T23) and susceptible lines (RJKB-T24). By contrast, 39 SA-induced genes specific to RJKB-T23 might function in the defense response to Foc. In this study, SA-induced genes were identified at the whole genome level, and the possibility, the defense response to Foc observed in a resistant line could be mediated by SA-induced genes, is suggested. These results will be useful for future research concerning the SA importance in Foc or other diseases resistance in B. rapa.

Genetics of Clubroot and Fusarium Wilt Disease Resistance in Brassica Vegetables: The Application of Marker Assisted Breeding for Disease Resistance.

The genus Brassica contains important vegetable crops, which serve as a source of oil seed, condiments, and forages. However, their production is hampered by various diseases such as clubroot and Fusarium wilt, especially in Brassica vegetables. Soil-borne diseases are difficult to manage by traditional methods. Host resistance is an important tool for minimizing disease and many types of resistance (R) genes have been identified. More than 20 major clubroot (CR) disease-related loci have been identified in Brassica vegetables and several CR-resistant genes have been isolated by map-based cloning. Fusarium wilt resistant genes in Brassica vegetables have also been isolated. These isolated R genes encode the toll-interleukin-1 receptor/nucleotide-binding site/leucine-rice-repeat (TIR-NBS-LRR) protein. DNA markers that are linked with disease resistance allele have been successfully applied to improve disease resistance through marker-assisted selection (MAS). In this review, we focused on the recent status of identifying clubroot and Fusarium wilt R genes and the feasibility of using MAS for developing disease resistance cultivars in Brassica vegetables.

In Arabidopsis thaliana Heterosis Level Varies among Individuals in an F1 Hybrid Population.

Heterosis or hybrid vigour is a phenomenon in which hybrid progeny exhibit superior yield and biomass to parental lines and has been used to breed F1 hybrid cultivars in many crops. A similar level of heterosis in all F1 individuals is expected as they are genetically identical. However, we found variation in rosette size in individual F1 plants from a cross between C24 and Columbia-0 accessions of Arabidopsis thaliana. Big-sized F1 plants had 26.1% larger leaf area in the first and second leaves than medium-sized F1 plants at 14 days after sowing in spite of the identical genetic background. We identified differentially expressed genes between big- and medium-sized F1 plants by microarray; genes involved in the category of stress response were overrepresented. We made transgenic plants overexpressing 21 genes, which were differentially expressed between the two size classes, and some lines had increased plant size at 14 or 21 days after sowing but not at all time points during development. Change of expression levels in stress-responsive genes among individual F1 plants could generate the variation in plant size of individual F1 plants in A. thaliana.

The histone modification H3 lysine 27 tri-methylation has conserved gene regulatory roles in the triplicated genome of Brassica rapa L.

Brassica rapa L. is an important vegetable and oilseed crop. We investigated the distribution of the histone mark tri-methylation of H3K27 (H3K27me3) in B. rapa and its role in the control of gene expression at two stages of development (2-day cotyledons and 14-day leaves) and among paralogs in the triplicated genome. H3K27me3 has a similar distribution in two inbred lines, while there was variation of H3K27me3 sites between tissues. Sites that are specific to 2-day cotyledons have increased transcriptional activity, and low levels of H3K27me3 in the gene body region. In 14-day leaves, levels of H3K27me3 were associated with decreased gene expression. In the triplicated genome, H3K27me3 is associated with paralogs that have tissue-specific expression. Even though B. rapa and Arabidopsis thaliana are not closely related within the Brassicaceae, there is conservation of H3K27me3-marked sites in the two species. Both B. rapa and A. thaliana require vernalization for floral initiation with FLC being the major controlling locus. In all four BrFLC paralogs, low-temperature treatment increases H3K27me3 at the proximal nucleation site reducing BrFLC expression. Following return to normal temperature growth conditions, H3K27me3 spreads along all four BrFLC paralogs providing stable repression of the gene.

The role of FRIGIDA and FLOWERING LOCUS C genes in flowering time of Brassica rapa leafy vegetables.

There is a wide variation of flowering time among lines of Brassica rapa L. Most B. rapa leafy (Chinese cabbage etc.) or root (turnip) vegetables require prolonged cold exposure for flowering, known as vernalization. Premature bolting caused by low temperature leads to a reduction in the yield/quality of these B. rapa vegetables. Therefore, high bolting resistance is an important breeding trait, and understanding the molecular mechanism of vernalization is necessary to achieve this goal. In this study, we demonstrated that BrFRIb functions as an activator of BrFLC in B. rapa. We showed a positive correlation between the steady state expression levels of the sum of the BrFLC paralogs and the days to flowering after four weeks of cold treatment, suggesting that this is an indicator of the vernalization requirement. We indicate that BrFLCs are repressed by the accumulation of H3K27me3 and that the spreading of H3K27me3 promotes stable FLC repression. However, there was no clear relationship between the level of H3K27me3 in the BrFLC and the vernalization requirement. We also showed that if there was a high vernalization requirement, the rate of repression of BrFLC1 expression following prolonged cold treatments was lower.

Allele specific DNA marker for fusarium resistance gene FocBo1 in Brassica oleracea.

The fusarium yellows resistance (YR) gene FocBo1 was previously identified and the DNA markers were developed to assist the breeding of YR cultivars in Brassica oleracea. However, the further analysis revealed discrepancies between the phenotypes and the genotypes predicted by those DNA markers in cabbage commercial cultivars. Since this discrepancy seemed to be due to unknown susceptible alleles of focbo1, we sequenced the gene in 19 accessions to determine the sequence variations between alleles and found that there were two resistant FocBo1 alleles and six susceptible alleles in the investigated population. The newly designed PCR markers detected three mutations in the susceptible alleles that generate premature termination codons. These were shown to accurately distinguish resistant and susceptible alleles in more than 200 accessions of B. oleracea inbred lines and cultivars. The study revealed that the locus is represented by 37.2% resistant and 62.8% susceptible alleles within seventy-eight commercial cultivars. Structural analysis of the gene revealed that a part of the allelic variation comes from intragenic recombination between alleles. Our results enable a more precise prediction of the phenotype by marker assisted selection, promoting the production of YR cultivars in B. oleracea.

Long noncoding RNAs in Brassica rapa L. following vernalization.

Brassica rapa L. is an important agricultural crop that requires a period of prolonged cold for flowering. This process is known as vernalization. Studies have shown that long noncoding RNAs (lncRNAs) play important roles in abiotic stress responses and several cold-responsive noncoding RNAs have been suggested to be involved in vernalization. We examined the transcriptome of the Chinese cabbage inbred line (B. rapa L. var. pekinensis) RJKB-T24, and identified 1,444 long intergenic noncoding RNAs (lincRNAs), 551 natural antisense transcripts (NATs), and 93 intronic noncoding RNAs (incRNAs); 549 of the 2,088 lncRNAs significantly altered their expression in response to four weeks of cold treatment. Most differentially expressed lncRNAs did not lead to a change of expression levels in mRNAs covering or near lncRNAs, suggesting that the transcriptional responses to four weeks of cold treatment in lncRNA and mRNA are independent. However, some differentially expressed mRNAs had NATs with expression altered in the same direction. These genes were categorized as having an abiotic stress response, suggesting that the paired-expression may play a role in the transcriptional response to vernalization or cold treatment. We also identified short-term cold treatment induced NATs in BrFLC and BrMAF genes, which are involved in vernalization. The lncRNAs we identified differed from those reported in Arabidopsis thaliana, suggesting the role of lncRNAs in vernalization differ between these two species.

Genome re-sequencing, SNP analysis, and genetic mapping of the parental lines of a commercial F1 hybrid cultivar of Chinese cabbage.

The genome-wide characterization of single nucleotide polymorphism (SNP) between cultivars or between inbred lines contributes to the creation of genetic markers that are important for plant breeding. Functional markers derived from polymorphisms within genes that affect phenotypic variation are especially valuable in plant breeding. Here, we report on the genome re-sequencing and analysis of the two parental inbred lines of the commercial F1 hybrid Chinese cabbage cultivar "W77". Through the genome-wide identification and classification of the SNPs and indels present in each parental line, we identified about 1,500 putative non-functional genes in each parent. We designed cleaved amplified polymorphic sequence (CAPS) markers using specific mutations found at Eco RI restriction sites in the parental lines and confirmed their Mendelian segregation by constructing a linkage map using 96 F2 plants derived from the F1 hybrid cultivar, "W77". Our results and data will be a useful genomic resource for future studies of gene function and metagenomic studies in Chinese cabbage.

Genome-wide characterization of DNA methylation, small RNA expression, and histone H3 lysine nine di-methylation in Brassica rapa L.

Epigenetic gene regulation is crucial to plant life and can involve dynamic interactions between various histone modifications, DNA methylation, and small RNAs. Detailed analysis of epigenome information is anticipated to reveal how the DNA sequence of the genome is translated into the plant's phenotype. The aim of this study was to map the DNA methylation state at the whole genome level and to clarify the relationship between DNA methylation and transcription, small RNA expression, and histone H3 lysine 9 di-methylation (H3K9me2) in Brassica rapa. We performed whole genome bisulfite sequencing, small RNA sequencing, and chromatin immunoprecipitation sequencing using H3K9me2 antibody in a Chinese cabbage inbred line, RJKB-T24, and examined the impact of epigenetic states on transcription. Cytosine methylation in DNA was analysed in different sequence contexts (CG, CHG, and CHH) (where H could be A, C, or T) and position (promoter, exon, intron, terminator, interspersed repeat regions), and the H3K9me2 and 24 nucleotide small interfering RNAs (24 nt-siRNA) were overlaid onto the B. rapa reference genome. The epigenome was compared with that of Arabidopsis thaliana and the relationship between the position of DNA methylation and gene expression, and the involvement of 24 nt siRNAs and H3K9me2 are discussed.

Comparison of transcriptome profiles by Fusarium oxysporum inoculation between Fusarium yellows resistant and susceptible lines in Brassica rapa L.

KEY MESSAGE: Resistant and susceptible lines in Brassica rapa have different immune responses against Fusarium oxysporum inoculation. Fusarium yellows caused by Fusarium oxysporum f. sp. conglutinans (Foc) is an important disease of Brassicaceae; however, the mechanism of how host plants respond to Foc is still unknown. By comparing with and without Foc inoculation in both resistant and susceptible lines of Chinese cabbage (Brassica rapa var. pekinensis), we identified differentially expressed genes (DEGs) between the bulked inoculated (6, 12, 24, and 72 h after inoculation (HAI)) and non-inoculated samples. Most of the DEGs were up-regulated by Foc inoculation. Quantitative real-time RT-PCR showed that most up-regulated genes increased their expression levels from 24 HAI. An independent transcriptome analysis at 24 and 72 HAI was performed in resistant and susceptible lines. GO analysis using up-regulated genes at 24 HAI indicated that Foc inoculation activated systemic acquired resistance (SAR) in resistant lines and tryptophan biosynthetic process and responses to chitin and ethylene in susceptible lines. By contrast, GO analysis using up-regulated genes at 72 HAI showed the overrepresentation of some categories for the defense response in susceptible lines but not in the resistant lines. We also compared DEGs between B. rapa and Arabidopsis thaliana after F. oxysporum inoculation at the same time point, and identified genes related to defense response that were up-regulated in the resistant lines of Chinese cabbage and A. thaliana. Particular genes that changed expression levels overlapped between the two species, suggesting that they are candidates for genes involved in the resistance mechanisms against F. oxysporum.

Role of DNA methylation in hybrid vigor in Arabidopsis thaliana.

Hybrid vigor or heterosis refers to the superior performance of F1 hybrid plants over their parents. Heterosis is particularly important in the production systems of major crops. Recent studies have suggested that epigenetic regulation such as DNA methylation is involved in heterosis, but the molecular mechanism of heterosis is still unclear. To address the epigenetic contribution to heterosis in Arabidopsis thaliana, we used mutant genes that have roles in DNA methylation. Hybrids between C24 and Columbia-0 (Col) without RNA polymerase IV (Pol IV) or methyltransferase I (MET1) function did not reduce the level of biomass heterosis (as evaluated by rosette diameter). Hybrids with a mutation in decrease in dna methylation 1 (ddm1) showed a decreased heterosis level. Vegetative heterosis in the ddm1 mutant hybrid was reduced but not eliminated; a complete reduction could result if there was a change in methylation at all loci critical for generating the level of heterosis, whereas if only a proportion of the loci have methylation changes there may only be a partial reduction in heterosis.