A genetic and cytogenetic map for the duck (Anas platyrhynchos)

A genetic linkage map for the duck (Anas platyrhynchos) was developed within a cross between two extreme Peking duck lines by linkage analysis of 155 polymorphic microsatellite markers, including 84 novel markers reported in this study. A total of 115 microsatellite markers were placed into 19 linkage groups. The sex-averaged map spans 1353.3 cM, with an average interval distance of 15.04 cM. The male map covers 1415 cM, whereas the female map covers only 1387.6 cM. All of the flanking sequences of the 155 polymorphic loci--44 monomorphic loci and a further 41 reported microsatellite loci for duck--were blasted against the chicken genomic sequence, and corresponding orthologs were found for 49. To integrate the genetic and cytogenetic map of the duck genome, 28 BAC clones were screened from a chicken BAC library using the specific PCR primers and localized to duck chromosomes by FISH, respectively. Of 28 BAC clones, 24 were detected definitely on duck chromosomes. Thus, 11 of 19 linkage groups were localized to 10 duck chromosomes. This genetic and cytogenetic map will be helpful for the mapping QTL in duck for breeding applications and for conducting genomic comparisons between chicken and duck.     

An extended anchored linkage map and virtual mapping for the American mink genome based on homology to human and dog

In this report we present an extended linkage map of the American mink (Neovison vison) consisting of 157 microsatellite markers and comprising at least one linkage group for each of the autosomes. Each linkage group has been assigned to a chromosome and oriented by fluorescence in situ hybridization (FISH) and/or by means of human/dog/mink comparative homology. The average interval between markers is 8.5 cM and the linkage groups collectively span 1340 cM. In addition, 217 and 275 mink microsatellites have been placed on human and dog genomes, respectively. In conjunction with the existing comparative human/dog/mink data, these assignments represent useful virtual maps for the American mink genome. Comparison of the current human/dog assembled sequential map with the existing Zoo-FISH-based human/dog/mink maps helped to refine the human/dog/mink comparative map. Furthermore, comparison of the human and dog genome assemblies revealed a number of large synteny blocks, some of which are corroborated by data from the mink linkage map.     

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

Background

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

Results

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. ‘SNP_only’ markers accounted for 89.25 % of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9 %, and their LOD scores varied from 3.22 to 4.04.

Conclusions

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Electronic supplementary material

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

Genetic linkage map of the house musk shrew, Suncus murinus, constructed with PCR-based and RFLP markers.

A genetic linkage map for Suncus murinus was previously constructed with 11 marker loci. In this study, we developed 172 new microsatellite and three RFLP markers, and re-constructed a new framework map by combining all markers. The new map comprises 42 markers that are distributed into 12 linkage groups, two of which are assigned to chromosomes, and spans 403.5 cM with an average inter-marker distance of 13.5 cM.     

Construction of a genetic map for Theileria parva: identification of hotspots of recombination

The tick-borne protozoan parasite Theileria parva is the causal agent of East Coast Fever (ECF), a severe lymphoproliferative disease of cattle in eastern, central and southern Africa. The life cycle of T. parva is predominantly haploid, with a brief diploid stage occurring in the tick vector that involves meiotic recombination. Resolved genetic studies of T. parva are currently constrained by the lack of a genome-wide high-definition genetic map of the parasite. We undertook a genetic cross of two cloned isolates of T. parva to construct such a map from 35 recombinant progeny, using a genome-wide panel of 79 variable number of tandem repeat markers. Progeny were established by in vitro cloning of cattle lymphocytes after infection with sporozoites prepared from Rhipicephalus appendiculatus ticks fed on a calf undergoing a dual infection with the two clonal parental stocks. The genetic map was determined by assigning individual markers to the four chromosome genome, whose physical length is approximately 8309 kilobasepairs (Kb). Segregation analysis of the markers among the progeny revealed a total genetic size of 1683.8 centiMorgans (cM), covering a physical distance of 7737.62 Kb (∼93% of the genome). The average genome-wide recombination rate observed for T. parva was relatively high, at 0.22 cM Kb⁻¹ per meiotic generation. Recombination hot-spots and cold-spots were identified for each of the chromosomes. A panel of 27 loci encoding determinants previously identified as immunorelevant or likely to be under selection were positioned on the linkage map. We believe this to be the first genetic linkage map for T. parva. This resource, with the availability of the genome sequence of T. parva, will promote improved understanding of the pathogen by facilitating the use of genetic analysis for identification of loci responsible for variable phenotypic traits exhibited by individual parasite stocks.     

A Linkage Map of the Canine Genome

A genetic linkage map of the canine genome has been developed by typing 150 microsatellite markers using 17 three-generation pedigrees, composed of 163 F₂individuals. One hundred and thirty-nine markers were linked to at least one other marker with a lod score ≥ 3.0, identifying 30 linkage groups. The largest chromosome had 9 markers spanning 106.1 cM. The average distance between markers was 14.03 cM, and the map covers an estimated 2073 cM. Eleven markers were informative on the mapping panel, but were unlinked to any other marker. These likely represent single markers located on small, distinct canine chromosomes. This map will be the initial resource for mapping canine traits of interest and serve as a foundation for development of a comprehensive canine genetic map.     

A Simple Sequence Repeat- and Single-Nucleotide Polymorphism-Based Genetic Linkage Map of the Brown Planthopper,Nilaparvata lugens

In this study, we developed the first genetic linkage map for the major rice insect pest, the brown planthopper (BPH, Nilaparvata lugens). The linkage map was constructed by integrating linkage data from two backcross populations derived from three inbred BPH strains. The consensus map consists of 474 simple sequence repeats, 43 single-nucleotide polymorphisms, and 1 sequence-tagged site, for a total of 518 markers at 472 unique positions in 17 linkage groups. The linkage groups cover 1093.9 cM, with an average distance of 2.3 cM between loci. The average number of marker loci per linkage group was 27.8. The sex-linkage group was identified by exploiting X-linked and Y-specific markers. Our linkage map and the newly developed markers used to create it constitute an essential resource and a useful framework for future genetic analyses in BPH.     

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

Background

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

Results

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

Conclusions

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-152) contains supplementary material, which is available to authorized users.     

Map-based analysis of the tenacious glume gene Tg-B1 of wild emmer and its role in wheat domestication

The domestication of wheat was instrumental in spawning the civilization of humankind, and it occurred through genetic mutations that gave rise to types with non-fragile rachises, soft glumes, and free-threshing seed. Wild emmer (Triticum turgidum ssp. dicoccoides), the tetraploid AB-genome progenitor of domesticated wheat has genes that confer tenacious glumes (Tg) that underwent genetic mutations to give rise to free-threshing wheat. Here, we evaluated disomic substitution lines involving chromosomes 2A and 2B of wild emmer accessions substituted for homologous chromosomes in tetraploid and hexaploid backgrounds. The results suggested that both chromosomes 2A and 2B of wild emmer possess genes that inhibit threshability. A population of recombinant inbred lines derived from the tetraploid durum wheat variety Langdon crossed with a Langdon — T. turgidum ssp. dicoccoides accession PI 481521 chromosome 2B disomic substitution line was used to develop a genetic linkage map of 2B, evaluate the genetics of threshability, and map the gene derived from PI 481521 that inhibited threshability. A 2BS linkage map comprised of 58 markers was developed, and markers delineated the gene to a 2.3 cM interval. Comparative analysis with maps containing the tenacious glume gene Tg-D1 on chromosome arm 2DS from Aegilops tauschii, the D genome progenitor of hexaploid wheat, revealed that the gene inhibiting threshability in wild emmer was homoeologous to Tg-D1 and therefore designated Tg-B1. Comparative analysis with rice and Brachypodium distachyon indicated a high level of divergence and poorly conserved colinearity, particularly near the Tg-B1 locus. These results provide a foundation for further studies involving Tg-B1, which, together with Tg-D1, had profound influences on wheat domestication.     

High density SNP and SSR-based genetic maps of two independent oil palm hybrids

Background

Oil palm is an important perennial oil crop with an extremely long selection cycle of 10 to 12 years. As such, any tool that speeds up its genetic improvement process, such as marker-assisted breeding is invaluable. Previously, genetic linkage maps based on AFLP, RFLP and SSR markers were developed and QTLs for fatty acid composition and yield components identified. High density genetic maps of crosses of different genetic backgrounds are indispensable tools for investigating oil palm genetics. They are also useful for comparative mapping analyses to identify markers closely linked to traits of interest.

Results

A 4.5 K customized oil palm SNP array was developed using the Illumina Infinium platform. The SNPs and 252 SSRs were genotyped on two mapping populations, an intraspecific cross with 87 palms and an interspecific cross with 108 palms. Parental maps with 16 linkage groups (LGs), were constructed for the three fruit forms of E. guineensis (dura, pisifera and tenera). Map resolution was further increased by integrating the dura and pisifera maps into an intraspecific integrated map with 1,331 markers spanning 1,867 cM. We also report the first map of a Colombian E. oleifera, comprising 10 LGs with 65 markers spanning 471 cM. Although not very dense due to the high level of homozygosity in E. oleifera, the LGs were successfully integrated with the LGs of the tenera map. Direct comparison between the parental maps identified 603 transferable markers polymorphic in at least two of the parents. Further analysis revealed a high degree of marker transferability covering 1,075 cM, between the intra- and interspecific integrated maps. The interspecific cross displayed higher segregation distortion than the intraspecific cross. However, inclusion of distorted markers in the genetic maps did not disrupt the marker order and no map expansion was observed.

Conclusions

The high density SNP and SSR-based genetic maps reported in this paper have greatly improved marker density and genome coverage in comparison with the first reference map based on AFLP and SSR markers. Therefore, it is foreseen that they will be more useful for fine mapping of QTLs and whole genome association mapping studies in oil palm.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-309) 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.     

Genetic and physical mapping of new EST-derived SSRs on the A and B genome chromosomes of wheat

The availability of genetic maps and phenotypic data of segregating populations allows to localize and map agronomically important genes, and to identify closely associated molecular markers to be used in marker-assisted selection and positional cloning. The objective of the present work was to develop a durum wheat intervarietal genetic and physical map based on genomic microsatellite or genomic simple sequence repeats (gSSR) markers and expressed sequence tag (EST)-derived microsatellite (EST-SSR) markers. A set of 122 new EST-SSR loci amplified by 100 primer pairs was genetically mapped on the wheat A and B genome chromosomes. The whole map also comprises 149 gSSR markers amplified by 120 primer pairs used as anchor chromosome loci, two morphological markers (Black colour, Bla1, and spike glaucousness, Ws) and two seed storage protein loci (Gli-A2 and Gli-B2). The majority of SSR markers tested (182) was chromosome-specific. Out of 275 loci 241 loci assembled in 25 linkage groups assigned to the chromosomes of the A and B genome and 34 remained unlinked. A higher percentage of markers (54.4%), localized on the B genome chromosomes, in comparison to 45.6% distributed on the A genome. The whole map covered 1,605 cM. The B genome accounted for 852.2 cM of genetic distance; the A genome basic map spanned 753.1 cM with a minimum length of 46.6 cM for chromosome 5A and a maximum of 156.2 cM for chromosome 3A and an average value of 114.5 cM. The primer sets that amplified two or more loci mapped to homoeologous as well as to non-homoeologous sites. Out of 241 genetically mapped loci 213 (88.4%) were physically mapped by using the nulli-tetrasomic, ditelosomic and a stock of 58 deletion lines dividing the A and B genome chromosomes in 94 bins. No discrepancies concerning marker order were observed but the cytogenetic maps revealed in some cases small genetic distance covered large physical regions. Putative function for mapped SSRs were assigned by searching against GenBank nonredundant database using TBLASTX algorithms.     

A comparative genetic map of the turkey genome

Genetic markers (microsatellites and SNPs) were used to create and compare maps of the turkey and chicken genomes. A physical map of the chicken genome was built by comparing sequences of turkey markers with the chicken whole-genome sequence by BLAST analysis. A genetic linkage map of the turkey genome (Meleagris gallopavo) was developed by segregation analysis of genetic markers within the University of Minnesota/Nicholas Turkey Breeding Farms (UMN/NTBF) resource population. This linkage map of the turkey genome includes 314 loci arranged into 29 linkage groups. An additional 40 markers are tentatively placed within linkage groups based on two-point LOD scores and 16 markers remain unlinked. Total map distance contained within linkage groups is 2,011 cM with the longest linkage group (47 loci) measuring 413.3 cM. Average marker interval over the 29 linkage groups was 6.4 cM. All but one turkey linkage group could be aligned with the physical map of the chicken genome. The present genetic map of the turkey provides a comparative framework for future genomic studies.     

Construction of a high density integrated genetic map for cucumber (Cucumis sativus L.)

The high-density consensus map was constructed based on the GY14 × PI 183967 map from an inter-subspecific cross and the extended S94 × S06 map from an intra-subspecific cross. The consensus map was composed of 1,369 loci, including 1,152 SSR loci, 192 SRAP loci, 21 SCAR loci and one STS locus as well as three gene loci of fruit external quality traits in seven chromosomes, and spanned 700.5 cM, of which 682.7 cM (97.5%) were covered by SSR markers. The average genetic distance and physical interval between loci were 0.51 cM and ~268 kbp, respectively. Additionally, the physical position of the sequence-associated markers aligned along the assembled cucumber genome sequence established a relationship between genetic maps and cucumber genome sequence and to a great extent validated the order of markers in individual maps and consensus map. This consensus map with a high marker density and well-ordered markers is a saturated and reliable linkage map for genetic analysis of cucumber or the Cucurbitaceae family of plants.     

A first-generation microsatellite linkage map of the Japanese quail

A linkage map of the Japanese quail (Coturnix japonica) genome was constructed based upon segregation analysis of 72 microsatellite loci in 433 F(2) progeny of 10 half-sib families obtained from a cross between two quail lines of different genetic origins. One line was selected for long duration of tonic immobility, a behavioural trait related to fearfulness, while the other was selected based on early egg production. Fifty-eight of the markers were resolved into 12 autosomal linkage groups and a Z chromosome-specific linkage group, while the remaining 14 markers were unlinked. The linkage groups range from 8 cM (two markers) to 206 cM (16 markers) and cover a total map distance of 576 cM with an average spacing of 10 cM between loci. Through comparative mapping with chicken (Gallus gallus) using orthologous markers, we were able to assign linkage groups CJA01, CJA02, CJA05, CJA06, CJA14 and CJA27 to chromosomes. This map, which is the first in quail based solely on microsatellites, is a major step towards the development of a quality molecular genetic map for this valuable species. It will provide an important framework for further genetic mapping and the identification of quantitative trait loci controlling egg production and fear-related behavioural traits in quail.     

Molecular marker analysis of Helianthus annuus L. 2. Construction of an RFLP linkage map for cultivated sunflower

A detailed linkage map of Helianthus annuus was constructed based on segregation at 234 RFLP loci, detected by 213 probes, in an F₂ population of 289 individuals (derived from a cross between the inbred lines HA89 and ZENB8). The genetic markers covered 1380 centiMorgans (cM) of the sunflower genome and were aranged in 17 linkage groups, corresponding to the haploid number of chromosomes in this species. One locus was found to be unlinked. Although the average interval size was 5.9 cM, there were a number of regions larger than 20 cM that were devoid of markers. Genotypic classes at 23 loci deviated significantly from the expected ratios (121 or 31), all showing a reduction in the ZENB8 homozygous class. The majority of these loci were found to map to four regions on linkage groups G, L and P.     

An integrated and comparative genetic map of the turkey genome

An integrated genetic linkage map was developed for the turkey (Meleagris gallopavo) that combines the genetic markers from the three previous mapping efforts. The UMN integrated map includes 613 loci arranged into 41 linkage groups. An additional 105 markers are tentatively placed within linkage groups based on two-point LOD scores and 19 markers remain unlinked. A total of 210 previously unmapped markers has been added to the UMN turkey genetic map. Markers from each of the 20 linkage groups identified in the Roslin map and the 22 linkage groups of the Nte map are incorporated into the new integrated map. Overall map distance contained within the 41 linkage groups is 3,365 cM (sex-averaged) with the largest linkage group (94 loci) measuring 533.1 cM. Average marker interval for the map was 7.86 cM. Sequences of markers included in the new map were compared to the chicken genome sequence by 'BLASTN'. Significant similarity scores were obtained for 95.6% of the turkey sequences encompassing an estimated 91% of the chicken genome. A physical map of the chicken genome based on positions of the turkey sequences was built and 36 of the 41 turkey linkage groups were aligned with the physical map, five linkage groups remain unassigned. Given the close similarities between the turkey and chicken genomes, the chicken genome sequence could serve as a scaffold for a genome sequencing effort in the turkey.     

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.     

Moderate-density molecular maps of Eucalyptus urophylla S. T. Blake and E. tereticornis Smith genomes based on RAPD markers

Moderate-density molecular maps were constructed for the genomes of Eucalyptus urophylla S. T. Blake and E. tereticornis Smith using RAPD markers and an interspecific cross between the two species. One hundred and eighty-three primers were employed to generate 245 and 264 parent-specific markers in E. urophylla and E. tereticornis, respectively, as well as 49 parent-shared markers. The normally segregating markers, including 208 (84.9%) specific to maternal E. urophylla, 175 (66.3%) to paternal E. tereticornis, and 48 shared by both parents, were used for framework map construction for each parental species. For maternal E. urophylla, the linkage map consisted of 23 linkage groups, 160 framework markers, and 60 accessory markers, defining a total map distance of 1504.6 cM and an average interval of 11.0 ± 8.07 cM. For paternal E. tereticornis, the linkage map contained 23 linkage groups, 126 framework markers, and 92 accessory markers, defining a total map distance of 1035.7 cM and an average interval of 10.1 ± 7.23 cM. Genome length was estimated at 1585.7 and 1507.5 cM for E. urophylla and E. tereticornis, respectively, indicating map coverage of 94.9 and 68.7% of the corresponding genomes. Construction of such maps will be valuable for quantitative trait loci (QTLs) detection, marker-assisted selection (MAS), comparative mapping, and whole genome based fingerprint characterization in Eucalyptus breeding programs.     

The First High-Density Genetic Map Construction in Tree Peony (Paeonia Sect. Moutan) using Genotyping by Specific-Locus Amplified Fragment Sequencing

Genetic linkage maps, permitting the elucidation of genome structure, are one of most powerful genomic tools to accelerate marker-assisted breeding. However, due to a lack of sufficient user-friendly molecular markers, no genetic linkage map has been developed for tree peonies (Paeonia Sect. Moutan), a group of important horticultural plants worldwide. Specific-locus amplified fragment sequencing (SLAF-seq) is a recent molecular marker development technology that enable the large-scale discovery and genotyping of sequence-based marker in genome-wide. In this study, we performed SLAF sequencing of an F1 population, derived from the cross P. ostti ‘FenDanBai’ × P. × suffruticosa ‘HongQiao’, to identify sufficient high-quality markers for the construction of high-density genetic linkage map in tree peonies. After SLAF sequencing, a total of 78 Gb sequencing data and 285,403,225 pair-end reads were generated. We detected 309,198 high-quality SLAFs from these data, of which 85,124 (27.5%) were polymorphic. Subsequently, 3518 of the polymorphic markers, which were successfully encoded in to Mendelian segregation types, and were in conformity with the criteria of high-quality markers, were defined as effective markers and used for genetic linkage mapping. Finally, we constructed an integrated genetic map, which comprised 1189 markers on the five linkage groups, and spanned 920.699 centiMorgans (cM) with an average inter-marker distance of 0.774 cM. There were 1115 ‘SNP-only’ markers, 18 ‘InDel-only’ markers, and 56 ‘SNP&InDel’ markers on the map. Among these markers, 450 (37.85%) showed significant segregation distortion (P < 0.05). In conclusion, this investigation reported the first large-scale marker development and high-density linkage map construction for tree peony. The results of this study will serve as a solid foundation not only for marker-assisted breeding, but also for genome sequence assembly for tree peony.