The molecular mechanism involved in BmNPV resistance was investigated using a genome wide microarray in midgut tissue of Indian silkworm Bombyx mori. In resistant race (Sarupat), 735 genes up-regulated and 589 genes down-regulated at 12 h post BmNPV infection. Similarly, in case of susceptible race (CSR-2), 2183 genes up-regulated and 2115 genes down-regulated. Among these, nine up-regulated and eight down-regulated genes were validated using real-time qPCR analysis. In Sarupat, vacuolar protein sorting associated, Xfin-like protein and carboxypeptidase E-like protein genes significantly up-regulated in infected midgut; prominently down-regulated genes were glutamate receptor ionotropic kainite 2-like, BTB/POZ domain and transferrin. Considerably up-regulated genes in the CSR-2 were peptidoglycan recognition protein S6 precursor and rapamycin while the conspicuous down-regulated genes were facilitated trehalose transporter and zinc transporter ZIP1-like gene. The up-regulation of genes in resistant race after BmNPV infection indicates their possible role in antiviral immune response.

Cecropins constitute one of the largest and most potent immune protein families found in insect species with diversified numbers and features. In view of the large number of cecropin proteins existing with much sequence variations among them, an overview of the multigene cecropin family in silkworm Bombyx mori was attempted in this study. Cecropin encodes an inducible 64 residue anti‐bacterial peptide and was clustered into two groups; first group viz. A and second group including B, D, E and Enbocin. Cecropin A consisted of two sub-groups located on chromosome number 6 of B.mori genome. Cecropin B consisted of six sub‐groups, cecropin D and E of one each and Enbocin of two. The second sub‐group formed in tandem array of multigene family locus over a length of 78.62 kb on chromosome number 26 in B.mori genome and was organized in positive as well as opposite orientation. The results indicated that cecropin B genes were organized in a close cluster with the intergenic sequence ranging from 1366 bp to 23526 bp. Interestingly a distantly related cecropin E was also located within the cecropin B multigene locus. Similarly distant members like cecropin D and Enbocin were also located in the 3’ region of cecropin B locus. The maximum intergenic region of 23526 bp observed between Cecropin D and Enbocin indicates that the two genes were distantly evolved. The phylogenetic analysis clearly indicates a positive correlation between the clusters and physical location on the chromosome, as the length of the intergenic region plays a major role to create newer cecropin families. EST database analysis suggests that most of the cecropin A members were expressed in the microbial fat body while, the cecropin B was equally expressed in fat body and other target tissues. The signal peptides were conserved in all the twelve paralogous gene sequences.

The activation of host response proteins against parasitic infection is dependent on the coregulation of immune gene expression. The infection of commercially important silkworm Bombyx mori through endoparasite Exorista bombycis enhanced host-response gene expression in integument early in the infection and was lowered asymptotically. Principal component analysis (PCA) showed heterogeneity while explaining ∼80 % variance among expression timings. PCA showed positive and negative correlation with gene expression and differentiated transcriptional timings, and revealed cross talk within the immune system. Pearson correlation analysis showed significant linear correlation (mean R (2) = >0.7; P < 0.004) between the expression of 16 pairs of genes in control, while the relation switched over to curvilinear due to parasitism. The genes showed pleiotropic interaction among them, with four genes each for prophenoloxidase activating enzyme (PPAE) and caspase. Besides, after parasitism, exclusive correlation of five gene pairs including PPAE-Spatzle pair (R (2) = 0.9; P < 0.011) was observed in the integument. In integument, the phenol oxidase (PO) activity showed a positive correlation with the tyrosine level (R (2) = 0.410; P < 0.002) and a curvilinear relation (R (2) = 0.745; P < 0.0002) with the expanding lysis area. The PO activity was positively correlated with BmToll expression and negatively correlated with paralytical peptide expression, revealing polygenic influence. Caspase expression was tightly regulated by signal genes in control integument, whereas they were deregulated after infection. Switchover from linear to curvilinear correlation and the appearance of new gene correlations in parasitized integument revealed deviation from gene coregulation, leading to impaired immune responses, characterized by lowered gene expression and varied phenotypic consequences.

The bipartite genome of an Indian isolate of Bombyx mori bidensovirus (BmBDV), one of the causative agents of the fatal silkworm disease 'Flacherie', was cloned and completely sequenced. Nucleotide sequence analysis of this Indian isolate of BmBDV revealed two viral DNA segments, VD1 and VD2 as well as a DNA polymerase motif which supports its taxonomical status as the type species of a new family of Bidnaviridae. The Indian isolate of BmBDV was found to have a total of six putative ORFs four of which were located on the VD1 with the other two being on the VD2 DNA segment. The VD1 DNA segment was found to code for three non-structural proteins including a viral DNA polymerase as well as one structural protein, while the VD2 DNA segment was found to code for one structural and one non-structural protein, similar to that of the Japanese and Zhenjiang isolates of BmBDV. A BmBDV ORF expression study was done through real time qPCR wherein the VD2 ORF 1 and 2 showed the maximum transcript levels. This is the first report of the genome characterization of an Indian isolate of BmBDV, infecting silkworm B. mori.

Attacin gene exists as paralogous conversion and is being used for identification of strain variations in insects based on the sequence variation. Hence, a study was undertaken to analyze the sequence variation of the attacin gene isoforms in the tasar silkworm Anthereae mylitta that exists in the form of different ecoraces depending upon the environment, food plant and location. Comparison of the previously reported attacin sequences with the DNA sequences of attacin A and B genes revealed six amino acid substitutions among the sequences of the ecoraces which however did not affect the functional domain of Attacin. The generated dendrogram clearly indicated unique branches for each ecorace with two separate gene clusters for attacin A and B. The Sarihan ecorace formed a separate sub-group under both the gene clusters. The present study also revealed the presence of Attacin_N Superfamily domain exclusively in Exon I separated from the Attacin_C Superfamily domain that was present in Exon II and part of Exon III, a prominent character of attacin gene. The phylogenetic reconstruction analysis of attacin gene in A.mylitta supported the common evolutionary origin of attacin genes belonging to the Lepidoteran and Dipteran families that formed two separate clusters.