Next Article in Journal
PD-L1 Expression in Mastocytosis
Next Article in Special Issue
The Direct Interaction between E93 and Kr-h1 Mediated Their Antagonistic Effect on Ovary Development of the Brown Planthopper
Previous Article in Journal
OsbZIP81, A Homologue of Arabidopsis VIP1, May Positively Regulate JA Levels by Directly Targetting the Genes in JA Signaling and Metabolism Pathway in Rice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
State Key Laboratory of Silkworm Genome Biology, Biological Science Research Center, Southwest University, Beibei, Chongqing 400715, China
2
Chongqing Key Laboratory of Sericultural Science, Southwest University, Chongqing 400715, China
3
Chongqing Engineering and Technology Research Center for Novel Silk Materials, Southwest University, Chongqing 400715, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(9), 2361; https://doi.org/10.3390/ijms20092361
Submission received: 15 April 2019 / Revised: 10 May 2019 / Accepted: 10 May 2019 / Published: 13 May 2019
(This article belongs to the Special Issue Molecular Ecology, Physiology and Biochemistry of Insects)

Abstract

:
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.

1. Introduction

Sex determination is a fundamental process but has been found to vary among different taxa. Sex determination mechanisms in the silkworm Bombyx mori differ from those in Drosophila melanogaster, especially the upstream of dsx which functions as a double-switch gene in the final step of the sex-determination cascade. In 2001, a double-switch gene, Bmdsx, was found at the final step in the sex-determination cascade of B. mori [1,2]. When binding to a particular sequence splicing repressor, the female-specific splice sites in Bmdsx pre-mRNA are inhibited, leading to the male-specific alternative splicing of Bmdsx [3,4]. A series of mutation analyses using an in vivo splicing assay system identified three distinct sequences (CE1, CE2, and CE3) in exon 4 to be exonic splicing silencers responsible for male-specific splicing of Bmdsx. The CE1 sequence was shown to bind to the nuclear protein BmPSI. Using UAA repeats to replace CE1 inhibited the binding of BmPSI to the RNA and caused female-specific splicing of Bmdsx in male cells [5]. The transgenic knockout strain of BmPSI had a defective gonad in males [6]. Then, the interaction between BmPSI and CE1 was shown to be essential for sex determination of B. mori by silencing the female-specific alternative splicing of Bmdsx. Hence, research on binding characteristics will contribute to the discovery of the mechanisms underlying these alternative splicing processes during sex determination in many species.
BmPSI has four KH_1 motifs. The hnRNP K-homology (KH) motif has been shown to participate in sequence-specific recognition of target RNA in many species [7,8]. For example, the human PSI homolog KSRP contains four KH domains and participates in the alternative splicing of neuron-specific c-src pre-mRNA [9,10]. Furthermore, Drosophila PSI regulates the splicing of P-element transposase pre-mRNA by binding a pseudo-splice site upstream of the authentic splice site using four KH motifs [11,12,13]. Though there have been a few previous studies regarding these special proteins with multiple KH motifs, there is still little known about the detail contribution of the KH_1 motif and the key amino acids, except for the Ile-Gly-X2-Gly-X2-Ile structure, which is conservative in the KH_1 domain [14,15]. BmPSI contains four KH_1 motifs for binding RNA and there is a certain interaction between BmPSI and CE1. Thus, it is a good model protein in which to explore the contribution of the KH motif to recognize and bind RNA.

2. Results

2.1. Four KH_1 Motifs Play Important Roles in the Interaction Between BmPSI and CE1, Especially the KH_1-1 Motif and KH_1-2 Motif

According to the sequence of BmPSI, the tertiary structure of BmPSI was predicted to contain four KH_1 motifs and two other motifs of unknown function. The KH_1 motif was found to be an RNA binding motif in vitro. Therefore, we speculated that these four KH_1 motifs play important roles in the interaction between BmPSI and CE1. We purified maltose binding protein (MBP), wild-type BmPSI, and four truncate BmPSI proteins that respectively truncated KH_1-1, KH_1-2, KH_1-3, and KH_1-4 (Figure 1a). The concentrations of these proteins were adjusted to be identical to that of a bicinchoninic acid assay (BCA), and final concentrations of the adjustment were confirmed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) (Figure 1b). After the concentrations of these six proteins were adjusted uniformly, they were used for an electrophoretic mobility shift assay (EMSA) (Figure 1c). As shown by the EMSA, MBP cannot bind to CE1, and the four mutant proteins almost lost their ability to interact with CE1, especially the ΔKH_1-1 and ΔKH_1-2, compared with the wild-type BmPSI, which was used as a positive control. In order to eliminate the effect of nonspecific binding, wild-type BmPSI, MBP- ΔKH_1-3, and MBP- ΔKH_1-4 proteins were used in a cold competition EMSA experiment (Figure 1d). When the concentration of the competitive probe was aggrandized, the signal intensity of the complex decreased sharply. These results showed that the interaction between CE1 and these three proteins are specific. Hence, we deduced that these four KH_1 motifs participate in the interaction between BmPSI and CE1. In addition, there was absolutely no binding signal found in the lines of MBP- ΔKH_1-1 and MBP- ΔKH_1-2, which was in contrast to the stable binding signal in the lines of MBP- ΔKH_1-3 and MBP- ΔKH_1-4. This result indicates that the KH_1-1 and KH_1-2 motifs play more important roles than the other two KH_1 motifs in the binding affinity of BmPSI.

2.2. Important Amino Acids in the KH_1-1 Motif of BmPSI Found in the Interaction between BmPSI and CE1

We have discussed how KH_1-1 and KH_1-2 play important roles in the interaction between BmPSI and CE1. Through the comparison of the four KH_1 motifs of BmPSI, seven potential key amino acids were predicted (Figure 2a). After cloning these amino acid mutant BmPSI proteins, the proteins were purified with the wild-type BmPSI for EMSA (Figure 2b). EMSA results indicated that five types of mutant proteins lost most of their binding ability compared with the wild-type BmPSI, especially I116G, L127G, and the mut IGGI.
To determine whether these three amino acid mutant proteins transform the secondary structure of BmPSI, circular dichroism (CD) spectroscopy was performed on the wild-type BmPSI and these three mutant proteins (Figure 3a). The results suggest that the three mutant proteins and the wild-type BmPSI have similar patterns. There are two minima at 208 and 222 nm in the result, which is considered to be the norm for α-helices [16]. This suggests that these amino acid mutations did not influence the interaction process by altering the secondary structure of BmPSI. To verify that these three amino acid mutants did lose their ability to bind CE1, isothermal titration calorimetry (ITC) was performed (Figure 3b) [17,18]. Firstly, the wild-type BmPSI was analyzed to explore the appropriate conditions for the experiment. Then, the detailed energetic profile of the interaction between BmPSI and CE1 was obtained (Table 1). The average K value of 1.74 × 107 (n =1.02) of the CE1/BmPSI system indicates that there is a strong interaction between BmPSI and CE1.
The three mutant proteins were analyzed individually with the same conditions grouped from the experiment of BmPSI and CE1. From the ITC results, we found that the ITC data of these three amino acid mutant proteins could not fit the model properly, indicating that these mutant proteins did lose their ability to interact with CE1.

2.3. The Combination of PSI and CE1 is Conserved in Spodoptera Litura

Among the three kinds of amino acid mutants, the two which were not IGGI were found close to the Ile-Gly-X2-Gly-X2-Ile structure, which was conserved in the KH_1 motifs. Hence, we considered that the two amino acids near the Ile-Gly-X2-Gly-X2-Ile structure strongly influenced the combination of PSI and RNA in not only B. mori but also other species. A phylogenic tree was performed with these sequences of homolog PSI in ten species (Figure 4a). It was shown that BmPSI, SlPSI, and Helicoverpa armigera PSI (HaPSI) were clustered into the same clade. From the conservative analysis of KH_1 motifs in these ten species, it was found that the KH_1-1 motif of SlPSI was similar to the KH_1-1 of BmPSI (Figure 4b). SlPSI was selected as the experimental subject to verify the importance of key amino acids for the combination of PSI and CE1. It was confirmed that the CE1 exists conservatively in the female-specific splicing of Sldsx in S. litura. Full length coding sequence (CDS) of SlPSI was cloned and the SlPSI was purified with Ni-NTA and gel filtration chromatography. Cold competition EMSA and SlPSI concentration gradient EMSA were used to show that the combination of PSI and CE1 was also conserved in S. litura (Figure 4c).

2.4. Identifying a Key Amino Acid for the Combination of SlPSI and CE1

Three types of amino acid mutations (L116I, L127G, and IGGI) were added into SlPSI. The mutant proteins were purified and the concentrations of mutant proteins were adjusted consistently for EMSA (Figure 5a). The results showed that mut IGGI can weakly bind to CE1, and that I116G lost most of its ability to interact with CE1 compared to SlPSI. The L127G also lost a lot of its binding capacity but this was not as conspicuous as with I116G. In order to eliminate the effect of nonspecific binding, the SlPSI mutL127G protein was used in a cold competition EMSA experiment (Figure 5b). When the concentration of the competitive probe was aggrandized, the signal intensity of the complex decreased sharply. This result showed that the interaction between CE1 and SlPSI mutL127G is specific. To summarize, we found two key amino acids in the KH-1-1 motif of SlPSI. One was the conservative key structure Ile-Gly-X2-Gly-X2-Ile of the KH_1 motif that has been reported previously [15] and the other was the amino acid I116, which is close to Ile-Gly-X2-Gly-X2-Ile. This is a new discovery of the key amino acid in proteins containing the KH_1 motif for binding its target RNA.

3. Discussion

BmPSI contains four KH_1 motifs and it is considered to be an informative model for exploring how these KH motifs contribute to RNA binding and recognition in these proteins with multiple KH motifs. The EMSA of BmPSI and four truncate BmPSI proteins showed that the four KH_1 motifs participate in the interaction between BmPSI and CE1, in which the KH_1-1 and KH_1-2 play more important roles in the combination process. Through sequence alignment, EMSA, CD spectroscopy, and ITC, we found three potentially key amino acids in the KH_1-1 motif of BmPSI. After verifying the interaction of SlPSI and CE1, these three kinds of amino acid mutations were cloned into SlPSI. Except for the L127G mutation, we found that the I116G and IGGI mutations lost most of their binding ability compared with the wild-type SlPSI.
Drosophila PSI requires multiple tandem KH domains for specific and high-affinity recognition of substrate RNA [11]; however, the KH_1-1 and KH_1-2 motifs of BmPSI play more important roles in the interaction of BmPSI and Bmdsx pre-mRNA compared with the other two KH_1 motifs. Differences in the binding ability of the L127 key site in BmPSI and SlPSI suggest that there are a lot of differences between these two relatively close lepidopterans. The binding affinity of mutant I116G sharply decreased in BmPSI and SlPSI, indicating that this is a key site in the KH_1 motif for the binding process. Furthermore, I116G is stably conserved in Lepidoptera, which indicates it is a key amino acid for binding affinity in proteins containing KH_1 motifs, just as EMSA confirmed it was. I116 is isoleucine, which is tightly close to the Ile-Gly-X2-Gly-X2-Ile structure. These discoveries provide further insight into how KH_1 motifs contribute to the combining of a protein and its target RNA sequence.
The isoforms of sex-specific alternative splicing of Bmdsx are translated into sex-specific BmDSX proteins which lead to differences in the expression of sex-specific genes, and, thus, sexual differentiation [19,20]. BmPSI is the most important part of the exonic splicing silencer complex for binding Bmdsx pre-mRNA with CE1 to promote male-specific splicing of Bmdsx [21,22]. This is the only mechanism known as being closely upstream of Bmdsx in the sex determination cascade of B. mori. The researchers could not generate a gender reversal strain through genetic modification of Bmdsx or the knockout of Bmpsi. Two possible reasons for this are that Bmdsx leads the expressions of many sex-specific genes and Bmpsi participates in not only sex determination but also other key pathways in B. mori. Therefore, a slight change in Bmdsx or Bmpsi would cause disorder in transgene stains. Thus, the key amino acid sites of the interaction between BmPSI and CE1 could be used to establish transgenic strains, in which the transgenic BmPSI loses its capacity to bind to CE1 but retains its other functions.
From an evolutionary perspective, the four conserved KH_1 motifs in SlPSI and CE1 in female specific splicing of Sldsx provide strong evidence of conserved structures for essential components of the sex determination process in Lepidoptera [23]. Furthermore, the conserved interaction of SlPSI and CE1 indicates functional conservation of essential components in the sex determination process in Lepidoptera. In addition, the discovery of this key site for the binding in both BmPSI and SlPSI not only contributes a lot to future research on sex determination of lepidopterans but also provides a theoretical basis for gender reversal in Bombyx mori [24,25].

4. Materials and Methods

4.1. Construction of Recombinant Expression Vectors

The pET His6 MBP TEV LIC cloning vector is an LIC N-terminal fusion vector for E. coli expression which contains a maltose binding protein. The wild type Bmpsi sequence was cloned into this vector behind a (His)6 affinity tag by LIC cloning. Four different KH_1 region truncations of Bmpsi were made with overlapping PCR. Seven point mutations in the KH_1-1 motif of Bmpsi were added with site-directed mutagenesis into the recombinant vector containing Bmpsi (Table 2).
The sequence of Slpsi was downloaded from NCBI and was predicted to be the S. litura far upstream element-binding protein (LOC 111350908). According to the sequence information on NCBI, we designed a pair of primers to clone the fragment of Slpsi into the pET His6 MBP TEV LIC cloning vector. Three kinds of point mutation in the KH_1-1 motif of Slpsi were added with site-directed mutagenesis into the recombinant vector (Table 3).

4.2. Preparation of Single Stranded RNA

Two kinds of CE1 single-stranded RNA were prepared for EMSA and ITC. One kind of CE1 RNA was labeled with Biotin on its 5′ end and 3′ end. The sequence of the single-stranded RNA was the same as the sequence of CE1 in the female-specific splicing of Bmdsx (5′-uuaauaauauaaguggugua-3′). The RNA probe was compounded by the Beijing Genomics Institute and purified with ion exchange HPLC. The RNA was dissolved in a buffer (20 mM Tris, 20 mM NaCl), which was used to purify the different kinds of protein.

4.3. The Overexpression and Purification of Different Kinds of Proteins

Transetta (DE3) Chemically Competent Cell was used to overexpress the wild-type BmPSI, the wild-type SlPSI, the four KH_1 region truncation proteins of BmPSI, the seven amino acid mutants of BmPSI proteins, and the three amino acid mutants of SlPSI proteins with the recombinant vectors discussed above. We overexpressed the maltose-binding protein with the original expression vector as a negative control. The purification of every kind of protein needed 2 L cells in LB media at 37 °C with an absorbance at 600 nm of 0.5. Then, we added Isopropyl β-D-Thiogalactoside (IPTG)to a final concentration of 1 mM and cultured the cells at 37 °C for another 4 h. After cell culturing, the cells were gathered and crushed with multi-gelation three times and ultra-sonicated (40%, 30 min). The supernatant was separated using centrifugation (15,000× g, 10 min). A series of buffers with different concentrations of imidazole were used to wash the HiTrap Chelating where the supernatant was loaded. The buffer containing the target protein was passed through gel filtration chromatography, which made it purer. Finally, we retrieved the pure target protein with the (His)6 -MBP-TEV component.

4.4. Bioinformatics Analysis

The sequences of the four KH_1 motifs in BmPSI were compared to find out the potential key amino acids which are conservative in KH_1-1and KH_1-2 but not KH_1-3 and KH_1-4 [26,27,28]. The PSI sequences of ten species were downloaded from NCBI to construct a phylogenic tree [29,30] and the sequences of the four KH_1 motifs from BmPSI were compared by multiple sequence alignment.

4.5. Electrophoretic Mobility Shift Assay

Electrophoretic mobility shift assay can be used to determine binding affinity, specificity, and stoichiometry of the RNA/protein interaction [31,32,33]. In this article, we used EMSA experiments in different ways which were designed to accommodate two different experimental objectives. A normal EMSA was used to verify the interaction between SlPSI and CE1 with cold competition EMSA and SlPSI concentration gradient EMSA. The other kind of EMSA was designed to compare the binding affinity with RNA of different proteins under the same conditions. At first, the concentrations of proteins in every group were adjusted consistently with BCA assay and SDS-PAGE. Then, the proteins were experimented on under exactly the same conditions. As a result, the signal intensity of the complex stood for the abilities of different proteins to interact with the RNA probe. One group of proteins included MBP, MBP-BmPSI, and four KH_1 region truncate BmPSI proteins, which were set to compare the binding abilities of different KH_1 motifs. Another group of proteins included MBP-BmPSI and the seven amino acid mutant BmPSI proteins, which were set to show the importance these key amino acids had in the interaction between BmPSI and CE1. The third group of proteins included MBP-SlPSI and the three amino acid mutant SlPSI proteins, which were set to detect the conservative property of the key amino acids we found in BmPSI for the combination of PSI and CE1 in other species.

4.6. CD Wavelength Scans

Circular dichroism spectroscopy is a technique which is widely used to detect secondary structure composition. The wild-type BmPSI protein and the three kinds of amino acid mutations of BmPSI, which were important to the interaction between BmPSI and CE1, were experimented on for these secondary structures with CD wavelength scans. Before the CD wavelength test, the concentrations of different proteins were adjusted consistently with a BCA protein assay kit and SDS-PAGE.

4.7. Isothermal Titration Calorimetry

Isothermal titration calorimetry is always used to estimate binding ability through the binding enthalpy (ΔH) [34]. An ITC isotherm of the CE1/BmPSI system was performed at 25 °C with an cell volume of 1.4 mL using a stirring rate of 307 rpm. The CE1 probe was injected 2 μL at a time into a solution of BmPSI at 210 s to allow complete equilibration of the CE1/BmPSI system between each injection. The most appropriate experimental conditions were grouped with the experiment between wild-type BmPSI and the CE1 probe. Then, the other amino acid mutants with the same concentration were experimented on under the same conditions. The differences in binding ability among these four proteins were determined from the curves constructed during these ITC experiments.

Author Contributions

Data curation, Y.W.; formal analysis, Y.W.; funding acquisition, X.-F.Z.; investigation, X.-F.Z.; methodology, Y.W. and X.-F.Z.; software, Y.W.; validation, Y.W.; writing—original draft, Y.W.; writing—review and editing, Q.Z., Q.-X.W., K.-X.W., and X.-F.Z.

Funding

This research was funded by the State Key Program of National Natural Science of China (grant no. 31530071) and the National Natural Science Foundation of China (grant no. 31272502).

Acknowledgments

We thank Lin Liu (South China Normal University) for providing the sequence of the female-specific splicing of Sldsx.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

dsxDoublesex
PSIP-element somatic inhibitor
CDCircular dichroism
ITCIsothermal titration calorimetry
EMSAElectrophoretic mobility shift assay

References

  1. Ohbayashi, F.; Suzuki, M.G.; Mita, K.; Okano, K.; Shimada, T. A homologue of the Drosophila doublesex gene is transcribed into sex-specific mRNA isoforms in the silkworm, Bombyx mori. Comp. Biochem. Phys. B 2001, 128, 145–158. [Google Scholar] [CrossRef]
  2. Suzuki, M.G.; Funaguma, S.; Kanda, T.; Tamura, T.; Shimada, T. Role of the male BmDSX protein in the sexual differentiation of Bombyx mori. Evol. Dev. 2005, 7, 58–68. [Google Scholar] [CrossRef]
  3. Funaguma, S.; Suzuki, M.G.; Tamura, T.; Shimada, T. The Bmdsx transgene including trimmed introns is sex-specifically spliced in tissues of the silkworm, Bombyx mori. J. Insect Sci. 2005, 5, 17–22. [Google Scholar] [CrossRef]
  4. Suzuki, M.G.; Ohbayashi, F.; Mita, K.; Shimada, T. The mechanism of sex-specific splicing at the doublesex gene is different between Drosophila melanogaster and Bombyx mori. Insect Biochem. Mol Biol. 2001, 31, 1201–1211. [Google Scholar] [CrossRef]
  5. Suzuki, M.G.; Imanishi, S.; Dohmae, N.; Nishimura, T.; Shimada, T.; Matsumoto, S. Establishment of a novel in vivo sex-specific splicing assay system to identify a trans-acting factor that negatively regulates splicing of Bombyx mori dsx female exons. J. Mol. Cell Biol. 2008, 28, 333–343. [Google Scholar] [CrossRef]
  6. Xu, J.; Chen, S.; Zeng, B.; James, A.A.; Tan, A.; Huang, Y. Bombyx mori P-element somatic inhibitor (BmPSI) is a key auxiliary factor for silkworm male sex determination. PLoS Genet. 2017, 13, e1006576. [Google Scholar] [CrossRef] [PubMed]
  7. Adinolfi, S.; Bagni, C.; Castiglione Morelli, M.A.; Fraternali, F.; Musco, G.; Pastore, A. Novel RNA-binding motif: The KH module. Biopolymers 1999, 51, 153–164. [Google Scholar] [CrossRef]
  8. Fraternali, F.; Amodeo, P.; Musco, G.; Nilges, M.; Pastore, A. Exploring protein interiors: The role of a buried histidine in the KH module fold. Proteins 1999, 34, 484–496. [Google Scholar] [CrossRef] [Green Version]
  9. Hall, M.P.; Huang, S.; Black, D.L. Differentiation-induced colocalization of the KH-type splicing regulatory protein with polypyrimidine tract binding protein and the c-src pre-mRNA. Mol. Biol. Cell 2004, 15, 774–786. [Google Scholar] [CrossRef] [PubMed]
  10. Art, J.; Besche, V.; Bros, M.; Li, H.; Handler, N.; Bauer, F.; Erker, T.; Behnke, F.; Monch, B.; Forstermann, U.; et al. Identification of the KH type splicing regulatory protein (KSRP) as a new important mediator of the anti-inflammatory effects of resveratrol. N.-S. Arch. Pharmacol. 2012, 385, 6. [Google Scholar]
  11. Chmiel, N.H.; Rio, D.C.; Doudna, J.A. Distinct contributions of KH domains to substrate binding affinity of Drosophila P-element somatic inhibitor protein. RNA 2006, 12, 283–291. [Google Scholar] [CrossRef]
  12. Guo, L.; Zaysteva, O.; Nie, Z.; Mitchell, N.C.; Amanda Lee, J.E.; Ware, T.; Parsons, L.; Luwor, R.; Poortinga, G.; Hannan, R.D.; et al. Defining the essential function of FBP/KSRP proteins: Drosophila Psi interacts with the mediator complex to modulate MYC transcription and tissue growth. Nucleic Acids Res. 2016, 44, 7646–7658. [Google Scholar] [CrossRef]
  13. Labourier, E.; Blanchette, M.; Feiger, J.W.; Adams, M.D.; Rio, D.C. The KH-type RNA-binding protein PSI is required for Drosophila viability, male fertility, and cellular mRNA processing. Gene Dev. 2002, 16, 72–84. [Google Scholar] [CrossRef]
  14. Makeyev, A.V.; Eastmond, D.L.; Liebhaber, S.A. Targeting a KH-domain protein with RNA decoys. RNA 2002, 8, 1160–1173. [Google Scholar] [CrossRef] [Green Version]
  15. Xu, H.X.; Xu, Y.J.; Liang, X.J.; Zeng, W.X. The biological function of heterogeneous nuclear ribonucleoprotein K (hnRNP K) and its roles in spermatogenisis. J. Agric. Biotechnol. 2015, 23, 661–670. [Google Scholar]
  16. Micsonai, A.; Wien, F.; Kernya, L.; Lee, Y.H.; Goto, Y.; Refregiers, M.; Kardos, J. Accurate secondary structure prediction and fold recognition for circular dichroism spectroscopy. Proc. Natl. Acad. Sci. USA 2015, 112, E3095–E3103. [Google Scholar] [CrossRef] [Green Version]
  17. Arroyo-Maya, I.J.; McClements, D.J. Application of ITC in foods: A powerful tool for understanding the gastrointestinal fate of lipophilic compounds. BBA-Gen. Subjects 2016, 1860, 1026–1035. [Google Scholar] [CrossRef]
  18. Johnson, R.A.; Manley, O.M.; Spuches, A.M.; Grossoehme, N.E. Dissecting ITC data of metal ions binding to ligands and proteins. BBA-Gen. Subjects 2016, 1860, 892–901. [Google Scholar] [CrossRef]
  19. Kawaoka, S.; Kadota, K.; Arai, Y.; Suzuki, Y.; Fujii, T.; Abe, H.; Yasukochi, Y.; Mita, K.; Sugano, S.; Shimizu, K.; et al. The silkworm W chromosome is a source of female-enriched piRNAs. RNA 2011, 17, 2144–2151. [Google Scholar] [CrossRef] [Green Version]
  20. Hara, K.; Fujii, T.; Suzuki, Y.; Sugano, S.; Shimada, T.; Katsuma, S.; Kawaoka, S. Altered expression of testis-specific genes, piRNAs, and transposons in the silkworm ovary masculinized by a W chromosome mutation. BMC Genomics 2012, 13, 119–128. [Google Scholar] [CrossRef]
  21. Suzuki, M.G.; Imanishi, S.; Dohmae, N.; Asanuma, M.; Matsumoto, S. Identification of a male-specific RNA binding protein that regulates sex-specific splicing of Bmdsx by increasing RNA binding activity of BmPSI. Mol. Cell Biol. 2010, 30, 5776–5786. [Google Scholar] [CrossRef]
  22. Kiuchi, T.; Koga, H.; Kawamoto, M.; Shoji, K.; Sakai, H.; Arai, Y.; Ishihara, G.; Kawaoka, S.; Sugano, S.; Shimada, T.; et al. A single female-specific piRNA is the primary determiner of sex in the silkworm. Nature 2014, 509, 633–636. [Google Scholar] [CrossRef]
  23. Du, Q.; Wen, L.; Zheng, S.C.; Bi, H.L.; Huang, Y.P.; Feng, Q.L.; Liu, L. Identification and functional characterization of doublesex gene in the testis of Spodoptera litura. Insect Sci. 2018. [Google Scholar] [CrossRef]
  24. Nagashim, E.; Hujieda, T.; Moriizum, S. On expression of several characters affected by difference of sexes in silkworm. Jpn. J. Genet. 1967, 42, 427. [Google Scholar]
  25. Shabdini, A.; Seidavi, A.R.; Gharahveysi, S.; Mirhosseini, S.Z. Effects of sex, generation and genetic group on economic traits in three Iranian commercial silkworm pure lines. J. Anim. Vet. Adv. 2011, 10, 76–80. [Google Scholar] [CrossRef]
  26. Amorim, A.R.; Neves, L.A.; Valencio, C.R.; Roberto, G.F.; Zafalon, G.F.D. An approach for COFFEE objective function to global DNA multiple sequence alignment. Comput. Biol. Chem. 2018, 75, 39–44. [Google Scholar] [CrossRef]
  27. Nute, M.; Saleh, E.; Warnow, T. Evaluating statistical multiple sequence alignment in comparison to other alignment methods on protein data sets. Syst. Biol. 2018, 68, 396–411. [Google Scholar] [CrossRef] [PubMed]
  28. Dong, D.; Su, W.; Shi, W.; Zou, Q.; Peng, S. VCSRA: A fast and accurate multiple sequence alignment algorithm with a high degree of parallelism. J. Genet. Genomics 2018, 45, 407–410. [Google Scholar] [CrossRef]
  29. Mirkov, T.E.; Wahlstrom, J.M.; Hagiwara, K.; Finardi-Filho, F.; Kjemtrup, S.; Chrispeels, M.J. Evolutionary relationships among proteins in the phytohemagglutinin-arcelin-alpha-amylase inhibitor family of the common bean and its relatives. Plant. Mol. Biol. 1994, 26, 1103–1113. [Google Scholar] [CrossRef]
  30. Lai, X.; Zhou, Y.; Xia, X.; Zhang, Q. Performance analysis of evolutionary algorithms for steiner tree problems. Evol. Comput. 2017, 25, 707–723. [Google Scholar] [CrossRef]
  31. Hellman, L.M.; Fried, M.G. Electrophoretic mobility shift assay (EMSA) for detecting protein-nucleic acid interactions. Nat. Protoc. 2007, 2, 1849–1861. [Google Scholar] [CrossRef] [PubMed]
  32. Jing, D.; Agnew, J.; Patton, W.F.; Hendrickson, J.; Beechem, J.M. A sensitive two-color electrophoretic mobility shift assay for detecting both nucleic acids and protein in gels. Proteomics 2003, 3, 1172–1180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fillebeen, C.; Wilkinson, N.; Pantopoulos, K. Electrophoretic mobility shift assay (EMSA) for the study of RNA-protein interactions: The IRE/IRP example. J. Vis. Exp. 2014, 94, e52230. [Google Scholar] [CrossRef] [PubMed]
  34. Jensen, L.B.; Mortensen, K.; Pavan, G.M.; Kasimova, M.R.; Jensen, D.K.; Gadzhyeva, V.; Nielsen, H.M.; Foged, C. Molecular characterization of the interaction between siRNA and PAMAM G7 dendrimers by SAXS, ITC, and molecular dynamics simulations. Biomacromolecules 2010, 11, 3571–3577. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Structure and binding ability of BmPSI and truncation proteins: (a) Comparison among BmPSI and four truncation proteins used for the binding experiment. “DUF” means “domain of unknown function”. (b) The dosage of these proteins. “NC” means “negative control”. (c) Binding ability of these proteins with CE1. Buffer and maltose binding protein (MBP) were set as the negative control and wild-type BmPSI was used as the positive control. (d) Cold competition electrophoretic mobility shift assay (EMSA) experiment between CE1 and BmPSI, MBP-ΔKH_1-3, and MBP-ΔKH_1-4 proteins.
Figure 1. Structure and binding ability of BmPSI and truncation proteins: (a) Comparison among BmPSI and four truncation proteins used for the binding experiment. “DUF” means “domain of unknown function”. (b) The dosage of these proteins. “NC” means “negative control”. (c) Binding ability of these proteins with CE1. Buffer and maltose binding protein (MBP) were set as the negative control and wild-type BmPSI was used as the positive control. (d) Cold competition electrophoretic mobility shift assay (EMSA) experiment between CE1 and BmPSI, MBP-ΔKH_1-3, and MBP-ΔKH_1-4 proteins.
Ijms 20 02361 g001
Figure 2. Important amino acids in the KH_1-1 motif of BmPSI: (a) Comparison among amino acid sequences of the KH_1 motifs in BmPSI. The black color represents that amino acids are completely identical. The gray color represents that partial amino acids are identical. Different colored boxes represent different potential key amino acid sites for the binding. (b) After the concentrations were adjusted uniformly, EMSA was performed on wild-type BmPSI and seven amino acid mutants of BmPSI.
Figure 2. Important amino acids in the KH_1-1 motif of BmPSI: (a) Comparison among amino acid sequences of the KH_1 motifs in BmPSI. The black color represents that amino acids are completely identical. The gray color represents that partial amino acids are identical. Different colored boxes represent different potential key amino acid sites for the binding. (b) After the concentrations were adjusted uniformly, EMSA was performed on wild-type BmPSI and seven amino acid mutants of BmPSI.
Ijms 20 02361 g002
Figure 3. Two amino acids in the KH_1-1 motif of BmPSI play important roles in the binding with CE1: (a) The secondary structures of wild-type BmPSI and three kind of amino acid mutant proteins detected with circular dichroism (CD) spectroscopy. (b) The RNA binding affinity of three kinds of amino acid mutant proteins and wild-type BmPSI were assessed, showing different exports of isothermal titration calorimetry (ITC) under the same conditions.
Figure 3. Two amino acids in the KH_1-1 motif of BmPSI play important roles in the binding with CE1: (a) The secondary structures of wild-type BmPSI and three kind of amino acid mutant proteins detected with circular dichroism (CD) spectroscopy. (b) The RNA binding affinity of three kinds of amino acid mutant proteins and wild-type BmPSI were assessed, showing different exports of isothermal titration calorimetry (ITC) under the same conditions.
Ijms 20 02361 g003
Figure 4. The combination of PSI and CE1 is conserved in Spodoptera litura: (a) A phylogenic tree based on the PSI of 10 species. (b) Conservation analysis of KH_1 motifs in these 10 species. The black color represents that amino acids are completely identical. The gray color represents that partial amino acids are identical. (c) SlPSI concentration gradient EMSA and cold competition EMSA.
Figure 4. The combination of PSI and CE1 is conserved in Spodoptera litura: (a) A phylogenic tree based on the PSI of 10 species. (b) Conservation analysis of KH_1 motifs in these 10 species. The black color represents that amino acids are completely identical. The gray color represents that partial amino acids are identical. (c) SlPSI concentration gradient EMSA and cold competition EMSA.
Ijms 20 02361 g004
Figure 5. EMSA of wild-type SlPSI and three kinds of potentially key amino acid mutant SlPSI proteins. (a) MBP-BmPSI was set as the positive control. The signal intensity of the complex stood for the abilities of different proteins to interact with the CE1 probe. (b) Cold competition EMSA experiment between CE1 and SlPSI mutL127G.
Figure 5. EMSA of wild-type SlPSI and three kinds of potentially key amino acid mutant SlPSI proteins. (a) MBP-BmPSI was set as the positive control. The signal intensity of the complex stood for the abilities of different proteins to interact with the CE1 probe. (b) Cold competition EMSA experiment between CE1 and SlPSI mutL127G.
Ijms 20 02361 g005
Table 1. The energetic profile of the interaction of the BmPSI with CE1 obtained from ITC.
Table 1. The energetic profile of the interaction of the BmPSI with CE1 obtained from ITC.
MoleculeΔH (cal/mol)ΔS (cal/mol/deg)nK (M−1)
BmPSI−3.3 × 104 ± 269.2−77.91.021.74 × 107
Table 2. Primer sequences of modifications on Bmpsi.
Table 2. Primer sequences of modifications on Bmpsi.
NamePrimers
Bmpsi5′-cgcggatccatgagtgattattcttctat-3′,
5′-ccgctcgagtcactgctggtggtcggagc-3′
ΔKH_1-15′-cgatcagagcatcaatattgtgaatcatcgaggaa-3′
5′-ttcctcgatgattcacaatattgatgctctgatcg-3′
ΔKH_1-25′-ccaggccctaatgcaatgctcctcgccaacaagga-3′
5′-tccttgttggcgaggagcattgcattagggcctgg-3′
ΔKH_1-35′-aacggactcgccaccactcttatatctagtgtcaa-3′
5′-ttgacactagatataagagtggtggcgagtccgtt-3′
ΔKH_1-45′-gaccgaccagagatgcgcaaagttggtgggcctgt-3′
5′-acaggcccaccaactttgcgcatctctggtcggtc-3′
I116G5′-gttggactaggaattggacgtggt-3′
5′-accacgtccaattcctagtccaac-3′
L127G5′-atcaccagagggcaagcagaatcc-3′
5′-ggattctgcttgccctctggtgat-3′
K134G,5′-tccggttgcgggatacaaatggca-3′
5′-tgccatttgtatcccgcaaccgga-3′
A161G5′-atacagagaggtaaagaattagtg-3′
5′-cactaattctttacctctctgtat-3′
L164G5′-gctaaagaaggagtgaaccaaatt-3′
5′-aatttggttcactccttctttagc-3′
V165G5′-aaagaattagggaaccaaattgtg-3′
5′-cacaatttggttccctaattcttt-3′
IGGI 117, 118, 121, 124 AAAA5′-ggactaatagctgcacgtggtgga-3′
5′-tccaccacgtgcagctattagtcc-3′,
5′-cacgtggtgcagaacaagccaccagactgca-3′
5′-tgcagtctggtggcttgttctgcaccacgtg-3′,
5′- gcagaacaagccaccagactgcaa-3′
5′-ttgcagtctggtggcttgttctgc-3′
Table 3. Primer sequences of modifications on Slpsi
Table 3. Primer sequences of modifications on Slpsi
NamePrimers
Slpsi5′- cgggatccatgagtgattattcttctatggc-3′
5′-ggaattcctattgatgatggtcggg-3′
I116G5′-atggttggactcggaattgggcgcggc-3′
5′-gccgcgcccaattccgagtccaaccat-3′
L127G5′-atcacacgcgggcaggcggagtca-3′
5′-tgactccgcctgcccgcgtgtgat-3′
IGGI 117, 118, 121, 124 AAAA5′-ggactcatagctgggcgcggcggtgag-3′
5′-ctcaccgccgcgcccagctatgagtcc-3′,
5′-ggactcatagctgcacgcggcggtgag-3′
5′-ctcaccgccgcgtgcagctatgagtcc-3′,
5′-gcacgcggcgcagagcaggccacacgc-3′
5′-gcgtgtggcctgctctgcgccgcgtgc-3′

Share and Cite

MDPI and ACS Style

Wang, Y.; Zhao, Q.; Wan, Q.-X.; Wang, K.-X.; Zha, X.-F. P-element Somatic Inhibitor Protein Binding a Target Sequence in dsx Pre-mRNA Conserved in Bombyx mori and Spodoptera litura. Int. J. Mol. Sci. 2019, 20, 2361. https://doi.org/10.3390/ijms20092361

AMA Style

Wang Y, Zhao Q, Wan Q-X, Wang K-X, Zha X-F. P-element Somatic Inhibitor Protein Binding a Target Sequence in dsx Pre-mRNA Conserved in Bombyx mori and Spodoptera litura. International Journal of Molecular Sciences. 2019; 20(9):2361. https://doi.org/10.3390/ijms20092361

Chicago/Turabian Style

Wang, Yao, Qin Zhao, Qiu-Xing Wan, Kai-Xuan Wang, and Xing-Fu Zha. 2019. "P-element Somatic Inhibitor Protein Binding a Target Sequence in dsx Pre-mRNA Conserved in Bombyx mori and Spodoptera litura" International Journal of Molecular Sciences 20, no. 9: 2361. https://doi.org/10.3390/ijms20092361

APA Style

Wang, Y., Zhao, Q., Wan, Q. -X., Wang, K. -X., & Zha, X. -F. (2019). P-element Somatic Inhibitor Protein Binding a Target Sequence in dsx Pre-mRNA Conserved in Bombyx mori and Spodoptera litura. International Journal of Molecular Sciences, 20(9), 2361. https://doi.org/10.3390/ijms20092361

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop