The Bombyx ovary-derived cell line endogenously expresses PIWI/PIWI-interacting RNA complexes

Genetic studies and large-scale sequencing experiments have revealed that the PIWI subfamily proteins and PIWI-interacting RNAs (piRNAs) play an important role in germ line development and transposon control. Biochemical studies in vitro have greatly contributed to the understanding of small interfering RNA (siRNA) and microRNA (miRNA) pathways. However, in vitro analyses of the piRNA pathway have been thus far quite challenging, because their expression is largely restricted to the germ line. Here we report that Bombyx mori ovary-derived cultured cell line, BmN4, endogenously expresses two PIWI subfamily proteins, silkworm Piwi (Siwi) and Ago3 (BmAgo3), and piRNAs associated with them. Siwi-bound piRNAs have a strong bias for uridine at their 5′ end and BmAgo3-bound piRNAs are enriched for adenine at position 10. In addition, Siwi preferentially binds antisense piRNAs, whereas BmAgo3 binds sense piRNAs. Moreover, we identified many pairs in which Siwi-bound antisense and BmAgo3-bound sense piRNAs are overlapped by precisely 10 nt at their 5′ ends. These signatures are known to be important for secondary piRNA biogenesis in other organisms. Taken together, BmN4 is a unique cell line in which both primary and secondary steps of piRNA biogenesis pathways are active. This cell line would provide useful tools for analysis of piRNA biogenesis and function.     

Metazoan Maelstrom is an RNA-binding protein that has evolved from an ancient nuclease active in protists

Piwi-interacting RNAs (piRNAs) guide Piwi argonautes to their transposon targets for silencing. The highly conserved protein Maelstrom is linked to both piRNA biogenesis and effector roles in this pathway. One defining feature of Maelstrom is the predicted MAEL domain of unknown molecular function. Here, we present the first crystal structure of the MAEL domain from Bombyx Maelstrom, which reveals a nuclease fold. The overall architecture resembles that found in Mg²⁺- or Mn²⁺-dependent DEDD nucleases, but a clear distinguishing feature is the presence of a structural Zn²⁺ ion coordinated by the conserved ECHC residues. Strikingly, metazoan Maelstrom orthologs across the animal kingdom lack the catalytic DEDD residues, and as we show for Bombyx Maelstrom are inactive as nucleases. However, a MAEL domain-containing protein from amoeba having both sequence motifs (DEDD and ECHC) is robustly active as an exoribonuclease. Finally, we show that the MAEL domain of Bombyx Maelstrom displays a strong affinity for single-stranded RNAs. Our studies suggest that the ancient MAEL nuclease domain evolved to function as an RNA-binding module in metazoan Maelstrom.     

The nuclear receptor gene nhr-25 plays multiple roles in the Caenorhabditis elegans heterochronic gene network to control the larva-to-adult transition

Developmental timing in the nematode Caenorhabditis elegans is controlled by heterochronic genes, mutations in which cause changes in the relative timing of developmental events. One of the heterochronic genes, let-7, encodes a microRNA that is highly evolutionarily conserved, suggesting that similar genetic pathways control developmental timing across phyla. Here we report that the nuclear receptor nhr-25, which belongs to the evolutionarily conserved fushi tarazu-factor 1/nuclear receptor NR5A subfamily, interacts with heterochronic genes that regulate the larva-to-adult transition in C. elegans. We identified nhr-25 as a regulator of apl-1, a homolog of the Alzheimer's amyloid precursor protein-like gene that is downstream of let-7 family microRNAs. NHR-25 controls not only apl-1 expression but also regulates developmental progression in the larva-to-adult transition. NHR-25 negatively regulates the expression of the adult-specific collagen gene col-19 in lateral epidermal seam cells. In contrast, NHR-25 positively regulates the larva-to-adult transition for other timed events in seam cells, such as cell fusion, cell division and alae formation. The genetic relationships between nhr-25 and other heterochronic genes are strikingly varied among several adult developmental events. We propose that nhr-25 has multiple roles in both promoting and inhibiting the C. elegans heterochronic gene pathway controlling adult differentiation programs.     

The influences of PRG-1 on the expression of small RNAs and mRNAs

Background

In metazoans, Piwi-related Argonaute proteins play important roles in maintaining germline integrity and fertility and have been linked to a class of germline-enriched small RNAs termed piRNAs. Caenorhabditis elegans encodes two Piwi family proteins called PRG-1 and PRG-2, and PRG-1 interacts with the C. elegans piRNAs (21U-RNAs). Previous studies found that mutation of prg-1 causes a marked reduction in the expression of 21U-RNAs, temperature-sensitive defects in fertility and other phenotypic defects.

Results

In this study, we wanted to systematically demonstrate the function of PRG-1 in the regulation of small RNAs and their targets. By analyzing small RNAs and mRNAs with and without a mutation in prg-1 during C. elegans development, we demonstrated that (1) mutation of prg-1 leads to a decrease in the expression of 21U-RNAs, and causes 35 ~ 40% of miRNAs to be down-regulated; (2) in C. elegans, approximately 3% (6% in L4) of protein-coding genes are differentially expressed after mutating prg-1, and 60 ~ 70% of these substantially altered protein-coding genes are up-regulated; (3) the target genes of the down-regulated miRNAs and the candidate target genes of the down-regulated 21U-RNAs are enriched in the up-regulated protein-coding genes; and (4) PRG-1 regulates protein-coding genes by down-regulating small RNAs (miRNAs and 21U-RNAs) that target genes that participate in the development of C. elegans.

Conclusions

In prg-1-mutated C. elegans, the expression of miRNAs and 21U-RNAs was reduced, and the protein-coding targets, which were associated with the development of C. elegans, were up-regulated. This may be the mechanism underlying PRG-1 function.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-321) contains supplementary material, which is available to authorized users.     

BmICE-2 is a novel pro-apoptotic caspase involved in apoptosis in the silkworm, Bombyx mori

In this study we identified a potential pro-apoptotic caspase gene, Bombyx mori(B. mori)ICE-2 (BmICE-2) which encoded a polypeptide of 284 amino acid residues, including a ¹⁶⁹QACRG¹⁷³ sequence which surrounded the catalytic site and contained a p20 and a p10 domain. BmICE-2 expressed in Escherichia coli (E. coli) exhibited high proteolytic activity for the synthetic human initiator caspase-9 substrates Ac-LEHD-pNA, but little activity towards the effector caspase-3 substrates Ac-DEVD-pNA. When BmICE-2 was transiently expressed in BmN-SWU1 silkworm B. mori cells, we found that the high proteolytic activity for Ac-LEHD-pNA triggered caspase-3-like protease activity resulting in spontaneous cleavage and apoptosis in these cells. This effect was not replicated in Spodoptera frugiperda 9 cells. In addition, spontaneous cleavage of endogenous BmICE-2 in BmN-SWU1 cells could be induced by actinomycin D. These results suggest that BmICE-2 may be a novel pro-apoptotic gene with caspase-9 activity which is involved apoptotic processes in BmN-SWU1 silkworm B. mori cells.     

Spatial expression of the mevalonate enzymes involved in juvenile hormone biosynthesis in the corpora allata in Bombyx mori

The developmental expressions of the mRNA of JH synthetic enzymes have been studied using homogenates of the corpora cardiaca-corpora allata (CC-CA) complexes in Bombyx mori [Kinjoh, T., Kaneko, Y., Itoyama, K., Mita, K., Hiruma, K., Shinoda, T., 2007. Control of juvenile hormone biosynthesis in Bombyx mori: cloning of the enzymes in the mevalonate pathway and assessment of their developmental expression in the corpora allata. Insect Biochemistry and Molecular Biology 37, 808-818]. The in situ hybridization analyses in the CC-CA complex showed that the distribution of the mRNAs of all the mevalonate enzymes and juvenile hormone (JH) acid O-methyltransferase occurred only in the CA cells, indicating that the fluctuations of the enzyme mRNA amounts in the CC-CA complexes were derived solely from the CA. In addition, the size of the CA and their nuclei was not associated with the JH synthetic activity by the CA until the pharate adult. Only female adult CA synthesized JH in B. mori, and the CA and the nuclei were significantly larger than those of male CA which do not synthesize JH.     

Organization of the Hox gene cluster of the silkworm, <Emphasis Type="Italic">Bombyx mori</Emphasis>: a split of the Hox cluster in a non-<Emphasis Type="Italic">Drosophila</Emphasis> insect

A bacterial artificial chromosome (BAC) contig was constructed by chromosome walking, starting from the Hox genes of the silkworm, Bombyx mori. Bombyx orthologues of the labial (lab) and zerknült (zen) genes were newly identified. The size of the BAC contig containing the Hox gene cluster—except the lab and Hox 2 genes—was estimated to be more than 2 Mb. The Bombyx Hox cluster was mapped to linkage group (LG) 6. The lab gene was mapped on the same LG, but far apart from the cluster. Fluorescence in situ hybridization analysis confirmed that the major Hox gene cluster and lab were at different locations on the same chromosome in B. mori.Edited by M. Akam     

The G protein-coupled receptors in the silkworm,Bombyx mori

G protein-coupled receptors (GPCRs) are the largest and most versatile family of transmembrane receptors in the cell, occupying the highest hierarchical positions in the regulation of many physiological processes. Although they have been extensively studied in a number of model insects, there have been few investigations of GPCRs in large Lepidopterans, such as Bombyx mori, an organism that provides a means to perform detailed tissue expression analyses, which may help to characterize GPCRs and their ligands. In addition, B. mori, also known as the silkworm, is an insect of substantial economic importance, due to its use in silk production and traditional medicines. In this work, we computationally identified 90 putative GPCRs in B. mori, 33 of which represent novel proteins. These GPCRs were annotated and compared with their homologs in Drosophila melanogaster and Anopheles gambiae. Phylogenetics analyses of the GPCRs from these three insects showed that GPCRs may easily duplicate or disappear during insect evolution, especially in the neuropeptide and protein hormone receptor subfamily. Interestingly, we observed a decrease in the quantity and diversity of the stress-tolerance gene, Methuselah, in B. mori, which may be related to its long history of domestication. Moreover, the presence of many Bombyx-specific GPCRs suggests that neither Drosophila nor Anopheles is good representatives for the GPCRs in the Class Insecta.     

A functional screen for genes involved in Xenopus pronephros development

In Xenopus, the pronephros is the functional larval kidney and consists of two identifiable components; the glomus, the pronephric tubules, which can be divided into four separate segments, based on marker gene expression. The simplicity of this organ, coupled with the fact that it displays the same basic organization and function as more complex mesonephros and metanephros, makes this an attractive model to study vertebrate kidney formation. In this study, we have performed a functional screen specifically to identify genes involved in pronephros development in Xenopus. Gain-of-function screens are performed by injecting mRNA pools made from a non-redundant X. tropicalis full-length plasmid cDNA library into X. laevis eggs, followed by sib-selection to identify the single clone that caused abnormal phenotypes in the pronephros. Out of 768 egg and gastrula stage cDNA clones, 31 genes, approximately 4% of the screened clones, affected pronephric marker expression examined by whole mount in situ hybridization or antibody staining. Most of the positive clones had clear expression patterns in pronephros and predicted/established functions highly likely to be involved in developmental processes. In order to carry out a more detailed study, we selected Sox7, Cpeb3, P53csv, Mecr and Dnajc15, which had highly specific expression patterns in the pronephric region. The over-expression of these five selected clones indicated that they caused pronephric abnormalities with different temporal and spatial effects. These results suggest that our strategy to identify novel genes involved in pronephros development was highly successful, and that this strategy is effective for the identification of novel genes involved in late developmental events.     

Molecular characterization and expression patterns of the actinin-associated LIM protein (ALP) subfamily genes in porcine skeletal muscle

The actinin-associated LIM protein (ALP) subfamily has important functions in cell signal transduction, cell proliferation, and integration of cytoskeletal architecture. To detect their functions in pig skeletal muscle, we cloned and characterized the pig ALP subfamily genes, drew their genomic structure maps, and detected their tissue expression patterns. We identified a new spliced variant of PDLIM3 in pig skeletal muscle and named it as PDLIM3-4, which was only expressed in the heart and skeletal muscle. Our results showed that PDLIM3-4 was expressed in adult pig skeletal muscle with the highest expression level, and both PDLIM3-4 isoform and PDLIM4 had different expression profiles during the prenatal and postnatal stages of skeletal muscle development among the three pig breeds. These studies provide useful information for further research on the functions of pig ALP subfamily genes in skeletal muscle development.     

Mitochondrial protein BmPAPI modulates the length of mature piRNAs

PIWI proteins and their associated PIWI-interacting RNAs (piRNAs) protect genome integrity by silencing transposons in animal germlines. The molecular mechanisms and components responsible for piRNA biogenesis remain elusive. PIWI proteins contain conserved symmetrical dimethylarginines (sDMAs) that are specifically targeted by TUDOR domain-containing proteins. Here we report that the sDMAs of PIWI proteins play crucial roles in PIWI localization and piRNA biogenesis in Bombyx mori-derived BmN4 cells, which harbor fully functional piRNA biogenesis machinery. Moreover, RNAi screenings for Bombyx genes encoding TUDOR domain-containing proteins identified BmPAPI, a Bombyx homolog of Drosophila PAPI, as a factor modulating the length of mature piRNAs. BmPAPI specifically recognized sDMAs and interacted with PIWI proteins at the surface of the mitochondrial outer membrane. BmPAPI depletion resulted in 3′-terminal extensions of mature piRNAs without affecting the piRNA quantity. These results reveal the BmPAPI-involved piRNA precursor processing mechanism on mitochondrial outer membrane scaffolds.     

Establishment and characterization of the Bombyx mandarina cell line

A new cell line, designated as NIAS-Boma-529b, was established from the larval fat bodies of Bombyx mandarina (B. mandarina), which is believed to be an ancestor of Bombyx mori (B. mori). This cell line has been cultured for approximately 150 passages during 2 years in an IPL-41 medium supplemented with 10% fetal bovine serum at a constant temperature of 26 °C. The morphology of this line includes adhesive round and spindle-shaped cells. Random-amplified polymorphic DNA analysis (RAPD) using 7 primers and a statistical analysis based on Nei’s genetic distance revealed that this cell line was closely related to B. mori-derived cell lines. An infection study also revealed that this cell line was susceptible to B. mori nucleopolyhedrovirus (BmNPV); however, it had no apparent susceptibility to Autographa californica NPV (AcNPV), which is closely related to BmNPV. Nevertheless, cells infected with AcNPV showed an extensive cytopathic effect (CPE), including a rough cell surface, rounding, nuclear expansion, and cell blebbing. These results suggest that this cell line can be useful to clarify the mechanism of host range determination of BmNPV and AcNPV.     

Non-uniqueness of factors constraint on the codon usage in Bombyx mori

Background

The analysis of codon usage is a good way to understand the genetic and evolutionary characteristics of an organism. However, there are only a few reports related with the codon usage of the domesticated silkworm, Bombyx mori (B. mori). Hence, the codon usage of B. mori was analyzed here to reveal the constraint factors and it could be helpful to improve the bioreactor based on B. mori.

Results

A total of 1,097 annotated mRNA sequences from B. mori were analyzed, revealing there is only a weak codon bias. It also shows that the gene expression level is related to the GC content, and the amino acids with higher general average hydropathicity (GRAVY) and aromaticity (Aromo). And the genes on the primary axis are strongly positively correlated with the GC content, and GC3s. Meanwhile, the effective number of codons (ENc) is strongly correlated with codon adaptation index (CAI), gene length, and Aromo values. However, the ENc values are correlated with the second axis, which indicates that the codon usage in B. mori is affected by not only mutation pressure and natural selection, but also nucleotide composition and the gene expression level. It is also associated with Aromo values, and gene length. Additionally, B. mori has a greater relative discrepancy in codon preferences with Drosophila melanogaster (D. melanogaster) or Saccharomyces cerevisiae (S. cerevisiae) than with Arabidopsis thaliana (A. thaliana), Escherichia coli (E. coli), or Caenorhabditis elegans (C. elegans).

Conclusions

The codon usage bias in B. mori is relatively weak, and many influence factors are found here, such as nucleotide composition, mutation pressure, natural selection, and expression level. Additionally, it is also associated with Aromo values, and gene length. Among them, natural selection might play a major role. Moreover, the “optimal codons” of B. mori are all encoded by G and C, which provides useful information for enhancing the gene expression in B. mori through codon optimization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1596-z) contains supplementary material, which is available to authorized users.     

Effects of egt gene transfer on the development of Bombyx mori

This study investigated the effects of gain of ecdysteroid UDP-glucosyltransferase (EGT) gene function mutation on the development of the silkworm, Bombyx mori. A novel piggyBac-derived plasmid containing the egt gene from B. mori nucleopolyhedrovirus (BmNPV) driven by a heat-shock protein (hsp) 23.7 promoter, with a neomycin-resistance gene (neo) controlled by the BmNPV ie-1 promoter and a green fluorescent protein gene (gfp) under the control of the B. mori actin 3 (A3) promoter was constructed. The vector was transferred into silkworm eggs by sperm-mediated gene transfer. Transgenic silkworms were produced after screening for neo and gfp genes and gene transfer was verified by polymerase chain reaction, dot-blot hybridization and western blotting. The hatching rate of G1 generation silkworm eggs was about 60% lower than that of normal silkworm eggs. The duration of the G1 generation larval period was extended, and the G2 generation pupal stage lasted four days longer than that in non-transgenic silkworms. The ecdysone blood level in G2 silkworms in the third instar molting stage was reduced by up to 90%. These results show that EGT suppressed transgenic silkworm molting, and that egt expression in egt-transgenic silkworms resulted in arrest of metamorphosis from pupae to moths.