Barley MLA3 recognizes the host-specificity effector Pwl2 from Magnaporthe oryzae.

Plant nucleotide-binding leucine-rich repeat (NLRs) immune receptors directly or indirectly recognize pathogen-secreted effector molecules to initiate plant defense. Recognition of multiple pathogens by a single NLR is rare and usually occurs via monitoring for changes to host proteins; few characterized NLRs have been shown to recognize multiple effectors. The barley (Hordeum vulgare) NLR gene Mildew locus a (Mla) has undergone functional diversification, and the proteins encoded by different Mla alleles recognize host-adapted isolates of barley powdery mildew (Blumeria graminis f. sp. hordei [Bgh]). Here, we show that Mla3 also confers resistance to the rice blast fungus Magnaporthe oryzae in a dosage-dependent manner. Using a forward genetic screen, we discovered that the recognized effector from M. oryzae is Pathogenicity toward Weeping Lovegrass 2 (Pwl2), a host range determinant factor that prevents M. oryzae from infecting weeping lovegrass (Eragrostis curvula). Mla3 has therefore convergently evolved the capacity to recognize effectors from diverse pathogens.

Bioengineering a plant NLR immune receptor with a robust binding interface toward a conserved fungal pathogen effector.

Bioengineering of plant immune receptors has emerged as a key strategy for generating novel disease resistance traits to counteract the expanding threat of plant pathogens to global food security. However, current approaches are limited by rapid evolution of plant pathogens in the field and may lack durability when deployed. Here, we show that the rice nucleotide-binding, leucine-rich repeat (NLR) immune receptor Pik-1 can be engineered to respond to a conserved family of effectors from the multihost blast fungus pathogen Magnaporthe oryzae. We switched the effector binding and response profile of the Pik NLR from its cognate rice blast effector AVR-Pik to the host-determining factor pathogenicity toward weeping lovegrass 2 (Pwl2) by installing a putative host target, OsHIPP43, in place of the native integrated heavy metal-associated domain (generating Pikm-1OsHIPP43). This chimeric receptor also responded to other PWL alleles from diverse blast isolates. The crystal structure of the Pwl2/OsHIPP43 complex revealed a multifaceted, robust interface that cannot be easily disrupted by mutagenesis, and may therefore provide durable, broad resistance to blast isolates carrying PWL effectors in the field. Our findings highlight how the host targets of pathogen effectors can be used to bioengineer recognition specificities that have more robust properties compared to naturally evolved disease resistance genes.

Disentangling the complex gene interaction networks between rice and the blast fungus identifies a new pathogen effector.

Studies focused solely on single organisms can fail to identify the networks underlying host-pathogen gene-for-gene interactions. Here, we integrate genetic analyses of rice (Oryza sativa, host) and rice blast fungus (Magnaporthe oryzae, pathogen) and uncover a new pathogen recognition specificity of the rice nucleotide-binding domain and leucine-rich repeat protein (NLR) immune receptor Pik, which mediates resistance to M. oryzae expressing the avirulence effector gene AVR-Pik. Rice Piks-1, encoded by an allele of Pik-1, recognizes a previously unidentified effector encoded by the M. oryzae avirulence gene AVR-Mgk1, which is found on a mini-chromosome. AVR-Mgk1 has no sequence similarity to known AVR-Pik effectors and is prone to deletion from the mini-chromosome mediated by repeated Inago2 retrotransposon sequences. AVR-Mgk1 is detected by Piks-1 and by other Pik-1 alleles known to recognize AVR-Pik effectors; recognition is mediated by AVR-Mgk1 binding to the integrated heavy metal-associated (HMA) domain of Piks-1 and other Pik-1 alleles. Our findings highlight how complex gene-for-gene interaction networks can be disentangled by applying forward genetics approaches simultaneously to the host and pathogen. We demonstrate dynamic coevolution between an NLR integrated domain and multiple families of effector proteins.

Calcium-dependent protein kinase 16 phosphorylates and activates the aquaporin PIP2;2 to regulate reversible flower opening in Gentiana scabra.

Flower opening is important for successful pollination in many plant species, and some species repeatedly open and close their flowers. This is thought to be due to turgor pressure changes caused by water influx/efflux, which depends on osmotic oscillations in the cells. In some ornamental plants, water-transporting aquaporins, also known as plasma membrane intrinsic proteins (PIPs), may play an important role in flower opening. However, the molecular mechanism(s) involved in corolla movement are largely unknown. Gentian (Gentiana spp.) flowers undergo reversible movement in response to temperature and light stimuli; using gentian as a model, we showed that the Gentiana scabra aquaporins GsPIP2;2 and GsPIP2;7 regulate repeated flower opening. In particular, phosphorylation of a C-terminal serine residue of GsPIP2;2 is important for its transport activity and relates closely to the flower re-opening rate. Furthermore, GsPIP2;2 is phosphorylated and activated by the calcium (Ca2+)-dependent protein kinase GsCPK16, which is activated by elevated cytosolic Ca2+ levels in response to temperature and light stimuli. We propose that GsCPK16-dependent phosphorylation and activation of GsPIP2;2 regulate gentian flower re-opening, with stimulus-induced Ca2+ signals acting as triggers.

A blast fungus zinc-finger fold effector binds to a hydrophobic pocket in host Exo70 proteins to modulate immune recognition in rice.

Exocytosis plays an important role in plant-microbe interactions, in both pathogenesis and symbiosis. Exo70 proteins are integral components of the exocyst, an octameric complex that mediates tethering of vesicles to membranes in eukaryotes. Although plant Exo70s are known to be targeted by pathogen effectors, the underpinning molecular mechanisms and the impact of this interaction on infection are poorly understood. Here, we show the molecular basis of the association between the effector AVR-Pii of the blast fungus Maganaporthe oryzae and rice Exo70 alleles OsExo70F2 and OsExo70F3, which is sensed by the immune receptor pair Pii via an integrated RIN4/NOI domain. The crystal structure of AVR-Pii in complex with OsExo70F2 reveals that the effector binds to a conserved hydrophobic pocket in Exo70, defining an effector/target binding interface. Structure-guided and random mutagenesis validates the importance of AVR-Pii residues at the Exo70 binding interface to sustain protein association and disease resistance in rice when challenged with fungal strains expressing effector mutants. Furthermore, the structure of AVR-Pii defines a zinc-finger effector fold (ZiF) distinct from the MAX (Magnaporthe Avrs and ToxB-like) fold previously described for a majority of characterized M. oryzae effectors. Our data suggest that blast fungus ZiF effectors bind a conserved Exo70 interface to manipulate plant exocytosis and that these effectors are also baited by plant immune receptors, pointing to new opportunities for engineering disease resistance.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution.

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications.

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Rice apoplastic CBM1-interacting protein counters blast pathogen invasion by binding conserved carbohydrate binding module 1 motif of fungal proteins.

When infecting plants, fungal pathogens secrete cell wall-degrading enzymes (CWDEs) that break down cellulose and hemicellulose, the primary components of plant cell walls. Some fungal CWDEs contain a unique domain, named the carbohydrate binding module (CBM), that facilitates their access to polysaccharides. However, little is known about how plants counteract pathogen degradation of their cell walls. Here, we show that the rice cysteine-rich repeat secretion protein OsRMC binds to and inhibits xylanase MoCel10A of the blast fungus pathogen Magnaporthe oryzae, interfering with its access to the rice cell wall and degradation of rice xylan. We found binding of OsRMC to various CBM1-containing enzymes, suggesting that it has a general role in inhibiting the action of CBM1. OsRMC is localized to the apoplast, and its expression is strongly induced in leaves infected with M. oryzae. Remarkably, knockdown and overexpression of OsRMC reduced and enhanced rice defense against M. oryzae, respectively, demonstrating that inhibition of CBM1-containing fungal enzymes by OsRMC is crucial for rice defense. We also identified additional CBM-interacting proteins (CBMIPs) from Arabidopsis thaliana and Setaria italica, indicating that a wide range of plants counteract pathogens through this mechanism.

Loss of PWT7, Located on a Supernumerary Chromosome, Is Associated with Parasitic Specialization of Pyricularia oryzae on Wheat.

Pyricularia oryzae, a blast fungus of gramineous plants, is composed of various host genus-specific pathotypes. The avirulence of an Avena isolate on wheat is conditioned by PWT3 and PWT4. We isolated the third avirulence gene from the Avena isolate and designated it as PWT7. PWT7 was effective as an avirulence gene only at the seedling stage or on leaves. PWT7 homologs were widely distributed in a subpopulation of the Eleusine pathotype and the Lolium pathotype but completely absent in the Triticum pathotype (the wheat blast fungus). The PWT7 homolog found in the Eleusine pathotype was one of the five genes involved in its avirulence on wheat. A comparative analysis of distribution of PWT7 and the other two genes previously identified in the Eleusine pathotype suggested that, in the course of parasitic specialization toward the wheat blast fungus, a common ancestor of the Eleusine, Lolium, Avena, and Triticum pathotypes first lost PWT6, secondly PWT7, and, finally, the function of PWT3. PWT7 or its homologs were located on core chromosomes in Setaria and Eleusine isolates but on supernumerary chromosomes in Lolium and Avena isolates. This is an example of interchromosomal translocations of effector genes between core and supernumerary chromosomes. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Characterization of FLOWERING LOCUS C 5 in Brassica rapa L.

Brassica rapa L., which includes Chinese cabbage, turnip, and pak choi, has more complex flowering time regulation than does Arabidopsis thaliana due to the presence of multiple paralogous flowering time genes. FLOWERING LOCUS C (FLC) is one of the key genes regulating the flowering time, and B. rapa has four FLC paralogs. BrFLC5 on the reference genome is deemed a pseudogene because of a mutation (from G to A) in the splice site of the third intron, but there are some accessions with a G nucleotide in the splice site. In this study, we genotyped 310 B. rapa accessions and found that 19 had homozygous and 81 had heterozygous putative functional BrFLC5 alleles. Accessions of turnip showed the highest proportion with a functional BrFLC5 allele. BrFLC5 acts as a floral repressor when overexpressed in A. thaliana. The BrFLC5 expression level varied in pre-vernalized plants, and this transcriptional variation was not associated with the G/A polymorphism in the third intron. Three accessions having a higher BrFLC5 expression in pre-vernalized plants had a 584-bp insertion in the promoter region. Many regions homologous to this 584-bp sequence are present in the B. rapa genome, and this 584-bp inserted region has tandem duplications of an AT-rich sequence in its central region. The possibility that a high expression of a functional BrFLC5 could contribute to producing premature bolting-resistant lines in B. rapa vegetables is discussed. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-023-01405-0.

Genome analyses reveal the hybrid origin of the staple crop white Guinea yam (Dioscorea rotundata).

White Guinea yam (Dioscorea rotundata) is an important staple tuber crop in West Africa. However, its origin remains unclear. In this study, we resequenced 336 accessions of white Guinea yam and compared them with the sequences of wild Dioscorea species using an improved reference genome sequence of D. rotundata In contrast to a previous study suggesting that D. rotundata originated from a subgroup of Dioscorea praehensilis, our results suggest a hybrid origin of white Guinea yam from crosses between the wild rainforest species D. praehensilis and the savannah-adapted species Dioscorea abyssinica We identified a greater genomic contribution from D. abyssinica in the sex chromosome of Guinea yam and extensive introgression around the SWEETIE gene. Our findings point to a complex domestication scenario for Guinea yam and highlight the importance of wild species as gene donors for improving this crop through molecular breeding.

V-primer: software for the efficient design of genome-wide InDel and SNP markers from multi-sample variant call format (VCF) genotyping data.

DNA markers are indispensable tools in genetics and genomics research as well as in crop breeding, particularly for marker-assisted selection. Recent advances in next-generation sequencing technology have made it easier to obtain genome sequences for various crop species, enabling the large-scale identification of DNA polymorphisms among varieties, which in turn has made DNA marker design more accessible. However, existing primer design software is not suitable for designing many types of genome-wide DNA markers from next-generation sequencing data. Here, we describe the development of V-primer, high-throughput software for designing insertion/deletion, cleaved amplified polymorphic sequence, and single-nucleotide polymorphism (SNP) markers. We validated the applicability of these markers in different crops. In addition, we performed multiplex PCR targeted amplicon sequencing using SNP markers designed with V-primer. Our results demonstrate that V-primer facilitates the efficient and accurate design of primers and is thus a useful tool for genetics, genomics, and crop breeding. V-primer is freely available at https://github.com/ncod3/vprimer.

Quantitative trait locus (QTL) analysis and fine-mapping for Fusarium oxysporum disease resistance in Raphanus sativus using GRAS-Di technology.

Fusarium wilt is a significant disease in radish, but the genetic mechanisms controlling yellows resistance (YR) are not well understood. This study aimed to identify YR-QTLs and to fine-map one of them using F2:3 populations developed from resistant and susceptible radish parents. In this study, two high-density genetic maps each containing shared co-dominant markers and either female or male dominant markers that spanned 988.6 and 1127.5 cM with average marker densities of 1.40 and 1.53 cM, respectively, were generated using Genotyping by Random Amplicon Sequencing-Direct (GRAS-Di) technology. We identified two YR-QTLs on chromosome R2 and R7, and designated the latter as ForRs1 as the major QTL. Fine mapping narrowed down the ForRs1 locus to a 195 kb region. Among the 16 predicted genes in the delimited region, 4 genes including two receptor-like protein and -kinase genes (RLP/RLK) were identified as prime candidates for ForRs1 based on the nucleotide sequence comparisons between the parents and their predicted functions. This study is the first to use a GRAS-Di for genetic map construction of cruciferous crops and fine map the YR-QTL on the R7 chromosome of radish. These findings will provide groundbreaking insights into radish YR breeding and understanding the genetics of YR mechanism.

Genome Analysis Revives a Forgotten Hybrid Crop Edo-dokoro in the Genus Dioscorea.

A rhizomatous Dioscorea crop 'Edo-dokoro' was described in old records of Japan, but its botanical identity has not been characterized. We found that Edo-dokoro is still produced by four farmers in Tohoku-machi of the Aomori prefecture, Japan. The rhizomes of Edo-dokoro are a delicacy to the local people and are sold in the markets. Morphological characters of Edo-dokoro suggest its hybrid origin between the two species, Dioscorea tokoro and Dioscorea tenuipes. Genome analysis revealed that Edo-dokoro likely originated by hybridization of a male D. tokoro to a female D. tenuipes, followed by a backcross with a male plant of D. tokoro. Edo-dokoro is a typical minor crop possibly maintained for more than 300 years but now almost forgotten by the public. We hypothesize that there are many such uncharacterized genetic heritages passed over generations by small-scale farmers that await serious scientific investigation for future use and improvement by using modern genomics information.

The complete chloroplast genome of Onobrychis gaubae (Fabaceae-Papilionoideae): comparative analysis with related IR-lacking clade species.

BACKGROUND: Plastome (Plastid genome) sequences provide valuable markers for surveying evolutionary relationships and population genetics of plant species. Papilionoideae (papilionoids) has different nucleotide and structural variations in plastomes, which makes it an ideal model for genome evolution studies. Therefore, by sequencing the complete chloroplast genome of Onobrychis gaubae in this study, the characteristics and evolutionary patterns of plastome variations in IR-loss clade were compared. RESULTS: In the present study, the complete plastid genome of O. gaubae, endemic to Iran, was sequenced using Illumina paired-end sequencing and was compared with previously known genomes of the IRLC species of legumes. The O. gaubae plastid genome was 122,688 bp in length and included a large single-copy (LSC) region of 81,486 bp, a small single-copy (SSC) region of 13,805 bp and one copy of the inverted repeat (IRb) of 29,100 bp. The genome encoded 110 genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes and possessed 83 simple sequence repeats (SSRs) and 50 repeated structures with the highest proportion in the LSC. Comparative analysis of the chloroplast genomes across IRLC revealed three hotspot genes (ycf1, ycf2, clpP) which could be used as DNA barcode regions. Moreover, seven hypervariable regions [trnL(UAA)-trnT(UGU), trnT(GGU)-trnE(UUC), ycf1, ycf2, ycf4, accD and clpP] were identified within Onobrychis, which could be used to distinguish the Onobrychis species. Phylogenetic analyses revealed that O. gaubae is closely related to Hedysarum. The complete O. gaubae genome is a valuable resource for investigating evolution of Onobrychis species and can be used to identify related species. CONCLUSIONS: Our results reveal that the plastomes of the IRLC are dynamic molecules and show multiple gene losses and inversions. The identified hypervariable regions could be used as molecular markers for resolving phylogenetic relationships and species identification and also provide new insights into plastome evolution across IRLC.

The rice OsERF101 transcription factor regulates the NLR Xa1-mediated immunity induced by perception of TAL effectors.

Plant nucleotide-binding leucine-rich repeat receptors (NLRs) initiate immune responses by recognizing pathogen effectors. The rice gene Xa1 encodes an NLR with an N-terminal BED domain, and recognizes transcription activator-like (TAL) effectors of Xanthomonas oryzae pv oryzae (Xoo). Our goal here was to elucidate the molecular mechanisms controlling the induction of immunity by Xa1. We used yeast two-hybrid assays to screen for host factors that interact with Xa1 and identified the AP2/ERF-type transcription factor OsERF101/OsRAP2.6. Molecular complementation assays were used to confirm the interactions among Xa1, OsERF101 and two TAL effectors. We created OsERF101-overexpressing and knockout mutant lines in rice and identified genes differentially regulated in these lines, many of which are predicted to be involved in the regulation of response to stimulus. Xa1 interacts in the nucleus with the TAL effectors and OsERF101 via the BED domain. Unexpectedly, both the overexpression and the knockout lines of OsERF101 displayed Xa1-dependent, enhanced resistance to an incompatible Xoo strain. Different sets of genes were up- or downregulated in the overexpression and knockout lines. Our results indicate that OsERF101 regulates the recognition of TAL effectors by Xa1, and functions as a positive regulator of Xa1-mediated immunity. Furthermore, an additional Xa1-mediated immune pathway is negatively regulated by OsERF101.

Comparative transcriptome analysis during tuberous stem formation in Kohlrabi (B. oleracea var. gongylodes) at early growth periods (seedling stages).

Tuberous stem of kohlrabi is an important agronomic trait, however, the molecular basis of tuberization is poorly understood. To elucidate the tuberization mechanism, we conducted a comparative transcriptomic analysis between kohlrabi and broccoli at 10 and 20 days after germination (DAG) as tuberous stem initiated between these time points. A total of 5580 and 2866 differentially expressed transcripts (DETs) were identified between genotypes (kohlrabi vs. broccoli) and growth stages (10 DAG vs. 20 DAG), respectively, and most of the DETs were down-regulated in kohlrabi. Gene ontology (GO) and KEGG pathway enrichment analyses showed that the DETs between genotypes are involved in cell wall loosening and expansion, cell cycle and division, carbohydrate metabolism, hormone transport, hormone signal transduction and in several transcription factors. The DETs identified in those categories may directly/indirectly relate to the initiation and development of tuberous stem in kohlrabi. In addition, the expression pattern of the hormone synthesis related DETs coincided with the endogenous IAA, IAAsp, GA, ABA, and tZ profiles in kohlrabi and broccoli seedlings, that were revealed in our phytohormone analysis. This is the first report on comparative transcriptome analysis for tuberous stem formation in kohlrabi at early growth periods. The resulting data could provide significant insights into the molecular mechanism underlying tuberous stem development in kohlrabi as well as in other tuberous organ forming crops.

De Novo Transcriptome Analysis Reveals Flowering-Related Genes That Potentially Contribute to Flowering-Time Control in the Japanese Cultivated Gentian Gentiana triflora.

Japanese cultivated gentians are perennial plants that flower in early summer to late autumn in Japan, depending on the cultivar. Several flowering-related genes, including GtFT1 and GtTFL1, are known to be involved in regulating flowering time, but many such genes remain unidentified. In this study, we obtained transcriptome profiling data using the Gentiana triflora cultivar 'Maciry', which typically flowers in late July. We conducted deep RNA sequencing analysis using gentian plants grown under natural field conditions for three months before flowering. To investigate diurnal changes, the plants were sampled at 4 h intervals over 24 h. Using these transcriptome data, we determined the expression profiles of leaves based on homology searches against the Flowering-Interactive Database of Arabidopsis. In particular, we focused on transcription factor genes, belonging to the BBX and MADS-box families, and analyzed their developmental and diurnal variation. The expression levels of representative BBX genes were also analyzed under long- and short-day conditions using in-vitro-grown seedlings, and the expression patterns of some BBX genes differed. Clustering analysis revealed that the transcription factor genes were coexpressed with GtFT1. Overall, these expression profiles will facilitate further analysis of the molecular mechanisms underlying the control of flowering time in gentians.

Post-transcriptional gene silencing of the chalcone synthase gene CHS causes corolla lobe-specific whiting of Japanese gentian.

MAIN CONCLUSION: Post-transcriptional gene silencing of the chalcone synthase gene CHS specifically suppresses anthocyanin biosynthesis in corolla lobes and is responsible for the formation of a stripe type bicolor in Japanese gentian. The flower of Japanese gentian is a bell-shaped corolla composed of lobes and plicae, which is painted uniformly blue. However, the gentian cultivar 'Hakuju' shows bicolor phenotype (blue-white stripe corolla), in which anthocyanin accumulation is suppressed only in corolla lobes. Expression analysis indicated that steady-state levels of chalcone synthase (CHS) transcripts were remarkably reduced in corolla lobes compared with plicae during petal pigmentation initiation. However, no significant difference in expression levels of other flavonoid biosynthetic structural and regulatory genes was detected in its lobes and plicae. On feeding naringenin in white lobes, anthocyanin accumulation was recovered. Northern blotting probed with CHS confirmed the abundant accumulation of small RNAs in corolla lobes. Likewise, small RNA-seq analysis indicated that short reads from its lobes were predominantly mapped onto the 2nd exon region of the CHS gene, whereas those from the plicae were scarcely mapped. Subsequent infection with the gentian ovary ringspot virus (GORV), which had an RNA-silencing activity, showed the recovery of partial pigmentation in lobes. Hence, these results strongly suggested that suppressing anthocyanin accumulation in the lobes of bicolored 'Hakuju' was attributed to the specific degradation of CHS mRNA in corolla lobes, which was through post-transcriptional gene silencing (PTGS). Herein, we revealed the molecular mechanism of strip bicolor formation in Japanese gentian, and showed that PTGS of CHS was also responsible for flower color pattern in a floricultural plant other than petunia and dahlia.

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.