High-throughput and genome-scale targeted mutagenesis using CRISPR in a nonmodel multicellular organism, Bombyx mori.

Large-scale genetic mutant libraries are powerful approaches to interrogating genotype-phenotype correlations and identifying genes responsible for certain environmental stimuli, both of which are the central goal of life science study. We produced the first large-scale CRISPR-Cas9-induced library in a nonmodel multicellular organism, Bombyx mori We developed a piggyBac-delivered binary genome editing strategy, which can simultaneously meet the requirements of mixed microinjection, efficient multipurpose genetic operation, and preservation of growth-defect lines. We constructed a single-guide RNA (sgRNA) plasmid library containing 92,917 sgRNAs targeting promoters and exons of 14,645 protein-coding genes, established 1726 transgenic sgRNA lines following microinjection of 66,650 embryos, and generated 300 mutant lines with diverse phenotypic changes. Phenomic characterization of mutant lines identified a large set of genes responsible for visual phenotypic or economically valuable trait changes. Next, we performed pooled context-specific positive screens for tolerance to environmental pollutant cadmium exposure, and identified KWMTBOMO12902 as a strong candidate gene for breeding applications in sericulture industry. Collectively, our results provide a novel and versatile approach for functional B. mori genomics, as well as a powerful resource for identifying the potential of key candidate genes for improving various economic traits. This study also shows the effectiveness, practicality, and convenience of large-scale mutant libraries in other nonmodel organisms.

Effects of poly (ADP-ribose) polymerase 1 (PARP1) on silk proteins in the silkworm, Bombyx mori.

Animal silk is economically important, while silk secretion is a complex and subtle mechanism regulated by many genes. We identified the poly (ADP-ribose) polymerase (PARP1) gene of the silkworm and successfully cloned its coding sequence (CDS) sequence. Using clustered regularly interspaced short palindromic repeat (CRISPR/Cas9) technology, we screened single guide RNA (sgRNA) with high knockout efficiency by cellular experiments and obtained PARP1 mutants by knocking out the PARP1 gene of the silkworm at the individual level. We found that the mutants mainly exhibited phenotypes such as smaller cocoon size and reduced cocoon shell rate than the wild type. We also detected the expression of silk protein genes in the mutant by quantitative real-time PCR (qPCR) and found that the expression of some silk protein genes was slightly down-regulated. Meanwhile, together with the results of transcriptomic analysis, we hypothesized that PARP1 may affect the synthesis of silk proteins, resulting in their failure to function properly. Our study may provide an important reference for future in-depth refinement of the molecular mechanism of silk protein expression in silk-producing animals, as well as a potential idea for future development of molecular breeding lines of silkworms to improve silk production.

Overview and Evolution of Insect Fibroin Heavy Chain (FibH).

The FibH gene, crucial for silk spinning in insects, encodes a protein that significantly influences silk fiber mechanics. Due to its large size and repetitive sequences, limited known sequences of insect FibH impede comprehensive understanding. Here, we analyzed 114 complete FibH gene sequences from Lepidoptera (71 moths, 24 butterflies) and 13 Trichoptera, revealing single-copy FibH in most species, with 2-3 copies in Hesperinae and Heteropterinae (subfamily of skippers). All FibH genes are structured with two exons and one intron (39-45 bp), with the second exon being notably longer. Moths exhibit higher GC content in FibH compared to butterflies and Trichoptera. The FibH composition varies among species, with moths and butterflies favoring Ala, Gly, Ser, Pro, Gln, and Asn, while Trichoptera FibH is enriched in Gly, Ser, and Arg, and has less Ala. Unique to Trichoptera FibH are Tyr, Val, Arg, and Trp, whereas Lepidoptera FibH is marked by polyAla (polyalanine), polySer (polyserine), and the hexapeptide GAGSGA. A phylogenetic analysis suggests that Lepidoptera FibH evolved from Trichoptera, with skipper FibH evolving from Papilionoidea. This study substantially expands the FibH repertoire, providing a foundation for the development of artificial silk.

Genome-wide CRISPR screening reveals genes essential for cell viability and resistance to abiotic and biotic stresses in Bombyx mori.

High-throughput genetic screens are powerful methods to interrogate gene function on a genome-wide scale and identify genes responsible to certain stresses. Here, we developed a piggyBac strategy to deliver pooled sgRNA libraries stably into cell lines. We used this strategy to conduct a screen based on genome-wide clustered regularly interspaced short palindromic repeat technology (CRISPR)-Cas9 in Bombyx mori cells. We first constructed a single guide RNA (sgRNA) library containing 94,000 sgRNAs, which targeted 16,571 protein-coding genes. We then generated knockout collections in BmE cells using the piggyBac transposon. We identified 1006 genes that are essential for cell viability under normal growth conditions. Of the identified genes, 82.4% (829 genes) were homologous to essential genes in seven animal species. We also identified 838 genes whose loss facilitated cell growth. Next, we performed context-specific positive screens for resistance to biotic or nonbiotic stresses using temperature and baculovirus separately, which identified several key genes and pathways from each screen. Collectively, our results provide a novel and versatile platform for functional annotations of B. mori genomes and deciphering key genes responsible for various conditions. This study also shows the effectiveness, practicality, and convenience of genome-wide CRISPR screens in nonmodel organisms.

Comparison of Long-Read Methods for Sequencing and Assembly of Lepidopteran Pest Genomes.

Lepidopteran species are mostly pests, causing serious annual economic losses. High-quality genome sequencing and assembly uncover the genetic foundation of pest occurrence and provide guidance for pest control measures. Long-read sequencing technology and assembly algorithm advances have improved the ability to timeously produce high-quality genomes. Lepidoptera includes a wide variety of insects with high genetic diversity and heterozygosity. Therefore, the selection of an appropriate sequencing and assembly strategy to obtain high-quality genomic information is urgently needed. This research used silkworm as a model to test genome sequencing and assembly through high-coverage datasets by de novo assemblies. We report the first nearly complete telomere-to-telomere reference genome of silkworm Bombyx mori (P50T strain) produced by Pacific Biosciences (PacBio) HiFi sequencing, and highly contiguous and complete genome assemblies of two other silkworm strains by Oxford Nanopore Technologies (ONT) or PacBio continuous long-reads (CLR) that were unrepresented in the public database. Assembly quality was evaluated by use of BUSCO, Inspector, and EagleC. It is necessary to choose an appropriate assembler for draft genome construction, especially for low-depth datasets. For PacBio CLR and ONT sequencing, NextDenovo is superior. For PacBio HiFi sequencing, hifiasm is better. Quality assessment is essential for genome assembly and can provide better and more accurate results. For chromosome-level high-quality genome construction, we recommend using 3D-DNA with EagleC evaluation. Our study references how to obtain and evaluate high-quality genome assemblies, and is a resource for biological control, comparative genomics, and evolutionary studies of Lepidopteran pests and related species.

Let-7 microRNA is a critical regulator in controlling the growth and function of silk gland in the silkworm.

The silk gland is characterized by high protein synthesis. However, the molecular mechanisms controlling silk gland growth and silk protein synthesis remain undetermined. Here we demonstrated that CRISPR/Cas9-based knockdown of let-7 or the whole cluster promoted endoreduplication and enlargement of the silk gland, accompanied by changing silk yield, whereas transgenic overexpression of let-7 led to atrophy and degeneration of the silk gland. Mechanistically, let-7 controls cell growth in the silk gland through coordinating nutrient metabolism processes and energy signalling pathways. Transgenic overexpression of pyruvate carboxylase, a novel target of let-7, resulted in enlargement of the silk glands, which is consistent with the abnormal phenotype of the let-7 knockdown. Overall, our data reveal a previously unknown miRNA-mediated regulation of silk gland growth and physiology and shed light on involvement of let-7 as a critical stabilizer and booster in carbohydrate metabolism, which may have important implications for understanding of the molecular mechanism and physiological function of specialized organs in other species.

High-Throughput Screening Identifies Two Novel Small Molecule Enhancers of Recombinant Protein Expression.

As a primary strategy for production of biological drugs, recombinant proteins produced by transient transfection of mammalian cells are essential for both basic research and industrial production. Here, we established a high-throughput screening platform for improving the expression levels of recombinant proteins. In total, 10,011 small molecule compounds were screened through our platform. After two rounds of screening, we identified two compounds, Apicidin and M-344, that significantly enhanced recombinant protein expression. Both of the selected compounds were histone deacetylase inhibitors, suggesting that the two small molecules increased the expression levels of recombinant proteins by promoting histone acetylation. Moreover, both molecules showed low cytotoxicity. Therefore, our findings suggest that these small molecules may have wide applications in the future.

Deep Insight into the Transcriptome of the Single Silk Gland of Bombyx mori.

The silk gland synthesizes and secretes a large amount of protein and stores liquid silk protein at an extremely high concentration. Interestingly, silk proteins and serine protease inhibitors are orderly arranged in the silk gland lumen and cocoon shells. Silk fiber formation and the spinning mechanism have not been fully elucidated. Therefore, we conducted a comparative transcriptome analysis of seven segments of the single silk gland to characterize internal changes in the silk gland during the 5th instar of mature larvae. In total, 3121 differentially expressed genes were identified in the seven segments. Genes highly expressed in the middle silk gland (MSG) were mainly involved in unsaturated fatty acid biosynthesis, fatty acid metabolism, apoptosis-fly, and lysosome pathways, whereas genes highly expressed in the posterior silk gland (PSG) were mainly involved in ribosome, proteasome, citrate cycle, and glycolysis/gluconeogenesis pathways. Thus, the MSG and PSG differ greatly in energy source use and function. Further, 773 gradually upregulated genes (from PSG to MSG) were involved in energy metabolism, silk protein synthesis, and secretion, suggesting that these genes play an important role in silk fiber formation. Our findings provide insights into the mechanism of silk protein synthesis and transport and silk fiber formation.

Bacillus bombysepticus α-Toxin Binding to G Protein-Coupled Receptor Kinase 2 Regulates cAMP/PKA Signaling Pathway to Induce Host Death.

Bacterial pathogens and their toxins target host receptors, leading to aberrant behavior or host death by changing signaling events through subversion of host intracellular cAMP level. This is an efficient and widespread mechanism of microbial pathogenesis. Previous studies describe toxins that increase cAMP in host cells, resulting in death through G protein-coupled receptor (GPCR) signaling pathways by influencing adenylyl cyclase or G protein activity. G protein-coupled receptor kinase 2 (GRK2) has a central role in regulation of GPCR desensitization. However, little information is available about the pathogenic mechanisms of toxins associated with GRK2. Here, we reported a new bacterial toxin-Bacillus bombysepticus (Bb) α-toxin that was lethal to host. We showed that Bb α-toxin interacted with BmGRK2. The data demonstrated that Bb α-toxin directly bound to BmGRK2 to promote death by affecting GPCR signaling pathways. This mechanism involved stimulation of Gαs, increase level of cAMP and activation of protein kinase A (PKA). Activated cAMP/PKA signal transduction altered downstream effectors that affected homeostasis and fundamental biological processes, disturbing the structural and functional integrity of cells, resulting in death. Preventing cAMP/PKA signaling transduction by inhibitions (NF449 or H-89) substantially reduced the pathogenicity of Bb α-toxin. The discovery of a toxin-induced host death specifically linked to GRK2 mediated signaling pathway suggested a new model for bacterial toxin action. Characterization of host genes whose expression and function are regulated by Bb α-toxin and GRK2 will offer a deeper understanding of the pathogenesis of infectious diseases caused by pathogens that elevate cAMP.

A Simple Method for the Cross-Section Area Determination of Single Profiled Fibers and Its Application.

One of the critical prerequisites for accurately measuring the mechanical properties of profiled fibers is the precise determination of their cross-sectional areas (CSAs). In this study, a new method is established for determining a single profiled fibers' CSA based on the frozen section method and digital photo, pixel-ratio method (FS-DP). FS-DP is used to obtain a transverse section of a fiber, by acquiring an image of the cross section using optical microscopy or scanning electron microscopy, and then calculating the CSA using Photoshop. Using FS-DP, it was found that the shape of a fiber of silk changes little in a range of 50 µm, but varies considerably over a range of 1 m, while the CSA of cocoon silk (900 m) first increases and then decreases. Mechanical property tests showed that the elongation, strength, elastic modulus, and toughness values of the cocoon silk are consistent with those reported previously. Additionally, FS-DP was also used to observe other profiled fibers. The application tests indicated that FS-DP can be used to quickly and accurately obtain the CSA of a single profiled fiber, and that it is suitable for the large-scale determination and analysis of the mechanical properties of profiled fibers.

Increasing the yield of middle silk gland expression system through transgenic knock-down of endogenous sericin-1.

Various genetically modified bioreactor systems have been developed to meet the increasing demands of recombinant proteins. Silk gland of Bombyx mori holds great potential to be a cost-effective bioreactor for commercial-scale production of recombinant proteins. However, the actual yields of proteins obtained from the current silk gland expression systems are too low for the proteins to be dissolved and purified in a large scale. Here, we proposed a strategy that reducing endogenous sericin proteins would increase the expression yield of foreign proteins. Using transgenic RNA interference, we successfully reduced the expression of BmSer1 to 50%. A total 26 transgenic lines expressing Discosoma sp. red fluorescent protein (DsRed) in the middle silk gland (MSG) under the control of BmSer1 promoter were established to analyze the expression of recombinant. qRT-PCR and western blotting showed that in BmSer1 knock-down lines, the expression of DsRed had significantly increased both at mRNA and protein levels. We did an additional analysis of DsRed/BmSer1 distribution in cocoon and effect of DsRed protein accumulation on the silk fiber formation process. This study describes not only a novel method to enhance recombinant protein expression in MSG bioreactor, but also a strategy to optimize other bioreactor systems.

Genetic breeding of silkworms: from traditional hybridization to molecular design.

Sericulture is one of the great inventions of the Chinese people and has become an important cultural feature of China. As China is the long-lasting center of silk production, genetic breeding of silkworm was highly developed historically, and has formed a comprehensive system for breeding and preservation of new varieties. However, silkworm breeding reached a bottleneck recently, because most of the traditional genetic resources have been utilized and silkworm strains have become homogeneous. Meanwhile, sericulture in China meets huge challenges in the 21st century. In recent years, with the development and rapid application of molecular biology, genomics, transgene and genome editing, silkworm genetic breeding has entered a new era. In this review, we summarize the development of silkworm genetic breeding, especially the progress and perspective of transgene and genome editing in genetic engineering of silkworms. We also discuss the future development of silkworm genetic breeding.

The C-terminus of DSX(F5) protein acts as a novel regulatory domain in Bombyx mori.

The doublesex gene regulates the somatic sexual development of Bombyx mori by alternatively splicing into sex-specific splice forms. In our previous study, the splice form Bmdsx (F7) , which encodes the BmDSX(F5) protein, was found to be expressed in a female-specific manner and to contain a novel C-terminus. In this study, we aimed to investigate the role of this C-terminus. Two transgenic lines, L1 and L2, were constructed to ectopically express Bmdsx (F7) in males. Phenotype and W chromosome-specific polymerase chain reaction (PCR) analysis showed that developmental abnormalities and sex reversal did not occur. Moreover, the sex ratio was also normal. Quantitative PCR revealed that the expression levels of SP1 and Vg were upregulated in the fat body of transgenic males. Additionally, the expression level of PBP was downregulated in the antenna of transgenic males. The results suggested that the C-terminus of BmDSX(F5) functioned as a regulatory domain during regulation of downstream target gene expression and that BmDSX(F5) participated in the sexual development of somatic cells together with other DSX proteins in B. mori.

CRISPR/Cas9-mediated targeted mutagenesis in Nicotiana tabacum.

Genome editing is one of the most powerful tools for revealing gene function and improving crop plants. Recently, RNA-guided genome editing using the type II clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein (Cas) system has been used as a powerful and efficient tool for genome editing in various organisms. Here, we report genome editing in tobacco (Nicotiana tabacum) mediated by the CRISPR/Cas9 system. Two genes, NtPDS and NtPDR6, were used for targeted mutagenesis. First, we examined the transient genome editing activity of this system in tobacco protoplasts, insertion and deletion (indel) mutations were observed with frequencies of 16.2-20.3% after transfecting guide RNA (gRNA) and the nuclease Cas9 in tobacco protoplasts. The two genes were also mutated using multiplexing gRNA at a time. Additionally, targeted deletions and inversions of a 1.8-kb fragment between two target sites in the NtPDS locus were demonstrated, while indel mutations were also detected at both the sites. Second, we obtained transgenic tobacco plants with NtPDS and NtPDR6 mutations induced by Cas9/gRNA. The mutation percentage was 81.8% for NtPDS gRNA4 and 87.5% for NtPDR6 gRNA2. Obvious phenotypes were observed, etiolated leaves for the psd mutant and more branches for the pdr6 mutant, indicating that highly efficient biallelic mutations occurred in both transgenic lines. No significant off-target mutations were obtained. Our results show that the CRISPR/Cas9 system is a useful tool for targeted mutagenesis of the tobacco genome.

Comparative proteomic analysis of silkworm fat body after knocking out fibroin heavy chain gene: a novel insight into cross-talk between tissues.

Cross-talk between tissues plays key roles in development of organisms; however, there are few researches on cross-talk between tissues in insects. Our previous studies showed that the pupal body weight was elevated after knocking out the fibroin heavy chain gene (BmFib-H), whereas the gene specifically expressed in silk glands of silkworm. Hence, the mutant is a good material for studying the cross-talk between tissues. It is considered that the fat body of silkworm during larval stage is used to store nutrients for pupal development. Herein, comparative proteomic of fat body on the 5th day of fifth instar was performed between BmFib-H gene knock-out Bombyx mori line (FGKO) and its wide-type Dazao. These results revealed that a single gene knock-out in silk gland triggered large-scale metabolic pathways changes in fat body. The levels of proteins involved in glycolysis/gluconeogenesis, pentose phosphate pathway, and glycine-serine biosynthetic pathway were down-regulated in the FGKO fat body. In contrast, the abundances of many proteins participating in protein synthesis, including ribosomal proteins, eukaryotic translation initiation factor, and elongation factor, were up-regulated. Moreover, the concentrations of glycogen and proteins in the FGKO fat body were greatly increased. These findings provided a novel insight into the cross-talk between silk gland and fat body in silkworm, and the presence of cross-talk between silk gland and fat body could regulate the redistribution of nutrients in the FGKO fat body leading to the increase of the pupal weight.

Remobilizing deleted piggyBac vector post-integration for transgene stability in silkworm.

Deletion of transposable elements post-genomic integration holds great promise for stability of the transgene in the host genome and has an essential role for the practical application of transgenic animals. In this study, a modified piggyBac vector that mediated deletion of the transposon sequence post-integration for transgene stability in the economically important silkworm Bombyx mori was constructed. The piggyBac vector architecture contains inversed terminal repeat sequences L1, L2 and R1, which can form L1/R1 and L2/R1 types of transposition cassettes. hsp70-PIG as the piggyBac transposase expression cassette for initial transposition, further remobilization and transgene stabilization test was transiently expressed in a helper vector or integrated into the modified vector to produce a transgenic silkworm. Shortening L2 increased the transformation frequency of L1/R1 into the silkworm genome compared to L2/R1. After the integration of L1/R1 into the genome, the remobilization of L2/R1 impaired the transposon structure and the resulting transgene linked with an impaired transposon was stable in the genome even in the presence of exogenously introduced transposase, whereas those flanked by the intact transposon were highly mobile in the genome. Our results demonstrated the feasibility of post-integration deletion of transposable elements to guarantee true transgene stabilization in silkworm. We suggest that the modified vector will be a useful resource for studies of transgenic silkworms and other piggyBac-transformed organisms.

Modifying the Mechanical Properties of Silk Fiber by Genetically Disrupting the Ionic Environment for Silk Formation.

Silks are widely used biomaterials, but there are still weaknesses in their mechanical properties. Here we report a method for improving the silk fiber mechanical properties by genetic disruption of the ionic environment for silk fiber formation. An anterior silk gland (ASG) specific promoter was identified and used for overexpressing ion-transporting protein in the ASG of silkworm. After isolation of the transgenic silkworms, we found that the metal ion content, conformation and mechanical properties of transgenic silk fibers changed accordingly. Notably, overexpressing endoplasmic reticulum Ca2+-ATPase in ASG decreased the calcium content of silks. As a consequence, silk fibers had more α-helix and ß-sheet conformations, and their tenacity and extension increased significantly. These findings represent the in vivo demonstration of a correlation between metal ion content in the spinning duct and the mechanical properties of silk fibers, thus providing a novel method for modifying silk fiber properties.

Multiplex genomic structure variation mediated by TALEN and ssODN.

BACKGROUND: Genomic structure variation (GSV) is widely distributed in various organisms and is an important contributor to human diversity and disease susceptibility. Efficient approaches to induce targeted genomic structure variation are crucial for both analytic and therapeutic studies of GSV. Here, we presented an efficient strategy to induce targeted GSV including chromosomal deletions, duplications and inversions in a precise manner. RESULTS: Utilizing Transcription Activator-Like Effector Nucleases (TALEN) designed to target two distinct sites, we demonstrated targeted deletions, duplications and inversions of an 8.9 Mb chromosomal segment, which is about one third of the entire chromosome. We developed a novel method by combining TALEN-induced GSV and single stranded oligodeoxynucleotide (ssODN) mediated gene modifications to reduce unwanted mutations occurring during the targeted GSV using TALEN or Zinc finger nuclease (ZFN). Furthermore, we showed that co-introduction of TALEN and ssODN generated unwanted complex structure variation other than the expected chromosomal deletion. CONCLUSIONS: We demonstrated the ability of TALEN to induce targeted GSV and provided an efficient strategy to perform GSV precisely. Furthermore, it is the first time to show that co-introduction of TALEN and ssODN generated unwanted complex structure variation. It is plausible to believe that the strategies developed in this study can be applied to other organisms, and will help understand the biological roles of GSV and therapeutic applications of TALEN and ssODN.

Overexpression and functional characterization of an Aspergillus niger phytase in the fat body of transgenic silkworm, Bombyx mori.

In a previous study, we isolated 1,119 bp of upstream promoter sequence from Bmlp3, a gene encoding a member of the silkworm 30 K storage protein family, and demonstrated that it was sufficient to direct fat body-specific expression of a reporter gene in a transgenic silkworm, thus highlighting the potential use of this promoter for both functional genomics research and biotechnology applications. To test whether the Bmlp3 promoter can be used to produce recombinant proteins in the fat body of silkworm pupae, we generated a transgenic line of Bombyx mori which harbors a codon-optimized Aspergillus niger phytase gene (phyA) under the control of the Bmlp3 promoter. Here we show that the Bmlp3 promoter drives high levels of phyA expression in the fat body, and that the recombinant phyA protein is highly active (99.05 and 54.80 U/g in fat body extracts and fresh pupa, respectively). We also show that the recombinant phyA has two optimum pH ranges (1.5-2.0 and 5.5-6.0), and two optimum temperatures (55 and 37 °C). The activity of recombinant phyA was lost after high-temperature drying, but treating with boiling water was less harmful, its residual activity was approximately 84% of the level observed in untreated samples. These results offer an opportunity not only for better utilization of large amounts of silkworm pupae generated during silk production, but also provide a novel method for mass production of low-cost recombinant phytase using transgenic silkworms.

High-Throughput Screening of PAM-Flexible Cas9 Variants for Expanded Genome Editing in the Silkworm (Bombyx mori).

Genome editing provides novel opportunities for the precise genome engineering of diverse organisms. Significant progress has been made in the development of genome-editing tools for Bombyx mori (B. mori) in recent years. Among these, CRISPR/Cas9, which is currently the most commonly used system in lepidopteran insects, recognizes NGG protospacer adjacent motif (PAM) sequences within the target locus. However, Cas9 lacks the ability to target all gene loci in B. mori, indicating the need for Cas9 variants with a larger editing range. In this study, we developed a high-throughput screening platform to validate Cas9 variants at all possible recognizable and editable PAM sites for target sequences in B. mori. This platform enabled us to identify PAM sites that can be recognized by both xCas9 3.7 and SpCas9-NG variants in B. mori and to assess their editing efficiency. Cas9 shows PAM sites every 13 base pairs in the genome, whereas xCas9 3.7 and SpCas9-NG have an average distance of 3.4 and 3.6 base pairs, respectively, between two specific targeting sites. Combining the two Cas9 variants could significantly expand the targeting range of the genome, accelerate research on the B. mori genome, and extend the high-throughput rapid screening platform to other insects, particularly those lacking suitable NGG PAM sequences.