A high-density SNP genetic linkage map for the silver-lipped pearl oyster, Pinctada maxima: a valuable resource for gene localisation and marker-assisted selection
© Jones et al.; licensee BioMed Central Ltd. 2013
Received: 18 July 2013
Accepted: 18 November 2013
Published: 20 November 2013
The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species extensively farmed for the highly sought "South Sea" pearls. Traditional breeding programs have been initiated for this species in order to select for improved pearl quality, but many economic traits under selection are complex, polygenic and confounded with environmental factors, limiting the accuracy of selection. The incorporation of a marker-assisted selection (MAS) breeding approach would greatly benefit pearl breeding programs by allowing the direct selection of genes responsible for pearl quality. However, before MAS can be incorporated, substantial genomic resources such as genetic linkage maps need to be generated. The construction of a high-density genetic linkage map for P. maxima is not only essential for unravelling the genomic architecture of complex pearl quality traits, but also provides indispensable information on the genome structure of pearl oysters.
A total of 1,189 informative genome-wide single nucleotide polymorphisms (SNPs) were incorporated into linkage map construction. The final linkage map consisted of 887 SNPs in 14 linkage groups, spans a total genetic distance of 831.7 centimorgans (cM), and covers an estimated 96% of the P. maxima genome. Assessment of sex-specific recombination across all linkage groups revealed limited overall heterochiasmy between the sexes (i.e. 1.15:1 F/M map length ratio). However, there were pronounced localised differences throughout the linkage groups, whereby male recombination was suppressed near the centromeres compared to female recombination, but inflated towards telomeric regions. Mean values of LD for adjacent SNP pairs suggest that a higher density of markers will be required for powerful genome-wide association studies. Finally, numerous nacre biomineralization genes were localised providing novel positional information for these genes.
This high-density SNP genetic map is the first comprehensive linkage map for any pearl oyster species. It provides an essential genomic tool facilitating studies investigating the genomic architecture of complex trait variation and identifying quantitative trait loci for economically important traits useful in genetic selection programs within the P. maxima pearling industry. Furthermore, this map provides a foundation for further research aiming to improve our understanding of the dynamic process of biomineralization, and pearl oyster evolution and synteny.
The silver-lipped pearl oyster, Pinctada maxima, is an important tropical aquaculture species that, along with P. margaritifera, produces almost 50% of marketed pearls worldwide by value . However, like most aquaculture industries, pearl culture is still in its infancy compared to terrestrial animal production systems and has yet to establish advanced selective breeding programs required for industry advancement. Profitability of the P. maxima industry is driven primarily by the grading of the five pearl quality traits: shape, size, colour, lustre and surface complexion. Large variation is observed during harvest for each of these traits, presenting the potential to increase industry profitability through selective breeding. Although traditional animal improvement methods have had some success improving traits which are easy to measure in candidates under selection (i.e. animal growth) , they are not particularly effective for complex pearl quality traits, which are generally hard to measure, expressed late in life and generally have low heritability [3, 4]. Promising developments in livestock genomics are opening up opportunities, allowing genomic information to be incorporated into breeding programs in order to increase the rate of genetic gain for complex commercial traits in oyster. The current impediment to the implementation of genomic approaches in mollusc breeding programs, however, is a significant lack of genomic resources such as genome-wide molecular markers, genomic maps and genome sequences [2, 5, 6].
A robust high-density genetic linkage map for P. maxima is a fundamental precursor to understanding the architecture and evolution of pearl oyster genomes, determining the genetic basis of complex phenotypic traits under natural and industrial settings, and identifying genes and quantitative trait loci (QTL) associated with bivalve shell biomineralization. Such resources are invaluable for the development and incorporation of marker-assisted selection (MAS) into breeding programs aiming to fast track improvements in pearl quality. Presently, no genetic maps are available for P. maxima, with information on this species’ genome largely limited to the general physical description of its chromosomes (i.e. N = 14, 10 submeta- or metacentric, and four telocentric chromosomes) .
Preliminary genetic linkage maps have been developed for only a few bivalves, including the edible oysters Crassostrea virginica , C. gigas [9, 10], Ostrea edulis  and one pearl oyster species, Pinctada fucata martensii . However, information from these maps is of limited use in P. maxima for molecular breeding studies, as they either consist of non-transferable markers [i.e. amplified fragment length polymorphisms (AFLPs)], have low marker density (100–200 markers), or the original species is phylogenetically too distant to be useful in a comparative genetic mapping approach .
Alongside the lack of genomic resources, several fundamental aspects of pearl oyster biology still remain unclear. For example, one of the most striking features of pearl oysters is that they are non-obligatory protandrous hermaphrodites (i.e. mature first as males and later change to females). In P. margaritifera for instance, individuals develop as males and remain so for the first two years of life before progressively changing to females reaching a sex ratio close to 1:1 at around 8 years old . Sex change is known to be largely driven by environmental factors such as stress. However, the genetic determinates of this unusual life history have yet to be investigated in detail. Genetic linkage maps may be implemented to unravel some of the genetic determinates of sex differentiation and sex change in oysters. Linkage maps are also highly desirable for evolutionary genetic research and comparative mapping which would improve our understanding of pearl oyster chromosome evolution and help identify homologous chromosomal segments involved in the genetic control of economical and adaptive traits for species in the genus Pinctada.
This study aimed to construct medium to high density sex-average and sex-specific genetic linkage maps for the silver-lipped pearl oyster, P. maxima, by utilising a recently developed single nucleotide polymorphism (SNP) array . Following robust linkage map construction, this study evaluates heterochiasmy between the sexes, extent of linkage disequilibrium (LD) across the genome, and the localization of important biomineralization genes. This comprehensive genetic resource allows for the first time the ability to obtain new insights into the biological and genomic architecture of this important marine species, including the identification of the genetic basis of complex phenotypic traits.
Reference mapping families and DNA extraction
SNP selection, genotyping and data integrity
A total of 1,189 informative P. maxima type I SNP markers (developed and validated in ) were selected for linkage map construction. Strict data integrity measures were implemented to ensure that only the most accurate SNPs were included in the analysis, as even a small proportion of genetic marker errors can dramatically affect the accuracy of genetic linkage maps . Briefly, SNPs were selected if they amplified successfully, returned clear genotype calling clusters, had a minor allele frequency (MAF) > 0.01, did not deviate from Hardy-Weinberg equilibrium (P value cut off), conformed to Mendelian inheritance (MI) patterns from parent to offspring, did not exhibit duplicated SNP clustering patterns and had a call rate > 90%. Of the available 1,189 SNPs, 1,167 were informative for the subset of 351 oysters belonging to the mapping families and included in linkage map analysis (Additional file 1). All SNPs have previously been annotated with gene identity and Gene Ontology (GO) terms as described in Jones et al. .
Map construction and genome coverage
To generate the most accurate genetic linkage maps two versions of CRI-MAP  were employed. Firstly, a modified version of CRI-MAP developed by Liu and Grosz  was utilised to calculate sex-average and sex-specific two-point recombination rates and logarithm of the odds (LOD) scores for all pairs of markers. Linkage groups (LGs) were then identified using AUTOGROUP  by identifying markers that co-segregate with pairwise LOD scores > 5. AUTOGROUP utilises an iterative process with a succession of parameters decreasing in stringency through five layers to ensure accurate groups are produced. Each layer consists of the following four parameters; the minimum threshold of LOD scores for linkage to be included, the minimum number of informative meiosis for a marker to be included, the maximum number of linkage groups a marker is allowed for having linkages, and the minimum threshold for the linkage ratio to which a marker shows qualified links to the best linkage group. By layer, the parameters were a) layer one: 100, 2.0, 2, 0.9, b) layer two: 50, 1.5, 5, 0.8, c) layer three: 20, 1.0, 8, 0.7, d) layer four 10, 0.5, 10, 0.6, and e) layer five: 5, 0.1, 15, 0.5. Layer five defines the cut off for a marker to be included in a linkage group. Linkage groups were numbered in order of decreasing number of markers placed within each linkage group during the AUTOGROUP phase. Once linkage groups of markers were established, predefined haplogroups of tightly linked loci (i.e. LOD ≥ 3.0 and theta ≤ 0.03) were identified using HAPLOGROUP. This produced a subset of primary (most informative) markers from each haplogroup that were selected for initial construction of a framework map.
The BUILD and FLIPS commands of the second version of CRI-MAP (2.503) modified by Jill Maddox and Ian Evans (unpublished data) were utilised to determine the marker order within each linkage group as it has been designed to deal with large datasets more efficiently. The linkage mapping strategy consisted of a hierarchal approach whereby markers were included if they could be assigned a position over the next most likely position firstly with a LOD score threshold cut off of LOD3 (the standard threshold for framework markers) which represents a 1:1000 chance of a marker being placed incorrectly. After placement of the LOD3 markers, the omitted non-framework markers and remaining secondary haplogroup markers were incorporated into the framework map using successive BUILD commands at the following decreasing LOD threshold cut offs; LOD2 (1:100 chance of incorrect marker placement), LOD1 (1:10 chance of incorrect marker placement) and finally the most likely position of remaining loci. For each BUILD at each LOD threshold cut-off, the marker order was verified using the FLIPS function with a moving window of five markers (FLIPS5). When a better marker order was established after FLIPS5, marker order was resolved and FLIPS5 was re-run until no further changes were apparent. CHROMPIC was then employed to ensure no incorrect double recombinants were present which may indicate incorrect marker positioning. Erroneous genotype calls were corrected and any markers identified with double recombinants were reanalysed with BUILD and FLIPS to determine if the double recombinants were real or the marker position was incorrect. Any markers with unresolved double recombinants were excluded and FLIPS5 was re-run to ensure the remaining marker order remained correct. The final map is referred to as the comprehensive map . Sex-specific maps were also constructed using the sex-average marker order and recalculating marker intervals based on separate male and female informative meiosis events. Final map distances were calculated using the option FIXED. The Kosambi mapping function  was used for all cM calculations and all maps were drawn using MapDraw version 2.2 .
To validate the map ordering of CRI-MAP, markers belonging to a large linkage group with a range of informative loci (LG8 - established by AUTOGROUP in CRI-MAP) were chosen to build an independent sex-average comprehensive linkage map with CarthaGène version 1.0 which incorporates an EM (expectation-maximization) algorithm and a local search technique to build a maximum likelihood map . The phasing function in TMAP version 1.1  was utilised to generate input files for CarthaGène which incorporated the eight reference mapping families (Figure 1). The map was built using the same hierarchical mapping LOD thresholds as outlined above (LOD3, LOD2, LOD1 and most likely position) using recurrent executions of "build", "polish" and "flips".
To calculate genome coverage of the linkage maps the observed and expected genome lengths need to be established. The observed genome length (Goa) was simply the addition of all observed linkage group lengths and the expected genome length (Ge) was calculated by multiplying the length (cM) of each linkage group by (m + 1)/(m - 1), where m is the number of loci in each linkage group see . The total expected genome length was the sum of Ge from all linkage groups. Genome coverage (Coa), was calculated by dividing Goa by Ge see .
Segregation distortion, defined as the deviation from Mendelian inheritance of co-dominant alleles, may be present as a result of gametic selection or post-zygotic selection. To determine if such biological processes are present, segregation distortion was investigated using log-likelihood ratio tests for goodness of fit to Mendelian expectations in the software suit LINKMFEX version 2.4 . Here, G-values were calculated for all markers across all mothers and fathers of each family and subsequently tested using the heterogeneity G-test as described in Sokal and Rohlf . For each marker, G Total (sum of G values across all parents) and G-Pooled [calculated from the sum of allele specific (A and B) and total numbers (N) of co-informative events] were calculated and compared to determine the direction of the distortion if present. Heterogeneity was then calculated by subtracting the Total G value from the Pooled G value .
Sex-specific and family-specific recombination heterogeneity
where, is the joint sex-specific recombination rate and represents the recombination rate when equal male and female recombination fractions are assumed. For each test, a false discovery rate (FDR) correction was applied to correct for multiple comparisons and minimise false positives .
To detect any differences in sex-specific recombination rates, ratios of female-to-male map distances were calculated (R = X f /X m ) for each interval and linkage group as well as over the entire map. In addition, standardised marker interval distances were calculated for each sex [standardised interval distance = 100 * (interval distance/total LG length)] and plotted against one another. For all linkage groups, distinct slopes were observed along the length of the linkage group. Breakpoints between the distinct slopes for each linkage group were assigned by visual inspection. Each slope was analysed using a simple linear regression of two continuous variables (female and male) as they represent biologically real differences. For regression analysis, data that produced each slope were grouped into three groups (left, middle and right) for all linkage groups except LG7 & LG9 where only two groups were produced (Additional file 2).
Here, represents the LOD scores maximum likelihood estimation (MLE) for the ith F1 reference family for a pair of markers, with being the total LOD score MLE of all ith reference families.
Extent of linkage disequilibrium
The extent of LD is an important consideration for association mapping as it indicates the relative size of chromosomal segments shared amongst individuals within a population, and thus determines the number of theoretical markers necessary to detect genetic associations to quantitative traits . Two commonly used estimates of LD, r 2  and D’ , were computed using GOLD software . The LD estimates were computed among all 1,167 SNPs using genotypic data on 995 oysters (the additional 660 oysters either have no pedigree information or belong to smaller families not suitable for linkage mapping). The extent of LD among SNPs, within and across the linkage groups, was estimated using position of SNPs on the current linkage map.
Genotyping, pedigrees and data integrity
The validation success of SNPs included on the custom genotyping array is detailed in Jones et al. . Strict data integrity on the SNPs based on a genotyped population consisting of 525 individuals produced a total of 1,189 SNPs suitable for linkage mapping . Of these, 1,167 produced polymorphic genotypes (MAF > 0.01) across the subset of 351 oysters belonging to the reference mapping families with an average genotyping call rate of > 99.4% (Additional file 1).
Summary statistics of the sex-average, female and male linkage maps of P. maxima
Map length (cM)
Expected genome length (Ge)
Average interval (cM)
No of intervals (Sex-Av.)
No. of SNPs
Female: Male ratio
0.99 (+/- 1.09)
1.34 (+/- 1.90)
1.30 (+/- 1.25)
1.20 (+/- 1.29)
1.63 (+/- 2.82)
1.18 (+/- 1.21)
1.14 (+/- 1.21)
2.00 (+/- 2.72)
1.67 (+/- 1.49)
1.41 (+/- 1.64)
2.56 (+/- 3.04)
2.54 (+/- 3.16)
1.28 (+/- 1.27)
1.96 (+/- 3.54)
1.81 (+/- 1.72)
1.68 (+/- 1.64)
3.20 (+/- 3.99)
1.58 (+/- 1.51)
1.94 (+/- 2.62)
6.26 (+/- 9.67)
1.88 (+/- 2.16)
2.14 (+/- 2.37)
4.03 (+/- 6.76)
1.64 (+/- 1.70)
2.36 (+/- 2.80)
5.40 (+/- 7.15)
2.81 (+/- 1.70)
2.61 (+/- 3.23)
4.41 (+/- 5.68)
2.57 (+/- 3.00)
2.30 (+/- 2.42)
3.94 (+/- 4.32)
2.08 (+/- 1.71)
2.42 (+/- 2.75)
3.46 (+/- 3.39)
3.28 (+/- 3.75)
3.46 (+/- 3.18)
6.43 (+/- 7.07)
2.60 (+/- 2.36)
3.07 (+/- 3.48)
5.04 (+/- 8.09)
3.50 (+/- 1.96)
Sex-specific and family-specific recombination heterogeneity
Investigations into family specific heterogeneity confirm that observed sex-specific recombination is truly caused by the sexes and not individual parental F1 individuals biasing the data. Only one interval in LG1 on the sex-average map (c10004 - c13798) returned significant recombination heterogeneity after FDR (χ 2 = 21.6, P = 0.0002, df = 4). This deviation was explained by deviations in only two of the eight families (130×148 and 131×145) providing evidence that the mapping parents are relatively homogeneous within the sexes for recombination differences.
Significant segregation distortions were detected in 121 of the 887 mapped SNPs and seven of the sixteen mapping parents following FDR correction (mean corrected alpha of 0.003) (Additional file 20). The majority (64.9%) of these distortions were localized to the mapping family 103×102 and to linkage groups 2, 4 and 10. As no significant family specific heterogeneity was detected for these distortions, they are not thought to be influencing calculations of mapping distances. However, to be conservative, only markers that did not cause conflicts in map position were mapped.
Biomineralization gene mapping
Positional information of biomineralization candidate genes can assist in determining which genes influence pearl quality traits by comparing their positions to QTLs. As described in Jones et al. , numerous SNPs were designed within known biomineralization genes. A total of nine SNPs designed within six biomineralization gene homologs were successfully mapped. These genes were Calreticulin, chitin synthase 1 (CS1), Lustrin A, N19, Pinctada fucata mantle gene (PFMG) complex and Pif177. Two SNPs from Lustin A were mapped, clustering together in a telomeric region of LG9 (Figure 3) and three SNPs designed within the PFMG complex were mapped to the centre of LG4 (Figure 2) along with the SNPs designed in Pif177 and CS1. A SNP from Calreticulin was mapped to the centre of LG1 and N19 was mapped to the end of LG1 (Figure 2).
Extent of linkage disequilibrium
Overall distribution of linkage disequilibrium (LD) estimates ( r 2 and D’ ) for all, non-syntenic and syntenic SNPs
Range of estimate
Number of SNP pairs
0 - 0.1
0.1 - 0.2
0.2 - 0.3
0.3 - 0.4
0.4 - 0.5
0.5 - 0.6
0.6 - 0.7
0.7 - 0.8
0.8 - 0.9
0.9 - 1
Mean (± SD) and median of r 2 and D’ linkage disequilibrium estimates over distance for all linkage groups
r 2 mean (± SD)
r 2 median
D’mean (± SD)
0.082 (± 0.138)
0.519 (± 0.321)
0 - 1 cM
0.075 (± 0.125)
0.494 (± 0.323)
1 - 2 cM
0.061 (± 0.094)
0.490 (± 0.311)
2 - 5 cM
0.058 (± 0.088)
0.472 (± 0.315)
5 - 10 cM
0.051 (± 0.078)
0.440 (± 0.299)
10 - 20 cM
0.042 (± 0.064)
0.410 (± 0.293)
20 - 50 cM
0.022 (± 0.034)
0.319 (± 0.265)
> 50 cM
0.020 (± 0.029)
0.307 (± 0.266)
The robust high-density genetic linkage map for P. maxima presented here is the most comprehensive map to date for any bivalve species. With the combination of physical mapping data, this map will contribute to a better understanding of the genome structure, function and evolution for P. maxima, and other species within the genus Pinctada. In addition, the identification of genetic associations and QTLs for commercial traits will be highly valuable to the pearling industry as this information will be directly used for genetic improvement of complex traits in farmed stock.
Genotyping, marker selection and pedigree validation
Missing genotypes or typing errors are known to interfere with the ordering of SNPs leading to incorrect estimation of map lengths . Therefore, thorough genotypic data integrity is vital for the generation of accurate maps. The majority of data integrity applied to this dataset has been described in Jones et al. . This previous work has provided a highly refined list of SNPs useful for linkage mapping analysis by excluding all SNPs with low polymorphism and removing erroneous genotype errors including deviations from HWE, Mendelian incompatibility, SNP duplication, low MAF and low call rates. The average number of informative meiosis events (83,377 phase known and 259,844 phase unknown for sex-average) for this subset of 1,167 SNPs over our eight families was high ensuring fine resolution throughout the map. However, in some cases, the order of closely linked markers (N = 397, recombination = 0) could not be determined even though the overall number of informative meiosis was relatively high. For these loci, a small proportion cannot be separated due to limitations in SNP discovery see , while others will require more informative meiosis events to eventually separate.
Parental relationships of mapping families utilised in this study were rigorously tested and confirmed using parentage analysis during previous projects, minimising the possibility of pedigree errors. This was further confirmed through testing for MI errors during map construction. No families were found with Mendelian Inheritance errors across many informative loci. Overall, the level of Mendelian inconsistencies was extremely low for a custom array consisting of novel SNPs. A few sporadic MI errors could be attributed to poor DNA quality in a few samples which were subsequently removed.
One of the problems in linkage mapping of oysters is that moderate distortions from expected Mendelian segregation are common . Markers that exhibit segregation distortion can potentially influence marker positions and linkage relationships, however, the presence of moderate segregation distortion has been reported to have little effect on the overall construction of linkage maps [20, 39], and maps have been successfully constructed in species exhibiting moderate segregation distortions . Additionally, the inclusion of distorted markers in mapping can be beneficial as they may have an association to genes that affect fitness and survival, particularly larval mortality, and they may also help with understanding the distribution of deleterious recessive genes throughout the genome.
The extent of segregation distortion throughout the P. maxima linkage map was investigated to determine their influence on marker order and mapping distances and to determine if markers exhibiting distortion clustered together. A total of 121 mapped SNP loci showed at least one significant distortion in a family after stringent FDR correction (average FDR alpha value of 0.0032). A high proportion of these distortions (79.8%) were localised to specific linkage groups (25.4% to LG1, 20.9% to LG2, 15.7% to LG4, 9.7% to LG3 and 8.2% to LG10) indicating a true biological phenomenon is in effect instead of random artefacts . Localised segregation distortion has been reported as a common feature in numerous bivalve species including, the Pacific oyster (Crassostrea gigas) [9, 39, 41], the European flat oyster (Ostrea edulis)  and Pinctada martensii . Since at least C. gigas is known to have a high genetic load , such distortions in these bivalves have been explained by zygotic viability selection due to the presence of deleterious recessive genes [9, 11, 12]. Segregation distortions reported here may also be attributed to the presence of deleterious recessive genes as has been observed in C. gigas, O. edulis and P. martensii, however, further research is warranted to confirm this.
Map construction and estimated genome size
Linkage map construction resulted in the generation of 14 linkage groups that correspond to the 14 haploid chromosomes of P. maxima . Approximately 76% of the SNPs (887 out of 1,167) were placed on the linkage map (Additional file 1). This comprehensive first-generation linkage map is a substantial resource and is a large improvement on any bivalve map to date with reference to the number of markers mapped (previous average number of markers mapped of 191) and genome coverage (previous average genome coverage of 80.4%) [8–12, 43]. The number of markers on this map (N = 887) more than doubles any previous attempt in bivalves and the predicted genome coverage (96%) is much higher than an average of 80% reported in previous bivalve maps. In addition, the distribution of inter-SNP spacing throughout the map demonstrates an even spread of markers throughout the genome with over 49% of the inter-marker distances being less than 1 cM (Median inter-marker distance of 1.03 cM).
The estimated genome size for P. maxima based on the sex-average linkage map is 865.6 cM (954.6 cM for the female map and 826.1 cM for the male map). This is significantly less than expected genome length reported for a previous linkage map for Pinctada martensii (1862.9 cM for the female map and 1838.4 for the male map) . Assuming similar recombination rates between the species, one possible explanation for this is that the inclusion of more markers refines positions and reduces the total cM of each linkage group as acknowledged by Shi et al. . Maps of low density are commonly longer than maps of high density and as more markers are added, map length decreases [8, 9, 44]. A low marker density in the previous P. martensii map is most likely the cause of the overestimation of genome size. The length of the P. maxima linkage groups in present study varied from 48.3 cM to 75.6 cM and exhibited a negative relationship with the number of markers mapped per linkage group (Table 1). As a result, linkage groups of smaller sizes showed similar recombination rates than those of larger sizes. This too may be a result of inflated map distances for linkage groups with fewer markers mapped.
Sex-specific maps and recombination rates
Sex-specific differences in recombination rates are not uncommon and have been reported in numerous vertebrate [45–48] and invertebrate [10, 12, 49, 50] species. In accordance to the Haldane rule, for organisms with a chromosomal mechanism of sex determination, recombination should be more frequent in the homogametic sex than in the heterogametic sex [32, 33, 51]. This observation has been termed heterochiasmy. However, many exceptions to this rule have been demonstrated including the tammar wallaby , the great reed warbler  and the saltwater crocodile . In addition, reports of sex-specific recombination in species without heteromorphic sex chromosomes have become apparent [47, 48, 52, 53]. Oysters are one taxa that lack specialised heteromorphic sex chromosomes [7, 54], but exhibit sex-specific recombination [this study, 8, 9]. Results observed here for P. maxima show that the male map (826.1 cM) is shorter than the female map (954.64 cM), suggesting a slight female bias in recombination with an overall ratio of female-to-male recombination of 1.15:1 and ratios reaching 2.12:1 in LG13. This is comparable to previous ratios of sex-specific recombination in oysters that range from 1.07:1–1.51:1 [8, 9, 12] and other aquaculture species (female-to-male ratios range from 1.2:1 – 3.25:1) [48, 52, 55, 56]. Such proliferation of studies that report female biased sex-specific recombination in species with no specialised sex chromosomes suggest that there must be another underling phenomenon of the timing, duration or biological features associated with meiosis that is responsible for the observed differences between the sexes.
Sex-specific recombination rates are also known to differ throughout regions within the genome . Dramatic localised sex-specific differences were detected throughout the P. maxima maps where male recombination rates were supressed relative to the female rates in areas proximal to centromeres, but elevated in regions distal to centromeres, with females showing the opposite pattern (see Figure 5 and Additional files 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15 and 16) [58, 59]. The expansion of the male genetic map in telemetric regions indicates that chiasmata would be found more frequently near the telomeres in meiosis in males compared to females. Similarly, chiasmata would be more common in centromeric regions during oogenesis. Such pronounced localised differences in recombination rates have not been previously reported in oysters , but are quite widespread throughout other aquaculture species including the rainbow trout (Oncorhynchus mykiss) , Atlantic salmon (Salmo salar)  and the zebrafish (Danio rerio) . This unusual pattern of sex-specific recombination is not well understood, however, several theories have been suggested [reviewed by Miles et al. . Briefly, sex-specific recombination could have been caused by a) differing environments in which the germ cells develop , b) temporal differences in initiation of meiosis between the sexes  and c) differences in the pairing and synapses of homologs at meiosis that cause different exchange patterns in oocytes and spermatocytes . However, further research is required to confirm these theories in P. maxima.
Nevertheless, it is remarkable that such strong sex-specific recombination patterns are present in a species without differentiated sex chromosomes, exhibits no sexual dimorphism and is a protandrous hermaphrodite (maturing first as a male and able to switch sex after 2 years of age). Definitely, the strong evidence of sex-specific recombination presented here may aid in identifying the mechanism behind sex-specific recombination, especially for species without differentiated sex chromosomes. To truly, elucidate the basis of sex-specific recombination, cytogenetic analysis of female and male meiosis would be required. The unusual life history of P. maxima (a protandrous hermaphrodite) may allow the estimation of female and male recombination rates in the same individual removing any effect of genetic background on such estimations. As suggested by Franch et al. , hermaphroditic species might play a crucial role in dissecting the contribution of sex-determining and sex-differentiating genes on meiotic recombination .
Biomineralization gene mapping
Nine mapped SNPs were designed within six homolog sequences of known candidate genes for biomineralization (Calreticulin, CS1, Lustrin A, N19, PFMG complex and Pif177). Clustering of these SNPs in four locations within the genome (telemetric region of LG9, centre of LG4, centre of LG1 and telemetric region of LG1) provides evidence that these regions may have a strong influence on biomineralization for P. maxima and signposts these regions for further investigation to determine true associations to biomineralization processes. Five of the nine biomineralization gene SNPs clustered within 27.2 cM near the centre of LG4 (total length of 71.9 cM). These SNPs represented three genes, the PFMG1, Pif177 and CS1 which are highly expressed in mantle tissue of pearl oysters and are known to be involved in nacre formation [63–65]. More specifically, PFMG1 and Pif177 (both initially described in P. fucata) are key calcium-binding proteins that specifically bind aragonite crystals and regulate nucleation and precipitation during nacre formation [64, 65]. PFMG1 and Pif177 have also previously been co-localised during EST clustering analysis for P. fucata sequences . The second gene, Pif177 consists of two proteins, Pif80 and Pif97 which are encoded by a single mRNA . Pif80, Pif97 and N16 (another nacre biomineralization protein) work in collaboration (along with chitin) to initiate aragonite crystallization and orientate the stacking of aragonite tablets in nacreous layers [65, 67]. Interestingly, the third protein clustering at this region, CS1, is also involved in the gene complex described above. CS1 is a key enzyme responsible for the deposition of chitin, a polysaccharide integral for calcium carbonate biomineral formation in mollusc shells . The co-localisation of these three major nacre biomineralization genes to central parts of LG4 provides strong evidence that this region is a hot spot for nacre biomineralization genes and would become a prime target for studies aiming to identify QTL for commercially valuable pearl quality traits.
Four additional SNPs designed within three genes (Lustrin A, Calreticulin and N19) were localised to three other regions throughout the linkage map. Two SNPs (Lustrin_A_c15856_1 and Lustrin_A_c15856_2) designed within a contig homologous to Lustrin A (c15856) were mapped to the same position (0 cM intermarker distance) close to a telomere of LG9. The SNP Calreticulin_c2420_1 (designed within a contig homologous to Calreticulin) was mapped to the centre of LG1 and N19_c591_1 (designed within a contig homologous to N19) was mapped to a telemetric region of LG1. Specific functions have been ascribed to each of these three genes, including; conferring elastic resilience to the molluscan shell and maintaining the structure and protein compounds of nacre for Lustrin A ; calcium binding, transport and storage during biomineralization for Calreticulin ; and finally, having a negative regulatory role in calcification for N19 . The localisation of these six biomineralization genes will increase the accuracy of identifying regions of interest for researchers interested in identifying genetic association to important nacre biomineralization genes and will also be important for comparative mapping studies investigating genome evolution and synteny.
Extent of linkage disequilibrium
Association studies aiming to identify genetic variations or quantitative trait nucleotide (QTN) that explain a large proportion of the phenotypic variance in a quantitative trait rely on the co-segregation of QTNs with the surrounding genetic markers or loci. If the marker and QTN are sufficiently close, this association will remain intact within the population over many generations . Such non-random association between loci is termed linkage disequilibrium (LD). The extent of LD is therefore important as it defines the density of genome-wide makers necessary for association analysis to detect markers associated with traits of commercial interest and are also in LD with QTNs. Generally, higher marker density is beneficial, although, if the extent of LD throughput the genome is high, fewer markers may be sufficient for association studies . Two estimates of the extent of LD were utilised in this study, D’ and r 2 . The D’ estimate of LD is suggested to be a good measure for the extent of LD in a population and variation in LD throughout the genome as it focuses on historical recombination. However, D’ is known to be more influenced by allelic frequency variation than the r 2 estimate . As such, r 2 is more useful in predicting the power of association mapping.
The LD estimates presented in this study are based on 995 oysters which include an additional 660 oysters to those utilised for linkage mapping. Additional animals from smaller families and unknown pedigree are particularly suitable for computing LD estimates. Estimates of LD among non-syntenic (on different chromosomes) SNP pairs represent background variation observed within the data. The mean estimate of LD among non-syntenic SNPs for P. maxima (mean r 2 of 0.020) is generally higher when compared to well characterised species (i.e. bovine with mean r 2 of 0.003; ). This may be due to the high relationship among animals in this population as compared to bovine. Therefore, the comparatively higher background LD estimates of non-syntenic SNPs in this study are not unexpected. For syntenic (on the same chromosome) SNP pairs greater than 50 cM apart, estimates of LD were similar to that of non-syntenic SNPs (r 2 and D’ of 0.014 and 0.307 respectively). This indicates that recombination between these long range SNP pairs is relatively high, and SNPs on distal ends of the chromosomes are behaving in a similar manner as non-syntenic SNPs.
Estimates of LD usually decline as map distance increases in most species. Here, the decline in LD over map distance is gradual for both D’ and r 2 estimates in P. maxima (Table 3 and Figure 8). However, the mean LD estimates among closely spaced markers are lower as compared to other well characterised species (e.g. bovine, human ) which suggests a low extent of LD within the current population. Limited studies of LD have been reported in invertebrates. The LD estimates reported here are contrary to what has been observed in another aquaculture species, the Pacific white shrimp (Litopenaeus vannamei). For L. vannamei, a steeper decline in LD with map distance suggests smaller LD blocks . In addition, estimates of r 2 for syntenic SNP pairs greater than 50 cM apart were higher than that observed in this study (0.15 compared to 0.014). This is likely due to a difference in the effective population size (N e ) between the two studies. The current study was based on multiple families derived from outbred populations (higher N e ), whereas, estimates of LD for L. vannamei were based on only 144 individuals from six family lines (lower N e ). However, the LD estimates of L. vannamei are probably more typical of aquaculture species in general as these are usually derived from limited numbers of stocks.
The low LD estimates for short range (0–1 cM) syntenic SNPs (r 2 and D’ of 0.083 and 0.519 respectively), and gradual decline in LD, suggests limited short range LD at the current marker density. To fully evaluate short range LD in this population, marker density needs to be increased. As a result, these LD estimates must be treated with caution. With a higher marker density, the decline of LD throughout the genome may be more pronounced once better estimates can be made between SNPs pairs at smaller map intervals. Furthermore, in this study, the extent of LD was compared against the linkage map, however, both linkage and LD maps are calculated using recombination rates. The extent of LD across a genome is better understood when presented against the physical map positions. Nevertheless, in the absence of a physical map, these results provide a preliminary estimate of broad patterns of LD observed within the oyster genome for this population. Even though the present SNP density will be useful for first-pass QTL and genome-wide association studies (GWAS), the low values of r 2 for most adjacent SNPs pairs suggest that density should be increased before fine-scale trait and LD mapping across the P. maxima genome is attempted.
This research developed a high-density genetic linkage map suitable for studies aiming to identify gene associations and QTLs for commercially important traits such as shell growth, pearl size, nacre colour and surface complexion in the silver-lipped pearl oyster. The genetic linkage map will be particularly useful for the mapping of QTLs in this species, especially since it is of high density, the mapped SNPs are genic, and numerous regions have been flagged with genes known to be involved in nacre biomineralization. The density of this linkage map would also be sufficient for preliminary GWAS analysis, however, higher density would be more appropriate considering the low extent of LD throughout the genome.
Finally, the transferability of mapped SNPs to species within the genus Pinctada has previously been shown to be high . For example, conversion rates of SNPs between species closely related to P. maxima (i.e. 61.3% in P. margaritifera and 58.5% in P. mazatlantica), illustrate the high utility for the P. maxima map in comparative mapping studies. When other genomic resources become available for bivalve species, comparative mapping studies utilising our linkage map will provide insights into many fundamental questions in the localization of genes, conservation of gene content and order, genome evolution and synteny in bivalves.
Amplified fragment length polymorphism
Chitin synthase 1
False discovery rate
Genome-wide association study
Logarithm of odds
Minor allele frequency
Pinctada fucata mantle gene
Quantitative trait loci
Quantitative trait nucleotide
Single nucleotide polymorphism.
This work was funded by Australia Research Council funding to D.R.J. Array genotyping was conducted at PathWest medical laboratories in Perth, Western. Atlas South Sea Pearls provided the oysters and the facilities necessary to breed, raise and maintain reference mapping families.
- Tisdell C, Poirine B: Economics of pearl farming. The pearl oyster. Edited by: Southgate PC, Lucas JS. 2008, Amsterdam, The Netherlands: Elsevier B V, 473-495.View ArticleGoogle Scholar
- Wada KT, Jerry DR: Population genetics and stock improvement. The pearl oyster. Edited by: Southgate PC, Lucas JS. 2008, Amsterdam, The Netherlands: Elsevier B V, 437-471.View ArticleGoogle Scholar
- Kvingedal R, Evans BS, Lind CE, Taylor JJU, Dupont-Nivet M, Jerry DR: Population and family growth response to different rearing location, heritability estimates and genotype × environment interaction in the silver-lip pearl oyster (Pinctada maxima). Aquaculture. 2010, 304 (1–4): 1-6.View ArticleGoogle Scholar
- Jerry DR, Kvingedal R, Lind CE, Evans BS, Taylor JJU, Safari AE: Donor-oyster derived heritability estimates and the effect of genotype x environment interaction on the production of pearl quality traits in the silver-lip pearl oyster, Pinctada maxima. Aquaculture. 2012, 338–341: 66-71.View ArticleGoogle Scholar
- Saavedra C, Bachere E: Bivalve genomics. Aquaculture. 2006, 256 (1–4): 1-14.View ArticleGoogle Scholar
- Curole JP, Hedgecock D: Bivalve genomics: complications, challenges, and future perspectives. Aquaculture genome technologies. 2007, Oxford, UK: Blackwell Publishing Ltd, 525-544.View ArticleGoogle Scholar
- Wada K, Komaru A: Karyotypes in five species of Pteriidae (Bivalvia: Pteriomorphia). Venus. 1985, 44 (3): 183-192.Google Scholar
- Yu Z, Guo X: Genetic linkage map of the eastern oyster Crassostrea virginica gmelin. Biol Bull. 2003, 204 (3): 327-338. 10.2307/1543603.View ArticlePubMedGoogle Scholar
- Li L, Guo X: AFLP-based genetic linkage maps of the pacific oyster Crassostrea gigas thunberg. Mar Biotechnol. 2004, 6 (1): 26-36. 10.1007/s10126-003-0001-0.View ArticlePubMedGoogle Scholar
- Hubert S, Hedgecock D: Linkage maps of microsatellite DNA markers for the pacific oyster Crassostrea gigas. Genetics. 2004, 168 (1): 351-362. 10.1534/genetics.104.027342.PubMed CentralView ArticlePubMedGoogle Scholar
- Lallias D, Beaumont AR, Haley CS, Boudry P, Heurtebise S, Lapègue S: A first-generation genetic linkage map of the european flat oyster Ostrea edulis (L.) based on AFLP and microsatellite markers. Anim Genet. 2007, 38 (6): 560-568. 10.1111/j.1365-2052.2007.01647.x.View ArticlePubMedGoogle Scholar
- Shi Y, Kui H, Guo X, Gu Z, Wang Y, Wang A: Genetic linkage map of the pearl oyster, Pinctada martensii (dunker). Aquacult Res. 2009, 41 (1): 35-44.View ArticleGoogle Scholar
- Cunha R, Blanc F, Bonhomme F, Arnaud-Haond S: Evolutionary patterns in pearl oysters of the genus Pinctada (Bivalvia: Pteriidae). Mar Biotechnol. 2011, 13 (2): 181-192. 10.1007/s10126-010-9278-y.View ArticlePubMedGoogle Scholar
- Chávez-Villalba J, Soyez C, Huvet A, Gueguen Y, Lo C, Moullac GL: Determination of gender in the pearl oyster Pinctada margaritifera. J Shellfish Res. 2011, 30 (2): 231-240. 10.2983/035.030.0206.View ArticleGoogle Scholar
- Jones DB, Jerry DR, Forêt S, Konovalov DA, Zenger KR: Genome-wide SNP validation and mantle tissue transcriptome analysis in the silver-lipped pearl oyster. Pinctada maxima. Mar Biotechnol. 2013, 15 (6): 647-658. 10.1007/s10126-013-9514-3.View ArticlePubMedGoogle Scholar
- Konovalov DA, Manning C, Henshaw MT: KINGROUP: a program for pedigree relationship reconstruction and kin group assignments using genetic markers. Mol Ecol Notes. 2004, 4 (4): 779-782. 10.1111/j.1471-8286.2004.00796.x.View ArticleGoogle Scholar
- Jones OR, Wang J: COLONY: a program for parentage and sibship inference from multilocus genotype data. Mol Ecol Resour. 2010, 10 (3): 551-555. 10.1111/j.1755-0998.2009.02787.x.View ArticlePubMedGoogle Scholar
- Garbe JR, Da YP: Pedigraph user manual version 2.4. 2008, Department of Animal Science, University of MinnesotaGoogle Scholar
- Sambrook J, Russell DW: Molecular cloning: a laboratory manual. 2001, New York: Cold spring harbor laboratory pressGoogle Scholar
- Hackett CA, Broadfoot LB: Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps. Heredity. 2003, 90 (1): 33-38. 10.1038/sj.hdy.6800173.View ArticlePubMedGoogle Scholar
- Green P, Falls K, Crooks S: Documentation for CRI-MAP. 1990, St Louis: Washington University, http://compgen.rutgers.edu/Crimap/Default.aspx,Google Scholar
- Liu X, Grosz M: Linkage mapping software user’s guide. 2006, Monsanto Company: St-LouisGoogle Scholar
- Keats BJB, Sherman SL, Morton NE, Robson EB, Buetow KH, Cartwright PE, Chakravarti A, Francke U, Green PP, Ott J: Guidelines for human linkage maps an international system for human linkage maps (ISLM, 1990). Ann Hum Genet. 1991, 55 (1): 1-6. 10.1111/j.1469-1809.1991.tb00392.x.View ArticlePubMedGoogle Scholar
- Kosambi D: The estimation of map distances from recombination values. Ann Eugen. 1944, 12: 172-175.View ArticleGoogle Scholar
- Liu RH, Meng JL: MapDraw: a microsoft excel macro for drawing genetic linkage maps based on given genetic linkage data. Yi Chuan. 2003, 25 (3): 317-321.PubMedGoogle Scholar
- de Givry S, Bouchez M, Chabrier P, Milan D, Schiex T: CarthaGène: multipopulation integrated genetic and radiation hybrid mapping. Bioinformatics. 2005, 21 (8): 1703-1704. 10.1093/bioinformatics/bti222.View ArticlePubMedGoogle Scholar
- Cartwright DA, Troggio M, Velasco R, Gutin A: Genetic mapping in the presence of genotyping errors. Genetics. 2007, 176 (4): 2521-2527. 10.1534/genetics.106.063982.PubMed CentralView ArticlePubMedGoogle Scholar
- Chakravarti A, Lasher LK, Reefer JE: A maximum likelihood method for estimating genome length using genetic linkage data. Genetics. 1991, 128 (1): 175-182.PubMed CentralPubMedGoogle Scholar
- Liao M, Zhang L, Yang G, Zhu M, Wang D, Wei Q, Zou G, Chen D: Development of silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) genetic maps using microsatellite and AFLP markers and a pseudo‒testcross strategy. Anim Genet. 2007, 38 (4): 364-370. 10.1111/j.1365-2052.2007.01617.x.View ArticlePubMedGoogle Scholar
- Danzmann R: LINKMFEX: LINKage analysis package for outcrossed families with male or female EXchange of the mapping parent, ver. 2.3. 2006, Ontario, Canada, http://www.uoguelph.ca/~rdanzman/software/LINKMFEX.zip,Google Scholar
- Sokal RR, Rohlf FJ: Biometry: the principles and practise of statistics in biological research. 2012, New York: W. H. Freeman and Co., 4Google Scholar
- Haldane J: Sex ratio and unisexual sterility in hybrid animals. J Genet. 1922, 12 (2): 101-109. 10.1007/BF02983075.View ArticleGoogle Scholar
- Ott J: Analysis of human genetic linkage (third edition). 1999, Baltimore: Johns Hopkins University PressGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B (Methodological). 1995, 57 (1): 289-300.Google Scholar
- Khatkar M, Nicholas F, Collins A, Zenger K, Cavanagh J, Barris W, Schnabel R, Taylor J, Raadsma H: Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC Genomics. 2008, 9 (1): 187-10.1186/1471-2164-9-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Hill WG, Robertson A: Linkage disequilibrium in finite populations. Theoret Appl Genet. 1968, 38 (6): 226-231. 10.1007/BF01245622.View ArticleGoogle Scholar
- Lewontin R: The interaction of selection and linkage. I. general considerations; heterotic models. Genetics. 1964, 49 (1): 49-PubMed CentralPubMedGoogle Scholar
- Abecasis GR, Cookson WOC: GOLD—Graphical overview of linkage disequilibrium. Bioinformatics. 2000, 16 (2): 182-183. 10.1093/bioinformatics/16.2.182.View ArticlePubMedGoogle Scholar
- Guo X, Li Q, Wang Q, Kong L: Genetic mapping and QTL analysis of growth-related traits in the pacific oyster. Mar Biotechnol. 2012, 14 (2): 218-226. 10.1007/s10126-011-9405-4.View ArticlePubMedGoogle Scholar
- Schwarz-Sommer Z, de Andrade Silva E, Berndtgen R, Lönnig W-E, Müller A, Nindl I, Stüber K, Wunder J, Saedler H, Gübitz T, et al: A linkage map of an F2 hybrid population of Antirrhinum majus and A. molle. Genetics. 2003, 163 (2): 699-710.PubMed CentralPubMedGoogle Scholar
- McGoldrick D, Hedgecock D: Fixation, segregation and linkage of allozyme loci in inbred families of the pacific oyster Crassostrea gigas (thunberg): implications for the causes of inbreeding depression. Genetics. 1997, 146 (1): 321-334.PubMed CentralPubMedGoogle Scholar
- Launey S, Hedgecock D: High genetic load in the pacific oyster Crassostrea gigas. Genetics. 2001, 159 (1): 255-265.PubMed CentralPubMedGoogle Scholar
- Li L, Xiang J, Liu X, Zhang Y, Dong B, Zhang X: Construction of AFLP-based genetic linkage map for zhikong scallop, Chlamys farreri Jones et Preston and mapping of sex-linked markers. Aquaculture. 2005, 245 (1–4): 63-73.View ArticleGoogle Scholar
- Wada H, Naruse K, Shimada A, Shima A: Genetic linkage map of a fish, the Japanese medaka Oryzias latipes. Mol Mar Biol Biotech. 1995, 4 (3): 269-274. 10.1007/BF02779020.View ArticleGoogle Scholar
- Zenger KR, McKenzie LM, Cooper DW: The first comprehensive genetic linkage map of a marsupial: the tammar wallaby (Macropus eugenii). Genetics. 2002, 162 (1): 321-330.PubMed CentralPubMedGoogle Scholar
- Hansson B, Åkesson M, Slate J, Pemberton JM: Linkage mapping reveals sex-dimorphic map distances in a passerine bird. Proc R Soc B Biol Sci. 2005, 272 (1578): 2289-2298. 10.1098/rspb.2005.3228.View ArticleGoogle Scholar
- Miles L, Isberg S, Glenn T, Lance S, Dalzell P, Thomson P, Moran C: A genetic linkage map for the saltwater crocodile (Crocodylus porosus). BMC Genomics. 2009, 10 (1): 339-10.1186/1471-2164-10-339.PubMed CentralView ArticlePubMedGoogle Scholar
- Lien S, Gidskehaug L, Moen T, Hayes B, Berg P, Davidson W, Omholt S, Kent M: A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns. BMC Genomics. 2011, 12 (1): 615-10.1186/1471-2164-12-615.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang Z, Parsch J: Positive correlation between evolutionary rate and recombination rate in drosophila genes with male-biased expression. Mol Biol Evol. 2005, 22 (10): 1945-1947. 10.1093/molbev/msi189.View ArticlePubMedGoogle Scholar
- Staelens J, Rombaut D, Vercauteren I, Argue B, Benzie J, Vuylsteke M: High-density linkage maps and sex-linked markers for the black tiger shrimp (Penaeus monodon). Genetics. 2008, 179 (2): 917-925. 10.1534/genetics.107.080150.PubMed CentralView ArticlePubMedGoogle Scholar
- Huxley JS: Sexual difference of linkage in Gammarus chevreuxi. J Genet. 1928, 20 (2): 145-156. 10.1007/BF02983136.View ArticleGoogle Scholar
- Singer A, Perlman H, Yan Y, Walker C, Corley-Smith G, Brandhorst B, Postlethwait J: Sex-specific recombination rates in zebrafish (Danio rerio). Genetics. 2002, 160 (2): 649-657.PubMed CentralPubMedGoogle Scholar
- Isberg SR, Johnston SM, Chen Y, Moran C: First evidence of higher female recombination in a species with temperature-dependent sex determination: the saltwater crocodile. J Hered. 2006, 97 (6): 599-602. 10.1093/jhered/esl035.View ArticlePubMedGoogle Scholar
- Ahmed M, Sparks AK: A preliminary study of chromosomes of two species of oysters (Ostrea lurida and Crassostrea gigas). J Fish Res Board Can. 1967, 24 (10): 2155-2159. 10.1139/f67-175.View ArticleGoogle Scholar
- Sakamoto T, Danzmann RG, Gharbi K, Howard P, Ozaki A, Khoo SK, Woram RA, Okamoto N, Ferguson MM, Holm L-E, et al: A microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) characterized by large sex-specific differences in recombination rates. Genetics. 2000, 155 (3): 1331-1345.PubMed CentralPubMedGoogle Scholar
- Franch R, Louro B, Tsalavouta M, Chatziplis D, Tsigenopoulos CS, Sarropoulou E, Antonello J, Magoulas A, Mylonas CC, Babbucci M, et al: A genetic linkage map of the hermaphrodite teleost fish Sparus aurata L. Genetics. 2006, 174 (2): 851-861. 10.1534/genetics.106.059014.PubMed CentralView ArticlePubMedGoogle Scholar
- Donis-Keller H, Green P, Helms C, Cartinhour S, Weiffenbach B, Stephens K, Keith TP, Bowden DW, Smith DR, Lander ES: A genetic linkage map of the human genome. Cell. 1987, 51 (2): 319-337. 10.1016/0092-8674(87)90158-9.View ArticlePubMedGoogle Scholar
- Berset-Brändli L, Jaquiéry J, Broquet T, Ulrich Y, Perrin N: Extreme heterochiasmy and nascent sex chromosomes in European tree frogs. Proc R Soc B Biol Sci. 2008, 275 (1642): 1577-1585. 10.1098/rspb.2008.0298.View ArticleGoogle Scholar
- Mank JE: The evolution of heterochiasmy: the role of sexual selection and sperm competition in determining sex-specific recombination rates in eutherian mammals. Genet Res. 2009, 91 (5): 355-10.1017/S0016672309990255.View ArticleGoogle Scholar
- Lynn A, Schrump S, Cherry J, Hassold T, Hunt P: Sex, not genotype, determines recombination levels in mice. Am J Hum Gen. 2005, 77 (4): 670-675. 10.1086/491718.View ArticleGoogle Scholar
- Hassold T, Hunt P: To err (meiotically) is human: the genesis of human aneuploidy. Nat Rev Genet. 2001, 2 (4): 280-291. 10.1038/35066065.View ArticlePubMedGoogle Scholar
- Tease C, Hultén MA: Inter-sex variation in synaptonemal complex lengths largely determine the different recombination rates in male and female germ cells. Cytogenet Genome Res. 2004, 107 (3–4): 208-215.View ArticlePubMedGoogle Scholar
- Weiss IM, Schönitzer V, Eichner N, Sumper M: The chitin synthase involved in marine bivalve mollusk shell formation contains a myosin domain. FEBS Lett. 2006, 580 (7): 1846-1852. 10.1016/j.febslet.2006.02.044.View ArticlePubMedGoogle Scholar
- Liu HL, Liu SF, Ge YJ, Liu J, Wang XY, Xie LP, Zhang RQ, Wang Z: Identification and characterization of a biomineralization related gene PFMG1 highly expressed in the mantle of Pinctada fucata. Biochemistry. 2007, 46 (3): 844-851. 10.1021/bi061881a.View ArticlePubMedGoogle Scholar
- Suzuki M, Saruwatari K, Kogure T, Yamamoto Y, Nishimura T, Kato T, Nagasawa H: An acidic matrix protein, Pif, is a key macromolecule for nacre formation. Science. 2009, 325 (5946): 1388-1390. 10.1126/science.1173793.View ArticlePubMedGoogle Scholar
- Kinoshita S, Wang N, Inoue H, Maeyama K, Okamoto K, Nagai K, Kondo H, Hirono I, Asakawa S, Watabe S: Deep sequencing of ESTs from nacreous and prismatic layer producing tissues and a screen for novel shell formation-related genes in the pearl oyster. PLoS One. 2011, 6 (6): e21238-10.1371/journal.pone.0021238.PubMed CentralView ArticlePubMedGoogle Scholar
- Metzler RA, Evans JS, Killian CE, Zhou D, Churchill TH, Appathurai NP, Coppersmith SN, Gilbert PUPA: Nacre protein fragment templates lamellar aragonite growth. J Am Chem Soc. 2010, 132 (18): 6329-6334. 10.1021/ja909735y.View ArticlePubMedGoogle Scholar
- Shen X, Belcher A, Hansma P, Stucky G, Morse D: Molecular cloning and characterization of lustrin A, a matrix protein from shell and pearl nacre of Haliotis rufescens. J Biol Chem. 1997, 272 (51): 32472-32481. 10.1074/jbc.272.51.32472.View ArticlePubMedGoogle Scholar
- Fan W, Hu Y, Li C, Xie L, Zhang R: Cloning, characterization, and expression analysis of calreticulin from pearl oyster Pinctada fucata. Tsinghua Sci Technol. 2008, 13 (4): 466-473. 10.1016/S1007-0214(08)70075-0.View ArticleGoogle Scholar
- Yano M, Nagai K, Morimoto K, Miyamoto H: A novel nacre protein N19 in the pearl oyster Pinctada fucata. Biochem Biophys Res Commun. 2007, 362 (1): 158-163. 10.1016/j.bbrc.2007.07.172.View ArticlePubMedGoogle Scholar
- Sonesson AK: Genomic selection for aquaculture: principles and procedures. Next generation sequencing and whole genome selection in aquaculture. Edited by: Liu Z. 2011, Ames, Iowa: Blackwell Publishing Ltd, 151-164.Google Scholar
- Du ZQ, Ciobanu DC, Onteru SK, Gorbach D, Mileham AJ, Jaramillo G, Rothschild MF: A gene-based SNP linkage map for pacific white shrimp. Litopenaeus vannamei. Anim Genet. 2010, 41 (3): 286-294. 10.1111/j.1365-2052.2009.02002.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.