Open access peer-reviewed chapter

Population Genetic Studies of Silkworm (Bombyx mori L.) and Phylogenetic Relationships

Written By

Teodora Atanasova Staykova and Panomir Ivanov Tzenov

Submitted: 23 June 2023 Reviewed: 26 June 2023 Published: 29 August 2023

DOI: 10.5772/intechopen.1002229

From the Edited Volume

Biodiversity and Ecology of Lepidoptera - Insights and Advances

Farzana Khan Perveen

Chapter metrics overview

76 Chapter Downloads

View Full Metrics

Abstract

The study aimed to evaluate the degree of genetic variability and phylogenetic relationships between 13 breeds of mulberry silkworm (Bombyx mori L.) from genetic resources of Bulgaria through isozyme polymorphism. PAGE was used. Among nine studied isoenzyme loci, by eight loci (Bes A, Bes B, Bes D, Bes E, Pgm A, Mdh A, Bph, and Alp A) we found intra-breed and inter-breed polymorphism. At the Hk locus, we found inter-breed polymorphism only. The number of alleles per polymorphic locus ranged from one to two. The degree of polymorphism ranged from 0% to 77.80%. Low levels of observed heterozygosity in comparison with the expected one have been calculated in all of breeds. The combined FIS value over all polymorphic loci was 0.3205, which reflects a substantial deficit of heterozygotes. The value of FST showed that 49.21% of the overall genetic diversity observed was among breeds. The dendrogram constructed manifested that the two breeds of Japanese origin (Daizo and Japanese 106) were genetically most distant from other breeds. The data for isoenzyme polymorphism and genetic structure of the tested breeds can be used for genetic improvement and to develop new hybrids for silk production.

Keywords

  • silkworm Bombyx mori L.
  • isoenzymes
  • polymorphism
  • population genetic parameters
  • phylogenetic relationships

1. Introduction

Until the end of the 80s of the last century, Bulgaria was one of the best producers of cocoons in Europe. For various reasons, mostly economic, this industry is now in decline. The favorable climate and the existing rich national traditions are prerequisites for its restoration and further development of sericulture in Bulgaria, because the demand for silk and its place in everyday life will continue to be highly valued due to its hygienic qualities and finesse. Bulgaria maintains a rich genetic resources of more than 250 breeds of different origins, which is also a prerequisite for the recovering of sericulture.

The selection process of the mulberry silkworm is related to the solution of some basic issues as selection of individuals with the highest productivity, reproductive ability, viability and resistance to diseases, as well as analyzing and evaluating the capabilities of the breed gene pool. The creation of new highly productive breeds requires evaluation of the promising features for selection, creation of lines with desired qualities and analysis of their gene pool, development of methods for creation of synthetic lines, and evaluation of their genotypic and phenotypic features.

A basic principle for improving breeds is the presence of genetic diversity. Genetic variability can be analyzed by different methods and markers. The method of electrophoresis provides an opportunity to analyze the genetic heterogeneity in populations by studying genetically determined protein polymorphism [1, 2, 3, 4, 5, 6]. The established variability in isozyme markers can be used to characterize the genetic heterogeneity and degree of polymorphism of breeds to study the intensity of gene flow and the origin of individual breeds [7]. Application of isoenzymes and other molecular markers helps to estimate genetic diversity much more accurately than that of morphological traits [8]. Isoenzyme analysis is useful for the study of intra- and inter-breed polymorphism of mulberry silkworm and determining the level of genetic variability and genetic relationships [4, 6, 9, 10, 11]. Isoenzymes like esterase, acid phosphatase (ACP), alkaline phosphatase (ALP), malate dehydrogenase (MDH), and phosphoglucomutase have been used by various researchers to study diversity in silkworm genotypes [8, 12].

Studies for detection of polymorphic enzymatic and nonenzymatic protein systems in the breeds of B. mori kept in Bulgaria were made [4, 6, 10, 11, 13, 14, 15]. Data on the level of polymorphism and heterozygosity in more than 50 breeds with different origin have been established. The level of polymorphism varied between 0 and 77.8%, and the degree of heterozygosity varied between 0 and 0.280 in different breeds. Phylogenetic relationships between breeds have been identified on the basis of the established genetic distance and similarity. All these data are important for the maintenance of biodiversity in the silkworm, as well as for the selection of this species. However, there is still no information about other breeds, as the genetic resources of the silkworm maintained in Bulgaria are rich and include a large number of breeds. All this motivates the present investigation which aim was to study the degree of genetic variability and phylogenetic relationships between 13 breeds of silkworms, introduced and created in Bulgaria, on the basis of isozyme polymorphism.

Advertisement

2. Isoenzyme polymorphism and population genetic characteristics of silkworm breeds from the genetic resources of Bulgaria and their phylogenetic relationships

2.1 Material and methods

The silkworm resources used in the present investigation include a total of 13 breeds with different geographical origin and phenotype characteristic. They were obtained from the Scientific Center of Sericulture in Vratsa at the Agricultural Academy in Bulgaria. Breed Vratza 16 was created in Bulgaria. All others have been introduced as follow: breeds AES-1 wh and AES-1 zb originated from Spain, Tg—from Italia, Japanese 106 and Daizo originated from Japan, Mir 5—from Egypt, Mziuri 1—from Georgia, Tahvon 106—from Notrh Korea, Ukrainian 19—from Ukraine, Sh 4—from China, Line 22—from Uzbekistan, and MNB—from Madagascar. Daizo is polyvoltine, while all other breeds are mono-bivoltine. Breeds Tg, Daizo, AES-1 wh, and AES-1 zb have color cocoons. All other breeds have white cocoons. All individuals were nourished at a standard regime of silkworm breeding.

Totally 493 larvae on the fifth day of the fifth instar were studied. Larvae were selected randomly from each breed and were submitted to electrophoretic analysis of hemolymph, silk glands, and midgut tissues.

The tissue extracts were prepared according to the procedure described earlier [1114, 16, 17]. The individual samples were studied by 7.5% polyacrilamide gel electrophoresis (PAGE) [18] for nonspecific esterases (EST, EC 3.1.1), malate dehydrogenase (MDH, EC 1.1.1.37), and acid phosphatase (ACP, EC 3.1.3.2)—from the hemolymph; hexokinase (HK, EC 2.7.1.1)—from the silk glands and alkaline phosphatase (ALP, 3.1.3.1)—from the midgut. The 6% PAGE was used to analyze phosphoglucomutase (PGM, EC 5.4.2.2) from the silk glands. The staining mixtures for the enzymatic activities tested were pointed previously [19].

The phenotypes of the discovered loci were recorded after the revelation of the isozyme activity regions. Allele frequencies, mean number of alleles per locus, proportion of polymorphic loci, observed (Ho) and expected (He) heterozygosity, deviation from the Hardy-Weinberg equilibrium, Nei’s genetic distance (D) [20], and Wright’s fixation index, FST [21] were calculated using BIOSYS-1 [22]. Dendrogram was constructed using Nei’s [20] genetic distance, by UPGMA [23] method using the PHYLIP [24] software package.

2.2 Results

Isoenzyme and allozyme polymorphism of nonspecific esterases and allozyme polymorphism of phosphoglucomutase, malate dehydrogenase, acid phosphatase, alkaline phosphatase, and hexokinase were detected by polyacrylamide gel electrophoresis (Figure 1).

Figure 1.

PAGE spectra of: a. Nonspecific esterases of hemolymph; b. Phosphoglucomutase of silk glands; c. Alkaline phosphatase of midgut; d. Acid phpsphatase of hemolymph; e. Hexokinase of silk glands; f. Malate dehydrogenase from hemolymph of Bombyx mori L.

The tested enzymes recorded a total of nine polymorphic loci with 26 alleles (Table 1). Breed specificity of gene pools with respect to allele content and allele frequencies was established.

Locus (alleles)Breeds
AES1 zbAES1 whTgMir5Japanese 106MNBTahvon 106Mziuri 1DaizoLine 22Sh 4Ukrainian 19Vratza 16
Bes A
A10.6220.6971.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.000
A00.3780.3030.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.000
Bes B
B10.0000.0000.0000.4240.0000.3710.0000.4620.0000.2500.1400.7500.891
B20.5140.6450.5000.2881.0000.4290.3460.5001.0000.2080.1350.2120.076
B30.4860.3550.5000.2880.0000.2000.6540.0380.0000.5420.7160.0380.033
Bes D
D10.5810.1580.5380.4700.0000.7710.8080.1410.0000.1810.1081.0001.000
D20.4190.8420.4620.5301.0000.2290.1920.1410.0000.6800.1490.0000.000
D30.0000.0000.0000.0000.0000.0000.0000.5000.0000.1390.5270.0000.000
D00.0000.0000.0000.0000.0000.0000.0000.2181.0000.0000.2160.0000.000
Bes E
E10.1350.0920.2000.0001.0001.0000.0770.0000.0000.1810.0000.0000.000
E20.0410.0920.0620.1820.0000.0000.1790.0001.0000.0000.0270.0000.000
E00.8240.8160.7380.8180.0000.0000.7441.0000.0000.8190.9731.0001.000
Pgm A
A10.0000.0000.0000.0150.0000.0000.0000.0131.0000.1670.0000.0000.000
A20.3380.4740.5500.6971.0001.0001.0000.5770.0000.3890.4730.7001.000
A30.6620.5260.4500.2880.0000.0000.0000.4170.0000.4440.5270.3000.000
Mdh A
A20.9461.0001.0001.0001.0001.0000.1431.0001.0001.0001.0000.9251.000
A30.0540.0000.0000.0000.0000.0000.8570.0000.0000.0000.0000.0750.000
Bph A
A0.1760.1970.0000.3641.0000.1000.1030.2051.0000.0000.1620.1880.250
B0.0270.0390.0000.6360.0000.5000.0000.2050.0000.3890.0000.0000.304
C0.0000.0000.0000.0000.0000.2710.0000.0000.0000.0000.0000.0000.000
D0.0000.0000.0000.0000.0000.1290.0000.0000.0000.0000.0000.0000.000
00.7970.7631.0000.0000.0000.0000.8970.5900.0000.6110.8380.8130.446
Alp A
A11.0000.5661.0001.0001.0001.0001.0001.0001.0001.0001.0001.0001.000
A00.0000.4340.0000.0000.0000.0000.0000.0000.0000.0000.0000.0000.000
Hk A
A10.0000.0000.0000.0001.0000.0000.0000.0000.0000.0000.0000.0000.000
A21.0001.0001.0001.0000.0001.0001.0001.0001.0001.0001.0001.0001.000

Table 1.

Allele frequencies in breeds tested.

In the gene pool of AES 1 zb and AES 1 wh breeds (Table 1), we found two alleles polymorphism with “null” allele of the blood esterase A locus—Bes A0 and A1. In all rest breeds, Bes A1 allele was fixed. Bes A2 allele, which was described earlier [6] in some Egyptian breeds, was not detected in the current sample of breeds. For the Bes B locus, three alleles were recorded (Bes B1, B2, and B3) in the gene pools of seven breeds—Mir 5, MNB, Mziuri 1, Line 22, Sh 4, Ukrainian 19, and Vratza 16. Polymorphism with two alleles (Bes B2 and B3) was established in AES 1 zb, AES 1 wh, Tg, and Tahvon 106. Bes B2 allele was fixed in the gene pool of Daizo and Japanese 106. Among the breeds with polymorphism on Bes B locus, the allele Bes B1 showed the highest frequency in Vratza 16, Bes B2—in AES 1 wh and Bes B3—in Sh 4. Polymorphism with four alleles (Bes D1, D2, D3, and D0) was found on the Bes D locus. Bes D1 allele was fixed in the gene pool of Ukrainian 19 and Vratza 16, Bes D2—in Japanese 106, Bes D0—in Daizo. All four alleles were presented in the gene pool of the breeds Mziuri 1 and Sh4. Bes D1, D2, and D3 alleles were presented in the gene pool of Line 22, Bes D1 and D2—in all the rest tested breeds. We obtained the highest frequency of the allele Bes D1 in Tahvon 106, of Bes D2—in AES 1wh, of Bes D3—in Sh 4 (except for breeds with fixed Bes D alleles). Mziury 1 and Sh 4 have similar frequencies of the allele Bes D0. Polymorphism with three alleles was found in Bes E locus (Bes E1, E2, and E0) in AES 1 zb, AES 1 wh, Tg, and Tahvon 106 breeds. Two of these alleles we obtained in Line 22 (Bes E1 and Eo), Mir 5, and Sh 4 (Bes E2 and Eo). Bes E1 allele was fixed in Japanese 106 and MNB gene pool, Bes E2—in Daizo and Bes E0—in Mziuri 1, Ukrainian 19, and Vratza 16. The “null” allele Bes E0 demonstrated the highest frequency in all polymorphic breeds tested (Table 1).

Among the studied breeds we found polymorphism at the phosphoglucomutase (Pgm) locus with three alleles—Pgm A1, A2, and A3 in Mir 5, Mziuri 1, and Line 22 breeds. Pgm A2 and A3 were presented in the gene pool of AES 1 zb, AES 1 wh, Tg, Sh 4, and Ukrainian 19 (Table 1). Monomorphism of phosphoglucomutase demonstrated Daizo (with fixed Pgm A1 allele), Japanese 106, MNB, Tahvon 106, and Vratza 16 (with fixed Pgm A2 allele). Pgm A1 allele had the highest frequency in Line 22, Pgm A2—in Ukrainian 19 and Pgm A3—in AES 1 zb.

We found two alleles at the malatedehydrogenase locus (Mdh) in breeds AES 1 zb, Tahvon 106, and Ukrainian 19 (Mdh A2 and A3). Among these three breeds Mdh A2 was the most common allele in AES 1 zb and Ukrainian 19, while Mdh A3 was the most common in Tahvon 106. Mdh A2 allele was fixed in the gene pools of the rest 10 breeds (Table 1). Mdh A1 allele, which was described earlier [6, 12] in other breeds, was not detected in the current sample of breeds.

Total of five alleles of the acid phosphatase locus (Bph) were found in tested breeds (Bph A, B, C, D, and the “null” allele Bph O) (Table 1). Four of them were presented in the gene pool of MNB breed (Bph A, B, C, and D), three—in AES 1 zb, AES 1 wh, Mziuri 1, and Vratza 16 (Bph A, B, and 0) and two—in Mir 5 (Bph A and B), Tahvon 106, Sh 4 and Ukrainian 19 (Bph A and O), and Line 22 (Bph B and O). Bph A allele was fixed in Japanese 106 and Daizo breeds, whereas Bph O allele was fixed in Tg. The “null” allele was the most common in eight of the tested breeds. Bph B was the most expressed allele in Mir 5 and MNB.

Two alleles at the alkaline phosphatase locus (Alph A1 and A0) with a higher allele frequency of Alph A1 were recorded in the breed AES 1wh, only (Table 1). Alph A1 allele was fixed in the gene pool of the rest 12 breeds.

We found inter-breed polymorphism with two alleles on the hexokinase (Hk A) locus. The Hk A1 allele was presented only in the gene pool of Japanese 106 breed, whereas Hk A2 allele was fixed in the gene pools of all the rest 12 breeds.

The number of alleles per locus calculated with BIOSYS-1 software package in the silkworm breeds analyzed using nine enzyme loci ranged from 1.0 (Japanese 106 and Daizo) to 2.0 (Mziuri 1, AES 1 zb, and AES 1wh) (Table 2). The degree of polymorphism (according to the criterion 0.99) was the highest for the two Spanish breeds AES 1 zb and AES 1 wh (77.80%), and the lowest—for the two Japanese breeds Japanese 106 and Daizo (0%). The observed heterozygosity (Ho) by polymorphic loci varied from 0.000 (for Japanese 106 and Daizo) to 0.222 (for AES 1 zb). The expected heterozygosity (He) was higher than the observed one (Ho) in all breeds with polymorphism. Significant differences (P ˂ 0.05) in genotype frequencies were seen at the most loci in breeds studied. Chi-square test (DF = 1÷3) showed that the deviations from the Hardy-Weinberg equilibrium were in result of excess of homozygotes and deficiency of heterozygotes.

BreedMean sample size per locusMean number of alleles per locusPercent polymorphic loci (P = 0.99)HoHe
Vratza 1646.0 ± 0.01.4 ± 0.322.20.085 ± 0.0740.095 ± 0.073
Ukrainian 1940.0 ± 0.01.6 ± 0.244.40.125 ± 0.0590.141 ± 0.062
Sh 437.0 ± 0.01.9 ± 0.444.40.144 ± 0.0570.215 ± 0.087
Line 2236.0 ± 0.01.9 ± 0.355.60.164 ± 0.0590.279 ± 0.094
Mziuri 139.0 ± 0.02.0 ± 0.444.40.157 ± 0.0670.255 ± 0.102
Tahvon 10639.0 ± 0.01.6 ± 0.244.40.097 ± 0.0420.153 ± 0.065
MNB35.0 ± 0.01.7 ± 0.433.30.156 ± 0.0850.185 ± 0.097
Japanese 10636.0 ± 0.01.0 ± 0.00.00.000 ± 0.0000.000 ± 0.000
Mir 533.0 ± 0.01.8 ± 0.355.60.189 ± 0.0740.264 ± 0.089
Tg40.0 ± 0.01.6 ± 0.244.40.119 ± 0.0630.214 ± 0.085
Daizo37.0 ± 0.01.0 ± 0.00.00.000 ± 0.0000.000 ± 0.000
AES 1 wh38.0 ± 0.02.0 ± 0.277.80.164 ± 0.0480.319 ± 0.066
AES 1 zb37.0 ± 0.02.0 ± 0.277.80.222 ± 0.0580.297 ± 0.070

Table 2.

Mean number of alleles per locus, proportion of polymorphic loci, observed (Ho) and expected heterozygosity (He).

The mean FST value over all loci, which is associated with the level of inter-breed differentiation, was 0.4921 and shows that 49.21% of the overall genetic diversity observed was among breeds (Table 3). 50.79% of genetic variations were within the breeds. The highest level of genetic diversity among breeds we found for the Hk A locus was 1.0000) and the lowest one for the Mdh A locus was 0.0571. For two loci (Bes D and Bes E), we established a level of inter-breed differentiation over than 50%. The heterozygosity in total populations FIT averaged to 0.6549 and shows that there was a deficit of heterozygotes in the tested breeds and correlates with the obtained lower level of heterozygosity observed compared to the expected one and with deviations from the Hardy-Weinberg equilibrium as well. The combined FIS value for all polymorphic loci was 0.3205, which also reflects a significant deficiency of heterozygotes.

LocusFISFSTFIT
Bes A0.67110.30840.7725
Bes B0.14550.37010.4618
Bes D0.47420.51510.7450
Bes E0.56430.66920.8559
Pgm A0.50230.41910.7109
Mdh A−0.07140.0571−0.0102
Bph A0.06180.45920.4925
Alp A0.30370.41460.5924
Hk A0.00001.00001.0000
Mean0.32050.49210.6549

Table 3.

F-statistics for all polymorphic loci studied.

The values of genetic distance [20] were calculated using the allele frequencies and ranged from 0.029 (between the breeds Vratza 16 and Ukrainian 19) to 0.730 (between Japanese 106 and Ukrainian 19).

Analysis of the results obtained from genetic distances and UPGMA dendrogram (Figure 2) revealed that all the 13 breeds were grouped into two major clusters. The first cluster included Japanese breeds Daizo and Japanese 106, while the second included the rest 11 breeds, which was distributed in two subgroups. The first of them included MNB breed. The second one included all others. This subgroup was distributed into two groups—four breeds (Mir 5, Tahvon 106, Ukrainian 19 and Vratza 16) were grouped to form one subgroup and six breeds (Sh4, Mziuri 1, AES 1 wh, Line 22, Tg and AES 1 zb) were grouped to form another subgroup.

Figure 2.

UPGMA dendrogram.

2.3 Discussion

The study of polymorphic enzymatic and nonenzymatic proteins in mulberry silkworm is important for the selection of this species. They could serve as a kind of “passport” of the parent breeds, on the basis of which it is possible to compile optimal variants of crossbreeding and predict the effect of heterosis. Studies on proteins and enzymes in the silkworm (B. mori L.) have long since begun. Breeds of different origins and geographical distribution, bred in different countries as Japan, China, Korea, India, Russia and others, have been studied. Different methods of analysis have been used. A wide range of results has been obtained, which reflects not only the application of different research methods but also the existing huge variety of breeds that are used in various breeding centers around the world. Significant and positive correlation has been established between some isoenzymes and some yield parameters [25, 26]. Successes in this area of research could significantly facilitate the production of desired phenotypes and introduce elements of rigorous planning in the selection process.

In this study, we indicated a total of 12 alleles of four esterase loci. Three of them were “null” alleles. In earlier studies, [6] have reported an allele Bes A2. We did not find this allele among the tested 13 breeds. “Null” alleles of the Bes A, D, and E loci were described in other breeds from Bulgarian germplasm resourses of silkworm [36, 12, 19, 27]. Polymorphism with five alleles was determined of the acid phosphatase. One of them was found as “null” type. Allozyme polymorphism with codominant alleles of this enzyme was reported earlier [1, 6, 10, 12, 28]. We found intra-breed polymorphism with three or two alleles of the phosphoglucomutase, malate dehydrogenase, and alkaline phosphatase, as well as inter-breed polymorphism with two alleles of the hexokinase. Some of the alleles of the polymorphic loci demonstrated breed specificity. For instance, Bes A0 allele was presented only in AES 1 zb and AES 1 wh. Alp A0 allele was presented in the gene pool of AES 1 wh, and Hk A1 was presented in Japanese 106 only.

The results based on population genetic analyses showed a certain degree of differentiation between the tested breeds. 50.79% of isoenzyme diversity is observed between breeds and 49.21% is maintained within breeds, which is in line with the diversity based of AFLP markers found in some Iranian breeds [29]. Larger proportion of genetic variations among B. mori strains (84.08%) and a relatively smaller within strains (15.92%) have been established by RAPID analysis [30]. The UPGMA dendrogram resolved 13 breeds into two main clusters. The two breeds of Japanese origin (Daizo and Japanese 106) form one cluster and all the others the other cluster. This grouping of silkworm breeds is probably related to their adaptation to the specific geographical conditions of the environment from which they originate, as well as to the founder effect in the introduced breeds. Japanese 106 and Ukrainian 19 were the most distant breeds, whereas Vratza 16 and Ukrainian 19 were the closest breeds.

Low level of heterozygosity among tested breeds has observed in this study. Heterozygote deficiencies probably results from low effective number of reproductive individuals, selection process, and inbreeding effect. Some authors [8] pointed that reduction in genetic diversity in silkworm might be mainly due to domestication, breeding systems, selection, genetic drift, and inbreeding. The effects of inbreeding can accumulate over many generations [31, 32]. Breeders use artificial selection for target characteristics which also leads to a reduction in genetic variations in the population. The import of breeds of different origins and their use in breeding programs would help maintain a higher level of genetic diversity, which is very important for selection of suitable parents required for successful development of improved breeds and hybrids of silkworm that have high adaptive potential [8]. In view of the differences found in the genetic structure of the studied breeds with different origins, the results obtained here would be useful for breeders in planning crossbreeding strategies to produce new hybrids and in the conservation programs of silkworm B. mori.

Advertisement

3. Conclusions

Our results complement the knowledge of the genetic variations among the silkworm breeds bred in Bulgaria. They confirmed that nonspecific esterases, acid phosphatase, and malate dehydrogenase from hemolymph, phosphoglucomutase and hexokinase from the silk glands, and alkaline phosphatase from the midgut are applicable to the study of genetic structure and phylogenetic relationships between breeds. The number of alleles, allelic, and genotypic frequencies at polymorphic loci show breed specificity. It is important to perform continuous evaluate the polymorphism degree of the breeds to avoid a marked increase of the homozygosity. This would result in the expression of deleterious genes that can cause high mortality or other adverse effects. The results obtained in the present study could help breeders in selecting parental pairs for crossbreeding and in determining the quality of parental forms in the early stages of development. The analysis of genetic structure on the basis of isoenzyme markers showed that the tested breeds of silkworm are genetically differentiated. Most of them have high degree of polymorphism and can be used for genetic improvement and to develop new hybrids for silk production.

Advertisement

Acknowledgments

This study was supported by Horizon ARACNE Project (Grant Agreement 101095188).

References

  1. 1. Eguchi M, Takahama Y, Ikeda M, Horii S. A novel variant of acid phosphatase isozyme from hemolymph of the silkworm, Bombyx mori. The Japanese Journal of Genetics. 1988;63(2):149-157
  2. 2. Marcato S, Trevisan M, Cappellozza L, Bisol PM. Variabilita di sistemi gene-enzima, in razze di Bombyx mori (Lepidoptera, Bombycidae) utilizzate per la produzione di seta. Estratto da Redia. 1990;73(2):595-608
  3. 3. He JL. Studies on the inheritance of a deletion type of the blood esterase isoenzymes in the silkworm, Bombyx mori L. Sericologia. 1995;35:17-24
  4. 4. Staykova T. Genetically-determined polymorphism of nonspecific esterases and phosphoglucomutase in eight introduced breeds of the silkworm, Bombyx mori, raised in Bulgaria. Journal of Insect Science. 2008;8(18):1-8
  5. 5. Moorthy SM, Chandrakanth N. Analysis of phenotypic diversity and protein polymorphism in some silkworm breeds of Bombyx mori L. Munis Entomology & Zoology. 2015;10(2):478-485
  6. 6. Staykova T, Ivanova E, Tzenov P, Vasileva Y, Pantaleeva D, Grekov D, et al. Genetic analysis of isoenzyme polymorphism in silkworm (Bombyx mori L.) (lepidoptera: Bombycidae) strains and phylogenetic relationships. Acta Zoologica Bulgarica. 2015;67(1):117-125
  7. 7. Kaidanov L. Genetics of Populations. High School, Russia: Moskow; 1996. p. 320
  8. 8. Bindroo B, Moorthy SM. Genetic divergence, implication of diversity, and conservation of silkworm, Bombyx mori. International Journal of Biodiversity. Hindawi publishing corporation. 2014:1-15
  9. 9. Egorova T, Nasirillaev Y. Polymorphic enzymes of silkworm and their using in breeding. In: State Committee of Science and Technics of Republic Uzbekistan. Tashkent, UZ: GFNTI Press; 1993. p. 120
  10. 10. Staykova T, Ivanova E, Tzenov P, Vasileva Y, Pantaleeva D, Petkov Z. Acid phosphatase as a marker for differentiation of silkworm (Bombyx mori) strains. Biotechnology & Biotechnological Equipment. 2010;24(2):379-384
  11. 11. Staykova T, Ivanova E, Grekov D, Avramova K. Genetic variability in silkworm (Bombyx mori L.) strains with different origin. Acta Zoologica Bulgarica. 2012;4:87-92
  12. 12. Staykova T, Tzenov P, Vasileva Y, Takova N, Ivanova E, Stoyanov I, et al. Isoenzyme polymorphism of silkworm Bombyx mori L. 1758 (lepidoptera: Bombicidae) breeds from germplasm resources of Bulgaria. Acta Zoologica Bulgarica. 2020;15:11-19
  13. 13. Staykova T. Electrophoretic analysis of haemolymph proteins during silkworm (Bombyx mori L.) ontogenesis. International Journal of Industrial Entomology. 2007;14(1):37-44
  14. 14. Staykova T, Grekov D, Panayotov M. Electrophoretic analysis of nonspecific esterases in silkworm (Bombyx mori L.) female genital organs and eggs. International Journal of Industrial Entomology. 2004;9(1):59-63
  15. 15. Staykova T, Ivanova E, Tzenov P, Vasileva Y, Pantaleeva D. Differentiation of silkworm (Bombyx mori L.) strains by isoenzime markers. Genetics and Breeding. 2009;38(1):47-55
  16. 16. Stoykova T, Popov P, Dimitrov B. Electrophoretic analysis of nonspecific haemolymph esterases during silkworm (Bombyx mori L.) ontogenesis. Sericologia. 2003;43(2):153-162
  17. 17. Staykova T. Inter- and intra-population genetic variability of introduced silkworm (Bombyx mori L.) strains raised in Bulgaria. Journal of Bioscience and Biotechnology. 2013;2(1):73-77
  18. 18. Davis BJ. Disc electrophoresis II. Method and application to human serum proteins. Annals of the New York Academy of Sciences. 1964;121(2):404-427
  19. 19. Staykova T, Tzenov P, Vasileva Y, Aziz S, Ivanova E, Stoyanov I, et al. Genetic characterization of silkworm (Bombyx mori L.) strains (lepidoptera: Bombicidae) with different geographical origin on the basis of isozyme markers. Acta Scientiarum – Biological Sciences. 2020;42(e4797):1-9
  20. 20. Nei M. Genetic distance between populations. The American Naturalist. 1972;106(949):283-292
  21. 21. Wriht S. The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution. 1965;19(3):395-420
  22. 22. Swofford D, Selander RB. BIOSYS-1: A Computer Program for the Analysis of Allelic Variation in Genetics. Rel. 1.0. Urbana, IL: Department of Genetics and Development University of Illinois at Urbana-Champaign; 1981
  23. 23. Sneath PHA, Sokal RR. Numerical Taxonomy: The Principles and Practice of Numerical Classification. San Francisco, CA: W.H. Freeman & Co.; 1973. p. 573
  24. 24. Felsenstein J. PHYLIP (Phylogeny Inference Package). Version 3.5c. Seattle, WA: University of Washington; 1993
  25. 25. Egorova T, Naletova E, Nasirillaev Y. Polymorphic system of silkworm haemolymph esterases as a criterion to make programs for parental specimens crossing. Biochemistry of Insects. Moskow. 1985:54-62
  26. 26. Maqbool A, Da H, Ahmad M, Zaffar G, Malik G, Mir SA. Haemolymph esterase and peroxidase spectra in bivoltine silkworm (Bombyx mori L.) genotypes and their utility. Journal of Cell and Tissue Research. 2015;15(3):5193-5199
  27. 27. Shabalina A. Esterase genetic polymorphism in haemo-lymph of larvae Bombyx mori L. Comptes rendus de l'Academie bulgare des Sciences. 1990;43(9):105-108
  28. 28. Eguchi M, Yoshitake N, Kai H. Types and inheritance of blood esterase in the silkworm, Bombyx mori L. Japan Journal of Genetics. 1965;40(1):15-19
  29. 29. Mirhoseini S, Dalirsefat S, Pourkheirandish M. Genetic characterization of Iranian native Bombyx mori strains using amplified fragment length polymorphism markers. Journal of Economic Entomology. 2007;100(3):939-945
  30. 30. Liu Y, Qin L, Li Y, Wang H, Xia RX, Qi YH, et al. Comparative genetic diversity and genetic structure of three Chinese silkworm species Bombyx mori L. (lepidoptera: Bombycidae), Antheraea pernyi Guérin-Meneville and Samia cynthia ricini Donovan (lepidoptera: Saturniidae). Neotropical Entomology. 2010;39(6):967-976
  31. 31. Lande R. Risk of population extinction from fixation of new deleterious mutations. Evolution. 1994;48(5):1460-1469
  32. 32. Whitlock MC. Selection, load and inbreeding depression in a large metapopulation. Genetics. 2002;160(3):1191-1202

Written By

Teodora Atanasova Staykova and Panomir Ivanov Tzenov

Submitted: 23 June 2023 Reviewed: 26 June 2023 Published: 29 August 2023