|A consensus microsatellite-based linkage map for the hermaphroditic bay scallop (Argopecten irradians) and its application in size-related QTL analysis.|
|Jump to Full Text|
|PMID: 23077533 Owner: NLM Status: MEDLINE|
|Bay scallop (Argopecten irradians) is one of the most economically important aquaculture species in China. In this study, we constructed a consensus microsatellite-based genetic linkage map with a mapping panel containing two hybrid backcross-like families involving two subspecies of bay scallop, A. i. irradians and A. i. concentricus. One hundred sixty-one microsatellite and one phenotypic (shell color) markers were mapped to 16 linkage groups (LGs), which corresponds to the haploid chromosome number of bay scallop. The sex-specific map was 779.2 cM and 781.6 cM long in female and male, respectively, whereas the sex-averaged map spanned 849.3 cM. The average resolution of integrated map was 5.9 cM/locus and the estimated coverage was 81.3%. The proportion of distorted markers occurred more in the hybrid parents, suggesting that the segregation distortion was possibly resulted from heterospecific interaction between genomes of two subspecies of bay scallop. The overall female-to-male recombination rate was 1.13:1 across all linked markers in common to both parents, and considerable differences in recombination also existed among different parents in both families. Four size-related traits, including shell length (SL), shell height (SH), shell width (SW) and total weight (TW) were measured for quantitative trait loci (QTL) analysis. Three significant and six suggestive QTL were detected on five LGs. Among the three significant QTL, two (qSW-10 and qTW-10, controlling SW and TW, respectively) were mapped on the same region near marker AiAD121 on LG10 and explained 20.5% and 27.7% of the phenotypic variance, while the third (qSH-7, controlling SH) was located on LG7 and accounted for 15.8% of the phenotypic variance. Six suggestive QTL were detected on four different LGs. The linkage map and size-related QTL obtained in this study may facilitate marker-assisted selection (MAS) in bay scallop.|
|Hongjun Li; Xiao Liu; Guofan Zhang|
Related Documents :
|16870283 - Chromosome visualisation in filamentous fungi.
6156113 - Association frequency and silver staining of nucleolus organizing regions in hyperthyro...
3536873 - Fluphenazine-resistant saccharomyces cerevisiae mutants defective in the cell division ...
|Type: Journal Article; Research Support, Non-U.S. Gov't Date: 2012-10-16|
|Title: PloS one Volume: 7 ISSN: 1932-6203 ISO Abbreviation: PLoS ONE Publication Date: 2012|
|Created Date: 2012-10-18 Completed Date: 2013-03-28 Revised Date: 2013-05-20|
Medline Journal Info:
|Nlm Unique ID: 101285081 Medline TA: PLoS One Country: United States|
|Languages: eng Pagination: e46926 Citation Subset: IM|
|Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China.|
|APA/MLA Format Download EndNote Download BibTex|
Pectinidae / anatomy & histology*, genetics*
Quantitative Trait Loci*
Journal ID (nlm-ta): PLoS One
Journal ID (iso-abbrev): PLoS ONE
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
Publisher: Public Library of Science, San Francisco, USA
Copyright: 2012 Li et al
Received Day: 1 Month: 5 Year: 2012
Accepted Day: 6 Month: 9 Year: 2012
collection publication date: Year: 2012
Electronic publication date: Day: 16 Month: 10 Year: 2012
Volume: 7 Issue: 10
E-location ID: e46926
PubMed Id: 23077533
Publisher Id: PONE-D-12-12475
|A Consensus Microsatellite-Based Linkage Map for the Hermaphroditic Bay Scallop (Argopecten irradians) and Its Application in Size-Related QTL Analysis Alternate Title:Linkage Map and QTL Analysis for Bay Scallop|
|Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China
|Biomedical Research Institute, United States of America
|Correspondence: * E-mail: firstname.lastname@example.org
[conflict] Competing Interests: The authors have declared that no competing interests exist.
Contributed by footnote: Conceived and designed the experiments: XL. Performed the experiments: HL. Analyzed the data: HL. Contributed reagents/materials/analysis tools: GZ. Wrote the paper: HL XL.
¤Current address: National Marine Environmental Monitoring Center, Dalian, China
Linkage maps have become important genetic tools which provide detailed information on genotype-phenotype relationships . Quantitative trait loci (QTL) analysis based on linkage map reveals the genetic basis and inheritance patterns of economically important traits, which is essential for marker-assisted selection (MAS) in future breeding programs. Genetic linkage maps have been published for some economically important aquaculture species such as the Asian seabass (Lates calcarifer) , the Atlantic salmon (Salmo salar) , the Japanese flounder (Paralichthys olivaceus) , the rainbow trout (Oncorhynchus mykiss) , the Pacific white shrimp (Litopenaeus vannamei)  and the tilapia (Oreochromis spp.) . Although marine bivalves represent a major portion in the aquaculture industry, only a few linkage maps have been constructed mainly using amplified fragment length polymorphisms (AFLPs). These AFLP-based linkage maps include the oysters (Crassostrea virginica, Crassostrea gigas and Ostrea edulis), the pearl oysters (Pinctada martensii and Pinctada fucata), the mussel Mytilus edulis and the scallops (Chlamys farreri,  and Argopecten irradians, ).
Many economically important traits of aquaculture species are continuous quantitative in phenotype, which are controlled by multiple loci across the whole genome. The identification and mapping of QTL controlling economic traits represent an important approach towards determining the molecular genetic basis of some complex traits, especially the size-related QTL in mollusks. Although QTL analysis in mollusk may have been limited by the lack of inbred lines, F2 populations and backcross families, the high fecundity and high levels of polymorphisms among different populations of mollusks may facilitate QTL mapping in these animals . A few QTL controlling growth and disease resistance traits have been located on some linkage maps in marine bivalve such as the oysters – and the scallops , , which will be useful for MAS in these aquaculture species.
Bay scallop (Argopecten irradians) is a simultaneous hermaphroditic marine bivalve, with uncontinuous distribution along the eastern coast of United States and Gulf of Mexico. It can be further divided into three geographical subspecies: A. i. irradians, A. i. concentricus and A. i. amphicostus. Since early 1980s, two subspecies, A. i. irradians and A. i. concentricus have been successively introduced to China and quickly gained popularity in Northern and Southern China respectively . The annual production of bay scallops in China increased steadily from over 200,000 tons in 1990s  to 600,000 tons in 2007 . Genetic improvement has contributed, at least partially, to the great success of bay scallop aquaculture in China . The efficiency of genetic improvement depends largely on the availability of genetic information such as linkage maps. Although we have developed sex-specific linkage maps for bay scallop mainly using AFLP markers , such dominant markers may not be easily shared among different mapping populations established by different laboratories. Among the currently available molecular markers, microsatellite markers are the most popular choice for linkage mapping. Microsatellites are short (1–6 bp) repetitive DNA sequences, which are highly abundant and evenly distributed throughout genomes , . As codominant markers, microsatellite data may transfer well among different mapping populations and will be necessary for future QTL studies among different studies. In this study, we present a consensus genetic map using microsatellite markers for bay scallop, which may provide a scaffold to enable integration with novel markers.
The mapping panel used in this study included two hybrid backcross-like families constructed with two geographical subspecies of bay scallop: the northern bay scallop A. i. irradians (denoted as N) and the southern bay scallop A. i. concentricus (denoted as S). The N stock was initially introduced from Massachusetts and Virginia (USA) in 1998 and 1999 , and has been cultured mainly in Bohai Bay in Northern China. The brood stocks of N stock used in this study have been selected for orange shell color for 3 generations . The brood stocks of S stock, whose shell color is black on the left valve and white on the right, were obtained from a hatchery stock cultivated in Guangdong province since its initial introduction from Florida (USA) in 1991 . Two hybrid families (♀N×♂S and ♀S×♂N) were produced by crossing the N and S stocks in a hatchery in Qingdao. Two backcross families (Figure 1) were constructed for mapping: (1) CC5: a crossbred family between a ♂ hybrid (♀S×♂N) and a ♀S, denoted as ♀S×♂(S×N); (2) CC10: a crossbred family between a ♂ hybrid (♀N×♂S) and a ♀S, denoted as ♀S×♂(N×S). As these crossbred families were not derived from inbred lines of Northern and Southern bay scallop, these crosses should be considered as “backcross-like”. The CC5 and CC10 families were raised under the same culture condition following the protocols described by Zheng et al. . Briefly, fertilized eggs of each family were raised separately in a 40 L polyethylene tank for hatching at 22–23°C with a salinity of 30–32‰. Stocking density was kept at 10 individuals/mL by adjusting seawater volume during larval culture. Spats were kept in large concrete tanks until they reached 500∼600 µm. Thereafter, scallops were placed in polyethylene bags and transferred to an outdoor nursery pond. Seventy-day-old juveniles were dispersed into lantern nets (10 layers per net) and hung on a long-line system for grow-out. The densities were adjusted monthly, from 200 juveniles per layer at an early stage to 30 adults per layer after about day 130. At the age of eight months, the offspring (N = 94 for each family) were sampled and their shell length (SL), shell height (SH), shell width (SW) and total weight (TW) were measured and recorded. The adductor muscles were dissected and stored in 80% ethanol. The shell color was determined according to the color system developed by Elek and Adamkewicz  in this study. As previous studies indicated that shell color in the bay scallop is under single-gene control where white color is recessive , , the genotypes for orange-colored and white-colored parent were recorded as heterozygote and homozygote respectively.
DNA from each individual of the two mapping families was extracted from 100 mg of adductor muscle using the traditional phenol/chloroform method. A total of 268 microsatellite markers evaluated and screened in this study were as follows: eight from Roberts et al. , 20 suffixed with AIMS from Zhan et al. , ; 28 from Wang et al.  and 212 from Li et al. –. Amplification of microsatellite DNA sequences by polymerase chain reaction (PCR) was performed in a 15 µL reaction system containing approximately 50 ng template DNA, 1×PCR buffer, 1.5 mM MgCl2, 200 µM dNTPs, 5 pmol of each primer and 0.5 U Taq DNA polymerase (Promega). PCR reaction profile was set as follows: denaturation for 5 min at 95°C, followed by 35 cycles of 30 s at 94°C, 30 s at the optimal primer-specific annealing temperature and 30 s at 72°C, with a final extension step at 72°C for 10 min. PCR products were separated using 12% non-denaturing polyacrylamide on a vertical PAGE system AE-6220 (ATTO, Japan). Gels were run at 300 V for 2–3 h depending on amplified fragment sizes, stained with ethidium bromide (EB) and visualized under ultraviolet light.
Polymorphic microsatellite markers derived from each parent were scored as co-dominant markers and a letter (A, B, C or D) was assigned for each allele. Prior to linkage analysis, segregation of alleles was examined for goodness-of-fit to the expected segregation ratio using the chi-square test at α = 0.05 level of confidence. The significant criteria were adjusted for multiple comparisons based on the number of linkage groups (LGs)  using the sequential Bonferroni correction . Segregating markers were categorized into four expected segregation types (including null-alleles): 1∶1∶1∶1 type (♀×♂: AB×CD and AB×AC), 1∶2∶1 type (AB×AB), 1∶1 ♀ type (AB×AA or CC), and 1∶1 ♂ type (AA or CC×AB). Informative markers including shell color were grouped into separate sets corresponding to each parent and analyzed independently to construct separate sex-specific framework maps using the double pseudotestcross mapping strategy. The software JoinMap 3.0  was used for map construction from each mapping family with the cross-pollinating (CP) population type, which represents a population resulting from a cross between two heterogeneous parents with linkage phases unknown. Clustering of markers was performed by calculating logarithm of odds (LOD) scores for recombination rate (θ) based on the G2 statistic test, which is not affected by segregation distortion , thus leading to less incidence of spurious linkage. A LOD threshold >3.5 and a recombination rate threshold <0.4 were set to obtain maps from each mapping population. After assigning the markers into respective LG, heterogeneity in recombination rate for each pair of markers was tested between sexes and families using a G test, which compared the observed numbers of recombinants and non-recombinants with the mean expected numbers.
After creating the individual linkage maps, LGs with common markers on individual maps were merged to create a composite map using “map integration” function of JoinMap. If individual maps produced no significant different recombination rate for all combinations of orthologous markers in certain groups, the integrated LGs were constructed using a LOD weighted average of recombination rate values. When marker order in certain LGs was disturbed between individuals, one parental LG was fixed as a framework and “fix order” command was used to identify the most likely marker order. For assessing the differences in recombination rate between sexes, comparisons of recombination differences between both parents in two mapping families were performed by analyzing all pairwise marker combinations using a two-way contingency G-test as implemented in JoinMap. Graphics of the LGs were constructed using MapChart 2.1 .
On the basis of the consensus linkage map, two methods were used to estimate the expected genome length. First, we calculated the average spacing s between markers by dividing the total observed map length of all LGs by the number of marker intervals (number of markers minus number of LGs). We estimated genome length (Ge1) by adding 2 s to the length of each LG to account for terminal chromosome regions . Second, we calculated estimated genome length (Ge2) by multiplying the length of each LG by (m+1)/(m−1), where m is the number of framework markers on each LG (method 4 in ). The average of the two estimates was used as the estimated genome length (Ge) to describe the genome length. The observed genome length (Go) was calculated as the total length considering all markers. The genome coverage was determined by Go/Ge.
Mean value, standard deviations, and the Pearson's correlation analysis of four size-related traits including shell length (SL), shell height (SH), shell width (SW) and total weight (TW) were analyzed with SPSS 12.0. QTL detection was performed by MapQTL 5.0 . Interval mapping (IM) and multiple QTL model (MQM) mapping were utilized to detect any significance associating size-related traits and marker loci. The LOD score significance thresholds were calculated by permutation tests in MapQTL at genome-wide (significant linkage, α<0.05, n = 1000) and chromosome-wide significant level (suggestive linkage, α<0.05, n = 1000). Interval mapping was performed initially using the default parameters of MapQTL. Map intervals exceeding the linkage-group-wide threshold from interval mapping were subsequently selected as cofactors for the MQM procedure. Kruskal-Wallis (KW) analysis (a single marker analysis) was also tested based on one-way analysis of variance (ANOVA), using four parental allelic combinations (m1f1, m1f2, m2f1, m2f2) of markers nearest each QTL peak. Association analysis was carried out to determine the QTL allele substitution effects.
Of a total of 268 microsatellite markers and one phenotypic marker (shell color) screened, 169 markers were informative in at least one family. The number of segregating markers was essentially similar between the two families, with 137 in family CC5 and 119 in family CC10. Within family CC5, 96 and 95 markers segregated in the female and male parents respectively, and 54 markers exhibited a segregation ratio of 1∶1∶1∶1, providing the most valuable co-dominant mapping information. For family CC10, the female and male parent contained 75 and 97 markers respectively, and 53 markers were informative in both parents. Among these markers, 40 were commonly shared with polymorphism across all the parents, and thus allowed integration of two individual genetic linkage maps as allelic bridges.
Null-allele segregation, without which unexpected offspring genotypes could not be explained, was estimated to be 25 loci (25/137, 18.2%) in family CC5 and 22 loci (22/119, 18.5%) in family CC10. Null alleles are not amplified by RCR, and may cause genotyping errors of heterozygote deficiency in population genetic studies . However, null alleles are not a problem for mapping in pedigreed populations . In family-segregation analyses, null alleles can be regarded as recessive alleles for mapping purpose . Thus, we included the informative null-allele segregation for map construction.
Chi-square tests were used to compare all observed progeny ratios of each marker against the genotypic ratios expected under a 1∶1 Mendelian segregation ratio. In family CC5, distorted segregation was observed at a total of 21 loci (15.3%), including 5 loci segregating in the female and 16 loci in the male parent (P<0.05). In family CC10, a total of 23 loci (19.3%), including 7 loci in the female and 16 loci in the male parent, showed significant segregation distortion (P<0.05). The distribution of distorted markers was asymmetrical between the two sexes, with higher incidence of segregation distortion in male parents than female. Deviations from Mendelian expectations implicated several homologous loci and LGs across families. It is noteworthy that four markers (AiAD154, AiAD040, AiDD105 and AiSD274) linked to each other on the homologous group between the two families. Apart from these observations, linked skewed markers appeared scattered across different LGs between families. After sequential Bonferroni corrections, only five deviations showed significant segregation distortions in two mapping families.
Based on the result of linkage grouping at LOD = 3.5, we constructed an individual sex-specific linkage map for each parent. In family CC5, the female and male maps included 96 and 95 makers on 16 and 18 LGs respectively. The female map covered 739.1 cM, with an average interval of 9.2 cM and the number of markers varied from 2 to 13 with an average of 6 markers per LG. The male map covered 696.0 cM, with an average interval of 9.0 cM. The number of loci per LG varied from 2 to 11 with a mean of 5.3 markers. In family CC10, the female and male maps had 75 and 97 markers on 17 and 16 LGs respectively. These two linkage maps covered 657.1 cM and 641.8 cM with an average of 4.4 and 6.0 markers per LG and an averaging spacing of 11.3 cM and 7.9 cM. The details of individual linkage maps, including length, number of markers, marker spacing and largest interval are shown in Table S1.
Bridge markers or homologous loci were used to identify co-linear regions between different individual maps. As shown in Figure 2, the integrated map contained 161 markers distributed over 16 LGs, which was consistent with the haploid chromosome number of the bay scallop . This map covered 849.3 cM with an average inter-marker spacing of 5.9 cM (Table 1). The length of LGs ranged from 24.3 cM to 77.9 cM and the number of markers varied from 4 to 16 per LG (Table 1). The marker for the phenotype of shell color was located on LG10, which was flanked by two microsatellite markers AiSD245 and AIMS013, with the nearest genetic distance of 0.5 cM.
The genome length of the integrated map estimated by two methods was 1038.1 cM and 1050.7 cM, respectively, with an average of 1044.4 cM. Based on the observed length of the frame map (849.3 cM) and the average expected genome length (1044.4 cM), the percentage of the genome covered by the integrated map was 81.3%.
The availability of the common microsatellite markers in the different individual maps allowed for a comparative analysis of recombination rate between sexes and families. Although slight difference in length was observed between female and male maps within the two families (Table S1), a higher was observed in female parents when only common informative markers between sexes were considered. The overall female-to-male recombination rate was 1.13∶1 across all linked markers that are common to both parents in two families. There were 141 and 178 comparable marker pairs in the “between-sexes within each family” category and “within-sexes between families” category, respectively. Of 141 tests of heterogeneity of recombination rate between sexes, 32 (22.7%) were significant at α = 0.05 level using G test. The proportion of significant difference in recombination rate was 20.2% (36/178). When the percentage of marker pairs with significant heterogeneity in recombination rate was compared between these two categories, there was no significant difference (Mann–Whitney test, P>0.05).
Descriptive statistics of the trait distributions were presented in Table S2. The Kolmogorov–Smirnov test indicated that the variances of SW and TW in the family CC5, and SW, SH and TW in family CC10 were coordinate with normal distribution, however, SL and SW in the family CC5, and SL in the family CC10 were not normally distributed. Phenotypic correlations estimated between any two traits were significant (P<0.01), with the Pearson correlation coefficients ranging from 0.772 to 0.941.
QTL affecting four size-related traits, including SL, SH, SW and TW were scanned on a genome scale. Permutation tests using MapQTL indicated that the genome-wide LOD significance thresholds were 3.1–3.3 for SL, SH, SW and TW, whereas the chromosome-wide LOD significance thresholds varied from 2.4 to 2.6. Three significant (genome-wide) and six suggestive (linkage-group-wide) QTL were detected on 6 LGs (LG5, LG6, LG7, LG10 and LG13) in the two mapping populations using the MQM method (Table 2 and Figure 3). The three significant QTL (qSH-7, qSW-10 and qTW-10) controlling shell height, shell width and total weight were located on LG7 and LG10, explaining 15.8%, 20.5% and 27.7% of the observed phenotypic variance, respectively. qSW-10 and qTW-10 distributed on the same region near marker AiAD121. Six suggestive QTL were detected on four different LGs. These QTL explained relatively low phenotypic variance, ranging from 8.9% to 15.3%. Considering QTL positions and confidence intervals, some QTL for SL, SH, SW and TW (qSL-7 and qSH-7, qSW-10 and qTW-10, qSL-13 and qTW-13) were mapped to the same locations of LG7, LG10 and LG13, whereas other QTL were mapped on different LGs. In both CC5 and CC10 families, common QTL were identified on LG10.
KW test was used for single marker to test associations between phenotypic traits and genotypes of marker nearing each QTL (Table 2). Except for qTW-6, all significant and suggestive QTL showed significance by KW analysis. Even though qTW-6-nearest maker AiAD079 showed no significance by KW analysis, other loci flanking this QTL showed significance. In several cases, the associated maker was not always the nearest-QTL one.
As with most other marine invertebrates, mapping resources such as purebred lines, near-isogenic lines and inbred lines were not available in A. irradians owing to the biological characteristics including one-year life cycle, high early mortality and inbreeding depression. Since A. irradians contains high levels of polymorphisms among local populations , we applied an F1 pseudo-testcross mapping strategy to take advantages of parental heterozygosity . Considering the co-dominance and easy transferability of microsatellites, two families were used for linkage map construction, which may increase the possibility of placing more microsatellite markers on a consensus map. In addition, the availability of orthologous markers also allowed for a direct comparison of recombination rates and gene orders among different parents. Although moderately dense linkage maps have been previously reported for bay scallop using AFLPs , , this integrated linkage map based on microsatellite DNA markers represents an important development for marine invertebrate genome research. When compared with the AFLP-based maps, the microsatellite-based linkage map has the advantage of high portability to other families or strains, which may allow more loci to be placed on a genome-wide anchored microsatellite map. Therefore, this study focused on obtaining a general order and distance among these markers rather than the fine resolution of order and distance. The microsatellite-based map may provide a good resource from which markers may be selected for future mapping projects in bay scallop and for comparative studies among other scallops.
Deviations from Mendelian segregation ratio have been detected in previous efforts to construct genetic linkage maps for marine bivalves , . For relative scarce linkage maps based on microsatellites, distorted segregation after sequential Bonferroni correction accounted for 30.45% and 29.96% of loci in two mapping families of Zhikong scallop . In this study, the frequency of distorted marker was 15.3% in family CC5 and 19.3% in family CC10. A similar result was also observed in the construction of linkage maps using AFLP markers for bay scallop (14.5% segregation distortion ). The level of segregation distortion varies greatly among species and likely reflects the unique characteristics of a genome. Many factors may give rise to segregation distortion. In the Pacific oyster, this phenomenon was explained by a large load of recessive deleterious mutations and low survival of identical-by-descent homozygotes, which was probably caused by deleterious recessive lethal genes . There existed evidence that selective mortality mainly occurred around the time of metamorphosis from the larval to the juvenile stage . High genetic load was confirmed by Hubert and Hedgecock , who constructed a Pacific oyster linkage map using 11-day-old larvae, which proved to be effective in reducing segregation distortion. In addition, segregation distortion may result from other factors, such as close linkage of markers to genes or chromosomal regions affecting gametogenesis, competition among gametes for preferential fertilization, genetic drift and chromosomal rearrangements , . In this study, it is noteworthy that the proportion of distorted markers was asymmetrical between the sexes, with more occurred in male (hybrid) parents. This result suggested the hypothesis that segregation distortion resulted from heterospecific interaction between genomes of two subspecies of bay scallop and semi-lethal effects may be more evident in hybrid individuals , . The masked recessive lethal genes in hybrid parents should be unveiled in backcrossing, which makes more segregation distortions observed in hybrid parents.
Great diversity in shell color has been observed in many marine mollusks, and the variations give considerable implications to breeders since they are associated with differences in survival, growth and other adaptive characteristics , . Previous studies indicated that variations in shell color were influenced by environmental factors , , but in many cases they were genetically determined. In this study, shell color inheritance pattern in CC5 and CC10 families showed a simple 1∶1 segregation, suggesting that shell color in bay scallop is under genetic control. This result is in agreement with other mollusks, such as the mussel M. edulis, the Pacific abalone Haliotis discus and the Zhikong scallop C. farreri, lending support to the single-gene two-allele model for genetic control of shell color. The fact that white× (white×orange) crosses produced about 50% white and 50% orange progeny implies that orange parent is heterozygous and white allele is recessive to orange. In our previous study, we have detected an association between orange color and higher quantitative traits in larval and adult stages of bay scallop . The microsatellite markers linked to shell color obtained in this study provide an opportunity to breed colorful lines with enhanced growth performance for bay scallop.
It is of interest to examine the difference of recombination rate between sexes and individuals for any species. Significant difference in sex recombination ratios has been reported in many aquatic animals . In this study, the recombination rate in female parents was higher than the males, with an overall female-to-male ratio of 1.13∶1. This is consistent with higher female recombination rates and longer map lengths observed in other marine mollusks, such as the oyster (Crassostrea gigas, , ) and the abalone (Haliotis rubra, F∶M = 1.45∶1 ). Difference in recombination rate between sexes may result from various complex factors . Since bay scallops are functional hermaphrodites that exhibit both male and female reproductive functions in a single individual, it should be possible to observe recombination in the same individual acting as a male or a female parent. Further investigation is needed to understand the sex specificity in recombination rate of the hermaphroditic bay scallop.
Family-specific recombination rates were also demonstrated in linkage maps based on multiple mapping populations. Significant differences in recombination and orders of marker were evident among same-sex parents of different families in C. gigas, which suggested that polymorphism for chromosomal rearrangements may exist in natural populations . Similarly, Zhan et al. attributed the heterogeneity in recombination rate to individual difference, not sexes. In this study, there was no significant difference in the percentage of marker pairs with significant heterogeneity in recombination rate between “within-sexes between families” and “between-sexes within each family” category, which suggests that heterogeneities in recombination rates would be mainly due to individual differences. Although chromosomal rearrangements , sampling variations associated with sample sizes and statistical treatments , and environmental differences between mapping families  may be possible causes for this common phenomenon, the underlying mechanism is still not well understood.
Accurate mapping techniques have been developed in recent years to detect QTL for economically important traits. The efficiency for QTL mapping depends on heritability of traits, average allelic substitution effect of the alleles involved, recombination distance between the QTL and associated markers, and sample size of the progeny . For most shellfish species, they are broadcast spawners with high fecundity. The large number of progeny produced in a single family can ensure the accuracy of linkage map construction and QTL mapping. In addition, most bivalve mollusks possess high levels of heterozygosity. Hence, pseudo-testcross strategy would be ideal for QTL analysis in marine bivalves , . However, QTL mapping study using crosses between heterozygous parents from a breeding population has limited QTL discrimination. As QTL allele-substitution is measured against a heterogeneous background, the detected QTL may be limited to loci with specific genetic configurations and/or strong effects. Moreover, QTL are environmentally influenced and it is difficult to dissect the genetic and environment effects , . As a result, QTL-associated markers need to be validated in different genetic backgrounds and environments for application in molecular breeding .
In this study, four size-related traits were highly correlated with each other, with the Pearson's correlation coefficient being greater than 0.772. Separate QTL analysis showed that LOD peaks controlling two size-related traits were located on the same region of LG7, LG10 and LG13, indicating either the linkage of two QTL or the presence of a single QTL with pleiotropic effects. Common QTL-underlying phenotypes may be the most likely explanation to the co-location of LOD peaks . However, given the current marker density and unavailability of selected lines, it is rather difficult to distinguish between a single gene (or gene cluster) with individual QTL peak. More markers are needed to identify the full spectrum of genetic components contributing to size differences and to narrow down the distance between a marker and a QTL peak.
Three significant and six suggestive QTL were mapped on five LGs in this study, while in Zhikong scallop, a total of eight size-related QTL were detected on nine different LGs using a genome scan . The proportion of variance explained by the significant QTL varied from 15.8% to 27.7%, which is similar to the results of 6.4–19.2% found in Zhikong scallop . Among the three significant QTL identified in this study, qSW-10 and qTW-10 located on LG10, explained high percentage of the observed phenotypic variance, suggesting that there may exist a few QTL/genome regions that control size-related traits of the bay scallop. Our findings agree well with some previous studies in crops and domestic animals that a large proportion of quantitative variation can be explained by the segregation of a few major QTL .
In this study, two mapping families produced considerably different results, with only QTL on LG10 shared between them. Similar results were also obtained from the QTL analysis with multiple families in plants and insects. In maize, two independent samples evaluated in four environments were used to detect 107 QTL, of which only 20 were in common . Previous studies indicated that family-specific QTL could be caused by statistical sampling and fixation of some QTL in some families , . Since bay scallops are highly fecund animals with intensive natural selection on early life stages, biased survival of progeny may contribute to the different QTL results of different families. Additionally, we used hatchery scallops as mapping parents, which may be fixed for homozygous QTL allele during aquaculture practice. No segregation occurred in mapping parents. This study demonstrates the need of using multiple families for QTL mapping analysis in bay scallop.
A one-way ANOVA revealed that allelic substitution of microsatellite near significant and suggestive QTL showed large effects on size-related traits. The genotypes at these microsatellite loci may be useful for growth improvement through MAS. It could be expected that only the closest marker to the QTL peak showed a significant association, but in several cases, the associated marker was not always the closest one. There could be multiple explanations for this observation, including low mapping information of the closest marker, large genome region in linkage disequilibrium with QTL and possible existence of a secondary segregating QTL . Future studies focusing on genotyping more markers would be valuable for improving the resolution of QTL mapping.
In this study, a consensus genetic linkage map was constructed for the bay scallop using microsatellite markers, providing a base-line map for this important aquaculture species. It offers a scaffold that can be applied for further genetic studies, such as increasing marker density, QTL analysis of economic traits, genome assembly and addressing the problems of scallop evolution. Size-related QTL analysis was performed in two backcross-like families and nine QTL were detected on five LGs. These results provide a good starting point for fine-mapping of the identified QTL for MAS and eventually identifying candidate genes responsible for growth-related traits for gene-assisted selection (GAS) in the bay scallop.
Statistics for individual sex-specific linkage maps of bay scallop (Argopecten irradians) in two reference families CC5 and CC10.
Click here for additional data file (pone.0046926.s001.doc)
Statistics of four size-related traits in two mapping families CC5 and CC10 (± standard deviation).
Click here for additional data file (pone.0046926.s002.doc)
|1.||Wu R, Ma C, Casella G (2007) Statistical genetics of quantitative traits: linkage, maps and QTL. Springer, New York, NY, USA.|
|2.||Wang C,, Bai Z,, He X,, Lin G,, Xia J,, et al. (Year: 2011) A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer. BMC Genomics12: 17421457569|
|3.||Lien S,, Gidskehaug L,, Moen T,, Hayes B,, Berg P,, et al. (Year: 2011) A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns. BMC Genomics12: 61522182215|
|4.||Castano-Sanchez C,, Fuji K,, Ozaki A,, Hasegawa O,, Sakamoto T,, et al. (Year: 2010) A second generation genetic linkage map of Japanese flounder (Paralichthys olivaceus). BMC Genomics11: 55420937088|
|5.||Guyomard R,, Boussaha M,, Krieg F,, Hervet C,, Quillet E, (Year: 2012) A synthetic rainbow trout linkage map provides new insights into the salmonid whole genome duplication and the conservation of synteny among teleosts. BMC Genomics13: 1522233149|
|6.||Du ZQ,, Ciobanu DC,, Onteru SK,, Gorbach D,, Mileham AJ,, et al. (Year: 2010) A gene-based SNP linkage map for pacific white shrimp, Litopenaeus vannamei. Animal Genetics41: 286–29419968647|
|7.||Lee BY,, Lee WJ,, Streelman JT,, Carleton KL,, Howe AE,, et al. (Year: 2005) A second-generation genetic linkage map of tilapia (Oreochromis spp.). Genetics170: 237–24415716505|
|8.||Yu Z,, Guo X, (Year: 2003) Genetic linkage map of the eastern oyster Crassostrea virginica Gmelin. The Biological Bulletin204: 327–33812807709|
|9.||Li L,, Guo X, (Year: 2004) AFLP-based genetic linkage maps of the Pacific oyster Crassostrea gigas Thunberg. Marine Biotechnology6: 26–3614564534|
|10.||Lallias D,, Beaumont AR,, Haley CS,, Boudry P,, Heurtebise S,, et al. (Year: 2007) A first-generation genetic linkage map of the European flat oyster Ostrea edulis (L.) based on AFLP and microsatellite markers. Animal Genetics38: 560–56817931405|
|11.||Shi Y,, Kui H,, Guo X,, Gu Z,, Wang Y,, et al. (Year: 2009) Genetic linkage map of the pearl oyster, Pinctada martensii (Dunker). Aquaculture Research41: 35–44|
|12.||Yu D,, Wang X,, Guo Y,, Huang G,, Gong S,, et al. (Year: 2007) Construction of genetic linkage maps in pearl oyster Pinctada fucata using AFLP markers. Journal of Fishery Sciences of China14: 361–368 (Chinese edition with English abstract).|
|13.||Lallias D,, Lapegue S,, Hecquet C,, Boudry P,, Beaumont AR, (Year: 2007) AFLP-based genetic linkage maps of the blue mussel (Mytilus edulis). Animal Genetics38: 340–34917559551|
|14.||Wang L,, Song L,, Chang Y,, Xu W,, Ni D,, et al. (Year: 2005) A preliminary genetic map of Zhikong scallop (Chlamys farreri Jones et Preston 1904). Aquaculture Research36: 643–653|
|15.||Li L,, Xiang J,, Liu X,, Zhang Y,, Dong B,, et al. (Year: 2005) Construction of AFLP-based genetic linkage map for Zhikong scallop, Chlamys farreri Jones et Preston and mapping of sex-linked markers. Aquaculture245: 63–73|
|16.||Qin Y,, Liu X,, Zhang H,, Zhang G,, Guo X, (Year: 2007) Genetic mapping of size-related quantitative trait loci (QTL) in the bay scallop (Argopecten irradians) using AFLP and microsatellite markers. Aquaculture272: 281–290|
|17.||Wang L,, Song L,, Zhang H,, Gao Q,, Guo X, (Year: 2007) Genetic linkage map of bay scallop, Argopecten irradians irradians (Lamarck 1819). Aquaculture Research38: 409–419|
|18.||Zhan A,, Hu J,, Hu X,, Hui M,, Wang M,, et al. (Year: 2009) Construction of microsatellite-based linkage maps and identification of size-related quantitative trait loci for Zhikong scallop (Chlamys farreri). Animal Genetics40: 821–83119515089|
|19.||Yu Z,, Guo X, (Year: 2006) Identification and mapping of disease resistance QTLs in the Eastern oyster, Crassostrea virginica Gmelin. Aquaculture254: 160–170|
|20.||Lallias D,, Gomez-Raya L,, Haley C,, Arzul I,, Heurtebise S,, et al. (Year: 2009) Combining two-Stage testing and interval mapping strategies to detect QTL for resistance to bonamiosis in the European flat oyster (Ostrea edulis). Marine Biotechnology11: 570–58419139958|
|21.||Sauvage C,, Boudry P,, De Koning DJ,, Haley CS,, Heurtebise S,, et al. (Year: 2010) QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (Crassostrea gigas). Animal Genetics41: 390–39920096029|
|22.||Waller TR, (Year: 1969) The Evolution of the Argopecten gibbus stock (Mollusca: Bivalvia) with emphasis on the tertiary and quaternary species of Eastern North America. Memoir (The Paleontological Society)43 (Supplement) 1–121|
|23.||Zhang F,, He Y,, Yang H, (Year: 2000) Introduction engineering of bay scallop and its comprehensive effects. Engineering Science2: 30–35 (Chinese edition with English abstract).|
|24.||Guo X,, Ford SE,, Zhang G, (Year: 1999) Molluscan aquaculture in China. Journal of Shellfish Research18: 19–31|
|25.||Guo X, (Year: 2009) Use and exchange of genetic resources in molluscan aquaculture. Reviews in Aquaculture1: 251–259|
|26.||Weber JL,, May PE, (Year: 1989) Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction. American Journal of Human Genetics44: 388–3962916582|
|27.||Li H,, Liu W,, Gao X,, Zhu D,, Wang J,, et al. (Year: 2011) Identification of host-defense genes and development of microsatellite markers from ESTs of hard clam Meretrix meretrix. Molecular Biology Reports38: 769–77520376702|
|28.||Li Y,, Liu W,, Xin G,, Gu Z,, Shao T,, et al. (Year: 2000) Introduction and F1 rearing of the American bay scallop Argopecten irradians irradians Lamarck. Heibei Fishery2: 29–32 (Chinese edition with English abstract).|
|29.||Xu F,, Zheng H,, Zhang H,, Liu X,, Zhang G, (Year: 2008) Comparison on the larval traits between strain “Zhongkehong” and common cultured population of bay scallop Argopecten irradians under different temperature. Journal of Fisheries of China32: 876–883 (Chinese edition with English abstract).|
|30.||Zhang F,, He Y,, Qi L,, Sun L,, Liu B, (Year: 1994) Introduction and spat-rearing of F1 generation of Argopecten irradians concentricus Say. Oceanologica et Limnologia Sinica25: 372–377 (Chinese edition with English abstract).|
|31.||Zheng H,, Zhang G,, Liu X,, Zhang F,, Guo X, (Year: 2004) Different responses to selection in two stocks of the bay scallop, Argopecten irradians irradians Lamarck (1819). Journal of Experimental Marine Biology and Ecology313: 213–223|
|32.||Elek JA,, Adamkewicz SL, (Year: 1990) Polymorphism for shell color in the Atlantic bay scallop Argopecten irradians irradians (Lamarck) (Mollusca: Bivalvia) on Martha's Vineyard island. American Malacological Bulletin7: 117–126|
|33.||Adamkewicz L,, Castagna M, (Year: 1988) Genetics of shell color and pattern in the bay scallop Argopecten irradians. Journal of Heredity79: 14–17|
|34.||Qin Y,, Liu X,, Zhang H,, Zhang G,, Guo X, (Year: 2007) Identification and mapping of amplified fragment length polymorphism markers linked to shell color in bay scallop, Argopecten irradians irradians (Lamarck, 1819). Marine Biotechnology9: 66–7317160637|
|35.||Roberts S,, Romano C,, Gerlanch G, (Year: 2005) Characterization of EST derived SSRs from the bay scallop, Argopecten irradians. Molecular Ecology Notes5: 567–568|
|36.||Zhan A,, Bao Z,, Wang X,, Hu J, (Year: 2005) Microsatellite markers derived from bay scallop (Argopecten irradians) expressed sequence tags. Fisheries Science71: 1341–1346|
|37.||Zhan A,, Hu J,, Wang X,, Lu W,, Hui M,, et al. (Year: 2006) A panel of polymorphic EST-derived microsatellite loci for the bay scallop (Argopecten irradians). Journal of Molluscan Studies72: 435–438|
|38.||Wang L (2005) The linkage mapping in Zhikong scallop Chlamys farreri and bay scallop Argopecten irradians Doctoral Dissertation of Institution of Oceanology, Chinese Academy of Sciences, Qingdao, China.|
|39.||Li H,, Liu X,, Hu J,, Bao Z,, Zhang G, (Year: 2007) A set of polymorphic microsatellite loci for the bay scallop, Argopecten irradians. Molecular Ecology Notes7: 422–424|
|40.||Li H (2009) Development of microsatellite DNA markers and construction of genetic linkage maps for bay scallop (Argopecten irradians). Doctoral dissertation. Institute of Oceanology, Chinese Academy of Sciences, Qingdao, China (Chinese edition with English abstract).|
|41.||Li H,, Liu X,, Du X,, Song R,, Zhang G,, et al. (Year: 2009) Development of microsatellite markers in bay scallop and their inheritance patterns in an F1 hybrid family. Marine Sciences33: 4–8 (Chinese edition with English abstract).|
|42.||Li H,, Liu X,, Zhang G, (Year: 2011) Development and linkage analysis of 104 new microsatellite markers for bay scallop (Argopecten irradians). Marine Biotechnology14: 1–921590517|
|43.||Woram RA,, McGowan C,, Stout JA,, Gharbi K,, Ferguson MM,, et al. (Year: 2004) A genetic linkage map for Arctic char (Salvelinus alpinus): evidence for higher recombination rates and segregation distortion in hybrid versus pure strain mapping parents. Genome47: 304–31515060583|
|44.||Rice RW, (Year: 1989) Analyzing tables of statistical test. Evolution43: 223–225|
|45.||Van Ooijen JW, Voorrips RE (2001) JoinMap 3.0, software for the calculation of genetic linkage maps. Plant Research International, Wageningen, the Netherlands.|
|46.||Stam P, (Year: 1993) Construction of integrated genetic linkage maps by means of a new computer package: Join Map. The Plant Journal3: 739–744|
|47.||Voorrips RE, (Year: 2002) MapChart: software for the graphical presentation of linkage maps and QTL. Journal of Heredity93: 77–7812011185|
|48.||Fishman L,, Kelly AJ,, Morgan E,, Willis JH, (Year: 2001) A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics159: 1701–171611779808|
|49.||Chakravarti A,, Lasher LK,, Reefer JE, (Year: 1991) A maximum likelihood method for estimating genome length using genetic linkage data. Genetics128: 175–1822060775|
|50.||Van Ooijen JW (2004) MapQTL version 5.0: software for the mapping of quantitative trait loci in experimental populations. Plant Research Int, Wageningen, the Netherlands.|
|51.||Shaw PW,, Pierce GJ,, Boyle PR, (Year: 1999) Subtle population structuring within a highly vagile marine invertebrate, the veined squid Loligo forbesi, demonstrated with microsatellite DNA markers. Molecular Ecology8: 407–417|
|52.||McGoldrick DJ,, Hedgecock D,, English L,, Baoprasertkul P,, Ward RD, (Year: 2000) The transmission of microsatellite alleles in Australian and North American stocks of the Pacific oyster (Crassostrea gigas): selection and null alleles. Journal of Shellfish Research19: 779–788|
|53.||Reece KS,, Ribeiro WL,, Gaffney PM,, Carnegie RB,, Allen SK Jr, (Year: 2004) Microsatellite marker development and analysis in the eastern oyster (Crassostrea virginica): confirmation of null alleles and non-Mendelian segregation ratios. Journal of Heredity95: 346–35215247315|
|54.||Wang Y,, Guo X, (Year: 2004) Chromosomal rearrangement in pectinidae revealed by rRNA loci and implications for bivalve evolution. Biology Bulletin207: 247–256|
|55.||Wang L,, Zhang H,, Song L,, Guo X, (Year: 2007) Loss of allele diversity in introduced populations of the hermaphroditic bay scallop Argopecten irradians. Aquaculture271: 252–259|
|56.||Grattapaglia D,, Sederoff R, (Year: 1994) Genetic linkage map of Eucalyptus grandis and Eucalyptus urophylla using a pseudotestcross: mapping strategy and RAPD markers. Genetics137: 1121–11377982566|
|57.||Launey S,, Hedgecock D, (Year: 2001) High genetic load in the Pacific oyster Crassostrea gigas. Genetics159: 255–26511560902|
|58.||Hubert S,, Hedgecock D, (Year: 2004) Linkage maps of microsatellite DNA markers for the Pacific oyster Crassostrea gigas. Genetics168: 351–36215454548|
|59.||Lyttle TW, (Year: 1991) Segregation distorters. Annual Review Genetics25: 511–557|
|60.||Pardo-Manuel de Villena F,, Sapienza C, (Year: 2001) Nonrandom segregation during meiosis: the unfairness of females. Mammalian Genome12: 331–33911331939|
|61.||Rieseberg LH,, Baird SJE,, Gardner KA, (Year: 2000) Hybridization, introgression, and linkage evolution. Plant Molecular Biology42: 205–22410688138|
|62.||Mitton JB, (Year: 1977) Shell color and pattern variation in Mytilus edulis and its adaptive significance. Chesapeake Science18: 387–390|
|63.||Gary FN, (Year: 1980) Genetics of shell color in Mytilus edulis L. and the association of growth rate with shell color. Journal of Experimental Marine Biology and Ecology47: 89–94|
|64.||Cole TJ, (Year: 1975) Inheritance of juvenile shell colour of the oyster drill, Urosalpix cinerea. Nature257: 794–7951186863|
|65.||Innes DJ,, Leslie EH, (Year: 1977) Inheritance of a shell-color polymorphism in the mussel. Journal of Heredity68: 203–204|
|66.||Liu X,, Wu F,, Zhao H,, Zhang G,, Guo X, (Year: 2009) A novel shell color variant of the Pacific abalone Haliotis discus hannai Ino subject to genetic control and dietary influence. Journal of Shellfish Research28: 419–424|
|67.||Moen T,, Hoyheim B,, Munck H,, Gomez-Raya L, (Year: 2004) A linkage map of Atlantic salmon (Salmo salar) reveals an uncommonly large difference in recombination rate between the sexes. Animal Genetics35: 81–9215025566|
|68.||Baranski M,, Loughnan S,, Austin CM,, Robinson N, (Year: 2006) A microsatellite linkage map of the blacklip abalone, Haliotis rubra. Animal Genetics37: 563–57017121601|
|69.||Lindahl KF, (Year: 1991) His and hers recombinational hotspots. Trends in Genetics7: 273–2761763424|
|70.||Liu B (1998) Statistical genomics: linkage, mapping, and QTL analysis. CRC Press, Boca Raton, FL.|
|71.||Beavis WD,, Grant D, (Year: 1991) A linkage map based on information from four F2 populations of maize (Zea mays L.). TAG Theoretical and Applied Genetics82: 636–644|
|72.||Mackay TFC, (Year: 1996) The nature of quantitative genetic variation revisited: Lessons from Drosophila bristles. BioEssays18: 113–1218851044|
|73.||Grattapaglia D,, Bertolucci FL,, Sederoff RR, (Year: 1995) Genetic mapping of QTLs controlling vegetative propagation in Eucalyptus grandis and E. urophylla using a pseudo-testcross strategy and RAPD markers. TAG Theoretical and Applied Genetics90: 933–947|
|74.||Melchinger AE,, Utz HF,, Schon CC, (Year: 1998) Quantitative trait locus (QTL) mapping using different testers and independent population samples in maize reveals low power of QTL detection and large bias in estimates of QTL effects. Genetics149: 383–4039584111|
|75.||Massault C,, Bovenhuis H,, Haley C,, de Koning DJ, (Year: 2008) QTL mapping designs for aquaculture. Aquaculture285: 23–29|
|76.||Wang CM,, Lo LC,, Feng F,, Zhu ZY,, Yue GH, (Year: 2008) Identification and verification of QTL associated with growth traits in two genetic backgrounds of Barramundi (Lates calcarifer). Animal Genetics39: 34–3918076743|
|77.||Sanchez-Molano E,, Cerna A,, Toro M,, Bouza C,, Hermida M,, et al. Detection of growth-related QTL in turbot (Scophthalmus maximus). BMC Genomics12: 47321958071|
Table 1 Number of markers and genetic length for sex-averaged linkage map of bay scallop.
|LG||Length of LG (cM)||No. markers||Averaged spacing(cM)|
Table 2 Size-related QTL identified in two families (CC5 and CC10) of bay scallop.
|Trait||QTL name||LG||Position (cM)||Confidence interval (cM)#||Nearest marker||LOD||PVE(%)||Phenotype means||Significance|
nt101SL, shell length; SH, shell height; TW, total weight; PVE, percentage of phenotypic variance explained by the QTL; association between four parental genotypes (f1m1, f2m1, f1m2, f2m2) of markers nearest each QTL peak and mean phenotypic values of each trait was tested by Kruskal-Wallis analysis.
#Confidence interval at 95% significance;
nt105*P<0.1, **P<0.05, *** P<0.01, ****P<0.005, ***** P<0.001.
Previous Document: Salivary immune responses to the 7-valent pneumococcal conjugate vaccine in the first 2 years of lif...
Next Document: Protein Tyrosine Phosphatase-Induced Hyperactivity Is a Conserved Strategy of a Subset of BaculoViru...