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

Hydrogen nanobubbles and oxidative windows: Investigating how varying hydrogen concentrations influence seed germination and stress defense

The hydrogen molecule can effectively regulate plant growth and development, improving plant resistance to abiotic stresses. However, studies regarding the optimal concentration of hydrogen and the associated mechanisms of action in organisms are lacking. This study showed that the maximum germination rate of radish seeds decreased from 90 % to 50 % under the stress of cadmium ions (Cd2+), and hydrogen nanobubble (NB) water significantly alleviated the stress effect of Cd2+ on radish seed germination. A hydrogen concentration of 0.8 ppm had the best effect, reducing Cd2+ accumulation in radish seeds by 63.23 % and increasing the maximum germination rate from 50 % to 65 %. At concentrations exceeding 1.2 ppm, the beneficial effect of hydrogen was weakened or even reversed. Consequently, we integrated the concept of the oxidative window into a REDOX balance model and demonstrated that an appropriate hydrogen concentration can effectively maintain the REDOX state within organisms. Transcriptome sequencing analysis revealed that hydrogen NB water modulated Cd2+ absorption and accumulation in seeds by regulating cell wall components, alleviating oxidative stress through oxidoreductase activity, and enhancing nutrient synthesis and metabolism. This collectively alleviated the inhibitory effect of Cd2+ on seed germination. This study is helpful for further understanding the effect of hydrogen concentration on the REDOX balance of seed germination, providing a theoretical basis for selecting hydrogen concentration to improve its effectiveness in agricultural fields.

Comparative transcriptome analysis reveals transcriptional regulation of anthocyanin biosynthesis in purple radish (Raphanus sativus L.)

Radish exhibits significant variation in color, particularly in sprouts, leaves, petals, fleshy roots, and other tissues, displaying a range of hues such as green, white, red, purple, and black. Although extensive research has been conducted on the color variation of radish, the underlying mechanism behind the variation in radish flower color remains unclear. To date, there is a lack of comprehensive research investigating the variation mechanism of radish sprouts, leaves, fleshy roots, and flower organs. This study aims to address this gap by utilizing transcriptome sequencing to acquire transcriptome data for white and purple radish flowers. Additionally, the published transcriptome data of sprouts, leaves, and fleshy roots were incorporated to conduct a systematic analysis of the regulatory mechanisms underlying anthocyanin biosynthesis in these four radish tissues. The comparative transcriptome analysis revealed differential expression of the anthocyanin biosynthetic pathway genes DFR, UGT78D2, TT12 and CPC in the four radish tissues. Additionally, the WGCNA results identified RsDFR.9c and RsUGT78D2.2c as hub genes responsible for regulating anthocyanin biosynthesis. By integrating the findings from the comparative transcriptome analysis, WGCNA, and anthocyanin biosynthetic pathway-related gene expression patterns, it is hypothesized that genes RsDFR.9c and RsUGT78D2.2c may serve as pivotal regulators of anthocyanins in the four radish tissues. Furthermore, the tissue-specific expression of the four copies of RsPAP1 is deemed crucial in governing anthocyanin synthesis and accumulation. Our results provide new insights into the molecular mechanism of anthocyanin biosynthesis and accumulation in different tissues of radish.

Whole-Genome Sequencing and Analysis of Tumour-Forming Radish (Raphanus sativus L.) Line

Spontaneous tumour formation in higher plants can occur in the absence of pathogen invasion, depending on the plant genotype. Spontaneous tumour formation on the taproots is consistently observed in certain inbred lines of radish (Raphanus sativus var. radicula Pers.). In this paper, using Oxford Nanopore and Illumina technologies, we have sequenced the genomes of two closely related radish inbred lines that differ in their ability to spontaneously form tumours. We identified a large number of single nucleotide variants (amino acid substitutions, insertions or deletions, SNVs) that are likely to be associated with the spontaneous tumour formation. Among the genes involved in the trait, we have identified those that regulate the cell cycle, meristem activity, gene expression, and metabolism and signalling of phytohormones. After identifying the SNVs, we performed Sanger sequencing of amplicons corresponding to SNV-containing regions to validate our results. We then checked for the presence of SNVs in other tumour lines of the radish genetic collection and found the ERF118 gene, which had the SNVs in the majority of tumour lines. Furthermore, we performed the identification of the CLAVATA3/ESR (CLE) and WUSCHEL (WOX) genes and, as a result, identified two unique radish CLE genes which probably encode proteins with multiple CLE domains. The results obtained provide a basis for investigating the mechanisms of plant tumour formation and also for future genetic and genomic studies of radish.

Degenerate oligonucleotide primer MIG-seq: an effective PCR-based method for high-throughput genotyping

Next-generation sequencing (NGS) library construction often involves using restriction enzymes to decrease genome complexity, enabling versatile polymorphism detection in plants. However, plant leaves frequently contain impurities, such as polyphenols, necessitating DNA purification before enzymatic reactions. To overcome this problem, we developed a PCR-based method for expeditious NGS library preparation, offering flexibility in number of detected polymorphisms. By substituting a segment of the simple sequence repeat sequence in the MIG-seq primer set (MIG-seq being a PCR method enabling library construction with low-quality DNA) with degenerate oligonucleotides, we introduced variability in detectable polymorphisms across various crops. This innovation, named degenerate oligonucleotide primer MIG-seq (dpMIG-seq), enabled a streamlined protocol for constructing dpMIG-seq libraries from unpurified DNA, which was implemented stably in several crop species, including fruit trees. Furthermore, dpMIG-seq facilitated efficient lineage selection in wheat and enabled linkage map construction and quantitative trait loci analysis in tomato, rice, and soybean without necessitating DNA concentration adjustments. These findings underscore the potential of the dpMIG-seq protocol for advancing genetic analyses across diverse plant species.

Total mRNA sequence dataset from Pectobacterium atrosepticum colonising potato or radish roots

Pectobacterium atrosepticum (Pba) is a gram-negative bacterium that causes blackleg and tuber soft rot of potato but can also asymptomatically colonise other (non-host) plant species. The aim of this study was to investigate the molecular processes and responses involved in Pba-host (potato) and Pba-non-host (radish) interactions, under laboratory conditions. To achieve this, we used total mRNA-sequencing to measure the gene expression patterns from all three species: Pba, potato and radish. We employed an end-point dual transcriptome approach. We used hydroponically grown potato (Solanum tuberosum var. Estima) and oil radish (Raphanus sativa var. Bento) roots inoculated with Pba SCRI1039 for 14 days compared to un-inoculated control plants or cultured bacteria. Total RNA was extracted from replicates of the two plant species and the bacterium using a Macherey-Nagel Nucleospin Plant RNA kit. The RNA from the 17 samples was then subjected to total mRNA-sequencing (paired-end) on an Illumina NovaSeq 6000™ sequencing platform. This gave between 39.2-58.1M reads per sample. The high-quality reads obtained were mapped to the corresponding reference genomes using Bowtie2 and the percentages of bacterium and plant transcripts calculated. This dataset constitutes the raw read fastq files and can be used to inform on genes active in plant rhizosphere-microbe interactions.

Current status of community resources and priorities for weed genomics research

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Sequence and epigenetic variations of R2R3-MYB transcription factors determine the diversity of taproot skin and flesh colors in different cultivated types of radish (Raphanus sativus L.)

KEY MESSAGE

This study found that three paralogous R2R3-MYB transcription factors exhibit functional divergence among different subspecies and cultivated types in radish. Cultivated radish taproots exhibit a wide range of color variations due to unique anthocyanin accumulation patterns in various tissues. This study investigated the universal principles of taproot color regulation that developed during domestication of different subspecies and cultivated types. The key candidate genes RsMYB1 and RsMYB2, which control anthocyanin accumulation in radish taproots, were identified using bulked segregant analysis in two genetic populations. We introduced the RsMYB1-RsF3'H-RsMYB1Met genetic model to elucidate the complex and unstable genetic regulation of taproot flesh color in Xinlimei radish. Furthermore, we analyzed the expression patterns of three R2R3-MYB transcription factors in lines with different taproot colors and investigated the relationship between RsMYB haplotypes and anthocyanin accumulation in a natural population of 56 germplasms. The results revealed that three paralogous RsMYBs underwent functional divergence during radish domestication, with RsMYB1 regulating the red flesh of Xinlimei radish, and RsMYB2 and RsMYB3 regulating the red skin of East Asian big long radish (R. sativus var. hortensis) and European small radish (R. sativus var. sativus), respectively. Moreover, RsMYB1-H1, RsMYB2-H10, and RsMYB3-H6 were identified as the primary haplotypes exerting regulatory functions on anthocyanin synthesis. These findings provide an understanding of the genetic mechanisms regulating anthocyanin synthesis in radish and offer a potential strategy for early prediction of color variations in breeding programs.

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Technology-enabled great leap in deciphering plant genomes

Plant genomes provide essential and vital basic resources for studying many aspects of plant biology and applications (for example, breeding). From 2000 to 2020, 1,144 genomes of 782 plant species were sequenced. In the past three years (2021-2023), 2,373 genomes of 1,031 plant species, including 793 newly sequenced species, have been assembled, representing a great leap. The 2,373 newly assembled genomes, of which 63 are telomere-to-telomere assemblies and 921 have been generated in pan-genome projects, cover the major phylogenetic clades. Substantial advances in read length, throughput, accuracy and cost-effectiveness have notably simplified the achievement of high-quality assemblies. Moreover, the development of multiple software tools using different algorithms offers the opportunity to generate more complete and complex assemblies. A database named N3: plants, genomes, technologies has been developed to accommodate the metadata associated with the 3,517 genomes that have been sequenced from 1,575 plant species since 2000. We also provide an outlook for emerging opportunities in plant genome sequencing.

A comprehensive evaluation of the potential of three next-generation short-read-based plant pan-genome construction strategies for the identification of novel non-reference sequence

Pan-genome studies are important for understanding plant evolution and guiding the breeding of crops by containing all genomic diversity of a certain species. Three short-read-based strategies for plant pan-genome construction include iterative individual, iteration pooling, and map-to-pan. Their performance is very different under various conditions, while comprehensive evaluations have yet to be conducted nowadays. Here, we evaluate the performance of these three pan-genome construction strategies for plants under different sequencing depths and sample sizes. Also, we indicate the influence of length and repeat content percentage of novel sequences on three pan-genome construction strategies. Besides, we compare the computational resource consumption among the three strategies. Our findings indicate that map-to-pan has the greatest recall but the lowest precision. In contrast, both two iterative strategies have superior precision but lower recall. Factors of sample numbers, novel sequence length, and the percentage of novel sequences' repeat content adversely affect the performance of all three strategies. Increased sequencing depth improves map-to-pan's performance, while not affecting the other two iterative strategies. For computational resource consumption, map-to-pan demands considerably more than the other two iterative strategies. Overall, the iterative strategy, especially the iterative pooling strategy, is optimal when the sequencing depth is less than 20X. Map-to-pan is preferable when the sequencing depth exceeds 20X despite its higher computational resource consumption.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

A graph-based pan-genome of Brassica oleracea provides new insights into its domestication and morphotype diversification

The domestication of Brassica oleracea has resulted in diverse morphological types with distinct patterns of organ development. Here we report a graph-based pan-genome of B. oleracea constructed from high-quality genome assemblies of different morphotypes. The pan-genome harbors over 200 structural variant hotspot regions enriched in auxin- and flowering-related genes. Population genomic analyses revealed that early domestication of B. oleracea focused on leaf or stem development. Gene flows resulting from agricultural practices and variety improvement were detected among different morphotypes. Selective-sweep and pan-genome analyses identified an auxin-responsive small auxin up-regulated RNA gene and a CLAVATA3/ESR-RELATED family gene as crucial players in leaf-stem differentiation during the early stage of B. oleracea domestication and the BoKAN1 gene as instrumental in shaping the leafy heads of cabbage and Brussels sprouts. Our pan-genome and functional analyses further revealed that variations in the BoFLC2 gene play key roles in the divergence of vernalization and flowering characteristics among different morphotypes, and variations in the first intron of BoFLC3 are involved in fine-tuning the flowering process in cauliflower. This study provides a comprehensive understanding of the pan-genome of B. oleracea and sheds light on the domestication and differential organ development of this globally important crop species.

Construction of SNP fingerprints and genetic diversity analysis of radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is a vegetable crop with economic value and ecological significance in the genus Radish, family Brassicaceae. In recent years, developed countries have attached great importance to the collection and conservation of radish germplasm resources and their research and utilization, but the lack of population genetic information and molecular markers has hindered the development of the genetic breeding of radish. In this study, we integrated the radish genomic data published in databases for the development of single-nucleotide polymorphism (SNP) markers, and obtained a dataset of 308 high-quality SNPs under strict selection criteria. With the support of Kompetitive Allele-Specific PCR (KASP) technology, we screened a set of 32 candidate core SNP marker sets to analyse the genetic diversity of the collected 356 radish varieties. The results showed that the mean values of polymorphism information content (PIC), minor allele frequency (MAF), gene diversity and heterozygosity of the 32 candidate core SNP markers were 0.32, 0.30, 0.40 and 0.25, respectively. Population structural analysis, principal component analysis and genetic evolutionary tree analysis indicated that the 356 radish materials were best classified into two taxa, and that the two taxa of the material were closely genetically exchanged. Finally, on the basis of 32 candidate core SNP markers we calculated 15 core markers using a computer algorithm to construct a fingerprint map of 356 radish varieties. Furthermore, we constructed a core germplasm population consisting of 71 radish materials using 32 candidate core markers. In this study, we developed SNP markers for radish cultivar identification and genetic diversity analysis, and constructed DNA fingerprints, providing a basis for the identification of radish germplasm resources and molecular marker-assisted breeding as well as genetic research.

Effect of structural variation in the promoter region of RsMYB1.1 on the skin color of radish taproot

Accumulation of anthocyanins in the taproot of radish is an agronomic trait beneficial for human health. Several genetic loci are related to a red skin or flesh color of radish, however, the functional divergence of candidate genes between non-red and red radishes has not been investigated. Here, we report that a novel genetic locus on the R2 chromosome, where RsMYB1.1 is located, is associated with the red color of the skin of radish taproot. A genome-wide association study (GWAS) of 66 non-red-skinned (nR) and 34 red-skinned (R) radish accessions identified three nonsynonymous single nucleotide polymorphisms (SNPs) in the third exon of RsMYB1.1. Although the genotypes of SNP loci differed between the nR and R radishes, no functional difference in the RsMYB1.1 proteins of nR and R radishes in their physical interaction with RsTT8 was detected by yeast-two hybrid assay or in anthocyanin accumulation in tobacco and radish leaves coexpressing RsMYB1.1 and RsTT8. By contrast, insertion- or deletion-based GWAS revealed that one large AT-rich low-complexity sequence of 1.3-2 kb was inserted in the promoter region of RsMYB1.1 in the nR radishes (RsMYB1.1nR), whereas the R radishes had no such insertion; this represents a presence/absence variation (PAV). This insertion sequence (RsIS) was radish specific and distributed among the nine chromosomes of Raphanus genomes. Despite the extremely low transcription level of RsMYB1.1nR in the nR radishes, the inactive RsMYB1.1nR promoter could be functionally restored by deletion of the RsIS. The results of a transient expression assay using radish root sections suggested that the RsIS negatively regulates the expression of RsMYB1.1nR, resulting in the downregulation of anthocyanin biosynthesis genes, including RsCHS, RsDFR, and RsANS, in the nR radishes. This work provides the first evidence of the involvement of PAV in an agronomic trait of radish.

RsPDR8, a member of ABCG subfamily, plays a positive role in regulating cadmium efflux and tolerance in radish (Raphanus sativus L.)

Radish (Raphanus sativus L.) is one of the most vital root vegetable crops worldwide. Cadmium (Cd), a non-essential and toxic heavy metal, can dramatically restrict radish taproot quality and safety. Although the Peiotrpic Drug Resistance (PDR) genes play crucial roles in heavy metal accumulation and transport in plants, the systematic identification and functional characterization of RsPDRs remain largely unexplored in radish. Herein, a total of 19 RsPDR genes were identified from the radish genome. A few RsPDRs, including RsPDR1, RsPDR8 and RsPDR12, showed significant differential expression under Cd and lead (Pb) stress in the 'NAU-YH' genotype. Interestingly, the plasma membrane-localized RsPDR8 exhibited significantly up-regulated expression and enhanced promoter activity under Cd exposure. Ectopic expression of RsPDR8 conferred Cd tolerance via reducing Cd accumulation in yeast cells. Moreover, the transient transformation of RsPDR8 revealed that it positively regulated Cd tolerance by promoting ROS scavenging and enhancing membrane permeability in radish. In addition, overexpression of RsPDR8 increased root elongation but deceased Cd accumulation compared with the WT plants in Arabidopsis, demonstrating that it could play a positive role in mediating Cd efflux and tolerance in plants. Together, these results would facilitate deciphering the molecular mechanism underlying RsPDR8-mediated Cd tolerance and detoxification in radish.

RsCDF3, a member of Cycling Dof Factors, positively regulates cold tolerance via auto-regulation and repressing two RsRbohs transcription in radish (Raphanus sativus L.)

Radish is one of the most economical root vegetable crops worldwide. Cold stress dramatically impedes radish taproot formation and development as well as reduces its yield and quality. Although the Cycling Dof Factors (CDFs) play crucial roles in plant growth, development and abiotic stress responses, how CDF TFs mediate the regulatory network of cold stress response remains largely unexplored in radish. Herein, a total of nine RsCDF genes were identified from the radish genome. Among them, the RsCDF3 exhibited obviously up-regulated expression under cold stress, especially at 12 h and 24 h. RsCDF3 was localized to the nucleus and displayed dramatic cold-induced promoter activity in tobacco leaves. Moreover, overexpression of RsCDF3 significantly enhanced cold tolerance of radish plants, whereas its knock-down plants exhibited the opposite phenotype. Interestingly, both in vitro and in vivo assays indicated that the RsCDF3 repressed the transcription of RsRbohA and RsRbohC via directly binding to their promoters, which contributed to maintaining the cellular homeostasis of reactive oxygen species (ROS) production and scavenging in radish. In addition, the RsCDF3 bound to its own promoter to mediate its transcription, thereby forming an autoregulatory feedback loop to cooperatively trigger RsRbohs-dependent cold tolerance. Together, we revealed a novel RsCDF3-RsRbohs module to promote the cold tolerance in radish plants. These findings would facilitate unveiling the molecular mechanism governing RsCDF3-mediated cold stress response in radish.

Genome and haplotype provide insights into the population differentiation and breeding improvement of Gossypium barbadense

INTRODUCTION

Sea-island cotton (Gossypium barbadense, Gb) is one of the major sources of high-grade natural fiber. Besides the common annual Gb cotton, perennial Gb cotton is also cultivated, but studies on perennial Gb cotton are rare.

OBJECTIVES

We aimed to make a systematic analysis of perennial sea-island cotton and lay a foundation for its utilization in breeding, and try to identify the representative structural variations (SVs) in sea-island cotton, and to reveal the population differentiation and adaptive improvement of sea-island cotton.

METHODS

Through genome assembly of one perennial Gb cotton accession (named Gb_M210936) and comparative genome analysis, variations during Gb cotton domestication were identified by comparing Gb_M210936 with annual Gb accession 3-79 and with wild allotetraploid cotton G. darwinii. Six perennial Gb accessions combining with the resequenced 1,129 cotton accessions were used to conduct population and genetic analysis. Large haplotype blocks (haploblocks), generated from interspecific introgressions and intraspecific inversions, were identified and were used to analyze their effects on population differentiation and agronomic traits of sea-island cotton.

RESULTS

One reference genome of perennial sea-island cotton was assembled. Representative SVs in sea-island cotton were identified, and 31 SVs were found to be associated with agronomic traits. Perennial Gb cotton had a closer kinship with the wild-to-landrace continuum Gb cotton from south America where Gb cotton is originally domesticated. Haploblocks were associated with agronomic traits improvement of sea-island cotton, promoted sea-island cotton differentiation into three subgroups, were suffered from breeding selection, and may drive Gb cotton to be adapted to central Asian.

CONCLUSION

Our study made up the lack of perennial Gb cotton genome, and clarified that exotic introgressions improved the traits of sea-island cotton, promoted the population differentiation, and drove sea-island cotton adaptive to central Asia, which will provide new insights for the genetic breeding improvement of sea-island cottons.

Pan-genome analysis of 13 Malus accessions reveals structural and sequence variations associated with fruit traits

Structural variations (SVs) and copy number variations (CNVs) contribute to trait variations in fleshy-fruited species. Here, we assemble 10 genomes of genetically diverse Malus accessions, including the ever-green cultivar 'Granny Smith' and the widely cultivated cultivar 'Red Fuji'. Combining with three previously reported genomes, we assemble the pan-genome of Malus species and identify 20,220 CNVs and 317,393 SVs. We also observe CNVs that are positively correlated with expression levels of the genes they are associated with. Furthermore, we show that the noncoding RNA generated from a 209 bp insertion in the intron of mitogen-activated protein kinase homology encoding gene, MMK2, regulates the gene expression and affects fruit coloration. Moreover, we identify overlapping SVs associated with fruit quality and biotic resistance. This pan-genome uncovers possible contributions of CNVs to gene expression and highlights the role of SVs in apple domestication and economically important traits.

Telomere-to-telomere assembly of cassava genome reveals the evolution of cassava and divergence of allelic expression

Cassava is a crucial crop that makes a significant contribution to ensuring human food security. However, high-quality telomere-to-telomere cassava genomes have not been available up to now, which has restricted the progress of haploid molecular breeding for cassava. In this study, we constructed two nearly complete haploid resolved genomes and an integrated, telomere-to-telomere gap-free reference genome of an excellent cassava variety, 'Xinxuan 048', thereby providing a new high-quality genomic resource. Furthermore, the evolutionary history of several species within the Euphorbiaceae family was revealed. Through comparative analysis of haploid genomes, it was found that two haploid genomes had extensive differences in linear structure, transcriptome features, and epigenetic characteristics. Genes located within the highly divergent regions and differentially expressed alleles are enriched in the functions of auxin response and the starch synthesis pathway. The high heterozygosity of cassava 'Xinxuan 048' leads to rapid trait segregation in the first selfed generation. This study provides a theoretical basis and genomic resource for molecular breeding of cassava haploids.

Genome assembly and analysis of Lactuca virosa: implications for lettuce breeding

Lettuce (Lactuca sativa L.) is a leafy vegetable crop with ongoing breeding efforts related to quality, resilience, and innovative production systems. To breed resilient and resistant lettuce in the future, valuable genetic variation found in close relatives could be further exploited. Lactuca virosa (2x = 2n = 18), a wild relative assigned to the tertiary lettuce gene pool, has a much larger genome (3.7 Gbp) than Lactuca sativa (2.5 Gbp). It has been used in interspecific crosses and is a donor to modern crisphead lettuce cultivars. Here, we present a de novo reference assembly of L. virosa with high continuity and complete gene space. This assembly facilitated comparisons to the genome of L. sativa and to that of the wild species L. saligna, a representative of the secondary lettuce gene pool. To assess the diversity in gene content, we classified the genes of the 3 Lactuca species as core, accessory, and unique. In addition, we identified 3 interspecific chromosomal inversions compared to L. sativa, which each may cause recombination suppression and thus hamper future introgression breeding. Using 3-way comparisons in both reference-based and reference-free manners, we show that the proliferation of long-terminal repeat elements has driven the genome expansion of L. virosa. Further, we performed a genome-wide comparison of immune genes, nucleotide-binding leucine-rich repeat, and receptor-like kinases among Lactuca spp. and indicated the evolutionary patterns and mechanisms behind their expansions. These genome analyses greatly facilitate the understanding of genetic variation in L. virosa, which is beneficial for the breeding of improved lettuce varieties.

PanKmer: k-mer-based and reference-free pangenome analysis

SUMMARY

Pangenomes are replacing single reference genomes as the definitive representation of DNA sequence within a species or clade. Pangenome analysis predominantly leverages graph-based methods that require computationally intensive multiple genome alignments, do not scale to highly complex eukaryotic genomes, limit their scope to identifying structural variants (SVs), or incur bias by relying on a reference genome. Here, we present PanKmer, a toolkit designed for reference-free analysis of pangenome datasets consisting of dozens to thousands of individual genomes. PanKmer decomposes a set of input genomes into a table of observed k-mers and their presence-absence values in each genome. These are stored in an efficient k-mer index data format that encodes SNPs, INDELs, and SVs. It also includes functions for downstream analysis of the k-mer index, such as calculating sequence similarity statistics between individuals at whole-genome or local scales. For example, k-mers can be "anchored" in any individual genome to quantify sequence variability or conservation at a specific locus. This facilitates workflows with various biological applications, e.g. identifying cases of hybridization between plant species. PanKmer provides researchers with a valuable and convenient means to explore the full scope of genetic variation in a population, without reference bias.

AVAILABILITY AND IMPLEMENTATION

PanKmer is implemented as a Python package with components written in Rust, released under a BSD license. The source code is available from the Python Package Index (PyPI) at https://pypi.org/project/pankmer/ as well as Gitlab at https://gitlab.com/salk-tm/pankmer. Full documentation is available at https://salk-tm.gitlab.io/pankmer/.

Identification and expression analysis of the SWEET genes in radish reveal their potential functions in reproductive organ development

BACKGROUND

Sugars produced by photosynthesis provide energy for biological activities and the skeletons for macromolecules; they also perform multiple physiological functions in plants. Sugar transport across plasma membranes mediated by the Sugar Will Eventually be Exported Transporter (SWEET) genes substantially affects these processes. However, the evolutionary dynamics and function of the SWEET genes are largely unknown in radish, an important Brassicaceae species.

METHODS AND RESULTS

Genome-wide identification and analysis of the RsSWEET genes from the recently updated radish reference genome was conducted using bioinformatics methods. The tissue-specific expression was analyzed using public RNA-seq data, and the expression levels in the bud, stamens, pistils, pericarps and seeds at 15 and 30 days after flowering (DAF) were determined by RT‒qPCR. Thirty-seven RsSWEET genes were identified and named according to their Arabidopsis homologous. They are unevenly distributed across the nine radish chromosomes and were further divided into four clades by phylogenetic analysis. There are 5-7 transmembrane domains and at least one MtN3_slv domain in the RsSWEETs. RNA-seq and RT‒qPCR revealed that the RsSWEETs exhibit higher expression levels in the reproductive organs, indicating that these genes might play vital roles in reproductive organ development. RsSWEET15.1 was found to be especially expressed in siliques according to the RNA-seq data, and the RT‒qPCR results further confirmed that it was most highly expressed levels in the seeds at 30 DAF, followed by the pericarp at 15 DAF, indicating that it is involved in seed growth and development.

CONCLUSIONS

This study suggests that the RsSWEET genes play vital roles in reproductive organ development and provides a theoretical basis for the future functional analysis of RsSWEETs in radish.

The role of intraspecific crop-weed hybridization in the evolution of weediness and invasiveness: Cultivated and weedy radish (Raphanus sativus) as a case study

PREMISE

The phenotype of hybrids between a crop and its wild or weed counterpart is usually intermediate and maladapted compared to that of their parents; however, hybridization has sometimes been associated with increased fitness, potentially leading to enhanced weediness and invasiveness. Since the ecological context and maternal genetic effects may affect hybrid fitness, they could influence the evolutionary outcomes of hybridization. Here, we evaluated the performance of first-generation crop-weed hybrids of Raphanus sativus and their parents in two contrasting ecological conditions.

METHODS

Using experimental hybridization and outdoor common garden experiments, we assessed differences in time to flowering, survival to maturity, plant biomass, and reproductive components between bidirectional crop-weed hybrids and their parents in agrestal (wheat cultivation, fertilization, weeding) and ruderal (human-disturbed, uncultivated area) conditions over 2 years.

RESULTS

Crop, weeds, and bidirectional hybrids overlapped at least partially during the flowering period, indicating a high probability of gene flow. Hybrids survived to maturity at rates at least as successful as their parents and had higher plant biomass and fecundity, which resulted in higher fitness compared to their parents in both environments, without any differences associated with the direction of the hybridization.

CONCLUSIONS

Intraspecific crop-weed hybridization, regardless of the cross direction, has the potential to promote weediness in weedy R. sativus in agrestal and ruderal environments, increasing the chances for introgression of crop alleles into weed populations. This is the first report of intraspecific crop-weed hybridization in R. sativus.

Application of Multi-Perspectives in Tea Breeding and the Main Directions

Tea plants are an economically important crop and conducting research on tea breeding contributes to enhancing the yield and quality of tea leaves as well as breeding traits that satisfy the requirements of the public. This study reviews the current status of tea plants germplasm resources and their utilization, which has provided genetic material for the application of multi-omics, including genomics and transcriptomics in breeding. Various molecular markers for breeding were designed based on multi-omics, and available approaches in the direction of high yield, quality and resistance in tea plants breeding are proposed. Additionally, future breeding of tea plants based on single-cellomics, pangenomics, plant-microbe interactions and epigenetics are proposed and provided as references. This study aims to provide inspiration and guidance for advancing the development of genetic breeding in tea plants, as well as providing implications for breeding research in other crops.

Pan-genome of Citrullus genus highlights the extent of presence/absence variation during domestication and selection

The rich genetic diversity in Citrullus lanatus and the other six species in the Citrullus genus provides important sources in watermelon breeding. Here, we present the Citrullus genus pan-genome based on the 400 Citrullus genus resequencing data, showing that 477 Mb contigs and 6249 protein-coding genes were absent in the Citrullus lanatus reference genome. In the Citrullus genus pan-genome, there are a total of 8795 (30.5%) genes that exhibit presence/absence variations (PAVs). Presence/absence variation (PAV) analysis showed that a lot of gene PAV were selected during the domestication and improvement, such as 53 favorable genes and 40 unfavorable genes were identified during the C. mucosospermus to C. lanatus landrace domestication. We also identified 661 resistance gene analogs (RGAs) in the Citrullus genus pan-genome, which contains 90 RGAs (89 variable and 1 core gene) located on the pangenome additional contigs. By gene PAV-based GWAS, 8 gene presence/absence variations were found associated with flesh color. Finally, based on the results of gene PAV selection analysis between watermelon populations with different fruit colors, we identified four non-reference candidate genes associated with carotenoid accumulation, which had a significantly higher frequency in the white flesh. These results will provide an important source for watermelon breeding.

A review of the pangenome: how it affects our understanding of genomic variation, selection and breeding in domestic animals?

As large-scale genomic studies have progressed, it has been revealed that a single reference genome pattern cannot represent genetic diversity at the species level. While domestic animals tend to have complex routes of origin and migration, suggesting a possible omission of some population-specific sequences in the current reference genome. Conversely, the pangenome is a collection of all DNA sequences of a species that contains sequences shared by all individuals (core genome) and is also able to display sequence information unique to each individual (variable genome). The progress of pangenome research in humans, plants and domestic animals has proved that the missing genetic components and the identification of large structural variants (SVs) can be explored through pangenomic studies. Many individual specific sequences have been shown to be related to biological adaptability, phenotype and important economic traits. The maturity of technologies and methods such as third-generation sequencing, Telomere-to-telomere genomes, graphic genomes, and reference-free assembly will further promote the development of pangenome. In the future, pangenome combined with long-read data and multi-omics will help to resolve large SVs and their relationship with the main economic traits of interest in domesticated animals, providing better insights into animal domestication, evolution and breeding. In this review, we mainly discuss how pangenome analysis reveals genetic variations in domestic animals (sheep, cattle, pigs, chickens) and their impacts on phenotypes and how this can contribute to the understanding of species diversity. Additionally, we also go through potential issues and the future perspectives of pangenome research in livestock and poultry.

Genome-wide characterization of RsHSP70 gene family reveals positive role of RsHSP70-20 gene in heat stress response in radish (Raphanus sativus L.)

Radish is an economical cool-season root vegetable crop worldwide. Heat shock protein 70 (HSP70) plays indispensable roles in plant growth, development and abiotic stress responses. Nevertheless, little information is available regarding the identification and functional characterization of HSP70 gene family in radish. Herein, a total of 34 RsHSP70 genes were identified at the radish genome level, among which nine and 25 RsHSP70s were classified into the HSP110/SSE and DnaK subfamilies, respectively. RNA-seq analysis revealed that some RsHSP70 genes had differential expression profile in radish leaf, root, stamen and pistil. A range of RsHSP70 genes exhibited differential expression under several abiotic stresses such as heat, salt and heavy metals. Intriguingly, the expression of four RsHSP70 genes (RsHSP70-7, RsHSP70-12, RsHSP70-20 and RsHSP70-22) was dramatically up-regulated under heat stress (HS). RT-qPCR and transient LUC reporter assay indicated that both the expression and promoter activity of RsHSP70-20 was strongly induced by HS. Notably, overexpression of RsHSP70-20 significantly enhanced thermotolerance by decreasing reactive oxygen species and promoting proline accumulation in radish, whereas its knock-down plants exhibited increased thermosensitivity, indicating that RsHSP70-20 positively regulate HS response in radish. These results would provide valuable information to decipher the molecular basis of RsHSP70-mediated thermotolerance in radish.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

The Disturbance of the Antioxidant System Results in Internal Blue Discoloration of Postharvest Cherry Radish (Raphanus sativus L. var. radculus pers) Roots

Internal blue discoloration in cherry radish (Raphanus sativus L. var. radculus pers) roots can appear after harvest. The antioxidant system and content of reactive oxygen species (ROS) will affect the blue discoloration. Currently, the reason for the blue discoloration is not yet clear. In order to reveal the mechanism of the blue discoloration of cherry radish, we selected the blue discolored cherry radish as the research object and the white cherry radish as the control. The difference in the antioxidant system between them were compared, including related enzymes and non-enzymatic antioxidants in this system. Meanwhile, the changes in the contents of 4-hydroxyglucobrassicin as a precursor substance and ROS were compared. The results showed that the activities of typical antioxidant enzymes decreased and the cycle of Glutathione peroxidase (GPX) and Ascorbic acid-Glutathione (ASA-GSH) was disturbed, leading to the reduction of antioxidant effect and the failure of timely and effective decomposition of superoxide anions (O₂^•-) and hydrogen peroxide (H₂O₂). In addition, the elevated level of O₂^•- and H₂O₂ led to the disorder of the antioxidant system, while the 4-hydroxybrassinoside was oxidized under the catalysis of peroxidase (POD) and eventually led to the internal blue discoloration in cherry radish. These results can provide a theoretical basis for solving the blue discoloration problem.

Combined widely targeted metabolomics and transcriptomics analysis reveals differentially accumulated metabolites and the underlying molecular bases in fleshy taproots of distinct radish genotypes

Radish is an important taproot crop with medicinal and edible uses that is cultivated worldwide. However, the differences in metabolites and the underlying molecular bases among different radish types remain largely unknown. In the present study, a combined analysis of liquid chromatography-electrospray ionization-tandem mass spectrometry (LC-ESI-MS/MS) and RNA-Seq data was conducted to uncover important differentially accumulated metabolites (DAMs) among radish accessions with green, white and red taproot flesh colours. A total of 657 metabolites were identified and 138 DAMs were commonly present in red vs. green and red vs. white accessions. Red accessions were rich in phenolic compounds, while green and white accessions had more amino acids. Additionally, 41 metabolites and 98 genes encoding 37 enzymes were enriched in the shikimate and phenolic biosynthesis pathways. CHS is the rate-limiting enzyme determining flavonoid differences among accessions. A total of 119 candidate genes might contribute to red accession-specific accumulated metabolites. Specifically, one gene cluster consisting of 16 genes, including one RsMYB1, RsMYBL2, RsTT8, RsDFR, RsANS, Rs4CL3, RsSCPL10, Rs3AT1 and RsSAP2 gene, two RsTT19 and RsWRKY44 genes and three RsUGT genes, might be involved in anthocyanin accumulation in red radish fleshy taproots. More importantly, an InDel marker was developed based on an RsMYB1 promoter mutation, and the accuracy reached 95.9% when it was used to select red-fleshed radishes. This study provides comprehensive insights into the metabolite differences and underlying molecular mechanisms in fleshy taproots among different radish genotypes and will be beneficial for the genetic improvement of radish nutritional quality.

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.

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 (F_ST) 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.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

Sinapis genomes provide insights into whole-genome triplication and divergence patterns within tribe Brassiceae

Sinapis alba and Sinapis arvensis are mustard crops within the Brassiceae tribe of the Brassicaceae family, and represent an important genetic resource for crop improvement. We performed the de novo assembly of Brassica nigra, S. alba, and S. arvensis, and conducted comparative genomics to investigate the pattern of genomic evolution since an ancient whole-genome triplication event. Both Sinapis species retained evidence of the Brassiceae whole-genome triplication approximately 20.5 million years ago (Mya), with subgenome dominance observed in gene density, gene expression, and selective constraint. While S. alba diverged from the ancestor of Brassica and Raphanus at approximately 12.5 Mya, the divergence time of S. arvensis and B. nigra was approximately 6.5 Mya. S. arvensis and B. nigra had greater collinearity compared with their relationship to either Brassica rapa or Brassica oleracea. Two chromosomes of S. alba (Sal03 and Sal08) were completely collinear with two ancestral chromosomes proposed in the Ancestral Crucifer Karyotype (ACK) genomic block model, the first time this has been observed in the Brassiceae. These results are consistent with S. alba representing a relatively ancient lineage of the species evolved from the common ancestor of tribe Brassiceae, and suggest that the phylogeny of the Brassica and Sinapis genera requires some revision. Our study provides new insights into the genome evolution and phylogenetic relationships of Brassiceae and provides genomic information for genetic improvement of these plants.

Characterization and acceleration of genome shuffling and ploidy reduction in synthetic allopolyploids by genome sequencing and editing

Polyploidy and the subsequent ploidy reduction and genome shuffling are the major driving forces of genome evolution. Here, we revealed short-term allopolyploid genome evolution by sequencing a synthetic intergeneric hybrid (Raphanobrassica, RRCC). In this allotetraploid, the genome deletion was quick, while rearrangement was slow. The core and high-frequency genes tended to be retained while the specific and low-frequency genes tended to be deleted in the hybrid. The large-fragment deletions were enriched in the heterochromatin region and probably derived from chromosome breaks. The intergeneric translocations were primarily of short fragments dependent on homoeology, indicating a gene conversion origin. To accelerate genome shuffling, we developed an efficient genome editing platform for Raphanobrassica. By editing Fanconi Anemia Complementation Group M (FANCM) genes, homoeologous recombination, chromosome deletion and secondary meiosis with additional ploidy reduction were accelerated. FANCM was shown to be a checkpoint of meiosis and controller of ploidy stability. By simultaneously editing FLIP genes, gene conversion was precisely introduced, and mosaic genes were produced around the target site. This intergeneric hybrid and genome editing platform not only provides models that facilitate experimental evolution research by speeding up genome shuffling and conversion but also accelerates plant breeding by enhancing intergeneric genetic exchange and creating new genes.

Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium

Phenotypic diversity and evolutionary innovation ultimately trace to variation in genomic sequence and rewiring of regulatory networks. Here, we constructed a pan-genome of the Gossypium genus using ten representative diploid genomes. We document the genomic evolutionary history and the impact of lineage-specific transposon amplification on differential genome composition. The pan-3D genome reveals evolutionary connections between transposon-driven genome size variation and both higher-order chromatin structure reorganization and the rewiring of chromatin interactome. We linked changes in chromatin structures to phenotypic differences in cotton fiber and identified regulatory variations that decode the genetic basis of fiber length, the latter enabled by sequencing 1,005 transcriptomes during fiber development. We showcase how pan-genomic, pan-3D genomic and genetic regulatory data serve as a resource for delineating the evolutionary basis of spinnable cotton fiber. Our work provides insights into the evolution of genome organization and regulation and will inform cotton improvement by enabling regulome-based approaches.

Integration of metabolome and transcriptome analyses reveals the mechanism of anthocyanin accumulation in purple radish leaves

Anthocyanins are natural pigments and play significant roles in multiple growth, development, and stress response processes in plants. The vegetables with high anthocyanin content have better colours, higher antioxidant activity than green vegetables and are potent antioxidants with health benefits. However, the mechanism of anthocyanin accumulation in purple and green leaves of Raphanus sativus (radish) is poorly understood and needs further investigation. In the present study, the pigment content in a green leaf cultivar "RA9" and a purple-leaf cultivar "MU17" was characterized and revealed that the MU17 had significantly increased accumulation of anthocyanins and reduced content of chlorophyll and carotenoid compared with that in RA9. Meanwhile, these two cultivars were subjected to a combination of metabolomic and transcriptome studies. A total of 52 massively content-changed metabolites and 3463 differentially expressed genes were discovered in MU17 compared with RA9. In addition, the content of significantly increased flavonoids (such as pelargonidin and cyanidin) was identified in MU17 compared to RA9 using an integrated analysis of metabolic and transcriptome data. Moreover, the quantitative real-time polymerase chain reaction results also confirmed the differences in the expression of genes related to pathways of flavonoids and anthocyanin metabolism in MU17 leaves. The present findings provide valuable information for anthocyanin metabolism and further genetic manipulation of anthocyanin biosynthesis in radish leaves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01245-w.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Gambierdiscus and Its Associated Toxins: A Minireview

Gambierdiscus is a dinoflagellate genus widely distributed throughout tropical and subtropical regions. Some members of this genus can produce a group of potent polycyclic polyether neurotoxins responsible for ciguatera fish poisoning (CFP), one of the most significant food-borne illnesses associated with fish consumption. Ciguatoxins and maitotoxins, the two major toxins produced by Gambierdiscus, act on voltage-gated channels and TRPA1 receptors, consequently leading to poisoning and even death in both humans and animals. Over the past few decades, the occurrence and geographic distribution of CFP have undergone a significant expansion due to intensive anthropogenic activities and global climate change, which results in more human illness, a greater public health impact, and larger economic losses. The global spread of CFP has led to Gambierdiscus and its toxins being considered an environmental and human health concern worldwide. In this review, we seek to provide an overview of recent advances in the field of Gambierdiscus and its associated toxins based on the existing literature combined with re-analyses of current data. The taxonomy, phylogenetics, geographic distribution, environmental regulation, toxin detection method, toxin biosynthesis, and pharmacology and toxicology of Gambierdiscus are summarized and discussed. We also highlight future perspectives on Gambierdiscus and its associated toxins.

Genome-wide identification of sugar transporter gene family in Brassicaceae crops and an expression analysis in the radish

BACKGROUND

Sugar not only is an important biomacromolecule that plays important roles in plant growth, development, and biotic and abiotic stress tolerance but also provides a skeleton for other macromolecules, such as proteins and nucleic acids. Sugar transporter proteins (STPs) play essential roles in plant sugar transport and ultimately affect the abovementioned life processes. However, the evolutionary dynamics of this important gene family in Brassicaceae crops are still largely unknown, and the functional differentiation of radish STP genes remains unclear.

RESULTS

In the present study, a comparative genomic study of STP genes in five representative Brassicaceae crops was conducted, and a total of 25, 25, 28, 36 and 49 STP genes were individually identified in Raphanus sativus (Rs), Brassica oleracea (Bo), B. rapa (Br), B. napus (Bn) and B. juncea (Bj), which were divided into four clades by phylogenetic analysis. The number of STP genes was no direct correlation with genome size and the total number of coding genes in Brassicaceae crops, and their physical and chemical properties showed no significant difference. Expression analysis showed that radish STP genes play vital roles not only in flower and seedpod development but also under heavy metal (cadmium, chromium and lead), NaCl and PEG-6000 stresses, Agrobacterium tumefaciens infection, and exogenous sugar treatment. RsSTP13.2 was significantly upregulated in the resistant radish cultivar by A. tumefaciens infection and induced by heavy metal, NaCl and PEG-6000 stress, indicating that it is involved in resistance to both biotic and abiotic stress in radish.

CONCLUSIONS

The present study provides insights into the evolutionary patterns of the STP gene family in Brassicaceae genomes and provides a theoretical basis for future functional analysis of STP genes in Brassicaceae crops.

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.

Pangenomes as a Resource to Accelerate Breeding of Under-Utilised Crop Species

Pangenomes are a rich resource to examine the genomic variation observed within a species or genera, supporting population genetics studies, with applications for the improvement of crop traits. Major crop species such as maize (Zea mays), rice (Oryza sativa), Brassica (Brassica spp.), and soybean (Glycine max) have had pangenomes constructed and released, and this has led to the discovery of valuable genes associated with disease resistance and yield components. However, pangenome data are not available for many less prominent crop species that are currently under-utilised. Despite many under-utilised species being important food sources in regional populations, the scarcity of genomic data for these species hinders their improvement. Here, we assess several under-utilised crops and review the pangenome approaches that could be used to build resources for their improvement. Many of these under-utilised crops are cultivated in arid or semi-arid environments, suggesting that novel genes related to drought tolerance may be identified and used for introgression into related major crop species. In addition, we discuss how previously collected data could be used to enrich pangenome functional analysis in genome-wide association studies (GWAS) based on studies in major crops. Considering the technological advances in genome sequencing, pangenome references for under-utilised species are becoming more obtainable, offering the opportunity to identify novel genes related to agro-morphological traits in these species.

SSR-Sequencing Reveals the Inter- and Intraspecific Genetic Variation and Phylogenetic Relationships among an Extensive Collection of Radish (Raphanus) Germplasm Resources

Simple Summary

Raphanus is an important genus of Brassicaceae and has undergone a lengthy evolutionary process. However, the inter- and intraspecific phylogenetic relationships and genetic diversity are not well understood. To elucidate these issues, we SSR-sequenced 939 wild, semi-wild and cultivated accessions, and discovered that Europe was the origin center of radishes with diverse European wild radishes, and Europe, South Asia and East Asia might be three independent domestication centers. There was considerable genetic differentiation within European cultivated radishes. European primitive cultivated radish exhibited gene flow with black radish/oil radish and rat-tail radish. Among Asian cultivated radishes, rat-tail radish was a sister to the clade of Chines big radish (including Japanese wild radish), suggesting that they may share the most recent common ancestry. Japanese wild radish had strong gene exchange with Japanese/Korea big radish, oil radish and rat-tail radish. American wild radish developed from natural hybridization between European wild radish and European small radish. All these demonstrated that European primitive cultivated type, American wild radish and Japanese wild radish might have played indispensable roles in radish evolution. Our study provides new perspectives into the origin, evolution and genetic diversity of Raphanus and facilitates the conservation and exploitation of radish germplasm resources.

Abstract

Raphanus has undergone a lengthy evolutionary process and has rich diversity. However, the inter- and intraspecific phylogenetic relationships and genetic diversity of this genus are not well understood. Through SSR-sequencing and multi-analysis of 939 wild, semi-wild and cultivated accessions, we discovered that the European wild radish (EWR) population is separated from cultivated radishes and has a higher genetic diversity. Frequent intraspecific genetic exchanges occurred in the whole cultivated radish (WCR) population; there was considerable genetic differentiation within the European cultivated radish (ECR) population, which could drive radish diversity formation. Among the ECR subpopulations, European primitive cultivated radishes (EPCRs) with higher genetic diversity are most closely related to the EWR population and exhibit a gene flow with rat-tail radishes (RTRs) and black radishes (BRs)/oil radishes (ORs). Among Asian cultivated radishes (ACRs), Chinese big radishes (CBRs) with a relatively high diversity are furthest from the EWR population, and most Japanese/Korean big radishes (JKBRs) are close to CBR accessions, except for a few old Japanese landraces that are closer to the EPCR. The CBR and JKBR accessions are independent of RTR accessions; however, phylogenetic analysis indicates that the RTR is sister to the clade of CBR (including JWR), which suggests that the RTR may share the most recent common ancestry with CBRs and JWRs. In addition, Japanese wild radishes (JWRs), (namely, R. sativus forma raphanistroides) are mainly scattered between CBRs and EPCRs in PCoA analysis. Moreover, JWRs have a strong gene exchange with the JKBR, OR and RTR subpopulations. American wild radishes (AWRs) are closely related to European wild and cultivated radishes, and have a gene flow with European small radishes (ESRs), suggesting that the AWR developed from natural hybridization between the EWR and the ESR. Overall, this demonstrates that Europe was the origin center of the radish, and that Europe, South Asia and East Asia appear to have been three independent domestication centers. The EPCR, AWR and JWR, as semi-wild populations, might have played indispensable transitional roles in radish evolution. Our study provides new perspectives into the origin, evolution and genetic diversity of Raphanus and facilitates the conservation and exploitation of radish germplasm resources.