Next Article in Journal
An Integrative Phenotype–Genotype Approach Using Phenotypic Characteristics from the UAE National Diabetes Study Identifies HSD17B12 as a Candidate Gene for Obesity and Type 2 Diabetes
Next Article in Special Issue
Identification of QTNs and Their Candidate Genes for 100-Seed Weight in Soybean (Glycine max L.) Using Multi-Locus Genome-Wide Association Studies
Previous Article in Journal
A Gold Standard, CRISPR/Cas9-Based Complementation Strategy Reliant on 24 Nucleotide Bookmark Sequences
Previous Article in Special Issue
Screening for Resistance to PVY in Australian Potato Germplasm
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 4
10.3390/genes11040459
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat

by
Lijun Zhang
*,
Mingchuan Ma
and
Longlong Liu
Crop Germplasms Resources Research Institute, Shanxi Academy of Agricultural Sciences, Taiyuan 030006, China
*
Author to whom correspondence should be addressed.
Genes 2020, 11(4), 459; https://doi.org/10.3390/genes11040459
Submission received: 26 February 2020 / Revised: 19 April 2020 / Accepted: 20 April 2020 / Published: 23 April 2020
(This article belongs to the Special Issue Selection Methods in Plant Breeding: From Visual Phenotyping to NGS)
Browse Figures
Versions Notes

Abstract

:
As a highly nutritious crop, Tartary buckwheat (Fagopyrum tartaricum) strongly adapts and grows in adverse environments and is widely grown in Asia. However, its flour contains a large proportion of the hull that adheres to the testa layer of the groats and is difficult to be removed in industrial processing. Fortunately, rice-Tartary, with the loose and non-adhering hull, provides potentiality of improving Tartary buckwheat that can dehull easily. Here, we performed high-throughput sequencing for two parents (Tartary buckwheat and rice-Tartary) and two pools (samples from the F2 population) and obtained 101 Gb raw sequencing data for further analysis. Sequencing reads were mapped to the reference genome of Tartary buckwheat, and a total of 633,256 unique SNPs and 270,181 unique indels were found in these four samples. Then, based on the Bulked Segregant Analysis (BSA), we identified a candidate genetic region, containing 45 impact SNPs/indels and 36 genes, that might underly non-adhering hull of rice-Tartary and should have value for breeding easy dehulling Tartary buckwheat.

1. Introduction

The genus Fagopyrum, belonging to the eudicot family Polygonaceae, is comprised of 27 species [1] in which two diploid species are cultivated as food crops currently: Tartary buckwheat (Fagopyrum tataricum) and common buckwheat (Fagopyrum esculentum). In recent years, Tartary buckwheat has received more attention because it outperforms common buckwheat in some respects. Tartary buckwheat has higher and more stable yield attributed to its self-low-seed abortion, compatibility, and tolerance to stress such as frost [2,3]. It is also a beneficial health crop which can produce some nutraceutical products and functional foods and shows a higher content of some nutritional components than common buckwheat, such as total vitamin B content, the flavonoid with high antioxidant activity (e.g., rutin) [4]. Rutin is known for its ability in aiding vitamin C usage, for its strong blood vessels, and for reducing high blood pressure and blood clots [5,6,7,8], and its content in Tartary buckwheat seeds is approximately 100 times (0.8–1.7%) higher than in common buckwheat (0.01%) [9]. However, the hull of Tartary buckwheat tightly adheres to the testa layer of the groats [2], and it is hard to remove hull completely. Therefore, Tartary buckwheat has to be consumed together with the flour inadvertently, which reduces the taste. Most normal mechanical dehulling in industrial processing, such as steam dehulling, is time-consuming, breaking seed, and only up to 70% of the hull can be removed [2]. Therefore, it is one of the main objectives of breeders to develop new easy dehulling Tartary buckwheat varieties to improve the taste feeling of the flour and to meet the market needs.
Rice-Tartary, also called Miqiao in China and bate Phapal in Nepal and India, is a particular Tartary buckwheat type with seeds similar to rice in several aspects including shape, size, and cooking methods [10,11] Unlike seeds of Tartary buckwheat which has three grooves on the hull, rice-Tartary has three length-wise openings or splits in its hull, which makes dehulling easier [12]. Through introgression genetic component of easy dehulling in rice-Tartary into the Tartary buckwheat genome, it is hoped to develop a new buckwheat type with loose or non-adhering hull for ease-dehull [2]. Finding the candidate genes which may relate to the construct of interest is the main task of effort in introgression, and specific genetic component of the acceptor parent (variety to be improved) could be replaced with that of donor parent (variety with desired trait performance). The easy dehulling gene in rice-Tartary has not been reported before, so identification of the genetic locus, cloning of genes, and developing linked markers related to the easy dehulling in rice-Tartary has significant implications for the breeding of Tartary buckwheat.
The traditional gene identification method (positional- or map-based cloning) is time and labor consuming and expensive. It involves developing the genetic population with segregating phenotype of target traits, such as F2, recombinant inbreed populations, near-isogeneic lines, and so on, which need several rounds of crossing and planting [2,13,14]. Avoiding development of genetic population, genome-wide association study (GWAS) exploiting phenotype variations within natural population reduces the time required for the identification of genetic loci of potential traits of interest [15]. However, this method depends on the collection of enough germplasm samples for phenotype variation, genotype at high-density genome-wide markers that are usually generated from chip arrays, or whole genome resequencing which obviously increases the cost [15].
Selective genotyping reduces cost and simplifies analytical processing through a focus on selected samples with extreme phenotypes for genotyping and analyzing [16,17]. Bulked segregant analysis (BSA) [18] and DNA pooling [19] further reduce the cost significantly by bulk samples in two tails of phenotype distribution respectively and analyses them as an integrated unit. Recently, BSA has been modified to locate the target genes using positive markers, so the putative markers do not need to be validated by genotyping entire populations [17,20,21]. As a result, the cost of genotyping has been dramatically reduced from two aspects: (1) testing on fewer selective samples and (2) focusing on fewer positive markers [17,22,23]. For example, considering a population consisting of 300 samples and bulking 30 extreme samples from each tail, BSA expends only 0.6% (~2/300) of total cost required for testing all samples.
In this study, we identified a genetic locus controlling easy de-hulling in rice-Tartary by combining the BSA and high-throughput sequencing to improve Tartary buckwheat with properties of easy dehulling after harvesting. An F2 population was generated using Tartary buckwheat and rice-Tartary as parents, from which 30 rice-Tartary-like lines with easy dehulling seeds and 30 Tartary-buckwheat-like lines with seeds hard to be dehulled were selected and pooled, respectively, forming two “sequencing” lines for genotyping. The identified genetic locus, SNPs or indels, linking the gene underlying the easy dehulling in rice-Tartary, valuable for breeding Tartary buckwheat could be dehulled easily.

2. Materials and Methods

2.1. Sample Collection and DNA Sequencing

In the present study, one Tartary buckwheat sample and one rice-Tartary sample were collected and used as parents to generate an F2 population. A total of 60 samples which contain 30 rice-Tartary-like lines (with easy dehulling seeds) and 30 Tartary-buckwheat-like lines (with hard dehulling seeds) were selected and pooled, respectively (Figure 1).
Genomic DNA was extracted from two parental samples and two pools, and then, four pair-end sequencing libraries with 400 bp insert size were constructed according to the factorial manual. All libraries were sequenced using the Illumina Hiseq (Illumina, San Diego, CA, USA) platform in pair-end model with 150 bp length and 150 sequencing reaction rounds.

2.2. Preprocessing of Raw Sequencing Data

To reduce the effect of sequencing errors, raw sequencing reads were preprocessed to acquire high-quality data. First, reads with adapter contamination were clipped with Adapter Removal (version 2) [24]. Second, low-quality bases were cut using the sliding window method for which window size was settled to 5 bases (PERL program v.5.18.2). The dynamic window slides from reads 5′ to 3′ end and stops at the base where its quality is less than 2 or the average quality of bases in a window is less than 20. Then, bases before that terminate window sliding were preserved as continuous high-quality bases. After that, pair of reads from the same sequencing template (PE) were discarded if one or two of them were shorter than or equal to 50 bases. The remaining reads were regarded as “high-quality clean reads” and were used in further analysis.

2.3. Variant Analysis of SNPs and Indels

The clean reads were mapped to the genomic sequence assembly of Tartary buckwheat [25] using BWA with default parameters [26]. Output results were saved in .bam file and sorted with Picard (v.1.107) (https://www.psc.edu/index.php/user-resources/software/picard). The concordance of PE reads reflected by their mapping were checked using “FixMateInformation” package in Picard. Bias caused by duplicates, in which multiple pairs of reads mapped to exact coordinates in reference, were processed to preserve the read pair with highest mapping score and to discard others, using “MarkeDuplicates” package in Picard. The coverage depth of a base in reference was measured. Based on the mapping results, the SNPs and indels of samples were called with the OTG-snpcaller [27] and UnifiedGenotyper in GATK [28], respectively.

2.4. Identification of Candidate Region by Delta SNP-Index

The SNPs were filtered according to (1) Fisher test of strand bias (FS) < 60; (2) Mapping Quality (MA) > 40; (3) Quality Depth (QD) > 4; and (4) Genotype Quality (GQ) > 20. Indels were filtered according to (1) FS ≤ 200 and (2) QD ≥ 4. SNP-index and delta SNP-index were performed using the method described by Takagi et al. [21] based on the depth information generated in variants calling. The average SNP-index of variants within a sliding window (window size of 1 MB and step size of 50 Kb) were calculated for dot plotting. Based on the results, SNPs with delta SNP-index higher than the threshold (confidence interval > 99%) were picked, and the regions flanking were regarded as the target region that might be underlying easy dehulling of rice-Tartary. The impact effect on gene function of SNPs and indels located in the target region were analyzed using ANNOVAR [29].

3. Results

3.1. DNA Sequencing for Two Parents and Two Pools

We generated approximately 671.62 M reads with a total of 101.42 G bases from four sequencing samples: (1) KF is the rice-Tartary sample as a female parent; (2) KM is the Tartary buckwheat sample as a male parent; (3) K21 is the pool with seeds similar to rice-Tartary from F2 population; and (4) K5 is the pool with seeds similar to Tartary buckwheat from F2 population (Table 1). The number of raw reads generated from those four samples ranges from 163.00 to 175.33 M with a total length ranging from 24.61 to 26.47 Gb. After removing adaptor contamination, low-quality base, and PE with one or two reads shorter than or equal to 50 bp, 624.28 M high-quality clean reads with a total length of 91.17 Gb were contained for further analysis. The clean reads from those four samples range from 151.50 to 164.15 M with total length range from 21.98 to 23.88 Gb.

3.2. Mapping High-Quality Clean Reads to the Reference Genome

The high-quality reads were mapped to the reference genome of Tartary buckwheat [25] using BWA [26]. There were 98.92%, 98.95%, 98.93%, and 98.95% high-quality clean reads for KF, K21, K5, and KM that could be mapped to 88.49%, 88.57%, 88.59%, and 88.39% of the genomic reference sequence, respectively (Table 2). More than ten high-quality clean reads from KF, K21, K5, and KM samples can cover 86.97%, 87.00%, 87.04%, and 86.62% of reference bases, respectively. The average coverage depth of reference contributed by high-quality clean reads from KF, K21, K5, and KM were 33.29%, 32.18%, 32.12%, and 31.40%, respectively (Figure 2 and Table 3).

3.3. SNPs and Indels between Samples and Reference Genome

Based on mapping results, SNPs and indels between samples and reference were called with the OTG-snpcaller [27] and UnifiedGenotyper in GATK [28], respectively. We observed 356,378, 581,386, 582,079, and 544,689 SNPs and 222,008, 261,703, 261,989, and 255,033 indels from KF, K21, K5, and KM, respectively (Table 4). The SNPs (544,689) and indels (255,033) in the KM sample (Tartary buckwheat) might be attributed to the difference between the materials used [25]. Because of the different phenotype of seed husk, more SNPs or indels were expected in KF (rice-Tatary); however, fewer SNPs (356,378) and indels (222,008) were observed in KF (rice-Tartary).
K21 and K5 showed more SNPs and indels than their parents (KF and KM). This may be explained in two aspects. First, K21 and K5 had bulked samples generated by pooling 30 extreme lines in two tails of phenotype distribution, while the KF and KM were single material. Second, progenies combine variations between their parents and the reference genome. There are some sites at which the base in reference is the same as the base in one parent but is different from the base in another parent. The genotype at those sites should be heterozygous in progenies. It is easy to infer that there should be more heterozygous genes in K21 and K5 than in parent samples KM. As expected, we do observe more heterotic sites in K21 and K5. There were 457,839 and 458,763 heterotic SNPs in K21 and K5, which were higher than the 158,972 and 306,428 in KF and KM. There were 74,637 and 74,970 heterotic indels in K21 and K5 that were also higher than the 22,871 and 44,379 in KF and KM.

3.4. Identification of Candidate Region Underlying Easy Dehulling in Rice-Tartary according to SNP-Index

SNP-index and delta SNP-index were calculated and plotted using the method described by Takagi et al. [21] which is based on depth information from SNP calling. As expected, the plot of SNP-index for K21 (rice-Tartary-like seeds) shows an obvious peak (Figure 3A) while that for K5 (Tartary Buckwheat-like seeds) presents the background (Figure 3B). The delta SNP-index, the difference between SNP-index for K5 and K21, magnified the peak at the beginning of chromosome with sequence identifier of “CM008279.1” (Figure 3C). The candidate region near the peak might contain a genetic locus underlying easy dehulling in rice-Tartary was defined which SNPs with delta SNP-index were higher than the threshold (confidence interval > 99%) (blue sector in Figure 4). The candidate region locates at the beginning of “CM008279.1”, ranging from 5,999,388 to 6,856,630 and spanning 857,243 bps. SNPs within this candidate target region were filtered according to their impact on gene function. Forty-four “nonsynonymous” SNPs and one “stop gain” SNP affecting 36 genes were observed. Seven SNPs with significant higher delta SNP-index are located in five genes with annotation of “ABC transporter-like (IPR003439) | AAA + ATPase domain (IPR003593) | ABC transporter, conserved site (IPR017871) | P-loop containing nucleoside triphosphate hydrolase (IPR027417)”; “Pentatricopeptide repeat (IPR002885) | Tetratricopeptide-like helical domain superfamily (IPR011990)”; and “Zinc finger, MIZ-type (IPR004181) | Zinc finger, RING/FYVE/PHD-type (IPR013083)”.

4. Discussion and Conclusions

Tartary buckwheat, as a highly nutritious crop, has attracted increasing attention worldwide, while the hull of Tartary buckwheat is difficult to be removed [1]. Easy dehulling is important for grain processing and high-quality flour production. Here, we identified a candidate region which might be related to easy dehulling in rice-Tartary and could be used for improving postprocessing properties of Tartary buckwheat. A previous study involving crossing Tartary buckwheat with rice-Tartary and performing a progeny testing had concluded that the easy dehulling of rice-Tartary was controlled by a single recessive homozygous gene [2], but the location of the gene has not been reported yet. This might be attributed to the difficulty in crossing Tartary buckwheat with rice-Tartary to develop genetic populations, which is a fundamental basis in the traditional approach for gene mapping, location, and cloning, such as map-based cloning.
Tartary buckwheat is the most difficult species to hybridize artificially in the genus of Fapogyrum because its flower is small (~2mm) and its pollen often dehisces before flower blowing [11]. Rice-Tartary has fewer flowers, and it is difficult to collect enough pollen for artificial pollination, which makes it unsuitable to be used as a male parent. Using rice-Tartary as a female parent is also a daunting task because its flower is close and the sepal is intact after cutting [30]. Therefore, the samples from the F2 population, which were generated by Tartary buckwheat and rice-Tartary, were precious for BSA. The BSA method has been successfully used to detect candidate gene for vital phenotypic characteristics in many other major crops, such as maize [31], barley [32], soybean [33,34], cucumber [35], tomato [36], and chickpea [37]. In this study, an obvious and significant peak in the SNP-index plot for K21 and the final delta SNP-index indicate a single region underlying the easy dehulling in rice-Tartary. This is consistent with results described in the previous publication suggesting a single recessive gene [2]. The SNPs strongly impacting function of genes and genes affected in candidate region encompassing the peak identified in this study have valuable significance for both developing easy dehulling varieties through trait introgression in breeding practice and gene cloning in academic research.
With the increasing interest in Tartary buckwheat products, more researchers focus on the hybridization of Tartary buckwheat. A genetic study involving seed shattering of Tartary buckwheat using intraspecific hybrids method has been reported by Fesenko [38]. In Japan, several laboratories also focus on improving the efficiency of emasculation methods such as treatments of hot water immersion, which has been utilized in crossing of rice [2]. Wang and Campbell report the hybridization between Tartary buckwheat and rice-Tartary [2]. Prospectively, the result acquired in this study, combined with the development of a simple and practicable and straightforward method of hybridization, will speed up Tartary breeding.
In conclusion, we finally identified a candidate genetic region containing 45 impact SNPs/indels and 36 associated genes through comparison and bioinformatic analysis. The candidate genetic region, SNPs/indels, and genes might underly the non-adhering hull phenotype of rice-Tartary and should have value for breeding easy dehulling Tartary buckwheat.

Author Contributions

L.Z. conceived and designed the study. L.L. supervised the study. L.Z., L.L. and M.M. contributed to acquisition, analysis, or interpretation of data. L.Z. wrote the draft report. All authors revised the report and approved the final version before submission. All authors have read and agreed to the published version of the manuscript.

Funding

This program was supported by the National Natural Science Foundation of China (31501372), by the China Agriculture Research System (CARS-07-A-2), and by Shanxi Academy of Agricultural Sciences (YGJPY2004 and YGC2019FZ2).

Conflicts of Interest

The authors declare that there are no conflicts of commercial interest related to this paper.

References

  1. Zhou, M.; Ivan, K.; Sun, H.W.; Nikhil, K.C.; Wieslander, G. Molecular Breeding and Nutritional Aspects of Buckwheat; Academic Press: Pittsburgh, PA, USA, 2016; pp. 203–207. [Google Scholar]
  2. Wang, Y.; Campbell, C.G. Tartary buckwheat breeding (Fagopyrum Tataricum L. Gaertn.) through hybridization with its Rice-Tartary type. Euphytica 2007, 156, 399–405. [Google Scholar] [CrossRef]
  3. Campbell, C.G. Buckwheat: Fagopyrum Esculentum Moench. Promoting the Conservation and Use of Underutilized and Neglected Crops (19); Institute of Plant Genetics and Crop Plant Research, Gatersleben, International Plant Genetic Resources Institute: Rome, Italy, 1994; 95p. [Google Scholar]
  4. Bonafaccia, G.; Marocchini, M.; Kreft, I. Composition and technological properties of the flour and bran from common and Tartary buckwheat. Food Chem. 2003, 80, 9–15. [Google Scholar] [CrossRef]
  5. Middleton, E.; Kandaswami, C.; Theoharides, T.C. The effects of plant flavonoids on mammalian cells: Implications for inflammation, heart disease, and cancer. Pharmacol. Rev. 2000, 52, 673–751. [Google Scholar] [PubMed]
  6. Suzuki, T.; Morishita, T.; Mukasa, Y.; Takigawa, S.; Yokota, S.; Ishiguro, K.; Noda, T. Breeding of ‘Manten-Kirari’, a non-bitter and trace-rutinosidase variety of Tartary buckwheat (Fagopyrum Tataricum Gaertn.). Breed. Sci. 2015, 64, 344–350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Panwar, A.; Nidhi, G.; Chauhan, R. Biosynthesis and accumulation of flavonoids in Fagopyrum spp. Eur. J. Plant Sci. Biotechnol. 2012, 6, 17–26. [Google Scholar]
  8. Nishimura, M.; Ohkawara, T.; Sato, Y.; Satoh, H.; Suzuki, T.; Ishiguro, K.; Noda, T.; Morishita, T.; Nishihira, J. Effectiveness of Rutin-Rich Tartary buckwheat (Fagopyrum Tataricum Gaertn.) ‘Manten-Kirari’ in body weight reduction related to its antioxidant properties: A randomised, double-blind, placebo-controlled study. J. Funct. Foods 2016, 26, 460–469. [Google Scholar] [CrossRef] [Green Version]
  9. Nina, F.; Rode, J.; Wang, Z.; Kreft, I.; Zhang, Z. Tartary Buckwheat (Fagopyrum Tataricum Gaertn.) as a Source of dietary rutin and quercitrin. J. Agric. Food Chem. 2003, 51, 6452–6455. [Google Scholar]
  10. Wang, T.Y.; Yang, K.L.; Lu, P.; Chen, W.P. Evaluation of Tartary buckwheat germplasm from Tibet and their phylogenic study. Buckwheat Trend 1996, 1, 14–21. (In Chinese) [Google Scholar]
  11. Hunt, H.V.; Shang, X.; Jones, M.K. Buckwheat: A crop from outside the major Chinese domestication centres? A review of the archaeobotanical, palynological and genetic evidence. Veg. Hist. Archaeobot. 2018, 27, 493–506. [Google Scholar] [CrossRef] [Green Version]
  12. Ohnishi, O. Geographical distribution of allozymes in natural populations of wild Tartary buckwheat. Fagopyrum 2000, 17, 29–34. [Google Scholar]
  13. Gadaleta, A.; Colasuonno, P.; Giove, S.L.; Blanco, A.; Giancaspro, A. Map-based cloning of QFhb.mgb-2A identifies a WAK2 gene responsible for Fusarium Head Blight resistance in wheat. Sci. Rep. 2013, 126, 1839–1850. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Li, J.; Wang, L.; Zhang, Q.; Liu, Y. Map-based cloning and expression analysis of BMR-6 in sorghum. J. Genet. 2015, 94, 445–452. [Google Scholar] [CrossRef] [PubMed]
  15. Liu, M.; Huang, Q.; Sun, W.; Ma, Z.; Huang, L.; Wu, Q.; Tang, Z.; Bu, T.; Li, C.; Chen, H. Genome-wide investigation of the heat shock transcription factor (Hsf) gene family in Tartary buckwheat (Fagopyrum tataricum). BMC Genom. 2019, 20, 871. [Google Scholar] [CrossRef] [PubMed]
  16. Darvasi, A.; Soller, M. Selective DNA pooling for determination of linkage between a molecular marker and a quantitative trait locus. Genetics 1994, 138, 1365–1373. [Google Scholar]
  17. Sun, Y.; Wang, J.; Crouch, J.H.; Xu, Y. Efficiency of selective genotyping for genetic analysis of complex traits and potential applications in crop improvement. Mol. Breed. 2010, 26, 493–511. [Google Scholar] [CrossRef]
  18. Michelmore, R.W.; Kesseli, R.V. Identification of markers linked to disease-resistance genes by bulked segregant analysis: A rapid method to detect markers in specific genomic regions by using segregating populations. Proc. Natl. Acad. Sci. USA 1991, 88, 9828–9832. [Google Scholar] [CrossRef] [Green Version]
  19. Giovannoni, J.J.; Wing, R.A.; Ganal, M.W.; Tanksley, S.D. Isolation of molecular markers from specific chromosomal intervals using DNA pools from existing mapping populations. Nucleic Acids Res. 1991, 19, 6553–6568. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Xu, Y.; Crouch, J.H. Marker-Assisted Selection in Plant Breeding: From Publications to Practice. Crop. Sci. 2008, 48, 391–407. [Google Scholar] [CrossRef] [Green Version]
  21. Takagi, H.; Abe, A.; Yoshida, K.; Kosugi, S.; Natsume, S.; Mitsuoka, C.; Uemura, A.; Utsushi, H.; Tamiru, M.; Takuno, S.; et al. QTL-Seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. Plant J. 2013, 74, 174–183. [Google Scholar] [CrossRef] [PubMed]
  22. Macgregor, S.; Zhao, Z.Z.; Henders, A.K.; Martin, N.G.; Montgomery, G.W.; Visscher, P.M. Highly cost-efficient genome-wide association studies using DNA pools and dense SNP arrays. Nucleic Acids Res. 2008, 36, e35. [Google Scholar] [CrossRef]
  23. Vikram, P.; Swamy, B.P.M.; Dixit, S.; Ahmed, H.U.; Cruz, M.T.S.; Singh, A.K.; Ye, G.; Kumar, A. Bulk segregant analysis: ‘An effective approach for mapping consistent-effect drought grain yield QTLs in rice’. Field Crop. Res. 2012, 134, 185–192. [Google Scholar] [CrossRef]
  24. Schubert, M.; Stinus, L.; Ludovic, O. AdapterRemoval v2: Rapid adapter trimming, identification, and read merging. BMC Res. Notes 2016, 9, 88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhang, L.; Li, X.; Ma, B.; Gao, Q.; Du, H.; Han, Y.; Li, Y.; Cao, Y.; Qi, M.; Zhu, Y.; et al. The Tartary buckwheat genome provides insights into rutin biosynthesis and abiotic stress tolerance. Mol. Plant 2017, 10, 1224–1237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Li, H.; Durbin, R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Zhu, P.; He, L.; Li, Y.; Huang, W.; Xi, F.; Lin, L.; Zhi, Q.; Zhang, W.; Tang, Y.T.; Geng, C.; et al. OTG-Snpcaller: An optimized pipeline based on TMAP and GATK for SNP calling from Ion Torrent data. PLoS ONE 2015, 10, e0138824. [Google Scholar] [CrossRef] [Green Version]
  28. Mckenna, A.; Hanna, M.; Banks, E.; Sivachenko, A.; Cibulskis, K.; Kernytsky, A.M.; Garimella, K.; Altshuler, D.; Gabriel, S.; Daly, M.J.; et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20, 1297–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Wang, K.; Mingyao, L.; Hakon, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  30. Mukasa, Y.; Suzuki, T.; Honda, Y. Suitability. of Rice-Tartary buckwheat for crossbreeding and for utilization of rutin. Jpn. Agric. Res. Q. 2009, 43, 199–206. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, S.; Yeh, C.-T.; Tang, H.M.; Nettleton, D.; Schnable, P.S. Gene mapping via bulked segregant RNA-Seq (BSR-Seq). PLoS ONE 2012, 7, e36406. [Google Scholar] [CrossRef] [Green Version]
  32. Mascher, M.; Matthias, J.; Kuon, J.E.; Axel, H.; Axel, A.; Sebastian, B.; Scholz, U.; Graner, A.; Stein, N. Mapping-by-sequencing accelerates forward genetics in barley. Genome Boil. 2014, 15, R78. [Google Scholar] [CrossRef] [Green Version]
  33. Campbell, B.W.; Hofstad, A.N.; Sreekanta, S.; Fu, F.; Kono, T.J.Y.; Orourke, J.A.; Vance, C.P.; Muehlbauer, G.J.; Stupar, R.M. Fast neutron-induced structural rearrangements at a soybean NAP1 locus result in gnarled trichomes. Theor. Appl. Genet. 2016, 129, 1725–1738. [Google Scholar] [CrossRef] [Green Version]
  34. Dobbels, A.A.; Michno, J.M.; Campbell, B.W.; Virdi, K.S.; Stec, A.O.; Muehlbauer, G.J.; Naeve, S.L.; Stupar, R.M. An induced chromosomal translocation in soybean disrupts a KASI ortholog and is associated with a high-sucrose and low-oil seed phenotype. Genes 2017, 7, 1215–1223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Lu, H.; Lin, T.; Klein, J.; Wang, S.; Qi, J.; Zhou, Q.; Sun, J.; Zhang, Z.; Weng, Y.; Huang, S. QTL-seq identifies an early flowering QTL located near flowering locus T in cucumber. Theor. Appl. Genet. 2014, 127, 1491–1499. [Google Scholar] [CrossRef]
  36. Illa-Berenguer, E.; van Houten, J.; Huang, Z.; Esther Van, D.K. Rapid and reliable identification of tomato fruit weight and locule number loci by QTL-seq. Theor. Appl. Genet. 2015, 128, 1329–1342. [Google Scholar] [CrossRef]
  37. Das, S.; Upadhyaya, H.D.; Bajaj, D.; Kujur, A.; Badoni, S.; Laxmi; Kumar, V.; Tripathi, S.; Gowda, C.L.L.; Sharma, S.; et al. Deploying QTL-Seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea. DNA Res. 2014, 22, 193–203. [Google Scholar] [CrossRef] [PubMed]
  38. Fesenko, I.V. Non-Shattering Accessions of Fagopyrum tataricum Gaertn. Carry recessive alleles at two loci affecting development of functional abscission layer. Fagopyrum 2006, 23, 7–10. [Google Scholar]
Genes 11 00459 g001 550
Figure 1. Phenotype of traits associated with the dehulling: (A) seed loading branch of Tartary buckwheat as male parent; (B) seed loading branch of rice-Tartary as female parent; (C) seeds of Tartary buckwheat as male parent; (D) seeds of rice-Tartary as female parent; (E) “Tartary-buckwheat-like” seeds in the F2 population; (F) “Tartary-buckwheat-like” seeds in the F2 population; (G) “rice-Tartary-like” seeds in the F2 population; and (H) “rice-Tartary-like” seeds in the F2 population.
Figure 1. Phenotype of traits associated with the dehulling: (A) seed loading branch of Tartary buckwheat as male parent; (B) seed loading branch of rice-Tartary as female parent; (C) seeds of Tartary buckwheat as male parent; (D) seeds of rice-Tartary as female parent; (E) “Tartary-buckwheat-like” seeds in the F2 population; (F) “Tartary-buckwheat-like” seeds in the F2 population; (G) “rice-Tartary-like” seeds in the F2 population; and (H) “rice-Tartary-like” seeds in the F2 population.
Genes 11 00459 g001
Genes 11 00459 g002 550
Figure 2. Coverage depth of reference contributed by clean reads. (A–D) Distribution of coverage depth of reference contributed by KF, K21, K5, and KM, respectively, and (E–H) cumulative distribution of coverage depth of reference contributed by KF, K21, K5, and KM, respectively.
Figure 2. Coverage depth of reference contributed by clean reads. (A–D) Distribution of coverage depth of reference contributed by KF, K21, K5, and KM, respectively, and (E–H) cumulative distribution of coverage depth of reference contributed by KF, K21, K5, and KM, respectively.
Genes 11 00459 g002
Genes 11 00459 g003 550
Figure 3. SNP-index plot: (A) the top panel is the plot of SNP-index for K21; (B) the middle panel is the plot of SNP-index for K5; (C) the bottom panel is the plot of delta SNP-index displaying the difference between A (top panel) and B (middle panel). The green and yellow horizontal lines represent thresholds for 95% and 99% confidence intervals, respectively.
Figure 3. SNP-index plot: (A) the top panel is the plot of SNP-index for K21; (B) the middle panel is the plot of SNP-index for K5; (C) the bottom panel is the plot of delta SNP-index displaying the difference between A (top panel) and B (middle panel). The green and yellow horizontal lines represent thresholds for 95% and 99% confidence intervals, respectively.
Genes 11 00459 g003
Genes 11 00459 g004 550
Figure 4. Candidate target region underlying easy dehulling of rice-Tartary: The panels from outside to inside were (Ⅰ) chromosomes; (II, III) the confidence intervals of 99% and 95% for delta SNP-index; (IV) Delta SNP-index; (V) SNP-index for K5; (VI) SNP-index for K21; and (VII) Candidate region.
Figure 4. Candidate target region underlying easy dehulling of rice-Tartary: The panels from outside to inside were (Ⅰ) chromosomes; (II, III) the confidence intervals of 99% and 95% for delta SNP-index; (IV) Delta SNP-index; (V) SNP-index for K5; (VI) SNP-index for K21; and (VII) Candidate region.
Genes 11 00459 g004
Table
Table 1. Summary of data generated by sequencing.
Table 1. Summary of data generated by sequencing.
SampleRaw ReadsClean High-Quality Reads
Read CountsBase-Count (bp)ReadsBasesGC ContentQ20 RateQ30 Rate
CountsPercentageCounts (bp)Percentage
KF175,326,80626,474,347,706164,150,52093.63%23,881,470,61790.21%37.73%98.77%95.68%
K21167,552,43225,300,417,232157,245,50893.85%22,887,256,30790.46%37.53%98.79%95.74%
K5165,742,38625,027,100,286154,379,48493.14%22,421,430,79089.59%37.62%98.74%95.57%
KM163,004,48424,613,677,084151,502,31292.94%21,978,979,07089.30%37.68%98.72%95.52%
Total671,626,108101,415,542,308624,277,824 91,169,136,784
KF: rice-Tartary sample as female parent; K21: pool of lines with “rice-type” seed from the F2 population; K5: pool of lines with “Tartary-type” seed from the F2 population; KM: Tartary buckwheat sample as male parent.
Table
Table 2. Summary of reads mapped to reference genome.
Table 2. Summary of reads mapped to reference genome.
SampleTotal ReadsAll Mapped ReadsReads with Multiple HitsUniquely Mapped Reads
CountsPercentageCountsPercentagesCountsPercentages
KF166,011,976164,221,55098.92%42,503,32525.60%121,718,22573.32%
K21158,977,495157,302,12898.95%41,559,58426.14%115,742,54472.80%
K5156,090,062154,425,74798.93%38,067,75824.39%116,357,98974.55%
KM153,020,562151,421,20698.95%39,323,84625.70%112,097,36073.26%
Table
Table 3. Summary of reference bases with different coverage depth.
Table 3. Summary of reference bases with different coverage depth.
SampleAverage DepthPercentage of Reference Bases with Different Coverage Depth
Coverage Depth ≥ 1Coverage Depth ≥ 4Coverage Depth ≥ 10
KF33.2988.49%87.90%86.97%
K2132.1888.57%88.00%87.00%
K532.1288.59%88.04%87.04%
KM31.4088.39%87.72%86.62%
Table
Table 4. Summary of SNPs and indel between sample and reference.
Table 4. Summary of SNPs and indel between sample and reference.
SampleSNPIndel
TotalHeterozygousHomozygousTotalHeterozygousHomozygous
KF356,378158,972197,406222,00822,871199,137
K21581,386457,839123,547261,70374,637187,066
K5582,079458,763123,316261,98974,970187,019
KM544,689306,428238,261255,03344,379210,654
Union633,256543,128352,990270,18189,724230,763
Union is the total number of unique SNPs found in all samples.

Share and Cite

MDPI and ACS Style

Zhang, L.; Ma, M.; Liu, L. Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat. Genes 2020, 11, 459. https://doi.org/10.3390/genes11040459

AMA Style

Zhang L, Ma M, Liu L. Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat. Genes. 2020; 11(4):459. https://doi.org/10.3390/genes11040459

Chicago/Turabian Style

Zhang, Lijun, Mingchuan Ma, and Longlong Liu. 2020. "Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat" Genes 11, no. 4: 459. https://doi.org/10.3390/genes11040459

APA Style

Zhang, L., Ma, M., & Liu, L. (2020). Identification of Genetic Locus Underlying Easy Dehulling in Rice-Tartary for Easy Postharvest Processing of Tartary Buckwheat. Genes, 11(4), 459. https://doi.org/10.3390/genes11040459

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop