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.

Genomic insight into the origin, domestication, dispersal, diversification and human selection of Tartary buckwheat.

BACKGROUND: Tartary buckwheat, Fagopyrum tataricum, is a pseudocereal crop with worldwide distribution and high nutritional value. However, the origin and domestication history of this crop remain to be elucidated. RESULTS: Here, by analyzing the population genomics of 567 accessions collected worldwide and reviewing historical documents, we find that Tartary buckwheat originated in the Himalayan region and then spread southwest possibly along with the migration of the Yi people, a minority in Southwestern China that has a long history of planting Tartary buckwheat. Along with the expansion of the Mongol Empire, Tartary buckwheat dispersed to Europe and ultimately to the rest of the world. The different natural growth environments resulted in adaptation, especially significant differences in salt tolerance between northern and southern Chinese Tartary buckwheat populations. By scanning for selective sweeps and using a genome-wide association study, we identify genes responsible for Tartary buckwheat domestication and differentiation, which we then experimentally validate. Comparative genomics and QTL analysis further shed light on the genetic foundation of the easily dehulled trait in a particular variety that was artificially selected by the Wa people, a minority group in Southwestern China known for cultivating Tartary buckwheat specifically for steaming as a staple food to prevent lysine deficiency. CONCLUSIONS: This study provides both comprehensive insights into the origin and domestication of, and a foundation for molecular breeding for, Tartary buckwheat.

Multi-omics identification of a key glycosyl hydrolase gene FtGH1 involved in rutin hydrolysis in Tartary buckwheat (Fagopyrum tataricum).

Rutin, a flavonoid rich in buckwheat, is important for human health and plant resistance to external stresses. The hydrolysis of rutin to quercetin underlies the bitter taste of Tartary buckwheat. In order to identify rutin hydrolysis genes, a 200 genotypes mini-core Tartary buckwheat germplasm resource was re-sequenced with 30-fold coverage depth. By combining the content of the intermediate metabolites of rutin metabolism with genome resequencing data, metabolite genome-wide association analyses (GWAS) eventually identified a glycosyl hydrolase gene FtGH1, which could hydrolyse rutin to quercetin. This function was validated both in Tartary buckwheat overexpression hairy roots and in vitro enzyme activity assays. Mutation of the two key active sites, which were determined by molecular docking and experimentally verified via overexpression in hairy roots and transient expression in tobacco leaves, exhibited abnormal subcellular localization, suggesting functional changes. Sequence analysis revealed that mutation of the FtGH1 promoter in accessions of two haplotypes might be necessary for enzymatic activity. Co-expression analysis and GWAS revealed that FtbHLH165 not only repressed FtGH1 expression, but also increased seed length. This work reveals a potential mechanism behind rutin metabolism, which should provide both theoretical support in the study of flavonoid metabolism and in the molecular breeding of Tartary buckwheat.

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.

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.

Transcriptome Analysis and HPLC Profiling of Flavonoid Biosynthesis in Citrus aurantium L. during Its Key Developmental Stages.

Citrus aurantium L. (sour orange) is a significant Chinese medicinal and fruit crop rich in flavonoids. However, the pathways and genes involved in flavonoid biosynthesis at the key developmental stages of Citrus aurantium L. are not fully understood. This study found that the total flavonoid concentration gradually decreased as the fruit developed. Additionally, it showed that neohesperidin was the main flavonoid in the early stages of sour orange fruit development. However, as the development stage progressed, naringin content increased rapidly and emerged as the main flavonoid component. From 27 cDNA libraries, RNA sequencing yielded 16.64 billion clean bases, including 8989 differentially expressed genes. We identified 74 flavonoid related unigenes mapped to the phenylalanine, tyrosine, and phenylpropanoid biosynthesis pathways. A total of 152 UDP-glucuronosyltransferase genes (UGTs) were identified from C. aurantium L. transcriptome database, in which 22 key flavonoid-correlated UGTs were divided into five main AtGT groups: E, G, I, L, M. We observed that the ethylene responsive factors (ERF) and myeloblastosis (MYB) family mainly regulated the key genes involved in flavonoid biosynthesis. Overall, our study generated extensive and detailed transcriptome data on the development of C. aurantium L. and characterized the flavonoid biosynthesis pattern during its fruit developmental stages. These results will benefit genetic modification or selection to increase the flavonoid content in sour oranges.

JA-induced FtBPM3 accumulation promotes FtERF-EAR3 degradation and rutin biosynthesis in Tartary buckwheat.

Buckwheat accumulates abundant flavonoids, which exhibit excellent health-promoting value. Flavonoids biosynthesis is mediated by a variety of phytohormones, among which jasmonates (JAs) induce numerous transcription factors, taking part in regulation of flavonoids biosynthesis genes. However, some transcriptional repressors appeared also induced by JAs. How these transcriptional repressors coordinately participate in JA signaling remains unclear. Here, we found that the disruption of the GCC-box in FtF3H promoter was associated with flavonoids accumulation in Tartary buckwheat. Further, our study illustrated that the nucleus-localized FtERF-EAR3 could inhibit FtF3H expression and flavonoids biosynthesis through binding the GCC-box in the promoter of FtF3H. The JA induced FtERF-EAR3 gene expression while facilitating FtERF-EAR3 protein degradation via the FtBPM3-dependent 26S proteasome pathway. Overall, these results illustrate a precise modulation mechanism of JA-responsive transcription suppressor participating in flavonoid biosynthesis, and will further help to improve the efficiency of flavonoids biosynthesis in Tartary buckwheat.

Physicochemical properties, content, composition and partial least squares models of A. trifoliata seeds oil.

Physicochemical properties, oil content, and fatty acids (FAs) composition are key for determining the value of oil crops. The aim of this study was to illustrate the potential of exploiting A. trifoliata as an edible oil crop, and establish a rapid measurement model for the A. trifoliata seeds oil (ASO) content and composition. In 130 A. trifoliata germplasms, the highest content of ASO was 51.27%, and unsaturated fatty acids (UFAs) mainly accounted for 74-78% of ASO. The partial least squares (PLS) model based on GC-MS and near-infrared spectroscopy was well-suited for the determination of ASO and UFA content; however, the PLS model for oleic acid (OA) and linoleic acid (LA) was not effective. The acid values and peroxide values for ASO also conformed to the Chinese food safety standards. Our findings will provide new insights and guidance for the use of A. trifoliata as oil crops..

Genetic Diversity of a Natural Population of Akebia trifoliata (Thunb.) Koidz and Extraction of a Core Collection Using Simple Sequence Repeat Markers.

Understand genetic diversity and genetic structure of germplasm is premise of germplasm conservation and utilization. And core collection can reduce the cost and difficulty of germplasm conservation. Akebia trifoliata (Thunb.) Koidz is an important medicinal, fruit and oil crop, particularly in China. In this study, 28 simple sequence repeat (SSR) markers were used to assess the genetic diversity and genetic structure of 955 A. trifoliata germplasms, determine their molecular identity and extract a core collection. The genetic diversity of the 955 germplasms was moderately polymorphic. The average number of alleles (Na), observed heterozygosity (H O ), expected heterozygosity (H E ), Shannon's information index (I∗), and polymorphic information content (PIC) were 3.71, 0.24, 0.46, 0.81, and 0.41, respectively. Four subpopulations were identified, indicating a weak genetic structure. A 955 germplasms could be completely distinguished by the characters of s28, s25, s74, s89, s68, s30, s13, s100, s72, s77, and s3. And each germplasm's molecular identity was made up of eleven characters. The core collection was composed of 164 germplasms (17.2% of 955 total germplasms in the population) and diversity parameters differed significantly from those of a random core collection. These results have implications for germplasm conservation. At the same time, based on the results, the 955 germplasm could be better used and managed.

Genome-wide development of insertion-deletion (InDel) markers for Cannabis and its uses in genetic structure analysis of Chinese germplasm and sex-linked marker identification.

BACKGROUND: Cannabis sativa L., a dioecious plant derived from China, demonstrates important medicinal properties and economic value worldwide. Cannabis properties have been usually harnessed depending on the sex of the plant. To analyse the genetic structure of Chinese Cannabis and identify sex-linked makers, genome-wide insertion-deletion (InDel) markers were designed and used. RESULTS: In this study, a genome-wide analysis of insertion-deletion (InDel) polymorphisms was performed based on the recent genome sequences. In total, 47,558 InDels were detected between the two varieties, and the length of InDels ranged from 4 bp to 87 bp. The most common InDels were tetranucleotides, followed by pentanucleotides. Chromosome 5 exhibited the highest number of InDels among the Cannabis chromosomes, while chromosome 10 exhibited the lowest number. Additionally, 31,802 non-redundant InDel markers were designed, and 84 primers evenly distributed in the Cannabis genome were chosen for polymorphism analysis. A total of 38 primers exhibited polymorphisms among three accessions, and of the polymorphism primers, 14 biallelic primers were further used to analyse the genetic structure. A total of 39 fragments were detected, and the PIC value ranged from 0.1209 to 0.6351. According to the InDel markers and the flowering time, the 115 Chinese germplasms were divided into two subgroups, mainly composed of cultivars obtained from the northernmost and southernmost regions, respectively. Additional two markers, "Cs-I1-10" and "Cs-I1-15", were found to amplify two bands (398 bp and 251 bp; 293 bp and 141 bp) in the male plants, while 389-bp or 293-bp bands were amplified in female plants. Using the two markers, the feminized and dioecious varieties could also be distinguished. CONCLUSION: Based on the findings obtained herein, we believe that this study will facilitate the genetic improvement and germplasm conservation of Cannabis in China, and the sex-linked InDel markers will provide accurate sex identification strategies for Cannabis breeding and production.

Comparative Analysis of Akebia trifoliata Fruit Softening at Different Flesh Ripening Stages Using Tandem Mass Tag Technology.

Owing to its medicinal and high nutritional values, Akebia trifoliata can be considered as a new type of medicinal and edible homologous resources, and it has begun to be widely cultivated in many areas of China. Over-softening of fruit would affect the sensorial quality, utilization rate, and consumer acceptance of the fruit postharvest. However, fruit softening has not been characterized and the molecular mechanism underlying A. trifoliata fruit softening during ripening remains unclear. A comparative proteomic analysis was performed on the fruit at three developmental stages using tandem mass tag technology. In total, 2,839 proteins and 302 differentially abundant proteins (DAPs) were identified. Bioinformatics analysis indicated that most DAPs were implicated in oxidoreductase activity, protein domain-specific binding and pyruvate metabolism. Moreover, 29 DAPs associated with cell wall metabolism, plant hormone, and stress and defense response pathways were validated using quantitative PCR. Notably, pectinesterase, pectate lyase, and ß-galactosidase, which are involved in cell wall degradation, as well as gibberellin regulated protein, cysteine protease, thaumatin-like protein and heat shock proteins which is involved in plant hormone, and stress and defense response, were significantly up-regulated in softening fruit compared with the levels in non-softening fruit. This indicated that they might play key roles in A. trifoliata fruit softening. Our findings will provide new insights into potential genes influencing fruit softening traits of A. trifoliata, which will help to develop strategies to improve fruit quality and reduce softening-related losses.

Integrative transcriptome and proteome analyses provide new insights into different stages of Akebia trifoliata fruit cracking during ripening.

BACKGROUND: Akebia trifoliata (Thunb.) Koidz may have applications as a new potential source of biofuels owing to its high seed count, seed oil content, and in-field yields. However, the pericarp of A. trifoliata cracks longitudinally during fruit ripening, which increases the incidence of pests and diseases and can lead to fruit decay and deterioration, resulting in significant losses in yield. Few studies have evaluated the mechanisms underlying A. trifoliata fruit cracking. RESULTS: In this study, by observing the cell wall structure of the pericarp, we found that the cell wall became thinner and looser and showed substantial breakdown in the pericarp of cracking fruit compared with that in non-cracking fruit. Moreover, integrative analyses of transcriptome and proteome profiles at different stages of fruit ripening demonstrated changes in the expression of various genes and proteins after cracking. Furthermore, the mRNA levels of 20 differentially expressed genes were analyzed, and parallel reaction monitoring analysis of 20 differentially expressed proteins involved in cell wall metabolism was conducted. Among the molecular targets, pectate lyases and pectinesterase, which are involved in pentose and glucuronate interconversion, and ß-galactosidase 2, which is involved in galactose metabolism, were significantly upregulated in cracking fruits than in non-cracking fruits. This suggested that they might play crucial roles in A. trifoliata fruit cracking. CONCLUSIONS: Our findings provided new insights into potential genes influencing the fruit cracking trait in A. trifoliata and established a basis for further research on the breeding of cracking-resistant varieties to increase seed yields for biorefineries.

Development of SSR markers via de novo transcriptome assembly in Akebia trifoliata (Thunb.) Koidz.

Owing to its high nutritive, economic, and medicinal values, Akebia trifoliata has received increased attention, making worthy of being used as a new fruit crop for further domestication and commercialization in China. However, molecular research of A. trifoliata has lagged as investigations of its genomic resources and molecular markers are rare. In this study, a cDNA library of A. trifoliata leaves was sequenced using the Illumina NovaSeq. 6000 sequencing system. In total, 101 417 transcripts, 63 757 unigenes, and 9494 simple sequence repeats were assembled and identified from the transcriptome datasets. The majority of the SSRs were di- and trinucleotide repeats. Length and number of SSR motifs ranged from 15 to 66, and 5 to 48 bp, respectively. Of which, the A/T mononucleotide motif and AG/TC and CT/GA dinucleotide motifs were the most abundant. Furthermore, 100 SSR primers were randomly selected to validate amplification and polymorphism, and 88 A. trifoliata accessions were definitively distinguished by 49 primers. With the Qinling mountains and Huaihe River line as the boundaries, the northern and southern accessions were clustered into different groups, but no clear geographical patterns (city or origin) were observed in the southern accessions. These newly identified molecular markers may provide a foundation for the genetic identification and diversity analysis and marker-assisted selection breeding in species of Akebia.

Comparative transcriptomics reveals the selection patterns of domesticated ramie.

Although domestication has dramatically altered the phenotype, physiology, and life history of ramie (Boehmeria nivea) plants, few studies have investigated the effects of domestication on the structure and expression pattern of genes in this fiber crop. To investigate the selective pattern and genetic relationships among a cultivated variety of ramie (BNZ: B. nivea, ZZ1) and four wild species, BNT (B. nivea var. tenacissima), BNN (B. nivea var. nipononivea), BNW (B. nivea var. nivea), and BAN (B. nivea var. viridula), in the section Tilocnide, we performed an RNA sequencing analysis of these ramie species. The de novo assembly of the "all-ramie" transcriptome yielded 119,114 unigenes with an average length of 633 bp, and a total of 7,084 orthologous gene pairs were identified. The phylogenetic tree showed that the cultivar BNZ clustered with BAN in one group, BNW was closely related to BNT, and BNN formed a separate group. Introgression analysis indicated that gene flow occurred from BNZ to BNN and BAN, and between BAN and BNN. Among these orthologs, 2,425 and 269 genes underwent significant purifying and positive selection, respectively. For these positively selected genes, oxidation-reduction process (GO:0055114) and stress response pathways (GO:0006950) were enriched, indicating that modulation of the cellular redox status was important during both ramie fiber evolution and improvement. Two genes related to the suppression of flowering and one gene annotated as a flowering-promoting factor were subjected to positive selection, probably caused by human manipulation. Additionally, five genes were homologs of those involved in abiotic stress tolerance and disease resistance, with higher expression levels in the cultivar BNZ than in the wild species. Collectively, the results of this study indicated that domestication has resulted in the upregulation of many genes involved in the abiotic and biotic stress responses, fiber yield, and plant growth of ramie.