- Research article
- Open Access
Marker development using SLAF-seq and whole-genome shotgun strategy to fine-map the semi-dwarf gene ari-e in barley
BMC Genomicsvolume 17, Article number: 911 (2016)
Barley semi-dwarf genes have been extensively explored and widely used in barley breeding programs. The semi-dwarf gene ari-e from Golden Promise is an important gene associated with some agronomic traits and salt tolerance. While ari-e has been mapped on barley chromosome 5H using traditional markers and next-generation sequencing technologies, it has not yet been finely located on this chromosome.
We integrated two methods to develop molecular markers for fine-mapping the semi-dwarf gene ari-e: (1) specific-length amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) to develop SNP markers, and (2) the whole-genome shotgun sequence to develop InDels. Both SNP and InDel markers were developed in the target region and used for fine-mapping the ari-e gene. Linkage analysis showed that ari-e co-segregated with marker InDel-17 and was delimited by two markers (InDel-16 and DGSNP21) spanning 6.8 cM in the doubled haploid (DH) Dash × VB9104 population. The genetic position of ari-e was further confirmed in the Hindmarsh × W1 DH population which was located between InDel-7 and InDel-17. As a result, the overlapping region of the two mapping populations flanked by InDel-16 and InDel-17 was defined as the candidate region spanning 0.58 Mb on the POPSEQ physical map.
The current study demonstrated the SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development as an efficient approach to map complex genomic region for isolation of functional gene. The ari-e gene was fine mapped from 10 Mb to 0.58 Mb interval.
Semi-dwarf genes were crucial in the green revolution because they conferred semi-dwarfness, which reduced lodging and increased crop yields [1–3]. Most of the green revolution semi-dwarf varieties resulted from a gene mutation in the gibberellic acid (GA) pathways, e.g. rice sd1 gene, wheat DELLA genes and maize d8 [1–5]. For barley, the world's fourth leading cereal crop, semi-dwarf genes such as sdw1/denso, uzu1 and ari-e have been extensively explored and widely accepted in barley breeding programs in different regions. Barley sdw1/denso varieties carry a mutation in GA 20 oxidase-2, a GA biosynthesis gene that is homologous to rice sd1 [6–9], which has been used to develop short-stature cultivars in western USA, Canada, Australia and European countries [10–14]. The uzu1 mutant had a significantly reduced response to brassinolide resulting from a missense mutation in the BR receptor gene HvBRI1 and has been widely used in East Asia [14–17]. The ari-e gene from Golden Promise is a commonly-used semi-dwarf gene in spring barley in Europe, especially Scotland. A new barley variety with the ari-e gene, Hindmarsh, has been grown in nearly half of Australia. Besides the magnitude of the effect of this gene on reducing plant height, ari-e has significant effects on some agronomic traits, such as grain size, grain composition, malting quality and yield [18, 19]. Golden Promise is a model cultivar for barley genetic transformation [20, 21]. Recently, some studies have reported that ari-e is the only semi-dwarf gene to improve salt tolerance [19, 22–24]. Genetic analysis showed that ari-e was a recessive semi-dwarf gene and was mapped on barley chromosome 5H as a quality character or quantitative trait locus (QTL) using traditional markers and next-generation sequencing technologies [19, 22, 23, 25, 26]. However, the ari-e gene is yet to be identified.
Map-based cloning is the most promising approach for molecular isolation of target genes. The first map-based cloning gene in barley was the powdery mildew resistance gene mlo . A growing number of genes have been isolated from barley based on genetic information. Eight genes are involved in disease resistance against powdery mildew (mlo, Rar1, mla and Ror2) [27–30], stem rust (Rpg1 and Rpg4/Rpg5) [31, 32], barley yellow mosaic virus (rym4/rym5)  and leaf strip (Rdg2a) . Besides resistance genes, some of the genes identified by map-based cloning contribute to complex agronomic traits such as photoperiod sensitivity (Ppd-H1) , vernalization requirement (Vrn3) , row type (Vrs1) , tolerance to boron toxicity (Bot1) , covered/naked caryopsis (nud) , anthocyanin pigmentation (ant2) , floret closing (Cly1) , awn length (lsk2) , increased number of internodes (mnd) , leaf color variation (HvSGRA)  and branch development (Cul4) . A larger number of molecular markers are the basis of map-based cloning for genes related to complex agronomic traits. There are thousands of molecular markers (RFLP, AFLP, SSR, STS, SNP and DArT) available for barley genomes (http://wheat.pw.usda.gov/GG2/index.shtml, http://bioinf.scri.ac.uk/barley_snpdb/index.html, http://www.diversityarrays.com/), allowing an efficient and rapid localization of genes at low resolution. Further narrowing of the target region and distinguishing which gene is responsible for the interested traits is limited by the size of the segregating population and molecular marker density.
Recently, the ongoing revolution in sequencing techniques has significantly reduced costs and advanced technologies to provide new strategies for discovering thousands of markers. Next-generation sequencing (NGS) facilitates the rapid identification of genomic variants for marker development in most species. For example, NGS identified 416,856 markers in wheat  and 11,805 SNP markers in cotton . In barley, 1,391 high confidence SNPs generated from genotyping by sequencing (GBS) were used to construct a linkage map of doubled haploid (DH) Golden Promise × Morex (GPMx) RILs population . Aligning genomic DNA sequences between two sequenced barley varieties, Morex and Barke, identified 436,640 InDel markers .
Bulked segregant analysis (BSA) combined with NGS technologies offers new strategies for marker discovery associated with genes/QTLs. Specific-length amplified fragment sequencing (SLAF-seq) is considered an efficient and high-resolution genotyping method , which has identified major QTLs in maize , soybean , cucumber , rice  and candidate gene isolation in cucumber  and barley . In this study, traditional marker assays and the combination of BSA with SLAF-seq were used to detect a genomic region in barley harboring the semi-dwarf gene ari-e. SNPs and InDel markers were developed by SLAF-seq and whole-genome shotgun sequencing technology, respectively, to fine-map ari-e. The current study exemplifies the use of SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development in the target region, which may pave the way for map-based cloning of the ari-e gene and unraveling the molecular mechanisms of the semi-dwarf phenotype.
Classical mapping of the ari-e gene
According to previous linkage mapping, ari-e is located on chromosome 5H between SSR marker Bmag337 and Bmag357 in Derkado × B83-12/21/5 DH population . Twenty-nine InDel markers in the target region developed by Zhou et al.  were tested in the parents of Dash and VB9104. Of these, six markers (InDel5072, InDel5078, InDel5086, InDel5131, InDel5135 and InDel5136) showed polymorphism. Based on the current barley assembly information, Hvces8 (MLOC_68431), which is homologous to rice cellulose synthase catalytic subunit genes 9 (OsCesA9), is also in the target region. After sequencing the Hvces8 gene in Dash and VB9104, there was 6 bp (GCATCG) insertion in its 5’untranslated region in VB9104. Therefore, the InDel marker of Hvces8 was designed. Two SSR markers were developed from the sequence of Morex_contig47526, one of which was polymorphic between the parents. As a result, seven InDel markers and one SSR marker (Table 1) were used to genotype 119 individuals from the Dash × VB9104 DH population and to construct the linkage map (Fig. 1a). The ari-e gene was mapped as a phenotypic marker and was located to a region on chromosome 5H, delimited by M47526 (1.4 cM) and InDelHvces8 (7.7 cM). The flanking markers (M47526 and InDelHvces8) interval corresponded to 15.6 Mb on the International Barley Sequencing Consortium (IBSC) physical map.
Analysis of SLAF-seq data and SLAF tags
A total of 63,634,776 fragments were procured, each with read lengths of ~100 bp (Table 2). Most of the bases (93.51 %) were high quality, with quality scores of at least 30 (Table 2). The SLAF numbers were 190,519 for Dash and 217,281 for VB9104. The average depth of the SLAF markers was 19.99-fold in Dash, 15.51-fold in VB9104, 30.19-fold in the semi-dwarf pool and 30.34-fold in the highpool. According to the results of SLAF positioning on the barley genome, 319,656 SLAF tags were anchored; the SLAF numbers and chromosome positions are shown in Table 3.
Polymorphic SNPs and association analysis
After aligning the sequence data to the barley reference, SNPs that differed from the reference sequence were identified. At the stage of SNP calling, each SNP supported by at least four reads in semi-dwarf or high pools were filtered out and 193,626 SNP were detected. Of these SNPs, 82,284 SNPs were polymorphic between Dash and VB9104 and were ultimately selected for association analysis. The statistics of marker numbers on each chromosome according to the positioning result are shown in Table 3. The SNP-index was calculated for each SNP. An average SNP-index of SNPs was calculated with 200 SNP-indexes located in a given genomic interval. SNP-index graphs were generated for the semi-dwarf (Fig. 2a) and high (Fig. 2b) pools by plotting the average SNP-index against the position of each sliding window in the barley genome assembly. After combining the SNP-index information into the semi-dwarf and highpools, the ∆(SNP-index) was calculated and plotted against the genome positions (Fig. 2c). According to the IBSC physical map, the candidate region on chromosome 5H between 279,152, 772 and 460,183,011 was obtained with 99 % confidence.
Develop SNP markers based on SLAF-seq
After aligning the marker sequence to the barley physical map, the preliminary mapping region of ari-e, defined by M47526 and InDelHvces8, corresponded to 15.6 Mb (from 387,773,791 to 403,357,465 bp on chromosome 5H of barley IBSC assembly). Combining the preliminary mapping and SLAF-seq technology suggested that the overlap of both results was delimited by M47526 and InDelHvces8 and was the candidate region of the semi-dwarf gene ari-e. To further evaluate the accuracy of SLAF genotyping data, SNPs which showed polymorphism between semi-dwarf and high pools, and were not segregated in each pool, were selected across the entire candidate region. Selected SNPs with one genotype derived from Dash and the other from VB9104 were designated as potential markers. Ultimately, eight SNPs and one InDel (DGSNP76) were identified after performing independent traditional Sanger sequencing in Dash and VB9104. Their primer sequences and physical map positions are shown in Table 4.
Develop InDel markers by whole-genome shotgun strategy
There were 2,029,068,470 clean reads with a read length of at least 100 bp after resequencing the barley cultivar, Hindmarsh. Moreover, 98.64 % (2,001,398,858 reads) were mapped to the recently released barley reference genome assembly by population sequencing (POPSEQ) methodology . About 5.7 million variations, including 412,597 small InDels, were identified in the whole genome. There were 50, 767 InDels on chromosome 5H. According to the POPSEQ assembly, the ari-e candidate region between M47526 and InDelHvces8 spans a total sequence length of approximately 17.5 Mb (478,393,185 to 495,938,660). Therefore, 17 InDels were selected to design primers in the ari-e targeted region (Table 5). The size comparison of PCR products on PAGE gels revealed insertion or deletion polymorphisms between Dash and VB9104. Ultimately, seven InDels (InDel-2, −4, −5, −6, −13, −16 and −17) showed size polymorphism between Dash and VB9104.
Mapping the developed markers in the Dash × VB9104 DH population
Fine-mapping with all of the polymorphic markers was performed in 119 Dash × VB9104 DH individuals (including lines from the semi-dwarf and high pools). After independent traditional Sanger sequencing, the markers derived from different SLAF sequences correlated well with the SLAF-seq genotyping information, suggesting that the markers mined from Dash, VB9104, semi-dwarf and high pools are reliable. InDel markers have unique amplicons and can be distinguished easily on 6 % polyacrylamide gel in the population. Linkage analysis of 119 Dash × VB9104 DH individuals showed that all markers were assigned to the target regions, and 15 of 16 markers were mapped between marker M47526 and InDelHvces8 intervals. Moreover, the InDel-17 marker co-segregated with ari-e. Accordingly, ari-e was delimited to a 6.8 cM region flanked by InDel-16 and DGSNP21 (Fig. 1b). The mapping results also confirmed the genotyping accuracy of SLAF-seq and the efficiency of marker discovery using whole-genome shotgun strategy.
Validation of markers on Hindmarsh × W1 DH population
To confirm the ari-e position, we conducted linkage analysis with 340 DH plants from the Hindmarsh × W1 cross. Five polymorphic markers including M47526, InDel-3, InDel-7, InDel-17 and InDel 5135 were applied to the segregating population for ari-e mapping. Linkage analysis identified the ari-e semi-dwarf gene delimited by two InDel markers InDel-7 (2.4 cM) and InDel-17 (1.4 cM), which were physically located in the region of 482,623,538 to 488,465,184 bp (POPSEQ) on chromosome 5H (Fig. 1c). The mapping position of ari-e in the Hindmarsh × W1 DH population was consistent with that in the Dash × VB9104 DH population and the overlapped region identified by the two populations was considered the ari-e candidate region, which allowed further narrowing of the ari-e locus to an 0.58 Mb interval between InDel-16 (487,890,066) and InDel-17 (488,465,184) on the barley POPSEQ physical map.
Analysis of the candidate region and gene annotation
In this 0.58 Mb region, only five genes were predicted according to the barley POPSEQ physical map, and the gene annotations are presented in Table 6. Mazzucotelli et al.  reported that E3 ubiquitin ligases play a major role in protein degradation and are involved in plant growth and development. Therefore, MLOC_66038 caught our attention, based on the gene annotation of barley which acts as the E3 ubiquitin-protein ligase, and seemed most likely to be responsible for the semi-dwarf phenotype. Sequence analyses of Hindmarsh, however, showed no difference in MLOC_66038 compared with the reference genome Morex. Alignment of the nucleotide sequences in the other four genes also revealed no variation between Hindmarsh and Morex.
Genetic variants of the target region were also identified and are listed in Additional file 1: Table S1. Most of the variants were in the intergenic region, with only a few located in upstream and downstream regions of the predicted genes. Such variants were at least 2.0 Kb and 1.5 Kb away from the upstream and downstream regions of the predicted genes, respectively (Additional file 1: Table S1). To test if the variants affect the expression pattern of the predicted genes, we designed five primer pairs (Additional file 2: Table S2) for qRT-PCR analysis in VB9104, Dash and Golden Promise. According to the qRT-PCR results, the relative expression levels of MLOC_15960, MLOC_72534 and HvLOC100842002 had not significantly changed in the three varieties (Additional file 3: Figure S1). However, HvLOC4347149 and MLOC_66038 exhibited higher expression levels in Golden Promise and lower expression levels in Dash compared with the VB9104 (Additional file 3: Figure S1), which indicates that both genes exhibit no consistent expression with ari-e. It seems that none of the five predicted genes was responsible for the ari-e semi-dwarf phenotype based on our genetic map and the recent release of the barley draft genome assembly.
Traditional methods for developing markers are tedious, time-consuming and expensive . Limited number of markers in the barley genomic region is a major obstacle in the fine-mapping of barley genes, which has encouraged the development of new marker systems that have a greater degree of polymorphism. SLAF-seq is a relatively new experimental method which is considered an efficient and high-resolution strategy for large-scale genotyping using an enhanced reduced representation library (RRL) sequencing method . It is an effective, low-cost technology for constructing high-density linkage maps, which has successfully facilitated the identification of major QTLs related to complex traits in plants [50, 51, 53]. Using the SLAF-seq approach, an SNP-based saturated genetic map was constructed and a major fruit length QTL, which explained 44.6 % of the phenotypic variance in cucumber, was detected . Combining SLAF-seq with BSA identified chromosome regions related to 1000 grain weight in rice and the gene that control inflorescence meristem function in maize [50, 53]. Xu et al.  reported that fruit flesh thickness in cucumber was fine-mapped by two SLAF markers in the 0.19 Mb genomic region. After anchoring the markers to the cucumber 9930 reference genome, 20 genes were predicted, and Csa2M058670.1 was identified as a good candidate gene due to its high expression level in thick fruit flesh varieties. In some cases, the interval identified by SLAF-seq is too large to isolate the target gene. The mapping region of barley HvSGRA (Stage Green-Revertible Albino Gene) defined by SSR and SLAF-seq, was insufficient to pin the candidate gene in barley . Therefore, new markers based on the barley physical map were designed to further fine-mapping and its candidate gene was isolated according to its annotation and a nonsense mutation. In the present study, 319,656 SLAF tags were developed by high-throughput sequencing, with 82,284 SNPs identified (Table 3). Associated analysis with the SNP-index obtained semi-dwarf gene ari-e related candidate regions on chromosome 5H with a size of 181.03 Mb (Fig. 2). The large candidate region of ari-e by SLAF-seq was possibly due to the low recombinants between semi-dwarf and high pools in the target region. After combining with primary mapping, the SNP markers converted from polymorphic SLAFs provide useful data for restricting the candidate-associated regions into 7.2 cM interval, delimiting an approximate 3.6 Mb region in 119 Dash × VB9104 DH lines. Therefore, SLAF-seq is a highly-efficient strategy for mapping the candidate gene, but further work needs to be done to refine the target interval to isolate the candidate, especially in barley with such a large haploid genome.
InDel markers are a gel-based molecular marker which has played a major role in genetic studies and QTL mapping in rice and Arabidopsis [58–60]. The recent release of barley draft assembly sequences provides a new approach to finding insertion/deletion polymorphisms in barley varieties at the DNA level. A total of 436,640 InDels were identified between two sequenced barley varieties, Morex and Barke, based on genome-wide alignment, and 1,140 InDel markers were integrated with the barley consensus map . Taking the same approach, we developed InDel markers in the target region; the newly-designed polymorphic InDel markers were linked with ari-e (Fig. 1b, c). Of these, InDel-17 marker co-segregated with ari-e and we narrowed the ari-e gene to the region flanked by InDel-16 and DGSNP21 in the Dash × VB9104 DH population (Fig. 1b). Most of the developed markers were co-segregated, possibly due to the small population size. Furthermore, we validated the mined markers to another mapping population (Hindmarsh × W1 DH population), and the ari-e gene was flanked by InDel-7 and InDel-17 (Fig. 1c). Combining both mapping results, ari-e can be narrowed between InDel-16 and InDel-17. As a result, alignment of the whole-genome shotgun sequence with the released barley genome sequence enabled the discovery of a large number of DNA markers. This study is an example of the efficient exploitation of targeted InDel markers developed by the whole-genome shotgun strategy in barley. Different alleles at the ari-e locus such as such as ari-e.1, ari-e.119, ari-e.156 and ari-e.228 described by Kucera et al. , were induced by ionizing radiations (neutrons and X-ray) and by chemical mutagens (propyl methanesulfonate and N-methyl-N-nitrosourea). Allelism to the ari-e mutant showed that ari-e.1 and Golden Promise mutant are alleles . The variety Golden Promise (a gamma-ray mutant of the variety Maythorpe) with ari-e is famous for its desirable agronomic traits such as earliness, short stiff straw, reduced awn length and better salt tolerance compared with other semi-dwarf varieties [23, 63]. Considering its importance to barley breeding, ari-e has been genetically characterized in numerous genetic backgrounds. Thomas et al.  used the phenotypic markers in both F2 populations, and a doubled monoploid population to genetically map the ari-e locus to chromosome 5H of barley near the centromere. In the AFLP and SSR marker map constructed for the Derkado × B83-12/21/5 DH population, the ari-e gene was flanked by SSR marker Bmag337 (3 cM) and Bmag357 (7 cM) . After genotyping with the 1,536-SNP Illumina GoldenGate oligonucleotide pool assay in a three-way barley cross, Malosetti et al.  identified that ari-e was close to SNP marker 2_1239 which has been included in the barley integrated map of Aghnoum et al.  in 5H bin 6. Based on Illumina’s oligo pool assays (OPAs) (the same SNP genotyping), ari-e is associated with SNP markers 1_1198 to 2_0449 in 5H bins 2 to 9 in BW043, which resulted from the introgression ari-e into the background of cv. ‘Bowman’ . All these results confirm that barley 5H bin 6 encompasses the major semi-dwarf gene ari-e. Using GBS data of GPMx RILs population, ari-e is closely linked with GBS marker MR_47526P1793R57, which is flanked by MR_335403P1239R45 and MR_1560792P1192F41 . These markers defined a 7.2 cM interval on the GPMx GBS map. According to the barley IBSC physical map, the interval defined at least 46 Mb and contained an estimated 397 genes . The same interval spanned a total sequence length of approximately 10 Mb on the barley POPSEQ map . The authors explained that both the relatively low recombination and lack of detected polymorphism resulted in the low-resolution map in the ari-e region of the GPMx RILs population. In our initial mapping, ari-e was defined at a 9.1 cM interval, which corresponds to 17.5 Mb on the barley POPSEQ physical map. Although the genetic distance was similar, the physical map interval of ari-e defined by M47526 and InDelHvces8 in the Dash × VB9104 DH population was larger than that delimited by MR_335403P1239R45 and MR_1560792P1192F41 in the GPMx RILs population. The target region decreased to 6.8 cM covering the physical map from 487,890,066 bp to 490,942,368 bp intervals in the Dash × VB9104 DH population. After we had validated the polymorphic markers to the Hindmarsh × W1 DH population, ari-e was flanked by InDel-7 and InDel-17 and the genomics DNA interval of 482,623,538 to 488,465,184 bp on chromosome 5H. Based on the barley physical map, its physical position and combining the mapping results of both populations, the candidate region of ari-e declined to an interval of 0.58 Mb, which represents significant progress toward identifying the functional gene compared to previous studies .
Ellis et al.  reported that recessive dwarf mutants in barley indicate a loss of function. It was confirmed in both barley recessive semi-dwarf mutants, uzu1 and sdw1/denso, which resulted from the loss of function of the BR receptor gene HvBRI1 and gibberellin biosynthesis gene HvGA20ox 2 , respectively [8, 15]. Genetic analysis showed that ari-e mutants also carried a recessive allele . We deduced that the semi-dwarfism of ari-e was possibly due to the mutation of one of the functional genes. According to the most current POPSEQ barley map , there are five predicted genes in the candidate region (Table 6). Considering the semi-dwarf phenotype of ari-e and the gene annotation of the five genes, MLOC_66038 predicted as an E3 ubiquitin-protein ligase may be a candidate gene for ari-e. E3 ubiquitin ligase members are involved in the regulation of some biological processes including hormonal control of vegetative growth, plant reproduction, light response, biotic and abiotic stress tolerance, and DNA repair . For example, both Arabidopsis sly1 and rice gid2 encode the F-box subunit of SCF E3 ubiquitin ligase that regulates GA responses and results in the dwarf phenotype [67, 68]. Unfortunately, no nucleotide sequence variants of MLOC_66038 or the other four genes were present between the ari-e variety Hindmarsh and the reference variety Morex.
In barley, semi-dwarf gene sdw1/denso affects not only plant height but also multiple agronomic traits, because sdw1/denso encodes GA 20-oxidase and is involved in gibberellin biosynthesis [6–8]. Therefore, sdw1/denso is expressed all of the tissues and affects many developmental processes through GA level [7, 8]. ari-e also reduces plant height and is related with lots of agronomic traits, such as erect juvenile growth habit, earliness, compact inflorescene, short awns, small seed size and cell size reduction in leaf blades [19, 24, 25]. This indicated that ari-e may be involved in different developmental processes. Furthermore, reduced expression of HvGA20ox 2 was identified in various organs in denso semi-dwarf mutant compared with tall variety . Rieu et al.  reported that AtGA20ox 1 (also named as semi-dwarf gene ga5) was showed different expression in different tissues, and at different developmental stages. The difference of plant height between ari-e semi-dwarf and tall varieties reaches the maximum at the heading time and the uniform leaf samples can easily be collected. Therefore, leaves of VB9104, Dash and Golden Promise were collected to detected relative expression levels of all predicted genes at the heading stage. But expression levels of the five genes were not consistent with ari-e. The unexpected results might be due to the limited resolution of the current barley physical map while the non-coding region may play a key role in inducing semi-dwarfness. Owing to the large size of the barley genome (5.1 Gb) and the high repeat content (80 %), the current physical map comprises 9,265 contigs with a cumulative size of 4.9 Gb representing 96 % of the physical length of the barley genome . Chromosomal assignments between the POPSEQ and IBSC maps agree in 97.6 % of the cases and discordant contig placements mostly occurred in the genetic centromere due to the severely reduced recombination frequency [55, 68]. Moreover, most of the map-based cloned barley genes have been located in the distal regions of the chromosome. Considering that ari-e occurred near the centromere of barley 5H , the repetitive nature of the barley genome [55, 70, 71] and limited resolution of the barley physical map, there may be other unanchored candidate genes for the ari-e semi-dwarf gene.
In summary, next-generation sequencing combined with the current barley physical map served as a hub for marker development in the desired region and the candidate gene will be isolated with the ongoing construction of the high-resolution barley physical map.
Traditional marker assays and SALF-seq in conjunction with BSA were conducted to map the semi-dwarf gene ari-e. SNPs and InDel markers were developed by SLAF-seq and whole-genome shotgun sequencing technology, respectively, to fine-map ari-e. The ari-e gene was mapped between two developed markers InDel-16 and InDel-17 with an interval of 0.58 Mb. The current study exemplifies the use of SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development in the target region, which may pave the way for map-based cloning of the ari-e gene and unraveling the molecular mechanisms of the semi-dwarf phenotype.
Two mapping populations and their parents were used for the genetic analysis and molecular mapping of the ari-e. The first population was comprised of 119 DH lines produced from a cross between Dash and VB9101, and the second DH population was derived from a cross of Hindmarsh × W1 with 340 lines. Hindmarsh with semi-dwarf phenotype, was selected from the Dash × VB9409 cross and used in whole-genome shotgun sequencing for marker development. Dash is a short, stiff straw cultivar with ari-e derived from Golden Promise. Thus, the semi-dwarf phenotype of Hindmarsh can be traced back to Golden Promise. Both VB9104 and W1 are mid-tall accessions. Leaves of VB9104, Dash and Golden Promise were harvested at the heading stage for RNA extraction. Each variety was collected two samples with four individual lines.
Young healthy leaves from the parents and DH individuals were collected, frozen in liquid nitrogen and used for DNA extraction. Genomic DNA was extracted from each plant using the cetyl-trimethyl-ammonium bromide (CTAB) method . DNA was quantified by a Nanodrop 2000 UV-vis spectrophotometer machine and diluted to 50 ng μl−1 with OD260/280 of 1.8–2.2. RNase A (ST579, Beyotime, China) was used to remove RNA contamination.
Classical mapping markers
After integrated linked markers of ari-e from previous mapping reports [19, 26, 64, 66] to the barley consensus map of chromosome 5H , we estimated that the genetic position from 45 cM to 60 cM as the target region of the ari-e gene. The 29 InDel markers in the ari-e gene region reported by Zhou et al.  were used to detect polymorphism between Dash and VB9104. Based on the current barley assembly information, Hvces8 (MLOC_68431) was identified to the target region and was sequenced to find polymorphism between Dash and VB9104. According to Liu et al. , Morex_contig47526 was most closely associated with ari-e. So it was screened for microsatellite DNA; two pairs of potential SSR markers were designed using Oligo Primer Analysis Software v.7. PCR amplified DNA from the parental lines and DH individuals as a template. The final volume was set to 20 μl containing 1 × Taq Mix (GST101, Bioteke, China), 0.3 μM of each primer and 200 ng of template DNA. The PCR reactions were performed as follows: denaturation at 94 °C for 3 min, followed by 35 cycles of 94 °C for 30 s, 60 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 5 min. PCR products were separated on 6 % polyacrylamide gel (acrylamide/bisacrylamide ratio of 37.5:1) in 0.5 × Tris-borate-ethylene-diaminetetra-acetic acid and run at room temperature for 2–4 h, stained with silver nitrate, and observed on white illumination. Size differences in polymorphisms were identified between Dash and VB9104. Primer sequences and amplified lengths of the polymorphism markers used in this study are listed in Table 1.
Construction of SLAF library for sequencing and analysis of SLAF-seq data
We selected 30 progeny of the Dash × VB9104 DH population with the lowest and highest stature and bulked their DNAs to make ‘semi-dwarf’ and ‘high’ pools, respectively. The parents and two pools were used for SLAF library construction and sequencing as described previously . DNA was digested with the restriction enzyme RsaI. Paired-end sequencing with a read length of 100 bp was conducted on the platform of Illumina Genome Analyzer II (Illumina Inc., San Diego, CA, USA). After sequencing, all reads were aligned to the reference sequence released by The International Barley Sequencing Consortium (IBSC) using BWA software . SNP calling was performed using GATK software (https://www.broadinstitute.org/gatk/guide/best-practices.php). We excluded SNPs which supported less than four reads in the two pools and showed no polymorphism between the parents because they may be false positives due to genomic repeat sequence, sequencing or alignment errors. Association mapping was conducted to identify candidate regions for plant height using the SNP-index (see [74, 75] for methods). The Δ(SNP-index) was obtained by subtracting the SNP-index of the semi-dwarf pool from that of the high pool. For each read depth, 99 % confidence intervals of Δ(SNP-index) were obtained following Takagi et al. . Default parameters were used with all software.
Markers development by SLAF-seq strategy and Dash × VB9104 DH population genotyping
To minimize the genetic interval for fine-mapping and to verify the accuracy of SLAF-seq, SNP flanking primers located in the traditional mapping position of ari-e in the Dash × VB9104 DH population were designed using Oligo Primer Analysis Software v.7 which ranged from 100 to 300 bp in length. PCR amplification was conducted as described above. Amplicons were sequenced in one direction using the specific PCR primers distal to the potential SNP position by Shanghai Sunny Biotechnology Co., Ltd. The Megalign program (DNAStar) was used for sequence alignment and to confirm SNP or InDel sites. The confirmed markers were genotyped in 119 Dash × VB9104 DH individuals following SNP marker detection with direct DNA sequencing.
InDel markers developed using the whole-genome shotgun strategy
In this study, the whole-genome shotgun strategy was selected to perform resequencing of the barley cultivar Hindmarsh with ari-e. A paired-end library with insertion size of about 350 bp was prepared and sequenced on the platform of Illumina HiSeqTM2000 in BGI-Shenzhen. After the raw reads had been produced, strict filtering was performed to ensure the reliability of further analysis by removing contaminated or low-quality reads. Clean reads were mapped to the reference genome sequence using the publicly-accepted aligner of BWA with default parameters. The reference used for reads mapping was population sequencing methodology (POPSEQ) . PCR duplication was removed using rmdup in the SAMtools software package. After alignment, only reads with a mapping quality score >40 were used to infer SNP and InDel variations using the classic pipeline of mpileup/bcftools in SAMtools. In this step, variations with less than five supporting reads or those located within 10 bp around both ends were removed. At the same time, the genetic effect was assessed for these variations using SnpEffsoftware based on their gene structure and locations. Moreover, the CDS sequences of these genes containing significant mutations were used to search homology proteins in databases of NR, Swiss-Prot, KEGG, COG and GO for gene function and involved pathway prediction.
Markers validation on the Hindmarsh × W1 DH population
Five markers showing polymorphism between Hindmarsh and W1 were used to genotype 340 individuals from the Hindmarsh × W1 DH population. PCR and genotyping methods are as described above.
Real-time quantitative RT-PCR of the predicted genes
RNA was extracted from the leaves of VB9104, Dash and Golden Promise collected at the heading stage using a Spin Column Plant Total RNA Purification Kit (SK1321, Sanggon Biotech (Shanghai) Co., Ltd). cDNA was prepared from 1 μg RNA using an AMV First Strand cDNA Synthesis Kit (SK2445, Sanggon Biotech (Shanghai) Co., Ltd). qPCR reactions were performed using SYBR Green (B639273, SG Fast qPCR Master Mix (High Rox), BBI) and the Applied Biosystems Step-one plus Real-time PCR System. The real-time PCR assays were performed in triplicate for each cDNA sample. To determine the transcription levels of the five genes in the candidate region, HvACTIN and HvGAPDH were employed as reference genes. Additional file 2: Table S2 lists the oligonucleotide sequences used for quantitative RT-PCR. For statistical analysis of qPCR data, cycle threshold (CT) values were used to determine ∆ CT values (∆ CT = CTtarget –CTreference), and expression levels of target genes relative to reference genes were determined as 2-∆ CT.
Bulked segregant analysis
Genotyping by sequencing
Golden Promise × Morex
International Barley Sequencing Consortium
Polyacrylamide gel electrophoresis
Polymerase chain reaction
Quantitative trait locus
Reduced representation library
Specific-length amplified fragment sequencing
Sasaki A, Ashikari M, Ueguchi-Tanaka M, Itoh H, Nishimura A, Swapan D, et al. A mutant gibberellin-synthesis gene in rice. Nature. 2002;416:701–2.
Spielmeyer W, Ellis MH, Chandler PM. Semidwarf (sd-1), “green revolution” rice, contains a defective gibberellin 20-oxidase gene. Proc Natl Acad Sci U S A. 2002;99(13):9043–8.
Hedden P. The genes of the Green Revolution. Trends Genet. 2003;19(1):5–9.
Peng J, Richards DE, Hartley NM, Murphy GP, Devos KM, Flintham JE, et al. “Green Revolution” genes encode mutant gibberellin response modulators. Nature. 1999;400:256–61.
Milach SCK, Federizzi LC. Dwarfing genes in plant improvement. Adv Agron. 2001;73:35–65.
Jia Q, Zhang J, Westcott S, Zhang X, Bellgard M, Lance R, et al. GA-20 oxidase as a candidate for the semidwarf gene sdw1/denso in barley. Funct Integr Genomics. 2009;9:255–62.
Jia Q, Zhang XQ, Westcott S, Broughton S, Cakir M, Yang J, et al. Expression level of a gibberellins 20-oxidase gene is associated with multiple agronomic and quality traits in barley. Theor Appl Genet. 2011;122:1451–60.
Jia Q, Li C, Shang Y, Zhu J, Hua W, Wang J, et al. Molecular characterization and functional analysis of barley semi-dwarf mutant Riso no. 9265. BMC Genomics. 2015;16:927.
Barboza L, Effgen S, Alonso-Blanco C, Kooke R, Keurentjes JJ, Koornneef M, et al. Arabidopsis semidwarfs evolved from independent mutations in GA20ox1, ortholog to green revolution dwarf alleles in rice and barley. Proc Natl Acad Sci U S A. 2013;110(39):15818–23.
Mickelson HR, Rasmusson DC. Genes for short stature in barley. Crop Sci. 1994;34:1180–3.
Fettell NA, Moody DB, Long N, Flood RG. Determinants of grain size in malting barley. In: Proceedings of the 10th Australian Barley Technical Symposium, 16–20 September 2001, Canberra, ACT, Australia. 2001.
Ivandic V, Malyshev S, Korzum V, Gramer A, Börner A. Comparative mapping of a gibberellic acid-insensitive dwarfing gene (Dwf2) on chromosome 4HS in barley. Theor Appl Genet. 1999;98:728–31.
Hellewell KB, Rasmusson DC, Gallo-Meagher M. Enhancing yield of semi dwarf barley. Crop Sci. 2000;40:352–8.
Zhang J. Inheritance of agronomic traits from the Chinese barley dwarfing gene donors ‘XiaoshanLixiahuang’ and ‘CangzhouLuodamai’. Plant Breed. 2000;119:523–4.
Chono M, Honda I, Zeniya H, Yoneyama K, Saisho D, Takeda K, et al. A semidwarf phenotype of barley uzu results from a nucleotide substitution in the gene encoding a putative brassinosteroid receptor. Plant Physiol. 2003;133:1209–19.
Saisho D, Tanno K, Chono M, Honda I, Kitano H, Takeda K. Spontaneous brassinolide-insensitive barley mutants ‘uzu’ adapted to east Asia. Breeding Sci. 2004;54:409–16.
Zhang J, Zhen L, Zhang CH. Analysis on the dwarfing genes in Zhepi 1 and Aizao 3: two dwarfing gene donors in barley breeding in China. Agri Sci China. 2006;5:643–7.
Thomas WTB, Powell W, Swanston JS. The effects of major genes on quantitatively varying characters in barley. 4. The GPert and denso loci and quality characters. Heredity. 1991;66:381–9.
Ellis RP, Forster BP, Gordon DC, Handley LL, Keith RP, Lawrence P. Phenotype/genotype associations for yield and salt tolerance in a barley mapping population segregating for two dwarfing genes. J Exp Bot. 2002;53:1163–76.
Wan Y, Lemaux PG. Generation of large numbers of independently transformed fertile barley plants. Plant Physiol. 1994;104:37–48.
Tingay S, McElroy D, Kalla R, Fieg S, Wang M, Thornton S, et al. Agrobacterium tumefaciens mediated barley transformation. Plant J. 1997;11:1369–76.
Forster BP, Pakniyat H, Macaylay M, Matheson W, Phillips MS, Thomas WTB, et al. Variation in the leaf sodium content of the Hordeum vulgare cultivar Maythorpe and its derived mutant cv. Golden Promise. Heredity. 1994;73:249–53.
Forster BP. Mutation genetics of salt tolerance in barley: An assessment of Golden Promise and other semi-dwarf mutants. Euphytica. 2001;120:317–28.
Walia H, Wilson C, Condamine P, Ismail AM, Xu J, Cui X, et al. Array-based genotyping and expression analysis of barley cv. Maythorpe and Golden Promise. BMC Genomics. 2007;8:87.
Thomas WTB, Powell W, Wood W. The chromosomal location of the dwarfing gene present in the spring barley variety Golden Promise. Heredity. 1984;53(1):177–83.
Liu H, Bayer M, Druka A, Russel JR, Hackett CA, Poland J, et al. An evaluation of genotyping by sequencing (GBS) to map the Breviatistatum-e (ari-e) locus in cultivated barley. BMC Genomics. 2014;15:104.
Buschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, et al. The barley mlo gene: a novel control element of plant pathogen resistance. Cell. 1997;88:695–705.
Shirasu K, Lahaye T, Tan MW, Zhou FS, Azevedo C, Schulze-Lefert P. A novel class of eukaryotic zinc-binding proteins is required for disease resistance signaling in barley and development in C. elegans. Cell. 1999;99:355–66.
Zhou FS, Kurth KC, Wei FS, Elliot C, Vale G, Yahiaoui N, et al. Cell-autonomous expression of barley Mla1 confers race-specific resistance to the powdery mildew fungus via a Rar1-independent signalling pathway. Plant Cell. 2001;13:337–50.
Collins NC, Thordal-Christensen H, Lipka V, Bau S, Kombrink E, Qiu JL, et al. SNARE-protein-mediated disease resistanceat the plant cell wall. Nature. 2003;425:973–7.
Brueggeman R, Rostoks N, Kudrna D, Kilian A, Han F, Chen J, et al. The barley stem rust-resistance gene Rpg1 is a novel disease-resistance gene with homology to receptor kinases. Proc Natl Acad Sci U S A. 2002;99(14):9328–33.
Brueggeman R, Druka A, Nirmala J, Cavileer T, Drader T, Rostoks N, et al. The stem rust resistance gene Rpg5 encodes a protein with nucleotide-binding-site, leucine-rich, and protein kinase domains. Proc Natl Acad Sci U S A. 2008;105(39):14970–5.
Stein N, Perovic D, Kumlehn J, Pellio B, Stracke S, Streng S, et al. The eukaryotic translation initiation factor 4E confers multiallelic recessive Bymovirus resistance in Hordeum vulgare (L.). Plant J. 2005;42(6):912–22.
Bulgarelli D, Biselli C, Collins NC, Consonni G, Stanca AM, Schulze-Lefert P, et al. The CC-NB-LRR-type Rdg2a resistance gene confers immunity to the seed-borne barley leaf stripe pathogen in the absence of hypersensitive cell death. PloS One. 2010;5(9):e12599.
Turner A, Beales J, Faure S, Dunford RP, Laurie DA. The pseudo-response regulator Ppd-H1 provides adaptation to photoperiod in barley. Science. 2005;310(5750):1031–4.
Yan L, Fu D, Li C, Blechl A, Tranquilli G, Bonafede M, et al. The wheat and barley vernalization gene VRN3 is an orthologue of FT. Proc Natl Acad Sci U S A. 2006;103(51):19581–6.
Komatsuda T, Pourkheirandish M, He C, Azhaguvel P, Kanamori H, Perovic D, et al. Six-rowed barley originated from a mutation in a homeodomain-leucine zipper I-class homeobox gene. Proc Natl Acad Sci U S A. 2007;104(4):1424–9.
Sutton T, Baumann U, Hayes J, Collins NC, Shi BJ, Schnurbusch ST, et al. Boron-toxicity tolerance in barley arising from efflux transporter amplification. Science. 2007;318:14469.
Taketa S, Amano S, Tsujino Y, Sato T, Saisho D, Kakeda K, et al. Barley grain with adhering hulls is controlled by an ERF family transcription factor gene regulating a lipid biosynthesis pathway. Proc Natl Acad Sci U S A. 2008;105(10):4062–7.
Cockram J, White J, Zuluaga DL, Smith D, Comadran J, Macaulay M, et al. Genome-wide association mapping to candidate polymorphism resolution in the unsequenced barley genome. Proc Natl Acad Sci U S A. 2010;107(50):21611–6.
Nair SK, Wang N, Turuspekov Y, Pourkheirandish M, Sinsuwongwat S, Chen G, et al. Cleistogamous flowering in barley arises from the suppression of microRNA-guided HvAP2 mRNA cleavage. Proc Natl Acad Sci U S A. 2010;107:490–5.
You T, Yamashita Y, Kanamori H, Matsumoto T, Lundqvist U, Sato K, et al. A SHORTINTERNODES (SHI) family transcription factor gene regulates awn elongation and pistil morphology in barley. J Exp Bot. 2012;63:5223–32.
Mascher M, Jost M, Kuon JE, Himmelbach A, Aßfalg A, Beier S, et al. Mapping-by-sequencing accelerates forward genetics in barley. Genome Biol. 2014;15(6):R78.
Qin D, Dong J, Xu F, Guo G, Ge S, Xu Q, et al. Characterization and fine mapping of a novel barley Stage Green-Revertible Albino Gene (HvSGRA) by bulked segregant analysis based on SSR assay and specific length amplified fragment sequencing. BMC Genomics. 2015;16:838.
Tavakol E, Okagaki R, Verderio G, Shariati JV, Hussien A, Bilgic H, et al. The barley Uniculme4 gene encodes a BLADE-ON-PETIOLE-like protein that controls tillering and leaf patterning. Plant Physiol. 2015;168(1):164–74.
Saintenac C, Jiang D, Wang S, Akhunov E. Sequence-based mapping of polyploidy wheat genome. G3. 2013;3:1105–14.
Chen W, Yao J, Chu L, Yuan Z, Li H, Zhang Y. Genetic mapping of the nulliplex-branch gene (gb_nb1) in cotton using next-generation sequencing. Theor Appl Genet. 2015;128(3):539–47.
Zhou G, Zhang Q, Tan C, Zhang X, Li CD. Development of the genome-wide InDel markers and their integration with SSR, DArT and SNP markers in single barley map. BMC Genomics. 2015;16:804.
Sun X, Liu D, Zhang X, Li W, Liu H, Hong W, et al. SLAF-seq: an efficient method of large-scale De Novo SNP discovery and genotyping using high-throughput sequencing. PloS One. 2013;8(3):e5870.
Xia C, Chen L, Rong T, Li R, Xiang Y, Wang P, et al. Identification of a new maize inflorescence meristem mutant and association analysis using SLAF-seq method. Euphytica. 2015;202(1):35–44.
Li B, Tian L, Zhang J, Huang L, Han F, Yan S, et al. Construction of a high-density genetic map based on large-scale markers developed by specific length amplified fragment sequencing (SLAF-seq) and its application to QTL analysis for isoflavone content in Glycine max. BMC Genomics. 2014;15(1):1086.
Wei Q, Wang Y, Qin X, Zhang Y, Zhang Z, Wang J, et al. An SNP-based saturated genetic map and QTL analysis of fruit-related traits in cucumber using specific-length amplified fragment (SLAF) sequencing. BMC Genomics. 2014;15(1):1158.
Xu F, Sun X, Chen Y, Huang Y, Tong C, Bao J. Rapid identification of major QTLs associated with rice grain weight and their utilization. PloS One. 2015;10(3):e0122206.
Xu X, Lu L, Zhu B, Xu Q, Qi X, Chen X. QTL mapping of cucumber fruit flesh thickness by SLAF-seq. Sci Rep. 2015;5:15829.
Mascher M, Muehlbauer GJ, Rokhsar DS, Chapman J, Schmutz J, Barry K, et al. Anchoring and ordering NGS contig assemblies by population sequencing (POPSEQ). Plant J. 2013;76:718–27.
Mazzucotelli E, Belloni S, Marone D, De Leonardis AM, Guerra D, Di Fonzo N, et al. The E3 ubiquitin ligase gene family in plants: regulation by degradation. Curr Genomics. 2006;7(8):509–22.
Li CD, Eckstein PE, Lu M, Rossnagel BG, Scoles GJ. Targeted development of a microsatellite marker associated with a true loose smut resistance gene in barley (Hordeumvulgare L.). Mol Breeding. 2001;8:235–42.
Hou XH, Li LC, Peng ZY, Wei BY, Tang SJ, Ding MY, et al. A platform of high-density INDEL/CAPS markers for map-based cloning in Arabidopsis. Plant J. 2010;63(5):880–8.
Păcurar DI, Păcurar ML, Street N, Bussell JD, Pop TI, Gutierrez L, et al. A collection of INDEL markers for map-based cloning in seven Arabidopsis accessions. J Exp Bot. 2012;63(7):2491–501.
Wu DH, Wu HP, Wang CS, Tseng HY, Hwu KK, Wu DH. Genome-wide InDel marker system for application in rice breeding and mapping studies. Euphytica. 2013;192(1):131–43.
Kucera J, Lundqvist U, Gustafsson Å. Inheritance of breviaristatum mutants in barley. Hereditas. 1975;80:263–78.
Franckowiak JD. Allelism tests among selected semi-dwarf barleys. Barley Genet Newsl. 1991;21:17–23.
Pakniyat H, Thomas WTB, Caligari PDS, Forster BP. Comparison of salt tolerance of GPert and non-GPert barleys. Plant Breed. 1997;116:189–91.
Malosetti M, van Eeuwijk FA, Boer MP, Casas AM, Elía M, Moralejo M, et al. Gene and QTL detection in a three-way barley cross under selection by a mixed model with kinship information using SNPs. Theor Appl Genet. 2011;122(8):1605–16.
Aghnoum R, Marcel TC, Johrde A, Pecchioni N, Schweizer P, Niks RE. Basal host resistance of barley to powdery mildew: connecting quantitative trait loci and candidate genes. Mol Plant Microbe Interact. 2010;23:91–102.
Druka A, Franckowiak J, Lundqvist U, Bonar N, Alexander J, Houston K, et al. Genetic dissection of barley morphology and development. Plant Physiol. 2011;155:617–27.
McGinnis KM, Thomas SG, Soule JD, Strader LC, Zale JM, Sun TP, et al. The Arabidopsis SLEEPY1 gene encodes a putative F-Box subunit of an SCF E3 ubiquitin ligase. Plant Cell. 2003;15(5):1120–30.
Gomi K, Sasaki A, Itoh H, Ueguchi-Tankaka M, Ashikari M, Kitano H, et al. GID2, an F-box subunit of the SCF E3 complex, specifically interacts with phosphorylated SLR1 protein and regulates the gibberellin-dependent degradation of SLR1 in rice. Plant J. 2004;37(4):626–34.
Rieu I, Ruiz-Rivero O, Fernandez-Garcia N, Griffiths J, Powers SJ, Gong F, et al. The gibberellin biosynthetic genes AtGA20ox1 and AtGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J. 2008;53:488–504.
Ariyadasa R, Mascher M, Nussbaumer T, Schulte D, Frenkel Z, Poursarebani N, et al. A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms. Plant Physiol. 2014;164(1):412–23.
Mayer K, Waugh R, Brown J, Schulman A, Langridge P, Platzer M, et al. A physical, genetic and functional sequence assembly of the barley genome. Nature. 2012;491:711–6.
Murray MG, Thompson WF. Rapid isolation of high molecular weight plant DNA. Nucl Acids Res. 1980;8(19):4321–5.
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.
Abe A, Kosugi S, Yoshida K, Natsume S, Takagi H, Kanzaki H, et al. Genome sequencing reveals agronomically important loci in rice using MutMap. Nature Biotech. 2012;30:174–8.
Takagi H, Abe A, Yoshida K, Kosugi S, Natsume S, Mitsuoka C, 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–83.
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, et al. MAPMAKER: An interactive computer package for construction primary genetic linkage maps of experimental and natural populations. Genomics. 1987;1:174–81.
Kosambi DD. The estimation of map distances from recombination values. Ann Eugen. 1944;12:172–5.
This work was financially supported by National Natural Science Foundation of China (31471495), China Agriculture Research System (CARS-05), Grain Research & Development Corporation of Australia, Science foundation of Zhejiang Sci-Tech University (16042063-Y) and Young Talent Project of Zhejiang Academy of Agricultural Sciences (2015R16R08E01).
Availability of data and material
The conclusions of this article is included within the article and its additional files.
QJ, JY and CL designed the experiments; QJ and JW genotyped markers of Dash/VB9104 population; JZ contributed to phenotype Dash/VB9104 population. CT performed bioinformatic analysis of whole-genome shotgun sequence data; XZ carried out whole-genome shotgun sequence data analysis and Hindmarsh/W1 phenotyping; HL and SW mapped the markers in Hindmarsh/W1; DM and SB developed Dash/VB9104 population and Hindmarsh/W1 population, respectively. QJ and CL wrote the paper. All authors have read, edited and approved the current version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate