Advances in Editing Silkworms (Bombyx mori) Genome by Using the CRISPR-Cas System.

CRISPR (clustered regularly interspaced short palindromic repeats)-Cas (CRISPR-associated) represents a powerful genome editing technology that revolutionized in a short period of time numerous natural sciences branches. Therefore, extraordinary progress was made in various fields, such as entomology or biotechnology. Bombyx mori is one of the most important insects, not only for the sericulture industry, but for numerous scientific areas. The silkworms play a key role as a model organism, but also as a bioreactor for the recombinant protein production. Nowadays, the CRISPR-Cas genome editing system is frequently used in order to perform gene analyses, to increase the resistance against certain pathogens or as an imaging tool in B. mori. Here, we provide an overview of various studies that made use of CRISPR-Cas for B. mori genome editing, with a focus on emphasizing the high applicability of this system in entomology and biological sciences.

1. Introduction

The life sciences research fields were revolutionized by the outstanding development of various genome editing tools. By using specific techniques of genome editing, the genomic DNA of every living organism can be submitted to guided changes, such as deletions, insertions, and sequence substitutions [1].

In recent years, several genome editing tools have been in the spotlight. Among them, there are three remarkable technologies, namely those relying on programmable nucleases (i.e., the transcription activator like effector nucleases (TALENs)), zinc finger nucleases (ZFNs), and clustered regularly interspaced short palindromic repeat - associated nucleases (CRISPR-Cas) [2,3]. Currently, by using engineered nucleases, remarkable advances are being made regarding the correction of genetic mutations, gene expression regulation, and the development of therapeutic agents; these approaches are also used for a better understanding of gene functions and the mechanisms underlying the development of certain genetic disorders or various diseases [4].

When it comes to genome editing, it is crucial to avoid off-target effects, but overall, the CRISPR-Cas system exhibits reliable results, owing to a great degree of fidelity [4,5]. Since its discovery in bacteria, the CRISPR-Cas system has been continuously exploited, representing an extremely versatile tool for the scientific community due to its reprogrammable feature. Currently, this system is used to edit the genomes of various organisms, such as bacteria, insects, plants, or human cells [6].

For more than 5000 years, humans have been involved in the domestication of one of the most economically important insects, namely Bombyx mori. Historically, B. mori gained prominence due to silk production, but now it is one of the most valued model organisms for life science branches. Just to mention a few recent achievements, B. mori has been used as a bioreactor, for recombinant proteins production [7,8,9,10], while its silk is used to produce extraordinary silk-based biomaterials that exhibit great importance for the medical field [11]. The silk proteome contains two major proteins, namely fibroin and sericin. The silk gland consists of three anatomical areas, namely the posterior (PSG), the middle (MSG), and the anterior (ASG) silk gland. Fibroin is synthesized in the PSG and sericin is assembled in the MSG [12].

In 2018, Ma et al. [13] excellently summarized the advances of using genome editing tools in silkworms. However, herein, we mainly focus on the potential of CRISPR-Cas technology in editing B. mori’s genome [13]. We navigate through the recent progresses in using the outstanding CRISPR-Cas system in B. mori and discuss the latest studies that utilized this approach in order to investigate the genes function, to regulate the gene expression and to enhance the resistance against B. mori nucleopolyhedrovirus (BmNPV). Additionally, the key role of B. mori in the scientific fields will be discussed.

2. The CRISPR-Cas System

CRISPR-Cas is one of the key methods employed by many molecular biology scientific laboratories. Since its first description [14], genome editing focused research was implemented by countless research groups [15,16,17].

2.1. The CRISPR-Cas Complex Role in the Immunity System

When investigating the iap gene product in the opportunistic pathogen Escherichia coli, Ishino et al. (1987) [18] observed an atypical structure, specifically the repetition of several homologous sequences. Later, this type of structure was observed in various bacterial, as well as archaeal strains [19,20]. Subsequently, these repetitive sequences were linked with exogenous genetic material, and following several years, their assembling mechanism and function were elucidated [21]. This type of sequence can be placed on the chromosomal DNA, but it can also be found on the plasmid DNA [22].

The scientists demonstrated that CRISPR-Cas, which is present in one-third of bacteria and nearly in all archaea, has a key role in host’s adaptive immunity. It protects the organism against various intruders, such as viruses, but it also offers protection against other mobile genetic elements, such as transposons or plasmids [23].

The CRISPR-Cas system structure includes three main components, i.e., the CRISPR arrays, the associated Cas proteins, and the leader nucleotide sequence. The first genetic component, the CRISPR locus, is characterized by identical repeats structures (21–37 bp) that are highly conserved and the spacer sequences that are acquired fragments of invader’s nucleic acid material. The CRISPR array is located downstream from cas genes. The latter encodes for Cas proteins that are crucial to the immune reaction [24].

Initially, only four distinct Cas proteins (1–4) were reported, due to the rapid evolution of biological sciences; currently, there numerous Cas proteins have been described [25,26], Cas1 being the most analyzed [27].

CRISPR-Cas has a great adaptability, with host-related specificities; thus, it exhibits a significant diversity. The varying feature is defined by the CRISPR array and the cas gene sequences. The classification of these types of systems is based on the signature Cas proteins. Currently, there are two major classes of CRISPR-Cas systems, each also di-vided in several groups [28]. Regarding the leader nucleotide sequence, it has been shown that it has a key role by carrying the essential promoter sequences for the transcription of CRISPR loci. Besides the promoter, the leader contains specific signals that are crucial for the adaptation stage from the first phase of CRISPR-Cas activation [29].

The adaptation is the first functional stage of the CRISPR-Cas mechanism, during which the foreign nucleic acid is recognized by several Cas proteins [30] and consequently integrated next to a leader sequence. Through this mechanism, in evolution, the spacers are arranged chronologically, and this feature helps bacteria and archaea to enhance their protection against the genetic material of the latest foreign encounter [31]. Each new acquired spacer is accompanied by a repeat sequence; therefore, the CRISPR array expands with every invasion [32].

The CRISPR array is transcribed in the second step, specifically in the biogenesis phase [33]. First of all, it is being transcribed into a precursor CRISPR RNA (crRNA). At the end of this phase, there are numerous mature crRNAs molecules, resulting from the action of RNase III that process the precursor crRNA. Each crRNA includes a spacer and a repeat sequence [31,34,35].

The last step of the system’s mechanism is the interference phase. It involves the degradation of the foreign nucleic acid, by targeting and cleaving it [36]. The products of the biogenesis phase, the crRNAs, act like guides for targeting the invader, which is then cleaved following a Cas proteins cascade that act like molecular scissors [37].

2.2. The CRISPR-Cas System as a Genome Editing Tool

When it comes to leading tools in genetic engineering, the CRISPR-Cas system can be considered the foremost instrument. After elucidating its function in various organisms, scientists aimed to exploit its versatility, in order to overcome the disadvantages of other available genome editing tools [38]. Even if scientific studies still report the use of ZFNs and TALENs techniques as editing tools, the CRISPR-Cas system is the most effective genome editing instrument, standing on top in regards to efficiency, cost-effectiveness, and the relative simplicity of use [39] (Table 1). Another considerable advantage of this system is represented by its capacity to simultaneously target multiple genes [40].

Table 1

Comparison between TALEN, ZFN, and CRISPR-Cas gene editing technologies.

Traits	TALEN	ZFN	CRISPR-Cas	References
Origin	Prokaryotic	Eukaryotic	Prokaryotic	[41]
Efficiency (%)	76	12	81	[42]
Specificity	Moderate	Low	High	[43,44,45,46]
Target site recognition	Any site	Any site	Pam motif (NGG) required	[43]
Multiplex potential	Low	Low	High	[43,45]
Processing time	Time consuming	Time consuming	Short	[45]
Methylation sensitive	Sensitive	Sensitive	Not sensitive	[42]
Engineering feasibility	Moderate/High	Moderate	Moderate/High	[42,45]
Dimerization required	Yes	Yes	No	[44]
Cost effectiveness	Moderate	No	Yes	[43,44]

Of the numerous CRISPR-Cas systems, CRISPR-Cas9 is currently the most used instrument in laboratories across the world [47]. The Cas9 nuclease is the signature protein of CRISPR-Cas II systems and it is responsible for double strand DNA breaks [27]. Three different methods to deliver the Cas9 endonuclease have been described. It can be directly delivered by microinjection into the embryos, while the other two delivery methods involve a plasmid that expresses the Cas9 enzyme, or a messenger RNA (mRNA) sequence that encodes it. Of the three techniques, in terms of genome engineering, the earliest mentioned is the best option due to its certain advantages. By directly delivering the protein, low immunogen effects were observed. Furthermore, the off-target activity is minimized compared with the other two methods [48]. The use of CRISPR-Cas9 is a simple but powerful genome editing tool, with various implementations, and their impact on new research trends has been reviewed elsewhere [49].

The CRISPR-Cas9 mechanism relies on the Cas9 nuclease and a guide sequence (gRNA). As the name implies, the gRNA has the role to guide the Cas9 nuclease to a target site in order to cleave the DNA. The key feature of gRNA is the extensive complementarity with the target sequence [50].The protospacer-adjacent motif (PAM) bordering the target complementary sequence has a key role, since in its absence, the CRISPR-Cas systems would degrade their own CRISPR loci. In order to perform a cleavage, the Cas9 protein scans for the PAM sequence. Even if the gRNA is complementary with the target sequence, the Cas9 endonuclease will not cleave it in the absence of PAM [51].

The central factor that influences the success of the gene editing process is the repair path of the double-strand breaks produced by Cas9. There are two main repair pathways: the homology-directed repair (HDR) and the nonhomologous end joining (NHEJ), respectively [52]. More often, NHEJ is exploited in order to acquire indels mutations, specially to obtain small deletions. These deletions are extremely useful for disclosing gene functions [53]. However, the HDR machinery is used not just to obtain knock-out or knock-down mutations, such as the expected output following NHEJ, but to generate target knock-ins. Therefore, by using HDR, exogenous sequences can be successfully integrated into the host’s genome. Currently, major efforts are being made in order to enhance the sequence replacement by using the HDR mechanism [54].

CRISPR-Cas9 is currently used in multiple research fields, such as agriculture (editing of various agricultural plant genomes or pest insect’s genome) [55,56,57], biotechnology, food industry, and medicine (modeling diseases using HeLa cells, deciphering HIV infection mechanisms, using various experimental models, such as Danio rerio to tackle cancer and neurological diseases, etc.) [58,59], just to mention a few (Figure 1) [60].

Figure 1

Schematic representation of the most important current applications of CRISPR-Cas9 in entomology, medicine, and agriculture. On top, a simplified description of CRISPR-Cas9 applicability in B. mori that is extensively described in the main text (created with BioRender.com, accessed on 2 December 2021).

2.3. CRISPR-Cas9 in Entomology

Being the most diverse and numerous category of organisms for decades [61], insects have been intensively studied. Countless studies have been performed due to insects’ key roles in ecology, agriculture, and medicine [62,63]. Considering this, numerous research groups aimed to use the CRISPR-Cas9 system to manipulate the insects’ genome. The first application of CRISPR-Cas9 was performed in Drosophila melanogaster [64] due to its strategic importance as arguably the main experimental model organism for life sciences [65]. Besides D. melanogaster, the researchers also used the CRISPR-Cas9 applicability on B. mori, Apis mellifera, Aedes aegypti, and Tribolium castaneum [66,67].

Gratz et al. (2013) [68] programmed CRISPR-Cas9 to edit Drosophila’s genome. The authors targeted the yellow gene, which is commonly used in various studies. First, they aimed to determine if this genome engineering tool could be efficient and could fulfill its role to induce breaks in the target sequence. By using the CRISPR-Cas9 system in Drosophila, not only was the yellow gene successfully knocked out, but the genome’s alterations were also germline transmitted. Subsequent to the deletion of the target gene, a donor sequence was designed. This sequence provided the template for the HDR repair pathway and its use was to test the accuracy of specific replacement of yellow gene with an exogenous sequence. These sequence replacements were transmitted to descendants as well. Their data showed that there was no off-target activity and it highlighted the feasibility of using the CRISPR-Cas9 technology in eukaryotes [68].

Aiming to further highlight the feasibility of choosing this system to perform genome alteration in Drosophila, Yu et al. (2013) [69] designed two gRNAs to induce mutations in two regions of the yellow gene. In addition, they targeted other six sequences, both euchromatic and heterochromatic loci. Remarkably, a definite mutation in ms(3)k81 was transmitted to descendants in a proportion of 100%. By successfully targeting heterochromatic loci, their result showed that the CRISPR-Cas9 system is efficient for altering the heterochromatin [69]. Drosophila have been used in numerous studies in order to examine the insecticide resistance [70,71,72]. In this direction, Douris et al. (2020) [56] notably summarized the progress in using CRISPR-Cas9 to explore the genetic basis of this mechanism.

The CRISPR-Cas9 technique was used to perform functional analysis concomitantly on two genes belonging to the cricket (Gryllus bimaculatus) [73]. G. bimaculatus is an important insect for experimental studies; for example, it plays an important role for evolutionary developmental studies and comparative biology, but it is also a relevant model organism for neurobiology and behavioral sciences [74]. The efficiency of inserting a donor sequence via a homology-independent technique was tested in two hox genes, namely Gb-Ubx and Gb-abd-A. After inserting the donor fragment into essential exons of both genes, their function was lost. Thus, functional investigations of hox genes could be carried out by using the knock-in/knock-out approaches [73].

Being one of the most important social insects [75] and as it plays a crucial role as a pollinator, the honeybee (Apis mellifera) has been intensively studied. It also plays a pivotal role in various therapeutic areas due to honey production. This natural product has extraordinary benefits for human health, exhibiting antioxidant, antiviral, and antibacterial effects [76]. Due to its special characteristics, the use of honey is not limited to humans, but this natural product is being used to improve certain features of other insects, such as silkworms [77]. There are numerous studies that detail functional analysis of A. mellifera genes by exploiting the CRISPR-Cas9 system. For instance, Hu et al. (2019) [67] reported the successful utilization of this system for knocking out the mrjp 1 gene from the honeybee genome. The CRISPR-Cas9 complex was delivered through microinjection and they tested two specific regions of embryos, for identifying the most convenient structure for delivering the gRNA and the Cas9 endonuclease. By microinjecting the construct at the dorsal posterior side, there was a low rate of successful manipulation (11.8%); however, when choosing the ventral cephalic side, the results showed a great rate of gene editing (93.3%). Trying to validate the previous results, the authors also targeted pax6. Based on the previously obtained results, they microinjected the CRISPR-Cas9 construct at the ventral cephalic side. The results showed an editing rate of 100% [67]. Targeting the same gene, mrjp 1, similar results have been obtained in another study [78]. Thus, functional analysis of A. mellifera genes can be effectively performed by using the CRISPR-Cas9 system.

Considering the same topic of gene function research, Nie et al. (2021) [79] used the CRISPR-Cas9 technology to determine if the yellow-y gene plays a crucial role in the process of cuticular melanin synthesis in A. mellifera. They targeted this gene due to its great potential for mutants screening, being a selectable marker. By disrupting it, the phenotype of worker cuticle has changed, mainly due to the black pigment decreasing, thus confirming the yellow-y gene critical role in melanin pigmentation. However, as future prospects, this could be a great genetic marker for upcoming genomic research [79].

Referring to A. mellifera sex determination, it is controlled by the heterozygosity at a particular locus that harbors the key complementary sex determiner (csd) gene. The bees that are heterozygous at this specific locus are females, while the males are homozygous or hemizygous [80]. In a recent study, Wang et al. (2021) [81] used the CRISPR-Cas9 tool in order to knock out the csd gene and, thus, eliminated the genetic difference between females and males. Subsequently, they aimed to observe the transcriptome difference between the two sexes in this particular genetic background. They also successfully induced target mutations in mutant haploid individuals. It was observed that the expression level of several male-biased genes was higher in the mutant males. On the other hand, the expression level of several specific female-biased genes was lower. Their data also confirmed that csd interacts with certain genes, such as fruitless, troponin T, and transformer-2 just to mention a few [81].

3. Bombyx mori

For numerous reasons, B. mori is one of the most studied insects, especially because it presents a real interest for the scientific community. It has been completely domesticated and plays the pivotal role in sericulture, being reared principally for silk production on a large scale [82]. Considering the great role of silk in the textile industry and its use as a biomaterial in medicine, major efforts are being made to enhance silk’s quality, but also to increase its quantity [12]. Moreover, there are numerous studies that describe the process of obtaining enhanced silk with appreciable properties by genetically manipulating the silkworms [83,84,85,86]. B. mori’s genome was manipulated in order to enhance silk’s properties. For instance, by knocking in the major ampullate silk protein from spiders into B. mori’s genome, a research group obtained silk with superior mechanical characteristics [84].

B. mori is an oligophagous insect and its main food source are the mulberry leaves, a nutritional preference that influences its biological and economical parameters. There is a major drawback when it comes to this source of nutrition, namely the limited availability of mulberry leaves. In order to be able to rear silkworms not only in spring and summer, artificial diets are currently being used for their nourishment [77,87].

Regardless of human medicine advances, certain microorganisms developed survival strategies and numerous infectious diseases remain a crucial problem worldwide [88]. In this context, B. mori is a reliable experimental model and exhibits the great advantage of a short development cycle that consists of four different stages. The first period of growth is the egg phase, followed by the larvae, pupa, and moth phases. The larvae phase plays a crucial role in silkworm’s development, being the longest phase and including five different stages. Another significant advantage of using B. mori is the fact that in the last instar larvae, the body size has nearly 5 cm; hence, it is facile to manipulate or to exploit it for various purposes. Moreover, this size also facilitates the body’s dissection; thus, multiple target tissues or organs can easily be obtained [89,90].

B. mori as a Model Organism

The level of attention received as a model experimental organism increased, since several groups made available B. mori’s genome data [91]. Various mutant strains have been described and the genetic analysis confirmed the numerous genetic traits. Another important aspect is that its manipulation is not associated with ethical concerns [90].

Being susceptible to various infectious agents, B. mori is one of the most used experimental models for drug screening, evaluation of different virulence factors, and identifying the pathogen genes responsible for its virulence [90]. Hitherto, several studies used B. mori to evaluate the effectiveness of antibiotics against certain human pathogens [92,93,94]. In a recent study, silkworms were used to examine the efficacy of three different glycopeptide antibiotics against Staphylococcus aureus infection. As a result, authors highlighted the great feasibility and efficacy of using B. mori to mimic bacterial infections in order to examine the therapeutic potential of antibiotics [95]. Moreover, recently, silkworms were used to evaluate the impact on several antibacterial compounds against Cutibacterium acnes [96]. As an experimental model organism, B. mori is currently being intensively used for various other purposes, as detailed in Table 2.

Table 2

Summary of various uses of B. mori as an experimental model organism.

Type of Model Organism	Brief Description	Purpose	References
Human disease model	Transgenic B. mori expressing hIR (human insulin receptor)	Drug evaluation for diabetes treatment	[97,98]
	Bacterial (Listeria monocytogenes) infection model in B. mori	Evaluating the interaction between host and pathogen; investigating the activity of vitamin A against microbial infections	[99]
	Fungal (C. albicans) infection model in silkworms	Assessing the C. albicans biofilm development	[100]
	Inducing deletions in the BmSpr gene, leading to sepiapterin reductase deficiency	Treatment options evaluations against sepiapterin reductase deficiency	[101]
Model for pesticide toxicity	Exposing silkworms to phoxim	Identifying specific biomarkers for phoxim stress; evaluating the toxicity reaction and the pretreatment with nanoparticulate titanium dioxide	[102]
	Inducing genotoxicity by feeding the silkworms with different doses of avermectin	Exploration of certain genes that are required for the DNA repairing mechanism	[103]
	Treating B. mori larvae with Fenvalerate-20EC	Evaluation of Fenvalerate-20EC impact on several digestive enzymes	[104]
Model for drugs toxicity	Injecting the silkworms with three different pharmacologically active agents (4-methyl umbelliferone, 7-ethoxycoumarine)	Evaluation of the metabolic pathway of these compounds	[92]
	Exposing the silkworms to fungal infections	Exploring pharmacokinetic parameters of an antifungal agent, Voriconazole	[105]
	Injecting cytotoxic drugs into B. mori larvae	Evaluation of cytotoxic drugs impact	[106]
Model for nanomaterials toxicity	Spreading silver nanoparticles on mulberry leaves	Toxicity evaluation of silver nanoparticles	[107,108,109,110]
	Injecting subcutaneously zinc oxide nanoparticles	Evaluation of zinc oxide nanoparticles toxicity, accumulation, and distribution	[111]
	Injecting in the dorsal vein different nanoparticles with great interest in various life science branches	Investigation of different silicon and carbon nanomaterials toxicity level against hemocytes	[112]

4. Applications of CRISPR-Cas in B. mori

The first communication of successful manipulation of B. mori’s genome by using the CRISPR-Cas9 tool, was reported by Wang et al. (2013) [113]. The authors targeted an essential gene [113], BmBlos2, that is orthologous to the Blos2 human gene [114]. Two sgRNAs (23-bp) were designed to induce mutations leading to loss of target gene function. Each complex formed by one sgRNA and the Cas9 nuclease was injected in the preblastoderm embryonic stage. Ordinarily, the larval integument is opaque, but when the BmBlos2 gene function is lost, the tegument becomes translucent. This effect could be viewed as a phenotypic marker for mutant’s detection. Of all individuals, 94% and 95.6%, respectively were successfully edited by using the two sgRNAs. This study highlighted the feasibility of using CRISPR-Cas9 not only in B. mori, but also in other lepidopteran insects. These findings are of great interest by revealing CRISPR-Cas9 system’s applicability in pest control approaches [113].

The multiplexable potential of CRISPR-Cas9 technology was highlighted by Liu et al. (2014) [115]. First, the BmBlos2 gene was targeted for site specific mutagenesis, in order to confirm the feasibility of using the CRISPR-Cas9 system in B. mori. Following this, other six genes were targeted to confirm the system’s multiplexable feature: tyrosine hydroxylase, red egg, yellow-e, kynureninase, ebony, and flugellos. Mutations were induced in each target gene, without evidence for the system’s off-target activity [115]. The multiplexable feature of CRISPR-Cas9 is facilitating the process of genome engineering by simultaneously and precisely inducing mutations in different sites. This property of CRISPR-Cas9 enables researchers to perform precise and elaborate target mutagenesis in a time-effective manner [116].

Other approaches aimed to knock out the BmKu70 gene by using CRISPR-Cas9 in order to target its second exon [117]. The BmKu70 gene is coding a highly conserved protein, Ku70, which plays a key role in numerous mechanisms: cell adhesion, apoptosis, the maintenance of telomeres length, etc. In addition, numerous studies reported that by inactivating the BmKu70 gene, the frequency of homologous repair is increased. In order to test this hypothesis, the authors performed a transient analysis in genetically manipulated embryos. They knocked in the Bm702 gene that is found on the Z chromosome. Their results confirmed that by knocking out BmKu70 the homologous repair frequency is expanded. These promising results are of great interest for fundamental research in B. mori. [117].

It has been demonstrated by Fujinaga et al. (2017) [118] that in B. mori, the insulin-like growth factor-like peptide (IGFLP), is closely related to the genital disc, particularly involved in its growth. Due to the lack of studies performed in vivo on this topic, the same research group aimed to confirm the role of IGFLP by inactivating it [119]. For this purpose, they used the CRISPR-Cas9 genome editing tool. The absence of IGFLP leads to smaller ovaries and a lower number of laid eggs compared with the wild type. On the other hand, the size of laid eggs and its development were not affected. These findings indicated that this hormone has no impact on B. mori’s fertility. The ecdysteroids play a crucial role in IGFLPs production, by inducing gene expression. Furthermore, it has been shown previously that ecdysteroids have a key role in ovary development. Therefore, the authors initially appraised that a low ecdysteroids titre caused reduced ovary weights. However, by analyzing the transgenic females, their data showed that the ecdysteroids titre was the same. This study reveals insights on IGFLPs impact on ovarian development [119].

In a recent study [120], a research group explored the effect of activating the BmFibH gene in the B. mori embryonic cells. For this purpose, they constructed a complex that involved the inactive form of Cas9 nuclease (dCas9) and a VPR activating domain, driven by a certain promoter. This activating domain consists of several different activators: VP64, p65, and Rta, respectively. When it comes to sgRNAs, three specific constructs were designed to target the promoter of the gene of interest. In order to confirm the success of BmFibH activation, first they determined its expression in normal embryonic cells. Their results showed that the BmFibH gene is strongly downregulated in untransformed cells. On the other hand, the transfected cells exhibited a higher BmFibH expression level. Moreover, their data showed that the activation of the target gene impacted the cellular stress responses [120]. Cui et al. (2018) [121] targeted the BmFibH gene in order to explore its role in the development of the silk gland. After the CRISPR-Cas9 construct was designed, a total number of 630 eggs were microinjected, but only 12.5% hatched. After analyzing the unhatched eggs, they observed that all embryos were genetically edited. By knocking out the BmFibH, severe changes were observed, such as naked pupa or thin cocoons. All individuals that exhibited naked pupae died. Moreover, by inactivating this gene, several other genes involved in the processes of degradation, such as the autophagy, were upregulated. These findings offer a better understanding of FibH protein’s role in silk gland development [121].

Keeping in mind the feasibility of expressing spider silk genes in B. mori [84], in order to obtain enhanced silk, Zhang et al. (2019) [122] used CRISPR-Cas9 to acquire high-performance fibers. By using this technique, the authors successfully knocked in spider silk genes in silkworms’ genome. Accordingly, they designed two types of systems, FibL-CRISPR-Cas9 and FibH-CRISPR-Cas9. To avoid disrupting the protein production, spider silk genes were knocked in into one of the introns of BmFibL or BmFibH. They demonstrated the feasibility of employing CRISPR-Cas9 in B. mori to obtain silk with enhanced mechanical properties at industrial scale. The strategy described in this study can be further used to obtain numerous exogenous proteins that exhibit great interest for medical applications and beyond [122].

The microRNAs (miRNAs) are key regulators when it comes to gene expression. Notably, they recognize by complementarity specific mRNAs and inactivate them [123]. It has been revealed that miR-2 is one of the most important miRNAs for wing morphogenesis in D. melanogaster and exhibits great influence on the Notch pathway. The BmAwd and BmFng genes are known positive regulators of this signaling pathway and are also potential miR-2 targets. On this topic, Ling et al. (2015) [124] used the CRISPR-Cas9 technology to investigate the function of miR-2 in B. mori. In the first phase of the study, the authors used the Gal4/UAS system to overexpress the miR-2, resulting in deformed adults’ wings. However, in the second phase of the research, the CRISPR-Cas9 system was used to knock out the two miR-2 target genes. The loss of function of BmFng and BmAwd also led to deformed wings. Both phases of the study confirmed that in silkworms, the miR-2 plays a crucial role in wing development [124].

In another study, Liu et al. (2020) [125] used the CRISPR-Cas9 system to explore the function of miR-34, another miRNA that exhibits a great impact on insect development. Firstly, they overexpressed miR-34 in a transgenic line constructed by using a pBac plasmid. The miR-34 overexpression negatively impacted the body size and wing morphology. Secondly, the CRISPR-Cas9 system was used to inactivate the miR-34 by using two different gRNAs. The miR-34 ablation led to larvae developmental delay. Using several bioinformatic tools, they predicted that the miR-34 target genes could be BmE74, BmCpg4, BmLcp, BmWcp11, and BmBrc-z2. Various analyses revealed that miR-34 target genes are just BmE74 and BmCpg4. While it is well known that BmE74 plays a key role in growth and morphogenesis, functional analysis of the second gene had to be performed. Given this, the CRISPR-Cas9 system was used to knock out BmCpg4 resulting in affected wings, thus highlighting the gene’s role in wing development [125].

Regarding the silkworms’ innate defense mechanisms, there are two major activation pathways involved in the expression of numerous evolutionarily conserved antimicrobial peptides (AMPs) genes, and the Toll and the Imd pathways [126,127]. When Gram-negative bacterial or fungal contamination occurs, the Toll pathway is activated. One of the key genes that is involved in the Toll pathway is BmCactus. Being a negative regulator, once the infection with the mentioned pathogens occurs, this gene is being phosphorylated and inactivated. Considering this, Park et al. (2018) [128] used the CRISPR-Cas9 technology to perform site-target mutagenesis targeting the BmCactus gene in a specific B. mori cell line. The authors designed six different gRNAs and transfected the CRISPR-Cas9 complex in the B. mori ovarian cell line by electroporation, but observed a very low survival rate of only 24%. Their data showed that all gRNAs determined site induced mutagenesis. By disrupting the BmCactus gene, the expression of several antimicrobial proteins (e.g., lysozyme and lebocin) was stimulated [128]. Keeping in mind the great importance of AMPs in clinical research [129], this study underlined the outstanding potential of B. mori usage in life science fields and the feasibility of using the advanced CRISPR-Cas9 genetic scissors to edit genomes for different purposes [128].

Ecdysteroids are steroid hormones that play a crucial role in the process of molting and metamorphosis in insects. The most important molting hormone is 20-hydroxyecdysone (20E) [130], but even if its biosynthesis processes have been intensively studied, the 20E metabolism has not been well-documented. However, there are several genes that are believed to be involved in the inactivation of 20E, but their biological functions have not been fully understood. Therefore, Li et al. (2015) [131] used B. mori to investigate the biological functions of one particular 20E inactivation enzyme, specifically the ecdysone oxidase (EO). Having a great impact on insects’ key processes, it is crucial to regulate the ecdysteroids concentration. Accordingly, the EO participates in ecdysteroids’ oxidation [132]. The CRISPR-Cas9 system was used to deplete the BmEo gene, and consequently, there was observed that the duration of the fifth instar larval phase was prolonged by 24 h [131].

Another key element that is involved in insects’ development is the juvenile hormone (JH) [133]. The central role in the JH degradation process is played by the juvenile hormone esterase (JHE). Zhang et al. (2017) [134] used CRISPR-Cas9 in B. mori to deplete the encoding gene for JHE, specifically BmJhe, in order to investigate its function. Their data showed that by knocking out this gene, the fourth and fifth instar stages were prolonged due the fact that JH metabolism was delayed. These findings are not only important for functional analysis, but also for the sericulture industry. For silk production, it is of great impact to expand the larval stages, because this leads to larger larvae, and thus, larger cocoons are produced. These findings highlight the feasibility of using genome editing tools for economic purposes [134].

Moreover, CRISPR-Cas was used in B. mori to perform epigenetic modifications. In this context, Liu et al. (2019) [135] explored the impact of methylation on silkworms’ development. This study is of great significance by providing a strategy for the investigation of DNA methylation importance in a target locus. Furthermore, their study represents a starting point for exploring the impact of DNA methylation status on different phenotypes in silkworms and beyond [135]. Furthermore, Xing et al. (2019) [136] used this technology for labeling endogenous regions in B. mori embryonic cell line and targeted the BmFibH gene. Using CRISPR-dCas9 as an imaging tool has extraordinary impact on performing fundamental research, but also on providing insights on insecticide resistance [136].

As for the detection of CRISPR-Cas induced mutations, there are a plethora of screening methods. On this topic, in a recent study, Brady et al. (2020) [137] described a new approach and provided a protocol for screening, characterizing, and stabilizing the mutant silkworm lines. The provided protocol involves several molecular methods that allow the recognition of the induced mutations on both autosomes and sex chromosomes [137].

In Table 3, several studies are described that used the CRISPR-Cas technology in order to genetically edit B. mori.

Table 3

CRISPR-Cas applications in B. mori.

Target Gene	Mutation Type	Delivery Approach	Objective	Gene Function	References
BmTim	Deletions	Plasmid	Functional gene analysis	Exhibit an impact on the embryo hatching process	[138]
BmApp	Deletions, insertions	mRNA	Functional gene analysis	Regulates wing development and cell mitosis	[139]
BmFoxo	Deletions	Plasmid	Analysis of BmFoxo and JH interaction	Involved in JH degradation	[140]
BmKmo	Deletions	Protein	Phenotypic analysis	Involved in the process of egg formation and eye coloring	[141]
BmYki	Deletions	Plasmid	Functional gene analysis	Involved in organ development and regeneration	[142]
BmTorso	Deletions	Plasmid	Functional gene analysis	Maintain the steroid hormones balance	[143]
BmEsp	Deletions	Plasmid	Functional gene analysis	Involved in a female’s reproducibility	[144]
BmTudor	Deletions, insertions	Plasmid	Investigating the frequency of homologous recombination	Included in the stress granule formation	[145]
BmIdgf	Deletions	mRNA	Analyzing the pigmentation mechanism	Plays a key role in the melanization mechanism	[146]
BmBngr-a2	Deletions	Plasmid	Exploring functional studies of certain ion transport peptides	Involved in water homeostasis	[147]
BmTctp	Deletions	Plasmid	Functional analysis	Involved in different cell process, such as growth, development, and proliferation	[148]
BmGr66	Deletions	Plasmid	A better understanding of the specific feeding preference	Involved in silkworms’ specific feeding preferences	[149]
BmOvo	Deletions	Plasmid	Functional analysis	Involved in germline sex determination and wing metamorphosis	[150]
BmPhyhd1	Deletions	Protein	Functional analysis	Exhibits a great impact on certain features of the epithelial cells	[151]
BmWnt1	Deletions	mRNA	Functional analysis	Involved in the embryogenesis	[152]
BmE75b	Deletions	mRNA	Functional analysis	Controls the developmental timing	[153]
BmOrco	Deletions	Plasmid	Exploration of adult mating behavior	Involved in silkworms’ olfactory system, being an odorant receptor co-receptor	[154]

Applicability of CRISPR-Cas in Anti-BmNPV Therapy

Viruses represent a major threat to numerous hosts, including humans or insects, being one of the most rapidly mutating biological entities. However, there is an urgent need to develop new methods to combat these pathogens. Being one of the most valued molecular tools, the CRISPR-Cas system is currently being used for developing antiviral strategies in various organisms [155,156,157]. Of all nucleases described in the specialized literature, it has been shown that different variants of Cas9, Cas12, and Cas13 exhibit the most promising results regarding the antiviral approaches [155].

In respect to insect virology, BmNPV was the first discovered pathogen (1998) [158]. In the sericulture industry, this baculovirus causes extensive economic losses; therefore, it has been intensively studied. Being part of Baculoviridae family, its infection leads to the most severe disease and only a few silkworm strains are resistant to this virus [159]. Currently, there are several traditional methods that help to enhance silkworms’ resistance to BmNPV, but they exhibit serious limitations [158,160,161].

However, CRISPR-Cas has been successfully used as an antiviral therapy in B. mori, especially against BmNPV. Chen et al. (2017) [162] selected two genes that are involved in baculovirus’ replication and propagation processes: immediate early-1 (ie-1) and me53, respectively. The authors designed two gRNA for each target gene. The engineered plasmid contained three expression cassettes. One cassette contained the Cas9 nuclease, the second one included the gRNAs and the last one harbored a selecting marker, specifically the enhanced green fluorescent protein (EGFP). Even if the silkworms’ viability and fecundity was not impacted, the transgenic homozygotes experienced a delay of the larval stage development. After performing viral inoculation in both wild type group and transgenic group, it was observed that the second category posed great viral resistance against BmNPV [162]. Likewise, another group of researchers targeted two other genes, ie-0 and ie-2, by using CRISPR-Cas9. After being inoculated with occlusion bodies of BmNPV, the survival rate reached 65% for the transgenic silkworms [163].

The multiplexable feature of the CRISPR-Cas9 technology was underlined in a study performed by Dong et al. (2019) [164]. The researchers targeted in BmNPV three different genes essential for viral replication: ie-1, the major envelope glycoprotein and the late expression factor-11. The success of this study revealed one promising strategy in inhibiting different B. mori by using the powerful CRISPR-Cas method.

5. Conclusions

B. mori is one of the most important domesticated insects due to its great potential as a biotechnological platform to produce recombinant proteins, but also because of its great success as an experimental model organism [92,110]. Due to its extraordinary prospects, a wide range of studies have been performed that accelerated fundamental research and beyond in silkworms. Being the most feasible technology in terms of genome editing, the CRISPR-Cas system is currently used in many laboratories specialized in medicine, agriculture, alimentary industry, and entomology research [34,165]. A hallmark of CRISPR-Cas is the relative ease of designing CRISPR-based experiments. As reviewed elsewhere, there are numerous bioinformatics tools that facilitate the guide RNA design, as well as the prediction and evaluation of editing results [166]. In our experience, a thorough guide RNA design can also be achieved by using standard sequence alignment tools and manual inspection of potential target regions. Considering the aforementioned, the utilization of the CRISPR-Cas system as a gene editing tool augmented the research in B. mori. Most of the studies have been focused on using the CRISPR-Cas system to perform functional gene analysis, to elucidate certain mechanisms [135,148], or to enhance silkworms’ resistance to BmNPV [164]. By reviewing the most remarkable work in this field, we provide deep insights that offer support for future research not only in B. mori, but also in other insect experimental models.

There are interesting future prospects when it comes to using the CRISPR-Cas technology in silkworms. Nowadays, by genetically manipulating the B. mori genome, major progress is being made for a better understanding of the process of fibroin and sericin synthesis, but also great advances are made to increase the understanding in the most important processes in silkworms. Notable great applications of CRISPR-Cas in B. mori refer to the development of enhanced silk fibers and the production of recombinant proteins that exhibit importance for various scientific fields. Even if compared with the other genome editing techniques, CRISPR-Cas9 exhibits lower off-target activity, although currently, it could not be confirmed that these unfavorable effects are completely eliminated. However, there is a major need to eliminate or at least reduce the off-target activity.