Genome sequencing reveals the genetic architecture of heterostyly and domestication history of common buckwheat.

Common buckwheat, Fagopyrum esculentum, is an orphan crop domesticated in southwest China that exhibits heterostylous self-incompatibility. Here we present chromosome-scale assemblies of a self-compatible F. esculentum accession and a self-compatible wild relative, Fagopyrum homotropicum, together with the resequencing of 104 wild and cultivated F. esculentum accessions. Using these genomic data, we report the roles of transposable elements and whole-genome duplications in the evolution of Fagopyrum. In addition, we show that (1) the breakdown of heterostyly occurs through the disruption of a hemizygous gene jointly regulating the style length and female compatibility and (2) southeast Tibet was involved in common buckwheat domestication. Moreover, we obtained mutants conferring the waxy phenotype for the first time in buckwheat. These findings demonstrate the utility of our F. esculentum assembly as a reference genome and promise to accelerate buckwheat research and breeding.

Genetic basis of maturity time is independent from that of flowering time and contributes to ecotype differentiation in common buckwheat (Fagopyrum esculentum Moench).

BACKGROUND: Common buckwheat is considered a quantitative short-day plant and is classified into the autumn (highly photoperiod sensitive), summer (weakly photoperiod sensitive), and intermediate ecotype. Understanding ecotype differentiation is essential for adaptive expansion and maximizing yield. The genetic analysis for ecotype has focused on photoperiod-dependent flowering time, whereas post-flowering traits such as seed set and maturity time might also regulate ecotype differentiation. RESULTS: A field experiment revealed that ecotype differentiation is mainly defined by the timing of seed set and maturation, whereas flowering time is less relevant. Thus, we focused on maturity time as a trait that defines the ecotype. To detect QTLs for maturity time, we developed two F2 populations derived from early × late-maturing accessions and intermediate × late-maturing accessions. Using genotyping by random amplicon sequencing-direct analysis, we generated a high-density linkage map. QTL analysis detected two major QTLs for maturity time, one in each F2 population. We also detected QTLs for flowering time at loci different from maturity time QTLs, which suggests that different genetic mechanisms regulate flowering and maturity. Association analysis showed that both QTLs for maturity time were significantly associated with variations in the trait across years. CONCLUSIONS: Maturity time appeared to be more suitable for explaining ecotype differentiation than flowering time, and different genetic mechanisms would regulate the timing of flowering and maturation. The QTLs and QTL-linked markers for maturity time detected here may be useful to extend the cultivation area and to fine-tune the growth period to maximize yield in buckwheat.

Targeted amplicon sequencing + next-generation sequencing-based bulked segregant analysis identified genetic loci associated with preharvest sprouting tolerance in common buckwheat (Fagopyrum esculentum).

BACKGROUND: Common buckwheat (2n = 2x = 16) is an outcrossing pseudocereal whose seeds contain abundant nutrients and potential antioxidants. As these beneficial compounds are damaged by preharvest sprouting (PHS) and PHS is likely to increase with global warming, it is important to find efficient ways to develop new PHS-tolerant lines. However, genetic loci and selection markers associated with PHS in buckwheat have not been reported. RESULTS: By next-generation sequencing (NGS) of whole-genome of parental lines, we developed a genome-wide set of 300 markers. By NGS- based bulked segregant analysis (NGS-BSA), we developed 100 markers linked to PHS tolerance. To confirm the effectiveness of marker development from NGS-BSA data, we developed 100 markers linked to the self-compatibility (SC) trait from previous NGS-BSA data. Using these markers, we developed genetic maps with AmpliSeq technology, which can quickly detect polymorphisms by amplicon-based multiplex targeted NGS, and performed quantitative trait locus (QTL) analysis for PHS tolerance in combination with NGS-BSA. QTL analysis detected two major and two minor QTLs for PHS tolerance in a segregating population developed from a cross between the PHS-tolerant 'Kyukei 29' and the self-compatible susceptible 'Kyukei SC7'. We found different major and minor QTLs in other segregating populations developed from the PHS-tolerant lines 'Kyukei 28' and 'NARO-FE-1'. Candidate markers linked to PHS developed by NGS-BSA were located near these QTL regions. We also investigated the effectiveness of markers linked to these QTLs for selection of PHS-tolerant lines among other segregating populations. CONCLUSIONS: We efficiently developed genetic maps using a method combined with AmpliSeq technology and NGS-BSA, and detected QTLs associated with preharvest sprouting tolerance in common buckwheat. This is the first report to identify QTLs for PHS tolerance in buckwheat. Our marker development system will accelerate genetic research and breeding in common buckwheat.

Genetic and genomic research for the development of an efficient breeding system in heterostylous self-incompatible common buckwheat (Fagopyrum esculentum).

Common buckwheat (Fagopyrum esculentum Moench; 2n = 2x = 16) is an annual crop that is cultivated widely around the world and contains an abundance of nutrients and bioactive compounds. However, the yield of buckwheat is low compared to that of other major crops, and it contains proteins that cause allergic reactions in some people. Much research has aimed to improve or eliminate these undesirable traits, and some major advances have recently been made. Here, we review recent advances in buckwheat breeding materials, tools, and methods, including the development of self-compatible lines, genetic maps, a buckwheat genome database, and an efficient breeding strategy. We also describe emerging breeding methods for high-value lines.

Identification of a gene encoding polygalacturonase expressed specifically in short styles in distylous common buckwheat (Fagopyrum esculentum).

Common buckwheat (Fagopyrum esculentum) is a heteromorphic self-incompatible (SI) species with two types of floral architecture: thrum (short style) and pin (long style). The floral morphology and intra-morph incompatibility are controlled by a single genetic locus, S. However, the molecular mechanisms underlying the heteromorphic self-incompatibility of common buckwheat remain unclear. To identify these mechanisms, we performed proteomic, quantitative reverse-transcription PCR, and linkage analyses. Comparison of protein profiles between the long and short styles revealed a protein unique to the short style. Amino-acid sequencing revealed that it was a truncated form of polygalacturonase (PG); we designated the gene encoding this protein FePG1. Phylogenetic analysis classified FePG1 into the same clade as PGs that function in pollen development and floral morphology. FePG1 expression was significantly higher in short styles than in long styles. It was expressed in flowers of a short-homostyle line but not in flowers of a long-homostyle line. Linkage analysis indicated that FePG1 was not linked to the S locus; it could be a factor downstream of this locus. Our finding of a gene putatively working under the regulation of the S locus provides useful information for elucidation of the mechanism of heteromorphic self-incompatibility.

Identification of a gene encoding glutathione S-transferase that is related to anthocyanin accumulation in buckwheat (Fagopyrum esculentum).

Anthocyanins are a group of flavonoids found in buckwheat (Fagopyrum esculentum) and many other plant species; however; little is known about their mechanisms of synthesis and regulation in buckwheat. We previously reported a spontaneous mutant buckwheat line that shows the green stem phenotype; this line does not accumulate anthocyanins but synthesizes flavonol and proanthocyanidin in the stem. Here, we used this line and lines developed by this line to search for genes related to anthocyanin accumulation in buckwheat. The lines with green stem showed flavonoid-3-O-glucosyltransferase activity against UDP-glucose, indicating that the flavonoid-3-O-glucosyltransferase gene was not controlling the green stem trait. We therefore searched the buckwheat genome database for a gene encoding glutathione S-transferase (GST), a flavonoid-binding protein that transports flavonoids to the vacuole, and identified a candidate gene, FeGST1. Expression analysis showed that FeGST1 was expressed in wild type buckwheat but not in the green stem lines. Linkage analysis with an F2 segregating population produced by crossing between the green stem line and a self-compatible line showed that FeGST1 segregated with stem color without any recombination. This indicates that the green stem trait could be caused by homozygous non-functional alleles of the FeGST1 locus.

Buckwheat R2R3 MYB transcription factor FeMYBF1 regulates flavonol biosynthesis.

Buckwheat (Fagopyrum esculentum) contains high amounts of flavonoids, especially flavonols (e.g., rutin), which are thought to be highly beneficial for human health. Little is known, however, about the regulation of flavonol synthesis in buckwheat. We identified a buckwheat gene encoding an R2R3 MYB transcription factor, and named this gene FeMYBF1. Analysis of the deduced amino acid sequence and phylogenetic analysis suggested that FeMYBF1 encodes an ortholog of the Arabidopsis flavonol regulators AtMYB11, AtMYB12 and AtMYB111. Expression of FeMYBF1 in a flavonol-deficient Arabidopsis triple mutant (myb11 myb12 myb111) restored flavonol synthesis. Constitutive expression of FeMYBF1 driven by the CaMV 35S promoter in Arabidopsis resulted in over-accumulation of flavonol glycosides and upregulation of the expression of AtFLS1. Transient expression assays showed that FeMYBF1 activated the promoter of the Arabidopsis gene encoding AtFLS1, and the promoters of buckwheat genes related to anthocyanin and proanthocyanidin synthesis such as dihydroflavonol 4-reductase (DFR) and leucoanthocyanidin dioxygenase (LDOX) in addition to genes encoding FLS. The results indicate that FeMYBF1 regulates flavonol synthesis and may have a role in synthesis of other flavonoid compounds, and also that buckwheat may have alternative pathway of flavonol synthesis through DFR and LDOX.

A new buckwheat dihydroflavonol 4-reductase (DFR), with a unique substrate binding structure, has altered substrate specificity.

BACKGROUND: Dihydroflavonol 4-reductase (DFR) is the key enzyme committed to anthocyanin and proanthocyanidin biosynthesis in the flavonoid biosynthetic pathway. DFR proteins can catalyse mainly the three substrates (dihydrokaempferol, dihydroquercetin, and dihydromyricetin), and show different substrate preferences. Although relationships between the substrate preference and amino acids in the region responsible for substrate specificity have been investigated in several plant species, the molecular basis of the substrate preference of DFR is not yet fully understood. RESULTS: By using degenerate primers in a PCR, we isolated two cDNA clones that encoded DFR in buckwheat (Fagopyrum esculentum). Based on sequence similarity, one cDNA clone (FeDFR1a) was identical to the FeDFR in DNA databases (DDBJ/Gen Bank/EMBL). The other cDNA clone, FeDFR2, had a similar sequence to FeDFR1a, but a different exon-intron structure. Linkage analysis in an F2 segregating population showed that the two loci were linked. Unlike common DFR proteins in other plant species, FeDFR2 contained a valine instead of the typical asparagine at the third position and an extra glycine between sites 6 and 7 in the region that determines substrate specificity, and showed less activity against dihydrokaempferol than did FeDFR1a with an asparagine at the third position. Our 3D model suggested that the third residue and its neighbouring residues contribute to substrate specificity. FeDFR1a was expressed in all organs that we investigated, whereas FeDFR2 was preferentially expressed in roots and seeds. CONCLUSIONS: We isolated two buckwheat cDNA clones of DFR genes. FeDFR2 has unique structural and functional features that differ from those of previously reported DFRs in other plants. The 3D model suggested that not only the amino acid at the third position but also its neighbouring residues that are involved in the formation of the substrate-binding pocket play important roles in determining substrate preferences. The unique characteristics of FeDFR2 would provide a useful tool for future studies on the substrate specificity and organ-specific expression of DFRs.

Isolation and characterization of genes encoding leucoanthocyanidin reductase (FeLAR) and anthocyanidin reductase (FeANR) in buckwheat (Fagopyrum esculentum).

Proanthocyanidins (PAs) are a major group of flavonoids synthesized via the phenylpropanoid biosynthesis pathway, however the pathway has not been fully characterized in buckwheat. Anthocyanidin reductase (ANR) and leucoanthocyanidin reductase (LAR) are involved in the last steps of PA biosynthesis. To isolate the genes for these enzymes from buckwheat we performed PCR using degenerate primers and obtained cDNAs of ANR and LAR, which we designated FeANR and FeLAR1. A search for homologs in a buckwheat genome database with both sequences returned two more LAR sequences, designated FeLAR2 and FeLAR3. Linkage analysis with an F2 segregating population indicated that the three LAR loci were not genetically linked. We detected high levels of PAs in roots and cotyledons of buckwheat seedlings and in buds and flowers of mature plants. FeANR and FeLAR1-3 were expressed in most organs but had different expression patterns. Our findings would be useful for breeding and further analysis of PA synthesis and its regulation in buckwheat.

Assembly of the draft genome of buckwheat and its applications in identifying agronomically useful genes.

Buckwheat (Fagopyrum esculentum Moench; 2n = 2x = 16) is a nutritionally dense annual crop widely grown in temperate zones. To accelerate molecular breeding programmes of this important crop, we generated a draft assembly of the buckwheat genome using short reads obtained by next-generation sequencing (NGS), and constructed the Buckwheat Genome DataBase. After assembling short reads, we determined 387,594 scaffolds as the draft genome sequence (FES_r1.0). The total length of FES_r1.0 was 1,177,687,305 bp, and the N50 of the scaffolds was 25,109 bp. Gene prediction analysis revealed 286,768 coding sequences (CDSs; FES_r1.0_cds) including those related to transposable elements. The total length of FES_r1.0_cds was 212,917,911 bp, and the N50 was 1,101 bp. Of these, the functions of 35,816 CDSs excluding those for transposable elements were annotated by BLAST analysis. To demonstrate the utility of the database, we conducted several test analyses using BLAST and keyword searches. Furthermore, we used the draft genome as a reference sequence for NGS-based markers, and successfully identified novel candidate genes controlling heteromorphic self-incompatibility of buckwheat. The database and draft genome sequence provide a valuable resource that can be used in efforts to develop buckwheat cultivars with superior agronomic traits.