A multi-population consensus genetic map reveals inconsistent marker order among maps likely attributed to structural variations in the apple genome

A high-density integrated map for grapevine based on three mapping populations genotyped by the Vitis18K SNP chip

We present a high-density integrated map for grapevine, allowing refinement and improved understanding of the grapevine genome, while demonstrating the applicability of the Vitis18K SNP chip for linkage mapping. The improvement of grapevine through biotechnology requires identification of the molecular bases of target traits by studying marker-trait associations. The Vitis18K SNP chip provides a useful genotyping tool for genome-wide marker analysis. Most linkage maps are based on single mapping populations, but an integrated map can increase marker density and show order conservation. Here we present an integrated map based on three mapping populations. The parents consist of the well-known wine cultivars 'Cabernet Sauvignon', 'Corvina' and 'Rhine Riesling', the lesser-known wine variety 'Deckrot', and a table grape selection, G1-7720. Three high-density population maps with an average inter-locus gap ranging from 0.74 to 0.99 cM were developed. These maps show high correlations (0.9965-0.9971) with the reference assembly, containing only 93 markers with large order discrepancies compared to expected physical positions, of which a third is consistent across multiple populations. Moreover, the genetic data aid the further refinement of the grapevine genome assembly, by anchoring 104 yet unanchored scaffolds. From these population maps, an integrated map was constructed which includes 6697 molecular markers and reduces the inter-locus gap distance to 0.60 cM, resulting in the densest integrated map for grapevine thus far. A small number of discrepancies, mainly of short distance, involve 88 markers that remain conflictual across maps. The integrated map shows similar collinearity to the reference assembly (0.9974) as the single maps. This high-density map increases our understanding of the grapevine genome and provides a useful tool for its further characterization and the dissection of complex traits.

Consensus genetic map and QTLs for morphological and agronomical traits in mungbean using DArT markers

Mungbean is an important but understudied food legume compared with other major grain crops. Genetic studies through development of high-through put markers, linkage map construction and QTL analysis can accelerate and improve the efficiency of mining for genes for breeding in this crop. This study used four mungbean F₅ recombinant inbred lines (RILs) from crosses of two wild types (ACC 1, ACC 87) and two cultivars (Berken, Kiloga) and DArT markers to construct individual and consensus linkage maps and to identify QTLs associated with 54 traits in mungbean. The number of polymorphic DArT markers identified among the four RIL populations varied from 1062 to 2013. The individual maps covered the lengths of 629.7-883.5 cM, comprising 672-981 DArT markers and 15-19 linkage groups (LG) with average distance between markers of 0.9-1.2 cM. The consensus map had the total length of 795.3 cM, comprising 1539 DArT markers and resolved 11 LGs with an average inter-marker distance of 0.65 cM. Sixty-two QTLs were identified for 39 traits across 10 LGs of the consensus map. Major QTLs were identified for two special traits, late flowering inherited from ACC 1 (6 QTLs, PVE of 11.2-29.9%) and perenniality inherited from ACC 87 (3 QTLs, PVE of 17.4-22.6%) in separate population analysis. Number of congruent QTLs across four mungbean populations and the consensus map was 18 for 13 traits. These results illustrated the high efficiency of DArT marker application in mungbean genetic dissection and suggested the future potential employment of identified QTLs for mungbean improvement.

A high-resolution consensus linkage map for barley based on GBS-derived genotypes

As genotyping-by-sequencing (GBS) is widely used in barley genetic studies, the translation of the physical position of GBS-derived SNPs into accurate genetic positions has become relevant. The main aim of this study was to develop a high-resolution consensus linkage map based on GBS-derived SNPs. The construction of this integrated map involved 11 bi-parental populations composed of 3743 segregating progenies. We adopted a uniform set of SNP-calling and filtering conditions to identify 50 875 distinct SNPs segregating in at least one population. These SNPs were grouped into 18 580 non-redundant SNPs (bins). The resulting consensus linkage map spanned 1050.1 cM, providing an average density of 17.7 bins and 48.4 SNPs per cM. The consensus map is characterized by the absence of large intervals devoid of marker coverage (significant gaps), the largest interval between bins was only 3.7 cM and the mean distance between adjacent bins was 0.06 cM. This high-resolution linkage map will contribute to several applications in genomic research, such as providing useful information on the recombination landscape for QTLs/genes identified via GWAS or ensuring a uniform distribution of SNPs when developing low-cost genotyping tools offering a limited number of markers.

Construction of a High-Density Genetic Map of Acca sellowiana (Berg.) Burret, an Outcrossing Species, Based on Two Connected Mapping Populations

Acca sellowiana, known as feijoa or pineapple guava, is a diploid, (2n = 2x = 22) outcrossing fruit tree species native to Uruguay and Brazil. The species stands out for its highly aromatic fruits, with nutraceutical and therapeutic value. Despite its promising agronomical value, genetic studies on this species are limited. Linkage genetic maps are valuable tools for genetic and genomic studies, and constitute essential tools in breeding programs to support the development of molecular breeding strategies. A high-density composite genetic linkage map of A. sellowiana was constructed using two genetically connected populations: H5 (TCO × BR, N = 160) and H6 (TCO × DP, N = 184). Genotyping by sequencing (GBS) approach was successfully applied for developing single nucleotide polymorphism (SNP) markers. A total of 4,921 SNP markers were identified using the reference genome of the closely related species Eucalyptus grandis, whereas other 4,656 SNPs were discovered using a de novo pipeline. The individual H5 and H6 maps comprised 1,236 and 1,302 markers distributed over the expected 11 linkage groups, respectively. These two maps spanned a map length of 1,593 and 1,572 cM, with an average inter-marker distance of 1.29 and 1.21 cM, respectively. A large proportion of markers were common to both maps and showed a high degree of collinearity. The composite map consisted of 1,897 SNPs markers with a total map length of 1,314 cM and an average inter-marker distance of 0.69. A novel approach for the construction of composite maps where the meiosis information of individuals of two connected populations is captured in a single estimator is described. A high-density, accurate composite map based on a consensus ordering of markers provides a valuable contribution for future genetic research and breeding efforts in A. sellowiana. A novel mapping approach based on an estimation of multipopulation recombination fraction described here may be applied in the construction of dense composite genetic maps for any other outcrossing diploid species.

From Genetic Maps to QTL Cloning: An Overview for Durum Wheat

Durum wheat is one of the most important cultivated cereal crops, providing nutrients to humans and domestic animals. Durum breeding programs prioritize the improvement of its main agronomic traits; however, the majority of these traits involve complex characteristics with a quantitative inheritance (quantitative trait loci, QTL). This can be solved with the use of genetic maps, new molecular markers, phenotyping data of segregating populations, and increased accessibility to sequences from next-generation sequencing (NGS) technologies. This allows for high-density genetic maps to be developed for localizing candidate loci within a few Kb in a complex genome, such as durum wheat. Here, we review the identified QTL, fine mapping, and cloning of QTL or candidate genes involved in the main traits regarding the quality and biotic and abiotic stresses of durum wheat. The current knowledge on the used molecular markers, sequence data, and how they changed the development of genetic maps and the characterization of QTL is summarized. A deeper understanding of the trait architecture useful in accelerating durum wheat breeding programs is envisioned.

Mapping QTLs for Alternaria blight in Linseed (Linum usitatissimum L.)

A SSR-based linkage map of linseed constructed based on 154 individual lines of F 2 mapping population derived from JRF-4 (disease-resistant) and Chambal (disease susceptible) genotypes. QTLs for Alternaria blight and other yield related traits identified. Out of 1720 SSRs, 216 SSRs were found polymorphic among the parents but due to segregation distortion 18 SSRs could not be used for linkage map construction. Total 191 SSRs were used to construct the linkage map and distributed in 15 linkage groups covering genome length of 1802.4 cM. A total of 10 QTLs were identified for 4 phenotypic traits including 4 QTLs for capsules/plant, 2 for capsule weight/plant, 2 for seed weight/plant and 2 for Alternaria blight resistance. This study laid a foundation for further validation and fine mapping with more advance and large set of marker for different QTL identification and marker-assisted selection in linseed.

Maximising recombination across macadamia populations to generate linkage maps for genome anchoring

The Proteaceae genus Macadamia has a recent history of domestication as a commercial nut crop. We aimed to establish the first sequence-based haploid-correlated reference genetic linkage maps for this primarily outcrossing perennial tree crop, with marker density suitable for genome anchoring. Four first generation populations were used to maximise the segregation patterns available within full-sib, biparental and self-pollinated progeny. This allowed us to combine segregation data from overlapping subsets of >4,000 informative sequence-tagged markers to increase the effective coverage of the karyotype represented by the recombinant crossover events detected. All maps had 14 linkage groups, corresponding to the Macadamia haploid chromosome number, and enabled the anchoring and orientation of sequence scaffolds to construct a pseudo-chromosomal genome assembly for macadamia. Comparison of individual maps indicated a high level of congruence, with minor discrepancies satisfactorily resolved within the integrated maps. The combined set of maps significantly improved marker density and the proportion (70%) of the genome sequence assembly anchored. Overall, increasing our understanding of the genetic landscape and genome for this nut crop represents a substantial advance in macadamia genetics and genomics. The set of maps, large number of sequence-based markers and the reconstructed genome provide a toolkit to underpin future breeding that should help to extend the macadamia industry as well as provide resources for the long term conservation of natural populations in eastern Australia of this unique genus.

Development of an integrated 200K SNP genotyping array and application for genetic mapping, genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost-effective manner. Genotyping arrays can provide fast, efficient and high-throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large-scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re-sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from 'Cuiguan' × 'Starkrimson' was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high-density SNP array presents an important tool for rapid and high-throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.

A reference high-density genetic map of greater yam (Dioscorea alata L.)

This study generated the first high-density genetic map for D. alata based on genotyping-by-sequencing and provides new insight on sex determination in yam. Greater yam (Dioscorea alata L.) is a major staple food in tropical and subtropical areas. This study aimed to produce the first reference genetic map of this dioecious species using genotyping-by-sequencing. In this high-density map combining information of two F1 outcrossed populations, 20 linkage groups were resolved as expected and 1579 polymorphic markers were ordered. The consensus map length was 2613.5 cM with an average SNP interval of 1.68 cM. An XX/XY sex determination system was identified on LG6 via the study of sex ratio, homology of parental linkage groups and the identification of a major QTL for sex determination. Homology with the sequenced D. rotundata is described, and the median physical distance between SNPs was estimated at 139.1 kb. The effects of segregation distortion and the presence of heteromorphic sex chromosomes are discussed. This D. alata linkage map associated with the available genomic resources will facilitate quantitative trait mapping, marker-assisted selection and evolutionary studies in the important yet scarcely studied yam species.

High-density multi-population consensus genetic linkage map for peach

Highly saturated genetic linkage maps are extremely helpful to breeders and are an essential prerequisite for many biological applications such as the identification of marker-trait associations, mapping quantitative trait loci (QTL), candidate gene identification, development of molecular markers for marker-assisted selection (MAS) and comparative genetic studies. Several high-density genetic maps, constructed using the 9K SNP peach array, are available for peach. However, each of these maps is based on a single mapping population and has limited use for QTL discovery and comparative studies. A consensus genetic linkage map developed from multiple populations provides not only a higher marker density and a greater genome coverage when compared to the individual maps, but also serves as a valuable tool for estimating genetic positions of unmapped markers. In this study, a previously developed linkage map from the cross between two peach cultivars 'Zin Dai' and 'Crimson Lady' (ZC2) was improved by genotyping additional progenies. In addition, a peach consensus map was developed based on the combination of the improved ZC2 genetic linkage map with three existing high-density genetic maps of peach and a reference map of Prunus. A total of 1,476 SNPs representing 351 unique marker positions were mapped across eight linkage groups on the ZC2 genetic map. The ZC2 linkage map spans 483.3 cM with an average distance between markers of 1.38 cM/marker. The MergeMap and LPmerge tools were used for the construction of a consensus map based on markers shared across five genetic linkage maps. The consensus linkage map contains a total of 3,092 molecular markers, consisting of 2,975 SNPs, 116 SSRs and 1 morphological marker associated with slow ripening in peach (SR). The consensus map provides valuable information on marker order and genetic position for QTL identification in peach and other genetic studies within Prunus and Rosaceae.

Highly dense linkage maps from 31 full-sibling families of turbot (Scophthalmus maximus) provide insights into recombination patterns and chromosome rearrangements throughout a newly refined genome assembly

Highly dense linkage maps enable positioning thousands of landmarks useful for anchoring the whole genome and for analysing genome properties. Turbot is the most important cultured flatfish worldwide and breeding programs in the fifth generation of selection are targeted to improve growth rate, obtain disease resistant broodstock and understand sex determination to control sex ratio. Using a Restriction-site Associated DNA approach, we genotyped 18,214 single nucleotide polymorphism in 1,268 turbot individuals from 31 full-sibling families. Individual linkage maps were combined to obtain a male, female and species consensus maps. The turbot consensus map contained 11,845 markers distributed across 22 linkage groups representing a total normalised length of 3,753.9 cM. The turbot genome was anchored to this map, and scaffolds representing 96% of the assembly were ordered and oriented to obtain the expected 22 megascaffolds according to its karyotype. Recombination rate was lower in males, especially around centromeres, and pairwise comparison of 44 individual maps suggested chromosome polymorphism at specific genomic regions. Genome comparison across flatfish provided new evidence on karyotype reorganisations occurring across the evolution of this fish group.

Effective BAC clone anchoring with genotyping-by-sequencing and Diversity Arrays Technology in a large genome cereal rye

Identification of bacterial artificial chromosome (BAC) clones containing specific sequences is a prerequisite for many applications, such as physical map anchoring or gene cloning. Existing BAC library screening strategies are either low-throughput or require a considerable initial input of resources for platform establishment. We describe a high-throughput, reliable, and cost-effective BAC library screening approach deploying genotyping platforms which are independent from the availability of sequence information: a genotyping-by-sequencing (GBS) method DArTSeq and the microarray-based Diversity Arrays Technology (DArT). The performance of these methods was tested in a very large and complex rye genome. The DArTseq approach delivered superior results: a several fold higher efficiency of addressing genetic markers to BAC clones and anchoring of BAC clones to genetic map and also a higher reliability. Considering the sequence independence of the platform, the DArTseq-based library screening can be proposed as an attractive method to speed up genomics research in resource poor species.

Genetic mapping and comparative genomics to inform restoration enhancement and culture of southern flounder, Paralichthys lethostigma

BACKGROUND

Southern flounder, Paralichthys lethostigma, historically support a substantial fishery along the Atlantic and Gulf coasts of the southern United States. Low year-class strengths over the past few years in the western Gulf of Mexico have raised concern that spawning stocks may be overfished. Current management of the resource includes releasing hatchery-raised juveniles to restock bays and estuaries; additionally, there is a growing interest in the potential for commercial aquaculture of the species. Currently, genomic resources for southern flounder do not exist. Here, we used two hatchery-reared families and double-digest, restriction-site-associated DNA (ddRAD) sequencing to create a reduced-representation genomic library consisting of several thousand single nucleotide polymorphisms (SNPs) located throughout the genome.

RESULTS

The relative position of each SNP-containing locus was determined to create a high-density genetic map spanning the 24 linkage groups of the southern flounder genome. The consensus map was used to identify regions of shared synteny between southern flounder and seven other fish species for which genome assemblies are available. Finally, syntenic blocks were used to localize genes identified from transcripts in European flounder as potentially being involved in ecotoxicological and osmoregulatory responses, as well as QTLs associated with growth and disease resistance in Japanese flounder, on the southern flounder linkage map.

CONCLUSIONS

The information provided by the linkage map will enrich restoration efforts by providing a foundation for interpreting spatial genetic variation within the species, ultimately furthering an understanding of the adaptive potential and resilience of southern flounder to future changes in local environmental conditions. Further, the map will facilitate the use of genetic markers to enhance restoration and commercial aquaculture.

Genotyping-by-sequencing application on diploid rose and a resulting high-density SNP-based consensus map

Roses, which have been cultivated for at least 5000 years, are one of the most important ornamental crops in the world. Because of the interspecific nature and high heterozygosity in commercial roses, the genetic resources available for rose are limited. To effectively identify markers associated with QTL controlling important traits, such as disease resistance, abundant markers along the genome and careful phenotyping are required. Utilizing genotyping by sequencing technology and the strawberry genome (Fragaria vesca v2.0.a1) as a reference, we generated thousands of informative single nucleotide polymorphism (SNP) markers. These SNPs along with known bridge simple sequence repeat (SSR) markers allowed us to create the first high-density integrated consensus map for diploid roses. Individual maps were first created for populations J06-20-14-3×"Little Chief" (J14-3×LC), J06-20-14-3×"Vineyard Song" (J14-3×VS) and "Old Blush"×"Red Fairy" (OB×RF) and these maps were linked with 824 SNPs and 13 SSR bridge markers. The anchor SSR markers were used to determine the numbering of the rose linkage groups. The diploid consensus map has seven linkage groups (LGs), a total length of 892.2 cM, and an average distance of 0.25 cM between 3527 markers. By combining three individual populations, the marker density and the reliability of the marker order in the consensus map was improved over a single population map. Extensive synteny between the strawberry and diploid rose genomes was observed. This consensus map will serve as the tool for the discovery of marker-trait associations in rose breeding using pedigree-based analysis. The high level of conservation observed between the strawberry and rose genomes will help further comparative studies within the Rosaceae family and may aid in the identification of candidate genes within QTL regions.

Construction of an ultra-high density consensus genetic map, and enhancement of the physical map from genome sequencing in Lupinus angustifolius

An ultra-high density genetic map containing 34,574 sequence-defined markers was developed in Lupinus angustifolius. Markers closely linked to nine genes of agronomic traits were identified. A physical map was improved to cover 560.5 Mb genome sequence. Lupin (Lupinus angustifolius L.) is a recently domesticated legume grain crop. In this study, we applied the restriction-site associated DNA sequencing (RADseq) method to genotype an F₉ recombinant inbred line population derived from a wild type × domesticated cultivar (W × D) cross. A high density linkage map was developed based on the W × D population. By integrating sequence-defined DNA markers reported in previous mapping studies, we established an ultra-high density consensus genetic map, which contains 34,574 markers consisting of 3508 loci covering 2399 cM on 20 linkage groups. The largest gap in the entire consensus map was 4.73 cM. The high density W × D map and the consensus map were used to develop an improved physical map, which covered 560.5 Mb of genome sequence data. The ultra-high density consensus linkage map, the improved physical map and the markers linked to genes of breeding interest reported in this study provide a common tool for genome sequence assembly, structural genomics, comparative genomics, functional genomics, QTL mapping, and molecular plant breeding in lupin.

High density SNP and DArT-based genetic linkage maps of two closely related oil palm populations

Oil palm (Elaeis guineensis Jacq.) is an outbreeding perennial tree crop with long breeding cycles, typically 12 years. Molecular marker technologies can greatly improve the breeding efficiency of oil palm. This study reports the first use of the DArTseq platform to genotype two closely related self-pollinated oil palm populations, namely AA0768 and AA0769 with 48 and 58 progeny respectively. Genetic maps were constructed using the DArT and SNP markers generated in combination with anchor SSR markers. Both maps consisted of 16 major independent linkage groups (2n = 2× = 32) with 1399 and 1466 mapped markers for the AA0768 and AA0769 populations, respectively, including the morphological trait "shell-thickness" (Sh). The map lengths were 1873.7 and 1720.6 cM with an average marker density of 1.34 and 1.17 cM, respectively. The integrated map was 1803.1 cM long with 2066 mapped markers and average marker density of 0.87 cM. A total of 82% of the DArTseq marker sequence tags identified a single site in the published genome sequence, suggesting preferential targeting of gene-rich regions by DArTseq markers. Map integration of higher density focused around the Sh region identified closely linked markers to the Sh, with D.15322 marker 0.24 cM away from the morphological trait and 5071 bp from the transcriptional start of the published SHELL gene. Identification of the Sh marker demonstrates the robustness of using the DArTseq platform to generate high density genetic maps of oil palm with good genome coverage. Both genetic maps and integrated maps will be useful for quantitative trait loci analysis of important yield traits as well as potentially assisting the anchoring of genetic maps to genomic sequences.

Development of CACTA transposon derived SCAR markers and their use in population structure analysis in Zea mays

Molecular marker technologies have proven to be an important breakthrough for genetic studies, construction of linkage maps and population genetics analysis. Transposable elements (TEs) constitute major fractions of repetitive sequences in plants and offer a wide range of possible areas to be explored as molecular markers. Sequence characterized amplified region (SCAR) marker development provides us with a simple and time saving alternative approach for marker development. We employed the CACTA-TD to develop SCARs and then integrated them into linkage map and used them for population structure and genetic diversity analysis of corn inbred population. A total of 108 dominant SCAR markers were designed out of which, 32 were successfully integrated in to the linkage map of maize RIL population and the remaining were added to a physical map for references to check the distribution throughout all chromosomes. Moreover, 76 polymorphic SCARs were used for diversity analysis of corn accessions being used in Korean corn breeding program. The overall average polymorphic information content (PIC) was 0.34, expected heterozygosity was 0.324 and Shannon's information index was 0.491 with a percentage of polymorphism of 98.67%. Further analysis by associating with desirable traits may also provide some accurate trait specific tagged SCAR markers. TE linked SCARs can provide an added level of polymorphism as well as improved discriminating ability and therefore can be useful in further breeding programs to develop high yielding germplasm.

Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives

There is a rapidly rising trend in the development and application of molecular marker assays for gene mapping and discovery in field crops and trees. Thus far, more than 50 SNP arrays and 15 different types of genotyping-by-sequencing (GBS) platforms have been developed in over 25 crop species and perennial trees. However, much less effort has been made on developing ultra-high-throughput and cost-effective genotyping platforms for applied breeding programs. In this review, we discuss the scientific bottlenecks in existing SNP arrays and GBS technologies and the strategies to develop targeted platforms for crop molecular breeding. We propose that future practical breeding platforms should adopt automated genotyping technologies, either array or sequencing based, target functional polymorphisms underpinning economic traits, and provide desirable prediction accuracy for quantitative traits, with universal applications under wide genetic backgrounds in crops. The development of such platforms faces serious challenges at both the technological level due to cost ineffectiveness, and the knowledge level due to large genotype-phenotype gaps in crop plants. It is expected that such genotyping platforms will be achieved in the next ten years in major crops in consideration of (a) rapid development in gene discovery of important traits, (b) deepened understanding of quantitative traits through new analytical models and population designs, (c) integration of multi-layer -omics data leading to identification of genes and pathways responsible for important breeding traits, and (d) improvement in cost effectiveness of large-scale genotyping. Crop breeding chips and genotyping platforms will provide unprecedented opportunities to accelerate the development of cultivars with desired yield potential, quality, and enhanced adaptation to mitigate the effects of climate change.

Integrated high-density consensus genetic map of Pyrus and anchoring of the 'Bartlett' v1.0 (Pyrus communis) genome

Genetic maps are essential tools for pear genetics and genomics research. In this study, we first constructed an integrated simple sequence repeat (SSR) and single nucleotide polymorphism (SNP)-based consensus genetic map for pear based on common SSR markers between nine published maps. A total of 5,085 markers, including 1,232 SSRs and 3,853 SNPs, were localized on a consensus map spanning 3,266.0 cM in total, with an average marker interval of 0.64 cM, which represents the highest density consensus map of pear to date. Using three sets of high-density SNP-based genetic maps with European pear genetic backgrounds, we anchored a total of 291.5 Mb of the ‘Bartlett’ v1.0 (Pyrus communis L.) genome scaffolds into 17 pseudo-chromosomes. This accounted for 50.5% of the genome assembly, which was a great improvement on the 29.7% achieved originally. Intra-genome and inter-genome synteny analyses of the new ‘Bartlett’ v1.1 genome assembly with the Asian pear ‘Dangshansuli’ (Pyrus bretschneideri Rehd.) and apple (Malus × domestica Borkh.) genomes uncovered four new segmental duplication regions. The integrated high-density SSR and SNP-based consensus genetic map provided new insights into the genetic structure patterns of pear and assisted in the genome assembly of ‘Bartlett’ through further exploration of different pear genetic maps.

High-quality de novo assembly of the apple genome and methylome dynamics of early fruit development

Using the latest sequencing and optical mapping technologies, we have produced a high-quality de novo assembly of the apple (Malus domestica Borkh.) genome. Repeat sequences, which represented over half of the assembly, provided an unprecedented opportunity to investigate the uncharacterized regions of a tree genome; we identified a new hyper-repetitive retrotransposon sequence that was over-represented in heterochromatic regions and estimated that a major burst of different transposable elements (TEs) occurred 21 million years ago. Notably, the timing of this TE burst coincided with the uplift of the Tian Shan mountains, which is thought to be the center of the location where the apple originated, suggesting that TEs and associated processes may have contributed to the diversification of the apple ancestor and possibly to its divergence from pear. Finally, genome-wide DNA methylation data suggest that epigenetic marks may contribute to agronomically relevant aspects, such as apple fruit development.

Genotyping-by-Sequencing Facilitates a High-Density Consensus Linkage Map for Aegilops umbellulata, a Wild Relative of Cultivated Wheat

High-density genetic maps are useful to precisely localize QTL or genes that might be used to improve traits of nutritional and/or economical importance in crops. However, high-density genetic maps are lacking for most wild relatives of crop species, including wheat. Aegilops umbellulata is a wild relative of wheat known for its potential as a source of biotic and abiotic stress resistance genes. In this work, we have developed a framework consensus genetic map using two biparental populations derived from accessions PI 298905, PI 542369, PI 5422375, and PI 554395. The framework map comprised 3009 genotype-by-sequence SNPs with a total map size of 948.72 cM. On average, there were three SNPs per centimorgan for each chromosome. Chromosome 1U was the shortest (66.5 cM), with only 81 SNPs, whereas the remaining chromosomes had between 391 and 591 SNP markers. A total of 2395 unmapped SNPs were added to the linkage maps through a recombination frequency approach, and increased the number of SNPs placed on the consensus map to a total of 5404 markers. Segregation distortion was disproportionally high for chromosome 1U for both populations used to construct component linkage maps, and thus segregation distortion could be one of the probable reasons for the exceptionally reduced linkage size for chromosome 1U. From comparative analysis, Ae. umbellulata chromosomes except 4U showed moderate to strong collinearity with corresponding homeologous chromosomes of hexaploid wheat and barley. The present consensus map may serve as a reference map in QTL mapping and validation projects, and also in genome assembly to develop a reference genome sequence for Ae. umbellulata.

Linkage Mapping and Comparative Genomics of Red Drum (Sciaenops ocellatus) Using Next-Generation Sequencing

Developments in next-generation sequencing allow genotyping of thousands of genetic markers across hundreds of individuals in a cost-effective manner. Because of this, it is now possible to rapidly produce dense genetic linkage maps for nonmodel species. Here, we report a dense genetic linkage map for red drum, a marine fish species of considerable economic importance in the southeastern United States and elsewhere. We used a prior microsatellite-based linkage map as a framework and incorporated 1794 haplotyped contigs derived from high-throughput, reduced representation DNA sequencing to produce a linkage map containing 1794 haplotyped restriction-site associated DNA (RAD) contigs, 437 anonymous microsatellites, and 44 expressed sequence-tag-linked microsatellites (EST-SSRs). A total of 274 candidate genes, identified from transcripts from a preliminary hydrocarbon exposure study, were localized to specific chromosomes, using a shared synteny approach. The linkage map will be a useful resource for red drum commercial and restoration aquaculture, and for better understanding and managing populations of red drum in the wild.

A High-Density Consensus Map of Common Wheat Integrating Four Mapping Populations Scanned by the 90K SNP Array

A high-density consensus map is a powerful tool for gene mapping, cloning and molecular marker-assisted selection in wheat breeding. The objective of this study was to construct a high-density, single nucleotide polymorphism (SNP)-based consensus map of common wheat (Triticum aestivum L.) by integrating genetic maps from four recombinant inbred line populations. The populations were each genotyped using the wheat 90K Infinium iSelect SNP assay. A total of 29,692 SNP markers were mapped on 21 linkage groups corresponding to 21 hexaploid wheat chromosomes, covering 2,906.86 cM, with an overall marker density of 10.21 markers/cM. Compared with the previous maps based on the wheat 90K SNP chip detected 22,736 (76.6%) of the SNPs with consistent chromosomal locations, whereas 1,974 (6.7%) showed different chromosomal locations, and 4,982 (16.8%) were newly mapped. Alignment of the present consensus map and the wheat expressed sequence tags (ESTs) Chromosome Bin Map enabled assignment of 1,221 SNP markers to specific chromosome bins and 819 ESTs were integrated into the consensus map. The marker orders of the consensus map were validated based on physical positions on the wheat genome with Spearman rank correlation coefficients ranging from 0.69 (4D) to 0.97 (1A, 4B, 5B, and 6A), and were also confirmed by comparison with genetic position on the previously 40K SNP consensus map with Spearman rank correlation coefficients ranging from 0.84 (6D) to 0.99 (6A). Chromosomal rearrangements reported previously were confirmed in the present consensus map and new putative rearrangements were identified. In addition, an integrated consensus map was developed through the combination of five published maps with ours, containing 52,607 molecular markers. The consensus map described here provided a high-density SNP marker map and a reliable order of SNPs, representing a step forward in mapping and validation of chromosomal locations of SNPs on the wheat 90K array. Moreover, it can be used as a reference for quantitative trait loci (QTL) mapping to facilitate exploitation of genes and QTL in wheat breeding.

QTL analysis of soft scald in two apple populations

The apple (Malus×domestica Borkh.) is one of the world's most widely grown and valuable fruit crops. With demand for apples year round, storability has emerged as an important consideration for apple breeding programs. Soft scald is a cold storage-related disorder that results in sunken, darkened tissue on the fruit surface. Apple breeders are keen to generate new cultivars that do not suffer from soft scald and can thus be marketed year round. Traditional breeding approaches are protracted and labor intensive, and therefore marker-assisted selection (MAS) is a valuable tool for breeders. To advance MAS for storage disorders in apple, we used genotyping-by-sequencing (GBS) to generate high-density genetic maps in two F1 apple populations, which were then used for quantitative trait locus (QTL) mapping of soft scald. In total, 900 million DNA sequence reads were generated, but after several data filtering steps, only 2% of reads were ultimately used to create two genetic maps that included 1918 and 2818 single-nucleotide polymorphisms. Two QTL associated with soft scald were identified in one of the bi-parental populations originating from parent 11W-12-11, an advanced breeding line. This study demonstrates the utility of next-generation DNA sequencing technologies for QTL mapping in F1 populations, and provides a basis for the advancement of MAS to improve storability of apples.

A high-density, multi-parental SNP genetic map on apple validates a new mapping approach for outcrossing species

Quantitative trait loci (QTL) mapping approaches rely on the correct ordering of molecular markers along the chromosomes, which can be obtained from genetic linkage maps or a reference genome sequence. For apple (Malus domestica Borkh), the genome sequence v1 and v2 could not meet this need; therefore, a novel approach was devised to develop a dense genetic linkage map, providing the most reliable marker-loci order for the highest possible number of markers. The approach was based on four strategies: (i) the use of multiple full-sib families, (ii) the reduction of missing information through the use of HaploBlocks and alternative calling procedures for single-nucleotide polymorphism (SNP) markers, (iii) the construction of a single backcross-type data set including all families, and (iv) a two-step map generation procedure based on the sequential inclusion of markers. The map comprises 15 417 SNP markers, clustered in 3 K HaploBlock markers spanning 1 267 cM, with an average distance between adjacent markers of 0.37 cM and a maximum distance of 3.29 cM. Moreover, chromosome 5 was oriented according to its homoeologous chromosome 10. This map was useful to improve the apple genome sequence, design the Axiom Apple 480 K SNP array and perform multifamily-based QTL studies. Its collinearity with the genome sequences v1 and v3 are reported. To our knowledge, this is the shortest published SNP map in apple, while including the largest number of markers, families and individuals. This result validates our methodology, proving its value for the construction of integrated linkage maps for any outbreeding species.

A dense SNP genetic map constructed using restriction site-associated DNA sequencing enables detection of QTLs controlling apple fruit quality

Background

Genetic map based quantitative trait locus (QTL) analysis is an important method for studying important horticultural traits in apple. To facilitate molecular breeding studies of fruit quality traits in apple, we aim to construct a high density map which was efficient for QTL mapping and possible to search for candidate genes directly in mapped QTLs regions.

Methods

A total of 1733 F1 seedlings derived from ‘Jonathan’ × ‘Golden Delicious’ was used for the map constructionand QTL analysis. The SNP markers were developed by restriction site-associated DNA sequencing (RADseq). Phenotyping data of fruit quality traits were calculated in 2008-2011. Once QTLs were mapped, candidate genes were searched for in the corresponding regions of the apple genome sequence underlying the QTLs. Then some of the candidate genes were validated using real-time PCR.

Results

A high-density genetic map with 3441 SNP markers from 297 individuals was generated. Of the 3441 markers, 2017 were mapped to ‘Jonathan’ with a length of 1343.4 cM and the average distance between markers was 0.67 cM, 1932 were mapped to ‘Golden Delicious’ with a length of 1516.0 cM and the average distance between markers was 0.78 cM. Twelve significant QTLs linked to the control of fruit weight, fruit firmness, sugar content and fruit acidity were mapped to seven linkage groups. Based on gene annotation, 80, 64 and 17 genes related to fruit weight, fruit firmness and fruit acidity, respectively, were analyzed.Among the 17 candidate genes associated with control of fruit acidity, changes in the expression of MDP0000582174 (MdMYB4) were in agreement with the pattern of changes in malic acid content in apple during ripening, and the relative expression of MDP0000239624 (MdME) was significantly correlated withfruit acidity.

Conclusions

We demonstrated the construction of a dense SNP genetic map in apple using next generation sequencing and that the increased resolution enabled the detection of narrow interval QTLs linked to the three fruit quality traits assessed. The candidate genes MDP0000582174 and MDP0000239624 were found to be related to fruit acidity regulation. We conclude that application of RADseq for genetic map construction improved the precision of QTL detection and should be utilized in future studies on the regulatory mechanisms of important fruit traits in apple.

Electronic supplementary material

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

High-density interspecific genetic maps of kiwifruit and the identification of sex-specific markers

Kiwifruit (Actinidia chinensis Planchon) is an important specialty fruit crop that suffers from narrow genetic diversity stemming from recent global commercialization and limited cultivar improvement. Here, we present high-density RAD-seq-based genetic maps using an interspecific F₁ cross between Actinidia rufa ‘MT570001’ and A. chinensis ‘Guihai No4’. The A. rufa (maternal) map consists of 2,426 single-nucleotide polymorphism (SNP) markers with a total length of 2,651 cM in 29 linkage groups (LGs) corresponding to the 29 chromosomes. The A. chinensis (paternal) map consists of 4,214 SNP markers over 3,142 cM in 29 LGs. Using these maps, we were able to anchor an additional 440 scaffolds from the kiwifruit draft genome assembly. Kiwifruit is functionally dioecious, which presents unique challenges for breeding and production. Three sex-specific simple sequence repeats (SSR) markers can be used to accurately sex type male and female kiwifruit in breeding programmes. The sex-determination region (SDR) in kiwifruit was narrowed to a 1-Mb subtelomeric region on chromosome 25. Localizing the SDR will expedite the discovery of genes controlling carpel abortion in males and pollen sterility in females.

Multiple QTLs Linked to Agro-Morphological and Physiological Traits Related to Drought Tolerance in Potato

Dissection of the genetic architecture of adaptation and abiotic stress-related traits is highly desirable for developing drought-tolerant potatoes and enhancing the resilience of existing cultivars, particularly as agricultural production in rain-fed areas may be reduced by up to 50 % by 2020. The "DMDD" potato progeny was developed at International Potato Center (CIP) by crossing the sequenced double monoploid line DM and a diploid cultivar of the Solanum tuberosum diploid Andigenum Goniocalyx group. Recently, a high-density integrated genetic map based on single nucleotide polymorphism (SNP), diversity array technology (DArT), simple sequence repeats (SSRs), and amplified fragment length polymorphism (AFLP) markers was also made available for this population. Two trials were conducted, in greenhouse and field, for drought tolerance with two treatments each, well-watered and terminal drought, in which watering was suspended 60 days after planting. The DMDD population was evaluated for agro-morphological and physiological traits before and after initiation of stress, at multiple time points. Two dense parental genetic maps were constructed using published genotypic data, and quantitative trait locus (QTL) analysis identified 45 genomic regions associated with nine traits in well-watered and terminal drought treatments and 26 potentially associated with drought stress. In this study, the strong influence of environmental factors besides water shortage on the expression of traits and QTLs reflects the multigenic control of traits related to drought tolerance. This is the first study to our knowledge in potato identifying QTLs for drought-related traits in field and greenhouse trials, giving new insights into genetic architecture of drought-related traits. Many of the QTLs identified have the potential to be used in potato breeding programs for enhanced drought tolerance.

High-density genetic linkage map construction and identification of fruit-related QTLs in pear using SNP and SSR markers

Pear (Pyrus spp) is an important fruit crop, grown in all temperate regions of the world, with global production ranked after grape and apples among deciduous tree crops. A high-density linkage map is a valuable tool for fine mapping quantitative trait loci (QTL) and map-based gene cloning. In this study, we firstly constructed a high-density linkage map of pear using SNPs integrated with SSRs, developed by the rapid and robust technology of restriction-associated DNA sequencing (RADseq). The linkage map consists of 3143 SNP markers and 98 SSRs, 3241 markers in total, spanning 2243.4 cM, with an average marker distance of 0.70 cM. Anchoring SSRs were able to anchor seventeen linkage groups to their corresponding chromosomes. Based on this high-density integrated pear linkage map and two years of fruit phenotyping, a total of 32 potential QTLs for 11 traits, including length of pedicel (LFP), single fruit weight (SFW), soluble solid content (SSC), transverse diameter (TD), vertical diameter (VD), calyx status (CS), flesh colour (FC), juice content (JC), number of seeds (NS), skin colour (SC), and skin smooth (SS), were identified and positioned on the genetic map. Among them, some important fruit-related traits have for the first time been identified, such as calyx status, length of pedicel, and flesh colour, and reliable localization of QTLs were verified repeatable. This high-density linkage map of pear is a worthy reference for mapping important fruit traits, QTL identification, and comparison and combination of different genetic maps.

A consensus framework map of durum wheat (Triticum durum Desf.) suitable for linkage disequilibrium analysis and genome-wide association mapping

BACKGROUND

Durum wheat (Triticum durum Desf.) is a tetraploid cereal grown in the medium to low-precipitation areas of the Mediterranean Basin, North America and South-West Asia. Genomics applications in durum wheat have the potential to boost exploitation of genetic resources and to advance understanding of the genetics of important complex traits (e.g. resilience to environmental and biotic stresses). A dense and accurate consensus map specific for T. durum will greatly facilitate genetic mapping, functional genomics and marker-assisted improvement.

RESULTS

High quality genotypic data from six core recombinant inbred line populations were used to obtain a consensus framework map of 598 simple sequence repeats (SSR) and Diversity Array Technology® (DArT) anchor markers (common across populations). Interpolation of unique markers from 14 maps allowed us to position a total of 2,575 markers in a consensus map of 2,463 cM. The T. durum A and B genomes were covered in their near totality based on the reference SSR hexaploid wheat map. The consensus locus order compared to those of the single component maps showed good correspondence, (average Spearman's rank correlation rho ρ value of 0.96). Differences in marker order and local recombination rate were observed between the durum and hexaploid wheat consensus maps. The consensus map was used to carry out a whole-genome search for genetic differentiation signatures and association to heading date in a panel of 183 accessions adapted to the Mediterranean areas. Linkage disequilibrium was found to decay below the r2 threshold=0.3 within 2.20 cM, on average. Strong molecular differentiations among sub-populations were mapped to 87 chromosome regions. A genome-wide association scan for heading date from 27 field trials in the Mediterranean Basin and in Mexico yielded 50 chromosome regions with evidences of association in multiple environments.

CONCLUSIONS

The consensus map presented here was used as a reference for genetic diversity and mapping analyses in T. durum, providing nearly complete genome coverage and even marker density. Markers previously mapped in hexaploid wheat constitute a strong link between the two species. The consensus map provides the basis for high-density single nucleotide polymorphic (SNP) marker implementation in durum wheat.

Fast and cost-effective genetic mapping in apple using next-generation sequencing

Next-generation DNA sequencing (NGS) produces vast amounts of DNA sequence data, but it is not specifically designed to generate data suitable for genetic mapping. Recently developed DNA library preparation methods for NGS have helped solve this problem, however, by combining the use of reduced representation libraries with DNA sample barcoding to generate genome-wide genotype data from a common set of genetic markers across a large number of samples. Here we use such a method, called genotyping-by-sequencing (GBS), to produce a data set for genetic mapping in an F1 population of apples (Malus × domestica) segregating for skin color. We show that GBS produces a relatively large, but extremely sparse, genotype matrix: over 270,000 SNPs were discovered but most SNPs have too much missing data across samples to be useful for genetic mapping. After filtering for genotype quality and missing data, only 6% of the 85 million DNA sequence reads contributed to useful genotype calls. Despite this limitation, using existing software and a set of simple heuristics, we generated a final genotype matrix containing 3967 SNPs from 89 DNA samples from a single lane of Illumina HiSeq and used it to create a saturated genetic linkage map and to identify a known QTL underlying apple skin color. We therefore demonstrate that GBS is a cost-effective method for generating genome-wide SNP data suitable for genetic mapping in a highly diverse and heterozygous agricultural species. We anticipate future improvements to the GBS analysis pipeline presented here that will enhance the utility of next-generation DNA sequence data for the purposes of genetic mapping across diverse species.

A high-density SNP Map of sunflower derived from RAD-sequencing facilitating fine-mapping of the rust resistance gene R12

A high-resolution genetic map of sunflower was constructed by integrating SNP data from three F₂ mapping populations (HA 89/RHA 464, B-line/RHA 464, and CR 29/RHA 468). The consensus map spanned a total length of 1443.84 cM, and consisted of 5,019 SNP markers derived from RAD tag sequencing and 118 publicly available SSR markers distributed in 17 linkage groups, corresponding to the haploid chromosome number of sunflower. The maximum interval between markers in the consensus map is 12.37 cM and the average distance is 0.28 cM between adjacent markers. Despite a few short-distance inversions in marker order, the consensus map showed high levels of collinearity among individual maps with an average Spearman's rank correlation coefficient of 0.972 across the genome. The order of the SSR markers on the consensus map was also in agreement with the order of the individual map and with previously published sunflower maps. Three individual and one consensus maps revealed the uneven distribution of markers across the genome. Additionally, we performed fine mapping and marker validation of the rust resistance gene R₁₂, providing closely linked SNP markers for marker-assisted selection of this gene in sunflower breeding programs. This high resolution consensus map will serve as a valuable tool to the sunflower community for studying marker-trait association of important agronomic traits, marker assisted breeding, map-based gene cloning, and comparative mapping.

YouGenMap: a web platform for dynamic multi-comparative mapping and visualization of genetic maps

Comparative genetic maps are used in examination of genome organization, detection of conserved gene order, and exploration of marker order variations. YouGenMap is an open-source web tool that offers dynamic comparative mapping capability of users' own genetic mapping between 2 or more map sets. Users' genetic map data and optional gene annotations are uploaded, either publically or privately, as long as they follow our template which is available in several standard file formats. Data is parsed and loaded into MySQL relational database to be displayed and compared against users' genetic maps or other public data available on YouGenMap. With the highly interactive GUIs, all public data on YouGenMap are maps available for visualization, comparison, search, filtration and download. YouGenMap web tool is available on the website (http://conifergdb.miamioh.edu/yougenmap) with the source-code repository at (http://sourceforge.net/projects/yougenmap/?source=directory).

A reference consensus genetic map for molecular markers and economically important traits in faba bean (Vicia faba L.)

BACKGROUND

Faba bean (Vicia faba L.) is among the earliest domesticated crops from the Near East. Today this legume is a key protein feed and food worldwide and continues to serve an important role in culinary traditions throughout Middle East, Mediterranean region, China and Ethiopia. Adapted to a wide range of soil types, the main faba bean breeding objectives are to improve yield, resistance to biotic and abiotic stresses, seed quality and other agronomic traits. Genomic approaches aimed at enhancing faba bean breeding programs require high-quality genetic linkage maps to facilitate quantitative trait locus analysis and gene tagging for use in a marker-assisted selection. The objective of this study was to construct a reference consensus map in faba bean by joining the information from the most relevant maps reported so far in this crop.

RESULTS

A combination of two approaches, increasing the number of anchor loci in diverse mapping populations and joining the corresponding genetic maps, was used to develop a reference consensus map in faba bean. The map was constructed from three main recombinant inbreed populations derived from four parental lines, incorporates 729 markers and is based on 69 common loci. It spans 4,602 cM with a range from 323 to 1041 loci in six main linkage groups or chromosomes, and an average marker density of one locus every 6 cM. Locus order is generally well maintained between the consensus map and the individual maps.

CONCLUSION

We have constructed a reliable and fairly dense consensus genetic linkage map that will serve as a basis for genomic approaches in faba bean research and breeding. The core map contains a larger number of markers than any previous individual map, covers existing gaps and achieves a wider coverage of the large faba bean genome as a whole. This tool can be used as a reference resource for studies in different genetic backgrounds, and provides a framework for transferring genetic information when using different marker technologies. Combined with syntenic approaches, the consensus map will increase marker density in selected genomic regions and will be useful for future faba bean molecular breeding applications.

Identification of genetic loci associated with fire blight resistance in Malus through combined use of QTL and association mapping

Fire blight, incited by the enterobacterium Erwinia amylovora, is a destructive disease of Rosaceae, particularly of apples and pears. There are reports on the molecular mechanisms underlying E. amylovora pathogenesis and how the host activates its resistance mechanism. The host's resistance mechanism is quantitatively controlled, although some major genes might also be involved. Thus far, quantitative trait loci (QTL) mapping and differential expression studies have been used to elucidate those genes and/or genomic regions underlying quantitative resistance present in the apple genome. In this study, an effort is undertaken to dissect the genetic basis of fire blight resistance in apple using both QTL and genome-wide association mapping. On the basis of an F1 pedigree of 'Coop 16' × 'Coop 17' and a genome-wide association study (GWAS) mapping population of Malus accessions (species, old and new cultivars and selections), new QTLs and associations have been identified. A total of three QTLs for resistance to fire blight, with above 95% significant logarithm of odds threshold value of 2.5, have been identified on linkage groups (LGs) 02, 06, and 15 of the apple genome with phenotypic variation explained values of 14.7, 20.1 and 17.4, respectively. Although elevated P-values with signals for marker-trait associations are observed for some LGs, these are not found to be significant. However, a total of 34 significant associations, with P-values ≥0.02, have been detected including 8 for lesion length at 7 days following inoculation (PL1), 14 for lesion length at 14 days following inoculation (PL2), and 12 for shoot length.