Draft sequences of the radish (Raphanus sativus L.) genome

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Glucose Oxidase of a Crucifer-Specialized Insect: A Potential Role in Suppressing Plant Defense via Modulating Antagonistic Plant Hormones

Glucose oxidase (GOX) is a representative compound found in most insect saliva that can suppress plant-defensive responses. However, little is known about the origin and role of GOX in the crucifer-specialized pest Plutella xylostella. In this study, we showed obvious regurgitation from the larval gut of P. xylostella and identified abundant peptides highly similar to known GOX. Three PxGOX genes were verified with PxGOX2 preferentially expressed in the gut. The heterologously expressed PxGOX2 confirmed its function to be a GOX, and it was detected in plant wounds together with the gut regurgitant. Further experiments revealed that PxGOX2 functioned as an effector and may suppress defensive responses in plant through the production of H2O2, which modulates levels of antagonistic salicylic acid and jasmonic acid. However, excessive H2O2 in the host plant may be neutralized by peroxidase, thus forming defensive feedback. Our findings provided new insights into understanding the GOX-mediated insect-plant interactions.

Allele mining through TILLING and EcoTILLING approaches in vegetable crops

MAIN CONCLUSION

The present review illustrates a comprehensive overview of the allele mining for genetic improvement in vegetable crops, and allele exploration methods and their utilization in various applications related to pre-breeding of economically important traits in vegetable crops. Vegetable crops have numerous wild descendants, ancestors and terrestrial races that could be exploited to develop high-yielding and climate-resilient varieties resistant/tolerant to biotic and abiotic stresses. To further boost the genetic potential of economic traits, the available genomic tools must be targeted and re-opened for exploitation of novel alleles from genetic stocks by the discovery of beneficial alleles from wild relatives and their introgression to cultivated types. This capability would be useful for giving plant breeders direct access to critical alleles that confer higher production, improve bioactive compounds, increase water and nutrient productivity as well as biotic and abiotic stress resilience. Allele mining is a new sophisticated technique for dissecting naturally occurring allelic variants in candidate genes that influence important traits which could be used for genetic improvement of vegetable crops. Target-induced local lesions in genomes (TILLINGs) is a sensitive mutation detection avenue in functional genomics, particularly wherein genome sequence information is limited or not available. Population exposure to chemical mutagens and the absence of selectivity lead to TILLING and EcoTILLING. EcoTILLING may lead to natural induction of SNPs and InDels. It is anticipated that as TILLING is used for vegetable crops improvement in the near future, indirect benefits will become apparent. Therefore, in this review we have highlighted the up-to-date information on allele mining for genetic enhancement in vegetable crops and methods of allele exploration and their use in pre-breeding for improvement of economic traits.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

S haplotype collection in Brassicaceae crops-an updated list of S haplotypes

Self-incompatibility is the system that inhibits pollen germination and pollen tube growth by self-pollen. This trait is important for the breeding of Brassica and Raphanus species. In these species, self-incompatibility is governed by the S locus, which contains three linked genes (a set called the S haplotype), i.e., S-locus receptor kinase, S-locus cysteine-rich protein/S-locus protein 11, and S-locus glycoprotein. A large number of S haplotypes have been identified in Brassica oleracea, B. rapa, and Raphanus sativus to date, and the nucleotide sequences of their many alleles have also been registered. In this state, it is important to avoid confusion between S haplotypes, i.e., an identical S haplotype with different names and a different S haplotype with an identical S haplotype number. To mitigate this issue, we herein constructed a list of S haplotypes that are easily accessible to the latest nucleotide sequences of S-haplotype genes, together with revisions to and an update of S haplotype information. Furthermore, the histories of the S-haplotype collection in the three species are reviewed, the importance of the collection of S haplotypes as a genetic resource is discussed, and the management of information on S haplotypes is proposed.

A comprehensive and conceptual overview of omics-based approaches for enhancing the resilience of vegetable crops against abiotic stresses

Abiotic stresses adversely affect the productivity and production of vegetable crops. The increasing number of crop genomes that have been sequenced or re-sequenced provides a set of computationally anticipated abiotic stress-related responsive genes on which further research may be focused. Knowledge of omics approaches and other advanced molecular tools have all been employed to understand the complex biology of these abiotic stresses. A vegetable can be defined as any component of a plant that is eaten for food. These plant parts may be celery stems, spinach leaves, radish roots, potato tubers, garlic bulbs, immature cauliflower flowers, cucumber fruits, and pea seeds. Abiotic stresses, such as deficient or excessive water, high temperature, cold, salinity, oxidative, heavy metals, and osmotic stress, are responsible for the adverse activity in plants and, ultimately major concern for decreasing yield in many vegetable crops. At the morphological level, altered leaf, shoot and root growth, altered life cycle duration and fewer or smaller organs can be observed. Likewise different physiological and biochemical/molecular processes are also affected in response to these abiotic stresses. In order to adapt and survive in a variety of stressful situations, plants have evolved physiological, biochemical, and molecular response mechanisms. A comprehensive understanding of the vegetable's response to different abiotic stresses and the identification of tolerant genotypes are essential to strengthening each vegetable's breeding program. The advances in genomics and next-generation sequencing have enabled the sequencing of many plant genomes over the last twenty years. A combination of modern genomics (MAS, GWAS, genomic selection, transgenic breeding, and gene editing), transcriptomics, and proteomics along with next-generation sequencing provides an array of new powerful approaches to the study of vegetable crops. This review examines the overall impact of major abiotic stresses on vegetables, adaptive mechanisms and functional genomic, transcriptomic, and proteomic processes used by researchers to minimize these challenges. The current status of genomics technologies for developing adaptable vegetable cultivars that will perform better in future climates is also examined.

Vegetable biology and breeding in the genomics era

Vegetable crops provide a rich source of essential nutrients for humanity and represent critical economic values to global rural societies. However, genetic studies of vegetable crops have lagged behind major food crops, such as rice, wheat and maize, thereby limiting the application of molecular breeding. In the past decades, genome sequencing technologies have been increasingly applied in genetic studies and breeding of vegetables. In this review, we recapitulate recent progress on reference genome construction, population genomics and the exploitation of multi-omics datasets in vegetable crops. These advances have enabled an in-depth understanding of their domestication and evolution, and facilitated the genetic dissection of numerous agronomic traits, which jointly expedites the exploitation of state-of-the-art biotechnologies in vegetable breeding. We further provide perspectives of further directions for vegetable genomics and indicate how the ever-increasing omics data could accelerate genetic, biological studies and breeding in vegetable crops.

Epistatic interactions among multiple copies of FLC genes with naturally occurring insertions correlate with flowering time variation in radish

Brassicaceae crops, which underwent whole-genome triplication during their evolution, have multiple copies of flowering-related genes. Interactions among multiple gene copies may be involved in flowering time regulation; however, this mechanism is poorly understood. In this study, we performed comprehensive, high-throughput RNA sequencing analysis to identify candidate genes involved in the extremely late-bolting (LB) trait in radish. Then, we examined the regulatory roles and interactions of radish FLOWERING LOCUS C (RsFLC) paralogs, the main flowering repressor candidates. Seven flowering integrator genes, five vernalization genes, nine photoperiodic/circadian clock genes and eight genes from other flowering pathways were differentially expressed in the early-bolting (EB) cultivar 'Aokubinagafuto' and LB radish cultivar 'Tokinashi' under different vernalization conditions. In the LB cultivar, RsFLC1 and RsFLC2 expression levels were maintained after 40 days of cold exposure. Bolting time was significantly correlated with the expression rates of RsFLC1 and RsFLC2. Using the EB × LB F2 population, we performed association analyses of genotypes with or without 1910- and 1627-bp insertions in the first introns of RsFLC1 and RsFLC2, respectively. The insertion alleles prevented the repression of their respective FLC genes under cold conditions. Interestingly, genotypes homozygous for RsFLC2 insertion alleles maintained high RsFLC1 and RsFLC3 expression levels under cold conditions, and two-way analysis of variance revealed that RsFLC1 and RsFLC3 expression was influenced by the RsFLC2 genotype. Our results indicate that insertions in the first introns of RsFLC1 and RsFLC2 contribute to the late-flowering trait in radish via different mechanisms. The RsFLC2 insertion allele conferred a strong delay in bolting by inhibiting the repression of all three RsFLC genes, suggesting that radish flowering time is determined by epistatic interactions among multiple FLC gene copies.

Differences in pseudogene evolution contributed to the contrasting flavors of turnip and Chiifu, two Brassica rapa subspecies

Pseudogenes are important resources for investigation of genome evolution and genomic diversity because they are nonfunctional but have regulatory effects that influence plant adaptation and diversification. However, few systematic comparative analyses of pseudogenes in closely related species have been conducted. Here, we present a turnip (Brassica rapa ssp. rapa) genome sequence and characterize pseudogenes among diploid Brassica species/subspecies. The results revealed that the number of pseudogenes was greatest in Brassica oleracea (CC genome), followed by B. rapa (AA genome) and then Brassica nigra (BB genome), implying that pseudogene differences emerged after species differentiation. In Brassica AA genomes, pseudogenes were distributed asymmetrically on chromosomes because of numerous chromosomal insertions/rearrangements, which contributed to the diversity among subspecies. Pseudogene differences among subspecies were reflected in the flavor-related glucosinolate (GSL) pathway. Specifically, turnip had the highest content of pungent substances, probably because of expansion of the methylthioalkylmalate synthase-encoding gene family in turnips; these genes were converted into pseudogenes in B. rapa ssp. pekinensis (Chiifu). RNA interference-based silencing of the gene encoding 2-oxoglutarate-dependent dioxygenase 2, which is also associated with flavor and anticancer substances in the GSL pathway, resulted in increased abundance of anticancer compounds and decreased pungency of turnip and Chiifu. These findings revealed that pseudogene differences between turnip and Chiifu influenced the evolution of flavor-associated GSL metabolism-related genes, ultimately resulting in the different flavors of turnip and Chiifu.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

CRISPR/Cas9-mediated genome editing of RsGL1a and RsGL1b in radish (Raphanus sativus L.)

The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) is a powerful tool widely used for genome editing in various organisms, including plants. It introduces and facilitates the study of rare genetic mutations in a short time and is a potent tool to assist in plant molecular breeding. Radish (Raphanus sativus L.) is an important Brassicaceae vegetable cultivated and consumed worldwide. However, the application of the CRISPR/Cas9 system is limited by the absence of an efficient transformation system in radish. This study aimed to establish a CRISPR/Cas9 system in radish employing the Agrobacterium-mediated genetic transformation system reported recently. For this purpose, we performed genome editing using the CRISPR/Cas9 system targeting the GLABRA1 (GL1) orthologs, RsGL1a and RsGL1b, that induces leaf trichome formation in radish. A Cas9/single guide RNA (sgRNA) vector with a common sgRNA corresponding to RsGL1a and RsGL1b was transferred. A total of eight T₀ plants were analyzed, of which six (editing efficiency 75%) had a mutated RsGL1a, five (62.5%) had a mutated RsGL1b, and five showed mutations in both RsGL1a and RsGL1b. Most mutations in T₀ plants were short (<3 bp) deletions or insertions, causing frameshift mutations that might produce non-functional proteins. Chimeric mutations were detected in several T₀ generation plants. In the T₁ generation, the hairless phenotype was observed only in plants with knockout mutations in both RsGL1a and RsGL1b. The majority of mutant alleles in T₀ plants, with the exception of the chimeric mutant plants detected, were stably inherited in the T₁ generation. In conclusion, we successfully knocked out RsGL1a and RsGL1b using the CRISPR/Cas9 system and demonstrated that both RsGL1a and RsGL1b independently contribute to the induction of leaf trichome formation in radish. In this study, genome-edited plants without T-DNA, which are useful as breeding material, were obtained. The findings prove the feasibility of genome editing in radish using a CRISPR/Cas9 system that could accelerate its molecular breeding to improve agronomically desirable traits.

Genome-wide characterization of homeodomain-leucine zipper genes reveals RsHDZ17 enhances the heat tolerance in radish (Raphanus sativus L.)

Homeodomain-leucine zipper (HD-Zip) transcription factors are involved in various biological processes of plant growth, development, and abiotic stress response. However, how they regulate heat stress (HS) response remains largely unclear in plants. In this study, a total of 83 RsHD-Zip genes were firstly identified from the genome of Raphanus sativus. RNA-Seq, RT-qPCR and promoter activity assays revealed that RsHDZ17 from HD-Zip Class I was highly expressed under heat, salt, and Cd stresses. RsHDZ17 is a nuclear protein with transcriptional activity at the C-terminus. Ectopic overexpression (OE) of RsHDZ17 in Arabidopsis thaliana enhanced the HS tolerance by improving the survival rate, photosynthesis capacity, and scavenging for reactive oxygen species (ROS). In addition, transient OE of RsHDZ17 in radish cotyledons impeded cell injury and augmented ROS scavenging under HS. Moreover, yeast one-hybrid, dual-luciferase assay, and electrophoretic mobility shift assay revealed that RsHDZ17 could bind to the promoter of HSFA1e. Collectively, these pieces of evidence demonstrate that RsHDZ17 could play a positive role in thermotolerance, partially through up-regulation of the expression of HSFA1e in plants. These results provide novel insights into the role of HD-Zips in radish and facilitate genetical engineering and development of heat-tolerant radish in breeding programs.

A Reference Genome Assembly of Hybrid-Derived California Wild Radish (Raphanus sativus × raphanistrum)

For agriculturally important plants, pollination and herbivory are 2 ecological factors that play into the success of crop yields. Each is also important in natural environments where invasive plants and their effect on species interactions may alter the native ecology. The California Wild Radish (Raphanus sativus × raphanistrum), a hybrid derived from an agriculturally important crop and a nonnative cultivar, is common in California. Remarkably, it has recently replaced wild populations of both progenitor species. Experiments on phenotypic variation for petal color and antiherbivore defenses suggest both pairs of polymorphisms are maintained as a result of pollinator- and herbivore-mediated natural selection. This species provides an opportunity to understand how natural selection shapes the evolution of ecologically important traits when traits are constrained by 2 opposing forces. Here we provide the genome assembly of the California Wild Radish displaying improvement to currently existing genomes for agronomically important crucifers. This genome sequence provides the tools to dissect the genomic architecture of traits related to herbivory and pollination using natural variation in the wild as well as the ability to infer demographic and selective history in the context of hybridization. Study systems like these will improve our understanding and predictions of evolutionary change for correlated traits.

An improved Raphanus sativus cv. WK10039 genome localizes centromeres, uncovers variation of DNA methylation and resolves arrangement of the ancestral Brassica genome blocks in radish chromosomes

This study presents an improved genome of Raphanus sativus cv. WK10039 uncovering centromeres and differentially methylated regions of radish chromosomes. Comprehensive genome comparison of radish and diploid Brassica species of U's triangle reveals that R. sativus arose from the Brassica B genome lineage and is a sibling species of B. nigra. Radish (Raphanus sativus L.) is a key root vegetable crop closely related to the Brassica crop species of the family Brassicaceae. We reported a draft genome of R. sativus cv. WK10039 (Rs1.0), which had 54.6 Mb gaps. To study the radish genome and explore previously unknown regions, we generated an improved genome assembly (Rs2.0) by long-read sequencing and high-resolution genome-wide mapping of chromatin interactions. Rs2.0 was 434.9 Mb in size with 0.27 Mb gaps, and the N50 scaffold length was 37.3 Mb (40-fold larger assembly compared to Rs1.0). Approximately 38% of Rs2.0 was comprised of repetitive sequences, and 52,768 protein-coding genes and 4845 non-protein-coding genes were predicted and annotated. The improved contiguity and coverage of Rs2.0, along with the detection of highly methylated regions, enabled localization of centromeres where R. sativus-specific centromere-associated repeats, full-length OTA and CRM LTR-Gypsy retrotransposons, hAT-Ac, CMC-EnSpm and Helitron DNA transposons, and sequences highly homologous to B. nigra centromere-specific CENH3-associated CL sequences were enriched. Whole-genome bisulfite sequencing combined with mRNA sequencing identified differential epigenetic marks in the radish genome related to tissue development. Synteny comparison and genomic distance analysis of radish and three diploid Brassica species of U's triangle suggested that the radish genome arose from the Brassica B genome lineage through unique rearrangement of the triplicated ancestral Brassica genome after splitting of the Brassica A/C and B genomes.

MIG-seq is an effective method for high-throughput genotyping in wheat (Triticum spp.)

MIG-seq (Multiplexed inter-simple sequence repeats genotyping by sequencing) has been developed as a low cost genotyping technology, although the number of polymorphisms obtained is assumed to be minimal, resulting in the low application of this technique to analyses of agricultural plants. We applied MIG-seq to 12 plant species that include various crops and investigated the relationship between genome size and the number of bases that can be stably sequenced. The genome size and the number of loci, which can be sequenced by MIG-seq, are positively correlated. This is due to the linkage between genome size and the number of simple sequence repeats (SSRs) through the genome. The applicability of MIG-seq to population structure analysis, linkage mapping, and quantitative trait loci (QTL) analysis in wheat, which has a relatively large genome, was further evaluated. The results of population structure analysis for tetraploid wheat showed the differences among collection sites and subspecies, which agreed with previous findings. Additionally, in wheat biparental mapping populations, over 3,000 SNPs/indels with low deficiency were detected using MIG-seq, and the QTL analysis was able to detect recognized flowering-related genes. These results revealed the effectiveness of MIG-seq for genomic analysis of agricultural plants with large genomes, including wheat.

Salinity responses and tolerance mechanisms in underground vegetable crops: an integrative review

The present review gives an insight into the salinity stress tolerance responses and mechanisms of underground vegetable crops. Phytoprotectants, agronomic practices, biofertilizers, and modern biotechnological approaches are crucial for salinity stress management. Underground vegetables are the source of healthy carbohydrates, resistant starch, antioxidants, vitamins, mineral, and nutrients which benefit human health. Soil salinity is a serious threat to agriculture that severely affects the growth, development, and productivity of underground vegetable crops. Salt stress induces several morphological, anatomical, physiological, and biochemical changes in crop plants which include reduction in plant height, leaf area, and biomass. Also, salinity stress impedes the growth of the underground organs, which ultimately reduces crop yield. Moreover, salt stress is detrimental to photosynthesis, membrane integrity, nutrient balance, and leaf water content. Salt tolerance mechanisms involve a complex interplay of several genes, transcription factors, and proteins that are involved in the salinity tolerance mechanism in underground crops. Besides, a coordinated interaction between several phytoprotectants, phytohormones, antioxidants, and microbes is needed. So far, a comprehensive review of salinity tolerance responses and mechanisms in underground vegetables is not available. This review aims to provide a comprehensive view of salt stress effects on underground vegetable crops at different levels of biological organization and discuss the underlying salt tolerance mechanisms. Also, the role of multi-omics in dissecting gene and protein regulatory networks involved in salt tolerance mechanisms is highlighted, which can potentially help in breeding salt-tolerant underground vegetable crops.

Construction of a high-density genetic map based on specific-locus amplified fragment sequencing and identification of loci controlling anthocyanin pigmentation in Yunnan red radish

Radish (Raphanus sativus L.) belongs to the family Brassicaceae. The Yunnan red radish variety contains fairly relatively large amounts of anthocyanins, making them important raw materials for producing edible red pigment. However, the genetic mechanism underlying this pigmentation has not been fully characterized. Herein, the radish inbred line YAAS-WR1 (white root-skin and white root-flesh) was crossed with the inbred line YAAS-RR1 (red root-skin and red root-flesh) to produce F1, F2, BC1P1, and BC1P2 populations. Genetic analyses revealed that the pigmented/non-pigmented (PiN) and purple/red (PR) traits were controlled by two genetic loci. The F2 population and the specific-locus amplified fragment sequencing (SLAF-seq) technique were used to construct a high-density genetic map (1230.16 cM), which contained 4032 markers distributed in nine linkage groups, with a mean distance between markers of 0.31 cM. Additionally, two QTL (QAC1 and QAC2) considerably affecting radish pigmentation were detected. A bioinformatics analysis of the QAC1 region identified 58 predicted protein-coding genes. Of these genes, RsF3'H, which is related to anthocyanin biosynthesis, was revealed as a likely candidate gene responsible for the PR trait. The results were further verified by analyzing gene structure and expression. Regarding QAC2, RsMYB1.3 was determined to be a likely candidate gene important for the PiN trait, with a 4-bp insertion in the first exon that introduced a premature termination codon in the YAAS-WR1 sequence. Assays demonstrated that RsMYB1.3 interacted with RsTT8 and activates RsTT8 and RsUFGT expression. These findings may help clarify the complex regulatory mechanism underlying radish anthocyanin synthesis. Furthermore, this study's results may be relevant for the molecular breeding of radish to improve the anthocyanin content and appearance of the taproots.

Pathogen-driven coevolution across the CBP60 plant immune regulator subfamilies confers resilience on the regulator module

Components of the plant immune signaling network need mechanisms that confer resilience against fast-evolving pathogen effectors that target them. Among eight Arabidopsis CaM-Binding Protein (CBP) 60 family members, AtCBP60g and AtSARD1 are partially functionally redundant, major positive immune regulators, and AtCBP60a is a negative immune regulator. We investigated possible resilience-conferring evolutionary mechanisms among the CBP60a, CBP60g and SARD1 immune regulatory subfamilies. Phylogenetic analysis was used to investigate the times of CBP60 subfamily neofunctionalization. Then, using the pairwise distance rank based on the newly developed analytical platform Protein Evolution Analysis in a Euclidean Space (PEAES), hypotheses of specific coevolutionary mechanisms that could confer resilience on the regulator module were tested. The immune regulator subfamilies diversified around the time of angiosperm divergence and have been evolving very quickly. We detected significant coevolutionary interactions across the immune regulator subfamilies in all of 12 diverse core eudicot species lineages tested. The coevolutionary interactions were consistent with the hypothesized coevolution mechanisms. Despite their unusually fast evolution, members across the CBP60 immune regulator subfamilies have influenced the evolution of each other long after their diversification in a way that could confer resilience on the immune regulator module against fast-evolving pathogen effectors.

Efficient Agrobacterium-mediated genetic transformation method using hypocotyl explants of radish (Raphanus sativus L.)

To investigate the gene function of radish (Raphanus sativus L.), several attempts have been made to generate genetically transformed radish. However, no efficient and relatively simple method for the genetic transformation of radish has been developed to date. In this study, we established an Agrobacterium-mediated genetic transformation method using the hypocotyl-derived explants of radish cultivar "Pirabikku". Primarily based on the Brassica transformation procedure, we optimized it for radish transformation. Using this system, the transformation efficiency of radish hypocotyl explants by Agrobacterium tumefaciens strain GV3101 harboring pIG121-Hm was 13.3%. The copy number of transfer DNA integrated into the genome was either one or two in the four independent transgenic plants. Two of the four plants exhibited male sterility and did not produce self-pollinated seeds. Examination of the expression of the β-glucuronidase (GUS) gene in T₁ plants from fertile T₀ plants showed that the GUS genes were inherited. The improvement in the genetic transformation in this study might pave the way for accelerated molecular breeding and genetic analysis of radish.

Pan-genome of Raphanus highlights genetic variation and introgression among domesticated, wild, and weedy radishes

Post-polyploid diploidization associated with descending dysploidy and interspecific introgression drives plant genome evolution by unclear mechanisms. Raphanus is an economically and ecologically important Brassiceae genus and model system for studying post-polyploidization genome evolution and introgression. Here, we report the de novo sequence assemblies for 11 genomes covering most of the typical sub-species and varieties of domesticated, wild and weedy radishes from East Asia, South Asia, Europe, and America. Divergence among the species, sub-species, and South/East Asian types coincided with Quaternary glaciations. A genus-level pan-genome was constructed with family-based, locus-based, and graph-based methods, and whole-genome comparisons revealed genetic variations ranging from single-nucleotide polymorphisms (SNPs) to inversions and translocations of whole ancestral karyotype (AK) blocks. Extensive gene flow occurred between wild, weedy, and domesticated radishes. High frequencies of genome reshuffling, biased retention, and large-fragment translocation have shaped the genomic diversity. Most variety-specific gene-rich blocks showed large structural variations. Extensive translocation and tandem duplication of dispensable genes were revealed in two large rearrangement-rich islands. Disease resistance genes mostly resided on specific and dispensable loci. Variations causing the loss of function of enzymes modulating gibberellin deactivation were identified and could play an important role in phenotype divergence and adaptive evolution. This study provides new insights into the genomic evolution underlying post-polyploid diploidization and lays the foundation for genetic improvement of radish crops, biological control of weeds, and protection of wild species' germplasms.

QTL mapping for Fusarium wilt resistance based on the whole-genome resequencing and their association with functional genes in Raphanus sativus

Two major QTL associated with resistance to Fusarium wilt (FW) were identified using whole-genome resequencing. Sequence variations and gene expression level differences suggest that TIR-NBS and LRR-RLK are candidate genes associated with FW-resistance. Fusarium wilt (FW) caused by Fusarium oxysporum f. sp. raphani is an important disease in radish, leading to severe decrease in yield and quality. YR4 as a novel genetic source to resistant to FW was confirmed through screening with five pathogen isolates. We have generated F₂ and F_2:3 populations segregated with FW resistance using YR4 and YR18 inbred lines. The disease symptom was evaluated in F_2:3 population (n = 180) in three independent studies over two years. We identified 4 QTL including the two major QTL (FoRsR7.1^59A and FoRsR9.3^59A). FoRsR7.1^59A and FoRsR9.3^59A were detected in three replicated experiments. FoRsR7.1^59A was delimited to the 2.18-Mb physical interval on chromosome R07, with a high LOD value (5.17-12.84) and explained phenotypic variation (9.34%-27.97%). The FoRsR9.3^59A represented relatively low LOD value (3.38-4.52) and explained phenotypic variation (6.24%-8.82%). On the basis of the re-sequencing data for the parental lines, we identified five putative resistance-related genes and 13 unknown genes with sequence variations at the gene and protein levels. A semi-quantitative RT-PCR analysis revealed that Rs382940 (TIR-NBS) and Rs382200 (RLK) were expressed only in 'YR4' from 0 to 6 days after the inoculation. Moreover, Rs382950 (TIR-NBS-LRR) was more highly expressed in 'YR4' from 3 to 6 days after the inoculation. These three genes might be important for FW-resistance in radish. We identified several markers based on these potential candidate genes. The marker set should be useful for breeding system to introduce the FW resistance loci from 'YR4' to improve tolerance to FW.

Genome-wide identification of AUX/IAA in radish and functional characterization of RsIAA33 gene during taproot thickening

Auxin/indole-3-acetic acid (Aux/IAA) genes encode short lived nuclear proteins that cooperated with auxin or auxin response factor (ARF), which are involved in plant growth and developmental processes. However, it's still ambiguous how the Aux/IAA genes regulate the process governing taproot thickening in radish. Herein, 65 Aux/IAA genes were identified from the radish genome. Gene duplication analysis showed that two pairs of tandem duplication and 17 (27%) segmental duplication events were identified among Aux/IAA family genes in radish. Transcriptomic analysis revealed that most of Aux/IAA genes (52/65) exhibited differential expression pattern in different root tissues, and six root-specific genes were highly expressed in root cortex, cambium, xylem, and root tip in radish. RT-qPCR analysis showed that the expression level of RsIAA33 was the highest at cortex splitting stage (CSS), and early expanding stage (ES). Furthermore, amiRNA-mediated gene silencing of RsIAA33 indicated that it could inhibit the reproductive growth, thus promoting taproot thickening and development. These results would provide valuable information for elucidating the molecular function of Aux/IAA genes involved in taproot thickening in radish.

Genome-wide transcriptome analysis of microspore abortion initiation in radish (Raphanus sativus L.)

The use of male sterile lines is one of the ideal means in hybrid seed production. Despite the widespread application of Ogura cytoplasmic male sterile (CMS) lines, the molecular mechanisms remain largely unknown. In this study, histological analyses of floral buds from a CMS line 40MA and its corresponding maintainer line 40MB were conducted, which indicate that microspore abortion was initiated shortly after the tetrad stage. RNA sequencing was performed to analyze the transcriptomes of floral buds from the tetrad stage and the early microspore stages of these two lines. More than 39 million clean reads were generated for each library, and the portions mapped to the reference genome were all above 70.60%. To further analyze the differentially expressed genes (DEGs), the samples were grouped into four pairs, of which the pair of 40MA and 40MB at the early microspore stage showed the most DEGs (5100 members). According to the abnormal appearance of the tapetum cells in 40MA, a series of tapetum development related genes were screened and analyzed. In addition, a total of 623 genes with differential expressions in the tetrad stage, but not in the early microspore stage between the two lines were filtered as the microspore abortion initiation related candidates. Twelve genes were selected to validate the sequencing result by quantitative RT-PCR. In this study, we identified a number of candidate genes involved in the initiation of microspore degeneration, which may provide a new perspective to unravel the molecular mechanism of Ogura CMS.

Molecular basis of heterosis and related breeding strategies reveal its importance in vegetable breeding

Heterosis has historically been exploited in plants; however, its underlying genetic mechanisms and molecular basis remain elusive. In recent years, due to advances in molecular biotechnology at the genome, transcriptome, proteome, and epigenome levels, the study of heterosis in vegetables has made significant progress. Here, we present an extensive literature review on the genetic and epigenetic regulation of heterosis in vegetables. We summarize six hypotheses to explain the mechanism by which genes regulate heterosis, improve upon a possible model of heterosis that is triggered by epigenetics, and analyze previous studies on quantitative trait locus effects and gene actions related to heterosis based on analyses of differential gene expression in vegetables. We also discuss the contributions of yield-related traits, including flower, fruit, and plant architecture traits, during heterosis development in vegetables (e.g., cabbage, cucumber, and tomato). More importantly, we propose a comprehensive breeding strategy based on heterosis studies in vegetables and crop plants. The description of the strategy details how to obtain F1 hybrids that exhibit heterosis based on heterosis prediction, how to obtain elite lines based on molecular biotechnology, and how to maintain heterosis by diploid seed breeding and the selection of hybrid simulation lines that are suitable for heterosis research and utilization in vegetables. Finally, we briefly provide suggestions and perspectives on the role of heterosis in the future of vegetable breeding.

Genome sequence and analysis of a Japanese radish (Raphanus sativus) cultivar named 'Sakurajima Daikon' possessing giant root

Aim

The complex genome of a Japanese radish (Raphanus sativus) cultivar named ‘Okute-Sakurajima’ with an extremely large edible round root was analysed to explore its genomic characteristics.

Methods and Results

Single-molecule real-time technology was used to obtain long sequence reads to cover 60× of the genome. De novo assembly generated 504.5 Mb contigs consisting of 1,437 sequences with the N50 value of 1.2 Mb and included 94.1% of the core eukaryotic genes. Nine pseudomolecules, comprising 69.3% of the assembled contigs, were generated along with a high-density SNP genetic map. The sequence data thus established revealed the presence of structural variations and rearrangements in the Brassicaceae genomes.

Conclusion and perspective

A total of 89,915 genes were identified in the ‘Okute-Sakurajima’ genome, 30,033 of which were newly found in this study. The genome information reported here will not only contribute to the establishment of a new resource for the radish genomics but also provide insights into the molecular mechanisms underlying formation of the giant root.

Identification of genome-wide single-nucleotide polymorphisms among geographically diverse radish accessions

Radish (Raphanus sativus L.) is cultivated around the world as a vegetable crop and exhibits diverse morphological and physiological features. DNA polymorphisms are responsible for differences in traits among cultivars. In this study, we determined genome-wide single-nucleotide polymorphisms (SNPs) among geographically diverse radish accessions using the double-digest restriction site-associated DNA sequencing (ddRAD-Seq) method. A total of 52,559 SNPs was identified in a collection of over 500 radish accessions (cultivated and wild) from East Asia, South and Southeast Asia, and the Occident and Near East. In addition, 2,624 SNP sites without missing data (referred to as common SNP sites) were identified among 510 accessions. Genetic diversity analyses, based on the common SNP sites, divided the cultivated radish accessions into four main groups, each derived from four geographical areas (Japan, East Asia, South and Southeast Asia, and the Occident and Near East). Furthermore, we discuss the origin of cultivated radish and its migration from the West to East Asia. SNP data generated in this work will facilitate further genetic studies on the radish breeding and production of DNA markers.

Genomics Armed With Diversity Leads the Way in Brassica Improvement in a Changing Global Environment

Meeting the needs of a growing world population in the face of imminent climate change is a challenge; breeding of vegetable and oilseed Brassica crops is part of the race in meeting these demands. Available genetic diversity constituting the foundation of breeding is essential in plant improvement. Elite varieties, land races, and crop wild species are important resources of useful variation and are available from existing genepools or genebanks. Conservation of diversity in genepools, genebanks, and even the wild is crucial in preventing the loss of variation for future breeding efforts. In addition, the identification of suitable parental lines and alleles is critical in ensuring the development of resilient Brassica crops. During the past two decades, an increasing number of high-quality nuclear and organellar Brassica genomes have been assembled. Whole-genome re-sequencing and the development of pan-genomes are overcoming the limitations of the single reference genome and provide the basis for further exploration. Genomic and complementary omic tools such as microarrays, transcriptomics, epigenetics, and reverse genetics facilitate the study of crop evolution, breeding histories, and the discovery of loci associated with highly sought-after agronomic traits. Furthermore, in genomic selection, predicted breeding values based on phenotype and genome-wide marker scores allow the preselection of promising genotypes, enhancing genetic gains and substantially quickening the breeding cycle. It is clear that genomics, armed with diversity, is set to lead the way in Brassica improvement; however, a multidisciplinary plant breeding approach that includes phenotype = genotype × environment × management interaction will ultimately ensure the selection of resilient Brassica varieties ready for climate change.

Resistance to Turnip Mosaic Virus in the Family Brassicaceae

Resistance to diseases caused by turnip mosaic virus (TuMV) in crop species of the family Brassicaceae has been studied extensively, especially in members of the genus Brassica. The variation in response observed on resistant and susceptible plants inoculated with different isolates of TuMV is due to a combination of the variation in the plant resistome and the variation in the virus genome. Here, we review the breadth of this variation, both at the level of variation in TuMV sequences, with one eye towards the phylogeny and evolution of the virus, and another eye towards the nature of the various responses observed in susceptible vs. different types of resistance responses. The analyses of the viral genomes allowed comparisons of pathotyped viruses on particular indicator hosts to produce clusters of host types, while the inclusion of phylogeny data and geographic location allowed the formation of the host/geographic cluster groups, the derivation of both of which are presented here. Various studies on resistance determination in particular brassica crops sometimes led to further genetic studies, in many cases to include the mapping of genes, and in some cases to the actual identification of the genes. In addition to summarizing the results from such studies done in brassica crops, as well as in radish and Arabidopsis (the latter as a potential source of candidate genes for brassica and radish), we also summarize work done using nonconventional approaches to obtaining resistance to TuMV.

Evolution of the Brassicaceae-specific MS5-Like family and neofunctionalization of the novel MALE STERILITY 5 gene essential for male fertility in Brassica napus

New genes (or lineage-specific genes) can facilitate functional innovations. MALE STERILITY 5 (MS5) in Brassica napus is a fertility-related new gene, which has two wild-type alleles (BnMS5a and BnMS5c ) and two mutant alleles (BnMS5b and BnMS5d ) that could induce male sterility. Here, we studied the history and functional evolution of MS5 homologs in plants by phylogenetic analysis and molecular genetic experiments. We identified 727 MS5 homologs and found that they define a Brassicaceae-specific gene family that has expanded partly via multiple tandem gene duplications and also probably transpositions. The MS5 in B. napus is inherited from a basic diploid ancestor of B. rapa. Molecular genetic experiments indicate that BnMS5a and BnMS5c are functionally distinct in B. napus and that BnMS5d can inhibit BnMS5a in B. napus in a dosage-dependent manner. The BnMS5a protein can move in coordination with meiotic telomeres and interact with the nuclear envelope protein SUN1, with a possible crucial role in meiotic chromosome behavior. In summary, BnMS5 belongs to a Brassicaceae-specific new gene family, and has gained a novel function that is essential for male fertility in B. napus through neofunctionalization that has likely occurred since the origin of B. rapa.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.

Genome-wide identification of the restorer-of-fertility-like (RFL) gene family in Brassica napus and expression analysis in Shaan2A cytoplasmic male sterility

Background

Cytoplasmic male sterility (CMS) is very important in hybrid breeding. The restorer-of-fertility (Rf) nuclear genes rescue the sterile phenotype. Most of the Rf genes encode pentatricopeptide repeat (PPR) proteins.

Results

We investigated the restorer-of-fertility-like (RFL) gene family in Brassica napus. A total of 53 BnRFL genes were identified. While most of the BnRFL genes were distributed on 10 of the 19 chromosomes, gene clusters were identified on chromosomes A9 and C8. The number of PPR motifs in the BnRFL proteins varied from 2 to 19, and the majority of BnRFL proteins harbored more than 10 PPR motifs. An interaction network analysis was performed to predict the interacting partners of RFL proteins. Tissue-specific expression and RNA-seq analyses between the restorer line KC01 and the sterile line Shaan2A indicated that BnRFL1, BnRFL5, BnRFL6, BnRFL8, BnRFL11, BnRFL13 and BnRFL42 located in gene clusters on chromosomes A9 and C8 were highly expressed in KC01.

Conclusions

In the present study, identification and gene expression analysis of RFL gene family in the CMS system were conducted, and seven BnRFL genes were identified as candidates for the restorer genes in Shaan2A CMS. Taken together, this method might provide new insight into the study of Rf genes in other CMS systems.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07163-z.

Resistance Gene Analogs in the Brassicaceae: Identification, Characterization, Distribution, and Evolution

The Brassicaceae consists of a wide range of species, including important Brassica crop species and the model plant Arabidopsis (Arabidopsis thaliana). Brassica spp. crop diseases impose significant yield losses annually. A major way to reduce susceptibility to disease is the selection in breeding for resistance gene analogs (RGAs). Nucleotide binding site-leucine rich repeats (NLRs), receptor-like kinases (RLKs), and receptor-like proteins (RLPs) are the main types of RGAs; they contain conserved domains and motifs and play specific roles in resistance to pathogens. Here, all classes of RGAs have been identified using annotation and assembly-based pipelines in all available genome annotations from the Brassicaceae, including multiple genome assemblies of the same species where available (total of 32 genomes). The number of RGAs, based on genome annotations, varies within and between species. In total 34,065 RGAs were identified, with the majority being RLKs (21,691), then NLRs (8,588) and RLPs (3,786). Analysis of the RGA protein sequences revealed a high level of sequence identity, whereby 99.43% of RGAs fell into several orthogroups. This study establishes a resource for the identification and characterization of RGAs in the Brassicaceae and provides a framework for further studies of RGAs for an ultimate goal of assisting breeders in improving resistance to plant disease.

Tumorous Stem Development of Brassica Juncea: A Complex Regulatory Network of Stem Formation and Identification of Key Genes in Glucosinolate Biosynthesis

Stem mustard is a stem variety of mustard, an important Brassica vegetable. The formation and development of the tumorous stem, which is the key organ for the direct yield and quality, is a complex biological process involving morphogenesis, material accumulation and gene regulation. In this study, we demonstrated through anatomical studies that stem swelling is mainly dependent on the increase in the number of cells and the volume of parenchyma cells in the cortex and pith. To further understand transcript and metabolic changes during stem swelling, we obtained 27,901 differentially expressed genes, of which 671 were specifically detected using transcriptome sequencing technology in all four stages of stem swelling. Functional annotation identified enrichment for genes involved in photosynthesis, energy metabolism, cell growth, sulfur metabolism and glucosinolate biosynthesis. Glucosinolates are a group of nitrogen- and sulfur-containing secondary metabolites, which largely exist in the Cruciferous vegetables. HPLC analysis of the contents and components of glucosinolates in four different stem development stages revealed eight glucosinolates, namely, three aliphatic glucosinolates (sinigrin, glucoalyssin and gluconapin), four indole glucosinolates (4-hydroxyglucobrassicin, glucobrassicin, 4-methoxyglucobrassicin and neoglucobrassicin) and one aromatic glucosinolate (gluconasturtiin). All these types of glucosinolates showed a significant downward trend during the stem swelling period. The content of aliphatic glucosinolates was the highest, with sinigrin being the main component. In addition, qPCR was used to validate the expression of nine genes involved in glucosinolate biosynthesis. Most of these genes were down-regulated during stem swelling in qPCR, which is consistent with transcriptome data. These data provide a basic resource for further molecular and genetic research on Brassica juncea.

Genome-wide sRNA and mRNA transcriptomic profiling insights into dynamic regulation of taproot thickening in radish (Raphanus sativus L.)

BACKGROUND

Taproot is the main edible organ and ultimately determines radish yield and quality. However, the precise molecular mechanism underlying taproot thickening awaits further investigation in radish. Here, RNA-seq was performed to identify critical genes involved in radish taproot thickening from three advanced inbred lines with different root size.

RESULTS

A total of 2606 differentially expressed genes (DEGs) were shared between 'NAU-DY' (large acicular) and 'NAU-YB' (medium obovate), which were significantly enriched in 'phenylpropanoid biosynthesis', 'glucosinolate biosynthesis', and 'starch and sucrose metabolism' pathway. Meanwhile, a total of 16 differentially expressed miRNAs (DEMs) were shared between 'NAU-DY' and 'NAU-YH' (small circular), whereas 12 miRNAs exhibited specific differential expression in 'NAU-DY'. Association analysis indicated that miR393a-bHLH77, miR167c-ARF8, and miR5658-APL might be key factors to biological phenomenon of taproot type variation, and a putative regulatory model of taproot thickening and development was proposed. Furthermore, several critical genes including SUS1, EXPB3, and CDC5 were characterized and profiled by RT-qPCR analysis.

CONCLUSION

This integrated study on the transcriptional and post-transcriptional profiles could provide new insights into comprehensive understanding of the molecular regulatory mechanism underlying taproot thickening in root vegetable crops.

Genetic characterization of a new radish introgression line carrying the restorer gene for Ogura CMS in Brassica napus

Creating a homologous restorer line for Ogura cytoplasmic male sterility (Ogu-CMS) in Brassica napus is meaningful for the wider application of Ogu-CMS system in rapeseed production. Previously, an independent development of a new Ogu-CMS restorer line (CLR650) was reported locally from crossing between Raphanobrassica (AACCRR, 2n = 56) and B. napus and a new version of Ogu CMS lines CLR6430 derived from CLR650 was characterized in this study. The results showed that the fertility restoration gene in CLR6430 presented a distorted segregation in different segregating populations. However, the majority of somatic cells from roots had a regular chromosome number (2n = 38) and no radish signal covered a whole chromosome was detected using GISH. Thirty-two specific markers derived from the introgressed radish fragments were developed based on the re-sequencing results. Unique radish insertions and differences between CLR6430 and R2000 were also identified through both radish-derived markers and PCR product sequences. Further investigations on the genetic behaviors, interactions between the fertility restoration and other traits and specific molecular markers to the introgression in CLR6430 were also conducted in this study. These results should provide the evidence of nucleotide differences between CLR6430 and R2000, and the specific markers will be helpful for breeding new Ogura restore lines in future.

Gene Regulatory Network Guided Investigations and Engineering of Storage Root Development in Root Crops

The plasticity of plant development relies on its ability to balance growth and stress resistance. To do this, plants have established highly coordinated gene regulatory networks (GRNs) of the transcription factors and signaling components involved in developmental processes and stress responses. In root crops, yields of storage roots are mainly determined by secondary growth driven by the vascular cambium. In relation to this, a dynamic yet intricate GRN should operate in the vascular cambium, in coordination with environmental changes. Despite the significance of root crops as food sources, GRNs wired to mediate secondary growth in the storage root have just begun to emerge, specifically with the study of the radish. Gene expression data available with regard to other important root crops are not detailed enough for us directly to infer underlying molecular mechanisms. Thus, in this review, we provide a general overview of the regulatory programs governing the development and functions of the vascular cambium in model systems, and the role of the vascular cambium on the growth and yield potential of the storage roots in root crops. We then undertake a reanalysis of recent gene expression data generated for major root crops and discuss common GRNs involved in the vascular cambium-driven secondary growth in storage roots using the wealth of information available in Arabidopsis. Finally, we propose future engineering schemes for improving root crop yields by modifying potential key nodes in GRNs.

Meiotic Chromosome Stability and Suppression of Crossover Between Non-homologous Chromosomes in xBrassicoraphanus, an Intergeneric Allotetraploid Derived From a Cross Between Brassica rapa and Raphanus sativus

Hybridization and polyploidization are major driving forces in plant evolution. Allopolyploids can be occasionally formed from a cross between distantly related species but often suffer from chromosome instability and infertility. xBrassicoraphanus is an intergeneric allotetraploid (AARR; 2n = 38) derived from a cross between Brassica rapa (AA; 2n = 20) and Raphanus sativus (RR; 2n = 18). xBrassicoraphanus is fertile and genetically stable, while retaining complete sets of both B. rapa and R. sativus chromosomes. Precise control of meiotic recombination is essential for the production of balanced gametes, and crossovers (COs) must occur exclusively between homologous chromosomes. Many interspecific hybrids have problems with meiotic division at early generations, in which interactions between non-homologous chromosomes often bring about aneuploidy and unbalanced gamete formation. We analyzed meiotic chromosome behaviors in pollen mother cells (PMCs) of allotetraploid and allodiploid F1 individuals of newly synthesized xBrassicoraphanus. Allotetraploid xBrassicoraphanus PMCs showed a normal diploid-like meiotic behavior. By contrast, allodiploid xBrassicoraphanus PMCs displayed abnormal segregation of chromosomes mainly due to the absence of homologous pairs. Notably, during early stages of meiosis I many of allodiploid xBrassicoraphanus chromosomes behave independently with few interactions between B. rapa and R. sativus chromosomes, forming many univalent chromosomes before segregation. Chromosomes were randomly assorted at later stages of meiosis, and tetrads with unequal numbers of chromosomes were formed at completion of meiosis. Immunolocalization of HEI10 protein mediating meiotic recombination revealed that COs were more frequent in synthetic allotetraploid xBrassicoraphanus than in allodiploid, but less than in the stabilized line. These findings suggest that structural dissimilarity between B. rapa and R. sativus chromosomes prevents non-homologous interactions between the parental chromosomes in allotetraploid xBrassicoraphanus, allowing normal diploid-like meiosis when homologous pairing partners are present. This study also suggests that CO suppression between non-homologous chromosomes is required for correct meiotic progression in newly synthesized allopolyploids, which is important for the formation of viable gametes and reproductive success in the hybrid progeny.

Identification of Conserved Gene-Regulatory Networks that Integrate Environmental Sensing and Growth in the Root Cambium

Cambium drives the lateral growth of stems and roots, contributing to diverse plant growth forms. The root crop is one of the outstanding examples of the cambium-driven growth. To understand its molecular basis, we used radish to generate a compendium of root-tissue- and stage-specific transcriptomes from two contrasting inbred lines during root growth. Expression patterns of key cambium regulators and hormone signaling components were validated. Clustering and gene ontology (GO) enrichment analyses of radish datasets followed by a comparative analysis against the newly established Arabidopsis early cambium data revealed evolutionary conserved stress-response transcription factors that may intimately control the cambium. Indeed, an in vivo network consisting of selected stress-response and cambium regulators indicated ERF-1 as a potential key checkpoint of cambial activities, explaining how cambium-driven growth is altered in response to environmental changes. The findings here provide valuable information about dynamic gene expression changes during cambium-driven root growth and have implications with regard to future engineering schemes, leading to better crop yields.

Gibberellin Promotes Bolting and Flowering via the Floral Integrators RsFT and RsSOC1-1 under Marginal Vernalization in Radish

Gibberellic acid (GA) is one of the factors that promotes flowering in radish (Raphanus Sativus L.), although the mechanism mediating GA activation of flowering has not been determined. To identify this mechanism in radish, we compared the effects of GA treatment on late-flowering (NH-JS1) and early-flowering (NH-JS2) radish lines. GA treatment promoted flowering in both lines, but not without vernalization. NH-JS2 plants displayed greater bolting and flowering pathway responses to GA treatment than NH-JS1. This variation was not due to differences in GA sensitivity in the two lines. We performed RNA-seq analysis to investigate GA-mediated changes in gene expression profiles in the two radish lines. We identified 313 upregulated, differentially expressed genes (DEGs) and 207 downregulated DEGs in NH-JS2 relative to NH-JS1 in response to GA. Of these, 21 and 8 genes were identified as flowering time and GA-responsive genes, respectively. The results of RNA-seq and quantitative PCR (qPCR) analyses indicated that RsFT and RsSOC1-1 expression levels increased after GA treatment in NH-JS2 plants but not in NH-JS1. These results identified the molecular mechanism underlying differences in the flowering-time genes of NH-JS1 and NH-JS2 after GA treatment under insufficient vernalization conditions.

Genome-based breeding approaches in major vegetable crops

Vegetable crops are major nutrient sources for humanity and have been well-cultivated since thousands of years of domestication. With the rapid development of next-generation sequencing and high-throughput genotyping technologies, the reference genome of more than 20 vegetables have been well-assembled and published. Resequencing approaches on large-scale germplasm resources have clarified the domestication and improvement of vegetable crops by human selection; its application on genetic mapping and quantitative trait locus analysis has led to the discovery of key genes and molecular markers linked to important traits in vegetables. Moreover, genome-based breeding has been utilized in many vegetable crops, including Solanaceae, Cucurbitaceae, Cruciferae, and other families, thereby promoting molecular breeding at a single-nucleotide level. Thus, genome-wide SNP markers have been widely used, and high-throughput genotyping techniques have become one of the most essential methods in vegetable breeding. With the popularization of gene editing technology research on vegetable crops, breeding efficiency can be rapidly increased, especially by combining the genomic and variomic information of vegetable crops. This review outlines the present genome-based breeding approaches used for major vegetable crops to provide insights into next-generation molecular breeding for the increasing global population.

Quantitative Trait Loci (QTLs) Associated with Microspore Culture in Raphanus sativus L. (Radish)

The radish is a highly self-incompatible plant, and consequently it is difficult to produce homozygous lines. Bud pollination in cross-fertilization plants should be done by opening immature pollen and attaching pollen to mature flowers. It accordingly takes a lot of time and effort to develop lines with fixed alleles. In the current study, a haploid breeding method has been applied to obtain homozygous plants in a short period of time by doubling chromosomes through the induction of a plant body in the haploid cells, in order to shorten the time to breed inbred lines. We constructed genetic maps with an F1 population derived by crossing parents that show a superior and inferior ability to regenerate microspores, respectively. Genetic maps were constructed from the maternal and parental maps, separately, using the two-way pseudo-testcross model. The phenotype of the regeneration rate was examined by microspore cultures and a quantitative trait loci (QTL) analysis was performed based on the regeneration rate. From the results of the culture of microspores in the F1 population, more than half of the group did not regenerate, and only a few showed a high regeneration rate. A total of five significant QTLs were detected in the F1 population, and five candidate genes were found based on the results. These candidate genes are divided into two classes, and appear to be related to either PRC2 subunits or auxin synthesis.

Characterization of the OFP Gene Family and its Putative Involvement of Tuberous Root Shape in Radish

The shape of the tuberous root, a very important quality trait, varies dramatically among radish cultivars. Ovate family proteins (OFPs) are plant-specific proteins that regulate multiple aspects of plant growth and development. To investigate the possible role of OFPs in radish tuberous root formation, 35 putative RsOFPs were identified from radish, and their expression patterns were detected during tuberous root development in six different radish cultivars. Phylogenetically, RsOFP2.3 clustered together with AtOFP1 and other members of this family that are known to regulate organ shape. Moreover, RsOFP2.3 expression was negatively correlated with tuberous root elongation after the cortex splitting stage, which made this gene the top candidate for the involvement of tuberous root shape. To further characterize the function of RsOFP2.3, it was ectopically expressed in Arabidopsis. RsOFP2.3 overexpression in Arabidopsis led to multiple phenotypical changes, especially the decreased length and increased width of the hypocotyl. Furthermore, RsOFP2.3 expression was induced by all the five classic plant hormones except ethylene, and it was most sensitive to exogenous gibberellic acid treatment. We also found that RsOFP2.3 was localized in the cytoplasm. Taken together, our results suggested the possible involvement for RsOFP2.3 in suppressing radish tuberous root elongation and that it encodes a functional protein which mainly inhibits the elongation of Arabidopsis aerial organs.

Identification and differential expression analysis of anthocyanin biosynthetic genes in root-skin color variants of radish (Raphanus sativus L.)

BACKGROUND

Taproot skin color is a major trait for assessing the commercial and nutritional quality of radish, and red-skinned radish is confirmed to improve consumer's interest and health. However, little is known about the molecular mechanisms responsible for controlling the formation of red-skinned radish.

OBJECTIVE

This study aimed to identify the differentially expressed anthocyanin biosynthetic genes between red- and white-skinned radishes and understand the molecular regulatory mechanism underlying red-skinned radish formation.

METHODS

Based on the published complete genome sequence of radish, the digital gene expression profiles of Yangzhouyuanbai (YB, white-skinned) and Sading (SD, red-skinned) were analyzed using Illumina sequencing.

RESULTS

A total of 3666 DEGs were identified in SD compared with YB. Interestingly, 46 genes encoded enzymes related to anthocyanin biosynthesis and 241 genes encoded transcription factors were identified. KEGG pathway analysis showed that the formation of red-skinned radish was mainly controlled by pelargonidin-derived anthocyanin biosynthetic pathway genes. This process included the upregulation of PAL, C4H, 4CL, CHS, CHI, F3H, DFR, LDOX, and UGT enzymes in SD. CHS genes were specifically expressed in SD, and it might be the key point for red pigment accumulation in red-skinned radish. Furthermore, MYB1/2/75, bHLH (TT8), and WD 40 showed higher expression in SD than in YB. Meanwhile, the corresponding low-abundance anthocyanin biosynthesis enzymes and upregulation of MYB4 might be the factors influencing the formation of white-skinned radish.

CONCLUSION

These findings provide new insights into the molecular mechanisms and regulatory network of anthocyanin biosynthesis in red-skinned radish.

Fine-mapping of a major QTL (Fwr1) for fusarium wilt resistance in radish

A major radish QTL (Fwr1) for fusarium wilt resistance was fine-mapped. Sequence and expression analyses suggest that a gene encoding a serine/arginine-rich protein kinase is a candidate gene for Fwr1. Fusarium wilt resistance locus 1 (Fwr1) is a major quantitative trait locus (QTL) mediating the resistance of radish inbred line 'B2' to Fusarium oxysporum, which is responsible for fusarium wilt. We previously detected Fwr1 on radish linkage group 3 (i.e., chromosome 5). In this study, a high-resolution genetic map of the Fwr1 locus was constructed by analyzing 354 recombinant F₂ plants derived from a cross between 'B2' and '835', the latter of which is susceptible to fusarium wilt. The Fwr1 QTL was fine-mapped to a 139.8-kb region between markers FM82 and FM87 in the middle part of chromosome 5. Fifteen candidate genes were predicted in this region based on a sequence comparison with the 'WK10039' radish reference genome. Additionally, we examined the time-course expression patterns of these predicted genes following an infection by the fusarium wilt pathogen. The ORF4 expression level was significantly higher in the resistant 'B2' plants than in the susceptible '835' plants. The ORF4 sequence was predicted to encode a serine/arginine-rich protein kinase and includes SNPs that result in nonsynonymous mutations, which may have important functional consequences. This study reveals a novel gene responsible for fusarium wilt resistance in radish. Further analyses of this gene may elucidate the molecular mechanisms underlying the fusarium wilt resistance of plants.

An ultra-high-density genetic map provides insights into genome synteny, recombination landscape and taproot skin colour in radish (Raphanus sativus L.)

High-density genetic map is a valuable tool for exploring novel genomic information, quantitative trait locus (QTL) mapping and gene discovery of economically agronomic traits in plant species. However, high-resolution genetic map applied to tag QTLs associated with important traits and to investigate genomic features underlying recombination landscape in radish (Raphanus sativus) remains largely unexplored. In this study, an ultra-high-density genetic map with 378 738 SNPs covering 1306.8 cM in nine radish linkage groups (LGs) was developed by a whole-genome sequencing-based approach. A total of 18 QTLs for 11 horticulture traits were detected. The map-based cloning data indicated that the R2R3-MYB transcription factor RsMYB90 was a crucial candidate gene determining the taproot skin colour. Comparative genomics analysis among radish, Brassica rapa and B. oleracea genome revealed several genomic rearrangements existed in the radish genome. The highly uneven distribution of recombination was observed across the nine radish chromosomes. Totally, 504 recombination hot regions (RHRs) were enriched near gene promoters and terminators. The recombination rate in RHRs was positively correlated with the density of SNPs and gene, and GC content, respectively. Functional annotation indicated that genes within RHRs were mainly involved in metabolic process and binding. Three QTLs for three traits were found in the RHRs. The results provide novel insights into the radish genome evolution and recombination landscape, and facilitate the development of effective strategies for molecular breeding by targeting and dissecting important traits in radish.

Genome- and Transcriptome-Wide Characterization of bZIP Gene Family Identifies Potential Members Involved in Abiotic Stress Response and Anthocyanin Biosynthesis in Radish (Raphanus sativus L.)

Basic leucine zipper (bZIP) transcription factors play crucial roles in various abiotic stress responses as well as anthocyanin accumulation. Anthocyanins are most abundant in colorful skin radish, which exhibit strong antioxidant activity that offers benefits for human health. Here, a total of 135 bZIP-encoding genes were identified from radish genome. Synteny analysis showed that 104 radish and 63 ArabidopsisbZIP genes were orthologous. Transcriptome analysis revealed that 10 RsbZIP genes exhibited high-expression levels in radish taproot (RPKM>10). Specifically, RsbZIP010 exhibited down-regulated expression under Cd, Cr and Pb stresses, whereas RsbZIP031 and RsbZIP059 showed significant down-regulation under heat and salt stresses, respectively. RT-qPCR analysis indicated that RsbZIP011 and RsbZIP102 were significantly up-regulated in the tissues of radish with high anthocyanin contents. Furthermore, the promoter sequences of 39 anthocyanin-related genes were found to contain G-box or ACE-box elements that could be recognized by bZIP family members. Taken together, several RsbZIPs might be served as critical regulators in radish taproot under Cd, Cr, Pb, heat and salt stresses. RsbZIP011 and RsbZIP102 were the potential participants in anthocyanin biosynthesis pathway of radish. These results facilitate further investigation on functional characterization of bZIP genes in response to abiotic stress and anthocyanin biosynthesis in radish.

Early Response of Radish to Heat Stress by Strand-Specific Transcriptome and miRNA Analysis

Radish is a crucial vegetable crop of the Brassicaceae family with many varieties and large cultivated area in China. Radish is a cool season crop, and there are only a few heat tolerant radish varieties in practical production with little information concerning the related genes in response to heat stress. In this work, some physiological parameter changes of young leaves under short-term heat stress were detected. Furthermore, we acquired 1802 differentially expressed mRNAs (including encoding some heat shock proteins, heat shock factor and heat shock-related transcription factors), 169 differentially expressed lncRNAs and three differentially expressed circRNAs (novel_circ_0000265, novel_circ_0000325 and novel_circ_0000315) through strand-specific RNA sequencing technology. We also found 10 differentially expressed miRNAs (ath-miR159b-3p, athmiR159c, ath-miR398a-3p, athmiR398b-3p, ath-miR165a-5p, ath-miR169g-3p, novel_86, novel_107, novel_21 and ath-miR171b-3p) by small RNA sequencing technology. Through function prediction and enrichment analysis, our results suggested that the significantly possible pathways/complexes related to heat stress in radish leaves were circadian rhythm-plant, photosynthesis-antenna proteins, photosynthesis, carbon fixation in photosynthetic organisms, arginine and proline metabolism, oxidative phosphorylation, peroxisome and plant hormone signal transduction. Besides, we identified one lncRNA-miRNA-mRNAs combination responsive to heat stress. These results will be helpful for further illustration of molecular regulation networks of how radish responds to heat stress.

Mobile RNAs and proteins: Prospects in storage organ development of tuber and root crops

Storage tuber and root crops make up a significant portion of the world's subsistence food supply. Because of their importance in food security, yield enhancement has become a priority. A major focus has been to understand the biology of belowground storage organ development. Considerable insights have been gained studying tuber development in potato. We now know that two mobile signals, a full-length mRNA, StBEL5, and a protein, StSP6A, play pivotal roles in regulating tuber development. Under favorable conditions, these signals move from leaves to a belowground modified stem (stolon) and regulate genes that activate tuberization. Overexpression of StBEL5 or StSP6A increases tuber yield even under non-inductive conditions. The mRNAs of two close homologs of StBEL5, StBEL11 and StBEL29, are also known to be mobile but act as repressors of tuberization. Polypyrimidine tract-binding proteins (PTBs) are RNA-binding proteins that facilitate the movement of these mRNAs. Considering their role in tuberization, it is possible that these mobile signals play a major role in storage root development as well. In this review, we explore the presence of these signals and their relevance in the development and yield potential of several important storage root crops.

What makes turnips: anatomy, physiology and transcriptome during early stages of its hypocotyl-tuber development

Brassica species are characterized by their tremendous intraspecific diversity, exemplified by leafy vegetables, oilseeds, and crops with enlarged inflorescences or above ground storage organs. In contrast to potato tubers that are edible storage organs storing energy as starch and are the vegetative propagation modules, the storage organs of turnips, grown from true seed, are swollen hypocotyls with varying degrees of root and stem that mainly store glucose and fructose. To highlight their anatomical origin, we use the term "hypocotyl-tuber" for these turnip vegetative storage organs. We combined cytological, physiological, genetic and transcriptomic approaches, aiming to identify the initial stages, molecular pathways and regulatory genes for hypocotyl-tuber induction in turnips (B. rapa subsp. rapa). We first studied the development of the hypocotyl zone of turnip and Pak choi and found that 16 days after sowing (DAS) morphological changes occurred in the xylem which indicated the early tuberization stage. Tissue culture experiments showed a clear effect of auxin on hypocotyl-tuber growth. Differentially expressed genes between 1 and 6 weeks after sowing in turnip hypocotyls, located in genomic regions involved in tuber initiation and/or tuber growth defined by QTL and selective sweeps for tuber formation, were identified as candidate genes that were studied in more detail for their role in hypocotyl-tuber formation. This included a Bra-FLOR1 paralogue with increased expression 16 DAS, when the hypocotyl starts swelling, suggesting dual roles for duplicated flowering time genes in flowering and hypocotyl-tuber induction. Bra-CYP735A2 was identified for its possible role in tuber growth via trans-zeatin. Weigthed Co-expression Network Analysis (WGCNA) identified 59 modules of co-expressed genes. Bra-FLOR1 and Bra-CYP735A2 were grouped in a module that included several genes involved in carbohydrate transport and metabolism, cell-wall growth, auxin regulation and secondary metabolism that serve as starting points to illuminate the transcriptional regulation of hypocotyl-tuber formation and development.