1. Introduction
With the publication of the silkworm (Bombyx mori) genome sketch in 2004 [1], a new era of silkworm research began, followed by a series of innovative research efforts that yielded genome fine maps [2], genetic variation maps [3], gene chips [4], and epigenetic maps [5], elevating silkworm research. In this context, transgenic silkworm technologies have been gradually established and have contributed immensely to the study of silkworm functional genomics [6,7,8]. CRISPR/Cas9 gene editing technology has become one of the important tools for studying silkworm functional genomics and utilizing it to improve and cultivate new high-quality silkworm strains is paramount.
The economic value of silkworms is primarily reflected in the application of silk proteins. The silk gland is the only organ in the silkworm capable of secreting silk proteins and is divided into the anterior (ASGs), middle (MSGs), and posterior (PSGs) silk glands according to their function and structure. The MSG mainly secretes sericin, which comprises sericin 1 (Ser1), 2 (Ser2), and 3 (Ser3) proteins according to their molecular weights [9]. The PSG is capable of secreting silk fibroin proteins, which mainly include silk fibroin heavy chain protein (Fib-H), silk fibroin light chain protein (Fib-L), and P25 protein [10,11,12]. The silk fibroin proteins secreted by the PSG are enveloped in the silk sericin protein secreted by the MSG, compressed by the ASG, and ejected from the body to form silk. The development of silk glands in the silkworm can be roughly divided into three stages: silk gland development in the embryonic and larval stages and silk gland degradation in the pupal stage [13,14,15]. In the embryonic stage, silk gland cells undergo approximately 10 mitotic divisions and differentiation of the MSG and PSG to form fully developed silk gland tissue approximately 8 days after egg laying [16]. Cyclins B and E and Scr have been shown to be involved in this process [17,18]. In the larval stage, the silk gland no longer undergoes mitosis and increases its DNA content only through intranuclear replication [19,20]; it has been shown that Let7 MicroRNAs (Let7) knockout promotes intranuclear replication in silk gland cells [21,22]. Apoptosis in the silk gland is thought to be regulated by ecdysteroids (insect growth and molting hormones) and involves both apoptosis and autophagy [23,24,25]. Silk gland development and silk protein synthesis require the involvement of a large number of genes. Although several studies have been conducted, the mechanisms by which the relevant genes regulate silk gland development or affect silk protein synthesis remain unexplored. Studying silk gland development and silk protein synthesis mechanisms will help us better understand the relationship between silk gland development and silk yield, provide theoretical guidance for the genetic improvement of high silk-yielding silkworm varieties, and facilitate the research and development of new materials.
Ecdysteroid kinase-like-1 (EcKL1) is a member of the ecdysteroid kinase-like family (EcKLs), which contains a Zinc finger C4 and HLH domain containing the kinase domain subfamily of choline kinases (CHK kinase-like) with a structural domain similar to ecdysteroid 22 kinase (EcKinase). Most insects contain 12–105 genes encoding this family. This family is believed to include kinases that act on other small molecule substrates and may function in detoxification processes [26]. Mutations in five EcKLs in the Drosophila One ancestral EcKL clade in the Drosophila (Dro5) branch significantly reduce developmental tolerance to caffeine and the fungal secondary metabolite tretinoin in Drosophila [27]. Overexpression of Juvenile hormone-inducible protein 26 (JHI-26) in the testes of male flies leads to paternal-associated lethality and reduces mating acceptability by females with whom they mate [28,29]. EcKLs genes do not seem to have received much attention in the silkworm, except for the identification of EcKinase, an enzyme responsible for phosphorylating ecdysteroids at C-22 in the ovary to form the physiologically inactive ecdysteroid 22-phosphates [30,31,32]. Other EcKL genes have remained uncharacterized until recently when researchers discovered that the silk gland of the silkworm is affected by Let-7, a gene that affects silk gland development. Let-7 can negatively regulate EcKinase by binding to it [21], suggesting that other EcKLs possessing structural domains similar to those of EcKinase have potential regulatory roles in the growth and development of silkworm silk glands.
This study aimed to determine the role of BmEcKL1 in silk gland development and silk protein synthesis in silkworms. We screened the potential target gene BmEcKL1 affecting silk gland development by analyzing the spatiotemporal expression of 14 EcKLs genes in the silkworm. We successfully knocked out BmEcKL1 in the MSG and PSG using CRISPR/Cas9 gene editing technology, demonstrating that BmEcKL1 is involved in regulating the growth and development of silkworm silk glands, and simultaneously affects silk protein synthesis to a certain extent.
2. Results
2.1. Expression Analysis of EcKL Genes in the Silkworm
EcKL genes seem to be rich in functionality, but very few studies focus on the silkworm. Based on the known information, it is difficult to select suitable research genes from numerous genes. Therefore, it is necessary to analyze the spatiotemporal expression of the EcKL family to help us select the genes for study. The expression profile during the entire silkworm period (Figure 1a) showed that most EcKL genes had low expression from the newly-hatched silkworm to the third instar, and expression suddenly increased from the third instar of sleeping silkworms until it reached a maximum during the wandering period. The expression profile of the 5th-instar day 3 tissues showed that the EcKL genes were highly expressed primarily in the head, and most were expressed in the malpighian-tubule, midgut, MSG, and PSG (Figure 1b). LOC10174263, LOC101743675, LOC101745199, and LOC101745257 were highly expressed in either the midgut or the malpighian-tubule, implying that they could be related to detoxification. In contrast, EcKL1 and EcKinase expression was lower in both tissues, but EcKL1 exhibited superior expression in the ASG, MSG, and PSG. We mapped the expression profiles of the EcKLs in the ASG, MSG, and PSG at the fifth instar and revealed that EcKL1, LOC101742632, and LOC101745257 were all expressed in the ASG (Figure 1c). In addition, the MSG full-period expression profile (Figure 1d) demonstrated high expression of EcKL1 and LOC101741944, while that of the PSG (Figure 1e) exhibited a more homogeneous expression pattern. Therefore, according to these results, we selected the EcKL1 gene, which had a good expression profile in the silk gland, as the target of the next study.
2.2. BmEcKL1 Knockout in MSG and PSG
The BmEcKL1 gene has a full CDS length of 1230 bp and contains two exons. Figure 2a shows the knockout vector structure diagram. The constructed vector was injected into silkworm eggs within 2 h after the eggs were laid, and the injected eggs were referred to as the G0 generation. The G0-generation silkworms were reared to adulthood and mated with wild-type silkworms (WT), and the laid eggs were referred to as the G1 generation. The eggs were identified by detecting green fluorescent protein (EGFP) expression in the eyes under the control of neuron-specific 3 × P3 promoter on the sixth day after the eggs of the G1 generation were laid. The G1 generation positive for green fluorescent eyes was reared to adulthood and rescreened (Figure 2c), and a positive gRNA strain with green fluorescent eyes was obtained and designated BmEcKL1-gRNA.
The BmEcKL1-gRNA G1 positive strains were then crossed with PSG-Cas9 and MSG-Cas9 [33,34] strains to produce F1 generations, which were screened for green and red fluorescent markers in the eyes of late-stage embryos and adults (Figure 2h,i); the positive F1 generations obtained were designated ΔBmEcKL1-PSG and ΔBmEcKL1-MSG, respectively. To verify whether mutation occurred in the BmEcKL1 gene sequence, the fragment containing the sgRNA target site was PCR amplified with site-specific primers. The results showed that various forms of DNA fragment deletions occurred at the sgRNA loci in ΔBmEcKL1-MSG- and ΔBmEcKL1-PSG-positive individuals (Figure 2d,e). Subsequently, we conducted experiments using mutants with deleted bases greater than 15 and not multiples of 3. Further analysis of the sequencing results showed that the individual knockout efficiency was 31% (Figure 2f), of which 28% were fragment deletions, 2% were fragment insertions, and 1% were base substitutions. Finally, the expression of BmEcKL1 at the transcriptional level was measured in the WT, ΔBmEcKL1-MSG, and ΔBmEcKL1-PSG strains using EcKL1 mutant primers (Figure 2g), which showed that the BmEcKL1 was successfully knocked out in both ΔBmEcKL1-MSG and ΔBmEcKL1-PSG individuals.
2.3. Deletion of BmEcKL1 Leads to MSG Abnormalities and Changes in Silk Protein Production
WT and ΔBmEcKL1-MSG individuals were mixed and reared to the fifth instar, and no significant differences in the changes in individual size and body weight were observed on day 3 of the fifth instar (5L3D; Figure 3a,e). Next, the larvae were dissected at different time points at the fifth instar stage, and the silk glands were assessed and visualized using a camera. At day 1 of the fifth instar (5L1D), the MSG size of the ΔBmEcKL1-MSG strain was almost the same as that of the WT (Figure 3b(i)) but slightly larger than that of the WT at 5L3D (Figure 3b(ii)). By 5L5D, the MSGs of the ΔBmEcKL1-MSG and WT clearly differed, with longer and larger ΔBmEcKL1-MSG M-MSGs and P-MSGs (Figure 3b(iii)). Compared with the WT, the ΔBmEcKL1-MSG P-MSG was significantly thicker and appeared roughly serrated on the surface by 5L6D (Figure 3b(iv)).
For a more comprehensive assessment, we measured the MSG length and weight at 5L6D. The MSG of the ΔBmEcKL1-MSG strain increased by approximately 10% (Figure 3f), while its weight increased by approximately 30% relative to those of the WT (Figure 3g). To understand whether the changes in the silk glands affect silk protein secretion, we determined cocoon and pupa phenotypes after clustering, and the results showed that the cocoons of the ΔBmEcKL1-MSG strain were slightly larger than those of the WT (Figure 3c). In addition, the full cocoon weight and layer weight increased to some extent (Figure 3h,j), but the difference in pupal size and weight between the two strains did not significantly differ (Figure 3d,i). Therefore, the cocoon layer ratio in the ΔBmEcKL1-MSG strain marginally increased (Figure 3k).
2.4. BmEcKL1 Knockout Leads to Enlarged PSG and Increased Silk Yield
WT and ΔBmEcKL1-PSG were mixed and reared to the fifth instar. We recorded the changes in individual size and body weight of the WT and ΔBmEcKL1-PSG strain at 5L3D, and the results showed that none of them appeared to be significantly altered (Figure 4a,e). To determine whether BmEcKL1 knockout affects PSG development, we dissected and visualized the silk glands at the fifth instar. The results showed that at 5L1D, the PSG of the ΔBmEcKL1-PSG strain had changed minimally compared with that of the WT (Figure 4b(i)), and from 5L3D onwards, the change was more apparent. The PSG of the ΔBmEcKL1-PSG strain was longer than that of the WT (Figure 4b(ii)), and by 5L5D, it had become longer and curlier (Figure 4b(iii)), becoming significantly longer and curlier by 5L6D (Figure 4b(iv)). Moreover, we measured the PSG length and weight at 5L6D of the ΔBmEcKL1-PSG strain. We observed that the former increased by approximately 20% (Figure 4f) and the latter by approximately 38% (Figure 4g) relative to those of the WT. Furthermore, we assessed the cocoon and pupa phenotypes and found that the ΔBmEcKL1-PSG cocoons were larger than those of the WT (Figure 4c), with an increase in full cocoon weight of approximately 10% (Figure 4h). The cocoon layer weight was also significantly increased (Figure 4j); nevertheless, pupae did not significantly differ (Figure 4d,i). Thus, the cocoon layer ratio of the ΔBmEcKL1-PSG strain was significantly increased (Figure 4k).
2.5. BmEcKL1 Knockout Alters Silk Protein Secretion
To investigate whether BmEcKL1 knockout affects the secreted silk protein, we analyzed the cocoons of the WT, ΔBmEcKL1-MSG, and ΔBmEcKL1-PSG strains by scanning electron microscopy (SEM). The results showed that the silk fibers of the ΔBmEcKL1-MSG and ΔBmEcKL1-PSG strains were thicker than those of the WT. The outer and inner layers of both the ΔBmEcKL1-MSG and ΔBmEcKL1-PSG cocoons were denser than those of the WT (Figure 5a,b), whereas those of the ΔBmEcKL1-PSG strain were less dense relative to those of the ΔBmEcKL1-MSG strain. This was corroborated by the observation of single filament fibers (Figure 5c).
The single filament fibers of the ΔBmEcKL1-MSG and ΔBmEcKL1-PSG strains were coarser than those of the WT, while those of the ΔBmEcKL1-PSG strain were coarser than those of the ΔBmEcKL1-MSG strain. In addition, the cross-sections of the ΔBmEcKL1-MSG and ΔBmEcKL1-PSG cocoons were thicker than those of the WT cocoons (Figure 5d). BmEcKL1 knockout appeared to enhance the performance of silk. To further investigate the reasons behind this phenomenon, we determined the expression of Fib-H, Fib-L, P25, sericin-1, and sericin-2, the main silk proteins secreted by silk glands, at the transcriptional level at 5L6D. The results showed no significant change in fibroin in the ΔBmEcKL1-MSG strain, whereas sericin-1 and sericin-2 were significantly increased (Figure 6a). In the ΔBmEcKL1-PSG strain, three fibroin proteins, Fib-H, Fib-L, and P25, were upregulated to varying degrees. However, the change in sericin was not significant (Figure 6b). Finally, we determined the transcriptional expression of genes related to silk gland development and hormone metabolism. The results showed that BmNvd, Bmsage, and BmSro expression increased in the ΔBmEcKL1-MSG strain, whereas BmRPL3 expression was downregulated, and BmSpo and Let7 did not significantly change. BmNvd and Bmsage expression was upregulated in the ΔBmEcKL1-PSG strain, whereas that of the others did not change significantly.
3. Discussion
In this study, we demonstrated that BmEcKL1 knockout affects silk gland development as well as silk protein synthesis. Using CRISPR/Cas9 technology, BmEcKL1 was knocked out in the MSG and PSG. On the one hand, BmEcKL1 knockout resulted in an enlarged and less smooth MSG and increased sericin protein production. However, although the cocoon size increased, the cocoon layer ratio only increased slightly. On the other hand, BmEcKL1 knockout resulted in an elongated and curlier PSG and higher fibroin protein production. Furthermore, the cocoon layer ratio considerably improved with the increase in cocoon size.
The silkworm is a valuable economic insect and model lepidopteran, occupying an important position in human history, culture, and economic life. Silkworms are characterized by silk ejection and cocooning, and the growth and development of the silk gland, the only silk-producing organ of silkworms, is closely related to silk protein synthesis. Silk glands are larger in high silk-producing strains compared with low silk-producing strains [35]. Transgenic techniques have revealed that Ras1 or Yorkie overexpression can increase silk gland size and improve silk production [36,37]. However, silk gland development is an extremely complex and synergistic process that involves hormone regulation, energy metabolism, material conversion, and other important processes requiring the participation of a large number of genes. The successive expression of genes, such as cyclins B and E and SGF1, regulates the formation and extension of the silk gland in the embryonic stage. The development of the silk gland in the larval stage is primarily affected by nuclear replication and silk gland cell size. Let7 knockout promotes the nuclear replication of silk gland cells, leading to silk gland enlargement and increased silk production. Collectively, both silk gland cell volume expansion and the intranuclear replication of DNA ultimately prepare the organ for the synthesis of large amounts of silk proteins at the fifth instar. In the present study, we knocked out BmEcKL1 in the MSG and PSG using CRISPR/Cas9 technology, which led to their expansion and improved silk production. Nevertheless, the ΔBmEcKL1 MSG strain was inferior to the ΔBmEcKL1 PSG strain in terms of the increase in silk yield. The difference between the two strains may be attributed to the fact that fibroin proteins account for a greater proportion of silk fibers than do sericin proteins, and the increase in fibroin can improve the cocoon layer ratio more under the same conditions.
The EcKLs have a highly conserved cluster, HxDhx3Nh3… D, known as the Brenner motif (where h denotes a large hydrophobic amino acid and x denotes any amino acid) [38,39], that is an important structure for this family of proteins to function, shared by several protein kinases [40], choline kinases [41,42], and aminoglycoside phosphotransferases [43]. EcKinase, as a representative of the EcKLs, has previously been shown to bind directly to ecdysone. Before thymus differentiation, free ecdysteroids are involved in controlling early events in embryonic development, including embryonic epidermogenesis and germ band lengthening, and are thought to be released from ovarian-derived ecdysteroid phosphates [44]. Considering how ecdysteroid phosphate is produced, researchers have suggested that in mature oocytes, EcKinase, which is responsible for the phosphorylation of ecdysteroid at the C-22 position, converts ecdysteroids into inactive phosphate. The larval EcKL expression profiles throughout the period exhibited some regularity. The analysis of BmEcKL1 expression showed that starting from the third instar, BmEcKL1 expression increased at the late stage of each instar, both at the early stage of molting and at the pupal transformation stage, which seemed to be similar to the periodic change of the 20-hydroxyecdysone (20E) titer. The transcript levels of silk gland development-related genes and 20E synthase were examined in the ΔBmEcKL1-MSG and ΔBmEcKL1-PSG strains. The results showed that 20E synthesis-related genes were elevated in both knockout lines, implying that BmEcKL1 might be directly or indirectly regulated by 20E, in addition to Bmsage, which is thought to control the number of silk gland cells through the cell cycle pathway. The expression of this gene in the knockout lines appeared to be altered by BmEcKL1 knockout; therefore, we hypothesized that BmEcKL1 knockout altered the energy metabolism pathway of the relevant hormones in the silkworm, causing changes in the relevant genes regulating silk gland development, such as Bmsage [45].
Silk is a natural animal protein, mainly comprising two types of proteins: fibroin and sericin. On the one hand, fibroin is a high molecular weight insoluble fibrous protein, and the current silk industry mainly utilizes it to produce raw silk and silk products. On the other hand, sericin is a low molecular weight soluble globular protein that has moisturizing, antioxidant, anti-bacterial, and UV-resistant functions, among others, and has high development and utilization value in the fields of cosmetics, medicine, food, and fiber-modification materials [46,47]. Fibroin and sericin account for 70–80% and 20–30%, respectively, in conventional varieties [48]. The cocoons of both ΔBmEcKL1-MSG and ΔBmEcKL1-PSG strains were enlarged to a certain extent, and SEM revealed that the silk fibers of both strains were coarser than those of the WT, with those of the ΔBmEcKL1-PSG strain being the coarsest. Thus, the ΔBmEcKL1-PSG strain silk fibers have superior toughness, whereas those of the ΔBmEcKL1-MSG strain have wider utility in more fields.
In conclusion, silk yield is paramount to the silk industry. Previously, high silk-yielding silkworm strains were obtained through traditional genetic breeding methods, which are labor- and time-intensive. However, the development and application of gene editing technology [49,50,51] have enabled researchers to produce silkworms with superior economic traits more conveniently, quickly, and efficiently by regulating silk gland development and silk protein synthesis. In this study, we demonstrated, that BmEcKL genes are involved in the regulation of silk gland development and silk protein synthesis, providing a new direction for future research on EcKLs in insects and a theoretical basis for the mechanism of silk gland development and silk protein synthesis mechanisms. Furthermore, our findings provide a reference for the selection and breeding of high silk-yielding silkworms. However, further studies are needed to fully elucidate how BmEcKL1 is involved in regulating silk gland development and silk protein synthesis. Meanwhile, whether knockout efficiency can be improved by designing multiple SgRNA tandem vectors, which will lead to greater changes in the silk glands and higher silk production, remains to be verified. Finally, whether other EcKLs genes also regulate silk gland development needs to be further investigated.
4. Materials and Methods
4.1. Bioinformatics Analysis
The full-length cDNA and protein sequences (ACCESSION: XM_012688665, XP_012544119) of BmEcKL1 were obtained from the NCBI database. BmEcKL1 was characterized and sequenced using ExPASy. Available online:
4.2. Silkworm Rearing
The multivoltine, non-lagging silkworm strain D9L was used in this study. The silkworms were reared at the Biological Research Center, Institute of Frontier and Interdisciplinary Sciences, Southwest University of China. The rearing of silkworms was divided into 10 zones, with 150 silkworms per zone. Silkworm larvae were reared on fresh mulberry leaves (Morus spp.) under photoperiod conditions of 12 h light/12 h dark at 25 ± 1 °C and 70 ± 5% humidity.
4.3. RNA Extraction and RT-qPCR Validation
Silkworms from a certain zone were randomly selected for sampling, and 5 silkworms from each period and tissue were selected and mixed as one sample. Sample RNA was extracted using TRIzol™ reagent (Invitrogen, Carlsbad, CA, USA), and reverse transcription was performed using a PrimeScript™ RT reagent Kit with gDNA Eraser (Takara Bio Inc., Shiga, Japan). All samples were stored at −80 °C. Three independent replicates were performed for each RNA extraction. Primers were designed using Primer Premier 5.0; all designed primers are listed in Supplementary Table S1. Quantitative PCR (qPCR) was performed using SYBR®Premix Ex Taq™II (Takara Bio Inc., Shiga, Japan) under the following conditions: 95 °C for 10 s, followed by 40 cycles of 95 °C for 5 s, and 60 °C for 30 s. B. mori sw22934 [55] was used as an internal reference gene. Three independent replicates were performed for each RT-qPCR experiment.
4.4. Vector Construction and Microinjection
The BmEcKL1 exon sequence was selected as the target region for knockout. SgRNA was designed using CCtop (available online:
4.5. Phenotype Statistics and Observation
The transgenic silkworms were screened using a fluorescence microscope (Olympus, Tokyo, Japan), and the genome was extracted as a template. PCR amplification was performed with BmEcKL1 target site detection primers (Table S1). The PCR fragment was subjected to agarose gel electrophoresis, purified, connected to T vector, transformed in E. coli, and single clones were sequenced by Sanger method using M13F/R primer. WT and transgenic silkworms were reared to the fifth instar for dissection, after which the silk glands were excised, washed with phosphate-buffered saline, and imaged using a fluorescence microscope (Olympus, Tokyo, Japan). The WT and transgenic silkworms were mixed and reared to the fifth instar, and their body weights were counted. The silkworms were reared to the moth stage, and cocoons were cut in advance to identify males and females.
4.6. Scanning Electron Microscopy
Cocoon morphology was assessed using a scanning electron microscope (JSM-5000; JEOL, Tokyo, Japan). The cocoon samples of each group were glued to conductive carbon tape and coated with a thin layer of platinum (NeoCoater MP-19020NCTR; JEOL). SEM observations were performed at 25 ± 1 °C.
4.7. Statistics and Reproducibility
All statistical analyses were performed using GraphPad Prism 8.0 (GraphPad Software, La Jolla, CA, USA), and the data are presented as the means ± SEM. The mean cycle threshold (Ct) value was converted to the relative expression level using the 2−ΔΔCt method [58,59]. Statistical analyses of the expression levels were performed using a two-tailed unpaired student t-test. Significant differences were defined as * p < 0.05, ** p < 0.01, *** p < 0.001, **** p < 0.0001.
Conceptualization, S.L. and Q.X.; methodology, S.L. and Q.X.; validation, S.L., J.L. and Y.S.; software, S.L.; formal analysis, S.L. and J.L.; investigation, S.L., J.L. and Y.S.; resources, X.H., P.L., F.W., G.S., P.Z. and Q.X.; data curation, S.L.; writing—original draft preparation, S.L. and Q.X.; writing—review and editing, S.L., J.L., Y.S., X.H., P.L., F.W., G.S., P.Z. and Q.X.; visualization, S.L.; supervision, P.Z. and Q.X.; project administration, S.L., X.H., P.L., F.W. and G.S.; funding acquisition, P.Z. and Q.X. All authors have read and agreed to the published version of the manuscript.
Not applicable.
Not applicable.
All of the data used in this study have been provided in the main text and the
The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this study.
Footnotes
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Enlarge this image.Figure 1. Spatiotemporal expression of EcKLs. EcKL1 exhibited superior expression in the ASG, MSG, and PSG. (a) Period expression profile of EcKLs during the entire larval stage. (b) Tissue expression of EcKLs in the day 3 fifth instar. (c) Fifth-instar temporal expression of EcKLs in the anterior silk gland. (d) Fifth-instar temporal expression of EcKLs in the middle silk gland. (e) Fifth-instar temporal expression of EcKLs in the posterior silk gland.
Enlarge this image.Figure 2. Construction of sgRNA strains and mutants. Relative expression means the expression of the experimental gene relative to the reference gene sw22934. (a) Schematic diagram of the sgRNA transgenic vector. (b) The designed sgRNA is located in the first exon. (c) Positive individuals selected at the G1 embryo and adult stages by screening for the EGFP marker that emit green light in the eyes. (d) Sequence alignment of sgRNA-targeted genomic regions in the ΔBmEcKL1-MSG strain exhibits multiple forms of mutations. (e) Sequence alignment of sgRNA-targeted genomic regions in the ΔBmEcKL1-PSG strain exhibits multiple forms of mutations. (f) Ratios of different mutation types. (g) Expression of BmEcKL1 at the transcriptional level. (h) Positive individuals of the ΔBmEcKL1-MSG strain were selected at the F1 embryo and adult stages by screening for green-emitting EGFP and red-emitting RFP markers in the eyes. (i) Positive individuals of the ΔBmEcKL1-PSG strain were selected at the F1 embryo and adult stages by screening for green-emitting EGFP and red-emitting RFP markers in the eyes. Significant differences were defined as *** p < 0.001, **** p < 0.0001.
Enlarge this image.Figure 3. Phenotypic statistics of the ΔBmEcKL1-MSG strain. For each set of experiments, 50 WT and 50 mutant silkworms were counted and performed three times in total. (a) Individual size differences between the WT and ΔBmEcKL1-MSG; there is no significant difference. (b) Silk gland phenotype of the ΔBmEcKL1-MSG at the fifth instar ((i) 5L1D; (ii) 5L3D; (iii) 5L5D; (iv) 5L6D); the mutant has a larger MSG. (c) The cocoon phenotype is slightly larger in the mutant pupae than in the WT. (d) The pupae phenotype; there is no significant difference. (e) WT and ΔBmEcKL1-MSG strains were weighed at the fifth-instar larva stage; there is no significant difference. (f) MSG length of the ΔBmEcKL1-MSG strain. (g) MSG weight of the ΔBmEcKL1-MSG strain. (h) Total cocoon weight of the ΔBmEcKL1-MSG strain. (i) Pupae weight of the ΔBmEcKL1-MSG strain. (j) Cocoon shell weight of the ΔBmEcKL1-MSG strain. (k) Cocoon shell ratio of the ΔBmEcKL1-MSG strain. Significant differences were defined as ** p < 0.01, *** p < 0.001, **** p < 0.0001. The meaning of ns is that there is no significant difference.
Enlarge this image.Figure 4. Phenotypic statistics of the ΔBmEcKL1-PSG strain. For each set of experiments, 50 WT and 50 mutant silkworms were counted and performed three times in total. (a) Individual size difference between the WT and ΔBmEcKL1-PSG; there is no significant difference. (b) Silk gland phenotype of ΔBmEcKL1-PSG at the fifth-instar larva stage ((i) 5L1D; (ii) 5L3D; (iii) 5L5D; (iv) 5L6D); the mutant has a longer PSG. (c) The cocoon phenotype is significantly larger in the mutant pupae than in the WT. (d) The pupae phenotype; there is no significant difference. (e) WT and ΔBmEcKL1-PSG strains were weighed at the fifth-instar larva stage; there is no significant difference. (f) PSG length of the ΔBmEcKL1-PSG strain. (g) PSG weight of the ΔBmEcKL1-PSG strain. (h) Total cocoon weight of the ΔBmEcKL1-PSG strain. (i) Pupa weight of the ΔBmEcKL1-PSG strain. (j) Cocoon shell weight of the ΔBmEcKL1-PSG strain. (k) The cocoon shell ratio of the ΔBmEcKL1-PSG strain. Significant differences were defined as ** p < 0.01, *** p < 0.001, **** p < 0.0001. The meaning of ns is that there is no significant difference.
Enlarge this image.Figure 5. Scanning electron microscopy images of the WT, ΔBmEcKL1-MSG, and ΔBmEcKL1-PSG cocoons. (a) The first column shows the outer surface; the mutant outer surface is denser than the WT. (b) The second column shows the inner surface; the mutant inner surface is denser than the WT. (c) The third column shows a single fiber; the mutant single fiber is coarser than the WT. (d) The fourth column shows the cross-section, the mutant has a thicker cross-section than the WT.
Enlarge this image.Figure 6. Transcriptional level analysis of related genes. For each set of experiments, 5 WT and 5 mutant silkworms were counted and performed three times in total. Relative expression means the expression of the experimental gene relative to the reference gene sw22934. (a) The transcription levels of major silk protein genes in the ΔBmEcKL1-MSG strain. (b) The transcription levels of major silk protein genes in the ΔBmEcKL1-PSG strain. (c) The transcription levels of genes involved in silk gland development in the ΔBmEcKL1-MSG strain. (d) The transcription levels of genes involved in silk gland development in the ΔBmEcKL1-PSG strain. Significant differences were defined as ** p < 0.01, *** p < 0.001. The meaning of ns is that there is no significant difference.
Supplementary Materials
The following supporting information can be downloaded at:
References
1. Xia, Q.Y.; Zhou, Z.Y.; Lu, C.; Cheng, D.J.; Dai, F.Y.; Li, B.; Zhao, P.; Zha, X.F.; Cheng, T.C.; Chai, C.L. et al. A Draft Sequence for the Genome of the Domesticated Silkworm (Bombyx mori). Science; 2004; 306, pp. 1937-1940. [DOI: https://dx.doi.org/10.1126/science.1102210]
2. Xia, Q.Y.; Wang, J.; Zhou, Z.Y.; Li, R.Q.; Fan, W.; Cheng, D.J.; Cheng, T.C.; Qin, J.J.; Duan, J.; Xu, H.F. et al. The Genome of a Lepidopteran Model Insect, the Silkworm (Bombyx mori). Insect Biochem. Mol. Biol.; 2008; 38, pp. 1036-1045. [DOI: https://dx.doi.org/10.1016/j.ibmb.2008.11.004]
3. Xia, Q.Y.; Guo, Y.R.; Zhang, Z.; Li, D.; Xuan, Z.L.; Li, Z.; Dai, F.Y.; Li, Y.R.; Cheng, D.J.; Li, R.Q. et al. Complete Resequencing of 40 Genomes Reveals Domestication Events and Genes in Silkworm (Bombyx mori). Science; 2009; 326, pp. 433-436. [DOI: https://dx.doi.org/10.1126/science.1176620] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19713493]
4. Xia, Q.Y.; Cheng, D.J.; Duan, J.; Wang, G.H.; Cheng, T.C.; Zha, X.F.; Liu, C.; Zhao, P.; Dai, F.Y.; Zhang, Z. et al. Microarray-Based Gene Expression Profiles in Multiple Tissues of the Domesticated Silkworm, Bombyx mori. Genome Biol.; 2007; 8, 13. [DOI: https://dx.doi.org/10.1186/gb-2007-8-8-r162] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17683582]
5. Xiang, H.; Zhu, J.D.; Chen, Q.; Dai, F.Y.; Li, X.; Li, M.W.; Zhang, H.Y.; Zhang, G.J.; Li, D.; Dong, Y. et al. Single Base-Resolution Methylome of the Silkworm Reveals a Sparse Epigenomic Map. Nat. Biotechnol.; 2010; 28, pp. 516-520. [DOI: https://dx.doi.org/10.1038/nbt.1626] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20436463]
6. Daimon, T.; Kiuchi, T.; Takasu, Y. Recent Progress in Genome Engineering Techniques in the Silkworm (Bombyx mori). Dev. Growth Differ.; 2014; 56, pp. 14-25. [DOI: https://dx.doi.org/10.1111/dgd.12096] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24175911]
7. Ando, T.; Fujiwara, H. Electroporation-Mediated Somatic Transgenesis for Rapid Functional Analysis in Insects. Development; 2013; 140, pp. 454-458. [DOI: https://dx.doi.org/10.1242/dev.085241] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23250219]
8. Kawamoto, M.; Jouraku, A.; Toyoda, A.; Yokoi, K.; Minakuchi, Y.; Katsuma, S.; Fujiyama, A.; Kiuchi, T.; Yamamoto, K.; Shimada, T. High-Quality Genome Assembly of the Silkworm, Bombyx mori. Insect Biochem. Mol. Biol.; 2019; 107, pp. 53-62. [DOI: https://dx.doi.org/10.1016/j.ibmb.2019.02.002] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30802494]
9. Couble, P.; Michaille, J.J.; Garel, A.; Couble, M.L.; Prudhomme, J.C. Developmental Switches of Sericin Messenger-Rna Splicing in Individual Cells of Bombyx mori Silkgland. Dev. Biol.; 1987; 124, pp. 431-440. [DOI: https://dx.doi.org/10.1016/0012-1606(87)90496-9]
10. Chevillard, M.; Couble, P.; Prudhomme, J.C. Complete Nucleotide-Sequence of the Gene Encoding the Bombyx mori Silk Protein-P25 and Predicted Amino-Acid-Sequence of the Protein. Nucleic Acids Res.; 1986; 14, pp. 6341-6342. [DOI: https://dx.doi.org/10.1093/nar/14.15.6341]
11. Sprague, K.U. Bombyx mori Silk Proteins—Characterization of Large Polypeptides. Biochemistry; 1975; 14, pp. 925-931. [DOI: https://dx.doi.org/10.1021/bi00676a008]
12. Yamaguchi, K.; Kikuchi, Y.; Takagi, T.; Kikuchi, A.; Oyama, F.; Shimura, K.; Mizuno, S. Primary Structure of the Silk Fibroin Light Chain Determined by Cdna Sequencing and Peptide Analysis. J. Mol. Biol.; 1989; 210, pp. 127-139. [DOI: https://dx.doi.org/10.1016/0022-2836(89)90295-7]
13. Edgar, B.A.; Orr-Weaver, T.L. Endoreplication Cell Cycles: More for Less. Cell; 2001; 105, pp. 297-306. [DOI: https://dx.doi.org/10.1016/S0092-8674(01)00334-8]
14. Perdrixgillot, S. DNA-Synthesis and Endomitoses in the Giant Nuclei of the Silkgland of Bombyx mori. Biochimie; 1979; 61, pp. 171-204. [DOI: https://dx.doi.org/10.1016/S0300-9084(79)80066-8] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/465570]
15. Matsui, H.; Kakei, M.; Iwami, M.; Sakurai, S. Hormonal Regulation of the Death Commitment in Programmed Cell Death of the Silkworm Anterior Silk Glands. J. Insect Physiol.; 2012; 58, pp. 1575-1581. [DOI: https://dx.doi.org/10.1016/j.jinsphys.2012.09.012]
16. Dhawan, S.; Gopinathan, K.P. Cell Cycle Events During the Development of the Silk Glands in the Mulberry Silkworm (Bombyx mori). Dev. Genes Evol.; 2003; 213, pp. 435-444. [DOI: https://dx.doi.org/10.1007/s00427-003-0343-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12883881]
17. Sudhakar, B.; Gopinathan, K.P. Expression of Cyclin E in Endomitotic Silk-Gland Cells from Mulberry Silkworm. Gene; 2000; 257, pp. 77-85. [DOI: https://dx.doi.org/10.1016/S0378-1119(00)00391-7]
18. Kokubo, H.; Ueno, K.; Amanai, K.; Suzuki, Y. Involvement of the Bombyx mori Scr Gene in Development of the Embryonic Silk Gland. Dev. Biol.; 1997; 186, pp. 46-57. [DOI: https://dx.doi.org/10.1006/dbio.1997.8578] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9188752]
19. Zhang, C.D.; Li, F.F.; Chen, X.Y.; Huang, M.H.; Zhang, J.; Cui, H.J.; Pan, M.H.; Lu, C. DNA Replication Events During Larval Silk Gland Development in the Silkworm, Bombyx mori. J. Insect Physiol.; 2012; 58, pp. 974-978. [DOI: https://dx.doi.org/10.1016/j.jinsphys.2012.04.017]
20. Chunduri, A.R.; Rajan, R.; Lima, A.; Ramamoorthy, S.; Mamillapalli, A. Dynamics of Nuclear Matrix Attachment Regions During 5th Instar Posterior Silk Gland Development in Bombyx mori. BMC Genom.; 2022; 23, 247. [DOI: https://dx.doi.org/10.1186/s12864-022-08446-3]
21. Wang, W.; Wang, X.R.; Luo, C.Y.; Pu, Q.; Yin, Q.; Xu, L.L.; Peng, X.Y.; Ma, S.Y.; Xia, Q.Y.; Liu, S.P. Let-7 Microrna Is a Critical Regulator in Controlling the Growth and Function of Silk Gland in the Silkworm. RNA Biol.; 2020; 17, pp. 703-717. [DOI: https://dx.doi.org/10.1080/15476286.2020.1726128]
22. Fu, Y.; Wang, Y.; Huang, Q.X.; Zhao, C.Y.; Li, X.M.; Kan, Y.C.; Li, D.D. Long Noncoding Rna Lncr17454 Regulates Metamorphosis of Silkworm through Let-7 Mirna Cluster. J. Insect Sci.; 2022; 22, 12. [DOI: https://dx.doi.org/10.1093/jisesa/ieac028] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35640247]
23. Li, Q.R.; Deng, X.J.; Yang, W.Y.; Huang, Z.J.; Tettamanti, G.; Cao, Y.; Feng, Q.L. Autophagy, Apoptosis, and Ecdysis-Related Gene Expression in the Silk Gland of the Silkworm (Bombyx mori) during Metamorphosis. Can. J. Zool.; 2010; 88, pp. 1169-1178. [DOI: https://dx.doi.org/10.1139/Z10-083]
24. Yi, H.Y.; Yang, W.Y.; Wu, W.M.; Li, X.X.; Deng, X.J.; Li, Q.R.; Cao, Y.; Zhong, Y.J.; Huang, Y.D. Bmcalpains Are Involved in Autophagy and Apoptosis During Metamorphosis and after Starvation in Bombyx mori. Insect Sci.; 2018; 25, pp. 379-388. [DOI: https://dx.doi.org/10.1111/1744-7917.12417]
25. Li, G.L.; Li, Y.Z.; He, C.H.; Wei, Y.T.; Cai, K.P.; Lu, Q.Y.; Liu, X.B.; Zhu, Y.Z.; Xu, K.Z. The Promoting Effects of Pyriproxyfen on Autophagy and Apoptosis in Silk Glands of Non-Target Insect Silkworm, Bombyx mori. Pestic. Biochem. Physiol.; 2023; 196, 105586. [DOI: https://dx.doi.org/10.1016/j.pestbp.2023.105586]
26. Scanlan, J.L.; Gledhill-Smith, R.S.; Battlay, P.; Robin, C. Genomic and Transcriptomic Analyses in Drosophila Suggest That the Ecdysteroid Kinase-Like (Eckl) Gene Family Encodes the ‘Detoxification-by-Phosphorylation’ Enzymes of Insects. Insect Biochem. Mol. Biol.; 2020; 123, 103429. [DOI: https://dx.doi.org/10.1016/j.ibmb.2020.103429] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32540344]
27. Scanlan, J.L.; Battlay, P.; Robin, C. Ecdysteroid Kinase-Like (Eckl) Paralogs Confer Developmental Tolerance to Caffeine in Drosophila melanogaster. Curr. Res. Insect Sci.; 2022; 2, 100030. [DOI: https://dx.doi.org/10.1016/j.cris.2022.100030]
28. Liu, C.; Wang, J.L.; Zheng, Y.; Xiong, E.J.; Li, J.J.; Yuan, L.L.; Yu, X.Q.; Wang, Y.F. Wolbachia-Induced Paternal Defect in Drosophila Is Likely by Interaction with the Juvenile Hormone Pathway. Insect Biochem. Mol. Biol.; 2014; 49, pp. 49-58. [DOI: https://dx.doi.org/10.1016/j.ibmb.2014.03.014]
29. Nozawa, M.; Minakuchi, Y.; Satomura, K.; Kondo, S.; Toyoda, A.; Tamura, K. Shared Evolutionary Trajectories of Three Independent Neo-Sex Chromosomes in Drosophila. Genome Res.; 2021; 31, pp. 2069-2079. [DOI: https://dx.doi.org/10.1101/gr.275503.121]
30. Sonobe, H.; Ohira, T.; Ieki, K.; Maeda, S.; Ito, Y.; Ajimura, M.; Mita, K.; Matsumoto, H.; Wilder, M.N. Purification, Kinetic Characterization, and Molecular Cloning of a Novel Enzyme, Ecdysteroid 22-Kinase. J. Biol. Chem.; 2006; 281, pp. 29513-29524. [DOI: https://dx.doi.org/10.1074/jbc.M604035200]
31. Fujinaga, D.; Gu, J.J.; Kawahara, H.; Ogihara, M.H.; Kojima, I.; Takeshima, M.; Kataoka, H. Twenty-Hydroxyecdysone Produced by Dephosphorylation and Ecdysteroidogenesis Regulates Early Embryonic Development in the Silkmoth, Bombyx mori. Insect Biochem. Mol. Biol.; 2020; 127, 103491. [DOI: https://dx.doi.org/10.1016/j.ibmb.2020.103491] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33096212]
32. Wang, C.F.; Zhang, Z.; Sun, W. Ecdysone Oxidase and 3-Dehydroecdysone-3β-Reductase Contribute to the Synthesis of Ecdysone During Early Embryonic Development of the Silkworm. Int. J. Biol. Sci.; 2018; 14, pp. 1472-1482. [DOI: https://dx.doi.org/10.7150/ijbs.26227] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30262999]
33. Liu, R.P.; Zeng, W.H.; Tan, T.T.; Chen, T.; Luo, Q.; Qu, D.W.; Tang, Y.Y.; Long, D.P.; Xu, H.F. Insights into the Regulatory Characteristics of Silkworm Fibroin Gene Promoters Using a Modified Gal4/Uas System. Transgenic Res.; 2019; 28, pp. 627-636. [DOI: https://dx.doi.org/10.1007/s11248-019-00175-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31564030]
34. Schoenfelder, S.; Fraser, P. Long-Range Enhancer-Promoter Contacts in Gene Expression Control. Nat. Rev. Genet.; 2019; 20, pp. 437-455. [DOI: https://dx.doi.org/10.1038/s41576-019-0128-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31086298]
35. Bellés, X. When Inordinate Tissue Growth Is Beneficial: Improving Silk Production by Increasing Silk Gland Size. Cell Res.; 2011; 21, pp. 862-863. [DOI: https://dx.doi.org/10.1038/cr.2011.68] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/21633374]
36. Ma, L.; Xu, H.F.; Zhu, J.Q.; Ma, S.Y.; Liu, Y.; Jiang, R.J.; Xia, Q.Y.; Li, S. Ras1 (Ca) Overexpression in the Posterior Silk Gland Improves Silk Yield. Cell Res.; 2011; 21, pp. 934-943. [DOI: https://dx.doi.org/10.1038/cr.2011.36]
37. Zhang, P.L.; Liu, S.M.; Song, H.S.; Zhang, G.Z.; Jia, Q.Q.; Li, S. Yorkie (Ca) Overexpression in the Posterior Silk Gland Improves Silk Yield in Bombyx mori. J. Insect Physiol.; 2017; 100, pp. 93-99. [DOI: https://dx.doi.org/10.1016/j.jinsphys.2017.06.001]
38. Brenner, S. Phosphotransferase Sequence Homology. Nature; 1987; 329, 21. [DOI: https://dx.doi.org/10.1038/329021a0]
39. Bairoch, A.; Claverie, J.M. Sequence Patterns in Protein-Kinases. Nature; 1988; 331, 22. [DOI: https://dx.doi.org/10.1038/331022a0]
40. Hanks, S.K.; Hunter, T. Protein Kinases. 6. The Eukaryotic Protein-Kinase Superfamily—Kinase (Catalytic) Domain-Structure and Classification. FASEB J.; 1995; 9, pp. 576-596. [DOI: https://dx.doi.org/10.1096/fasebj.9.8.7768349]
41. Peisach, D.; Gee, P.; Kent, C.; Xu, Z.H. The Crystal Structure of Choline Kinase Reveals a Eukaryotic Protein Kinase Fold. Structure; 2003; 11, pp. 703-713. [DOI: https://dx.doi.org/10.1016/S0969-2126(03)00094-7] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/12791258]
42. Yuan, C.; Kent, C. Identification of Critical Residues of Choline Kinase A2 from Caenorhabditis Elegans. J. Biol. Chem.; 2004; 279, pp. 17801-17809. [DOI: https://dx.doi.org/10.1074/jbc.M401382200] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/14960577]
43. Burk, D.L.; Hon, W.C.; Leung, A.K.W.; Berghuis, A.M. Structural Analyses of Nucleotide Binding to an Aminoglycoside Phosphotransferase. Biochemistry; 2001; 40, pp. 8756-8764. [DOI: https://dx.doi.org/10.1021/bi010504p] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11467935]
44. Beydon, P.; Girault, J.P.; Lafont, R. Ecdysone Metabolism in Pieris-Brassicae During the Feeding Last Larval Instar. Arch. Insect Biochem. Physiol.; 1987; 4, pp. 139-149. [DOI: https://dx.doi.org/10.1002/arch.940040207]
45. Hou, S.H.; Sun, Y.; Wu, Y.C.; Cheng, T.C.; Liu, C. Bmsage Is Involved in the Determination of Cell Number in the Silk Gland of Bombyx mori. Insect Biochem. Mol. Biol.; 2019; 113, 103205. [DOI: https://dx.doi.org/10.1016/j.ibmb.2019.103205] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31421207]
46. Kato, N.; Sato, S.; Yamanaka, A.; Yamada, H.; Fuwa, N.; Nomura, M. Silk Protein, Sericin, Inhibits Lipid Peroxidation and Tyrosinase Activity. Biosci. Biotechnol. Biochem.; 1998; 62, pp. 145-147. [DOI: https://dx.doi.org/10.1271/bbb.62.145] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9501526]
47. Zhang, Y.Q. Applications of Natural Silk Protein Sericin in Biomaterials. Biotechnol. Adv.; 2002; 20, pp. 91-100. [DOI: https://dx.doi.org/10.1016/S0734-9750(02)00003-4]
48. Inoue, S.; Tanaka, K.; Arisaka, F.; Kimura, S.; Ohtomo, K.; Mizuno, S. Silk Fibroin of Bombyx mori Is Secreted, Assembling a High Molecular Mass Elementary Unit Consisting of H-Chain, L-Chain, and P25, with a 6:6:1 Molar Ratio. J. Biol. Chem.; 2000; 275, pp. 40517-40528. [DOI: https://dx.doi.org/10.1074/jbc.M006897200]
49. Wang, Y.Q.; Li, Z.Q.; Xu, J.; Zeng, B.S.; Ling, L.; You, L.; Chen, Y.Z.; Huang, Y.P.; Tan, A.J. The Crispr/Cas System Mediates Efficient Genome Engineering in Bombyx mori. Cell Res.; 2013; 23, pp. 1414-1416. [DOI: https://dx.doi.org/10.1038/cr.2013.146]
50. Takasu, Y.; Kobayashi, I.; Beumer, K.; Uchino, K.; Sezutsu, H.; Sajwan, S.; Carroll, D.; Tamura, T.; Zurovec, M. Targeted Mutagenesis in the Silkworm (Bombyx mori) Using Zinc Finger Nuclease Mrna Injection. Insect Biochem. Mol. Biol.; 2010; 40, pp. 759-765. [DOI: https://dx.doi.org/10.1016/j.ibmb.2010.07.012]
51. Gudmunds, E.; Wheat, C.W.; Khila, A.; Husby, A. Functional Genomic Tools for Emerging Model Species. Trends Ecol. Evol.; 2022; 37, pp. 1104-1115. [DOI: https://dx.doi.org/10.1016/j.tree.2022.07.004] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35914975]
52. Expasy. Available online: https://www.expasy.org/ (accessed on 10 February 2023).
53. Interpro. Available online: http://www.ebi.ac.uk/interpro/ (accessed on 15 February 2023).
54. Smart. Available online: http://smart.embl-heidelberg.de/ (accessed on 21 February 2023).
55. Wang, G.H.; Xia, Q.Y.; Cheng, D.J.; Duan, J.; Zhao, P.; Chen, J.; Zhu, L. Reference Genes Identified in the Silkworm Bombyx mori During Metamorphism Based on Oligonucleotide Microarray and Confirmed by Qrt-Pcr. Insect Sci.; 2008; 15, pp. 405-413. [DOI: https://dx.doi.org/10.1111/j.1744-7917.2008.00227.x]
56. Cctop. Available online: https://www.cos.uni-heidelberg.de/en (accessed on 16 May 2022).
57. Yusa, K. Piggybac Transposon. Microbiol. Spectr.; 2015; 3, MDNA3-0028-2014. [DOI: https://dx.doi.org/10.1128/microbiolspec.MDNA3-0028-2014]
58. Heid, C.A.; Stevens, J.; Livak, K.J.; Williams, P.M. Real Time Quantitative Pcr. Genome Res.; 1996; 6, pp. 986-994. [DOI: https://dx.doi.org/10.1101/gr.6.10.986] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/8908518]
59. Winer, J.; Jung, C.K.S.; Shackel, I.; Williams, P.M. Development and Validation of Real-Time Quantitative Reverse Transcriptase-Polymerase Chain Reaction for Monitoring Gene Expression in Cardiac Myocytes In Vitro. Anal. Biochem.; 1999; 270, pp. 41-49. [DOI: https://dx.doi.org/10.1006/abio.1999.4085] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10328763]
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Details
; Shen, Guanwang 2 ; Zhao, Ping 2 ; Xia, Qingyou 2
1 Integrative Science Center of Germplasm Creation in Western China (CHONGQING) Science City, Biological Science Research Center, Southwest University, Chongqing 400716, China;
2 Integrative Science Center of Germplasm Creation in Western China (CHONGQING) Science City, Biological Science Research Center, Southwest University, Chongqing 400716, China;