Construction of integrated genetic linkage maps of the tiger shrimp (Penaeus monodon) using microsatellite and AFLP markers

The linkage maps of male and female tiger shrimp (P. monodon) were constructed based on 256 microsatellite and 85 amplified fragment length polymorphism (AFLP) markers. Microsatellite markers obtained from clone sequences of partial genomic libraries, tandem repeat sequences from databases and previous publications and fosmid end sequences were employed. Of 670 microsatellite and 158 AFLP markers tested for polymorphism, 341 (256 microsatellite and 85 AFLP markers) were used for genotyping with three F₁ mapping panels, each comprising two parents and more than 100 progeny. Chi‐square goodness‐of‐fit test (χ²) revealed that only 19 microsatellite and 28 AFLP markers showed a highly significant segregation distortion (P < 0.005). Linkage analysis with a LOD score of 4.5 revealed 43 and 46 linkage groups in male and female linkage maps respectively. The male map consisted of 176 microsatellite and 49 AFLP markers spaced every ～11.2 cM, with an observed genome length of 2033.4 cM. The female map consisted of 171 microsatellite and 36 AFLP markers spaced every ～13.8 cM, with an observed genome length of 2182 cM. Both maps shared 136 microsatellite markers, and the alignment between them indicated 38 homologous pairs of linkage groups including the linkage group representing the sex chromosome. The karyotype of P. monodon is also presented. The tentative assignment of the 44 pairs of P. monodon haploid chromosomes showed the composition of forty metacentric, one submetacentric and three acrocentric chromosomes. Our maps provided a solid foundation for gene and QTL mapping in the tiger shrimp.     

Construction of a high-density,high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.     

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.     

Construction of Genetic Linkage Maps and Comparative Genome Analysis of Catfish Using Gene-Associated Markers

A genetic linkage map of the channel catfish genome (N = 29) was constructed using EST-based microsatellite and single nucleotide polymorphism (SNP) markers in an interspecific reference family. A total of 413 microsatellites and 125 SNP markers were polymorphic in the reference family. Linkage analysis using JoinMap 4.0 allowed mapping of 331 markers (259 microsatellites and 72 SNPs) to 29 linkage groups. Each linkage group contained 3–18 markers. The largest linkage group contained 18 markers and spanned 131.2 cM, while the smallest linkage group contained 14 markers and spanned only 7.9 cM. The linkage map covered a genetic distance of 1811 cM with an average marker interval of 6.0 cM. Sex-specific maps were also constructed; the recombination rate for females was 1.6 times higher than that for males. Putative conserved syntenies between catfish and zebrafish, medaka, and Tetraodon were established, but the overall levels of genome rearrangements were high among the teleost genomes. This study represents a first-generation linkage map constructed by using EST-derived microsatellites and SNPs, laying a framework for large-scale comparative genome analysis in catfish. The conserved syntenies identified here between the catfish and the three model fish species should facilitate structural genome analysis and evolutionary studies, but more importantly should facilitate functional inference of catfish genes. Given that determination of gene functions is difficult in nonmodel species such as catfish, functional genome analysis will have to rely heavily on the establishment of orthologies from model species.     

Bovine chromosome 4 workshop: consensus and comprehensive linkage maps

The report of the bovine chromosome 4 (BTA4) workshop is presented. Six laboratories contributed a total of 30,168 informative meioses from 62 loci. Twenty-two loci were typed by at least two independent laboratories and were used to construct a consensus linkage map of BTA4. The remaining 40 loci were subsequently incorporated into a comprehensive map. The sex-averaged consensus map covered 131.4 cM. The female map was 124.3 cM in length, while the male map was 134.3 cM. The comprehensive sex-averaged map spanned 141.6 cM. The length of the female and male comprehensive maps were 123.1 cM and 156.4 cM, respectively. Average genetic distance between loci was 6 and 2.3 cM for the consensus and comprehensive linkage maps, respectively.     

Genetic linkage maps of barfin flounder (<Emphasis Type="Italic">Verasper moseri</Emphasis>) and spotted halibut (<Emphasis Type="Italic">Verasper variegatus</Emphasis>) based on AFLP and microsatellite markers

Barfin flounder (Verasper moseri) and spotted halibut (Verasper variegatus) are two economically important marine fish species for aquaculture in China, Korea and Japan. Construction of genetic linkage maps is an interesting issue for molecular marker-assisted selection (MAS) and for better understanding the genome structure. In the present study, we constructed genetic linkage maps for both fish species using AFLP and microsatellite markers based on an interspecific F₁ hybrid family (female V. moseri and male V. variegatus). The female genetic map comprised 98 markers (58 AFLP markers and 40 microsatellite markers), distributing in 27 linkage groups, and spanning 637 cM with an average resolution of 8.9 cM. Whereas the male genetic map consisted of 86 markers (48 AFLP and 38 microsatellite markers) in 24 linkage groups, covering a length of 625 cM with an average marker spacing of 10 cM. The expected genome length was 1,128 cM in female and 1,115 cM in male, and the estimated coverage of genome was 56% for both genetic maps. Moreover, five microsatellite markers were observed to be common to both genetic maps. This is the first time to report the genetic linkage maps of V. moseri and V. variegatus that could serve as the basis for genetic improvement and selective breeding, candidate genes cloning, and genome structure research.     

Identification of Sex-determining Loci in Pacific White Shrimp <Emphasis Type="Italic">Litopeneaus vannamei</Emphasis> Using Linkage and Association Analysis

The Pacific white shrimp Litopenaeus vannamei is a predominant aquaculture shrimp species in the world. Like other animals, the L. vannamei exhibited sexual dimorphism in growth trait. Mapping of the sex-determining locus will be very helpful to clarify the sex determination system and further benefit the shrimp aquaculture industry towards the production of mono-sex stocks. Based on the data used for high-density linkage map construction, linkage-mapping analysis was conducted. The sex determination region was mapped in linkage group (LG) 18. A large region from 0 to 21.205 cM in LG18 showed significant association with sex. However, none of the markers in this region showed complete association with sex in the other populations. So an association analysis was designed using the female parent, pool of female progenies, male parent, and pool of male progenies. Markers were de novo developed and those showing significant differences between female and male pools were identified. Among them, three sex-associated markers including one fully associated marker were identified. Integration of linkage and association analysis showed that the sex determination region was fine-mapped in a small region along LG18. The identified sex-associated marker can be used for the sex detection of this species at genetic level. The fine-mapped sex-determining region will contribute to the mapping of sex-determining gene and help to clarify sex determination system for L. vannamei.     

A Genome-Wide Search for Linkage of Estimated Glomerular Filtration Rate (eGFR) in the Family Investigation of Nephropathy and Diabetes (FIND)

Objective

Estimated glomerular filtration rate (eGFR), a measure of kidney function, is heritable, suggesting that genes influence renal function. Genes that influence eGFR have been identified through genome-wide association studies. However, family-based linkage approaches may identify loci that explain a larger proportion of the heritability. This study used genome-wide linkage and association scans to identify quantitative trait loci (QTL) that influence eGFR.

Methods

Genome-wide linkage and sparse association scans of eGFR were performed in families ascertained by probands with advanced diabetic nephropathy (DN) from the multi-ethnic Family Investigation of Nephropathy and Diabetes (FIND) study. This study included 954 African Americans (AA), 781 American Indians (AI), 614 European Americans (EA) and 1,611 Mexican Americans (MA). A total of 3,960 FIND participants were genotyped for 6,000 single nucleotide polymorphisms (SNPs) using the Illumina Linkage IVb panel. GFR was estimated by the Modification of Diet in Renal Disease (MDRD) formula.

Results

The non-parametric linkage analysis, accounting for the effects of diabetes duration and BMI, identified the strongest evidence for linkage of eGFR on chromosome 20q11 (log of the odds [LOD] = 3.34; P = 4.4×10⁻⁵) in MA and chromosome 15q12 (LOD = 2.84; P = 1.5×10⁻⁴) in EA. In all subjects, the strongest linkage signal for eGFR was detected on chromosome 10p12 (P = 5.5×10⁻⁴) at 44 cM near marker rs1339048. A subsequent association scan in both ancestry-specific groups and the entire population identified several SNPs significantly associated with eGFR across the genome.

Conclusion

The present study describes the localization of QTL influencing eGFR on 20q11 in MA, 15q21 in EA and 10p12 in the combined ethnic groups participating in the FIND study. Identification of causal genes/variants influencing eGFR, within these linkage and association loci, will open new avenues for functional analyses and development of novel diagnostic markers for DN.     

A Preliminary Genetic Linkage Map of the Pacific Abalone Haliotis discus hannai Ino

Preliminary genetic linkage maps were constructed for the Pacific abalone (Haliotis discus hannai Ino) using amplified fragment length polymorphism (AFLP), randomly amplified polymorphic DNA (RAPD), and microsatellite markers segregating in a F₁ family. Nine microsatellite loci, 41 RAPD, and 2688 AFLP markers were genotyped in the parents and 86 progeny of the mapping family. Among the 2738 markers, 384 (including 365 AFLP markers, 10 RAPD markers, and 9 microsatellite loci) were polymorphic and segregated in one or both parents: 241 in the female and 146 in the male. The majority of these markers, 232 in the female and 134 in the male, segregated according to the expected 1:1 Mendelian ratio (α = 0.05). Two genetic linkage maps were constructed using markers segregating in the female or the male parent. The female framework map consisted of 119 markers in 22 linkage groups, covering 1773.6 cM with an average intermarker space of 18.3 cM. The male framework map contained 94 markers in 19 linkage groups, spanning 1365.9 cM with an average intermarker space of 18.2 cM. The sex determination locus was mapped to the male map but not to the female map, suggesting a XY-male determination mechanism. Distorted markers showing excess of homozygotes were mapped in clusters, probably because of their linkage to a gene that is incompatible between two parental populations.     

Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori, based on bacterial artificial chromosome end sequences

We have developed a linkage map for the silkworm Bombyx mori based on single nucleotide polymorphisms (SNPs) between strains p50T and C108T initially found on regions corresponding to the end sequences of bacterial artificial chromosome (BAC) clones. Using 190 segregants from a backcross of a p50T female x an F1 (p50T x C108T) male, we analyzed segregation patterns of 534 SNPs between p50T and C108T, detected among 3840 PCR amplicons, each associated with a p50T BAC end sequence. This enabled us to construct a linkage map composed of 534 SNP markers spanning 1305 cM in total length distributed over the expected 28 linkage groups. Of the 534 BACs whose ends harbored the SNPs used to construct the linkage map, 89 were associated with 107 different ESTs. Since each of the SNP markers is directly linked to a specific genomic BAC clone and to whole-genome sequence data, and some of them are also linked to EST data, the SNP linkage map will be a powerful tool for investigating silkworm genome properties, mutation mapping, and map-based cloning of genes of industrial and agricultural interest.     

Development and Evaluation of SoySNP50K,a High-Density Genotyping Array for Soybean

The objective of this research was to identify single nucleotide polymorphisms (SNPs) and to develop an Illumina Infinium BeadChip that contained over 50,000 SNPs from soybean (Glycine max L. Merr.). A total of 498,921,777 reads 35–45bp in length were obtained from DNA sequence analysis of reduced representation libraries from several soybean accessions which included six cultivated and two wild soybean (G. soja Sieb. et Zucc.) genotypes. These reads were mapped to the soybean whole genome sequence and 209,903 SNPs were identified. After applying several filters, a total of 146,161 of the 209,903 SNPs were determined to be ideal candidates for Illumina Infinium II BeadChip design. To equalize the distance between selected SNPs, increase assay success rate, and minimize the number of SNPs with low minor allele frequency, an iteration algorithm based on a selection index was developed and used to select 60,800 SNPs for Infinium BeadChip design. Of the 60,800 SNPs, 50,701 were targeted to euchromatic regions and 10,000 to heterochromatic regions of the 20 soybean chromosomes. In addition, 99 SNPs were targeted to unanchored sequence scaffolds. Of the 60,800 SNPs, a total of 52,041 passed Illumina’s manufacturing phase to produce the SoySNP50K iSelect BeadChip. Validation of the SoySNP50K chip with 96 landrace genotypes, 96 elite cultivars and 96 wild soybean accessions showed that 47,337 SNPs were polymorphic and generated successful SNP allele calls. In addition, 40,841 of the 47,337 SNPs (86%) had minor allele frequencies ≥10% among the landraces, elite cultivars and the wild soybean accessions. A total of 620 and 42 candidate regions which may be associated with domestication and recent selection were identified, respectively. The SoySNP50K iSelect SNP beadchip will be a powerful tool for characterizing soybean genetic diversity and linkage disequilibrium, and for constructing high resolution linkage maps to improve the soybean whole genome sequence assembly.     

Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome

Here, we present an adaptation of restriction‐site‐associated DNA sequencing (RAD‐seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white‐eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18–25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired‐end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost‐effective way.     

Comparison of microsatellites, single-nucleotide polymorphisms (SNPs) and composite markers derived from SNPs in linkage analysis

There is growing evidence that a map of dense single-nucleotide polymorphisms (SNPs) can outperform a map of sparse microsatellites for linkage analysis. There is also argument as to whether a clustered SNP map can outperform an evenly spaced SNP map. Using Genetic Analysis Workshop 14 simulated data, we compared for linkage analysis microsatellites, SNPs, and composite markers derived from SNPs. We encoded the composite markers in a two-step approach, in which the maximum identity length contrast method was employed to allow for recombination between loci. A SNP map 2.3 times as dense as a microsatellite map (approximately 2.9 cM compared to approximately 6.7 cM apart) provided slightly less information content (approximately 0.83 compared to approximately 0.89). Most inheritance information could be extracted when the SNPs were spaced < 1 cM apart. Comparing the linkage results on using SNPs or composite markers derived from them based on both 3 cM and 0.3 cM resolution maps, we showed that the inter-SNP distance should be kept small (< 1 cM), and that for multipoint linkage analysis the original markers and the derived composite markers had similar power; but for single point linkage analysis the resulting composite markers lead to more power. Considering all factors, such as information content, flexibility of analysis method, map errors, and genotyping errors, a map of clustered SNPs can be an efficient design for a genome-wide linkage scan.     

Construction of a genetic linkage map in Fenneropenaeus chinensis using SNP markers

Genetic linkage maps of Fenneropenaeus chinensis were constructed using a “double pseudo-testcross” strategy with 200 single nucleotide polymorphisms (SNPs) markers. This study represents the first SNP genetic linkage map for F. chinensis. The parents and F ₁ progeny of 100 individuals were used as mapping populations. 21 genetic linkage groups in the male and female maps were identified. The male linkage map was composed of 115 loci and spanned 879.7 cM, with an average intermarker spacing of 9.4 cM, while the female map was composed of 119 loci and spanned 876.2 cM, with an average intermarker spacing of 8.9 cM. The estimated coverage of the linkage maps was 51.94% for the male and 53.77% for the female, based on two estimates of genome length. The integrated map contains 180 markers distributed in 16 linkage groups, and spans 899.3 cM with an average marker interval of 5.2 cM. This SNP genetic map lays the foundation for future shrimp genomics and genetic breeding studies, especially the discovery of gene or regions for economically important traits in Chinese shrimp.     

High-density linkage maps and sex-linked markers for the black tiger shrimp (Penaeus monodon)

We report on the construction of sex-specific high-density linkage maps and identification of sex-linked markers for the black tiger shrimp (Penaeus monodon). Overall, we identified 44 male and 43 female linkage groups (2n = 88) from the analysis of 2,306 AFLP markers segregating in three full-sib families, covering 2,378 and 2,362 cM, respectively. Twenty-one putatively homologous linkage groups, including the sex-linkage groups, were identified between the female and male linkage maps. Six sex-linked AFLP marker alleles were inherited from female parents in the three families, suggesting that the P. monodon adopts a WZ-ZZ sex-determining system. Two sex-linked AFLP markers, one of which we converted into an allele-specific assay, confirmed their association with sex in a panel of 52 genetically unrelated animals.     

Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus

Background

Kelp (Saccharina japonica) has been intensively cultured in China for almost a century. Its genetic improvement is comparable with that of rice. However, the development of its molecular tools is extremely limited, thus its genes, genetics and genomics. Kelp performs an alternative life cycle during which sporophyte generation alternates with gametophyte generation. The gametophytes of kelp can be cloned and crossed. Due to these characteristics, kelp may serve as a reference for the biological and genetic studies of Volvox, mosses and ferns.

Results

We constructed a high density single nucleotide polymorphism (SNP) linkage map for kelp by restriction site associated DNA (RAD) sequencing. In total, 4,994 SNP-containing physical (tag-defined) RAD loci were mapped on 31 linkage groups. The map expanded a total genetic distance of 1,782.75 cM, covering 98.66% of the expected (1,806.94 cM). The length of RAD tags (85 bp) was extended to 400–500 bp with Miseq method, offering us an easiness of developing SNP chips and shifting SNP genotyping to a high throughput track. The number of linkage groups was in accordance with the documented with cytological methods. In addition, we identified a set of microsatellites (99 in total) from the extended RAD tags. A gametophyte sex determining locus was mapped on linkage group 2 in a window about 9.0 cM in width, which was 2.66 cM up to marker_40567 and 6.42 cM down to marker_23595.

Conclusions

A high density SNP linkage map was constructed for kelp, an intensively cultured brown alga in China. The RAD tags were also extended so that a SNP chip could be developed. In addition, a set of microsatellites were identified among mapped loci, and a gametophyte sex determining locus was mapped. This map will facilitate the genetic studies of kelp including for example the evaluation of germplasm and the decipherment of the genetic bases of economic traits.

Electronic supplementary material

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

Construction of a High-Density Microsatellite Genetic Linkage Map and Mapping of Sexual and Growth-Related Traits in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

High-density genetic linkage maps of half-smooth tongue sole were developed with 1007 microsatellite markers, two SCAR markers and an F1 family containing 94. The female map was composed of 828 markers in 21 linkage groups, covering a total of 1447.3 cM, with an average interval 1.83 cM between markers. The male map consisted of 794 markers in 21 linkage groups, spanning 1497.5 cM, with an average interval of 1.96 cM. The female and male maps had 812 and 785 unique positions, respectively. The genome length of half-smooth tongue sole was estimated to be 1527.7 cM for the females and 1582.1 cM for the males. Based on estimations of the map lengths, the female and male maps covered 94.74 and 94.65% of the genome, respectively. The consensus map was composed of 1007 microsatellite markers and two SCAR markers in 21 linkage groups, covering a total of 1624 cM with an average interval of 1.67 cM. Furthermore, 159 sex-linked SSR markers were identified. Five sex-linked microsatellite markers were confirmed in their association with sex in a large number of individuals selected from different families. These sex-linked markers were mapped on the female map LG1f with zero recombination. Two QTLs that were identified for body weight, designated as We-1 and We-2, accounted for 26.39% and 10.60% of the phenotypic variation. Two QTLs for body width, designated Wi-1 and Wi-2, were mapped in LG4f and accounted for 14.33% and 12.83% of the phenotypic variation, respectively. Seven sex-related loci were mapped in LG1f, LG14f and LG1m by CIM, accounting for 12.5–25.2% of the trait variation. The results should prove to be very useful for improving growth traits using molecular MAS.     

Linkage maps for the Pacific abalone (genus Haliotis) based on microsatellite DNA markers

This study presents linkage maps for the Pacific abalone (Haliotis discus hannai) based on 180 microsatellite DNA markers. Linkage mapping was performed using three F1 outbred families, and a composite linkage map for each sex was generated by incorporating map information from the multiple families. A total of 160 markers are placed on the consolidated female map and 167 markers on the male map. The numbers of linkage groups in the composite female and male maps are 19 and 18, respectively; however, by aligning the two maps, 18 linkage groups are formed, which are consistent with the haploid chromosome number of H. discus hannai. The female map spans 888.1 cM (Kosambi) with an average spacing of 6.3 cM; the male map spans 702.4 cM with an average spacing of 4.7 cM. However, we encountered several linkage groups that show a high level of heterogeneity in recombination rate between families even within the same sex, which reduces the precision of the consolidated maps. Nevertheless, we suggest that the composite maps are of significant potential use as a scaffold to further extend the coverage of the H. discus hannai genome with additional markers.