Rewiring of the seed metabolome during Tartary buckwheat domestication

Application of GWAS and mGWAS in Livestock and Poultry Breeding

In recent years, genome-wide association studies (GWAS) and metabolome genome-wide association studies (mGWAS) have emerged as crucial methods for investigating complex traits in animals and plants. These have played pivotal roles in research on livestock and poultry breeding, facilitating a deeper understanding of genetic diversity, the relationship between genes, and genetic bases in livestock and poultry. This article provides a review of the applications of GWAS and mGWAS in animal genetic breeding, aiming to offer reference and inspiration for relevant researchers, promote innovation in animal genetic improvement and breeding methods, and contribute to the sustainable development of animal husbandry.

Genetic analysis of yield components in buckwheat using high-throughput sequencing analysis and wild resource populations

Fagopyrum tataricum, an important medicinal and edible crop, possesses significant agricultural and economic value. However, the development of buckwheat varieties and yields has been hindered by the delayed breeding progress despite the abundant material resources in China. Current research indicates that quantitative trait loci (QTLs) play a crucial role in controlling plant seed type and yield. To address these limitations, this study constructed recombinant inbred lines (RILs) utilizing both cultivated species and wild buckwheat as raw materials. In total, 84,521 Single Nucleotide Polymorphism (SNP) markers were identified through Genotyping-by-Sequencing (GBS) technology, and high-resolution and high-density SNP genetic maps were developed, which had significant value for QTL mapping, gene cloning and comparative mapping of buckwheat. In this study, we successfully identified 5 QTLs related to thousand grain weight (TGW), 9 for grain length (GL), and 1 for grain width (GW) by combining seed type and TGW data from 202 RIL populations in four different environments, within which one co-located QTL for TGW were discovered on the first chromosome. Transcriptome analysis during different grain development stages revealed 59 significant expression differences between the two materials, which can serve as candidate genes for further investigation into the regulation of grain weight and yield enhancement. The mapped major loci controlling TGW, GL and GW will be valuable for gene cloning and reveal the mechanism underlying grain development and marker-assisted selection in Tartary buckwheat.

Cratoxylum formosum ssp. pruniflorum induces gastric cancer cell apoptosis and pyroptosis through the elevation of ROS and cell cycle arrest

Cratoxylum formosum ssp. pruniflorum (CF), a traditional medicinal plant in Southern China, is widely recognized as a popular medicinal and tea plant traditionally utilized by diverse linguistic groups in the region for the treatment of gastrointestinal ailments. The objective of this study was to explore the active components and mechanisms of CF against gastric cancer (GC). The chemical ingredients of CF were obtained by using UPLC-MS/MS-based metabolomics. MGC-803 and HGC-27 cells were employed to investigate the direct anti-GC effect. The potential targets and signaling pathway of CF were identified through network pharmacology and proteomics, followed by subsequent experimental validation. Through UPLC-MS/MS metabolomics analysis, a total of 197 chemical ingredients were identified in CF leaves. Network pharmacology and proteomics techniques revealed 25 potential targets for GC, with a protein-protein interaction (PPI) network highlighting 12 cores targets, including CTNNB1, CDK2, et al. Furthermore, seven key CF ingredients - vismione B, feruloylcholine, α-amyrin, vanillic acid, galangin, cinnamic acid, and caffeic acid - were found to mediate anti-GC effects through pathways such as reactive oxygen species (ROS) and cell cycle signaling pathway. In vitro experiments demonstrated that CF significantly inhibited the proliferation and migration of GC cells, increased intracellular reactive oxygen species (ROS), malondialdehyde (MDA) and lactate dehydrogenase (LDH) levels, arrested the cell cycle at the S-phase, induced apoptosis and pyroptosis, and upregulated expression of apoptosis proteins (Bax, Bax/Bcl-2, cleaved-Caspase-3/Caspase-3), and pyroptosis proteins (GSDMD-N/GSDMD and GSDME-N/GSDME), while downregulating expression of cell cycle proteins (CDK2 and cyclin A1) as well as necroptosis proteins (RIP1 and MLKL). Collectively, these findings reveal CF's therapeutic potential against GC by the augmentation of ROS production, cell cycle arrest, promotion of apoptosis, and pyroptosis, offering valuable evidence for the development and utilization of CF in clinical settings.

Population comparative genomics discovers gene gain and loss during grapevine domestication

Plant domestication are evolutionary experiments conducted by early farmers since thousands years ago, during which the crop wild progenitors are artificially selected for desired agronomic traits along with dramatic genomic variation in the course of moderate to severe bottlenecks. However, previous investigations are mainly focused on small-effect variants, while changes in gene contents are rarely investigated due to the lack of population-level assemblies for both the crop and its wild relatives. Here, we applied comparative genomic analyses to discover gene gain and loss during grapevine domestication using long-read assemblies of representative population samples for both domesticated grapevines (V. vinifera ssp. vinifera) and their wild progenitors (V. vinifera ssp. sylvestris). Only ∼7% of gene families were shared by 16 Vitis genomes while ∼8% of gene families were specific to each accession, suggesting dramatic variations of gene contents in grapevine genomes. Compared to wild progenitors, the domesticated accessions exhibited an increased presence of genes associated with asexual reproduction, while the wild progenitors showcased a higher abundance of genes related to pollination, revealing the transition from sexual reproduction to clonal propagation during domestication processes. Moreover, the domesticated accessions harbored fewer disease-resistance genes than wild progenitors. The SVs occurred frequently in aroma and disease-resistance related genes between domesticated grapevines and wild progenitors, indicating the rapid diversification of these genes during domestication. Our study provides insights and resources for biological studies and breeding programs in grapevine.

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.

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.

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.

Agro-Morphological and Molecular Characterization Reveal Deep Insights in Promising Genetic Diversity and Marker-Trait Associations in Fagopyrum esculentum and Fagopyrum tataricum

Characterisation of genetic diversity is critical to adequately exploit the potential of germplasm collections and identify important traits for breeding programs and sustainable crop improvement. Here, we characterised the phenotypic and genetic diversity of a global collection of the two cultivated buckwheat species Fagopyrum esculentum and Fagopyrum tataricum (190 and 51 accessions, respectively) using 37 agro-morphological traits and 24 SSR markers. A wide range of variation was observed in both species for most of the traits analysed. The two species differed significantly in most traits, with traits related to seeds and flowering contributing most to differentiation. The accessions of each species were divided into three major phenoclusters with no clear geographic clustering. At the molecular level, the polymorphic SSR markers were highly informative, with an average polymorphic information content (PIC) of over 0.65 in both species. Genetic diversity, as determined by Nei's expected heterozygosity (He), was high (He = 0.77 and He = 0.66, respectively) and differed significantly between species (p = 0.03) but was homogeneously distributed between regions, confirming the lack of genetic structure as determined by clustering approaches. The weak genetic structure revealed by the phenotypic and SSR data and the low fixation indices in both species suggested frequent seed exchange and extensive cultivation and selection. In addition, 93 and 140 significant (p < 0.05) marker-trait associations (MTAs) were identified in both species using a general linear model and a mixed linear model, most of which explained >20% of the phenotypic variation in associated traits. Core collections of 23 and 13 phenotypically and genetically diverse accessions, respectively, were developed for F. esculentum and F. tataricum. Overall, the data analysed provided deep insights into the agro-morphological and genetic diversity and genetic relationships among F. esculentum and F. tataricum accessions and pointed to future directions for genomics-based breeding programs and germplasm management.

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.

Metabolic-GWAS provides insights into genetic architecture of seed metabolome in buckwheat

BACKGROUND

Buckwheat (Fagopyrum spp.), belonging to the Polygonaceae family, is an ancient pseudo-cereal with high nutritional and nutraceutical properties. Buckwheat proteins are gluten-free and show balanced amino acid and micronutrient profiles, with higher content of health-promoting bioactive flavonoids that make it a golden crop of the future. Plant metabolome is increasingly gaining importance as a crucial component to understand the connection between plant physiology and environment and as a potential link between the genome and phenome. However, the genetic architecture governing the metabolome and thus, the phenome is not well understood. Here, we aim to obtain a deeper insight into the genetic architecture of seed metabolome in buckwheat by integrating high throughput metabolomics and genotyping-by-sequencing applying an array of bioinformatics tools for data analysis.

RESULTS

High throughput metabolomic analysis identified 24 metabolites in seed endosperm of 130 diverse buckwheat genotypes. The genotyping-by-sequencing (GBS) of these genotypes revealed 3,728,028 SNPs. The Genome Association and Prediction Integrated Tool (GAPIT) assisted in the identification of 27 SNPs/QTLs linked to 18 metabolites. Candidate genes were identified near 100 Kb of QTLs, providing insights into several metabolic and biosynthetic pathways.

CONCLUSIONS

We established the metabolome inventory of 130 germplasm lines of buckwheat, identified QTLs through marker trait association and positions of potential candidate genes. This will pave the way for future dissection of complex economic traits in buckwheat.

A metabolome genome-wide association study implicates histidine N-pi-methyltransferase as a key enzyme in N-methylhistidine biosynthesis in Arabidopsis thaliana

A genome-wide association study (GWAS), which uses information on single nucleotide polymorphisms (SNPs) from many accessions, has become a powerful approach to gene identification. A metabolome GWAS (mGWAS), which relies on phenotypic information based on metabolite accumulation, can identify genes that contribute to primary and secondary metabolite contents. In this study, we carried out a mGWAS using seed metabolomic data from Arabidopsis thaliana accessions obtained by liquid chromatography-mass spectrometry to identify SNPs highly associated with the contents of metabolites such as glucosinolates. These SNPs were present in genes known to be involved in glucosinolate biosynthesis, thus confirming the effectiveness of our analysis. We subsequently focused on SNPs detected in an unknown methyltransferase gene associated with N-methylhistidine content. Knockout and overexpression of A. thaliana lines of this gene had significantly decreased and increased N-methylhistidine contents, respectively. We confirmed that the overexpressing line exclusively accumulated histidine methylated at the pi position, not at the tau position. Our findings suggest that the identified methyltransferase gene encodes a key enzyme for N-methylhistidine biosynthesis in A. thaliana.

The SlDOG1 Affect Biosynthesis of Steroidal Glycoalkaloids by Regulating GAME Expression in Tomato

Steroidal alkaloids (SAs) and steroidal glycoalkaloids (SGAs) are common constituents of plant species belonging to the Solanaceae family. However, the molecular mechanism regulating the formation of SAs and SGAs remains unknown. Here, genome-wide association mapping was used to elucidate SA and SGA regulation in tomatoes: a SlGAME5-like glycosyltransferase (Solyc10g085240) and the transcription factor SlDOG1 (Solyc10g085210) were significantly associated with steroidal alkaloid composition. In this study, it was found that rSlGAME5-like can catalyze a variety of substrates for glycosidation and can catalyze SA and flavonol pathways to form O-glucoside and O-galactoside in vitro. The overexpression of SlGAME5-like promoted α-tomatine, hydroxytomatine, and flavonol glycoside accumulation in tomatoes. Furthermore, assessments of natural variation combined with functional analyses identified SlDOG1 as a major determinant of tomato SGA content, which also promoted SA and SGA accumulation via the regulation of GAME gene expression. This study provides new insights into the regulatory mechanisms underlying SGA production in tomatoes.