Comparison of buckwheat genomes reveals the genetic basis of metabolomic divergence and ecotype differentiation

Evolutionary analysis of MADS-box genes in buckwheat species and functional study of FdMADS28 in flavonoid metabolism

The MADS-box gene family is a transcription factor family that is widely expressed in plants. It controls secondary metabolic processes in plants and encourages the development of tissues like roots and flowers. However, the phylogenetic analysis and evolutionary model of MADS-box genes in Fagopyrum species has not been reported yet. This study identified the MADS-box genes of three buckwheat species at the whole genome level, and conducted systematic evolution and physicochemical analysis. The results showed that these genes can be divided into four subfamilies, with fragment duplication being the main way for the gene family expansion. During the domestication process from golden buckwheat to tartary buckwheat and the common buckwheat, the Ka/Ks ratio indicated that most members of the family experienced strong purification selection pressure, and with individual gene pairs experiencing positive selection. In addition, we combined the expression profile data of the MADS genes, mGWAS data, and WGCNA data to mine genes FdMADS28/48/50 that may be related to flavonoid metabolism. The results also showed that overexpression of FdMADS28 could increase rutin content by decreasing Kaempferol pathway content in hairy roots, and increase the resistance and growth of hairy roots to PEG and NaCl. This study systematically analyzed the evolutionary relationship of MADS-box genes in the buckwheat species, and elaborated on the expression patterns of MADS genes in different tissues under biotic and abiotic stresses, laying an important theoretical foundation for further elucidating their role in flavonoid metabolism.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

Exploration of unique starch physicochemical properties of novel buckwheat lines created by crossing Golden buckwheat and Tatary buckwheat

Buckwheat is considered as a healthy cereal food, and it is essential to cultivate new buckwheat lines with good starch physicochemical properties for both consumers and food producers. Six novel buckwheat (Duoku, Dk) were generated by crossing of Golden buckwheat and Tatary buckwheat, and their kernel appearance properties and starch physicochemical properties were analyzed together with one domestic line (Cimiqiao) and one wild line (Yeku). The results showed that Dk samples had better appearance properties than two control samples. The Dk samples showed lower amylose content, similar amylopectin molecular structure and chain length distributions, and larger starch granules compared with Cimiqiao. The digestion results showed that two Dk samples: Dk6 & Dk9 had high resistant starch content; while the other two Dk samples: Dk37 & Dk38 had a steady glucose releasing rate. The Dk samples also showed high gelatinization temperature, indicating they were good raw materials for producing glass noodle. This study proved that Dk buckwheat had unique starch physicochemical properties, and could be used as new food materials in the future.

Comparative and population genomics of buckwheat species reveal key determinants of flavor and fertility

Common buckwheat (Fagopyrum esculentum) is an ancient crop with a world-wide distribution. Due to its excellent nutritional quality and high economic and ecological value, common buckwheat is becoming increasingly important throughout the world. The availability of a high-quality reference genome sequence and population genomic data will accelerate the breeding of common buckwheat, but the high heterozygosity due to the outcrossing nature has greatly hindered the genome assembly. Here we report the assembly of a chromosome-scale high-quality reference genome of F. esculentum var. homotropicum, a homozygous self-pollinating variant of common buckwheat. Comparative genomics revealed that two cultivated buckwheat species, common buckwheat (F. esculentum) and Tartary buckwheat (F. tataricum), underwent metabolomic divergence and ecotype differentiation. The expansion of several gene families in common buckwheat, including FhFAR genes, is associated with its wider distribution than Tartary buckwheat. Copy number variation of genes involved in the metabolism of flavonoids is associated with the difference of rutin content between common and Tartary buckwheat. Furthermore, we present a comprehensive atlas of genomic variation based on whole-genome resequencing of 572 accessions of common buckwheat. Population and evolutionary genomics reveal genetic variation associated with environmental adaptability and floral development between Chinese and non-Chinese cultivated groups. Genome-wide association analyses of multi-year agronomic traits with the content of flavonoids revealed that Fh05G014970 is a potential major regulator of flowering period, a key agronomic trait controlling the yield of outcrossing crops, and that Fh06G015130 is a crucial gene underlying flavor-associated flavonoids. Intriguingly, we found that the gene translocation and sequence variation of FhS-ELF3 contribute to the homomorphic self-compatibility of common buckwheat. Collectively, our results elucidate the genetic basis of speciation, ecological adaptation, fertility, and unique flavor of common buckwheat, and provide new resources for future genomics-assisted breeding of this economically important crop.

A chromosome-scale Rhubarb (Rheum tanguticum) genome assembly provides insights into the evolution of anthraquinone biosynthesis

Rhubarb is the collective name for various perennial plants from the genus Rheum L. and the Polygonaceae family. They are one of the most ancient, commonly used, and important herbs in traditional Chinese medicine. Rhubarb is a major source of anthraquinones, but how they are synthesized remains largely unknown. Here, we generate a genome sequence assembly of one important medicinal rhubarb R. tanguticum at the chromosome level, with 2.76 Gb assembled into 11 chromosomes. The genome is shaped by two recent whole-genome duplication events and recent bursts of retrotransposons. Metabolic analyses show that the major anthraquinones are mainly synthesized in its roots. Transcriptomic analysis reveals a co-expression module with a high correlation to anthraquinone biosynthesis that includes key chalcone synthase genes. One CHS, four CYP450 and two BGL genes involved in secondary metabolism show significantly upregulated expression levels in roots compared with other tissues and clustered in the co-expression module, which implies that they may also act as candidate genes for anthraquinone biosynthesis. This study provides valuable insights into the genetic bases of anthraquinone biosynthesis that will facilitate improved breeding practices and agronomic properties for rhubarb in the future.

Multiomics analysis reveals the molecular mechanisms underlying virulence in Rhizoctonia and jasmonic acid-mediated resistance in Tartary buckwheat (Fagopyrum tataricum)

Rhizoctonia solani is a devastating soil-borne pathogen that seriously threatens the cultivation of economically important crops. Multiple strains with a very broad host range have been identified, but only 1 (AG1-IA, which causes rice sheath blight disease) has been examined in detail. Here, we analyzed AG4-HGI 3 originally isolated from Tartary buckwheat (Fagopyrum tataricum), but with a host range comparable to AG1-IA. Genome comparison reveals abundant pathogenicity genes in this strain. We used multiomic approaches to improve the efficiency of screening for disease resistance genes. Transcriptomes of the plant-fungi interaction identified differentially expressed genes associated with virulence in Rhizoctonia and resistance in Tartary buckwheat. Integration with jasmonate-mediated transcriptome and metabolome changes revealed a negative regulator of jasmonate signaling, cytochrome P450 (FtCYP94C1), as increasing disease resistance probably via accumulation of resistance-related flavonoids. The integration of resistance data for 320 Tartary buckwheat accessions identified a gene homolog to aspartic proteinase (FtASP), with peak expression following R. solani inoculation. FtASP exhibits no proteinase activity but functions as an antibacterial peptide that slows fungal growth. This work reveals a potential mechanism behind pathogen virulence and host resistance, which should accelerate the molecular breeding of resistant varieties in economically essential crops.

Succession of endophytic fungi and rhizosphere soil fungi and their correlation with secondary metabolites in Fagopyrum dibotrys

Golden buckwheat (Fagopyrum dibotrys, also known as F. acutatum) is a traditional edible herbal medicinal plant with a large number of secondary metabolites and is considered to be a source of therapeutic compounds. Different ecological environments have a significant impact on their compound content and medicinal effects. However, little is known about the interactions between soil physicochemical properties, the rhizosphere, endophytic fungal communities, and secondary metabolites in F. dibotrys. In this study, the rhizosphere soil and endophytic fungal communities of F. dibotrys in five different ecological regions in China were identified based on high-throughput sequencing methods. The correlations between soil physicochemical properties, active components (total saponins, total flavonoids, proanthocyanidin, and epicatechin), and endophytic and rhizosphere soil fungi of F. dibotrys were analyzed. The results showed that soil pH, soil N, OM, and P were significantly correlated with the active components of F. dibotrys. Among them, epicatechin, proanthocyanidin, and total saponins were significantly positively correlated with soil pH, while proanthocyanidin content was significantly positively correlated with STN, SAN, and OM in soil, and total flavone content was significantly positively correlated with P in soil. In soil microbes, Mortierella, Trechispora, Exophiala, Ascomycota_unclassified, Auricularia, Plectosphaerella, Mycena, Fungi_unclassified, Agaricomycetes_unclassified, Coprinellus, and Pseudaleuria were significantly related to key secondary metabolites of F. dibotrys. Diaporthe and Meripilaceae_unclassified were significantly related to key secondary metabolites in the rhizome. This study presents a new opportunity to deeply understand soil-plant-fungal symbioses and secondary metabolites in F. dibotrys, as well as provides a scientific basis for using biological fertilization strategies to improve the quality of F. dibotrys.

Recent Advances in Assembly of Complex Plant Genomes

Over the past 20 years, tremendous advances in sequencing technologies and computational algorithms have spurred plant genomic research into a thriving era with hundreds of genomes decoded already, ranging from those of nonvascular plants to those of flowering plants. However, complex plant genome assembly is still challenging and remains difficult to fully resolve with conventional sequencing and assembly methods due to high heterozygosity, highly repetitive sequences, or high ploidy characteristics of complex genomes. Herein, we summarize the challenges of and advances in complex plant genome assembly, including feasible experimental strategies, upgrades to sequencing technology, existing assembly methods, and different phasing algorithms. Moreover, we list actual cases of complex genome projects for readers to refer to and draw upon to solve future problems related to complex genomes. Finally, we expect that the accurate, gapless, telomere-to-telomere, and fully phased assembly of complex plant genomes could soon become routine.

Genomic analyses provide insights into the polyploidization-driven herbicide adaptation in Leptochloa weeds

Polyploidy confers a selective advantage under stress conditions; however, whether polyploidization mediates enhanced herbicide adaptation remains largely unknown. Tetraploid Leptochloa chinensis is a notorious weed in the rice ecosystem, causing severe yield loss in rice. In China, L. chinensis has only one sister species, the diploid L. panicea, whose damage is rarely reported. To gain insights into the effects of polyploidization on herbicide adaptation, we first assembled a high-quality genome of L. panicea and identified genome structure variations with L. chinensis. Moreover, we identified herbicide-resistance genes specifically expanded in L. chinensis, which may confer a greater herbicide adaptability in L. chinensis. Analysis of gene retention and loss showed that five herbicide target-site genes and several herbicide nontarget-site resistance gene families were retained during polyploidization. Notably, we identified three pairs of polyploidization-retained genes including LcABCC8, LcCYP76C1 and LcCYP76C4 that may enhance herbicide resistance. More importantly, we found that both copies of LcCYP76C4 were under herbicide selection during the spread of L. chinensis in China. Furthermore, we identified another gene potentially involved in herbicide resistance, LcCYP709B2, which is also retained during polyploidization and under selection. This study provides insights into the genomic basis of the enhanced herbicide adaptability of Leptochloa weeds during polyploidization and provides guidance for the precise and efficient control of polyploidy weeds.

The pan-plastome of tartary buckwheat (fagopyrum tataricum): key insights into genetic diversity and the history of lineage divergence

BACKGROUND

Tartary buckwheat (Fagopyrum tataricum) is an important food and medicine crop plant, which has been cultivated for 4000 years. A nuclear genome has been generated for this species, while an intraspecific pan-plastome has yet to be produced. As such a detailed understanding of the maternal genealogy of Tartary buckwheat has not been thoroughly investigated.

RESULTS

In this study, we de novo assembled 513 complete plastomes of Fagopyrum and compared with 8 complete plastomes of Fagopyrum downloaded from the NCBI database to construct a pan-plastome for F. tartaricum and resolve genomic variation. The complete plastomes of the 513 newly assembled Fagopyrum plastome sizes ranged from 159,253 bp to 159,576 bp with total GC contents ranged from 37.76 to 37.97%. These plastomes all maintained the typical quadripartite structure, consisting of a pair of inverted repeat regions (IRA and IRB) separated by a large single copy region (LSC) and a small single copy region (SSC). Although the structure and gene content of the Fagopyrum plastomes are conserved, numerous nucleotide variations were detected from which population structure could be resolved. The nucleotide variants were most abundant in the non-coding regions of the genome and of those the intergenic regions had the most. Mutational hotspots were primarily found in the LSC regions. The complete 521 Fagopyrum plastomes were divided into five genetic clusters, among which 509 Tartary buckwheat plastomes were divided into three genetic clusters (Ft-I/Ft-II/Ft-III). The genetic diversity in the Tartary buckwheat genetic clusters was the greatest in Ft-III, and the genetic distance between Ft-I and Ft-II was the largest. Based on the results of population structure and genetic diversity analysis, Ft-III was further subdivided into three subgroups Ft-IIIa, Ft-IIIb, and Ft-IIIc. Divergence time estimation indicated that the genera Fagopyrum and Rheum (rhubarb) shared a common ancestor about 48 million years ago (mya) and that intraspecies divergence in Tartary buckwheat began around 0.42 mya.

CONCLUSIONS

The resolution of pan-plastome diversity in Tartary buckwheat provides an important resource for future projects such as marker-assisted breeding and germplasm preservation.

Six Underutilized Grain Crops for Food and Nutrition in China

Underutilized grain crops are an essential part of the food system that supports humankind. A number of these crops can be found in China, such as barley, buckwheat, broomcorn millet, foxtail millet, oat, and sorghum, which have characteristics such as containing more nutritional elements, being resistant to biotic and abiotic stresses, and having strong adaptability to poor environments. The diversity of these crops provides options for farmers’ livelihoods and healthy food for the population. Although some mentioned crops such as barley, oat, and sorghum are not underutilized crops globally, they could be considered underutilized in China as they were more important in the past and could be revitalized for food and nutrition in the future. This paper reviews current progress in research and development in the areas of germplasm resource conservation, variety improvement, cultivation technologies, processing, and the nutrition and benefits of six underutilized grain crops in China. It is concluded that underutilized grain crops could play a critical role in food and nutritional security in China.

Identification of Medicinal Compounds of Fagopyri Dibotryis Rhizome from Different Origins and Its Varieties Using UPLC-MS/MS-Based Metabolomics

Fagopyrum dibotrys, being native to southwest China, is widely distributed in Yunnan, Guizhou Provinces and Chongqing City. However, the quality of medicinal materials growing in different origins varies greatly, and cannot meet the market demand for high-quality F. dibotrys. In this study, 648 metabolites were identified, and phenolic compounds of F. dibotrys from different origins were clearly separated by principal component analysis (PCA). Our results suggested that the medicinal differences of F. dibotrys from different origins can be elucidated via the variations in the abundance of the phenolic and flavonoid compounds. We found that the epicatechin, total flavonoids and total tannin content in Yunnan Qujing (YQ) and Yunnan Kunming (YK) were higher than those in Chongqing Shizhu (CS), Chongqing Fuling (CF) and Guizhou Bijie (GB), suggesting that Yunnan Province can be considered as one of the areas that produce high-quality medicinal materials. Additionally, 1,6-di-O-galloyl-β-D-glucose, 2,3-di-O-galloyl-D-glucose and gallic acid could be used as ideal marker compounds for the quality control of F. dibotrys from different origins caused by metabolites, and the F. dibotrys planted in Yunnan Province is well worth exploiting.