High-density SNP linkage map construction and QTL mapping for flavonoid-related traits in a tea plant (Camellia sinensis) using 2b-RAD sequencing

BACKGROUND

Flavonoids are important components that confer upon tea plants a unique flavour and health functions. However, the traditional breeding method for selecting a cultivar with a high or unique flavonoid content is time consuming and labour intensive. High-density genetic map construction associated with quantitative trait locus (QTL) mapping provides an effective way to facilitate trait improvement in plant breeding. In this study, an F1 population (LJ43×BHZ) was genotyped using 2b-restriction site-associated DNA (2b-RAD) sequencing to obtain massive single nucleotide polymorphism (SNP) markers to construct a high-density genetic map for a tea plant. Furthermore, QTLs related to flavonoids were identified using our new genetic map.

RESULTS

A total of 13,446 polymorphic SNP markers were developed using 2b-RAD sequencing, and 4,463 of these markers were available for constructing the genetic linkage map. A 1,678.52-cM high-density map at an average interval of 0.40 cM with 4,217 markers, including 427 frameset simple sequence repeats (SSRs) and 3,800 novel SNPs, mapped into 15 linkage groups was successfully constructed. After QTL analysis, a total of 27 QTLs related to flavonoids or caffeine content (CAF) were mapped to 8 different linkage groups, LG01, LG03, LG06, LG08, LG10, LG11, LG12, and LG13, with an LOD from 3.14 to 39.54, constituting 7.5% to 42.8% of the phenotypic variation.

CONCLUSIONS

To our knowledge, the highest density genetic map ever reported was constructed since the largest mapping population of tea plants was adopted in present study. Moreover, novel QTLs related to flavonoids and CAF were identified based on the new high-density genetic map. In addition, two markers were located in candidate genes that may be involved in flavonoid metabolism. The present study provides valuable information for gene discovery, marker-assisted selection breeding and map-based cloning for functional genes that are related to flavonoid content in tea plants.

Construction of a highly saturated Genetic Map for Vitis by Next-generation Restriction Site-associated DNA Sequencing

BACKGROUND

High-saturate molecular linkage maps are an important tool in studies on plant molecular biology and assisted breeding. Development of a large set of single nucleotide polymorphisms (SNPs) via next-generation sequencing (NGS)-based methods, restriction-site associated DNA sequencing (RAD-seq), and the generation of a highly saturated genetic map help improve fine mapping of quantitative trait loci (QTL).

RESULTS

We generated a highly saturated genetic map to identify significant traits in two elite grape cultivars and 176 F₁ plants. In total, 1,426,967 high-quality restriction site-associated DNA tags were detected; 51,365, 23,683, and 70,061 markers were assessed in 19 linkage groups (LGs) for the maternal, paternal, and integrated maps, respectively. Our map was highly saturated in terms of marker density and average "Gap ≤ 5 cM" percentage.

CONCLUSIONS

In this study, RAD-seq of 176 F₁ plants and their parents yielded 8,481,484 SNPs and 1,646,131 InDel markers, of which 65,229 and 4832, respectively, were used to construct a highly saturated genetic map for grapevine. This map is expected to facilitate genetic studies on grapevine, including an evaluation of grapevine and deciphering the genetic basis of economically and agronomically important traits. Our findings provide basic essential genetic data the grapevine genetic research community, which will lead to improvements in grapevine breeding.

A high density SLAF-SNP genetic map and QTL detection for fibre quality traits in Gossypium hirsutum

Background

Upland Cotton (Gossypium hirsutum) is a very important cash crop known for its high quality natural fiber. Recent advances in sequencing technologies provide powerful tools with which to explore the cotton genome for single nucleotide polymorphism marker identification and high density genetic map construction toward more reliable quantitative trait locus mapping.

Results

In the present study, a RIL population was developed by crossing a Chinese high fiber quality cultivar (Yumian 1) and an American high fiber quality line (CA3084), with distinct genetic backgrounds. Specific locus amplified fragment sequencing (SLAF-seq) technology was used to discover SNPs, and a genetic map containing 6254 SNPs was constructed, covering 3141.72 cM with an average distance of 0.5 cM between markers. A total of 95 QTL were detected for fiber quality traits in three environments, explaining 5.5-24.6% of the phenotypic variance. Fifty-five QTL found in multiple environments were considered stable QTL. Nine of the stable QTL were found in all three environments. We identified 14 QTL clusters on 13 chromosomes, each containing one or more stable QTL.

Conclusion

A high-density genetic map of Gossypium hirsutum developed by using specific locus amplified fragment sequencing technology provides detailed mapping of fiber quality QTL, and identification of ‘stable QTL’ found in multiple environments. A marker-rich genetic map provides a foundation for fine mapping, candidate gene identification and marker-assisted selection of favorable alleles at stable QTL in breeding programs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5294-5) contains supplementary material, which is available to authorized users.

QTL mapping of melon fruit quality traits using a high-density GBS-based genetic map

BACKGROUND

Melon shows a broad diversity in fruit morphology and quality, which is still underexploited in breeding programs. The knowledge of the genetic basis of fruit quality traits is important for identifying new alleles that may be introduced in elite material by highly efficient molecular breeding tools.

RESULTS

In order to identify QTLs controlling fruit quality, a recombinant inbred line population was developed using two commercial cultivars as parental lines: "Védrantais", from the cantalupensis group, and "Piel de Sapo", from the inodorus group. Both have desirable quality traits for the market, but their fruits differ in traits such as rind and flesh color, sugar content, ripening behavior, size and shape. We used a genotyping-by-sequencing strategy to construct a dense genetic map, which included around five thousand variants distributed in 824 bins. The RIL population was phenotyped for quality and morphology traits, and we mapped 33 stable QTLs involved in sugar and carotenoid content, fruit and seed morphology and major loci controlling external color of immature fruit and mottled rind. The median confidence interval of the QTLs was 942 kb, suggesting that the high density of the genetic map helped in increasing the mapping resolution. Some of these intervals contained less than a hundred annotated genes, and an integrative strategy combining gene expression and resequencing data enabled identification of candidate genes for some of these traits.

CONCLUSION

Several QTLs controlling fruit quality traits in melon were identified and delimited to narrow genomic intervals, using a RIL population and a GBS-based genetic map.

High-density genetic map construction and mapping of the homologous transformation sterility gene (hts) in wheat using GBS markers

BACKGROUND

Homologous transformation sterility-1 (HTS-1) is a novel wheat mutant that exhibits pistillody, the transformation of stamens into pistils or pistil-like structures. More extreme phenotypes of this mutation can have six pistils or pistil-like structures without any stamens in a floret. Thus, HTS-1 is highly valuable for studies of wheat hybrid breeding and flower development. Previous studies have shown that two major genes (Pis1 and hts) control pistillody in HTS-1. The Pis1 gene controls the three-pistil trait in the three-pistil wheat mutant and has been mapped on chromosome 2D, but the hts gene has not been mapped or identified. To do so, we crossed HTS-1 with CM28TP (three-pistil mutant) and constructed a high-density linkage map with the F₂ population (200 individuals).

RESULTS

The map covered 2779.96 cM, and the genetic distance per chromosome ranged from 37.59 cM to 318.95 cM. The average distance between markers was 1.04 cM. We then mapped hts between GBS-SNP markers 4A_109 and 4A_119, separated by 2.0 cM and 5.2 Mb. To find the candidate genes, the hts region was enlarged to 7.2 Mb, encompassing 752 protein-coding genes. We identified TaWin1 as a possible candidate gene after comparing the 752 genes with 206 common differentially expressed genes between pistillody stamens (PS) versus normal stamens (S) and pistils (P) versus S. Real-time PCR indicated that TaWin1 was highly expressed in HTS-1 during the pistil-and-stamen-differentiating stage, at levels approximately 120 times greater than those in CM28TP. Further analysis indicated that TaWin1 was mainly expressed in HTS-1 PS, supporting its status as a candidate gene of hts. Thus, TaWin1 overexpression probably leads to the transformation of stamens into pistils in wheat.

CONCLUSIONS

The results of this study provide a foundation for further research on stamen and pistil development, with implications for wheat-hybrid breeding programs.

Genetic imprints of domestication for disease resistance, oil quality, and yield component traits in groundnut (Arachis hypogaea L.)

Ploidy difference between wild Arachis species and cultivated genotypes hinder transfer of useful alleles for agronomically important traits. To overcome this genetic barrier, two synthetic tetraploids, viz., ISATGR 1212 (A. duranensis ICG 8123 × A. ipaensis ICG 8206) and ISATGR 265-5A (A. kempff-mercadoi ICG 8164 × A. hoehnei ICG 8190), were used to generate two advanced backcross (AB) populations. The AB-populations, namely, AB-pop1 (ICGV 91114 × ISATGR 1212) and AB-pop2, (ICGV 87846 × ISATGR 265-5A) were genotyped with DArT and SSR markers. Genetic maps were constructed for AB-pop1 and AB-pop2 populations with 258 loci (1415.7 cM map length and map density of 5.5 cM/loci) and 1043 loci (1500.8 cM map length with map density of 1.4 cM/loci), respectively. Genetic analysis identified large number of wild segments in the population and provided a good source of diversity in these populations. Phenotyping of these two populations identified several introgression lines with good agronomic, oil quality, and disease resistance traits. Quantitative trait locus (QTL) analysis showed that the wild genomic segments contributed favourable alleles for foliar disease resistance while cultivated genomic segments mostly contributed favourable alleles for oil quality and yield component traits. These populations, after achieving higher stability, will be useful resource for genetic mapping and QTL discovery for wild species segments in addition to using population progenies in breeding program for diversifying the gene pool of cultivated groundnut.

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri 'DangshanSuli' v1.0 genome

Background

Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far.

Results

To anchor the scaffolds of P. bretschneideri ‘DangshanSuli’ (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, ‘Mantianhong’ (MTH) × ‘Hongxiangsu’ (HXS) and ‘Yuluxiang’ (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross ‘Bayuehong’ (BYH) × ‘Dangshansuli’ (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously.

Conclusions

The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5224-6) contains supplementary material, which is available to authorized users.

High-density genetic map construction and quantitative trait loci identification for growth traits in (Taxodium distichum var. distichum × T. mucronatum) × T. mucronatum

BACKGROUND

'Zhongshanshan' is the general designation for the superior interspecific hybrid clones of Taxodium species, which is widely grown for economic and ecological purposes in southern China. Growth is the priority objective in 'Zhongshanshan' tree improvement. A high-density linkage map is vital to efficiently identify key quantitative trait loci (QTLs) that affect growth.

RESULTS

In total, 403.16 Gb of data, containing 2016,336 paired-end reads, was obtained after preprocessing. The average sequencing depth was 28.49 in T. distichum var. distichum, 25.18 in T. mucronatum, and 11.12 in each progeny. In total, 524,662 high-quality SLAFs were detected, of which 249,619 were polymorphic, and 6166 of the polymorphic markers met the requirements for use in constructing a genetic map. The final map harbored 6156 SLAF markers on 11 linkage groups, and was 1137.86 cM in length, with an average distance of 0.18 cM between adjacent markers. Separate QTL analyses of traits in different years by CIM detected 7 QTLs. While combining multiple-year data, 13 QTLs were detected by ICIM. 5 QTLs were repeatedly detected by the two methods, and among them, 3 significant QTLs (q6-2, q4-2 and q2-1) were detected in at least two traits. Bioinformatic analysis discoveried a gene annotated as a leucine-rich repeat receptor-like kinase gene within q4-2.

CONCLUSIONS

This map is the most saturated one constructed in a Taxodiaceae species to date, and would provide useful information for future comparative mapping, genome assembly, and marker-assisted selection.

High-resolution genetic linkage map of European pear (Pyrus communis) and QTL fine-mapping of vegetative budbreak time

BACKGROUND

Genomic analysis technologies can promote efficient fruit tree breeding. Genotyping by sequencing (GBS) enables generating efficient data for high-quality genetic map construction and QTL analysis in a relatively accessible way. Furthermore, High-resolution genetic map construction and accurate QTL detection can significantly narrow down the putative candidate genes associated with important plant traits.

RESULTS

We genotyped 162 offspring in the F1 'Spadona' x 'Harrow Sweet' pear population using GBS. An additional 21 pear accessions, including the F1 population's parents, from our germplasm collection were subjected to GBS to examine diverse genetic backgrounds that are associated to agriculturally relevant traits and to enhance the power of SNP calling. A standard SNP calling pipeline identified 206,971 SNPs with Asian pear ('Suli') as the reference genome and 148,622 SNPs with the European genome ('Bartlett'). These results enabled constructing a genetic map, after further stringent SNP filtering, consisting of 2036 markers on 17 linkage groups with a length of 1433 cM and an average marker interval of 0.7 cM. We aligned 1030 scaffolds covering a total size of 165.5 Mbp (29%) of the European pear genome to the 17 linkage groups. For high-resolution QTL analysis covering the whole genome, we used phenotyping for vegetative budbreak time in the F1 population. New QTLs associated to vegetative budbreak time were detected on linkage groups 5, 13 and 15. A major QTL on linkage group 8 and an additional QTL on linkage group 9 were confirmed. Due to the significant genotype-by-environment (GxE) effect, we were able to identify novel interaction QTLs on linkage groups 5, 8, 9 and 17. Phenotype-genotype association analysis in the pear accessions for main genotype effect was conducted to support the QTLs detected in the F1 population. Significant markers were detected on every linkage group to which main genotype effect QTLs were mapped.

CONCLUSIONS

This is the first vegetative budbreak study of European pear that makes use of high-resolution genetic mapping. These results provide tools for marker-assisted selection and accurate QTL analysis in pear, and specifically at vegetative budbreak, considering the significant GxE and phenotype-plasticity effects.

Interspecific genetic maps in Miscanthus floridulus and M. sacchariflorus accelerate detection of QTLs associated with plant height and inflorescence

Miscanthus is recognized as a promising lignocellulosic crop for the production of bioethanol and bioproducts worldwide. To facilitate the identification of agronomical important traits and establish genetics knowledge, two genetic maps were developed from a controlled interspecific cross between M. floridulus and M. sacchariflorus. A total of 650 SSR markers were mapped in M. floridulus, spanning 19 linkage groups and 2053.31 cM with an average interval of 3.25 cM. The map of M. sacchariflorus comprised 495 SSR markers in 19 linkage groups covering 1684.86 cM with an average interval of 3.54 cM. The estimation on genome length indicated that the genome coverage of parental genetic maps were 93.87% and 89.91%, respectively. Eighty-eight bi-parental common markers were allowed to connect the two maps, and six pairs of syntenic linkage groups were recognized. Furthermore, quantitative trait loci (QTL) mapping of three agronomic traits, namely, plant height (PH), heading time (HT), and flowering time (FT), demonstrated that a total of 66 QTLs were identified in four consecutive years using interval mapping and multiple-QTL model. The LOD value of these QTLs ranged from 2.51 to 10.60, and the phenotypic variation explained varied from 9.50 to 37.10%. QTL cluster in syntenic groups MF19/MS7 contained six stable QTLs associated with PH, HT, and FT. In conclusion, we report for the first time the genetic mapping of biomass traits in M. floridulus and M. sacchariflorus. These results will be a valuable genetic resource, facilitating the discovery of essential genes and breeding of Miscanthus.

Construction of a high-density genetic map and mapping of QTLs for soybean (Glycine max) agronomic and seed quality traits by specific length amplified fragment sequencing

BACKGROUND

Soybean is not only an important oil crop, but also an important source of edible protein and industrial raw material. Yield-traits and quality-traits are increasingly attracting the attention of breeders. Therefore, fine mapping the QTLs associated with yield-traits and quality-traits of soybean would be helpful for soybean breeders. In the present study, a high-density linkage map was constructed to identify the QTLs for the yield-traits and quality-traits, using specific length amplified fragment sequencing (SLAF-seq).

RESULTS

SLAF-seq was performed to screen SLAF markers with 149 F8:11 individuals from a cross between a semi wild soybean, 'Huapidou', and a cultivated soybean, 'Qihuang26', which generated 400.91 M paired-end reads. In total, 53,132 polymorphic SLAF markers were obtained. The genetic linkage map was constructed by 5111 SLAF markers with segregation type of aa×bb. The final map, containing 20 linkage groups (LGs), was 2909.46 cM in length with an average distance of 0.57 cM between adjacent markers. The average coverage for each SLAF marker on the map was 81.26-fold in the male parent, 45.79-fold in the female parent, and 19.84-fold average in each F8:11 individual. According to the high-density map, 35 QTLs for plant height (PH), 100-seeds weight (SW), oil content in seeds (Oil) and protein content in seeds (Protein) were found to be distributed on 17 chromosomes, and 14 novel QTLs were identified for the first time. The physical distance of 11 QTLs was shorter than 100 Kb, suggesting a direct opportunity to find candidate genes. Furthermore, three pairs of epistatic QTLs associated with Protein involving 6 loci on 5 chromosomes were identified. Moreover, 13, 14, 7 and 9 genes, which showed tissue-specific expression patterns, might be associated with PH, SW, Oil and Protein, respectively.

CONCLUSIONS

With SLAF-sequencing, some novel QTLs and important QTLs for both yield-related and quality traits were identified based on a new, high-density linkage map. Moreover, 43 genes with tissue-specific expression patterns were regarded as potential genes in further study. Our findings might be beneficial to molecular marker-assisted breeding, and could provide detailed information for accurate QTL localization.

Genetic mapping of ovary colour and quantitative trait loci for carotenoid content in the fruit of Cucurbita maxima Duchesne

The high content of carotenoids, sugars, dry matter, vitamins and minerals makes the fruit of winter squash (Cucurbita maxima Duchesne) a valuable fresh-market vegetable and an interesting material for the food industry. Due to their nutritional value, long shelf-life and health protective properties, winter squash fruits have gained increased interest from researchers in recent years. Despite these advantages, the genetic and genomic resources available for C. maxima are still limited. The aim of this study was to use the genetic mapping approach to map the ovary colour locus and to identify the quantitative trait loci (QTLs) for high carotenoid content and flesh colour. An F₆ recombinant inbred line (RIL) mapping population was developed and used for evaluations of ovary colour, carotenoid content and fruit flesh colour. SSR markers and DArTseq genotyping-by-sequencing were used to construct an advanced genetic map that consisted of 1824 molecular markers distributed across linkage groups corresponding to 20 chromosomes of C. maxima. Total map length was 2208 cM and the average distance between markers was 1.21 cM. The locus affecting ovary colour was mapped at the end of chromosome 14. The identified QTLs for carotenoid content in the fruit and fruit flesh colour shared locations on chromosomes 2, 4 and 14. QTLs on chromosomes 2 and 4 were the most meaningful. A correlation was clearly confirmed between fruit flesh colour as described by the chroma value and carotenoid content in the fruit. A high-density genetic map of C. maxima with mapped loci for important fruit quality traits is a valuable resource for winter squash improvement programmes.

Construction of a high-density genetic map and fine QTL mapping for growth and nutritional traits of Crassostrea gigas

BACKGROUND

Both growth and nutritional traits are important economic traits of Crassostrea gigas (C. gigas) in industry. But few work has been done to study the genetic architecture of nutritional traits of the oyster. In this study, we constructed a high-density genetic map of C. gigas to help assemble the genome sequence onto chromosomes, meanwhile explore the genetic basis for nutritional traits via quantitative trait loci (QTL) mapping.

RESULTS

The constructed genetic map contained 5024 evenly distributed markers, with an average marker interval of 0.68 cM, thus representing the densest genetic map produced for the oyster. According to the high collinearity between the consensus map and the oyster genome, 1574 scaffold (about 70%) of the genome sequence of C. gigas were successfully anchored to 10 linkage groups (LGs) of the consensus map. Using this high-qualified genetic map, we then conducted QTL analysis for growth and nutritional traits, the latter of which includes glycogen, amino acid (AA), and fatty acid (FA). Overall, 41 QTLs were detected for 17 traits. In addition, six candidate genes identified in the QTL interval showed significant correlation with the traits on transcriptional levels. These genes include growth-related genes AMY and BMP1, AA metabolism related genes PLSCR and GR, and FA metabolism regulation genes DYRK and ADAMTS.

CONCLUSION

Using the constructed high-qualified linkage map, this study not only assembled nearly 70% of the oyster genome sequence onto chromosomes, but also identified valuable markers and candidate genes for growth and nutritional traits, especially for AA and FA that undergone few studies before. These findings will facilitate genome assembly and molecular breeding of important economic traits in C. gigas.

High-density genetic map construction and quantitative trait loci analysis of the stony hard phenotype in peach based on restriction-site associated DNA sequencing

BACKGROUND

Peach (Prunus persica) is an important fruit crop that generally softens rapidly after harvest resulting in a short shelf-life. By contrast, stony hard (SH) peach fruit does not soften and hardly produces ethylene. To explore the candidate genes responsible for the SH phenotype, a high-density genetic map was constructed by restriction-site associated DNA sequencing technology.

RESULTS

In the present study, the linkage map consisted of 1310 single nucleotide polymorphism markers, spanning 454.2 cM, with an average marker distance of 0.347 cM. The single nucleotide polymorphisms were able to anchor eight linkage groups to their corresponding chromosomes. Based on this high-density integrated peach linkage map and two years of fruit phenotyping, two potential quantitative trait loci for the SH trait were identified and positioned on the genetic map. Additionally, Prupe.6G150900.1, a key gene in abscisic acid (ABA) biosynthesis, displayed a differential expression profile identical to the ABA accumulation pattern: mRNA transcripts were maintained at a high level during storage of SH peaches but occurred at low levels in melting fruit.

CONCLUSION

Thus Prupe.6G150900.1 might play a crucial role in the SH phenotype of peach in which ABA negatively regulates ethylene production. Also, this high-density linkage map of peach will contribute to the mapping of important fruit traits and quantitative trait loci identification.

Genetic linkage map and QTL identification for adventitious rooting traits in red gum eucalypts

The eucalypt species, Eucalyptus tereticornis and Eucalyptus camaldulensis, show tolerance to drought and salinity conditions, respectively, and are widely cultivated in arid and semiarid regions of tropical countries. In this study, genetic linkage map was developed for interspecific cross E. tereticornis × E. camaldulensis using pseudo-testcross strategy with simple sequence repeats (SSRs), intersimple sequence repeats (ISSRs), and sequence-related amplified polymorphism (SRAP) markers. The consensus genetic map comprised totally 283 markers with 84 SSRs, 94 ISSRs, and 105 SRAP markers on 11 linkage groups spanning 1163.4 cM genetic distance. Blasting the SSR sequences against E. grandis sequences allowed an alignment of 64% and the average ratio of genetic-to-physical distance was 1.7 Mbp/cM, which strengths the evidence that high amount of synteny and colinearity exists among eucalypts genome. Blast searches also revealed that 37% of SSRs had homologies with genes, which could potentially be used in the variety of downstream applications including candidate gene polymorphism. Quantitative trait loci (QTL) analysis for adventitious rooting traits revealed six QTL for rooting percent and root length on five chromosomes with interval and composite interval mapping. All the QTL explained 12.0-14.7% of the phenotypic variance, showing the involvement of major effect QTL on adventitious rooting traits. Increasing the density of markers would facilitate the detection of more number of small-effect QTL and also underpinning the genes involved in rooting process.

A high-resolution genetic linkage map and QTL fine mapping for growth-related traits and sex in the Yangtze River common carp (Cyprinus carpio haematopterus)

Background

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection for economic traits in aquaculture species. The Yangtze River common carp (Cyprinus carpio haematopterus) is one of the most important aquacultured strains in China. However, quite limited genetics and genomics resources have been developed for genetic improvement of economic traits in such strain.

Results

A high-resolution genetic linkage map was constructed by using 7820 2b-RAD (2b-restriction site-associated DNA) and 295 microsatellite markers in a F2 family of the Yangtze River common carp (C. c. haematopterus). The length of the map was 4586.56 cM with an average marker interval of 0.57 cM. Comparative genome mapping revealed that a high proportion (70%) of markers with disagreed chromosome location was observed between C. c. haematopterus and another common carp strain (subspecies) C. c. carpio. A clear 2:1 relationship was observed between C. c. haematopterus linkage groups (LGs) and zebrafish (Danio rerio) chromosomes. Based on the genetic map, 21 QTLs for growth-related traits were detected on 12 LGs, and contributed values of phenotypic variance explained (PVE) ranging from 16.3 to 38.6%, with LOD scores ranging from 4.02 to 11.13. A genome-wide significant QTL (LOD = 10.83) and three chromosome-wide significant QTLs (mean LOD = 4.84) for sex were mapped on LG50 and LG24, respectively. A 1.4 cM confidence interval of QTL for all growth-related traits showed conserved synteny with a 2.06 M segment on chromosome 14 of D. rerio. Five potential candidate genes were identified by blast search in this genomic region, including a well-studied multi-functional growth related gene, Apelin.

Conclusions

We mapped a set of suggestive and significant QTLs for growth-related traits and sex based on a high-density genetic linkage map using SNP and microsatellite markers for Yangtze River common carp. Several candidate growth genes were also identified from the QTL regions by comparative mapping. This genetic map would provide a basis for genome assembly and comparative genomics studies, and those QTL-derived candidate genes and genetic markers are useful genomic resources for marker-assisted selection (MAS) of growth-related traits in the Yangtze River common carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4613-1) contains supplementary material, which is available to authorized users.

Mapping of a major quantitative trait locus for bakanae disease resistance in rice by genome resequencing

Bakanae disease (BD) has emerged as a serious threat in almost all rice cultivation regions worldwide. Nampyeong is a Korean japonica rice variety known to be resistant to BD. In this study, quantitative trait locus (QTL) mapping was performed with F2 and F3 plants derived from a cross between the Nampyeong variety and a susceptible Korean japonica line, DongjinAD. First, resequencing of Nampyeong and DongjinAD was performed, which identified 171,035 single nucleotide polymorphisms (SNPs) between the two parental varieties. Using these SNPs, 161 cleaved amplified polymorphic sequence (CAPS) markers and six derived CAPS markers were developed; then, a genetic map was constructed from the genotypes of 180 plants from the DongjinAD/Nampyeong F2 plants. The total length of the constructed genetic map was 1386 cM, with an average interval of 8.9 cM between markers. The BD mortality rates of each F3 family were measured by testing 40 F3 progenies using in vitro seedling screening method. QTL analysis based on the genetic map and mortality rate data revealed a major QTL, qFfR1, on rice chromosome 1. qFfR1 was located at 89.8 cM with a logarithm of the odds (LOD) score of 22.7. Further, there were three markers at this point: JNS01033, JNS01037, and JNS01041. A total of 15 genes were identified with annotations related to defense against plant diseases among the 179 genes in the qFfR1 interval at 95% probability, thereby providing potential candidate genes for qFfR1. qFfR1 and its closely linked markers will be useful in breeding rice varieties resistant to BD.

Gene-centromere mapping in meiotic gynogenetic European seabass

Background

Fully isogenic lines in fish can be developed using “mitotic” gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families.

Results

From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly.

Conclusion

Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3826-z) contains supplementary material, which is available to authorized users.

Genetic complexity of miscanthus cell wall composition and biomass quality for biofuels

BACKGROUND

Miscanthus sinensis is a high yielding perennial grass species with great potential as a bioenergy feedstock. One of the challenges that currently impedes commercial cellulosic biofuel production is the technical difficulty to efficiently convert lignocellulosic biomass into biofuel. The development of feedstocks with better biomass quality will improve conversion efficiency and the sustainability of the value-chain. Progress in the genetic improvement of biomass quality may be substantially expedited by the development of genetic markers associated to quality traits, which can be used in a marker-assisted selection program.

RESULTS

To this end, a mapping population was developed by crossing two parents of contrasting cell wall composition. The performance of 182 F1 offspring individuals along with the parents was evaluated in a field trial with a randomized block design with three replicates. Plants were phenotyped for cell wall composition and conversion efficiency characters in the second and third growth season after establishment. A new SNP-based genetic map for M. sinensis was built using a genotyping-by-sequencing (GBS) approach, which resulted in 464 short-sequence uniparental markers that formed 16 linkage groups in the male map and 17 linkage groups in the female map. A total of 86 QTLs for a variety of biomass quality characteristics were identified, 20 of which were detected in both growth seasons. Twenty QTLs were directly associated to different conversion efficiency characters. Marker sequences were aligned to the sorghum reference genome to facilitate cross-species comparisons. Analyses revealed that for some traits previously identified QTLs in sorghum occurred in homologous regions on the same chromosome.

CONCLUSION

In this work we report for the first time the genetic mapping of cell wall composition and bioconversion traits in the bioenergy crop miscanthus. These results are a first step towards the development of marker-assisted selection programs in miscanthus to improve biomass quality and facilitate its use as feedstock for biofuel production.

High-density genetic map construction and QTLs identification for plant height in white jute (Corchorus capsularis L.) using specific locus amplified fragment (SLAF) sequencing

Background

Genetic mapping and quantitative trait locus (QTL) detection are powerful methodologies in plant improvement and breeding. White jute (Corchorus capsularis L.) is an important industrial raw material fiber crop because of its elite characteristics. However, construction of a high-density genetic map and identification of QTLs has been limited in white jute due to a lack of sufficient molecular markers. The specific locus amplified fragment sequencing (SLAF-seq) strategy combines locus-specific amplification and high-throughput sequencing to carry out de novo single nuclear polymorphism (SNP) discovery and large-scale genotyping. In this study, SLAF-seq was employed to obtain sufficient markers to construct a high-density genetic map for white jute. Moreover, with the development of abundant markers, genetic dissection of fiber yield traits such as plant height was also possible. Here, we present QTLs associated with plant height that were identified using our newly constructed genetic linkage groups.

Results

An F₈ population consisting of 100 lines was developed. In total, 69,446 high-quality SLAFs were detected of which 5,074 SLAFs were polymorphic; 913 polymorphic markers were used for the construction of a genetic map. The average coverage for each SLAF marker was 43-fold in the parents, and 9.8-fold in each F₈ individual. A linkage map was constructed that contained 913 SLAFs on 11 linkage groups (LGs) covering 1621.4 cM with an average density of 1.61 cM per locus. Among the 11 LGs, LG1 was the largest with 210 markers, a length of 406.34 cM, and an average distance of 1.93 cM between adjacent markers. LG11 was the smallest with only 25 markers, a length of 29.66 cM, and an average distance of 1.19 cM between adjacent markers. ‘SNP_only’ markers accounted for 85.54% and were the predominant markers on the map. QTL mapping based on the F₈ phenotypes detected 11 plant height QTLs including one major effect QTL across two cultivation locations, with each QTL accounting for 4.14–15.63% of the phenotypic variance.

Conclusions

To our knowledge, the linkage map constructed here is the densest one available to date for white jute. This analysis also identified the first QTL in white jute. The results will provide an important platform for gene/QTL mapping, sequence assembly, genome comparisons, and marker-assisted selection breeding for white jute.

Construction of a high-density genetic map and the X/Y sex-determining gene mapping in spinach based on large-scale markers developed by specific-locus amplified fragment sequencing (SLAF-seq)

Background

Cultivated spinach (Spinacia oleracea L.) is one of the most widely cultivated types of leafy vegetable in the world, and it has a high nutritional value. Spinach is also an ideal plant for investigating the mechanism of sex determination because it is a dioecious species with separate male and female plants. Some reports on the sex labeling and localization of spinach in the study of molecular markers have surfaced. However, there have only been two reports completed on the genetic map of spinach. The lack of rich and reliable molecular markers and the shortage of high-density linkage maps are important constraints in spinach research work. In this study, a high-density genetic map of spinach based on the Specific-locus Amplified Fragment Sequencing (SLAF-seq) technique was constructed; the sex-determining gene was also finely mapped.

Results

Through bio-information analysis, 50.75 Gb of data in total was obtained, including 207.58 million paired-end reads. Finally, 145,456 high-quality SLAF markers were obtained, with 27,800 polymorphic markers and 4080 SLAF markers were finally mapped onto the genetic map after linkage analysis. The map spanned 1,125.97 cM with an average distance of 0.31 cM between the adjacent marker loci. It was divided into 6 linkage groups corresponding to the number of spinach chromosomes. Besides, the combination of Bulked Segregation Analysis (BSA) with SLAF-seq technology(super-BSA) was employed to generate the linkage markers with the sex-determining gene. Combined with the high-density genetic map of spinach, the sex-determining gene X/Y was located at the position of the linkage group (LG) 4 (66.98 cM–69.72 cM and 75.48 cM–92.96 cM), which may be the ideal region for the sex-determining gene.

Conclusions

A high-density genetic map of spinach based on the SLAF-seq technique was constructed with a backcross (BC₁) population (which is the highest density genetic map of spinach reported at present). At the same time, the sex-determining gene X/Y was mapped to LG4 with super-BSA. This map will offer a suitable basis for further study of spinach, such as gene mapping, map-based cloning of Specific genes, quantitative trait locus (QTL) mapping and marker-assisted selection (MAS). It will also provide an efficient reference for studies on the mechanism of sex determination in other dioecious plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3659-9) contains supplementary material, which is available to authorized users.

Identification of SSR markers closely linked to the yellow seed coat color gene in heading Chinese cabbage (Brassica rapa L. ssp. pekinensis)

Research on the yellow-seeded variety of heading Chinese cabbage will aid in broadening its germplasm resources and lay a foundation for AA genome research in Brassica crops. Here, an F₂ segregating population of 1575 individuals was constructed from two inbred lines (brown-seeded ‘92S105’ and yellow-seeded ‘91-125’). This population was used to identify the linkage molecular markers of the yellow seed coat trait using simple sequence repeat (SSR) techniques combined with a bulk segregant analysis (BSA). Of the 144 SSR primer pairs on the A01-A10 chromosomes from the Brassica database (http://brassicadb.org/brad/), two pairs located on the A06 chromosome showed polymorphic bands between the bulk DNA pools of eight brown-seeded and eight yellow-seeded F₂ progeny. Based on the genome sequence, 454 SSR markers were designed to A06 to detect these polymorphic bands and were synthesized. Six SSR markers linked to the seed coat color gene were successfully selected for fine linkage genetic map construction, in which the two closest flanking markers, SSR449a and SSR317, mapped the Brsc-ye gene to a 40.2 kb region with distances of 0.07 and 0.06 cM, respectively. The molecular markers obtained in this report will assist in the marker-assisted selection and breeding of yellow-seeded lines in Brassica rapa L. and other close species.

Summary: Genetic mapping of the yellow seed coat gene Brsc-ye to a 40.2 kb region of the Chinese cabbage (Brassica rapa ssp. pekinensis) genome will allow marker-assisted selection and breeding of yellow-seeded lines.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Genome-wide development and deployment of informative intron-spanning and intron-length polymorphism markers for genomics-assisted breeding applications in chickpea

The discovery and large-scale genotyping of informative gene-based markers is essential for rapid delineation of genes/QTLs governing stress tolerance and yield component traits in order to drive genetic enhancement in chickpea. A genome-wide 119169 and 110491 ISM (intron-spanning markers) from 23129 desi and 20386 kabuli protein-coding genes and 7454 in silico InDel (insertion-deletion) (1-45-bp)-based ILP (intron-length polymorphism) markers from 3283 genes were developed that were structurally and functionally annotated on eight chromosomes and unanchored scaffolds of chickpea. A much higher amplification efficiency (83%) and intra-specific polymorphic potential (86%) detected by these markers than that of other sequence-based genetic markers among desi and kabuli chickpea accessions was apparent even by a cost-effective agarose gel-based assay. The genome-wide physically mapped 1718 ILP markers assayed a wider level of functional genetic diversity (19-81%) and well-defined phylogenetics among domesticated chickpea accessions. The gene-derived 1424 ILP markers were anchored on a high-density (inter-marker distance: 0.65cM) desi intra-specific genetic linkage map/functional transcript map (ICC 4958×ICC 2263) of chickpea. This reference genetic map identified six major genomic regions harbouring six robust QTLs mapped on five chromosomes, which explained 11-23% seed weight trait variation (7.6-10.5 LOD) in chickpea. The integration of high-resolution QTL mapping with differential expression profiling detected six including one potential serine carboxypeptidase gene with ILP markers (linked tightly to the major seed weight QTLs) exhibiting seed-specific expression as well as pronounced up-regulation especially in seeds of high (ICC 4958) as compared to low (ICC 2263) seed weight mapping parental accessions. The marker information generated in the present study was made publicly accessible through a user-friendly web-resource, "Chickpea ISM-ILP Marker Database". The designing of multiple ISM and ILP markers (2-5 markers/gene) from an individual gene (transcription factor) with numerous aforementioned desirable genetic attributes can widen the user-preference to select suitable primer combination for simultaneous large-scale assaying of functional allelic variation, natural allelic diversity, molecular mapping and expression profiling of genes among chickpea accessions. This will essentially accelerate the identification of functionally relevant molecular tags regulating vital agronomic traits for genomics-assisted crop improvement by optimal resource expenses in chickpea.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping

Background

Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola.

Results

Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed.

Conclusions

We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3079-2) contains supplementary material, which is available to authorized users.