BioMed Research International

BioMed Research International / 2015 / Article

Research Article | Open Access

Volume 2015 |Article ID 161287 | 11 pages |

Identification and Analysis of the SET-Domain Family in Silkworm, Bombyx mori

Academic Editor: Yujungeorge Zheng
Received14 Jul 2015
Revised10 Sep 2015
Accepted13 Sep 2015
Published19 Oct 2015


As an important economic insect, Bombyx mori is also a useful model organism for lepidopteran insect. SET-domain-containing proteins belong to a group of enzymes named after a common domain that utilizes the cofactor S-adenosyl-L-methionine (SAM) to achieve methylation of its substrates. Many SET-domain-containing proteins have been shown to display catalytic activity towards particular lysine residues on histones, but emerging evidence also indicates that various nonhistone proteins are specifically targeted by this clade of enzymes. To explore their diverse functions of SET-domain superfamily in insect, we identified, cloned, and analyzed the SET-domains proteins in silkworm, Bombyx mori. Firstly, 24 genes containing SET domain from silkworm genome were characterized and 17 of them belonged to six subfamilies of SUV39, SET1, SET2, SUV4-20, EZ, and SMYD. Secondly, SET domains of silkworm SET-domain family were intraspecifically and interspecifically conserved, especially for the catalytic core “NHSC” motif, substrate binding site, and catalytic site in the SET domain. Lastly, further analyses indicated that silkworm SET-domain gene BmSu(var)3-9 owned different characterization and expression profiles compared to other invertebrates. Overall, our results provide a new insight into the functional and evolutionary features of SET-domain family.

1. Introduction

SET-domain superfamily includes all but one of the methyltransferases that methylate specific site of histone lysine (K) residues involved in epigenetic regulation. This family is characterized by the highly conserved SET (Su(var)3-9, E(z), Trithorax) domain that consists of approximately 130 amino acids and is responsible for the catalytic activity of these methyltransferases. Through transferring a methyl group from S-adenosyl-L-methionine (AdoMet) to the amino group of a lysine residue on different sites of histone proteins, SET-domain methyltransferases function in making different histone methylation marks. Histone methylation is very important for the chromatin modification and regulation of gene expression [14], which plays a crucial role in the animal development and a number of other biological processes, such as heterochromatin establishment, transcription regulation, parental imprinting, and cell fate destination. Researches also suggest that SET-domain proteins are closely related to many human diseases [510].

Based on the substrate specificity towards histones, SET-domain-containing lysine methyltransferases can be classified into subgroups like KMT1 (H3K9), KMT2 (H3K4), KMT3 (H3K36), KMT4 (H3K79), KMT5 (H4K20), and KMT6 (H3K27) methyltransferases [11]. SET-domain proteins can be divided into seven main subfamilies (SUV39, SET1, SET2, EZ, SUV4-20, SMYD, and RIZ families) as well as some unclassified members like SET7/9 and SET8. Members in each subfamily have very high similarity not only in the SET-domain amino acid sequence but also in the flank motifs of SET domain.

Members of this family have been studied extensively for their function in modifying histones by methylating them directly and thus changing the mode of chromatin to regulate the binding of cofactors in many research models, such as human, mouse, Drosophila, Arabidopsis, and yeast. On the other hand, more and more recent studies give sufficient evidence of the fact that SET-domain-containing proteins also regulate many nonhistone substrates, including some transcriptional factors and tumor suppressors [1217]. Epigenetic regulation of SET-domain family remains to be studied in lepidopteron insects. Based on silkworm genome database, we identified, cloned, and analyzed the SET-domains proteins in silkworm, Bombyx mori. Here, we provide an overview of the common and unique features of silkworm SET-domain family members, which will provide more information and important reference to the whole SET-domain proteins.

2. Materials and Methods

2.1. Silkworm

The silkworm strain, Dazao (p50), used in this study is maintained by the State Key Laboratory of Silkworm Genome Biology. The silkworm larvae were reared with fresh mulberry leaves under 25°C, with a 12 h/12 h photoperiod.

2.2. Datasets

We combined the updated Silkworm 9x genomic sequencing database, silkworm EST database, CDS database, and silkworm predicted protein database (all found at SilkDB, Gene sequences of mammals, C. elegans, D. melanogaster, and some other insects were downloaded from GenBank (

2.3. Screening of SET-Domain Family Genes

Protein sequences of SET-domain family members from other species were used to query the silkworm database with -value less than 0.1. The hits in the screening were furthermore confirmed by blast in NCBI protein database. Besides, the conserved SET-domain sequences of other species were also used as query sequence to blast in silkworm database. Moreover, we used online protein domain prediction program SMART ( and Pfam ( to validate the SET domain in screening hits from silkworm database.

2.4. RNA Extraction

Different tissues of silkworm larvae were collected and stored in liquid nitrogen until use. Trizol reagent (Invitrogen, USA) was used to extract the total RNA of silkworm tissues. RNA concentration was calculated by spectrophotometer. RNA samples were digested by RNase-free DNase I (TaKaRa, Japan) to get rid of genomic DNA contamination. 2 μg of each RNA sample was used to synthesize the first strand of cDNA by M-MLV Reverse Transcriptase following the manufacturer’s instructions (Promega, USA).

2.5. Verification of Identified Genes and Expression Analysis

Primers were designed for most of the identified silkworm SET-domain family members to clone them from silkworm cDNA. Silkworm cytoplasmic gene Actin3 (forward primer: 5′-AAC ACC CCG TCC TGC TCA CTG-3′; reverse primer: 5′-GGG CGA GAC GTG TGA TTT CCT-3′) was used as an internal control. PCR amplification was performed using cDNA of deferent silkworm tissues to examine their expression profile. The 20 μL PCR reaction volume is as follows: initial denaturation at 94°C for 3 min, followed by 30 cycles of 30 s at 94°C, annealing at temperatures (usually set at 55°C) for 45 s and 1 min extension at 72°C, and extension at 72°C for 10 min. The PCR products were analyzed by 1% agarose gels.

2.6. Phylogenetic Analysis

SET-domain amino acid sequences of silkworm SET-domain family candidates were aligned to each other and also with the representative SET domains of other species by the program ClustalX. The sequences of silkworm SET-domain-containing proteins were applied to construct phylogenetic trees by neighbor-joining algorithm (1000 bootstrap replicates) with the program MEGA4.0.

3. Results and Discussion

3.1. Identification of Silkworm SET-Domain Family

We have identified the silkworm SET-domain family genes from silkworm database SilkDB for the first time. We have found 25 genes containing SET domain from silkworm genome by screening. Referring to the classification method of SET-domain family members in other species which is based on the SET-domain sequence and the feature of its flank motifs or domains, we were able to characterize 17 of them into six subfamilies of SUV39, SET1, SET2, SUV4-20, EZ, and SMYD (see Table 1). Silkworm SET-domain genes are mainly located on chromosomes 1, 3, 4, 15, 16, and 23, and 12 of them have ESTs in the database. We did not find any homologous genes of mammal RIZ subfamily in silkworm database; this may be the species difference like the SET proteins SET1, SET2 can be found in yeast while homolog genes of Su(var)3-9 and EZH1/2 are missing.

SubfamilyGene nameGene holonameName in SilkDBConserved domain
(addition to SET)
ORF length (bp)Chromosome numberScaffoldESTsProbe-valueMethylase activity in other species

SUV39BmG9a-like/Ehmt2Euchromatic histone-lysine N-methyltransferase 2BGIBMGA007949Ankyrin repeats, Pre-SET, and ALAD416115nscaf28887sw00703, sw17591,
sw15556, and sw17724
H3K9, H3K27
BmSetdb1SET-domain, bifurcated 1BGIBMGA011407Tudor, MBD, Pre-SET, and Post-SET289523 Nscaf3026Nosw21024, sw184430.0H3K9
BmSu(var)3-9Suppressor of variegation 3-9BGIBMGA011500, BGIBMGA011499CHROMO, Pre-SET, and Post-SET178223 Nscaf302712sw11427, sw044710.0H3K9

SET1BmMllMyeloid/lymphoid or mixed-lineage leukemiaBGIBMGA010221PHD, FYRN, FYRC, and Post-SET126577nscaf29862sw12266, sw149740.0H3K4
BmSet1SET-domain-containing 1BGIBMGA012978Post-SET408016 Nscaf3058Nosw14237, sw14238,
sw16463, and sw04801
BmTrrTrithorax-relatedBGIBMGA002503, BGIBMGA002502PHD, FYRN, and Post-SET9nscaf25114sw05846, sw04238,
sw06661, and sw10095
BmMll5Myeloid/lymphoid or mixed-lineage leukemia 5BGIBMGA004310PHD20nscaf2789Nosw00989, sw16148H3K4
BmSET2SET-domain-containing 2BGIBMGA003106AWS, Post-SET36334nscaf2589Nosw11671, sw11672H3K36

SET2BmNsd1Nuclear receptor binding SET-domain protein 1BGIBMGA002246PWWP, PHD, AWS, and Post-SET660026nscaf23303sw11742H3K36, H4K20
BmAsh1Absent, small, or homeotic discs 1BGIBMGA001497AWS, Post-SET, BROMO, PHD, and BAH875421nscaf21362sw115620.0H3K4, H3K9, and H4K20

SUV4-20BmSu(var)4-20Suppressor of variegation 4-20 BGIBMGA000919209413 Nscaf18984sw13536H4K20

EZBmE(z)Enhancer of zesteBGIBMGA014476SANT, CxCScaffold7076sw15702, sw162170.0H3K27

SMYDBGIBMGA014048zf-MYND1782 Nscaf31151sw10812
BGIBMGA011073115523 Nscaf3015Nosw20741
BGIBMGA00210910711nscaf2210Nosw18112, sw02496
BGIBMGA008839zf-MYND3nscaf2927Nosw17374, sw17507

OthersBmSet8/PR-set7SET-domain-containing (lysine methyltransferase) 8BGIBMGA01285391216 Nscaf30581H4K20
BmSetmarSET-domain and mariner transposase fusion geneBGIBMGA00755781015nscaf28872

3.2. Cloning and Bioinformatics Analysis

According to the gene sequence from silkworm genome database, we designed primers to clone and verify whether these genes are true hits. We have cloned most members of the silkworm SET-domain family (see Table 1). Through bioinformatics analysis, we found that SET domains of silkworm SET-domain family are intraspecifically and interspecifically conserved, especially for the catalytic core “NHSC” motif, substrate binding site, and catalytic site in the SET domain (Figure 1).

3.3. Phylogenetic Analysis of Silkworm SET-Domain Subfamilies

We have selected the representative subfamily to make phylogenetic analysis; we found that the domain architecture of SUV4-20 subfamily and SETDB1 and Su(var)3-9 of SUV39 subfamily is highly conserved to other species. The three members of silkworm SET2 subfamily are clustered to their homolog of other species, respectively.

3.3.1. SUV4-20 Subfamily

SUV4-20 subfamily methylates H4K20 in other species. Human, mouse, and Xenopus laevis have two SUV4-20 members, while in Drosophila and several other insects there is only one member. In silkworm, we identified only one SUV4-20 member like other insects. Phylogenetic analysis of SUV4-20 subfamily suggests that BmSuv4-20 clusters with Suv4-20 proteins of Drosophila, Anopheles gambiae, and other insects (Figure 2). Besides, SUV4-20 proteins of vertebrates cluster together. Among vertebrates, Suv4-20h2 (Suv4-20 homolog 2) of human and Rattus norvegicus are shorter than Suv4-20h1 (Suv4-20 homolog 1) and cluster together. BmSuv4-20 has the same domain structure as SUV4-20 proteins of other species that there is only one domain the SET located on the N-terminal of them (all shown in Figure 2).

3.3.2. SUV39 Subfamily

SUV39 methyltransferase subfamily functions in specifically methylating H3K9. In mammals, this subfamily consists of six main members, SUV39H1, SUV39H2, G9a, GLP1 (G9a-like protein 1/EuHMT1), SETDB1 (ESET), and CLLL8 (SETDB2). The SET domains of this subfamily are located on the C-terminal of the proteins which are flanked by Pre-SET and Post-SET domain. The Pre-SET domain close to SET domain is a special character of this subfamily’s proteins which assists SET domain in carrying out the enzyme activity (Figure 3). We identified three members of this subfamily in silkworm, BmG9a-like, BmSetdb1, and BmSu(var)3-9 (actually, we found two variants of Su(var)3-9 in silkworm; see details in Study of Silkworm SET-Domain Gene BmSu(var)3-9). We have not identified the highly similar homolog proteins of G9a and SETDB2 in silkworm. There is one protein that we name BmG9a-like, because it has very similar domain structure to G9a of other species including ANK (ankyrin repeats) and Pre-SET domain, while Post-SET is replaced by an ALARD domain. Since BmG9a-like shows very low sequence similarity to G9a of other species, it is alone on the phylogenetic tree. Based on all of that, we speculate that BmG9a-like is a unique member existing in silkworm. However, more evidence should be provided. SETDB1 has very special SET domain which is bifurcated by a 100–300 amino acids’ insert. We also aligned the SET-domain sequence of BmSetdb1 with other species. SET sequence of SETDB1 is highly between silkworm and other species while the insert sequences differ from one species to another (Figure 4).

3.3.3. SET2 Subfamily

SET2 subfamily of mammals mainly includes 5 members, NSD1, NSD2, NSD3, SET2 (HIF1/HYPB), and ASH1. NSD1, NSD2, and NSD3 share the conserved domains such as SET, AWS (Associated With SET), Post-SET, PWWP (a domain containing highly conserved Pro-Trp-Trp-Pro motif), PHD, and Ring finger. In silkworm, we just identified one homolog of the NSD proteins, BmNSD1, which has five of the above mentioned conserved domains except the Ring finger. However, so far, there is no homolog protein of NSD identified in Drosophila. HIF1 contains four conserved domains, SET, Post-SET, AWS, and WW. BmHIF1 identified from silkworm database have all domains but the WW domain, which may be because the sequence of BmHIF1 in silkworm is not complete, which is also confirmed by the relatively shorter sequence of BmHIF1 compared to Drosophila and other insects. ASH1 has four conserved domains except SET domain, AT hook, BROMO, BAH, and PHD. BmASH1 has all the domains except AT hook. Phylogenetic analysis of SET2 subfamily shows that members of silkworm SET2 subfamily cluster together with their homolog proteins in other species separately (Figure 5).

3.3.4. Expression Profile of Silkworm SET-Domain Genes

The expression profiles of silkworm SET-domain family members show that they have widely high expression level in gonad (testis and ovary) except BGIBMGA002076, which may be related to the gonad’s function to propagate the genetic information (Figure 6).

3.3.5. Study of Silkworm SET-Domain Gene BmSu(var)3-9

Through amplification of homolog gene of Su(var)3-9 in silkworm, we found that silkworm has two transcript isoforms. Using RACE technology to obtain the full-length cDNA sequence showed that the situation is different from other species; the two transcript isoforms have a different sequence of 846 bp at the 5′ end which belongs to the 5UTR of the longer spliceosome, and they encode a protein of 317 and 593 amino acids separately. Embryonic expression analysis of different period and different organs of three-day, fifth-instar larvae of silkworm displays that Su(var)3-9 is highly expressed during the 1–9 days of development of embryo, whereas it is relatively low expressed in the larvae organs (Figures 7 and 8).

4. Conclusion

SET-domain superfamily is a group of histone lysine methyltransferases which form different methylation marks (mono-, di-, and trimethylation, also known as me1, me2, and me3) on histone by transferring the methyl from S-adenosyl-L-methionine (AdoMet) to lysine residues of histone proteins, with production of S-adenosyl-L-homocysteine (AdoHcy). The core enzymatic sites of the methyltransferases are located in the SET domain. Since the first SET-domain protein Drosophila DmSu(var)3-9 had been identified at the end of the last century, researchers confirmed the methyltransferase activity of Su(var)3-9 protein in mammals for the first time at the beginning of the 20th century. SET-domain proteins are widely identified from the lower eukaryote yeast to higher human. Human has more than 100 SET-domain-containing proteins. The great SET-domain family is linked to its great function in epigenetic regulation. SET-domain-containing proteins play important roles in the development process of human, Drosophila, and Arabidopsis. There are very few reports on the study of histone methylation in silkworm. Currently, the only one SET-domain protein reported is BmE(z) which is the main member of PcG (Polycomb group genes) proteins. The methylation of H3K27 by BmE(z) has been confirmed in silkworm by knocking down BmE(z) which leads to the reduction of H3K27me3 level [18]. Besides, demethyltransferase has been reported in silkworm, BmLid. The function of demethylating H3K4me2, H3K4me3, H3K9me2, H3K9me3, and H3K27me3 of BmLid has been validated and it is supposed to have wider catalytic substrate than its homolog protein in mammals and Drosophila [19].

In our study, from the silkworm database, 24 genes of SET-domain family were identified and 17 out of them were included in SET1, SET2, SUV39, SUV4-20, EZ, and SMYD and six subfamilies. Except for three from SUV39 subfamily members, four from SET1 subfamily members, three from SET2 subfamily members, one from SUV4-20 subfamily members, one from EZ subfamily members, and five from SMYD subfamily members, the remaining five were independent members and two out of them were the homologous genes. Although the RIZ subfamily had not been identified in the silkworm, the certain members had already covered the common members of other species. Comparative analysis for identifying the SET-domain family members showed that they were highly conserved in each species, including the substrate binding site and catalytic site. Evolutionary analysis of the representative subfamily showed that the structure of the SUV4-20 family was highly consistent in the large number of species including the silkworm. Identified in the silkworm, Bombyx mori SUV39 subfamily members SETDB1 and su(VaR)3-9 domain composition are highly conserved among various species, especially SETDB1 unique for sequence insertion and splitting the set domain in the silkworm, Bombyx mori and other species are extremely conservative. In addition, G9a-like gene BGIBMGA007949 in the silkworm possessed the conservative composition and structure of G9a and group specificity. BGIBMGA001497, BGIBMGA003106, and BGIBMGA002246 in the silkworm were three respective homologous genes of SET2 subfamily. Interestingly, Su(VaR)3-9 in embryos from day 1 to day 9 and third day of the fifth larval stage were highly expressed, suggesting an important role possibly in the embryonic differentiation process.

In summary, based on the silkworm genome database, we identified, cloned, and analyzed the SET-domains proteins in silkworm, Bombyx mori. We intend to study the common and unique features of silkworm SET-domain family members, which will provide more information and important reference to the whole SET-domain proteins.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Authors’ Contribution

Chunqin Zheng and Hongjuan Cui contributed to the study design. Chunqin Zheng and Hailong Zhao performed the research and conducted the data analysis. Chunqin Zheng and Hailong Zhao wrote the paper. Hailong Zhao and Chunqin Zheng contributed equally to this work.


This work was supported by the National Basic Research Program of China (no. 2012cb114603), the Research Fund for the Doctoral Program of Higher Education of China (20130182110003), the Natural Science Foundation of Chongqing (cstc2013jcyjys0007), and the National Natural Science Foundation of China (no. 31172268).


  1. C. L. Peterson and M.-A. Laniel, “Histones and histone modifications,” Current Biology, vol. 14, no. 14, pp. R546–R551, 2004. View at: Publisher Site | Google Scholar
  2. J. Mellor, “It takes a PHD to read the histone code,” Cell, vol. 126, no. 1, pp. 22–24, 2006. View at: Publisher Site | Google Scholar
  3. D. Rodenhiser and M. Mann, “Epigenetics and human disease: translating basic biology into clinical applications,” Canadian Medical Association Journal, vol. 174, no. 3, pp. 341–348, 2006. View at: Publisher Site | Google Scholar
  4. B. D. Strahl and C. D. Allis, “The language of covalent histone modifications,” Nature, vol. 403, no. 6765, pp. 41–45, 2000. View at: Publisher Site | Google Scholar
  5. M. Jaskelioff and C. L. Peterson, “Chromatin and transcription: histones continue to make their marks,” Nature Cell Biology, vol. 5, no. 5, pp. 395–399, 2003. View at: Publisher Site | Google Scholar
  6. C. Martin and Y. Zhang, “The diverse functions of histone lysine methylation,” Nature Reviews Molecular Cell Biology, vol. 6, no. 11, pp. 838–849, 2005. View at: Publisher Site | Google Scholar
  7. R. Hamamoto, Y. Furukawa, M. Morita et al., “SMYD3 encodes a histone methyltransferase involved in the proliferation of cancer cells,” Nature Cell Biology, vol. 6, no. 8, pp. 731–740, 2004. View at: Publisher Site | Google Scholar
  8. M. F. Fraga, E. Ballestar, A. Villar-Garea et al., “Loss of acetylation at Lys16 and trimethylation at Lys20 of histone H4 is a common hallmark of human cancer,” Nature Genetics, vol. 37, no. 4, pp. 391–400, 2005. View at: Publisher Site | Google Scholar
  9. S. Hake, A. Xiao, and C. Allis, “Linking the epigenetic ‘language’ of covalent histone modifications to cancer,” British Journal of Cancer, vol. 90, no. 4, pp. 761–769, 2004. View at: Publisher Site | Google Scholar
  10. S. B. Hake, A. Xiao, and C. D. Allis, “Linking the epigenetic ‘language’ of covalent histone modifications to cancer,” British Journal of Cancer, vol. 90, no. 4, pp. 761–769, 2004. View at: Publisher Site | Google Scholar
  11. M. Mohan, H.-M. Herz, and A. Shilatifard, “Snapshot: chromatin lysine methylase complexes: transcriptional regulation and epigenetics,” Cell, vol. 149, no. 2, pp. 498–498.e1, 2012. View at: Publisher Site | Google Scholar
  12. K. Zhang, W. Lin, J. A. Latham et al., “The Set1 methyltransferase opposes Ipl1 Aurora kinase functions in chromosome segregation,” Cell, vol. 122, no. 5, pp. 723–734, 2005. View at: Publisher Site | Google Scholar
  13. R. Van Duyne, R. Easley, W. Wu et al., “Lysine methylation of HIV-1 Tat regulates transcriptional activity of the viral LTR,” Retrovirology, vol. 5, article 40, 2008. View at: Publisher Site | Google Scholar
  14. A. He, X. Shen, Q. Ma et al., “PRC2 directly methylates GATA4 and represses its transcriptional activity,” Genes & Development, vol. 26, no. 1, pp. 37–42, 2012. View at: Publisher Site | Google Scholar
  15. M. J. Oudhoff, S. A. Freeman, A. L. Couzens et al., “Control of the hippo pathway by Set7-dependent methylation of Yap,” Developmental Cell, vol. 26, no. 2, pp. 188–194, 2013. View at: Publisher Site | Google Scholar
  16. G. K. Dhami, H. Liu, M. Galka et al., “Dynamic methylation of numb by set8 regulates its binding to p53 and apoptosis,” Molecular Cell, vol. 50, no. 4, pp. 565–576, 2013. View at: Publisher Site | Google Scholar
  17. D. R. Calnan, A. E. Webb, J. L. White et al., “Methylation by Set9 modulates FoxO3 stability and transcriptional activity,” Aging, vol. 4, no. 7, pp. 462–479, 2012. View at: Google Scholar
  18. Z. Li, T. Tatsuke, K. Sakashita et al., “Identification and characterization of Polycomb group genes in the silkworm, Bombyx mori,” Molecular Biology Reports, vol. 39, no. 5, pp. 5575–5588, 2012. View at: Publisher Site | Google Scholar
  19. B. Zhou, X. Yang, J. Jiang et al., “Silkworm (Bombyx mori) BmLid is a histone lysine demethylase with a broader specificity than its homolog in Drosophila and mammals,” Cell Research, vol. 20, no. 9, pp. 1079–1082, 2010. View at: Publisher Site | Google Scholar

Copyright © 2015 Hailong Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

1079 Views | 463 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.