Construction of High Density Sweet Cherry (Prunus avium L.) Linkage Maps Using Microsatellite Markers and SNPs Detected by Genotyping-by-Sequencing (GBS)

Linkage maps are valuable tools in genetic and genomic studies. For sweet cherry, linkage maps have been constructed using mainly microsatellite markers (SSRs) and, recently, using single nucleotide polymorphism markers (SNPs) from a cherry 6K SNP array. Genotyping-by-sequencing (GBS), a new methodology based on high-throughput sequencing, holds great promise for identification of high number of SNPs and construction of high density linkage maps. In this study, GBS was used to identify SNPs from an intra-specific sweet cherry cross. A total of 8,476 high quality SNPs were selected for mapping. The physical position for each SNP was determined using the peach genome, Peach v1.0, as reference, and a homogeneous distribution of markers along the eight peach scaffolds was obtained. On average, 65.6% of the SNPs were present in genic regions and 49.8% were located in exonic regions. In addition to the SNPs, a group of SSRs was also used for construction of linkage maps. Parental and consensus high density maps were constructed by genotyping 166 siblings from a ‘Rainier’ x ‘Rivedel’ (Ra x Ri) cross. Using Ra x Ri population, 462, 489 and 985 markers were mapped into eight linkage groups in ‘Rainier’, ‘Rivedel’ and the Ra x Ri map, respectively, with 80% of mapped SNPs located in genic regions. Obtained maps spanned 549.5, 582.6 and 731.3 cM for ‘Rainier’, ‘Rivedel’ and consensus maps, respectively, with an average distance of 1.2 cM between adjacent markers for both ‘Rainier’ and ‘Rivedel’ maps and of 0.7 cM for Ra x Ri map. High synteny and co-linearity was observed between obtained maps and with Peach v1.0. These new high density linkage maps provide valuable information on the sweet cherry genome, and serve as the basis for identification of QTLs and genes relevant for the breeding of the species.     

High-throughput genotyping in citrus accessions using an SNP genotyping array

We developed a 384 multiplexed SNP array, named CitSGA-1, for the genotyping of Citrus cultivars, and evaluated the performance and reliability of the genotyping. SNPs were surveyed by direct sequence comparison of the sequence tagged site (STS) fragment amplified from genomic DNA of cultivars representing the genetic diversity of citrus breeding in Japan. Among 1497 SNPs candidates, 384 SNPs for a high-throughput genotyping array were selected based on physical parameters of Illumina’s bead array criteria. The assay using CitSGA-1 was applied to a hybrid population of 88 progeny and 103 citrus accessions for breeding in Japan, which resulted in 73,726 SNP calls. A total of 351 SNPs (91 %) could call different genotypes among the DNA samples, resulting in a success rate for the assay comparable to previously reported rates for other plant species. To confirm the reliability of SNP genotype calls, parentage analysis was applied, and it indicated that the number of reliable SNPs and corresponding STSs were 276 and 213, respectively. The multiplexed SNP genotyping array reported here will be useful for the efficient construction of linkage map, for the detection of markers for marker-assisted breeding, and for the identification of cultivars.     

An EST-derived SNP and SSR genetic linkage map of cassava (Manihot esculenta Crantz)

Cassava (Manihot esculenta Crantz) is one of the most important food security crops in the tropics and increasingly being adopted for agro-industrial processing. Genetic improvement of cassava can be enhanced through marker-assisted breeding. For this, appropriate genomic tools are required to dissect the genetic architecture of economically important traits. Here, a genome-wide SNP-based genetic map of cassava anchored in SSRs is presented. An outbreeder full-sib (F1) family was genotyped on two independent SNP assay platforms: an array of 1,536 SNPs on Illumina's GoldenGate platform was used to genotype a first batch of 60 F1. Of the 1,358 successfully converted SNPs, 600 which were polymorphic in at least one of the parents and was subsequently converted to KBiosciences' KASPar assay platform for genotyping 70 additional F1. High-precision genotyping of 163 informative SSRs using capillary electrophoresis was also carried out. Linkage analysis resulted in a final linkage map of 1,837 centi-Morgans (cM) containing 568 markers (434 SNPs and 134 SSRs) distributed across 19 linkage groups. The average distance between adjacent markers was 3.4?cM. About 94.2% of the mapped SNPs and SSRs have also been localized on scaffolds of version 4.1 assembly of the cassava draft genome sequence. This more saturated genetic linkage map of cassava that combines SSR and SNP markers should find several applications in the improvement of cassava including aligning scaffolds of the cassava genome sequence, genetic analyses of important agro-morphological traits, studying the linkage disequilibrium landscape and comparative genomics.     

High-throughput single nucleotide polymorphism (SNP) identification and mapping in the sesame (<Emphasis Type="Italic">Sesamum indicum</Emphasis> L.) genome with genotyping by sequencing (GBS) analysis

Sesame (Sesamum indicum L. syn. Sesamum orientale L.) is considered to be the first oil seed crop known to man. Despite its versatile use as an oil seed and a leafy vegetable, sesame is a neglected crop and has not been a subject of molecular genetic research until the last decade. There is thus limited knowledge regarding genome-specific molecular markers that are indispensible for germplasm enhancement, gene identification, and marker-assisted breeding in sesame. In this study, we employed a genotyping by sequencing (GBS) approach to a sesame recombinant inbred line (RIL) population for high-throughput single nucleotide polymorphism (SNP) identification and genotyping. A total of 15,521 SNPs were identified with 14,786 SNPs (95.26 %) located along sesame genome assembly pseudomolecules. By incorporating sesame-specific simple sequence repeat (SSR) markers developed in our previous work, 230.73 megabases (99 %) of sequence from the genome assembly were saturated with markers. This large number of markers will be available for sesame geneticists as a resource for candidate polymorphisms located along the physical chromosomes of sesame. Defining SNP loci in genome assembly sequences provides the flexibility to utilize any genotyping strategy to survey any sesame population. SNPs selected through a high stringency filtering protocol (770 SNPs) for improved map accuracy were used in conjunction with SSR markers (50 SSRs) in linkage analysis, resulting in 13 linkage groups that encompass a total genetic distance of 914 cM with 432 markers (420 SNPs, 12 SSRs). The genetic linkage map constitutes the basis for future work that will involve quantitative trait locus (QTL) analyses of metabolic and agronomic traits in the segregating RIL population.     

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.     

Development and preliminary evaluation of a 90?K Axiom? SNP array for the allo-octoploid cultivated strawberry Fragaria × ananassa

Background

A high-throughput genotyping platform is needed to enable marker-assisted breeding in the allo-octoploid cultivated strawberry Fragaria × ananassa. Short-read sequences from one diploid and 19 octoploid accessions were aligned to the diploid Fragaria vesca ‘Hawaii 4’ reference genome to identify single nucleotide polymorphisms (SNPs) and indels for incorporation into a 90 K Affymetrix® Axiom® array. We report the development and preliminary evaluation of this array.

Results

About 36 million sequence variants were identified in a 19 member, octoploid germplasm panel. Strategies and filtering pipelines were developed to identify and incorporate markers of several types: di-allelic SNPs (66.6%), multi-allelic SNPs (1.8%), indels (10.1%), and ploidy-reducing “haploSNPs” (11.7%). The remaining SNPs included those discovered in the diploid progenitor F. iinumae (3.9%), and speculative “codon-based” SNPs (5.9%). In genotyping 306 octoploid accessions, SNPs were assigned to six classes with Affymetrix’s “SNPolisher” R package. The highest quality classes, PolyHigh Resolution (PHR), No Minor Homozygote (NMH), and Off-Target Variant (OTV) comprised 25%, 38%, and 1% of array markers, respectively. These markers were suitable for genetic studies as demonstrated in the full-sib family ‘Holiday’ × ‘Korona’ with the generation of a genetic linkage map consisting of 6,594 PHR SNPs evenly distributed across 28 chromosomes with an average density of approximately one marker per 0.5 cM, thus exceeding our goal of one marker per cM.

Conclusions

The Affymetrix IStraw90 Axiom array is the first high-throughput genotyping platform for cultivated strawberry and is commercially available to the worldwide scientific community. The array’s high success rate is likely driven by the presence of naturally occurring variation in ploidy level within the nominally octoploid genome, and by effectiveness of the employed array design and ploidy-reducing strategies. This array enables genetic analyses including generation of high-density linkage maps, identification of quantitative trait loci for economically important traits, and genome-wide association studies, thus providing a basis for marker-assisted breeding in this high value crop.

Electronic supplementary material

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

A multi-array multi-SNP genotyping algorithm for Affymetrix SNP microarrays

MOTIVATION: Modern strategies for mapping disease loci require efficient genotyping of a large number of known polymorphic sites in the genome. The sensitive and high-throughput nature of hybridization-based DNA microarray technology provides an ideal platform for such an application by interrogating up to hundreds of thousands of single nucleotide polymorphisms (SNPs) in a single assay. Similar to the development of expression arrays, these genotyping arrays pose many data analytic challenges that are often platform specific. Affymetrix SNP arrays, e.g. use multiple sets of short oligonucleotide probes for each known SNP, and require effective statistical methods to combine these probe intensities in order to generate reliable and accurate genotype calls. RESULTS: We developed an integrated multi-SNP, multi-array genotype calling algorithm for Affymetrix SNP arrays, MAMS, that combines single-array multi-SNP (SAMS) and multi-array, single-SNP (MASS) calls to improve the accuracy of genotype calls, without the need for training data or computation-intensive normalization procedures as in other multi-array methods. The algorithm uses resampling techniques and model-based clustering to derive single array based genotype calls, which are subsequently refined by competitive genotype calls based on (MASS) clustering. The resampling scheme caps computation for single-array analysis and hence is readily scalable, important in view of expanding numbers of SNPs per array. The MASS update is designed to improve calls for atypical SNPs, harboring allele-imbalanced binding affinities, that are difficult to genotype without information from other arrays. Using a publicly available data set of HapMap samples from Affymetrix, and independent calls by alternative genotyping methods from the HapMap project, we show that our approach performs competitively to existing methods. AVAILABILITY: R functions are available upon request from the authors.     

Genotyping by Sequencing for SNP-Based Linkage Map Construction and QTL Analysis of Chilling Requirement and Bloom Date in Peach [Prunus persica (L.) Batsch]

Low-cost, high throughput genotyping methods are crucial to marker discovery and marker-assisted breeding efforts, but have not been available for many ‘specialty crops’ such as fruit and nut trees. Here we apply the Genotyping-By-Sequencing (GBS) method developed for cereals to the discovery of single nucleotide polymorphisms (SNPs) in a peach F₂ mapping population. Peach is a genetic and genomic model within the Rosaceae and will provide a template for the use of this method with other members of this family. Our F₂ mapping population of 57 genotypes segregates for bloom time (BD) and chilling requirement (CR) and we have extensively phenotyped this population. The population derives from a selfed F₁ progeny of a cross between ‘Hakuho’ (high CR) and ‘UFGold’ (low CR). We were able to successfully employ GBS and the TASSEL GBS pipeline without modification of the original methodology using the ApeKI restriction enzyme and multiplexing at an equivalent of 96 samples per Illumina HiSeq 2000 lane. We obtained hundreds of SNP markers which were then used to construct a genetic linkage map and identify quantitative trait loci (QTL) for BD and CR.     

Evaluation of SNP Data from the Malus Infinium Array Identifies Challenges for Genetic Analysis of Complex Genomes of Polyploid Origin

High throughput arrays for the simultaneous genotyping of thousands of single-nucleotide polymorphisms (SNPs) have made the rapid genetic characterisation of plant genomes and the development of saturated linkage maps a realistic prospect for many plant species of agronomic importance. However, the correct calling of SNP genotypes in divergent polyploid genomes using array technology can be problematic due to paralogy, and to divergence in probe sequences causing changes in probe binding efficiencies. An Illumina Infinium II whole-genome genotyping array was recently developed for the cultivated apple and used to develop a molecular linkage map for an apple rootstock progeny (M432), but a large proportion of segregating SNPs were not mapped in the progeny, due to unexpected genotype clustering patterns. To investigate the causes of this unexpected clustering we performed BLAST analysis of all probe sequences against the ‘Golden Delicious’ genome sequence and discovered evidence for paralogous annealing sites and probe sequence divergence for a high proportion of probes contained on the array. Following visual re-evaluation of the genotyping data generated for 8,788 SNPs for the M432 progeny using the array, we manually re-scored genotypes at 818 loci and mapped a further 797 markers to the M432 linkage map. The newly mapped markers included the majority of those that could not be mapped previously, as well as loci that were previously scored as monomorphic, but which segregated due to divergence leading to heterozygosity in probe annealing sites. An evaluation of the 8,788 probes in a diverse collection of Malus germplasm showed that more than half the probes returned genotype clustering patterns that were difficult or impossible to interpret reliably, highlighting implications for the use of the array in genome-wide association studies.     

A Multi-Population Consensus Genetic Map Reveals Inconsistent Marker Order among Maps Likely Attributed to Structural Variations in the Apple Genome

Genetic maps serve as frameworks for determining the genetic architecture of quantitative traits, assessing structure of a genome, as well as aid in pursuing association mapping and comparative genetic studies. In this study, a dense genetic map was constructed using a high-throughput 1,536 EST-derived SNP GoldenGate genotyping platform and a global consensus map established by combining the new genetic map with four existing reliable genetic maps of apple. The consensus map identified markers with both major and minor conflicts in positioning across all five maps. These major inconsistencies among marker positions were attributed either to structural variations within the apple genome, or among mapping populations, or genotyping technical errors. These also highlighted problems in assembly and anchorage of the reference draft apple genome sequence in regions with known segmental duplications. Markers common across all five apple genetic maps resulted in successful positioning of 2875 markers, consisting of 2033 SNPs and 843 SSRs as well as other specific markers, on the global consensus map. These markers were distributed across all 17 linkage groups, with an average of 169±33 marker per linkage group and with an average distance of 0.70±0.14 cM between markers. The total length of the consensus map was 1991.38 cM with an average length of 117.14±24.43 cM per linkage group. A total of 569 SNPs were mapped onto the genetic map, consisting of 140 recombinant individuals, from our recently developed apple Oligonucleotide pool assays (OPA). The new functional SNPs, along with the dense consensus genetic map, will be useful for high resolution QTL mapping of important traits in apple and for pursuing comparative genetic studies in Rosaceae.     

Transcriptome-based single nucleotide polymorphism markers for genome mapping in Japanese pear (Pyrus pyrifolia Nakai)

The development of single nucleotide polymorphism (SNP) markers in Japanese pear (Pyrus pyrifolia Nakai) offers the opportunity to use DNA markers for marker-assisted selection in breeding programs because of their high abundance, codominant inheritance, and potential for automated high-throughput analysis. We developed a 1,536-SNP bead array without a reference genome sequence from more than 44,000 base changes on the basis of a large-scale expressed sequence tag (EST) analysis combined with 454 genome sequencing data of Japanese pear ‘Housui’. Among the 1,536 SNPs on the array, 756 SNPs were genotyped, and 609 SNP loci were mapped to linkage groups on a genetic linkage map of ‘Housui’, based on progeny of an interspecific cross between European pear (Pyrus communis L.) ‘Bartlett’ and ‘Housui’. The newly constructed genetic linkage map consists of 951 loci, comprising 609 new SNPs, 110 pear genomic simple sequence repeats (SSRs), 25 pear EST–SSRs, 127 apple SSRs, 61 pear SNPs identified by the “potential intron polymorphism” method, and 19 other loci. The map covers 22 linkage groups spanning 1341.9 cM with an average distance of 1.41 cM between markers and is anchored to reference genetic linkage maps of European pears and apples. A total of 514 contigs containing mapped SNP loci showed significant similarity to known proteins by functional annotation analysis.     

Development of a molecular linkage map of pearl millet integrating DArT and SSR markers

Pearl millet is an important component of food security in the semi-arid tropics and is assuming greater importance in the context of changing climate and increasing demand for highly nutritious food and feed. Molecular tools have been developed and applied for pearl millet on a limited scale. However, the existing tool kit needs to be strengthened further for its routine use in applied breeding programs. Here, we report enrichment of the pearl millet molecular linkage map by exploiting low-cost and high-throughput Diversity Arrays Technology (DArT) markers. Genomic representation from 95 diverse genotypes was used to develop a DArT array with circa 7,000 clones following PstI/BanII complexity reduction. This array was used to genotype a set of 24 diverse pearl millet inbreds and 574 polymorphic DArT markers were identified. The genetic relationships among the inbred lines as revealed by DArT genotyping were in complete agreement with the available pedigree data. Further, a mapping population of 140 F₇ Recombinant Inbred Lines (RILs) from cross H 77/833-2 × PRLT 2/89-33 was genotyped and an improved linkage map was constructed by integrating DArT and SSR marker data. This map contains 321 loci (258 DArTs and 63 SSRs) and spans 1148 cM with an average adjacent-marker interval length of 3.7 cM. The length of individual linkage groups (LGs) ranged from 78 cM (LG 3) to 370 cM (LG 2). This better-saturated map provides improved genome coverage and will be useful for genetic analyses of important quantitative traits. This DArT platform will also permit cost-effective background selection in marker-assisted backcrossing programs as well as facilitate comparative genomics and genome organization studies once DNA sequences of polymorphic DArT clones are available.     

The USDA Barley Core Collection: Genetic Diversity,Population Structure,and Potential for Genome-Wide Association Studies

New sources of genetic diversity must be incorporated into plant breeding programs if they are to continue increasing grain yield and quality, and tolerance to abiotic and biotic stresses. Germplasm collections provide a source of genetic and phenotypic diversity, but characterization of these resources is required to increase their utility for breeding programs. We used a barley SNP iSelect platform with 7,842 SNPs to genotype 2,417 barley accessions sampled from the USDA National Small Grains Collection of 33,176 accessions. Most of the accessions in this core collection are categorized as landraces or cultivars/breeding lines and were obtained from more than 100 countries. Both STRUCTURE and principal component analysis identified five major subpopulations within the core collection, mainly differentiated by geographical origin and spike row number (an inflorescence architecture trait). Different patterns of linkage disequilibrium (LD) were found across the barley genome and many regions of high LD contained traits involved in domestication and breeding selection. The genotype data were used to define ‘mini-core’ sets of accessions capturing the majority of the allelic diversity present in the core collection. These ‘mini-core’ sets can be used for evaluating traits that are difficult or expensive to score. Genome-wide association studies (GWAS) of ‘hull cover’, ‘spike row number’, and ‘heading date’ demonstrate the utility of the core collection for locating genetic factors determining important phenotypes. The GWAS results were referenced to a new barley consensus map containing 5,665 SNPs. Our results demonstrate that GWAS and high-density SNP genotyping are effective tools for plant breeders interested in accessing genetic diversity in large germplasm collections.     

Construction of High-Density Linkage Maps of Populus deltoides × P. simonii Using Restriction-Site Associated DNA Sequencing

Although numerous linkage maps have been constructed in the genus Populus, they are typically sparse and thus have limited applications due to low throughput of traditional molecular markers. Restriction-site associated DNA sequencing (RADSeq) technology allows us to identify a large number of single nucleotide polymorphisms (SNP) across genomes of many individuals in a fast and cost-effective way, and makes it possible to construct high-density genetic linkage maps. We performed RADSeq for 299 progeny and their two parents in an F₁ hybrid population generated by crossing the female Populus deltoides ‘I-69’ and male Populus simonii ‘L3’. A total of 2,545 high quality SNP markers were obtained and two parent-specific linkage maps were constructed. The female genetic map contained 1601 SNPs and 20 linkage groups, spanning 4,249.12 cM of the genome with an average distance of 2.69 cM between adjacent markers, while the male map consisted of 940 SNPs and also 20 linkage groups with a total length of 3,816.24 cM and an average marker interval distance of 4.15 cM. Finally, our analysis revealed that synteny and collinearity are highly conserved between the parental linkage maps and the reference genome of P. trichocarpa. We demonstrated that RAD sequencing is a powerful technique capable of rapidly generating a large number of SNPs for constructing genetic maps in outbred forest trees. The high-quality linkage maps constructed here provided reliable genetic resources to facilitate locating quantitative trait loci (QTLs) that control growth and wood quality traits in the hybrid population.     

Genetic linkage maps of Populus nigra L. including AFLPs, SSRs, SNPs, and sex trait

Black poplar (Populus nigra L.) is a tree of ecological and economic interest. A better knowledge of P. nigra genome is needed for an effective protection and use of its genetic resources. The main objective of this study is the construction of a highly informative genetic map of P. nigra species including genes of adaptive and economic interest. Two genotypes originated from contrasted natural Italian populations were crossed to generate a F₁ mapping pedigree of 165 individuals. Amplification fragment length polymorphism (AFLP), simple sequence repeat (SSR), and single nucleotide polymorphism (SNP) markers were used to genotype 92 F₁ individuals, and the pseudo-test-cross strategy was applied for linkage analysis. The female parent map included 368 markers (274 AFLPs, 91 SSRs, and 3 SNPs) and spanned 2,104 cM with 20 linkage groups, and the male parent map, including 317 markers (205 AFLPs, 106 SSRs, 5 SNPs, and sex trait), spanned 2,453 cM with 23 main linkage groups. The sex, as morphological trait, was mapped on the linkage group XIX of the male parent map. The generated maps are among the most informative in SSRs when compared to the Populus maps published so far and allow a complete alignment with the 19 haploid chromosomes of Populus sequence genome. These genetic maps provide informative tools for a better understanding of P. nigra genome structure and genetic improvement of this ecologically and economically important European tree species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Simple-sequence repeat markers used in merging linkage maps of melon (Cucumis melo L.)

A set of 118 simple sequence repeat (SSR) markers has been developed in melon from two different sources: genomic libraries (gSSR) and expressed sequence-tag (EST) databases (EST-SSR). Forty-nine percent of the markers showed polymorphism between the Piel de Sapo (PS) and PI161375 melon genotypes used as parents for the mapping populations. Similar polymorphism levels were found in gSSR (51.2%) and EST-SSR (45.5%). Two populations, F₂ and a set of double haploid lines (DHLs), developed from the same parent genotypes were used for map construction. Twenty-three SSRs and 79 restriction fragment length polymorphisms (RFLPs), evenly distributed through the melon genome, were used to anchor the maps of both populations. Ten cucumber SSRs, 41 gSSRs, 16 EST-SSR, three single nucleotide polymorphism (SNP) markers, and the Nsv locus were added in the DHL population. The maps developed in the F₂ and DHL populations were co-linear, with similar lengths, except in linkage groups G1, G9, and G10. There was segregation distortion in a higher proportion of markers in the DHL population compared with the F₂, probably caused by selection during the construction of DHLs through in vitro culture. After map merging, a composite genetic map was obtained including 327 transferable markers: 226 RFLPs, 97 SSRs, three SNPs, and the Nsv locus. The map length is 1,021 cM, distributed in 12 linkage groups, and map density is 3.11 cM/marker. SSR markers alone cover nearly 80% of the map length. This map is proposed as a basis for a framework melon map to be merged with other maps and as an anchor point for map comparison between species of the Cucurbitaceae family.Electronic Supplementary Material Supplementary material is available for this article at     

A ddRAD Based Linkage Map of the Cultivated Strawberry,Fragaria xananassa

The cultivated strawberry (Fragaria ×ananassa Duch.) is an allo-octoploid considered difficult to disentangle genetically due to its four relatively similar sub-genomic chromosome sets. This has been alleviated by the recent release of the strawberry IStraw90 whole genome genotyping array. However, array resolution relies on the genotypes used in the array construction and may be of limited general use. SNP detection based on reduced genomic sequencing approaches has the potential of providing better coverage in cases where the studied genotypes are only distantly related from the SNP array’s construction foundation. Here we have used double digest restriction-associated DNA sequencing (ddRAD) to identify SNPs in a 145 seedling F₁ hybrid population raised from the cross between the cultivars Sonata (♀) and Babette (♂). A linkage map containing 907 markers which spanned 1,581.5 cM across 31 linkage groups representing the 28 chromosomes of the species. Comparing the physical span of the SNP markers with the F. vesca genome sequence, the linkage groups resolved covered 79% of the estimated 830 Mb of the F. ×ananassa genome. Here, we have developed the first linkage map for F. ×ananassa using ddRAD and show that this technique and other related techniques are useful tools for linkage map development and downstream genetic studies in the octoploid strawberry.     

Identification of Pyrus Single Nucleotide Polymorphisms (SNPs) and Evaluation for Genetic Mapping in European Pear and Interspecific Pyrus Hybrids

We have used new generation sequencing (NGS) technologies to identify single nucleotide polymorphism (SNP) markers from three European pear (Pyrus communis L.) cultivars and subsequently developed a subset of 1096 pear SNPs into high throughput markers by combining them with the set of 7692 apple SNPs on the IRSC apple Infinium® II 8K array. We then evaluated this apple and pear Infinium® II 9K SNP array for large-scale genotyping in pear across several species, using both pear and apple SNPs. The segregating populations employed for array validation included a segregating population of European pear (‘Old Home’בLouise Bon Jersey’) and four interspecific breeding families derived from Asian (P. pyrifolia Nakai and P. bretschneideri Rehd.) and European pear pedigrees. In total, we mapped 857 polymorphic pear markers to construct the first SNP-based genetic maps for pear, comprising 78% of the total pear SNPs included in the array. In addition, 1031 SNP markers derived from apple (13% of the total apple SNPs included in the array) were polymorphic and were mapped in one or more of the pear populations. These results are the first to demonstrate SNP transferability across the genera Malus and Pyrus. Our construction of high density SNP-based and gene-based genetic maps in pear represents an important step towards the identification of chromosomal regions associated with a range of horticultural characters, such as pest and disease resistance, orchard yield and fruit quality.     

Development of a high-throughput SNP array for sea cucumber (Apostichopus japonicus) and its application in genomic selection with MCP regularized deep neural networks

High-throughput single nucleotide polymorphism (SNP) genotyping assays are powerful tools for genetic studies and genomic breeding applications for many species. Though large numbers of SNPs have been identified in sea cucumber (Apostichopus japonicus), but, as yet, no high-throughput genotyping platform is available for this species. In this study, we designed and developed a high-throughput 24 K SNP genotyping array named HaishenSNP24K for A. japonicus, based on the multi-objective-local optimization (MOLO) algorithm and HD-Marker genotyping method. The SNP array exhibited a relatively high genotyping call rate (> 96%), genotyping accuracy (>95%) and exhibited highly polymorphic in sea cucumber populations. In addition, we also assessed its application in genomic selection (GS). Deep neural networks (DNN) that can capture the complicated interactions of genes have been proposed as a promising tool in GS for SNP-based genomic prediction of complex traits in animal breeding. To overcome the problem of over-fitting when using the HaishenSNP24K array as high-dimensional DNN input, we developed minmax concave penalty (MCP) regularization for sparse deep neural networks (DNN-MCP) that finds an optimal sparse structure of a DNN by minimizing the square error subject to the non-convex penalty MCP on the parameters (weights and biases). Compared to two linear models, namely RR-GBLUP and Bayes B, and the nonlinear model DNN, DNN-MCP has greatly improved the genomic prediction ability for three quantitative traits (e.g., wet weight, dry weight and survival time) in the sea cucumber population. To the best of our knowledge, this is the first work to develop a high-throughput SNP array for A. japonicus and a new model DNN-MCP for genomic prediction of complex traits in GS. The present results provide evidence that supports the HaishenSNP24K array with DNN-MCP will be valuable for genetic studies and molecular breeding in A. japonicus.     

High-resolution, high-throughput SNP mapping in Drosophila melanogaster

Single nucleotide polymorphisms (SNPs) are useful markers for genetic mapping experiments in model organisms. Here we report the establishment of a high-density SNP map and high-throughput genotyping assays for Drosophila melanogaster. Our map comprises 27,367 SNPs in common laboratory Drosophila stocks. These SNPs were clustered within 2,238 amplifiable markers at an average density of 1 marker every 50.3 kb, or 6.3 genes. We have also constructed a set of 62 Drosophila stocks, each of which facilitates the generation of recombinants within a defined genetic interval of 1-2 Mb. For flexible, high-throughput SNP genotyping, we used fluorescent tag-array mini-sequencing (TAMS) assays. We designed and validated TAMS assays for 293 SNPs at an average resolution of 391.3 kb, and demonstrated the utility of these tools by rapidly mapping 14 mutations that disrupt embryonic muscle patterning. These resources enable high-resolution high-throughput genetic mapping in Drosophila.