Next Article in Journal
Protein Phosphatase 2A in the Regulation of Wnt Signaling, Stem Cells, and Cancer
Previous Article in Journal
Relating Stool Microbial Metabolite Levels, Inflammatory Markers and Dietary Behaviors to Screening Colonoscopy Findings in a Racially/Ethnically Diverse Patient Population
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 9
Issue 3
10.3390/genes9030120
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes)

by
Xiang Zhang
1,
Misaki Mizukoshi
1,
Hong Zhang
1,
Engkong Tan
1,
Yoji Igarashi
1,
Yutaka Suzuki
2,
Susumu Mitsuyama
1,
Shigeharu Kinoshita
1,
Kazuyoshi Saito
3,
Shugo Watabe
4 and
Shuichi Asakawa
1,*
1
Department of Aquatic Bioscience, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Bunkyo, Tokyo 113-8657, Japan
2
Department of Medical Genome Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Kashiwa, Chiba 277-8562, Japan
3
Akita Prefectural Fisheries Promotion Center, Oga, Akita 010-0531, Japan
4
School of Marine Bioscience, Kitasato University, Sagamihara, Kanagawa 252-0373, Japan
*
Author to whom correspondence should be addressed.
Genes 2018, 9(3), 120; https://doi.org/10.3390/genes9030120
Submission received: 20 December 2017 / Revised: 15 February 2018 / Accepted: 20 February 2018 / Published: 26 February 2018
(This article belongs to the Section Animal Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Next-generation sequencing enables genome-wide genotyping of a large population and further facilitates the construction of a genetic linkage map. Low-coverage whole-genome sequencing has been employed for genetic linkage map construction in several species. However, this strategy generally requires available high-quality reference genomes and/or designed inbred pedigree lines, which restrict the scope of application for non-model and unsequenced species. Here, using torafugu (Takifugu rubripes) as a test model, we propose a new strategy for ultrahigh-density genetic linkage map construction using low-coverage whole-genome sequencing of a haploid/doubled haploid (H/DH) population without above requirements. Low-coverage (≈1×) whole-genome sequencing data of 165 DH individuals were used for de novo assembly and further performed single nucleotide polymorphisms (SNPs) calling, resulting in the identification of 1,070,601 SNPs. Based on SNP genotypes and de novo assembly, genotypes were associated with short DNA segments and an ultrahigh-density linkage map was constructed containing information of 802,277 SNPs in 3090 unique positions. Comparative analyses showed near-perfect concordance between the present linkage map and the latest published torafugu genome (FUGU5). This strategy would facilitate ultrahigh-density linkage map construction in various sexually reproducing organisms for which H/DH populations can be generated.

1. Introduction

A genetic linkage map is a powerful tool in genetic and genomic research. It lays a strong foundation for comparative genomics and provides vital clues toward understanding genome evolution and divergence [1,2,3,4]. Moreover, it facilitates genotype–phenotype association mapping and enables investigating the genetics of complex phenotypic traits [5,6,7]. It also contributes toward characterization of genome structure and serves as the backbone for anchoring unplaced/misplaced scaffolds for chromosome-scale assembly [8,9,10,11].
Owing to the rapid development of next-generation sequencing (NGS) in the last decade, the ability to simultaneously sequence a large number of individuals in a multiplex manner has now become possible, so that an entire population can be rapidly genotyped for linkage mapping [12]. Since 2010, low-coverage whole-genome resequencing has been employed for construction of the genetic linkage maps of rice [13], shiitake mushroom [14], and safflower [15]. However, all of these cases relied on the availability of high-quality reference genome sequences and/or designed inbred pedigree lines to carry out the linkage mapping prior to resequencing of the mapping population. This requirement currently limits the wide application of low-coverage whole-genome resequencing in non-model organisms, especially for those with unexplored genomes.
To resolve this issue, we hypothesized that an ultrahigh-density genetic linkage map could be constructed using low-coverage whole-genome sequencing of haploid/doubled-haploid (H/DH) population without requiring a high-quality reference genome and the laborious establishment of inbred lines. H/DH individuals can be generated via natural or artificial uniparental reproduction, which is found in a wide range of species in several kingdoms. Indeed, in recent years, the H/DH population has been exploited as an ideal population type for genetic linkage map construction, particularly in plants [16,17,18,19,20] and teleosts [21,22,23,24] due to their well-developed H/DH technologies [25,26]. The main advantage of H/DH individuals for genotype sequencing is that a relatively low sequencing coverage is sufficient without loss of accuracy compared to the coverage necessary for sequencing more common diploid individuals owing to the presence of heterozygous single nucleotide polymorphisms (SNPs). To test this strategy, we established a DH population of the wild pufferfish or torafugu (Takifugu rubripes) through gynogenesis.
Torafugu is a popular species with economic importance in the waters of East Asia, and has emerged as an ideal model in genomic studies owing to its compact genome [27]. In fact, the torafugu genome is considered to be one of the smallest (≈400 Mb) among vertebrates and is approximately eight times smaller than the human genome [27]. Another advantage of torafugu as a model for genetic analysis is its similarity to mammals, including a shared body plan and physiological systems. Thus, the compact genome can favor the discovery of genes and gene regulatory regions with clear counterparts in the human genome, and the torafugu genome can further serve as a reference to understand the structure, function, and evolution of vertebrate genomes [28,29]. In the most recent fifth version of the torafugu genome assembly (FUGU5), 72% of the scaffolds have been located and oriented after integration with the torafugu genetic linkage map, comprising 1220 microsatellite markers; however, the remaining 14% have been located but not oriented, and the other 14% have not yet been assigned [30]. Therefore, the construction of a higher-density genetic linkage map of torafugu is needed to be able to expand the contiguity and improve the quality of the genome assembly.
Based on the proposed strategy, an ultrahigh-density linkage map of torafugu was constructed. The accuracy of the obtained linkage map was validated with comparison to the published genome FUGU5. The proposed strategy represents a cost-effective and less complex tool for genetic linkage map construction and can be widely applied in a wide diversity of sexual organisms, especially non-model and unsequenced species, for which H/DH populations can be generated.

2. Materials and Methods

2.1. Preparation of a Doubled-Haploid Population of Torafugu

A wild female torafugu was purchased from a market in Akita Prefecture and was subjected to mito-gynogenesis for generating a DH population according to the process described in detail in our previous paper [31]. In brief, mature oocytes were fertilized with inactive sperm of a male torafugu (from the same market) that had been pretreated with ultraviolet radiation at 40 mJ/cm2. After fertilization for 3 h, the eggs were subjected to 45 min of cold-shock treatment at 0.6 °C, followed by incubation in aerated tanks with fresh seawater at 18.0 °C. Several days after artificial insemination, hundreds of eggs were observed to contain embryonic bodies, which were selected for further analysis.

2.2. Whole-Genome Sequencing

Genomic DNA was extracted from each of the selected 192 eggs using Agencourt DNAdvance Kit (Beckman Coulter Genomics, Beverly, MA, USA) after homogenization. An average of 125 ng DNA was obtained from each sample. DNA libraries of these individuals were prepared and barcoded according to the Nextera DNA Library Prep Reference Guide (Illumina, San Diego, CA, USA) and were then subjected to sequencing in two lanes of Illumina HiSeq 2000 system. A total of 74.43 Gb of sequencing data, consisting of 2 × 100-bp paired-end reads with an average insert size of 230 bp, were obtained from 192 samples of the generated DH torafugu population. Potential remnants of adapter sequences were removed, low-quality bases with a Phred quality score below 20 were trimmed, and the 23 samples with very low sequencing coverage were removed. After these processes, a total of 71.32 Gb of sequencing data from 169 samples were ultimately reserved and applied to further analysis.

2.3. De Novo Assembly and SNP Calling

The obtained sequencing data were utilized to perform de novo assembly on the SOAPdenovo2 [32] assembler under a k-mer value of 58. The sequencing data of each sample were mapped to the obtained de novo assembly using Burrows–Wheeler Aligner (BWA) [33], followed by SNP calling using SAMtools with default parameters [34]. Four samples were identified as partial diploids due to the existence of many heterologous SNPs. Considering that DNA polymorphisms would have effects on de novo assembly, the sequencing data (1.74 Gb) from these four samples were removed. Thus, a total of 69.58 Gb sequencing data from 165 samples were reserved. According to the torafugu genome size (approximately 400 Mb), the total sequencing data coverage was estimated at 174, whereas the average coverage for each sample was 1.05 ± 0.76. The remaining sequencing data (69.58 Gb) were utilized to perform a second round of de novo assembly, and SNP calling was then performed for each sample with the same parameters.

2.4. Marker Genotype Scoring

Owing to the low-coverage (≈1×) sequencing, most of the SNPs in each sample were detected once or less, which led to a large quantity of missed genotypes and insufficient data for genotyping calibration. As shown in Figure 1, we arbitrarily assigned a phase to every allele of each SNP (with one type scored as “A”, the other scored as “B”, and unknown scored as “-”) because of the unknown linkage phase. Subsequently, a low-call-rate SNP dataset of unknown phase was obtained. However, based on the contig/scaffold of the above de novo assembly, the genotype of adjacent SNPs could be testified and/or compensated by each other. The SNPs located on short DNA segments were phased at the segment level and merged together as a new genetic marker termed the short segment genotype (SSG). As shown in Figure 2, when recombination occurred outside of a short DNA segment (Seg1 in Sample 1, Seg1 and Seg2 in Sample 2), in most cases, the converted SSG possessed the same genotypes at the SNP sites within it. By contrast, when recombination occurred on a DNA segment such as Seg2 in Sample 1, the converted SSG was regarded to possess the same genotype as the genotype of SNPs dominating within the segment. In that case, the site of crossover was temporarily located to the terminal end of the DNA segment. Practically, this temporal operation would hardly influence the formation of the linkage map because the length of the DNA segment is much shorter than that of the entire chromosome. In addition, two crossover events hardly ever occur within such a short DNA segment. The maximum segment length of the SSGs was set to 8 kb, so that the SSGs (>0.9 call rate) could harbor as much genome-wide SNPs information as possible while maintaining the length of segments as short as possible (see details in Figure S1). A high-call-rate SSG dataset was generated and subjected to 0.9-call rate filtering. Each SSG was selected when it met the segregation ratio of around 1:1 in the population.

2.5. Construction of the Genetic Linkage Map

We used a modified version of an approach designed for phase-unknown genetic linkage mapping in ants to construct a genetic linkage map using the phase-unknown SSGs, described in brief as follows [35]. In the first step, the complete SSGs array was doubled for each dataset, and the genotype score of every doubled SSG was switched (i.e., convert “A” to “B” and “B” to “A”) so that a new dataset containing every possible linkage phase for each SSG in each individual was generated. In the second step, linkage groups were generated using the linkage analysis program MSTmap [36]. Since the new dataset was doubled, each linkage group was represented twice. In step three, one of the two identical linkage groups was discarded. Step four involved the detection of genotyping errors and conversion to unknown phase, considering that double recombination is a very unlikely event among close markers. In the fifth and final step, the accuracy of the order of SSGs in each linkage group was checked with the R package ‘ASMap’ [37].

2.6. Comparative Analyses between the Linkage Map and FUGU5

The sequences of the SSGs located on the genetic linkage map obtained with the strategy outlined above were aligned to the published torafugu genome (FUGU5) using BLAST (version 2.2.29) with an e-value cut-off of 1 × 10−100 and identity of at least 95%. Genomic synteny was visualized using CIRCOS 0.69 software [38]. The flanking sequences of all mapped SNPs (35 bp on either side) were also subjected to a BLAST search against the published genome with an e-value cut-off of 1 × 10−26.

3. Results

3.1. De Novo Assembly

In total, 69.58 Gb of sequencing data from 165 samples of the generated DH torafugu population were used to perform de novo assembly. The total sequencing data coverage was 174, whereas the average coverage for each sample was approximately 1. After performing de novo assembly using the sequencing data, a relative low-quality assembly of a total size of 356.59 Mb and N50 size of 22,235 bp was generated, which was composed of 54,127 scaffolds with the length ranging from 200 to 264,568 bp.

3.2. Genetic Markers Genotyping

After SNP calling from the sequencing data of each DH individual, a total of 1,070,601 SNPs were discovered in the population using the above de novo assembly as reference. Despite the existence of a large quantity of missed genotypes and insufficient data for genotyping calibration due to the low-coverage (≈1×) sequencing of each sample, the genotypes of adjacent SNPs could be testified/compensated by each other based on the above de novo assembly. Therefore, as shown in Figure 1, the DNA sequences of the assembly were sliced into short segments. The information of SNPs located on each segment was combined and a genotype was assigned to each segment. This low-call-rate SNP dataset was then converted into a high-call-rate SSG dataset. After 0.9-call rate filtering, 37,398 SSGs containing information of 833,594 SNPs were retained and could be used for construction of the genetic linkage map (see details in Supplementary Materials).

3.3. Ultrahigh-Density Genetic Linkage Map Construction for Torafugu

An ultrahigh-density genetic linkage map (Figure 3) was constructed using the high-call-rate SSG dataset. As summarized in Table 1, the map consists of 37,343 SSGs in 3090 unique positions, containing the information of 802,277 SNPs (74.9% of total SNPs). The genetic linkage map contained 22 linkage groups, consistent with the number of chromosomes of the torafugu haploid genome. The genetic distances ranged from 62.75 cM of linkage group 10 (LG10) to 198.25 cM of LG1, with a total length of 2319.65 cM. Based on the unique marker positions, the estimated marker intervals ranged from 0.70 cM/marker in LG22 to 0.79 cM/marker in LG1, with an average marker interval of 0.75 cM/marker on the genetic linkage map. As shown in the heat map in Figure 4, the recombination fractions were considerably low between adjacent markers of each linkage group, indicating a low recombination frequency between them, whereas the LOD (logarithm of odds) scores between adjacent markers of each linkage group were high, indicating strong linkage between them. The 22 linkage groups appeared to be distinctly clustered.

3.4. Comparative Analyses

The sequence information of the 37,343 SSGs of the linkage map obtained with the proposed strategy was subjected to BLASTN analyses against the latest published genome FUGU5. Overall, 31,822 SSGs could be mapped to the 22 chromosomes of FUGU5. As shown in Figure 5, the Circos plot indicated near-perfect concordance between the genetic and physical position of each matched SSG, and 5521 SSGs could be mapped to the 1583 (65.97 Mb) unassembled scaffolds. The genetic positions of these SSGs are also highlighted in Figure 3, indicating the regions where FUGU5 can be improved. Furthermore, 180 of these scaffolds (28.2 Mb) contained more than one unique genetic position, suggesting that they might be located on the chromosomes with direction. The flanking sequences of the SNPs contained in this genetic linkage map were also well aligned to FUGU5, and 532,424 SNPs mapped to the 22 chromosomes of the genome. The plot shown in Figure 6 reflects more detailed collinearity between the orders of the SNPs of each linkage group and each chromosome. The results also indicated possible mis-assembled regions or large segmental polymorphisms in chromosome 2, 3, 4, 5, 6, 7, 11, 12, 17, 19, and 20 of FUGU5.

4. Discussion

We successfully developed an effective strategy for the construction of an ultrahigh-density genetic linkage map of torafugu based on low-coverage (≈1×) whole-genome sequencing of each individual of a DH population generated through mito-gynogenesis. The sequencing data were used for de novo assembly and further SNP calling to generate a low-call-rate SNP dataset with unknown phase. Based on the relatively low-quality de novo assembly, an SSG was designed as a high-call-rate genetic marker to assign a genotype to a short DNA segment after combing the information of its constituent low-call-rate SNPs. The high-call-rate SSG dataset enabled the construction of an ultrahigh-density genetic linkage map containing most of the information of SNPs (sub-million in this case) of the mapping population. The accuracy of the present linkage map was verified by subsequent analyses of recombination fractions and assessment of LOD scores for all marker pairs, along with comparative analyses between the linkage map and FUGU5. In addition, integration with the present linkage map allowed for validation and further refinement of FUGU5. Based on these indicators, an improved genome assembly of torafugu will be achieved in our future work.
In the whole-genome assembly project, de novo assembly and linkage map construction are independent works contributing to chromosome-scale assembly. However, in our case, both of them can be achieved from the low-coverage whole-genome sequencing of the mapping population. Thus, for a non-sequenced species, de novo assembly, linkage map construction, and further chromosome-scale assembly can be efficiently completed by adopting our strategy.
The strategy of the present study was developed for genetic linkage map construction based on an H/DH dataset with phase-unknown format. This strategy would be ideally implemented for various types of sexually reproducing organisms that could be used to generate large numbers of H/DH individuals, especially plants, teleosts, and fungi, without requiring other complex crossing schemes and designed inbred pedigree lines to identify the linkage phase. Notably, our strategy based on SSG could be applied to the construction of an ultrahigh-density genetic linkage map using only single gamete cells, the ubiquitously existing haploids, combined with single-cell sequencing technology, which further extends the application range. One of the most important advantages of our strategy is that it does not require the parental genetic phases and/or a high-quality reference genome, which are necessary for existing single-gamete sequencing strategies [39,40,41,42] to achieve high-quality genetic mapping.
The lack of a requirement of a high-quality reference genome for low-coverage whole-genome sequencing expands the application of our proposed strategy to a wide range of non-model and non-sequenced species. Moreover, our approach has an advantage of simplicity, in that whole-genome sequencing was applied to each individual sample, whereas existing techniques such as specific-locus amplified fragment sequencing (SLAF-seq) [43], reduced-representation libraries (RRLs) [44], and restriction-site-associated DNA sequencing (RAD-seq) [45] demand complicated processes for sequencing a small target portion of the whole genome. We successfully captured most of the SNPs of the population, which would allow for thoroughly characterizing complex genomes, whereas SLAF-seq, RAD-seq, and RRLs are only able to call a small portion of SNPs of the population. The low-coverage sequencing also makes the present strategy cost-effective. Collectively, these advantages demonstrate the potential of the present strategy as a good candidate for facilitating ultrahigh-density linkage map construction of various sexually reproducing organisms that can generate an H/DH population, and even can produce gametes combined with the single-cell sequencing technology platform.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/9/3/120/s1. Figure S1: Detailed information of the conversion of SNPs dataset to SSG dataset.

Acknowledgments

This work was supported by JSPS KAKENHI, Grant Numbers 15H02461, 24248034 and 25640099.

Author Contributions

S.A. conceive and designed the entire studies. X.Z. conceived and designed the data processing strategy and developed concept of SSG. S.A. and X.Z. wrote the manuscript. S.W. supervised generation of gynogenetic torafugu. H.Z., S.K., K.S. and S.A. generated the gynogenetic individuals. M.M., E.T., Y.I., Y.S., S.M. and S.A. contributed sequencing and SNP typing. All authors agreed on the final version of the submitted paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burt, D.W. Comparative mapping in farm animals. Brief Funct. Genomic. Proteomic. 2002, 1, 159–168. [Google Scholar] [CrossRef] [PubMed]
  2. Choi, H.K.; Mun, J.H.; Kim, D.J.; Zhu, H.; Baek, J.M.; Mudge, J.; Roe, B.; Ellis, N.; Doyle, J.; Kiss, G.B.; et al. Estimating genome conservation between crop and model legume species. Proc. Natl. Acad. Sci. USA 2004, 101, 15289–15294. [Google Scholar] [CrossRef] [PubMed]
  3. Naruse, K.; Fukamachi, S.; Mitani, H.; Kondo, M.; Matsuoka, T.; Kondo, S.; Hanamura, N.; Morita, Y.; Hasegawa, K.; Nishigaki, R.; et al. A detailed linkage map of medaka, Oryzias latipes: Comparative genomics and genome evolution. Genetics 2000, 154, 1773–1784. [Google Scholar] [PubMed]
  4. Hollenbeck, C.M.; Portnoy, D.S.; Gold, J.R. A genetic linkage map of red drum (Sciaenops ocellatus) and comparison of chromosomal syntenies with four other fish species. Aquaculture 2015, 435, 265–274. [Google Scholar] [CrossRef]
  5. Kujur, A.; Upadhyaya, H.D.; Shree, T.; Bajaj, D.; Das, S.; Saxena, M.S.; Badoni, S.; Kumar, V.; Tripathi, S.; Gowda, C.L.; et al. Ultra-high density intra-specific genetic linkage maps accelerate identification of functionally relevant molecular tags governing important agronomic traits in chickpea. Sci. Rep. 2015, 5, 9468. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Sun, C.; Niu, Y.; Ye, X.; Dong, J.; Hu, W.; Zeng, Q.; Chen, Z.; Tian, Y.; Zhang, J.; Lu, M. Construction of a high-density linkage map and mapping of sex determination and growth-related loci in the mandarin fish (Siniperca chuatsi). BMC Genomics 2017, 18, 446. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, J.; Li, L.T.; Li, M.; Khan, M.A.; Li, X.G.; Chen, H.; Yin, H.; Zhang, S.L. High-density genetic linkage map construction and identification of fruit-related QTL in pear using SNP and SSR markers. J. Exp. Bot. 2014, 65, 5771–5781. [Google Scholar] [CrossRef] [PubMed]
  8. Howe, K.; Clark, M.D.; Torroja, C.F.; Torrance, J.; Berthelot, C.; Muffato, M.; Collins, J.E.; Humphray, S.; McLaren, K.; Matthews, L.; et al. The zebrafish reference genome sequence and its relationship to the human genome. Nature 2013, 496, 498–503. [Google Scholar] [CrossRef] [PubMed]
  9. International Barley Genome Sequencing Consortium; Mayer, K.F.; Waugh, R.; Brown, J.W.; Schulman, A.; Langridge, P.; Platzer, M.; Fincher, G.B.; Muehlbauer, G.J.; Sato, K.; et al. A physical, genetic and functional sequence assembly of the barley genome. Nature 2012, 491, 711–716. [Google Scholar] [CrossRef] [PubMed]
  10. Shao, C.; Bao, B.; Xie, Z.; Chen, X.; Li, B.; Jia, X.; Yao, Q.; Orti, G.; Li, W.; Li, X.; et al. The genome and transcriptome of Japanese flounder provide insights into flatfish asymmetry. Nat. Genet. 2017, 49, 119–124. [Google Scholar] [CrossRef] [PubMed]
  11. Wang, Y.; Lu, Y.; Zhang, Y.; Ning, Z.; Li, Y.; Zhao, Q.; Lu, H.; Huang, R.; Xia, X.; Feng, Q.; et al. The draft genome of the grass carp (Ctenopharyngodon idellus) provides insights into its evolution and vegetarian adaptation. Nat. Genet. 2015, 47, 625–631. [Google Scholar] [CrossRef] [PubMed]
  12. Davey, J.W.; Hohenlohe, P.A.; Etter, P.D.; Boone, J.Q.; Catchen, J.M.; Blaxter, M.L. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat. Rev. Genet. 2011, 12, 499–510. [Google Scholar] [CrossRef] [PubMed]
  13. Xie, W.; Feng, Q.; Yu, H.; Huang, X.; Zhao, Q.; Xing, Y.; Yu, S.; Han, B.; Zhang, Q. Parent-independent genotyping for constructing an ultrahigh-density linkage map based on population sequencing. Proc. Natl. Acad. Sci. USA 2010, 107, 10578–10583. [Google Scholar] [CrossRef] [PubMed]
  14. Au, C.H.; Cheung, M.K.; Wong, M.C.; Chu, A.K.; Law, P.T.; Kwan, H.S. Rapid genotyping by low-coverage resequencing to construct genetic linkage maps of fungi: A case study in Lentinula edodes. BMC Res. Notes 2013, 6, 307. [Google Scholar] [CrossRef] [PubMed]
  15. Bowers, J.E.; Pearl, S.A.; Burke, J.M. Genetic mapping of millions of SNPs in safflower (Carthamus tinctorius L.) via whole-genome resequencing. G3-Genes Genom. Genet. 2016, 6, 2203–2211. [Google Scholar] [CrossRef] [PubMed]
  16. Lefebvre, V.; Palloix, A.; Caranta, C.; Pochard, E. Construction of an intraspecific integrated linkage map of pepper using molecular markers and doubled haploid progenies. Genome 1995, 38, 112–121. [Google Scholar] [CrossRef] [PubMed]
  17. Semagn, K.; Bjornstad, A.; Skinnes, H.; Maroy, A.G.; Tarkegne, Y.; William, M. Distribution of DArT, AFLP, and SSR markers in a genetic linkage map of a doubled-haploid hexaploid wheat population. Genome 2006, 49, 545–555. [Google Scholar] [CrossRef] [PubMed]
  18. Zhang, J.; Guo, W.; Zhang, T. Molecular linkage map of allotetraploid cotton (Gossypium hirsutum L. x Gossypium barbadense L.) with a haploid population. Theor. Appl. Genet. 2002, 105, 1166–1174. [Google Scholar] [PubMed]
  19. Xu, Y.; Huang, L.; Ji, D.; Chen, C.; Zheng, H.; Xie, C. Construction of a dense genetic linkage map and mapping quantitative trait loci for economic traits of a doubled haploid population of Pyropia haitanensis (Bangiales, Rhodophyta). BMC Plant Biol. 2015, 15, 228. [Google Scholar] [CrossRef] [PubMed]
  20. Delourme, R.; Falentin, C.; Fomeju, B.F.; Boillot, M.; Lassalle, G.; Andre, I.; Duarte, J.; Gauthier, V.; Lucante, N.; Marty, A.; et al. High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L. BMC Genomics 2013, 14, 120. [Google Scholar] [CrossRef] [PubMed]
  21. Young, W.P.; Wheeler, P.A.; Coryell, V.H.; Keim, P.; Thorgaard, G.H. A detailed linkage map of rainbow trout produced using doubled haploids. Genetics 1998, 148, 839–850. [Google Scholar] [PubMed]
  22. Brieuc, M.S.; Waters, C.D.; Seeb, J.E.; Naish, K.A. A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event. G3-Genes Genom. Genet. 2014, 4, 447–460. [Google Scholar] [CrossRef] [PubMed]
  23. Kelly, P.D.; Chu, F.; Woods, I.G.; Ngo-Hazelett, P.; Cardozo, T.; Huang, H.; Kimm, F.; Liao, L.Y.; Yan, Y.L.; Zhou, Y.Y.; et al. Genetic linkage mapping of zebrafish genes and ESTs. Genome Res. 2000, 10, 558–567. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, X.W.; Liang, L.Q. A genetic linkage map of common carp (Cyprinus carpio L.) and mapping of a locus associated with cold tolerance. Aquaculture 2004, 238, 165–172. [Google Scholar] [CrossRef]
  25. Komen, H.; Thorgaard, G.H. Androgenesis, gynogenesis and the production of clones in fishes: A review. Aquaculture 2007, 269, 150–173. [Google Scholar] [CrossRef]
  26. Dunwell, J.M. Haploids in flowering plants: Origins and exploitation. Plant Biotechnol. J. 2010, 8, 377–424. [Google Scholar] [CrossRef] [PubMed]
  27. Brenner, S.; Elgar, G.; Sandford, R.; Macrae, A.; Venkatesh, B.; Aparicio, S. Characterization of the pufferfish (Fugu) genome as a compact model vertebrate genome. Nature 1993, 366, 265–268. [Google Scholar] [CrossRef] [PubMed]
  28. Grutzner, F.; Lutjens, G.; Rovira, C.; Barnes, D.W.; Ropers, H.H.; Haaf, T. Classical and molecular cytogenetics of the pufferfish Tetraodon nigroviridis. Chromosome Res. 1999, 7, 655–662. [Google Scholar] [CrossRef] [PubMed]
  29. Hedges, S.B.; Kumar, S. Vertebrate genomes compared. Science 2002, 297, 1283–1285. [Google Scholar] [CrossRef] [PubMed]
  30. Kai, W.; Kikuchi, K.; Tohari, S.; Chew, A.K.; Tay, A.; Fujiwara, A.; Hosoya, S.; Suetake, H.; Naruse, K.; Brenner, S.; et al. Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals. Genome Biol. Evol. 2011, 3, 424–442. [Google Scholar] [CrossRef] [PubMed]
  31. Zhang, H.; Hirose, Y.; Yoshino, R.; Kinoshita, S.; Saito, K.; Tan, E.K.; Suzuki, Y.; Shimizu, A.; Okano, H.; Kudoh, J.; et al. Assessment of homozygosity levels in the mito-gynogenetic torafugu (Takifugu rubripes) by genome-wide SNP analyses. Aquaculture 2013, 380, 114–119. [Google Scholar] [CrossRef]
  32. Luo, R.B.; Liu, B.H.; Xie, Y.L.; Li, Z.Y.; Huang, W.H.; Yuan, J.Y.; He, G.Z.; Chen, Y.X.; Pan, Q.; Liu, Y.J.; et al. SOAPdenovo2: An empirically improved memory-efficient short-read de novo assembler. Gigascience 2012, 1. [Google Scholar] [CrossRef] [PubMed]
  33. Li, H.; Durbin, R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed]
  34. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef] [PubMed]
  35. Gadau, J. Phase-unknown linkage mapping in ants. Cold Spring Harb. Protoc. 2009, pdb.prot5251. [Google Scholar] [CrossRef] [PubMed]
  36. Wu, Y.H.; Bhat, P.R.; Close, T.J.; Lonardi, S. Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008, 4, e1000212. [Google Scholar] [CrossRef] [PubMed]
  37. Taylor, J.; Butler, D. R package ASMap: Efficient genetic linkage map construction and diagnosis. J. Stat. Softw. 2017, 79, 1–29. [Google Scholar] [CrossRef]
  38. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [PubMed]
  39. Lu, S.J.; Zong, C.H.; Fan, W.; Yang, M.Y.; Li, J.S.; Chapman, A.R.; Zhu, P.; Hu, X.S.; Xu, L.Y.; Yan, L.Y.; et al. Probing meiotic recombination and aneuploidy of single sperm cells by whole-genome sequencing. Science 2012, 338, 1627–1630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Xu, S.; Ackerman, M.S.; Long, H.A.; Bright, L.; Spitze, K.; Ramsdell, J.S.; Thomas, W.K.; Lynch, M. A male-specific genetic map of the microcrustacean Daphnia pulex based on single-sperm whole-genome sequencing. Genetics 2015, 201, 31–38. [Google Scholar] [CrossRef] [PubMed]
  41. Ottolini, C.S.; Newnham, L.J.; Capalbo, A.; Natesan, S.A.; Joshi, H.A.; Cimadomo, D.; Griffin, D.K.; Sage, K.; Summers, M.C.; Thornhill, A.R.; et al. Genome-wide maps of recombination and chromosome segregation in human oocytes and embryos show selection for maternal recombination rates. Nat. Genet. 2015, 47, 727–735. [Google Scholar] [CrossRef] [PubMed]
  42. Ottolini, C.S.; Capalbo, A.; Newnham, L.; Cimadomo, D.; Natesan, S.A.; Hoffmann, E.R.; Ubaldi, F.M.; Rienzi, L.; Handyside, A.H. Generation of meiomaps of genome-wide recombination and chromosome segregation in human oocytes. Nat. Protoc. 2016, 11, 1229–1243. [Google Scholar] [CrossRef] [PubMed]
  43. Sun, X.; Liu, D.; Zhang, X.; Li, W.; Liu, H.; Hong, W.; Jiang, C.; Guan, N.; Ma, C.; Zeng, H.; et al. SLAF-seq: An efficient method of large-scale de novo SNP discovery and genotyping using high-throughput sequencing. PLoS ONE 2013, 8, e58700. [Google Scholar] [CrossRef] [PubMed]
  44. Altshuler, D.; Pollara, V.J.; Cowles, C.R.; Van Etten, W.J.; Baldwin, J.; Linton, L.; Lander, E.S. An SNP map of the human genome generated by reduced representation shotgun sequencing. Nature 2000, 407, 513–516. [Google Scholar] [PubMed]
  45. Baird, N.A.; Etter, P.D.; Atwood, T.S.; Currey, M.C.; Shiver, A.L.; Lewis, Z.A.; Selker, E.U.; Cresko, W.A.; Johnson, E.A. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 2008, 3, e3376. [Google Scholar] [CrossRef] [PubMed]
Genes 09 00120 g001 550
Figure 1. Flowchart of high-call-rate markers genotyping. (a) Low-coverage whole-genome sequencing of the doubled-haploid (DH) mapping population was performed to generate a library of short gun reads (in blue) for each sample (S1‒Sn). (b) De novo assembly was performed using the total whole-genome sequencing data to generate scaffolds (in red) containing a series of segments (in green). (c) Single nucleotide polymorphisms (SNP) calling and genotyping of each sample was carried out to construct a low-call-rate SNPs array (SNPs × Samples) with unknown phase. The SNPs (SNP1‒5) located on one short segment (Seg1) were phased at the segment level and merged together to assign the short segment genotype (SSG1). The high-call-rate SSGs array (SSGs × Samples) was then constructed with unknown phase at the chromatid level.
Figure 1. Flowchart of high-call-rate markers genotyping. (a) Low-coverage whole-genome sequencing of the doubled-haploid (DH) mapping population was performed to generate a library of short gun reads (in blue) for each sample (S1‒Sn). (b) De novo assembly was performed using the total whole-genome sequencing data to generate scaffolds (in red) containing a series of segments (in green). (c) Single nucleotide polymorphisms (SNP) calling and genotyping of each sample was carried out to construct a low-call-rate SNPs array (SNPs × Samples) with unknown phase. The SNPs (SNP1‒5) located on one short segment (Seg1) were phased at the segment level and merged together to assign the short segment genotype (SSG1). The high-call-rate SSGs array (SSGs × Samples) was then constructed with unknown phase at the chromatid level.
Genes 09 00120 g001
Genes 09 00120 g002 550
Figure 2. Schematic illustration of the conversion of SSGs under different circumstances. The maternal chromatids possess linkage phase score A (blue) and B (red). In Sample 1, recombination occurs on Seg2 but outside of Seg1. In Sample 2, recombination occurs outside of both Seg1 and Seg2, which is the most common case.
Figure 2. Schematic illustration of the conversion of SSGs under different circumstances. The maternal chromatids possess linkage phase score A (blue) and B (red). In Sample 1, recombination occurs on Seg2 but outside of Seg1. In Sample 2, recombination occurs outside of both Seg1 and Seg2, which is the most common case.
Genes 09 00120 g002
Genes 09 00120 g003 550
Figure 3. Ideograms of the genetic linkage map of torafugu. The genetic position of each SSG is illustrated in each linkage group. SSGs are indicated as a blue line if the sequences could be aligned to the published genome FUGU5 using BLASTN (version 2.2.29), and are otherwise shown as a pink line. LG: linkage group; the y-axis represents the genetic position (cM).
Figure 3. Ideograms of the genetic linkage map of torafugu. The genetic position of each SSG is illustrated in each linkage group. SSGs are indicated as a blue line if the sequences could be aligned to the published genome FUGU5 using BLASTN (version 2.2.29), and are otherwise shown as a pink line. LG: linkage group; the y-axis represents the genetic position (cM).
Genes 09 00120 g003
Genes 09 00120 g004 550
Figure 4. Heat map of recombination fractions (upper triangle) and logarithm of odds (LOD) scores (lower triangle) for all pairs of markers (SSGs) in the 22 linkage groups. LG represents linkage group.
Figure 4. Heat map of recombination fractions (upper triangle) and logarithm of odds (LOD) scores (lower triangle) for all pairs of markers (SSGs) in the 22 linkage groups. LG represents linkage group.
Genes 09 00120 g004
Genes 09 00120 g005 550
Figure 5. Comparison between the present genetic linkage map and FUGU5. Links connect the genetic and physical locations of each mapped SSG. Grey blocks on the left represent the 22 linkage groups of the present linkage map, whereas colorful blocks on the right represent the 22 chromosomes of FUGU5. LG: linkage group; chr: chromosome.
Figure 5. Comparison between the present genetic linkage map and FUGU5. Links connect the genetic and physical locations of each mapped SSG. Grey blocks on the left represent the 22 linkage groups of the present linkage map, whereas colorful blocks on the right represent the 22 chromosomes of FUGU5. LG: linkage group; chr: chromosome.
Genes 09 00120 g005
Genes 09 00120 g006 550
Figure 6. Concordance of SNP positions on the present genetic linkage map with those on FUGU5. The x-axis represents the physical position (Mb) of SNPs on FUGU5, whereas the y-axis represents the genetic position (cM) of SNPs on the linkage map. Red arrows indicate potentially misplaced or segmentally polymorphic regions, whereas green arrows indicate potentially mis-oriented or segmentally polymorphic regions in FUGU5.
Figure 6. Concordance of SNP positions on the present genetic linkage map with those on FUGU5. The x-axis represents the physical position (Mb) of SNPs on FUGU5, whereas the y-axis represents the genetic position (cM) of SNPs on the linkage map. Red arrows indicate potentially misplaced or segmentally polymorphic regions, whereas green arrows indicate potentially mis-oriented or segmentally polymorphic regions in FUGU5.
Genes 09 00120 g006
Table
Table 1. Summary of the genetic linkage map of torafugu.
Table 1. Summary of the genetic linkage map of torafugu.
Linkage GroupGenetic Distance (cM)SSGs NumbersSNPs NumbersUnique PositionsMarker Interval
LG1198.25299764,1192520.79
LG296.95158834,3691280.76
LG3112.13187639,2801470.76
LG495.15172235,5351320.72
LG596.42145431,0431360.71
LG681.72143131,6161100.74
LG7128.42181937,2251660.77
LG8132.44205544,3021740.76
LG9104.88167034,9751430.73
LG1062.75139630,320820.77
LG11104.43154333,1861430.73
LG1279.40140830,6361100.72
LG13122.28179738,7391570.78
LG1497.47177037,7501280.76
LG1592.43160535,5581190.78
LG1664.71141730,875890.73
LG17110.91172738,6201500.74
LG1867.20115726,774890.76
LG19112.95171435,3121540.73
LG20120.46183038,9891520.79
LG21133.31189340,9751800.74
LG22104.98147432,0791490.70
Total2319.6537,343802,27730900.75 1
1 Mean value of marker intervals.

Share and Cite

MDPI and ACS Style

Zhang, X.; Mizukoshi, M.; Zhang, H.; Tan, E.; Igarashi, Y.; Suzuki, Y.; Mitsuyama, S.; Kinoshita, S.; Saito, K.; Watabe, S.; et al. Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes). Genes 2018, 9, 120. https://doi.org/10.3390/genes9030120

AMA Style

Zhang X, Mizukoshi M, Zhang H, Tan E, Igarashi Y, Suzuki Y, Mitsuyama S, Kinoshita S, Saito K, Watabe S, et al. Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes). Genes. 2018; 9(3):120. https://doi.org/10.3390/genes9030120

Chicago/Turabian Style

Zhang, Xiang, Misaki Mizukoshi, Hong Zhang, Engkong Tan, Yoji Igarashi, Yutaka Suzuki, Susumu Mitsuyama, Shigeharu Kinoshita, Kazuyoshi Saito, Shugo Watabe, and et al. 2018. "Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes)" Genes 9, no. 3: 120. https://doi.org/10.3390/genes9030120

APA Style

Zhang, X., Mizukoshi, M., Zhang, H., Tan, E., Igarashi, Y., Suzuki, Y., Mitsuyama, S., Kinoshita, S., Saito, K., Watabe, S., & Asakawa, S. (2018). Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes). Genes, 9(3), 120. https://doi.org/10.3390/genes9030120

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop