Translator Disclaimer
1 October 2010 Amplified Fragment Length Polymorphism Mapping of Quantitative Trait Loci for Economically Important Traits in the Silkworm, Bombyx mori
Author Affiliations +

Cocoon related characteristics are economically important traits in the silkworm, Bombyx mori L. (Lepidoptera: Bombycidae). In this study a genetic linkage map was developed that identified QTL controlling the cocoon weight, cocoon shell weight, and cocoon shell percentage using 161 amplified fragment length polymorphism (AFLP) markers. Twenty PstI/TaqI primer combinations were employed to genotype 78 F2 progenies derived from a cross between P107 Japanese inbred line and Khorasan Lemon Iranian native strain. Among polymorphic markers, 159 AFLP markers were assigned to 24 linkage groups at the LOD threshold of 2.5 that varied in length from 4 to 299 cM. The total length of the linkage map was 2747 cM, giving an average marker resolution of 19.31 cM. A total of 21 AFLP markers were identified that were distributed over the ten linkage groups linked to the three studied traits using the composite interval mapping method. The explained variation rate by QTL controlling cocoon weight, cocoon shell weight, and cocoon shell percentage ranged from 0.02% to 64.85%, 0.2% to 49.11%, and 0.04% to 84.20%, respectively. These QTL controlled by different actions as well as under dominance, additive, partial dominance, dominance, and over dominance.


The silkworm, Bombyx mori L. (Lepidoptera: Bombycidae), domesticated for silk production for about 5000 years, is an agriculturally important insect and comprises a large number of geographical races and inbred lines that show substantial variation in their qualitative and quantitative traits (Mirhoseini et al. 2007). With the establishment of stable transformation (Yamao et al. 1999; Tamura et al. 2000), B. mori has shown the potential to produce pharmaceutically important proteins in high yield (Tomita et al. 2003), opening up new applications for sericulture in medical, agricultural, and industrial fields (Yamamoto et al. 2006). Currently, it is the major economic resource for nearly 30 million families in countries such as China, India, Vietnam, and Thailand (Miao et al. 2005). In order to make sericulture economically viable, genes affecting growth rate, yield, fiber quality, and virus resistance can be tagged with molecular markers for rapid construction of genetically improved strains. Taking the exclusive investigational advantages of this organism into account, as well as its economic importance, an International Consortium on Lepidopteran Genomics was formed a few years ago to support international cooperation to sequence the genome of B. mori and to initiate comparative genomics of other economically important Lepidoptera (Nagaraju and Goldsmith 2002).

Most traits in nature and of importance to agriculture are quantitatively inherited and therefore are difficult to study due to the complex nature of their inheritance. However, recent advances of genomic technologies have led to revolutionary means for unraveling the secrets of genetic variation in quantitative traits. Genomic technologies allow the molecular characterization of polymorphic markers throughout the entire genome that are then used to identify and map the genes or quantitative trait loci (QTL) underlying a quantitative trait based on linkage analysis (Wu et al. 2007).

A complete linkage map is necessary to efficiently carry out molecular-based analyses such as molecular marker-assisted selection, quantitative trait loci (QTL) mapping of agronomically important traits, prediction of heterosis, and comprehensive investigations of genomic evolution between lineages (Tan et al. 2001). Presently, genome studies in B. mori have generated genetic linkage maps based on morphological markers (Doira et al. 1992) and molecular markers including RFLP (Goldsmith 1991; Shi et al. 1995; Nguu et al. 2005), RAPD (Promboon et al. 1995; Yasukochi 1998; Li et al. 2000), SADF and RAPD (He et al. 2001), AFLP (Tan et al. 2001; Lu et al. 2004; Sima et al. 2006), microsatellites (Miao et al. 2005), and SNP (Yamamoto et al. 2006).

Projects have been initiated to find molecular markers that are tightly linked to traits relevant for sericulture, with the related goals of developing tools for marker assisted selection and positional cloning. RAPD or cDNA markers have been associated with the four known densonucleovirus nonsusceptibility loci, nsd-1 (Abe et al. 1998), nsd-2 (Abe et al. 2000), Nid-1 (Kadono-Okuda et al. 2003), and nsd-Z (Li et al. 2001). Two large contigs on chromosome 17 that encompass cDNAs closely linked to Nid-1 and nsd-2 have been isolated and sequenced (Kadono-Okuda et al. 2003) and are being examined for candidate genes in susceptible and nonsusceptible strains (Goldsmith et al. 2005). A similar strategy was used to screen for RAPD markers linked to resistance to NPV, a potentially devastating pathogen (Yao et al. 2003), and fluoride resistance (Chen et al. 2003). Progress has also been made in assigning RAPDs (Chatterjee and Pradeep 2003), inter-simple sequence repeat markers (Chatterjee and Mohandas 2003), and AFLP (Lu et al. 2004; Li et al. 2006; Sima et al. 2006) to QTL for characters such as larval growth rate and pupal and cocoon weight. A suite of additional fingerprinting tools has been developed for these applications (Nagaraju and Goldsmith 2002). In addition, a collection of about 8500 expressed sequence tags (ESTs) is now available in GenBank and provides an additional source of important anchors in the ongoing Bombyx genome study (Nguu et al. 2005).

Although for genome mapping, the ideal genetic marker is codominant, multiallelic, and hypervariable (i.e., segregates in almost every family), some dominant markers are also very useful and powerful in particular situations (Wu et al. 2007). The amplified fragment length polymerphism (AFLP) technique (Zabeau and Vos 1992; Vos et al. 1995) has demonstrated to be a convenient and reliable tool to generate highly polymorphic molecular markers that greatly facilitate building linkage maps (Qi et al. 1997; Waugh et al. 1997). AFLP markers do allow one to construct linkage maps with wide genome coverage without engaging in extensive sequencing or marker development programs. AFLP markers are also faster than individual codominant marker types because a single polymerase chain reaction (PCR) can derive multiple loci simultaneously (Erickson et al. 2004). Because of these features, AFLP has been widely employed for genetic mapping in various organisms.

In this study, significant molecular markers and a large segregating population size were employed to detect QTL linked to economically important traits relevant to the B. mori cocoon and to better identify the genome regions of these QTL. Since high-resolution QTL mapping is critical for positional cloning and gene isolation (Zhong et al. 2006), a high resolution AFLP-based genetic linkage map and the results of QTL mapping for economically important cocoon traits are reported.

Materials and Methods

Insect materials and crosses

One F2 segregating family from mating between a Japanese inbred line (P107) as female parent and an Iranian native strain (Khorasan Lemon) as male parent were used in the study. These inbred line and strain exhibit high phenotype diversity for economically important characters such as whole cocoon weight, cocoon shell weight, and cocoon shell percentage, suggesting that considerable polymorphism exists at the DNA level (Dalirsefat and Mirhoseini, 2007). Indeed, the highest and the least quantities of these traits corresponded to P107 and Khorasan Lemon, respectively. These inbred line and strain have undergone a high degree of inbreeding and are relatively homozygous. A total number of 78 progenies, including 39 males and 39 females from F2 population, were used to construct the genetic linkage map and QTL detection. The parents and F1 progenies were used to establish the segregation pattern of the molecular markers. The crossing experiments were done in the Iran Silkworm Research Center (ISRC) located in Rasht, Iran.

AFLP analysis

Genomic DNAs were isolated individually from all the parents, F1, and F2 populations, in the moth stage following the phenol/chloroform method (Suzuki et al. 1972) and as modified by Nagaraja and Nagaraju (1995). DNAs were quantified using a known standard (DNA lambda, Roche, on agarose gels.

All individuals were subjected to genotyping with AFLP markers according to Vos et al. (1995) with some modifications. Briefly, genomic DNA was double digested with PstI and TaqI restriction enzymes, which produce polymorphic DNA fragments in B. mori (Tan et al. 2001; Mirhoseini et al. 2007). The DNA fragments were ligated with PstI and TaqI adaptors, generating template DNA for PCR amplification. Two primers were designed on the basis of adaptor sequences and restriction site sequences to use in PCR amplification. Selective nucleotide sequences were added to the 3′ end of each primer. PCR amplification was conducted in two steps: a pre-amplification and a selective amplification. For the selective amplification, a total of 81 primer combinations obtained from two sets of PstI and TaqI selective primers (Table 1) were screened. Among them, 20 primer pairs that produced fragments with clear dominance inheritance patterns and reproducibility were used for the linkage analysis. Polymorphism screening of AFLP products was conducted on a 6% Polyacrylamide gel using SequiGen 38×30 cm gel apparatus (BioRad Laboratories Inc., Bands were detected by the silver staining procedure (Promega,, Technical manual No.023), and gel images were scanned and saved as jpeg files for scoring and further analysis.

Linkage analysis and map construction Using genotype information of 81 AFLP primer combinations, 20 primer combinations that produced clearly readable and polymorphic fragments among parents were employed to analyze linkage mapping. Twenty polymorphic primer combinations generating 161 polymorphic AFLP fragments with a clear dominance inheritance pattern were employed to construct the linkage map and detect QTL; that is, the suitable fragments must show complete dominance expression in one parent and complete recessive expression in the other, and all F1 individuals must be heterozygous. The AFLP fragments were scored based on 0 and 1 and then converted to A, B, C, and D letters according to the Map manager QTX (Manly et al. 2001) instruction manual. The data were analyzed using the Kosambi map function (Kosambi 1944) of Map manager QTX (Manly et al. 2001) to develop a linkage map for the population. By genotyping 78 progenies from the F2 population using 161 polymorphic bands, a genotypic data matrix in a dimension of 78 × 161 was constructed and used for linkage mapping. First recombination rates among markers were evaluated, and then recombination rates converted to the map distance based on centiMorgan using the Kosambi map function (Kosambi 1944). Computer software QTL cartographer version 2.5 (Wang et al. 2007) was used to determine the QTL positions, the expected additive and dominance effects, and the phenotypic variance explained by individual QTL. The LOD threshold value for declaring the presence of a QTL was determined by a permutation test (n = 1000) (Churchill and Doerge 1994). Genome-wide threshold levels were used to declare significant QTL based at the 5% significance level. Average levels of dominance (h) were estimated using the ratio dominance/additive effects (Stuber et al. 1987).

All molecular experiments were established in the genomics laboratory of the Agricultural Biotechnology Research Institute of North Region (Rasht) under the supervision of the Agricultural Biotechnology Research Institute of Iran (ABRII).


Linkage map construction

Among the 81 AFLP primer combinations screened, approximately one-third of the primer combinations (n = 28) produced polymorphic fragments between the P107 inbred line and the Khorasan Lemon native strain. Twenty pairs of AFLP primer combinations were selected for segregation analysis on the F2 population based on reproducibility and the degree of polymorphism. Only the polymorphic fragments that segregated in a dominant manner and could be scored unambiguously were used for linkage map construction. An example of AFLP gel electrophoresis and polymorphism screening related to the Ptat-Ttac primer combination is shown in Figure 1.

Table 1.

Adapters and primers used in AFLP analysis.


Twenty PstI/TaqI primer combinations produced 845 clearly detected bands, of which 161 qualified polymorphic fragments showing good agreement of 3:1 segregation (for a dominant marker, the segregation ratio is 3:1 in the F2 population) were analyzed for linkage mapping. The frequency of polymorphic AFLP markers derived from the clearly detected bands in the P107 × Khorasan Lemon cross in the silkworm was 19.35%. This frequency was close to that obtained in the Dazao × C100 cross of the silkworm (25.7%) (Lu et al. 2004) but it was dramatically lower than in the no. 782 × od100 cross (60.7%) of the silkworm (Tan et al. 2001).

91 fragments of 161 polymorphic fragments (56.52%) were detected in the male parent (Khorasan Lemon strain), and 70 fragments (43.48%) were observed in the female parent (P107 inbred line). On average, each primer combination generated 8.05 polymorphic fragments that could be used for linkage mapping. The number of polymorphic bands produced using the 20 primer combinations ranged from 3 bands (7.32%) corresponding to P33-T32 to 14 bands (23.33%) for P22-T31 (Table 2).

The linkage map generated from the P107 × Khorasan Lemon cross contained 159 AFLP markers (two markers were unlinked) that were assigned to 24 linkage groups at the LOD threshold of 2.5 (Figure 2). Average distance between markers was 19.31 cM. The total recombination distance over 24 linkage groups was 2747 cM, which was longer than previous estimates in B. mori i.e. 1800 cM for the dense RAPD map (Yasukochi 1998), 1868.10 cM and 2677.50 M for the AFLP maps in two F2 subgroups (Sima et al. 2006), and 1305 cM for SNPs based linkage map (Yamamoto et al. 2006). However, it was shorter than 6512 cM (Tan et al. 2001) and 3676.7 cM (Li et al. 2006) for the AFLP maps and 3431.9 cM (Miao et al. 2005) for the SSR markers reported in backcrossed populations of B. mori. Miao et al. (2005) suggested that although many conditions influence map length, including differences in mating strategy and strains used, the distribution of markers is a possible causative aspect, and increased marker density should converge on a more realistic map length value. As Tan and Ma (1998) demonstrated theoretically, with additional markers typed, the map length may increase when marker density is not saturated or may decrease when marker density is in a saturation state (Tan et al. 2001). For example, Causse et al. (1994) constructed a rice map with 762 markers covering 4026.3 cM, whereas Harushima et al. (1998) obtained a 2275-marker genetic map of rice covering 1521.6 cM. This may explain why the length of our AFLP map is more than that of the B. mori linkage map studies mentioned above except for the maps of Tan et al. (2001), Miao et al. 2005, and Li et al. (2006).

Figure 1.

An example of AFLP gel electrophoresis and polymorphism screening corresponding to the Ptat- Ttac primer combination. GM-107 and GF-KL represent female and male parents, respectively. M stands for standard molecular size marker. Polymorphic bands are shown by arrows. High quality figures are available online.


Considering that the estimated genome size of B. mori is 530 Mbp (Gage 1974), the average physical distance per recombination distance is about 193 kb/cM. It seems that the AFLP markers did not exhibit significant clustering near centromeres or the distal region of chromosomes, suggesting that they provide good coverage of the genome (!”#$% et al. 2006, Figure 2).

Table 2.

The observed polymorphisms of twenty TaqI and PstI primer combinations used in the parents and F2 population.


Phenotypic values

The average phenotypic values of cocoon weight, cocoon shell weight, and cocoon shell percentage traits corresponding to each parent and F2 offspring are shown in Table 3. An extremely high significant difference (p < 0.01) for these traits was revealed as a result of comparing the mean phenotypic values between parents using the t-test. The mean cocoon weight in line P107 as female parent and Khorasan Lemon as male parent was 1.479 and 1.404 g, respectively. The mean cocoon shell weight in line P107 was 0.324 g, approximately 0.113 g more than that in Khorasan Lemon (0.211 g). In addition, the cocoon shell percentage in the line P107 was estimated to be 22.23%, which was seven percent more than that in Khorasan Lemon (15.23%) as the male parent. Except for cocoon weight value in the F2 population, which was higher than both parent values, both the mean cocoon shell weight and cocoon shell percentage traits in the F2 generation were closely equal to the mean parent values (Table 3).

Figure 2.

AFLP Linkage map of the silkworm (Bombyx mori) based on 76 F2 males derived from the cross between P107 and Khorasan Lemon. The numbers on the left side of each linkage group are genetic distances in Kosambi centiMorgans. AFLP markers are designated by the PstI and TaqI primer names on the right side of each linkage group. High quality figures are available online.


The frequency distribution of phenotypic values related to each trait studied in the F2 offspring is illustrated in Figure 3. As it is shown, a number of F2 offspring demonstrated out of range parent values, especially in the line P107 for the three studied traits. The least and the highest phenotypic values of cocoon weight in the F2 population were 1.13 g and 1.83 g, respectively, and the highest value was about 0.35 g more than that in the line P107 (1.479 g). The least and the highest phenotypic values of cocoon shell weight in the F2 population were 0.22 g and 0.37 g, respectively; consequently, a number of the F2 offspring had almost 0.046 g more than that in the line P107 phenotypic value (0.324 g). In addition, the least and the highest phenotypic values of cocoon shell percentage in the F2 progenies were 15.72% and 22.6%, respectively; the highest value was nearly 0.4% higher than that in the line P107 (22.2%) (Figure 3).

QTL analysis

A total of 21 different loci, including QTL, controlling cocoon weight, cocoon shell weight, and cocoon shell percentage traits were detected in the linkage map using a composite interval mapping method at the LOD threshold of 2.5 (Table 4). The selected LOD score plots at the threshold of 3 for linkage groups with the identified QTL provided a basis for identifying the molecular markers most closely linked to the QTL (Figure 4).

In particular, 12 QTL controlling cocoon weight were identified on the LG1, LG5, LG6, LG8, LG16, LG17, and LG19. The additive effects of these QTL ranged from -0.1581 (cw8) to +0.0887 (cw19b), and their dominance effects ranged from -0.3852 (cw19a) to +0.3881 (cw1a). Two QTL for cocoon shell weight were identified and located on the LG16 and LG22 with additive effects ranging from -0.0164 (cshw16c) to +0.0459 (cshw22b) and dominance effects ranging from -0.0314 (cshw16c) to +0.037 (cshw22b). Finally, fourteen QTL were identified for cocoon shell percentage located on the LG8, LG9, LG16, LG19, and LG23. The additive effects of these QTL ranged from -0.0347 (cshp23b) to +0.0118 (cshp8), and their dominance effects ranged from -0.0361 (cshp8) to + 0.0352 (cshp19a) (Table 4).

Table 3.

Mean phenotypic values of traits studied in the parents and F2 population.


Figure 3.

Frequency distribution of phenotypic values related to the studied traits in the F2 segregating Bombyx mori population derived from mating between the P107 inbred line and the Khorasan Lemon native strain (the mean phenotypic values of the parents for the three traits are shown in Table 3). High quality figures are available online.


Table 4.

Genetic parameters as estimated by composite interval mapping of the QTL affecting the cocoon related traits in the silkworm.


The explained variation rate by QTL controlling cocoon weight, cocoon shell weight, and cocoon shell percentage ranged from 0.02% to 64.85%, 0.2% to 49.11%, and 0.04% to 84.20%, respectively (Table 4).


In the present study, an AFLP-based linkage map containing 159 AFLP markers in a total length of 2747 cM and an average marker resolution of 19.31 cM was developed for B. mori. Using this map, a total of 21 AFLP markers linked to cocoon weight, cocoon shell weight, and cocoon shell percentage were identified using a composite interval mapping method (Table 4). Recently, 11 QTL (Lu et al. 2004) and 40 QTL (Li et al. 2006) for whole cocoon weight, cocoon shell weight, ratio of shell weight and weight of pupae have been reported, and Javadi Taklimi (2006) accounted 5 QTL controlling ratios of shell weight in a backcrossed population (BC1) of B. mori. These differences may be due to sample size, number and type of primer combinations used, and crosses established. Though apart from employing different primer pairs, Lu et al. (2004) and Li et al. (2006) used 44 BC1 progenies, and Javadi Taklimi (2006) applied only seven polymorphic primer combinations. Practically, factors such as the number of molecular markers used, types of crosses, sample size of segregating population, number of genes controlling the traits, and existence of gene interaction may influence the statistical power of QTL mapping (Zhong et al. 2006).

Figure 4.

Composite interval mapping of the cocoon related traits in Bombyx mori. Significance thresholds are indicated by horizontal lines, with LOD = 3.0 (genome wide P < 0.001) as determined by 1000 permutations of the mapping data (Churchill and Doerge 1994). High quality figures are available online.


Figure 4.


Figure 4.


Figure 4.


Figure 4.


It seems that some QTL had pliotropic effects on the traits. One QTL controlling both cocoon weight and cocoon shell weight traits located on the LG16 between neighboring markers Ptc/Taca7-Ptat/Ttac3 was observed. In addition, three QTL were observed controlling both cocoon weight and cocoon shell percentage traits located on the LG8, LG16, and LG19 between neighboring markers Paac/Ttg8-Ptat/Ttg7, Paac/Tat8-Ptat/Ttac2 and Paac/Tagc4-Paac/Tat7.

A total of 19.35% of clearly readable and qualified AFLP bands were polymorphic between the P107 inbred line and the Khorasan Lemon native strain of B. mori. A higher level (61%) of polymorphic AFLP marker has been reported by Tan et al. (2001) in a single backcross (no. 782 and od100) family of B. mori. To explain this approach, they discussed several factors:

  1. Employing two distinct B. mori strains in the present study, P107 and Khorasan Lemon are two examples of distinct silkworm inbred lines and strains. The former is from the Japanese bivoltine system and the latter is from the Iranian native monovoltine system.

  2. Detecting high levels of polymorphisms by the AFLP technique (Huys et al. 1996; Latorra et al. 1996; Mackill et al. 1996; Wan et al. 1999)

  3. A large fraction of the silkworm genome consists of families of transposable elements such as Bm1, BMC1 (a member of the LINE1 family), mariner, mariner-like elements (Bmmar1), long terminal repeat transposons (LTRs), non-long terminal transposons (nonLTRs), and others (Ueda et al. 1986; Herrer and Wang 1991; Xiong and Eickbush 1993; Xiong et al. 1993; Robertson and Asplund 1996; Tomita et al. 1997; Shimizu et al. 2000; Wang et al. 2000).

Among the 20 pairs of AFLP primer combinations applied in this study, an average of 9.2 polymorphic AFLP markers per primer combination for linkage analysis and QTL mapping was recognized. This rate was considerably lower than two other AFLP linkage and QTL mapping studies on silkworm with 35.7 (Tan et al. 2001) and 36.4 (Lu et al. 2004) fragments per primer. This may be due to the degree of differences between parental lines and strains and the primer combinations used.

The present AFLP map consisted of 24 linkage groups, whereas the haploid genome of B. mori has 28 chromosomes. As Promboon et al. (1995), Young et al. (1998), and He (1998) reported, this may be due to nonequivalence between the number of linkage groups and the number of chromosomes. In the RFLP based linkage map by Goldsmith (1991), 15 linkage groups were reported. However, by using morphological (Doira et al. 1992), RAPD (Yasukochi et al. 1998), RFLP (Nguu et al. 2005), AFLP (Sima et al. 2006), and SNP (Yamamoto et al. 2006) markers, 28 linkage groups and using SSR (Miao et al. 2005) markers, 29 linkage groups have been recognized in B. mori. It has also been shown that the large number of chromosomes in the haploid B. mori genome (n = 28), typical of Lepidoptera, makes it difficult to construct maps without missing some chromosomes (Yasukochi 1998).

The whole cocoon weight, cocoon shell weight, and cocoon shell percentage are the major economic traits in B. mori that are controlled by a polygene (Li et al. 2006). In the present study, a single F2 population derived from a cross between P107 Japanese inbred line and Khorasan Lemon Iranian native strain was used as the mapping population. Among the 21 QTL for the traits studied, one QTL had dominance effect, 13 QTL had under dominance or recessive effects, and seven QTL had over dominance effects (Table 4). Li et al. (2006) in a backcrossed population (BC1) derived from a cross between C100 and Dazao detected 40 QTL for whole cocoon weight and related traits, of which 19 were additive effect QTL and 21 were reduced effects QTL.

In summary, 159 AFLP markers were employed to construct a linkage map for B. mori with an average marker resolution of 19.31 cM. We identified 21 QTL (n = 21) using the composite interval mapping method that affects whole cocoon weight and related traits. The effects of these QTL were under dominance, dominance, and over dominance. Since AFLP amplification is highly reproducible, the development of an AFLP linkage map and subsequently the identification of strain-specific markers for tracking allele frequency changes and quantitative trait loci (QTL) analysis for economically important traits provides an invaluable tool for improving B. mori breeds, strains, and hybrids in order to enhance the silk production.


This work was supported by funds from Agricultural Biotechnology Research Institute of Iran (ABRII) under the supervision of Agricultural Research and Education Organization (AREO) (an organization of the Ministry Jihad-e Agriculture) and Agricultural Biotechnology Research Institute of North Region (Rasht). We thank Mr. Moeineddin Mavvajpour (the director) and Dr. Alireza Seydavi (the adviser) from ISRC (Iran Silkworm Research Centre) for establishing silkworm crosses and phenotypic information and for their invaluable advice. The technical assistance of the staff of ISRC is gratefully acknowledged. We are also grateful to Mohammad Naserani for contribution of data and typing of phenotypic records in Excel, for editing and inputing data in the linkage map and QTL software for further analysis.



Agricultural Biotechnology Research Institute of Iran;


Amplified Fragment Length Polymorphism;


Agricultural Research and Education Organization;




Composite Interval Mapping;




Expressed Sequence Tags;


Interval Mapping;


Iran Silkworm Research Center;


Linkage Group;


Logarithm of Odds;


Likelihood Ratio Test;


Nuclear Polyhidrosis Virus;


Polymerase Chain Reaction;


Quantitative Trait Loci;


Random Amplified Polymorphic DNA;


Restriction Fragment Length Polymorphism;


Selective Amplification of DNA Fragments;


Single Nucleotide Polymorphism



H Abe , T Harada , M Kanehara , T Shimada , F Ohbayashi , T Oshiki . 1998. Genetic mapping of RAPD markers linked to the densonucleosis refractoriness gene, nsd-1 in the silkworm, Bombyx mori. Genes and Genetic Systems 73: 237–42. Google Scholar


H Abe , T Sugasaki , M Kanehara , T Shimada , ST Gomi , F Ohbayashi . 2000. Identification and genetic mapping of RAPD markers linked to the densonucleosis refractoriness gene, nsd-2, in the silkworm, Bombyx mori. Genes and Genetic Systems 75: 93–96. Google Scholar


MA Causse , TM Fulton , YG Cho , SN Ahn , J Chumwongse , K Wu , J Xiao , Z Yu , PC Ronald , SE Hanington , G Second , SR McCouch , SD Tanksley . 1994. Saturated molecular map of the rice genome based on an interspecific backcross population. Genetics 138: 1251–1274. Google Scholar


SN Chatterjee , AR Pradeep . 2003. Molecular markers (RAPD) associated with growth, yield and origin of the silkworm, Bombyx mori L. in India. Genetika 39: 1612–24. Google Scholar


SN Chatterjee , TP Mohandas . 2003. Identification of ISSR markers associated with productivity traits in silkworm, Bombyx mori L. Genome 46: 438–47. Google Scholar


K Chen , Q Yao , M Li , Y Wang . 2003. Screening of RAPD markers for fluoride resistance in Bombyx mori L. International Journal of Industrial Entomology 7: 11–14. Google Scholar


GA Churchill , RW Doerge . 1994. Empirical threshold values for quantitative trait mapping. Genetics 138: 963–971. Google Scholar


SB Dalirsefat , SZ Mirhoseini . 2007. Assessing genetic diversity in Iranian native silkworm (Bombyx mori L.) strains and Japanese commercial lines using AFLP markers. Iranian Journal of Biotechnology 5(1): 25–33. Google Scholar


H Doira , H Fujii , Y Kawaguchi , H Kihara , Y Banno . 1992. Important Genetic Resources. In: Genetic Stocks and Mutations of Bombyx mori. p 73. Isseido Press. Google Scholar


DL Erickson , BC Fenster , HK StenØien , D Price . 2004. Quantitative trait locus analyses and the study of evolutionary process. Molecular Ecology 13: 2505–2522. Google Scholar


LP Gage . 1974. The Bombyx mori genome: Analysis by DNA reassociation kinetics. Chromosoma 45: 27–42. Google Scholar


MR Goldsmith . 1991. The Bombyx mori genome-mapping project. Sericologia 31: 145–155. Google Scholar


MR Goldsmith , T Shimada , H Abe . 2005. The genetics and genomics of the silkworm, Bombyx mori. Annual Review of Entomology 50: 71–100. Google Scholar


Y Harushima , M Yano , A Shomura , M Sato , T Shimano , Y Kuboki , T Yamamoto , SY Lin , BA Antonio , A Parco , H Kajiya , N Huang , K Yamamoto , Y Nagamura , N Kurata , GS Khush , T Sasaki . 1998. High-density rice genetic linkage map with 2275 markers using a single F2 population. Genetics 148: 479–494. Google Scholar


NJ He . 1998. Construction of the molecular linkage map of the silkworm (Bombyx mori L.). Ph.D. Thesis, Southwest Agricultural University, Chongqing, China. Google Scholar


LJ He , C Lu , B Li , ZY Zhou , ZH Xiang . 2001. The construction of linkage map of Bombyx mori by combination of SADF and RAPD. Acta Entomologia Sinica 44(4): 476–482. (in Chinese) Google Scholar


RJ Herrer , J Wang . 1991. Evidence for a relationship between the Bombyx mori middle repetitive Bm1 sequence family and U1 snRNA. Genetica 84: 31–37. Google Scholar


G Huys , R Coopman , P Janssen , K Kersters . 1996. High-resolution genotypic analysis of the genus Aeromonas by AFLP fingerprinting. International Journal of Systematic Bacteriology 46: 572–580. Google Scholar


SA Javadi Taklimi . 2006. Identification QTL controlling cocoon shell percentage in silkworm. M.Sc. dissertation, Faculty of Agricultural Science, University of Guilan. Google Scholar


K Kadono-Okuda , DO Ogoyi , J Nohata , S Sasanuma , M Sasanuma , R Eguchi , W Hara , K Yamamoto , K Mita . 2003. Sequence analysis of BAC contigs covering a densovirus-nonsusceptibility gene, nsd-2 in Bombyx mori. Lepidoptera Abstracts. Presented at the Sixth International Workshop Molecular Biology and Genetics of Lepidoptera Kolymbari, Crete. Google Scholar


DD Kosambi . 1944. The estimation of map distance from recombination values. Annals of Eugenics 12: 172–175. Google Scholar


D Latorra , MS Schanfield . 1996. Analysis of human specificity in AFLP systems APOB, PAH, and D1S80. Forensic Science International 83: 15–25. Google Scholar


B Li , C Lu , ZY Zhou . 2000. The Construction of RAPD linkage map of Bombyx mori. Acta Genetica Sinica 27(2): 127–132. (in Chinese) Google Scholar


B Li , C Lu , AC Zhao , ZH Xiang . 2006. Multiple interval mapping for whole cocoon weight and related economically important traits QTL in silkworm (Bombyx mori). Agricultural Sciences in China 5(10): 798–804. Google Scholar


M Li , Q Yao , C Hou , L Cheng , K Chen . 2001. Studies on RAPD markers linked to the densonucleosis refractoriness gene, nsd-Z in the silkworm, Bombyx mori. Sericologia 41: 409–415. Google Scholar


C Lu , B Li , AC Zhao , ZH Xiang . 2004. QTL mapping of economically important traits in Silkworm (Bombyx mori), Science in China Series C Life Sciences 47(5): 477–484. Google Scholar


DJ Mackill , Z Zhang , ED Redona , PM Colowit . 1996. Level of polymorphism and genetic mapping of AFLP markers in rice. Genome 39: 969–977. Google Scholar


KF Manly , JR Cudmore , JM Meer . 2001. Map Manager QTX: Cross-platform software for genetic mapping. Mammalian Genome 12: 930–932. Google Scholar


XX Miao , SJ Xu , MH Lia , MW Lia , JH Huang , FY Dai , SW Marino , DR Mills , P Zeng , K Mita , SH Jia , Y Zhang , WB Liu , H Xiang , QH Guo , AY Xu , XY Kong , HX Lin , YZ Shi , G Lu , X Zhang , W Huang , Y Yasukochi , T Sugasaki , T Shimada , J Nagaraju , ZH Xiang , SY Wang , MR Goldsmith , C Lu , GP Zhao , YP Huang . 2005. Simple sequence repeat-based consensus linkage map of Bombyx mori. Proceedings of the National Academy of Sciences U.S.A. 102: 16303–16308. Google Scholar


SZ Mirhoseini , SB Dalirsefat , M Pour Khairandish . 2007. Genetic characterization of Iranian native Bombyx mori strains by using amplified fragment length polymorphism markers. Journal of Economic Entomology 100(3): 939–945. Google Scholar


GM Nagaraja , J Nagaraju . 1995. Genome fingerprinting in silkworm, Bombyx mori using random arbitrary primers. Electrophoresis 16: 1633–1638. Google Scholar


J Nagaraju , MR Goldsmith . 2002. Silkworm genomics-progress and prospects. Current Science 83: 415–425. Google Scholar


EK Nguu , K Kadono-Okada , K Mase , E Kosegawa , W Hara . 2005. Molecular linkage map for the silkworm, Bombyx mori based on restriction fragment length polymorphism of cDNA clones. Journal of Insect Biotechnology and Sericology 74: 5–13. Google Scholar


A Promboon , T Shimada , H Fujiwara , M Kobayashi . 1995. Linkage map of RAPDs in the silkworm, Bombyx mori. Genetics Research 66: 1–7. Google Scholar


X Qi , P Stam , P Lindhout . 1997. Use of the locus specific AFLP markers to construct a high density molecular map in barley. Theoretical and Applied Genetics 96: 376–384. Google Scholar


HM Robertson , ML Asplund . 1996. Bmmar1: A basal lineage of the mariner family of transposable elements in the silkworm moth, Bombyx mori. Insect Biochemistry and Molecular Biology 26: 945–954. Google Scholar


J Shi , DG Heckel , MR Goldsmith . 1995. A genetic linkage map for the domesticated silkworm, Bombyx mori based on restriction fragment length polymorphisms. Genetics Research 66: 109–126. Google Scholar


K Shimizu , M Kamba , H Sonobe , T Kanda , AG Klinakis , C Savakis , T Tamura . 2000. Extrachromosomal transposition of the transposable element minos occurs in embryos of the silkworm Bombyx mori. Insect Molecular Biology 9: 277–281. Google Scholar


YH Sima , B Li , DX Chen , DB Sun , AC Zhao , L Zhang , C Lu , SM He , ZH Xiang . 2006. Construction and analysis of an AFLP molecular linkage map of the silkworm (Bombyx mori). Chinese Journal of Agricultural Biotechnology 3(1): 25–31. Google Scholar


CW Stuber , MD Edwards , JF Wendel . 1987. Molecular markers facilitated investigations of quantitative trait loci in maize. II. Factors influencing yield and its component traits. Crop Science 27: 639–648. Google Scholar


Y Suzuki , LP Gage , DD Brown . 1972. The genes for fibroin in Bombyx mori. Journal of Molecular Biology 70: 637–649. Google Scholar


T Tamura , C Thibert , C Royer , T Kanda , A Eappen , M Kamba , N Kômoto , JL Thomas , B Mauchamp , G Chavancy , P Shirk , M Fraser , JC Prudhomme , P Couble . 2000. Germline transformation of the silkworm, Bombyx mori L. using a piggyBac transposon-derived vector. Nature Biotechnology 18: 81–84. Google Scholar


YD Tan , RL Ma . 1998. Estimates of lengths of genome and chromosomes of rice using molecular markers. Journal of Biomathematics 13:1022–1027. (in Chinese) Google Scholar


YD Tan , CL Wan , YF Zhu , C Lu , Z Xiang , HW Deng . 2001. An amplified fragment length polymorphism map of the silkworm, Genetics 157(3): 1277–1284. Google Scholar


M Tomita , H Munetsuna , T Sato , T Adachi , R Hino , M Hayashi , K Shimizu , N Nakamura , T Tamura , K Yoshizato . 2003. Transgenic silkworms produce recombinant human type III procollagen in cocoons. Nature Biotechnology 21: 52–56. Google Scholar


S Tomita , BH Sohn , T Tamura . 1997. Cloning and characterization of a mariner-like element in the silkworm, Bombyx mori. Genes and Genetic Systems 72: 219–228. Google Scholar


H Ueda , S Mizuno , K Shimura . 1986. Transposable genetic element found in the 5'-flanking region of the fibroin H-chain gene in a genomic clone from the silkworm Bombyx mori. Journal of Molecular Biology 190: 319–327. Google Scholar


P Vos , R Hogers , M Bleeker , M Reijans , TVD Lee , M Hornes , A Frijters , J Pot , J Peleman , M Kuiper , M Zabeau . 1995. AFLP: A new technique for DNA fingerprinting. Nucleic Acids Research 23: 4407–4414. Google Scholar


CL Wan , YF Zhu , YD Tan , C Lu . 1999. Application of AFLP markers to detection of genetic polymorphic loci in the silkworm (Bombyx mori, L.). Biotechnology 9: 4–9. (in Chinese) Google Scholar


S Wang , CJ Basten , ZB Zeng . 2007. Windows QTL cartographer Version 2.5. Statistical Genetics North Carolina State University, USA. Google Scholar


W Wang , L Swevers , K Iatrou . 2000. Mariner (Mos1) transposase and genomic integration of foreign gene sequences in Bombyx mori cells. Insect Molecular Biology 9: 145–155. Google Scholar


R Waugh , N Bonar , E Baird , B Thomas , A Graner , WTB Thos , W Powell . 1997. Homology of AFLP products in three mapping populations of barley. Molecular Genetics and Genomics 255: 311–321. Google Scholar


R Wu , CX Ma , G Casella . 2007. Statistical Genetics of Quantitative Traits-Linkage, Maps, and QTL. Springer. Google Scholar


Y Xiong , TH Eickbush . 1993. Dong, a non-long terminal repeat (non-LTR) retrotransposable element from Bombyx mori. Nucleic Acids Research 21: 1318. Google Scholar


Y Xiong , WD Burke , TH Eickbush . 1993. Pao, a highly divergent retrotransposable element from Bombyx mori containing long terminal repeats with tandem copies of the putative R region. Nucleic Acids Research 21: 2117–2123. Google Scholar


K Yamamoto , J Narukawa , K Kadono-Okuda , J Nohata , M Sasanuma , Y Suetsugu , Y Banno , H Fujii , MR Goldsmith , K Mita . 2006. Construction of a single nucleotide polymorphism linkage map for the silkworm, Bombyx mori based on bacterial artificial chromosome end sequences. Genetics 173: 151–161. Google Scholar


M Yamao , N Katayama , H Nakazawa , M Yamakawa , Y Hayashi , S Hara , K Kamei , H Mori . 1999. Gene targeting in the silkworm by use of a baculovirus. Genes and Development 13: 511–516. Google Scholar


Q Yao , MW Li , Y Wang , WB Wang , J Lu , Y Dong , KP Chen . 2003. Screening of molecular markers for NPV resistance in Bombyx mori L. (Lep., Bombycidae). Journal of Applied Entomology 127(3): 134–136. Google Scholar


Y Yasukochi . 1998. A dense genetic map of the silkworm, Bombyx mori covering all chromosomes based on 1018 molecular markers. Genetics 150: 1513–1525. Google Scholar


WP Young , PA Wheeler , VH Coryell , P Keim , GH Thorgaard . 1998. A detailed linkage map of rainbow trout produced using doubled haploids. Genetics 148: 839–850. Google Scholar


M Zabeau , P Vos . 1992. Selective restriction fragment amplification: A general method for DNA fingerprinting, European Patent Application no. 92402629.7, publication no. 0 534 858 A1. Google Scholar


D Zhong , DM Menge , EA Temu , H Chen , G Yan . 2006. Amplified fragment length polymorphism mapping of quantitative trait loci for malaria parasite susceptibility in the Yellow Fever mosquito Aedes aegypti. Genetics 173: 1337–1345. Google Scholar
This is an open access paper. We use the Creative Commons Attribution 3.0 license that permits unrestricted use, provided that the paper is properly attributed.
Seyed Z Mirhoseini, Babak Rabiei, Payam Potki, and Seyed B Dalirsefat "Amplified Fragment Length Polymorphism Mapping of Quantitative Trait Loci for Economically Important Traits in the Silkworm, Bombyx mori," Journal of Insect Science 10(153), 1-21, (1 October 2010).
Received: 7 April 2009; Accepted: 1 May 2010; Published: 1 October 2010

Back to Top