Sequence-based ultra-dense genetic and physical maps reveal structural variations of allopolyploid cotton genomes

BackgroundSNPs are the most abundant polymorphism type, and have been explored in many crop genomic studies, including rice and maize. SNP discovery in allotetraploid cotton genomes has lagged behind that of other crops due to their complexity and polyploidy. In this study, genome-wide SNPs are detected systematically using next-generation sequencing and efficient SNP genotyping methods, and used to construct a linkage map and characterize the structural variations in polyploid cotton genomes.ResultsWe construct an ultra-dense inter-specific genetic map comprising 4,999,048 SNP loci distributed unevenly in 26 allotetraploid cotton linkage groups and covering 4,042 cM. The map is used to order tetraploid cotton genome scaffolds for accurate assembly of G. hirsutum acc. TM-1. Recombination rates and hotspots are identified across the cotton genome by comparing the assembled draft sequence and the genetic map. Using this map, genome rearrangements and centromeric regions are identified in tetraploid cotton by combining information from the publicly-available G. raimondii genome with fluorescent in situ hybridization analysis.ConclusionsWe report the genotype-by-sequencing method used to identify millions of SNPs between G. hirsutum and G. barbadense. We construct and use an ultra-dense SNP map to correct sequence mis-assemblies, merge scaffolds into pseudomolecules corresponding to chromosomes, detect genome rearrangements, and identify centromeric regions in allotetraploid cottons. We find that the centromeric retro-element sequence of tetraploid cotton derived from the D subgenome progenitor might have invaded the A subgenome centromeres after allotetrapolyploid formation. This study serves as a valuable genomic resource for genetic research and breeding of cotton.

Electronic supplementary material

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

A set of EST-SNPs for map saturation and cultivar identification in melon

Background

There are few genomic tools available in melon (Cucumis melo L.), a member of the Cucurbitaceae, despite its importance as a crop. Among these tools, genetic maps have been constructed mainly using marker types such as simple sequence repeats (SSR), restriction fragment length polymorphisms (RFLP) and amplified fragment length polymorphisms (AFLP) in different mapping populations. There is a growing need for saturating the genetic map with single nucleotide polymorphisms (SNP), more amenable for high throughput analysis, especially if these markers are located in gene coding regions, to provide functional markers. Expressed sequence tags (ESTs) from melon are available in public databases, and resequencing ESTs or validating SNPs detected in silico are excellent ways to discover SNPs.

Results

EST-based SNPs were discovered after resequencing ESTs between the parental lines of the PI 161375 (SC) × 'Piel de sapo' (PS) genetic map or using in silico SNP information from EST databases. In total 200 EST-based SNPs were mapped in the melon genetic map using a bin-mapping strategy, increasing the map density to 2.35 cM/marker. A subset of 45 SNPs was used to study variation in a panel of 48 melon accessions covering a wide range of the genetic diversity of the species. SNP analysis correctly reflected the genetic relationships compared with other marker systems, being able to distinguish all the accessions and cultivars.

Conclusion

This is the first example of a genetic map in a cucurbit species that includes a major set of SNP markers discovered using ESTs. The PI 161375 × 'Piel de sapo' melon genetic map has around 700 markers, of which more than 500 are gene-based markers (SNP, RFLP and SSR). This genetic map will be a central tool for the construction of the melon physical map, the step prior to sequencing the complete genome. Using the set of SNP markers, it was possible to define the genetic relationships within a collection of forty-eight melon accessions as efficiently as with SSR markers, and these markers may also be useful for cultivar identification in Occidental melon varieties.     

Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq)

Background

Cultivated peanut, or groundnut (Arachis hypogaea L.), is an important oilseed crop with an allotetraploid genome (AABB, 2n = 4x = 40). In recent years, many efforts have been made to construct linkage maps in cultivated peanut, but almost all of these maps were constructed using low-throughput molecular markers, and most show a low density, directly influencing the value of their applications. With advances in next-generation sequencing (NGS) technology, the construction of high-density genetic maps has become more achievable in a cost-effective and rapid manner. The objective of this study was to establish a high-density single nucleotide polymorphism (SNP)-based genetic map for cultivated peanut by analyzing next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq) reads.

Results

We constructed reduced representation libraries (RRLs) for two A. hypogaea lines and 166 of their recombinant inbred line (RIL) progenies using the ddRADseq technique. Approximately 175 gigabases of data containing 952,679,665 paired-end reads were obtained following Solexa sequencing. Mining this dataset, 53,257 SNPs were detected between the parents, of which 14,663 SNPs were also detected in the population, and 1,765 of the obtained polymorphic markers met the requirements for use in the construction of a genetic map. Among 50 randomly selected in silico SNPs, 47 were able to be successfully validated. One linkage map was constructed, which was comprised of 1,685 marker loci, including 1,621 SNPs and 64 simple sequence repeat (SSR) markers. The map displayed a distribution of the markers into 20 linkage groups (LGs A01–A10 and B01–B10), spanning a distance of 1,446.7 cM. The alignment of the LGs from this map was shown in comparison with a previously integrated consensus map from peanut.

Conclusions

This study showed that the ddRAD library combined with NGS allowed the rapid discovery of a large number of SNPs in the cultivated peanut. The first high density SNP-based linkage map for A. hypogaea was generated that can serve as a reference map for cultivated Arachis species and will be useful in genetic mapping. Our results contribute to the available molecular marker resources and to the assembly of a reference genome sequence for the peanut.

Electronic supplementary material

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

Fine mapping QMi-C11 a major QTL controlling root-knot nematodes resistance in Upland cotton

The identification and utilization of a high-level of host plant resistance is the most effective and economical approach to control root-knot nematode (Meloidogyne incognita). In an earlier study, we identified a major quantitative trait locus (QTL) for resistance to root-knot nematode in the M-120 RNR Upland cotton line (Gossypium hirsutum L.) of the Auburn 623 RNR source. The QTL is located in a 12.9-cM interval flanked by the two SSR markers CIR069 and CIR316 on the distal segment of chromosome 11. To construct a fine map around the target region, a bulked segregation analysis was performed using two DNA pools consisting of five individuals, with each being homozygous for the two parental alleles. From a survey of 1,152 AFLP primer combinations, 9 AFLP markers closely linked to the target region were identified. By screening an additional 1,221 F₂ individuals developed from the initial mapping population, the Mi-C11 locus was delimited to a 3.6-cM interval flanked by the SSR marker CIR069 and the AFLP marker E14M27-375. These results further elucidate the genetic fine structure of the Mi-C11 locus and provide the basis for map-based isolation of the nematode resistance gene in M-120 RNR.     

Construction of a SNP and SSR linkage map in autotetraploid blueberry using genotyping by sequencing

The construction of the first genetic map in autotetraploid blueberry has been made possible by the development of new SNP markers developed using genotyping by sequencing in a mapping population created from a cross between two key highbush blueberry cultivars, Draper × Jewel (Vaccinium corymbosum). The novel SNP markers were supplemented with existing SSR markers to enable the alignment of parental maps.  In total, 1794 single nucleotide polymorphic (SNP) markers and 233 simple sequence repeat (SSR) markers exhibited segregation patterns consistent with a random chromosomal segregation model for meiosis in an autotetraploid. Of these, 700 SNPs and 85 SSRs were utilized for construction of the ‘Draper’ genetic map, and 450 SNPs and 86 SSRs for the ‘Jewel’ map.  The ‘Draper’ map comprises 12  linkage groups (LG), associated with the haploid chromosome number for blueberry, and totals 1621 cM while the ‘Jewel’ map comprises 20 linkage groups totalling 1610 cM. Tentative alignments of the two parental maps have been made on the basis of shared SSR alleles and linkages to double-simplex markers segregating in both parents. Tentative alignments of the two parental maps have been made on the basis of shared SSR alleles and linkages to double-simplex markers segregating in both parents.     

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin–streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N₅₀ of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F₂ population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome (‘A’ or ‘D’) each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.     

Construction of an SSR and RAD Marker-Based Genetic Linkage Map for Carnation (<Emphasis Type="Italic">Dianthus caryophyllus</Emphasis> L.)

A high-density genetic map, an essential tool for comparative genomic studies and quantitative trait locus fine mapping, can also facilitate genome sequence assembly. The sequence-based marker technology known as restriction site-associated DNA (RAD) enables synchronous, single nucleotide polymorphism marker discovery, and genotyping using massively parallel sequencing. We constructed a high-density linkage map for carnation (Dianthus caryophyllus L.) based on simple sequence repeat (SSR) markers in combination with RAD markers developed by double-digest RAD sequencing (ddRAD-seq). A total of 2404 (285 SSR and 2119 RAD) markers could be assigned to 15 linkage groups spanning 971.5 cM, with an average marker interval of 0.4 cM. The total length of scaffolds with identified map positions was 95.6 Mb, which is equivalent to 15.4 % of the estimated genome size. The generated map is the first SSR and RAD marker-based high-density linkage map reported for carnation. The ddRAD-seq pipeline developed in this study should also help accelerate genetic and genomics analyses and molecular breeding of carnation and other non-model crops.     

A SNP haplotype associated with a gene resistant to Xanthomonas axonopodis pv. malvacearum in upland cotton (Gossypium hirsutum L.)

An F_4:5 population of 285 families with each tracing back to a different F₂ plant, derived from a cotton bacterial blight resistant line ‘DeltaOpal’ and a susceptible line ‘DP388’, was artificially inoculated with bacterial blight race 18 (Xanthomonas axonopodis pv. malvacearum) to assay their resistance or susceptibility to the disease. The segregation in the F_4:5 population indicates that the resistance was conditioned by a single dominant gene designated B _12. Simple sequence repeat (SSR) markers identified as putatively linked to the resistance gene by bulked segregant analysis were confirmed on the entire F_4:5 population. Three SSR markers, CIR246, BNL3545 and BNL3644 on chromosome 14, were found closely linked to B ₁₂ . The association between CIR246 and B ₁₂ was validated among 354 plants of 16 diverse varieties. Based on Monsanto SSR/single nucleotide polymorphism (SNP) consensus map, SNP markers closely linked to CIR246 were used to screen ‘DeltaOpal’ and ‘DP388’ for polymorphism. The polymorphic SNP markers were run on the F_4:5 population and the four SNP markers spanning 3.4 cM were found to flank the resistance gene on chromosome 14. The linkage between B ₁₂ and the 4-SNP marker haplotype was validated using 18 elite cotton lines. This 4-SNP marker haplotype can be used for marker assisted selection for bacterial blight resistance breeding programs or for screening germplasm collections for this locus rapidly.     

Construction of a high-density SNP genetic map in flue-cured tobacco based on SLAF-seq

Tobacco (Nicotiana tabacum L., 2n = 48) is an important agronomic crop and model plant. Flue-cured tobacco is the most important type and accounts for approximately 80 % of tobacco production worldwide. The low genetic diversity of flue-cured tobacco impedes the construction of a high-density genetic linkage map using simple sequence repeat (SSR) markers and warrants the exploitation of single nucleotide polymorphic (SNP) markers from genomic regions. In this article, initially using specific locus-amplified fragment sequencing, we discovered 10,891 SNPs that were subsequently used as molecular markers for genetic map construction. Combined with SSR markers, a final high-density genetic map was generated containing 4215 SNPs and 194 SSRs distributed on 24 linkage groups (LGs). The genetic map was 2662.43 cM in length, with an average distance of 0.60 cM between adjacent markers. Furthermore, by mapping the SNP markers to the ancestral genomes of Nicotiana tomentosiformis and Nicotiana sylvestris, a large number of genome rearrangements were identified as occurring after the polyploidization event. Finally, using this novel integrated map and mapping population, two major quantitative trait loci (QTLs) were identified for flue-curing and mapped to the LG6 of tobacco. This is the first report of SNP markers and a SNP-based linkage map being developed in tobacco. The high-density genetic map and QTLs related to tobacco curing will support gene/QTL fine mapping, genome sequence assembly and molecular breeding in tobacco.     

QTL mapping for bacterial wilt resistance in peanut (<Emphasis Type="Italic">Arachis hypogaea</Emphasis> L.)

Bacterial wilt (BW) caused by Ralstonia solanacearum is a serious, global, disease of peanut (Arachis hypogaea L.), but it is especially destructive in China. Identification of DNA markers linked to the resistance to this disease will help peanut breeders efficiently develop resistant cultivars through molecular breeding. A F₂ population, from a cross between disease-resistant and disease-susceptible cultivars, was used to detect quantitative trait loci (QTL) associated with the resistance to this disease in the cultivated peanut. Genome-wide SNPs were identified from restriction-site-associated DNA sequencing tags using next-generation DNA sequencing technology. SNPs linked to disease resistance were determined in two bulks of 30 resistant and 30 susceptible plants along with two parental plants using bulk segregant analysis. Polymorphic SSR and SNP markers were utilized for construction of a linkage map and for performing the QTL analysis, and a moderately dense linkage map was constructed in the F₂ population. Two QTL (qBW-1 and qBW-2) detected for resistance to BW disease were located in the linkage groups LG1 and LG10 and account for 21 and 12 % of the bacterial wilt phenotypic variance. To confirm these QTL, the F₈ RIL population with 223 plants was utilized for genotyping and phenotyping plants by year and location as compared to the F₂ population. The QTL qBW-1 was consistent in the location of LG1 in the F₈ population though the QTL qBW-2 could not be clarified due to fewer markers used and mapped in LG10. The QTL qBW-1, including four linked SNP markers and one SSR marker within 14.4-cM interval in the F₈, was closely related to a disease resistance gene homolog and was considered as a candidate gene for resistance to BW. QTL identified in this study would be useful to conduct marker-assisted selection and may permit cloning of resistance genes. Our study shows that bulk segregant analysis of genome-wide SNPs is a useful approach to expedite the identification of genetic markers linked to disease resistance traits in the allotetraploidy species peanut.     

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.     

Development of gene-targeted SNP markers for genomic mapping in Pacific abalone Haliotis discus hannai Ino

Useful and novel DNA markers are needed for aquaculture genetics and breeding. In this study, we report the discovery and development of gene-targeted single nucleotide polymorphisms (SNPs) for genomic mapping in the Pacific abalone Haliotis discus hannai Ino. Single EST or EST-contigs from 66 genes that had positive BLASTx matches (E-value ≤ 1e-8) were used for polymerase chain reaction (PCR) amplification. PCR products from the two parents of one mapping family were directly sequenced, and 83 SNP loci were found from 17 genes. Allele-specific PCR (AS-PCR) was developed and optimized for genotyping of 11 SNP loci in 120 progeny of the mapping family. Nine of the loci conformed to the expected Mendelian ratio of 1:1 based on the χ² test (P > 0.05) and could potentially be used for linkage map construction. Our data also indicate that the sequencing of two parents may be a practical strategy for the discovery of informative SNPs for linkage mapping in a particular mapping population.     

Molecular linkage map of allotetraploid cotton (<Emphasis Type="Italic">Gossypium hirsutum</Emphasis> L. × <Emphasis Type="Italic">Gossypium barbadense</Emphasis> L.) with a haploid population

In the present study, a haploid population from the cross of the two cultivated allotetraploid cottons, Gossypium hirsutum L. and Gossypium barbadense L., was developed by means of Vsg, a virescently marked semigamous line of Sea island cotton, and some target haploids were successfully doubled with colchicine. A molecular linkage map was constructed with 58 doubled and haploid plants. Among the total of 624 marker loci (510 SSRs and 114 RAPDs), 489 loci were assembled into 43 linkage groups and covered 3,314.5 centi-Morgans (cM). Using the monosomic and telodisomic genetic stocks, the linkage groups of the present map were associated with chromosomes of the allotetraploid genome, and some of the unassociated groups were connected to corresponding A or D subgenomes. Through the analysis of the assignment of the duplicated SSR loci in the chromosomes or the linkage groups, ten pairs of possible homoeologous chromosome (or linkage group) regions were identified. Among them, the pairs of Chrs. 1 and 15, Chrs. 4 and 22, and Chrs. 10 and 20 had already been determined as homoeologous by classical genetic and cytogenetic research, and the pair of Chrs. 9 and 23 had also been identified by the ISH method of molecular cytogenetics. But, from present research, it was assumed that Chrs. 5 and 18 might be a new pair of homoeologous chromosomes of the allotetraploid cotton genome detected by molecular mapping of the cotton genome.     

Enrichment of an intraspecific genetic map of upland cotton by developing markers using parental RAD sequencing

RAD sequencing was performed using DH962 and Jimian5 as upland cotton mapping parents. Sequencing data for DH962 and Jimian5 were assembled into the genome sequences of ≈55.27 and ≈57.06 Mb, respectively. Analysing genome sequences of the two parents, 1,323 SSR, 3,838 insertion/deletion (InDel), and 9,366 single-nucleotide polymorphism (SNP) primer pairs were developed. All of the SSRs, 121 InDels, 441 SNPs, and other 6,747 primer pairs were screened in the two parents, and a total of 535 new polymorphic loci were identified. A genetic map including 1,013 loci was constructed using these results and 506 loci previously published for this population. Twenty-seven new QTLs for yield and fibre quality were identified, indicating that the efficiency of QTL detection was greatly improved by the increase in map density. Comparative genomics showed there to be considerable homology and collinearity between the A_T and A₂ genomes and between the D_T and D₅ genomes, although there were a few exchanges and introgressions among the chromosomes of the A₂ genome. Here, the development of markers using parental RAD sequencing was effective, and a high-density intraspecific genetic map was constructed. This map can be used for molecular marker-assisted selection in cotton.     

Cot-based sampling of genomes for polymorphic low-copy DNA

DNA polymorphisms are powerful tools for many evolutionary and genomic studies in plants including molecular breeding. Single nucleotide polymorphisms (SNPs) are the most elemental DNA marker for genomic studies, but even with advances in DNA sequencing technology, SNP discovery remains costly and computationally demanding, especially in large genomes that are rich in repetitive DNA such as those of many plants. Here we report a method using DNA renaturation kinetics (Cot techniques), sequencing, and BLAST-based screening to identify low-copy, non-coding DNA sequences that were subsequently found to be relatively rich in polymorphisms. A total of of 63 such fragments isolated from a diploid D genome cotton species (Gossypium raimondii) revealed a higher frequency of polymorphisms than that observed for cotton expressed sequence tags or hypomethylated (PstI-susceptible) genomic DNA. While microsatellite-derived loci show still higher polymorphism rates, they often fall in repetitive elements and their sequence analysis is often complicated by alignment difficulties. The potential applications of Cot-filtered noncoding (CFNC) DNA in development of DNA markers are discussed.