Genome-wide analysis of the mulberry (Morus abla L.) GH9 gene family and the functional characterization of MaGH9B6 during the development of the abscission zone.

Plant glycoside hydrolase family 9 genes (GH9s) are widely distributed in plants and involved in a variety of cellular and physiological processes. In the current study, nine GH9 genes were identified in the mulberry and were divided into two subfamilies based on the phylogenetic analysis. Conserved motifs and gene structure analysis suggested that the evolution of the two subfamilies is relatively conserved and the glycoside hydrolase domain almost occupy the entire coding region of the GH9s gene. Only segmental duplication has played a role in the expansion of gene family. Collinearity analysis showed that mulberry GH9s had the closest relationship with poplar GH9s. MaGH9B1, MaGH9B6, MaGH9B5, and MaGH9B3 were detected to have transcript accumulation in the stalk of easy-to drop mature fruit drop, suggesting that these could play a role in mulberry fruit drop. Multiple cis-acting elements related to plant hormones and abiotic stress responses were found in the mulberry GH9 promoter regions and showed different activities under exogenous abscisic acid (ABA) and 2,4- dichlorophenoxyacetic acid (2,4-D) stresses. We found that the lignin content in the fruit stalk decreased with the formation of the abscission zone (AZ), which could indirectly reflect the formation process of the AZ. These results provide a theoretical basis for further research on the role of GH9s in mulberry abscission.

MaABI5 and MaABF1 transcription factors regulate the expression of MaJOINTLESS during fruit abscission in mulberry (Morus alba L.).

Mulberry holds significant economic value. However, during the ripening stage of its fruit, the phenomenon of abscission, resulting in heavy fruit drop, can severely impact the yield. The formation of off-zone structures is a critical factor in the fruit abscission process, and this process is regulated by multiple transcription factors. One such key gene that plays a significant role in the development of the off-zone in the model plant tomato is JOINTLESS, which promotes the expression of abscission-related genes and regulates the differentiation of abscission zone tissue cells. However, there is a lack of information about fruit abscission mechanism in mulberry. Here, we analyzed the MaJOINTLESS promoter and identified the upstream regulators MaABF1 and MaABI5. These two regulators showed binding with MaJOINTLESS promoter MaABF1 (the ABA Binding Factor/ABA-Responsive Element Binding Proteins) activated the expression of MaJOINTLESS, while MaABI5 (ABSCISIC ACID-INSENSITIVE 5) inhibited the expression of MaJOINTLESS. Finally, the differentially expressed genes (DEGs) were analyzed by transcriptome sequencing to investigate the expression and synergistic relationship of endogenous genes in mulberry during abscission. GO classification and KEGG pathway enrichment analysis showed that most of the DEGs were concentrated in MAPK signaling pathway, flavonoid biosynthesis, citric acid cycle, phytohormone signaling, amino acid biosynthesis, and glycolysis. These results provide a theoretical basis for subsequent in-depth study of physiological fruit abscission in mulberry.

Specific Expression of Antimicrobial Peptides from the Black Soldier Fly in the Midgut of Silkworms (Bombyx mori) Regulates Silkworm Immunity.

Antimicrobial peptides are molecules with strong antimicrobial activity and are of substantial interest for the immunization of insects. As a type of dipteran insect that can turn organic waste into animal feed, the black soldier fly (BSF) can "turn waste into treasure". In this study, we investigated the antimicrobial activity of the antimicrobial peptide genes, HiCG13551 and Hidiptericin-1, of BSF in silkworms, by overexpressing the genes specifically in the midgut. Changes in the mRNA levels of the transgenic silkworms after infection with Staphylococcus aureus were evaluated using transcriptome sequencing. The results showed that Hidiptericin-1 had stronger antimicrobial activity than HiCG13551. KEGG enrichment analysis showed that the differentially expressed genes in the transgenic overexpressed Hidiptericin-1 silkworm lines from the D9L strain were mainly enriched in the starch and sucrose metabolism, pantothenate and CoA biosynthesis, drug metabolism (other enzymes), biotin metabolism, platinum drug resistance, galactose metabolism, and pancreatic secretion pathways. In addition, immune-related genes were up-regulated in this transgenic silkworm strain. Our study may provide new insights for future immune studies on insects.

Metabonomic Analysis of Silkworm Midgut Reveals Differences between the Physiological Effects of an Artificial and Mulberry Leaf Diet.

Bombyx mori is a model lepidopteran insect of great economic value. Mulberry leaves are its only natural food source. The development of artificial diets can not only resolve the seasonal shortage of mulberry leaves but also enable changes to be made to the feed composition according to need. Metabolomic differences between the midguts of male and female silkworms fed either on fresh mulberry leaves or an artificial diet were studied using liquid chromatography-mass spectrography (LC-MS/MS) analysis. A total of 758 differential metabolites were identified. Our analysis showed that they were mainly involved in disease resistance and immunity, silk quality, and silkworm growth and development. These experimental results provide insights into the formulation of optimized artificial feed for silkworms.

Effects of Serine/Arginine Enriched Protein BmUP on the Development of Male Silkworm Reproductive Organs.

Serine/arginine-rich proteins are a class of highly conserved splicing factor proteins involved in constitutive and alternative splicing. We screened a low molecular weight serine/arginine rich protein from silkworms and named it BmUP. Temporal and spatial expression analysis indicated that the BmUP gene was specifically expressed in the silkworm testis, and the highest expression occurred in the pre-pupa stage from the fifth instar to the moth stages. Here, we generated BmUP knockout individuals with the CRISPR/Cas9 system. Both the internal and external genitalia of knockout individuals were abnormal in knockout compared with wild-type male silkworms. In transgenic silkworms overexpressing BmUP, male silkworms showed a phenotype similar to that of the knockout individuals, whereas female individuals showed no significant differences from the wild type. In addition, by conducting promoter analysis, we identified Bmachi, a transcription factor that regulates the BmUP gene. Gel migration experiments revealed that BmAchi specifically binds the BmUP promoter. Quantitative real-time PCR showed that an increase in Bmachi expression up-regulated the expression of BmUP. In contrast, when the expression of Bmachi decreased, the expression of BmUP also downregulated in the experimental group compared with the control group. These results provide new insights for studying the effects of serine/arginine-rich proteins on the development of silkworm genitals.

Midgut and Head Transcriptomic Analysis of Silkworms Reveals the Physiological Effects of Artificial Diets.

Silkworms, a model lepidopteran insect, have a very simple diet. Artificial diets as an alternative nutrient source for silkworms are gradually being developed. To understand the effects of various nutrients on the growth and development of silkworms, we studied the transcriptomic differences in the midgut and head tissues of male and female silkworms fed either fresh mulberry leaves or artificial diets. In the artificial diet group, compared with the control group (fed mulberry leaves), 923 and 619 differentially expressed genes (DEGs) were identified from the midgut, and 2969 and 3427 DEGs were identified from the head, in female and male silkworms. According to our analysis, the DEGs were mainly involved in the digestion and absorption of nutrients and silkworm innate immunity. These experimental results provide insights into the effects of different foods, such as artificial diets or fresh mulberry leaves, on silkworms.

Long Non-Coding RNA Bmdsx-AS1 Effects on Male External Genital Development in Silkworm.

Long non-coding RNAs (lncRNAs) have been suggested to play important roles in some biological processes. However, the detailed mechanisms are not fully understood. We previously identified an antisense lncRNA, Bmdsx-AS1, that is involved in pre-mRNA splicing of the sex-determining gene Bmdsx in the silkworm. In this study, we analyzed the changes in the male external genitalia of transgenic overexpressed Bmdsx-AS1 silkworm lines and analyzed downstream and upstream responses. We found that Bmdsx-AS1 transgenic silkworms, compared with wild type, showed more claspers in the male external genitalia. Quantitative real-time PCR (qPCR) results indicated that overexpression of Bmdsx-AS1 decreased the expression of genes in the EGFR signaling pathway. Knockdown of Bmdsx-AS1 increased the activity of the EGFR pathway. Through promoter prediction, promoter truncation and electrophoretic mobility shift assay (EMSA) analyses, we found that the protein encoded by the Hox gene BmAbd-B specifically binds to the promoter of Bmdsx-AS1. Moreover, overexpression of BmAbd-B in the silkworm BmE cell line indicated that BmAbd-B negatively regulates the mRNA expression of Bmdsx-AS1. Our study provides insights into the regulatory mechanism of the lncRNA in the silkworm.

Spliceosomal Protein Gene BmSPX Regulates Reproductive Organ Development in Bombyx mori.

Sex determination and differentiation are nearly universal to all eukaryotic organisms, encompassing diverse systems and mechanisms. Here, we identified a spliceosomal protein gene BmSPX involved in sex determination of the lepidopeteran insect, Bombyx mori. In a transgenic silkworm line that overexpressed the BmSPX gene, transgenic silkworm males exhibited differences in their external genitalia compared to wild-type males, but normal internal genitalia. Additionally, transgenic silkworm females exhibited a developmental disorder of the reproductive organs. Upregulation of BmSPX significantly increased the expression levels of sex-determining genes (BmMasc and BmIMP) and reduced the female-type splice isoform of Bmdsx, which is a key switch gene downstream of the sex-determination pathway. Additionally, co-immunoprecipitation assays confirmed an interaction between the BmSPX protein and BmPSI, an upstream regulatory factor of Bmdsx. Quantitative real-time PCR showed that BmSPX over-expression upregulated the expression of the Hox gene abdominal-B (Adb-B), which is required for specification of the posterior abdomen, external genitalia, and gonads of insects, as well as the genes in the Receptor Tyrosine Kinase (RTK) signaling pathway. In conclusion, our study suggested the involvement of BmSPX, identified as a novel regulatory factor, in the sex-determination pathway and regulation of reproductive organ development in silkworms.

A Novel Splice Variant of the Masculinizing Gene Masc with piRNA-Cleavage-Site Defect Functions in Female External Genital Development in the Silkworm, Bombyx mori.

In the silkworm, the sex-determination primary signal Fem controls sex differentiation by specific binding of Fem-derived piRNA to the cleavage site in Masc mRNA, thus inhibiting Masc protein production in the female. In this study, we identified a novel splicing isoform of Masc, named Masc-S, which lacks the intact sequence of the cleavage site, encoding a C-terminal truncated protein. Results of RT-PCR showed that Masc-S was expressed in both sexes. Over-expression of Masc-S and Masc in female-specific cell lines showed that Masc-S could be translated against the Fem-piRNA cut. By RNA-protein pull-down, LC/MS/MS, and EMSA, we identified a protein BmEXU that specifically binds to an exclusive RNA sequence in Masc compared to Masc-S. Knockdown of Masc-S resulted in abnormal morphology in female external genital and increased expression of the Hox gene Abd-B, which similarly occurred by Bmexu RNAi. These results suggest that the splice variant Masc-S against Fem-piRNA plays an important role in female external genital development, of which function is opposite to that of full-length Masc. Our study provides new insights into the regulatory mechanism of sex determination in the silkworm.

An antisense lncRNA functions in alternative splicing of Bmdsx in the silkworm, Bombyx mori.

Long non-coding RNAs (lncRNAs) are recently thought to play important roles in some physiological processes. In this study, we identified a lncRNA (named as Bmdsx-AS1) which is the antisense transcript and locates on the position of the crucial sex-determining gene Bmdsx in the silkworm. Quantitative real-time PCR and Fluorescence in situ hybridization demonstrated that Bmdsx-AS1 was highly expressed in the silkworm testis. After knock-down of Bmdsx-AS1, the splicing pattern of Bmdsx pre-mRNA was altered in male silkworm. Transgenic overexpression of Bmdsx-AS1 indicated male-specific splicing form of Bmdsx arose in female silkworm. Furthermore, by using RNA-protein pull down, LC-MS/MS and EMSA, we found a splicing factor hnRNPH binding specifically to Bmdsx-AS1. Co-Immunoprecipitation suggested that hnRNPH interacted with BmPSI, one of the upstream regulating factors of Bmdsx. Thus, our results suggested that the antisense lncRNA Bmdsx-AS1 was involved in alternative splicing of Bmdsx in the silkworm.

P-element Somatic Inhibitor Protein Binding a Target Sequence in dsx Pre-mRNA Conserved in Bombyx mori and Spodoptera litura.

Bombyx mori doublesex (Bmdsx) functions as a double-switch gene in the final step of the sex-determination cascade in the silkworm Bombyx mori. The P-element somatic inhibitor (PSI) protein in B. mori interacts with Bmdsx pre-mRNA in CE1 as an exonic splicing silencer to promote male-specific splicing of Bmdsx. However, the character of the interaction between BmPSI and Bmdsx pre-mRNA remains unclear. Electrophoretic mobility shift assay (EMSA) results showed that the four KH_1 motifs in BmPSI are all essential for the binding, especially the former two KH_1 motifs. Three active sites (I116, L127, and IGGI) in the KH_1 motif were found to be necessary for the binding through EMSA, circular dichroism (CD) spectroscopy, and isothermal titration calorimetry (ITC). The PSI homologous protein in S. litura (SlPSI) was purified and the binding of SlPSI and CE1 was verified. Compared with BmPSI, the mutant SlPSI proteins of I116 and IGGI lost their ability to bind to CE1. In conclusion, the binding of PSI and dsx pre-mRNA are generally conserved in both B. mori and S. litura. These findings provide clues for sex determination in Lepidoptera.

Zinc finger protein rotund deficiency affects development of the thoracic leg in Bombyx mori.

The insect limb develops from the imaginal disc or larval leg during metamorphosis. The molecular mechanisms involved in the development from the larval to the adult leg are poorly understood. Herein, we cloned the full length of a zinc finger gene rotund from Bombyx mori (Bmrn), which contained a 1419 bp open reading frame, and encoded a 473 amino acid protein. Reverse transcription polymerase chain reaction and Western blot analyses demonstrated that Bmrn was expressed at higher levels in the epidermis than in other tissues tested, and it showed a very high expression level during metamorphosis. Knock-down of Bmrn produced defects in the tarsus and pretarsus, including the fusion and reduction of tarsomeres, and the developmental arrest of pretarsus. Our data showed that Bmrn is involved in the formation of the tarsus and pretarsus, whereas its homologous gene in Drosophila has been shown to affect three tarsal segments (t2-t4), suggesting that the remodeling of the leg has involved changes in the patterning of gene regulation during evolution.

Nuclear hormone receptor BmFTZ-F1 is involved in regulating the fibroin heavy chain gene in the silkworm, Bombyx mori.

BACKGROUND: The synthesis of silk protein is controlled by hormones. The expression of the nuclear hormone Bmftz-f1 in the posterior silk gland (PSG) is induced by 20-hydroxyecdysone in vivo and in vitro. However, whether Bmftz-f1 regulates silk protein expression is unknown. METHODS: In our study, western blotting and quantitative polymerase chain reactions were conducted to detect the expression of FTZ-F1 in the PSG. Electrophoretic mobility shift, chromatin immunoprecipitation, far-western blotting, bimolecular fluorescence complementation, and dual luciferase reporter assays were performed to investigate the effect of FTZ-F1 on the fibH promoter. RESULTS: (1) The expression of the hormone receptor BmFTZ-F1 was opposite to that of fibH. It was highly expressed in the PSG during the fourth molting stage and the beginning of the fifth instar, and then its expression decreased gradually until it disappeared at the end of the fifth instar and the wandering stage. (2) We identified a FTZ-F1 response element 390bp upstream of the transcription initiation site of the fibH promoter. (3) BmFTZ-F1 interacted with the basic helix-loop-helix transcription factor Bmdimm. (4) BmFTZ-F1 down-regulated fibH promoter activity and counteracted the effect of Bmdimm on fibH expression. CONCLUSIONS: Integrating these results, we conclude that BmFTZ-F1 regulates the transcription of fibH by binding to the FTZ-F1 response element in the fibH promoter and counteracts the effect of Bmdimm on fibH expression. GENERAL SIGNIFICANCE: These findings provide new insights into the mechanism of regulation of the silk protein gene.

Identification of MBF2 family genes in Bombyx mori and their expression in different tissues and stages and in response to Bacillus bombysepticus infection and starvation.

The Multiprotein bridge factor 2 (MBF2) gene was first identified as a co-activator involved in BmFTZ-F1-mediated activation of the Fushi tarazu gene. Herein, nine homologous genes of MBF2 gene are identified. Evolutionary analysis showed that this gene family is insect-specific and that the family members are closely related to response to pathogens (REPAT) genes. Tissue distribution analysis revealed that these genes could be expressed in a tissue-specific manner. Developmental profiles analysis showed that the MBF2 gene family members were highly expressed in the different stages. Analysis of the expression patterns of nine MBF2 family genes showed that Bacillus bombysepticus treatment induced the up-regulation of several MBF2 family genes, including MBF2-4, -7, -9, -8. Furthermore, we found the MBF2 family genes were modulated by starvation and the expression of these genes recovered upon re-feeding, except for MBF2-5, -9. These findings suggested roles for these proteins in insect defense against pathogens and nutrient metabolism, which has an important guiding significance for designing pest control strategies.

Detection and characterization of Wolbachia infection in silkworm.

Wolbachia naturally infects a wide variety of arthropods, where it plays important roles in host reproduction. It was previously reported that Wolbachia did not infect silkworm. By means of PCR and sequencing we found in this study that Wolbachia is indeed present in silkworm. Phylogenetic analysis indicates that Wolbachia infection in silkworm may have occurred via transfer from parasitic wasps. Furthermore, Southern blotting results suggest a lateral transfer of the wsp gene into the genomes of some wild silkworms. By antibiotic treatments, we found that tetracycline and ciprofloxacin can eliminate Wolbachia in the silkworm and Wolbachia is important to ovary development of silkworm. These results provide clues towards a more comprehensive understanding of the interaction between Wolbachia and silkworm and possibly other lepidopteran insects.

Ectopic expression of the male BmDSX affects formation of the chitin plate in female Bombyx mori.

Mating structures are involved in successful copulation, intromission, and/or insemination. These structures enable tight coupling between external genitalia of two sexes. During Bombyx mori copulation, the double harpagones in the external genitalia of males clasp the female chitin plate, which is derived from the larval eighth abdominal segment; abnormal development of the female chitin plate affects copulation. We report that ERK phosphorylation (p-ERK) and expression of Abdominal-B (Abd-B) in the posterior abdomen of the female adult is lower than in the male. Ectopic expression of the male-specific spliced form of B. mori doublesex (Bmdsx(M)) in females, however, up-regulates Abd-B and spitz (spi) expression, increasing EGFR signaling activity, and thus forming an abnormal chitin plate and reduced female copulation. These findings indicate that Bmdsx affects the development of the eighth abdominal segment by regulating the activity of EGFR signaling and the expression of Abd-B, resulting in an extra eighth abdominal segment (A8) in males versus the loss of this segment in adult females.

New insights into the genomic organization and splicing of the doublesex gene, a terminal regulator of sexual differentiation in the silkworm Bombyx mori.

Sex-determination mechanisms differ among organisms. The primary mechanism is diverse, whereas the terminal regulator is relatively-conserved. We analyzed the transcripts of the Bombyx mori doublesex gene (Bmdsx), and reported novel results concerning the genomic organization and expression of Bmdsx. Bmdsx consists of nine exons and eight introns, of which two exons are novel and have not been reported previously. Bmdsx transcripts are spliced to generate seventeen alternatively-spliced forms and eleven putative trans-spliced variants. Thirteen of the alternatively-spliced forms and five of the putative trans-spliced forms are reported here for the first time. Sequence analysis predicts that ten female-specific, six male-specific splice forms and one splice form found in males and females will result in four female-specific, two male-specific Dsx proteins and one Dsx protein common to males and females. The Dsx proteins are expected to be functional and regulate downstream target genes. Some of the predicted Dsx proteins are described here for the first time. Therefore the expression of the dsx gene in B. mori results in a variety of cis- and trans-spliced transcripts and multiple Dsx proteins. These findings show that in B. mori there is a complicated pattern of dsx splicing, and that the regulation of splicing and sex-specific functions of lepidopteran dsx have evolved complexity.

Novel female-specific trans-spliced and alternative splice forms of dsx in the silkworm Bombyx mori.

The Bombyx mori doublesex gene (Bmdsx) plays an important role in somatic sexual development. Its pre-mRNA splices in a sex-specific manner to generate two female-specific and one male-specific splice forms. The present study investigated six novel dsx variants generated by trans-splicing between female dsx transcripts and two additional novel genes, dsr1 and dsr2. Expression analysis indicated that Bmdsx-dsr1 represented splicing noise, whereas dsr2, which trans-spliced with dsx to generate five variants, regulated the expression of the female-specific B. mori dsx transcript Bmdsx(F)s. We unexpectedly found a novel exon 2n insertion during Bmdsx transcription, which did not influence the validity of the novel protein, BmDSX(F3). Ectopic expression of BmDSX(F3) repressed the pheromone-binding protein gene and the testis-specific gene A2 in males, and activated of the storage protein 1 gene. Our findings suggest that trans-splicing is a novel regulatory function of Bmdsx, which participates in female sexual development by regulating the expression of three BmDSX(F) proteins.

Structural insights into small RNA sorting and mRNA target binding by Arabidopsis Argonaute Mid domains.

The RISC-associated Argonaute (Ago) proteins play the catalytic role for RISC-mediated gene regulation by selecting small RNAs and subsequent targeting and cleavage of complementary mRNAs. Ago Mid domains are proposed to play essential roles in small RNA sorting. Here, we report the crystal structures of Arabidopsis Ago1 Mid domain and its chimera mutant with part of Ago1 replaced by Ago4. The structures demonstrate that a single amino insertion in the nucleotide specificity loop of AtAgo1 will change the nucleotide binding preference of AtAgo1 from "5'-U" to "5'-A". Moreover, we identify a long positively charged groove located along the "5'-end-nucleotide specificity loop" and occupied by several sulfate ions with the distance of 9-11Å distance, indicating a putative mRNA target binding groove.

Proteomic analysis of larval integument, trachea and adult scale from the silkworm, Bombyx mori.

Multidimensional LC-tandem MS was used to investigate the protein compositions of three tissues of silkworm, Bombyx mori. A total of 162, 259, and 175 peptides from silkworm larval integument and trachea, and adult scale obtained by database search were matched to 48, 51, and 40 proteins, respectively. Forty-one cuticular proteins were identified from three tissues and covered all five cuticular protein families of silkworm. In the adult scale, all seven cuticular proteins were identified for the first time in the final pellet after SDS extraction. The majority of cuticular proteins were found in each tissue differentially, suggesting that tissue-specific cuticular proteins were involved in the building of the specialized tissues. Seventy-three non-cuticular proteins were also identified in this analysis mainly including muscular proteins, proteinases, inhibitors, transport proteins, and redox-related proteins.