Crystals in hand

A genetic association of the ENPP1 gene with primary hand osteoarthritis was recently reported in this journal. ENPP1 encodes an enzyme that regulates soft tissue calcification. The study as it stands is far from complete because the actual causal variant(s) within ENPP1 has not been identified and no functional study on the activity of the enzyme in hand osteoarthritis was presented. Nevertheless, the study stimulates interest and will encourage others in the field to test ENPP1 as a possible osteoarthritis susceptibility gene in their cohorts. The genetic basis of osteoarthritis is slowly being uncovered, and this report constitutes another interesting find.     

Genetic variations in <Emphasis Type="Italic">KIFC1</Emphasis> and the risk of aspirin exacerbated respiratory disease in a Korean population: an association analysis

Modest effects of genes in various pathways are significant in the etiology of complex human diseases, including aspirin exacerbated respiratory disease (AERD). By functioning as a relevant component of respiratory processes, the human kinesin family member C1 (KIFC1) is hypothesized to play a role in AERD pathogenesis. A case–control analysis was carried out by comparing the genotype distribution of six KIFC1 single-nucleotide polymorphisms between 93 AERD cases and 96 aspirin-tolerant asthma controls in a Korean population. After controlling for confounds, logistic and regression models via various modes of genetic inheritance facilitated the association analysis. Initial results revealed significant association at 0.05 level of significance between several KIFC1 variations and AERD (P = 0.01–0.05, OR = 1.81–1.90) as well as fall rate of forced expiratory volume in the 1st second, an important diagnostic marker of airways constriction (P = 0.04–0.05). However, the signals were not deemed significant after multiple testing corrections (P ^corr > 0.05). Although the results do not support a major role of KIFC1 in AERD pathogenesis in a Korean asthma cohort, further replication and validation studies are required to clarify the current findings.     

A tagging SNP in <Emphasis Type="Italic">ALOX5AP</Emphasis> and risk of stroke: a haplotype-based analysis among eastern Chinese Han population

A genome-wide approach found significant association of two at-risk haplotypes (HapA, HapB) in the ALOX5AP gene with myocardial infarction and stroke. To date, it is still controversial whether ALOX5AP gene polymorphisms are risk factors for stroke. The aim of the present study is to investigate the association between the ALOX5AP gene polymorphism and the risk for stroke in Eastern Chinese Han population with a haplotype-based analysis. We conducted a comprehensive association study of 507 stroke patients and 510 healthy controls to assess the association between the ALOX5AP tagging single-nucleotide polymorphisms (tSNPs) and stroke risk. Genotyping was performed using the PCR–RFLP assay. In the single-locus analysis, we found that the rs9579646 AG genotype was associated with a marginally decreased risk for stroke (adjusted odds ratio, 0.65; 95% confidence interval, 0.45–0.96), compared with the AA genotype. Haplotype-based association analysis of block 2 involving rs10507391 and rs12429692 revealed that the decreased risk of stroke was significantly associated with haplotype AA (OR, 0.66; 95% CI, 0.46–0.95). These results suggested that the genetic variants in ALOX5AP might modulate the risk of stroke in Eastern Chinese Han population. The frequencies of single-marker alleles and haplotypes showed remarkable differences from those in other populations.     

Genetic variation in <Emphasis Type="Italic">IGFBP2</Emphasis> and <Emphasis Type="Italic">IGFBP5</Emphasis> is associated with breast cancer in populations of African descent

The insulin-like growth factor (IGF) signaling pathway is thought to play a major role in the etiology of breast cancer. Although incidence rates of breast cancer overall are lower in African Americans than in Caucasians, African-American women have a higher incidence under age 40 years, are diagnosed with more advanced disease, and have poorer prognosis. We investigated the association of breast cancer and genetic variants in genes in the IGF signaling pathway in a population-based case–control study of African-American women. We found significant associations at a locus encompassing parts of the IGFBP2 and IGFBP5 genes on chromosome 2q35, which we then replicated in a case–control study of Nigerian women. Based on those initial findings, we genotyped a total of 34 single nucleotide polymorphisms (SNPs) across the region in both study populations. Statistically significant associations with breast cancer were observed across approximately 50 kb of DNA sequence encompassing three exons in the 3′ end of IGFBP2 and three exons in the 3′ end of IGFBP5. SNPs were associated with breast cancer risk with P values as low as P = 0.0038 and P = 0.01 in African-Americans and Nigerians, respectively. This study is the first to report associations between genetic variants in IGFBP2 and IGFBP5 and breast cancer risk.     

Association of <Emphasis Type="Italic">ADAM33</Emphasis> gene polymorphisms with adult-onset asthma and its severity in an Indian adult population

ADAM33, a member of the ADAM (a disintegrin and metalloprotease) gene family, is an asthma susceptibility gene originally identified by positional cloning. In the present study, we investigated the possible association of five single-nucleotide polymorphisms (SNPs) in the ADAM33 (rs511898, rs528557, rs44707, rs597980 and rs2787094) with adult-onset asthma in an Indian population. The study included 175 patients with mild intermittent (n = 44), mild persistent (n = 108) or moderate persistent (n = 23) subgroups of asthma, and 253 nonasthmatic control individuals. SNPs were genotyped with the help of restriction fragment length polymorphism polymerase chain reaction (RFLP-PCR) method, and data were analysed using chi-square test and logistic regression model. Bonferroni’s correction for multiple comparisons was applied for each hypothesis. Genotypes and allele frequencies of SNPs rs511898 and rs528557 were significantly associated with adult-onset asthma (P = 0.010-<0.001). A significant association of the homozygous mutant genotype and mutant alleles of SNPs rs2787094, rs44707 and rs597980 with the asthma was also observed (P = 0.020-<0.001). A positive association between asthma and haplotypes AGCCT, GGCCT, AGACT, GCAGT, GGACT, ACCCC and AGACC were also found (P = 0.036-<0.001, OR = 2.07–8.49). Haplotypes AGCGT, GCAGC, ACAGC, ACAGT, GGAGC and GGCGT appear to protect against asthma (P = 0.013-<0.0001, OR = 0.34–0.10). Our data suggest that ADAM33 gene polymorphisms serve as genetic risk factors for asthma in Indian adult population.     

Copy number variations of CBF genes at the Fr‐A2 locus are essential components of winter hardiness in wheat

Winter hardiness is important for the adaptation of wheat to the harsh winter conditions in temperate regions and is thus also an important breeding goal. Here, we employed a panel of 407 European winter wheat cultivars to dissect the genetic architecture of winter hardiness. We show that copy number variation (CNV) of CBF (C‐repeat Binding Factor) genes at the Fr‐A2 locus is the essential component for winter survival, with CBF‐A14 CNV being the most likely causal polymorphism, accounting for 24.3% of the genotypic variance. Genome‐wide association mapping identified several markers in the Fr‐A2 chromosomal region, which even after accounting for the effects of CBF‐A14 copy number explained approximately 15% of the genotypic variance. This suggests that additional, as yet undiscovered, polymorphisms are present at the Fr‐A2 locus. Furthermore, CNV of Vrn‐A1 explained an additional 3.0% of the genotypic variance. The allele frequencies of all loci associated with winter hardiness were found to show geographic patterns consistent with their role in adaptation. Collectively, our results from the candidate gene analysis, association mapping and genome‐wide prediction show that winter hardiness in wheat is a quantitative trait, but with a major contribution of the Fr‐A2 locus.     

Association between <Emphasis Type="Italic">GSTM1</Emphasis>, <Emphasis Type="Italic">GSTT1</Emphasis>, and <Emphasis Type="Italic">GSTP1</Emphasis> polymorphisms and lung cancer risk in a Turkish population

Several studies focused on investigating genetic polymorphisms in order to estimate genetic contribution to lung cancer often showed conflicting results. In this study, we investigated the role of GSTM1, GSTT1, GSTP1 exon 5 and exon 6 polymorphisms on developing lung cancer and histological subtypes in 213 lung cancer patients and 231 controls. GSTM1 null, GSTT1 null, and GSTP1 exon 5 variant genotypes did not show a significant risk for developing lung cancer overall. Significant association was noted between GSTP1 exon 6 variant genotypes and overall lung cancer risk (OR 2.17, 95% CI 1.25–3.78; P = 0.006). These results show that GSTP1 exon 6 polymorphism might be an important factor in determining lung cancer susceptibility in a Turkish population.     

<Emphasis Type="Italic">MUTYH</Emphasis> Tyr165Cys, <Emphasis Type="Italic">OGG1</Emphasis> Ser326Cys and <Emphasis Type="Italic">XPD</Emphasis> Lys751Gln polymorphisms and head neck cancer susceptibility: a case control study

In the present study we investigated the association between three polymorphisms of the MUTYH (Tyr165Cys, rs34612342), the OGG1 (Ser326Cys, rs1052133) and the XPD (Lys751Gln, rs13181) genes with head and neck cancer risk. Genotypes were determined in DNA from peripheral blood lymphocytes of 265 patients with head and neck squamous cell carcinoma (HNSCC) as well as 280 cancer-free controls by PCR-restriction fragment lenght polymorphisms. We found an association between HNSCC and the Ser326Cys (OR 1.69; 95% CI 1.19–2.45) as well as Cys326Cys (OR 4.56; 95% CI 2.07–10.05) variants of the OGG1 gene. The gene–gene interaction between MUTYH and OGG1 as well as OGG1 and XPD polymorphic variants may contribute to higher prevalence of HNSCC. We also found an association between Ser326Cys and Cys326Cys variants of OGG1 gene and smoking status in HNSCC patients (OR 1.97; 95% CI 1.25–3.11), (OR 3.54; 95% CI 1.39–9.04), respectively. Moreover, we also observed a protective association between Tyr165Cys variant of the MUTYH gene and non-smoking status in HNSCC (OR 0.34; 95% CI 0.17–0.66). We also found a link between gene–gene interaction (MUTYH and OGG1 or OGG1 and XPD) and smoking (ORs 2.17–4.20 and 2.18–5.23) or non-smoking status (ORs 0.11 and 7.61) in HNSCC patients, respectively. In conclusion our data showed that the Ser326Cys polymorphism of the OGG1 gene may modify the risk of HNSCC associated with smoking. Finally we suggested that this polymorphism might be used as predictive factor for head and neck cancer in Polish population.     

Genetic dissection of sex determinism, inflorescence morphology and downy mildew resistance in grapevine

A genetic linkage map of grapevine was constructed using a pseudo-testcross strategy based upon 138 individuals derived from a cross of Vitis vinifera Cabernet Sauvignon × Vitis riparia Gloire de Montpellier. A total of 212 DNA markers including 199 single sequence repeats (SSRs), 11 single strand conformation polymorphisms (SSCPs) and two morphological markers were mapped onto 19 linkage groups (LG) which covered 1,249 cM with an average of 6.7 cM between markers. The position of SSR loci in the maps presented here is consistent with the genome sequence. Quantitative traits loci (QTLs) for several traits of inflorescence and flower morphology, and downy mildew resistance were investigated. Two novel QTLs for downy mildew resistance were mapped on linkage groups 9 and 12, they explain 26.0–34.4 and 28.9–31.5% of total variance, respectively. QTLs for inflorescence morphology with a large effect (14–70% of total variance explained) were detected close to the Sex locus on LG 2. The gene of the enzyme 1-aminocyclopropane-1-carboxylic acid synthase, involved in melon male organ development and located in the confidence interval of all QTLs detected on the LG 2, could be considered as a putative candidate gene for the control of sexual traits in grapevine. Co-localisations were found between four QTLs, detected on linkage groups 1, 14, 17 and 18, and the position of the floral organ development genes GIBBERELLIN INSENSITIVE1, FRUITFULL, LEAFY and AGAMOUS. Our results demonstrate that the sex determinism locus also determines both flower and inflorescence morphological traits.     

Isolation and characterization of nine microsatellite loci for the takin (<Emphasis Type="Italic">Budorcas taxicolor</Emphasis>)

We report on the isolation and characterization of nine microsatellite markers in the takin (Budorcas taxicolor) from genomic DNA-enriched libraries. Twenty-eight microsatellites were screened from the libraries, and nine of the screened microsatellites were polymorphic. The number of observed alleles for each locus in 28 individuals ranged from two to seven, and the expected and observed heterozygosity was 0.105–0.758 and 0.071–0.821, respectively. Four loci (TK01, TK02, TK04 and TK08) of nine deviated from Hardy-Weinberg expectation and no significant linkage association was found among all these loci. These microsatellite markers provide useful tool for population genetic studies of the takin.     

Linkage disequilibrium based association mapping of fiber quality traits in <Emphasis Type="Italic">G. hirsutum</Emphasis> L. variety germplasm

Cotton is the world’s leading cash crop, but it lags behind other major crops for marker-assisted breeding due to limited polymorphisms and a genetic bottleneck through historic domestication. This underlies a need for characterization, tagging, and utilization of existing natural polymorphisms in cotton germplasm collections. Here we report genetic diversity, population characteristics, the extent of linkage disequilibrium (LD), and association mapping of fiber quality traits using 202 microsatellite marker primer pairs in 335 G. hirsutum germplasm grown in two diverse environments, Uzbekistan and Mexico. At the significance threshold (r ² ≥ 0.1), a genome-wide average of LD extended up to genetic distance of 25 cM in assayed cotton variety accessions. Genome wide LD at r ² ≥ 0.2 was reduced to ~5–6 cM, providing evidence of the potential for association mapping of agronomically important traits in cotton. Results suggest linkage, selection, inbreeding, population stratification, and genetic drift as the potential LD-generating factors in cotton. In two environments, an average of ~20 SSR markers was associated with each main fiber quality traits using a unified mixed liner model (MLM) incorporating population structure and kinship. These MLM-derived significant associations were confirmed in general linear model and structured association test, accounting for population structure and permutation-based multiple testing. Several common markers, showing the significant associations in both Uzbekistan and Mexican environments, were determined. Between 7 and 43% of the MLM-derived significant associations were supported by a minimum Bayes factor at ‘moderate to strong’ and ‘strong to very strong’ evidence levels, suggesting their usefulness for marker-assisted breeding programs and overall effectiveness of association mapping using cotton germplasm resources. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genetic analysis of candidate genes modifying the age-at-onset in Huntington’s disease

The expansion of a polymorphic CAG repeat in the HD gene encoding huntingtin has been identified as the major cause of Huntington’s disease (HD) and determines 42–73% of the variance in the age-at-onset of the disease. Polymorphisms in huntingtin interacting or associated genes are thought to modify the course of the disease. To identify genetic modifiers influencing the age at disease onset, we searched for polymorphic markers in the GRIK2, TBP, BDNF, HIP1 and ZDHHC17 genes and analysed seven of them by association studies in 980 independent European HD patients. Screening for unknown sequence variations we found besides several silent variations three polymorphisms in the ZDHHC17 gene. These and polymorphisms in the GRIK2, TBP and BDNF genes were analysed with respect to their association with the HD age-at-onset. Although some of the factors have been defined as genetic modifier factors in previous studies, none of the genes encoding GRIK2, TBP, BDNF and ZDHHC17 could be identified as a genetic modifier for HD.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Characterization of phytoene synthase 1 gene (<Emphasis Type="Italic">Psy1</Emphasis>) located on common wheat chromosome 7A and development of a functional marker

Phytoene synthase (Psy), a critical enzyme in the carotenoid biosynthetic pathway, demonstrated high association with the yellow pigment (YP) content in wheat grain. Characterization of Psy genes and the development of functional markers for them are of importance for marker-assisted selection in wheat breeding. In this study, the full-length genomic DNA sequence of a Psy gene (Psy-A1) located on chromosome 7A, was characterized by in silico cloning and experimental validation. The cloned Psy-A1 comprises six exons and five introns, 4,175 bp in total, and an ORF of 1,284 bp. A co-dominant marker, YP7A, was developed based on polymorphisms of two haplotypes of Psy-A1, yielding 194 and 231-bp fragments in cultivars with high and low YP content, respectively. The marker YP7A was mapped on chromosome 7AL using an RIL population from cross PH82-2/Neixing 188, and a set of Chinese Spring nullisomic–tetrasomic lines and ditelosomic line 7AS. Psy-A1, co-segregating with the STS marker YP7A, was linked to SSR marker Xwmc809 on chromosome 7AL with a genetic distance of 5.8 cM, and explained 20–28% of the phenotypic variance for YP content across three environments. A total of 217 Chinese wheat cultivars and advanced lines were used to validate the association between the polymorphic band pattern and grain YP content. The results showed that the functional marker YP7A was closely related to grain YP content and, therefore, could be used in wheat breeding programs targeting of YP content for various wheat-based products. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Development of tetranucleotide microsatellite loci for the finless porpoise (<Emphasis Type="Italic">Neophocaena phocaenoides</Emphasis>)

The finless porpoise, Neophocaena phocaenoides, is endemic to the coastal waters of the Indo-Pacific, ranging from the Persian Gulf to Japan. Nine tetranucleotide microsatellite loci were isolated from the finless porpoise (Neophocaena phocaenoides). Polymorphism of each locus was assessed in 39 unrelated individuals from the Yellow Sea and the South China Sea of the Chinese waters. The number of alleles per locus varied from 2 to 11. The ranges of observed and expected heterozygosity were 0.154–0.795 and 0.146–0.839, respectively. Cross-species amplification of these loci was tested in other cetacean species. These microsatellite markers described here will be suitable for population genetic studies of finless porpoises and other cetacean species.     

The association of polymorphisms in DNA base excision repair genes <Emphasis Type="Italic">XRCC1</Emphasis>, <Emphasis Type="Italic">OGG1</Emphasis> and <Emphasis Type="Italic">MUTYH</Emphasis> with the risk of childhood acute lymphoblastic leukemia

The aim of this study was to evaluate the association of polymorphisms in genes encoding three key proteins of DNA base excision repair (BER): the OGG1 Ser326Cys, the MUTYH Tyr165Cys and the XRCC1 Arg399Gln with the risk of childhood acute lymphoblastic leukemia (ALL). Our study included 97 children patients with ALL (mean age 5.4 ± 2.5) and 131 healthy children (mean age 6.2 ± 2.8) used as controls. Genetic polymorphisms in BER pathway genes were examined using PCR and restriction fragment length polymorphism (RFLP). We have demonstrated that the OGG1 Cys/Cys genotype increases the risk of ALL (OR 5.36) whereas the Ser/Ser genotype variant strongly reduces the risk of this cancer among Polish children (OR 0.45). Although we did not observe the differences in single nucleotide polymorphisms (SNPs) in MUTYH and XRCC1 genes between control group and children with ALL, we have shown that the combined genotypes of examined genes can modulate the risk of childhood ALL in Polish population. We found that the combined genotype Arg/Gln–Cys/Cys of XRCC1/OGG1 (OR 3.83) as well as the Cys/Cys–Tyr/Tyr of OGG1/MUTYH (OR 6.75) increases the risk of ALL. In contrast, the combined genotype Arg/Arg–Ser/Ser of XRCC1/OGG1 (OR 0.40) as well as the Ser/Ser–Tyr/Tyr of OGG1/MUTYH (OR 0.43) played a protective role against this malignant disease. In conclusion, we suggest that polymorphisms of BER genes may be used as an important predictive factor for acute lymphoblastic leukemia in children.     

Estimation of population structure in coastal Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco var. menziesii] using allozyme and microsatellite markers

Characterizing population structure using neutral markers is an important first step in association genetic studies in order to avoid false associations between phenotypes and genotypes that may arise from non-selective demographic factors. Population structure was studied in a wide sample of ∼1,300 coastal Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco var. menziesii] trees from Washington and Oregon. This sample is being used for association mapping between cold hardiness and phenology phenotypes and single-nucleotide polymorphisms in adaptive-trait candidate genes. All trees were genotyped for 25 allozyme and six simple sequence repeat (SSR) markers using individual megagametophytes. Population structure analysis was done separately for allozyme and SSR markers, as well as for both datasets combined. The parameter of genetic differentiation (θ or F _ST) was standardized to take into account high within-population variation in the SSR loci and to allow comparison with allozyme loci. Genetic distance between populations was positively and significantly correlated with geographic distance, and weak but significant clinal variation was found for a few alleles. Although the STRUCTURE simulation analysis inferred the same number of populations as used in this study and as based on previous analysis of quantitative adaptive trait variation, clustering among populations was not significant. In general, results indicated weak differentiation among populations for both allozyme and SSR loci (θ _s = 0.006–0.059). The lack of pronounced population subdivision in the studied area should facilitate association mapping in this experimental population, but we recommend taking the STRUCTURE analysis and population assignments for individual trees (Q-matrix) into account in association mapping.     

Program death 1 (<Emphasis Type="Italic">PD1</Emphasis>) haplotyping in patients with breast carcinoma

Located on chromosome 2q37.3, the programmed death 1 (PD1) gene encodes for PD-1 (also known as CD279), a negative co-stimulator in the immune system. PD-1 renders potent inhibitory effects on T and B lymphocytes as well as monocyte responses. Expression of PD-1 ligands by tumor cells has been reported to contribute in immune system evasion. We aimed, in current study, to investigate the association of two single nucleotide polymorphisms in PD1 gene, +7146 G to A (PD-1.3) and +7785 C to T (PD-1.5 or +872), with susceptibility and/or progression of breast carcinoma. Four hundred forty-three women with breast cancer and 328 age-sex match healthy donors were recruited in present study. Genotyping was performed using Nested polymerase chain reaction-restriction fragment length polymorphisms. Arlequin software package was used to check for the Hardy–Weinberg equilibration and to determine the haplotypes. Results revealed no significant differences in the frequencies of genotypes and alleles at PD-1.3 (P = 0.252 and 0.279 for genotypes and alleles, respectively) and PD-1.5 positions (P = 0.522 and 0.278 for genotypes and alleles, respectively). Four haplotypes were observed among populations with no differences in the frequency between patients and controls. Our results also revealed no association between PD1 genotypes and tumor stage, tumor size, tumor grade, lymph node involvement, vascular invasion, distant metastasis, and Nottingham prognostic index. Present data do not confirm association of PD-1.3 (+7146) G/A and PD-1.5 (+7785 or +872) C/T genetic markers with susceptibility of Iranians to breast cancer.     

A genetic linkage map for an apricot (<Emphasis Type="Italic">Prunus armeniaca</Emphasis> L.) BC<Subscript>1</Subscript> population mapping <Emphasis Type="Italic">plum pox virus</Emphasis> resistance

Plum pox virus (sharka; PPV) can cause severe crop loss in economically important Prunus species such as peach, plum, apricot, and cherry. Of these species, certain apricot cultivars (‘Stark Early Orange’, ‘Goldrich’, ‘Harlayne’) display significant levels of resistance to the disease and are the genetic substrate for studies of several xlaboratories working cooperatively to genetically characterize and mark the resistance locus or loci for marker-assisted breeding. The goals of the work presented in this communication are the characterization of the genetics of PPV resistance in ‘Stark Early Orange’ and the development of co-dominant molecular markers for marker-assisted selection (MAS) in PPV resistance breeding. We present the first genetic linkage map for an apricot backcross population of ‘Stark Early Orange’ and the susceptible cultivar ‘Vestar’ that segregates for resistance to PPV. This map is comprised of 357 loci (330 amplified fragment length polymorphisms (AFLPs), 26 simple sequence repeats (SSRs), and 1 morphological marker for PPV resistance) assigned to eight linkage groups. Twenty-two of the mapped SSRs are shared in common with genetic reference map for Prunus (T × E; Joobeur et al. 1998) and anchor our apricot map to the general Prunus map. A PPV resistance locus was mapped in linkage group 1 and four AFLP markers segregating with the PPV resistance trait, identified through bulk segregant analysis, facilitated the development of SSRs in this region. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Lalli, D.A. and Salava, J. contributed equally to this work.     

A first set of gene-derived polymorphic microsatellite markers in half-smooth tongue sole (<Emphasis Type="Italic">Cynoglossus semilaevis</Emphasis>)

Half-smooth tongue sole (Cynoglossus semilaevis) is a commercially important marine fish species in China. A set of type I microsatellite markers were identified through bioinformatic mining of the GenBank database. Thirteen of these markers showed polymorphisms through genotyping a sample of 30 individuals. A total of 47 alleles were detected, with an average of 3.6 alleles per locus. The number of alleles, observed and expected heterozygosity per locus ranged from two to five, from 0.14 to 0.93 and from 0.27 to 0.77, respectively. Three loci significantly deviated from Hardy-Weinberg equilibrium after Bonferroni correction (P < 0.0038) and no significant linkage disequilibrum between pairs of loci was found. The markers identified in this study will contribute to construction of genetic linkage maps and mapping of quantitative trait loci (QTL) of C. semilaevis.