Host-response proteins identified from B. mori Accession number Total length (amino acids) Number of fully conserved substitutions Number of conserved substitutions Number of  semiconserved substitutions Total substitutions (%) Amino acid replaced at position* Conserved domains at which substitution occurred** Functions of the conserved domains**
Heat shock protein 70 ACL36370 658 354 123 48 525 (79.78%) ---- ----  
Chaperonin subunit 4 delta NP_001040107 537 263 101 28 392
(72.99%)
p.Thr217Ala p.Ile283Met p.Ala436Ser T Complex protein 1
(cd03338)
To sequester non-native proteins inside their central cavity and promote folding
SERPIN-4 ACZ81437 410 27 47 18 92 (22.44%) p.Ile143Leu p.Leu297Met SERPIN
(cd00172)
Serine protease inhibitors which regulate coagulation cascades.
Paralytic peptide precursor NP_001036883 131 21 5 17 43 (32.82%) ----- ----- ------
Prophenoloxidase activating Enzyme AAL31707 386 79 40 21 140
36.27%
p.Val229Ala Trypsin – like Serine protease
(cd00190)
Serine protease
Prophenoloxidase-2 AAG09303 693 185 142 44 371 (53.53%) p.Ile98Leu p.Trp214Tyr p.Asp281Ala p.Asp306Glu p.Asp314Glu p.Phe582Tyr Hemocyanin
(Pfam00372)
Hemocyanins and insect larval storage proteins
Aromatic-L-amino-acid decarboxylase (dopa decarboxylase) NP_001037174 478 230 96 20 346 (72.38%) p.Lys205Arg p.Arg207Lys p.Ala250Thr p.Lys359Glu p.Arg410Lys DOPA decarboxylase family
(cd06450)
Catalyzes decarboxylation of tyrosine.
Spatzle-1 precursor NP_001108066 277 11 10 05 26 (9.38%) ----- -- --
18 Wheeler precursor (Toll7-2) NP_001116821 1295 185 78 34 297 (22.93%) p.Leu240Met p.Leu293lIe p.Ile327Val p.Gln492Leu p.Ala1044Val p.Glu1137Asp Amino acids up to the position 492 are in Leucine rich repeats (Pfam13855).
Alanine at 1044 and Glutamine at 1137 are in the Toll/interleukin-1 receptor (TIR) homology domain (Pfam01582)
TIR is an intracellular signalling domain that mediates protein-protein interactions between the Toll-like receptors and signal-transduction components
Embryonic polarity protein dorsal isoform B NP_001036896 572 99 50 25 174 (30.42%) p.Ile90Leu
p.Arg128Lys
p.Ile301Val
Isoleucine and arginine are substituted in the N-terminal sub-domain of the Rel homology domain (RHD; cd07887) whereas Isoleucine at 301 is substituted in the IPT domain of Dorsal (cd01177). RHD is transcription factors, which play roles in mediating innate immunity.
Immunoglobulin-like fold, Plexins, Transcription factors (IPT) domain is part of NF kappa B transcription factor activated by stressful conditions
Nuclear factor NF-kappa-B p110 subunit isoform 1 (Relish/BmRel) NP_001095935 937 75 58 45 178 (18.99%) p.Asp134Glu
p.Asn194Asp
p.Asn224Ser
p.Thr230Ser
p.Val670Ala
Substitution upto 224th position is in N-terminal sub-domain of the Rel homology domain (RHD) (cd07884). Valine substitution is at 670 in the Transient Receptor Potential Ca2+ Channel (TIGR00870) or in the ankyrin repeats (COG0666) RHD domains are transcription factors that mediate innate immunity.
TIGR00870 are members of the family that mediate entry of extracellular Ca2+ into cells in response to depletion of intracellular Ca2+ .
Ankyrin repeats associated with signal transduction.
Cactus BAI67121 326 24 31 130 185 (56.75%) p.Ser23Thr Ankyrin repeat protein (PHA02874) Associated with signal transduction
Notch homolog NP_001157370 2463 173 47 25 245 (9.95%) p.Leu622Val p.Asp742Asn p.Thr847Ser Calcium binding EGF – like domain (smart00179; cd00054) Calcium binding site
Chitinase isoform 3 precursor NP_001166832 544 185 78 34 297 (54.59%) p.Arg173Lys p.Pro404Ala p.Thr410Ala Substitution at 173 is in the catalytically inactive chitolectin-chitotriosidase (cd02872). The domain includes a large number of catalytically inactive chitinase-like lectins
Beta-N-acetylglucosaminidase BAJ20189 633 76 66 25 167 (26.38%) p.Leu140Ile p.Ala343Thr p.Trp351Phe The substitutions are in the domain  beta-acetyl hexosaminidase and glycosyl hydrolase (cd06562; Pfam14845 & 00728; COG3525) Beta-N-acetylhexosaminidases catalyze the removal of beta-1,4-linked N-acetyl-D-hexosamine residues.
Glycosyl hydrolase (GH) family 20 is a catalytic domain with a triose phosphate isomerase (TIM) barrel fold that catalyse GH reactions
Autophagy 5-like NP_001135959 264 48 68 16 132 (50%) p.Gly211Asn p.Leu219Met Autophagy protein 5 (Pfam04106) Apg5 is directly required for the import of aminopeptidase I via the cytoplasm-to-vacuole targeting pathway.
Apoptosis-inducing factor NP_001189462 597 85 101 47 233 (39.03%) p.Gly93Ala p.Phe408Val p.Pro442Ala p.Phe449Tyr Phenylalanine substitution at 408th position is in pyridine nucleotide – disulphide oxidoreductase (Pfam07992).
Substitution at 442 and 449 are in the domain Apoptosis inducing factor (AIF; Pfam14721).
Class I and class II oxidoreductases and also NADH oxidases and peroxidases.
AIF is a bifunctional mitochondrial flavoprotein critical for energy metabolism and induction of caspase-independent apoptosis
Caspase Nc NP_001182396 438 21 34 26 81 (18.49) p.Leu297Phe p.Asp359Asn Caspase domain, interleukin-1 beta converting enzyme (ICE) homologues (cd00032’ smart 00115; pfam006560 Cysteine-dependent aspartate-directed proteases that mediate programmed cell death (apoptosis)
Trehalase-2 BAE45249 642 168 99 47 314 (48.91%) p.Asn148Ser p.Val161Phe p.Phe199Leu p.Val391Leu Trehalase (pfam 01204) It is known to recycle trehalose to glucose and protects of  proteins and membranes against a variety of stresses
*p.Thr215Ala represents threonine at 215th position is substituted by alanine; **Accession number and details from NCBI Conserved domains
Table 3: Key to amino acid substitutions and conservations in the host- response proteins of B. mori when aligned with the similar proteins from representative species of other insect orders presented in supplementary Table 2.