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.

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.

Gene flow signature in the S-allele region of cultivated buckwheat.

BACKGROUND: Buckwheat (Fagopyrum esculentum Moench.) is an annual crop that originated in southern China. The nutritious seeds are used in cooking much like cereal grains. Buckwheat is an outcrossing species with heteromorphic self-incompatibility due to its dimorphic (i.e., short- and long-styled) flowers and intra-morph infertility. The floral morphology and intra-morph incompatibility are both determined by a single S locus. Plants with short-styled flowers are heterozygous (S/s) and plants with long-styled flowers are homozygous recessive (s/s) at this locus, and the S/S genotype is not found. Recently, we built a draft genome assembly of buckwheat and identified the 5.4-Mb-long S-allele region harbored by short-styled plants. In this study, the first report on the genome-wide diversity of buckwheat, we used a genotyping-by-sequencing (GBS) dataset to evaluate the genome-wide nucleotide diversity within cultivated buckwheat landraces worldwide. We also investigated the utility of the S-allele region for phylogenetic analysis of buckwheat. RESULTS: Buckwheat showed high nucleotide diversity (0.0065), comparable to that of other outcrossing plants, based on a genome-wide simple nucleotide polymorphism (SNP) analysis. Phylogenetic analyses based on genome-wide SNPs showed that cultivated buckwheat comprises two groups, Asian and European, and revealed lower nucleotide diversity in the European group (0.0055) and low differentiation between the Asian and European groups. The nucleotide diversity (0.0039) estimated from SNPs in the S-allele region is lower than that in genome-wide SNPs. Phylogenetic analysis based on this region detected three diverged groups, S-1, S-2, and S-3. CONCLUSION: The SNPs detected using the GBS dataset were effective for elucidating the evolutionary history of buckwheat, and led to the following conclusions: (1) the low nucleotide diversity of the entire genome in the European group and low differentiation between the Asian and European groups suggested genetic bottlenecks associated with dispersion from Asia to Europe, and/or recent intensified cultivation and selection in Europe; and (2) the high diversification in the S-allele region was indicative of gene flows from wild to cultivated buckwheat, suggesting that cultivated buckwheat may have multiple origins.

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.

Diversification of 13S globulins, allergenic seed storage proteins, of common buckwheat.

The α polypeptide of the 13S globulin subunit of common buckwheat is the counterpart of the major allergenic ß polypeptide. Trypsin digestibility varies between variants of the α polypeptide with and without a tandem repeat insert. To evaluate the intra-species diversity of 13S globulin, the comprehensive screening of a genomic DNA library was performed, resulting in the isolation of 14 and 3 genes for Met-poor and Met-rich subunits, respectively. Although most tandem repeat units were 45 bp in length, the two-repeat gene Glb2B and all one-repeat genes contained an additional 3 bp. Secondary structure predictions and polyacrylamide gel electrophoresis demonstrated that the sense strand of Glb2B-CCG, the additional 3 bp-deletion clone of Glb2B, formed a more rigid secondary structure than that of the wild-type. Thus, the large intra-species variation of 13S globulin revealed in this study and its diversification might be attributable to the unique nature of the tandem repeat sequences.

S-LOCUS EARLY FLOWERING 3 is exclusively present in the genomes of short-styled buckwheat plants that exhibit heteromorphic self-incompatibility.

The different forms of flowers in a species have attracted the attention of many evolutionary biologists, including Charles Darwin. In Fagopyrum esculentum (common buckwheat), the occurrence of dimorphic flowers, namely short-styled and long-styled flowers, is associated with a type of self-incompatibility (SI) called heteromorphic SI. The floral morphology and intra-morph incompatibility are both determined by a single genetic locus named the S-locus. Plants with short-styled flowers are heterozygous (S/s) and plants with long-styled flowers are homozygous recessive (s/s) at the S-locus. Despite recent progress in our understanding of the molecular basis of flower development and plant SI systems, the molecular mechanisms underlying heteromorphic SI remain unresolved. By examining differentially expressed genes from the styles of the two floral morphs, we identified a gene that is expressed only in short-styled plants. The novel gene identified was completely linked to the S-locus in a linkage analysis of 1,373 plants and had homology to EARLY FLOWERING 3. We named this gene S-LOCUS EARLY FLOWERING 3 (S-ELF3). In an ion-beam-induced mutant that harbored a deletion in the genomic region spanning S-ELF3, a phenotype shift from short-styled flowers to long-styled flowers was observed. Furthermore, S-ELF3 was present in the genome of short-styled plants and absent from that of long-styled plants both in world-wide landraces of buckwheat and in two distantly related Fagopyrum species that exhibit heteromorphic SI. Moreover, independent disruptions of S-ELF3 were detected in a recently emerged self-compatible Fagopyrum species and a self-compatible line of buckwheat. The nonessential role of S-ELF3 in the survival of individuals and the prolonged evolutionary presence only in the genomes of short-styled plants exhibiting heteromorphic SI suggests that S-ELF3 is a suitable candidate gene for the control of the short-styled phenotype of buckwheat plants.

Construction of a BAC library for buckwheat genome research - an application to positional cloning of agriculturally valuable traits.

We have constructed a BAC library for common buckwheat Fagopyrum esculentum Moench. The library includes 142,005 clones with an average insert size of approximately 76 kb, equivalent to approximately a 7 to approximately 8-fold coverage of the genome. Polymerase chain reaction based screening of the library with AGAMOUS and FLORICAULA/LEAFY primers, has identified 7 and 9 BACs, respectively, which are consistent with the genome coverage. This library represents the first large insert genomic library for F. esculentum and it can be served as a genetic resource facilitating agricultural, pharmacological, physiological, and evolutionary studies of the species. To demonstrate the utilization of the library for characterizing agriculturally valuable traits, we developed a sequence tagged site marker tightly linked to the dwarf E locus as well as to the self-incompatibility complex locus and screened the library to initiate positional cloning of the causative genes.

Original birthplace of cultivated common buckwheat inferred from genetic relationships among cultivated populations and natural populations of wild common buckwheat revealed by AFLP analysis.

The genetic relationships among seven cultivated populations and eight natural populations of wild common buckwheat were analyzed using amplified fragment length polymorphism (AFLP). The genetic distance was estimated for each pair of the 15 populations based on the AFLP data, and a phylogenetic tree was constructed using the neighbor-joining (NJ) method based on the genetic distance. All the cultivated populations were grouped in a cluster. The natural populations were grouped into two clusters composed of (1) the Sanjiang group (three populations from eastern Tibet and one population from Adong village of Yunnan province) and (2) two populations from Yunnan province and two populations from Sichuan province. The Sanjiang group is more closely related to cultivated populations. These results indicate that the direct ancestor of common buckwheat was natural populations of wild common buckwheat from the Sanjiang area.

Amplified fragment length polymorphism linkage analysis of common buckwheat (Fagopyrum esculentum) and its wild self-pollinated relative Fagopyrum homotropicum.

Common buckwheat (Fagopyrum esculentum) (2n = 2x = 16) and Fagopyrum homotropicum (2n = 2x = 16) were mated in an interspecific cross and amplified fragment length polymorphism (AFLP) linkage maps were constructed by analyzing segregation in the F2 population. Six hundred and sixty-nine bands were identified using 20 AFLP primer combinations, of which 462 (69%) segregated in the F2 population. The map of F. esculentum has eight linkage groups with 223 markers covering a total of 508.3 cM. The map of F. homotropicum has eight linkage groups with 211 markers covering 548.9 cM. There was one to one correspondence of the esculentum and homotropicum linkage groups. Three morphological markers, distylous self-incompatibility, shattering habit, and winged seed, were located on the AFLP map. Distylous self-incompatibility and shattering habit are tightly linked to each other (1.3 cM) and are located near the center of linkage group 1. Winged seed is located on linkage group 4.