Chemical composition analysis and transcriptomics reveal the R2R3-MYB genes and phenol oxidases regulating the melanin formation in black radish.

Melanins are dark-brown to black-colored biomacromolecules which have been thoroughly studied in animals and microorganisms. However, the biochemical and molecular basis of plant melanins are poorly understood. We first characterized melanin from the black radish (Raphanus sativus var. niger) 'HLB' through spectroscopic techniques. p-Coumaric acid was identified as the main precursor of radish melanin. Moreover, a joint analysis of transcriptome and coexpression network was performed for the two radish accessions with black and white cortexes, 'HLB' and '55'. A set of R2R3-type RsMYBs and enzyme-coding genes exhibited a coexpression pattern, and were strongly correlated with melanin formation in radish. Transient overexpression of two phenol oxidases RsLAC7 (laccase 7) or RsPOD22-1 (peroxidase 22-1) resulted in a deeper brown color around the infiltration sites and a significant increase in the total phenol content. Furthermore, co-injection of the transcriptional activator RsMYB48/RsMYB97 with RsLAC7 and/or RsPOD22-1, markedly increased the yield of black extracts. Spectroscopic analyses revealed that these extracts are similar to the melanin found in 'HLB'. Our findings advance the understanding of structural information and the transcriptional regulatory mechanism underlying melanin formation in radish.

Identification and fine mapping of Brmmd1 gene controlling recessive genic male sterility in Brassica rapa L.

Recessive genic male sterility (RGMS) provides an effective approach for the commercial exploitation of heterosis, especially in Brassica crops. Although some artificial RGMS mutants have been reported in B. rapa, no causal genes derived from these natural mutants have been identified so far. In this study, a spontaneous RGMS mutant Bcajh97-01A derived from the 'Aijiaohuang' line traced back to the 1980 s was identified. Genetic analysis revealed that the RGMS trait was controlled by a single locus in the Bcajh97-01A/B system. Bulk segregant analysis (BSA) in combination with linkage analysis was employed to delimit the causal gene to an approximate 129 kb interval on chromosome A02. The integrated information of transcriptional levels and the predicted genes in the target region indicated that the Brmmd1 (BraA02g017420) encoding a PHD-containing nuclear protein was the most likely candidate gene. A 374 bp miniature inverted-repeat transposable element (MITE) was inserted into the first exon to prematurely stop the Brmmd1 gene translation, thus blocking the normal expression of this gene at the tetrad stage in the Bcajh97-01A. Additionally, a co-segregating structure variation (SV) marker was developed to rapidly screen the RGMS progenies from Bcajh97-01A/B system. Our findings reveal that BraA02g017420 is the causal gene responsible for the RGMS trait. This study lays a foundation for marker-assisted selection and further molecular mechanism exploration of pollen development in B. rapa.

Genome-Wide Identification of the CDPK Gene Family and Their Involvement in Taproot Cracking in Radish.

Taproot cracking, a severe and common physiological disorder, markedly reduces radish yield and commercial value. Calcium-dependent protein kinase (CDPK) plays a pivotal role in various plant developmental processes; however, its function in radish taproot cracking remains largely unknown. Here, 37 RsCDPK gene members were identified from the long-read radish genome "QZ-16". Phylogenetic analysis revealed that the CDPK members in radish, tomato, and Arabidopsis were clustered into four groups. Additionally, synteny analysis identified 13 segmental duplication events in the RsCDPK genes. Analysis of paraffin-embedded sections showed that the density and arrangement of fleshy taproot cortex cells are important factors that affect radish cracking. Transcriptome sequencing of the fleshy taproot cortex revealed 5755 differentially expressed genes (DEGs) (3252 upregulated and 2503 downregulated) between non-cracking radish "HongYun" and cracking radish "505". These DEGs were significantly enriched in plant hormone signal transduction, phenylpropanoid biosynthesis, and plant-pathogen interaction KEGG pathways. Furthermore, when comparing the 37 RsCDPK gene family members and RNA-seq DEGs, we identified six RsCDPK genes related to taproot cracking in radish. Soybean hairy root transformation experiments showed that RsCDPK21 significantly and positively regulates root length development. These findings provide valuable insights into the relationship between radish taproot cracking and RsCDPK gene function.

Dynamic profiling of intact glucosinolates in radish by combining UHPLC-HRMS/MS and UHPLC-QqQ-MS/MS.

Glucosinolates (GSLs) and their degradation products in radish confer plant defense, promote human health, and generate pungent flavor. However, the intact GSLs in radish have not been investigated comprehensively yet. Here, an accurate qualitative and quantitative analyses of 15 intact GSLs from radish, including four major GSLs of glucoraphasatin (GRH), glucoerucin (GER), glucoraphenin (GRE), and 4-methoxyglucobrassicin (4MGBS), were conducted using UHPLC-HRMS/MS in combination with UHPLC-QqQ-MS/MS. Simultaneously, three isomers of hexyl GSL, 3-methylpentyl GSL, and 4-methylpentyl GSL were identified in radish. The highest content of GSLs was up to 232.46 µmol/g DW at the 42 DAG stage in the 'SQY' taproot, with an approximately 184.49-fold increase compared to the lowest content in another sample. That the GSLs content in the taproots of two radishes fluctuated in a similar pattern throughout the five vegetative growth stages according to the metabolic profiling, whereas the GSLs content in the '55' leaf steadily decreased over the same period. Additionally, the proposed biosynthetic pathways of radish-specific GSLs were elucidated in this study. Our findings will provide an abundance of qualitative and quantitative data on intact GSLs, as well as a method for detecting GSLs, thus providing direction for the scientific progress and practical utilization of GSLs in radish.

Genetic Diversity Analysis and Core Germplasm Collection Construction of Radish Cultivars Based on Structure Variation Markers.

Radish is an economically important root vegetable worldwide. In this study, the 217 cultivated radish accessions were collected and genotyped. To detect the genotypes of these accessions, a total of 24 structure variation (SV) markers distributed on nine chromosomes were employed to analyze genetic diversity and construct a core germplasm collection of radish. The results of polymorphism information content (PIC) indicated a good polymorphism of these SV markers. Population structure analysis and principal component analysis (PCA) results showed that the 217 radish accessions fell into three main populations (P1, P2, and P3). Genetic diversity analysis showed that these populations were highly associated with geographical distribution. The values of the fixation index (FST) indicated a high genetic diversity between P2 and P3, and a moderate genetic diversity between P1 and P2, and P1 and P3. Furthermore, the 43 core germplasm were exploited for creating cytoplasmic male sterility (CMS) lines and cultivating new radish varieties. The high genetic diversity of 217 radish germplasms will not only provide valuable resources for future genetic mapping and functional genomic research, but also facilitate core germplasm utilization and the molecular breeding of radish.

Loss-of-function of SAWTOOTH 1 affects leaf dorsiventrality genes to promote leafy heads in lettuce.

The mechanisms underlying leafy heads in vegetables are poorly understood. Here, we cloned a quantitative trait locus (QTL) controlling leafy heads in lettuce (Lactuca sativa). The QTL encodes a transcription factor, SAWTOOTH 1 (LsSAW1), which has a BEL1-like homeodomain and is a homolog of Arabidopsis thaliana. A 1-bp deletion in Lssaw1 contributes to the development of leafy heads. Laser-capture microdissection and RNA-sequencing showed that LsSAW1 regulates leaf dorsiventrality and loss-of-function of Lssaw1 downregulates the expression of many adaxial genes but upregulates abaxial genes. LsSAW1 binds to the promoter region of the adaxial gene ASYMMETRIC LEAVES 1 (LsAS1) to upregulate its expression. Overexpression of LsAS1 compromised the effects of Lssaw1 on heading. LsSAW1 also binds to the promoter region of the abaxial gene YABBY 1 (LsYAB1), but downregulates its expression. Overexpression of LsYAB1 led to bending leaves in LsSAW1 genotypes. LsSAW1 directly interacts with KNOTTED 1 (LsKN1), which is necessary for leafy heads in lettuce. RNA-seq data showed that LsSAW1 and LsKN1 exert antagonistic effects on the expression of thousands of genes. LsSAW1 compromises the ability of LsKN1 to repress LsAS1. Our results suggest that downregulation or loss-of-function of adaxial genes and upregulation of abaxial genes allow for the development of leafy heads.

Epidermal Patterning Factor 2-like (McEPFL2): A Putative Candidate for the Continuous Ridge (cr) Fruit Skin Locus in Bitter Gourd (Momordica charantia L.).

Bitter gourd (Momordica charantia L.) is an economically important vegetable and medicinal crop in many Asian countries. Limited work has been conducted in understanding the genetic basis of horticulturally important traits in bitter gourd. Bitter gourd is consumed primarily for its young, immature fruit, and fruit appearance plays an important role in market acceptability. One such trait is the ridges on the fruit skin. In the present study, molecular mapping of a locus underlying fruit ridge continuity was conducted. Genetic analysis in segregating populations, derived from the crosses between two inbred lines Y1 with continuous ridges (CR) and Z-1-4 with discontinuous ridges (DCR), suggested that CR was controlled by a single recessive gene (cr). High-throughput genome sequencing of CR and DCR bulks combined with high-resolution genetic mapping in an F2 population delimited cr into a 108 kb region with 16 predicted genes. Sequence variation analysis and expression profiling supported the epidermal patterning factor 2-like (McEPFL2) gene as the best candidate of the cr locus. A 1 bp deletion in the first exon of McEPFL2 in Y1 which would result in a truncated McEPFL2 protein may be the causal polymorphism for the phenotypic difference between Y1 and Z-1-4. The association of this 1 bp deletion with CR was further supported by gDNA sequencing of McEPFL2 among 31 bitter gourd accessions. This work provides a foundation for understanding the genetic and molecular control of fruit epidermal pattering and development, which also facilitates marker-assisted selection in bitter melon breeding.

Identification of Novel Locus RsCr6 Related to Clubroot Resistance in Radish (Raphanus sativus L.).

Clubroot is a devastating disease that causes substantial yield loss worldwide. However, the inheritance and molecular mechanisms of clubroot resistance during pathogen infection in radish remain largely unclear. In this study, we investigated the inheritance of clubroot resistance in the F2 population derived from crossing clubroot-resistant (CR) and clubroot-susceptible inbred lines "GLX" and "XNQ," respectively. Genetic analysis revealed that a single dominant gene controlled the clubroot resistance of "GLX" with a Mendelian ratio of resistance and susceptibility of nearly 3:1. Bulked segregant analysis combined with whole-genome resequencing (BSA-seq) was performed to detect the target region of RsCr6 on chromosome Rs8. Linkage analysis revealed that the RsCr6 locus was located between two markers, HB321 and HB331, with an interval of approximately 92 kb. Based on the outcomes of transcriptome analysis, in the RsCr6 locus, the R120263140 and R120263070 genes with a possible relation to clubroot resistance were considered candidate genes. In addition, three core breeding materials containing the two reported quantitative trait loci (QTLs) and our novel locus RsCr6 targeting clubroot resistance were obtained using marker-assisted selection (MAS) technology. This study reveals a novel locus responsible for clubroot resistance in radishes. Further analysis of new genes may reveal the molecular mechanisms underlying the clubroot resistance of plants and provide a theoretical basis for radish resistance breeding.

Upregulation of a KN1 homolog by transposon insertion promotes leafy head development in lettuce.

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

Fine mapping of a candidate gene for cool-temperature-induced albinism in ornamental kale.

BACKGROUND: The symptoms of cool-temperature-induced chlorosis (CTIC) are widely existed in higher plants. Although many studies have shown that the genetic mechanism of CTIC is generally controlled by recessive genes in model plants, the dominant inheritance of albinism has not been reported thus far. Here, two CTIC mutants, Red Kamome and White Kamome, were utilized to analyse the inheritance of the albino trait in ornamental kale. The objective of this investigation is to fine-map the target locus and identify the most likely candidate genes for albinism. RESULTS: Genetic analysis revealed that the albinism in the inner leaves of ornamental kale followed semi-dominant inheritance and was controlled by a single locus in two segregating populations. BSR-seq in combination with linkage analysis was employed to fine-map the causal gene, named AK (Albino Kale), to an approximate 60 kb interval on chromosome C03. Transcriptome data from two extreme pools indicated that the differentially expressed gene of Bol015404, which encodes a cytochrome P450 protein, was the candidate gene. The Bol015404 gene was demonstrated to be upregulated in the albino leaves of ornamental kale by qPCR. Additionally, the critical temperature for the albinism was determined between 10 °C and 16 °C by gradient test. CONCLUSIONS: Using two independent segregating populations, the albino mutants were shown to be controlled by one semi-dominant gene, AK, in ornamental kale. The Bol015404 gene was co-segregated with albinism phenotypes, suggesting this unknown function P450 gene as the most likely candidate gene. The albino trait appeared caused by the low temperatures rather than photoperiod. Our results lay a solid foundation on the genetic control of albinism in ornamental kale.

Independent activation of the BoMYB2 gene leading to purple traits in Brassica oleracea.

KEY MESSAGE: Transposon insertion and point mutation independently activated the BoMYB2 gene in three purple cultivars of Brassica oleracea including kale, kohlrabi, and cabbage. Several varieties of B. oleracea have both green and purple cultivars. In this study, the causal genes for the purple traits in kale, kohlrabi and cabbage were cloned using map-based cloning approach. The purple traits in all three varieties were mapped to the same locus as the BoMYB2 gene in cauliflower. Surprisingly, the insertion of Harbinger transposon of BoMYB2 in cauliflower was not found in purple kale, kohlrabi and cabbage. Sequencing of the BoMYB2 gene in purple kale and purple kohlrabi discovered a 7606 bp CACTA-like transposon in its promoter region. Transient assay and promoter activity study showed that the insertion upregulated the expression of the BoMYB2 gene. On the other hand, the activation of the BoMYB2 gene in purple cabbage was caused by point mutation and/or 1-bp insertion in its promoter region. Sequence analysis of the BoMYB2 gene in different varieties suggested that the activating events most likely occurred independently after the divergence of cabbage, cauliflower, and kale/kohlrabi. Our results not only contribute to a better understanding of anthocyanin inheritance in B. oleracea, but also provide useful information for future hybrid breeding of purple cultivars through combination of different functional alleles of the BoMYB2 gene.

RNA sequencing provides insights into the evolution of lettuce and the regulation of flavonoid biosynthesis.

Different horticultural types of lettuce exhibit tremendous morphological variation. However, the molecular basis for domestication and divergence among the different horticultural types of lettuce remains unknown. Here, we report the RNA sequencing of 240 lettuce accessions sampled from the major horticultural types and wild relatives, generating 1.1 million single-nucleotide polymorphisms (SNPs). Demographic modeling indicates that there was a single domestication event for lettuce. We identify a list of regions as putative selective sweeps that occurred during domestication and divergence, respectively. Genome-wide association studies (GWAS) identify 5311 expression quantitative trait loci (eQTL) regulating the expression of 4105 genes, including nine eQTLs regulating genes associated with flavonoid biosynthesis. GWAS for leaf color detects six candidate loci responsible for the variation of anthocyanins in lettuce leaves. Our study provides a comprehensive understanding of the domestication and the accumulation of anthocyanins in lettuce, and will facilitate the breeding of cultivars with improved nutritional value.

GC content fluctuation around plant small RNA-generating sites.

GC content of small RNA-generating sites and their flanking sequences in Arabidopsis thaliana and rice was systematically analyzed in silico. High GC content fluctuation (GCF) is observed at the borders of sRNA sites, while the GCF within sRNA sites is low. Furthermore, the GC content along sequences of some miniature inverted-repeat transposable element (MITE) families coincides with the abundance of MITE-derived small RNAs. The GCF within tasiRNA clusters is negatively correlated with its phasing score. We conclude that high GC content and low GCF could increase the expression of small RNA. Our results provide further insights on small RNA expression, which may be applied to improve the silencing efficiency of RNAi and virus-induced gene silencing.