Silkworm genebanks assume paramount importance as the reservoirs of biodiversity and source of alleles that can be easily retrieved for genetic enhancement of popular breeds. More than 4000 Bombyx mori L (Lepidoptera: Bombycidae) strains are currently available and these strains are maintained through continuous sibling mating. This repeated sibling mating makes the populations of each strain more homozygous, but leads to loss of unique and valuable genes through the process of inbreeding depression. Hence, it is essential to maintain a minimal degree of heterozygosity within the population of each silkworm strain, especially in the traditional geographic strains, to avoid such loss. As a result, accurate estimation of genetic diversity is becoming more important in silkworm genetic resources conservation. Application of molecular markers help estimate genetic diversity much more accurately than that of morphological traits. Since a minimal amount of heterozygosity in each silkworm strain is essential for better conservation by avoiding inbreeding depression, this article overviews both theoretical and practical importance of heterozygosity together with impacts of inbreeding depression and the merits and demerits of neutral molecular markers for measurements of both heterozygosity and inbreeding depression in the silkworm Bombyx mori.

Improvement of high yielding, disease resistant silkworm strains became imminent to increase production of silk, which is a major revenue earner for sericulturists. Since environment interacts with phenotype, conventional breeding did not result in commendable yield improvement in synthetic strains of silkworm, Bombyx mori. Identification of DNA markers associated with different economically important biomass traits and its introgression could assist molecular breeding and expression of stabilized high yielding characters, but genetic basis of most quantitative traits in silkworm is poorly understood due to its polygenic control. Correlation analysis (R = 0.9) revealed significant interrelation among biomass traits viz., larval duration (TLD), larval weight (LWT), cocoon weight (CWT), shell weight (SWT), shell ratio (SR) and floss content. PCR using inter simple sequence repeat (ISSR) primers revealed 92% polymorphism among 14 tropical and temperate strains of B. mori, with average diversity index of 0.747. Stepwise multiple regression analysis (MRA) selected 35 ISSR markers positively or negatively correlated with different biomass traits, illustrated polygenic control. ISSR marker 830.81050bp was significantly associated with LWT, CWT, SWT, SR and floss content, indicated its pleiotropic role. Two ISSR markers, 835.51950bp and 825.9710bp showed significant association with floss content and TLD. These markers were segregated in F2 generation and Chi-square test confirmed (χ2 = ~45; P < 0.05) its genetic contribution to the associated biomass traits. Strains, with both positively and negatively correlated markers, had intermediate mean value for biomass traits (eg. SWT = 0.17 ± 0.014 g in GNM and Moria) indicated interaction of loci in natural populations. Low yielding Indian strains grouped together by Hierarchical clustering. Chinese and Japanese strains were distributed in the periphery of ALSCAL matrix indicated convergence of genetic characters in Indian strains. Average genetic distance between Chinese strains and Indian strains (0.193) significantly (P < 0.01) varied from that between Chinese and Japanese strains. Interaction of loci and allelic substitutions induced phenotypic plasticity in temperate B. mori populations on tropic adaptation in India. These outcomes show possibility to combine favorable alleles at different QTL to increase larval, cocoon and shell weight.

Deforestation and exploitation has led to the fragmentation of habitats and scattering of populations of the economically important eri silkworm, Samia cynthia ricini, in north-east India. Genetic analysis of 15 eri populations, using ISSR markers, showed 98% inter-population, and 23% to 58% intra-population polymorphism. Nei’s genetic distance between populations increased significantly with altitude (R2 = 0.71) and geographic distance (R2 = 0.78). On the dendrogram, the lower and upper Assam populations were clustered separately, with intermediate grouping of those from Barpathar and Chuchuyimlang, consistent with geographical distribution. The Nei’s gene diversity index was 0.350 in total populations and 0.121 in subpopulations. The genetic differentiation estimate (Gst) was 0.276 among scattered populations. Neutrality tests showed deviation of 118 loci from Hardy-Weinberg equilibrium. The number of loci that deviated from neutrality increased with altitude (R2 = 0.63). Test of linkage disequilibrium showed greater contribution of variance among eri subpopulations to total variance. D’2IS exceeded D’2ST, showed significant contribution of random genetic drift to the increase in variance of disequilibrium in subpopulations. In the Lakhimpur population, the peripheral part was separated from the core by a genetic distance of 0.260. Patchy habitats promoted low genetic variability, high linkage disequilibrium and colonization by new subpopulations. Increased gene flow and habitat-area expansion are required to maintain higher genetic variability and conservation of the original S. c. ricini gene pool.