Survey of long terminal repeat retrotransposons of domesticated silkworm (Bombyx mori)

Long terminal retrotransposons are major components of eukaryotic transposable elements. We have surveyed the long terminal repeats (LTR) retrotransposons of domesticated silkworm (Bombyx mori) by mining the data produced by Bombyx mori Genome Sequencing Project. At least 29 separate families of LTR retrotransposons are identified in this survey, comprising of 11.8% of the complete sequence. Families of domesticated silkworm LTR retrotransposons can be mainly classified into three groups: gypsy-like, copia-like, Pao-Bel. Fourteen families identified consist of gypsy-like elements, four families consist of copia-like elements and seven families consist of Pao-Bel elements. In addition to the three groups of LTR retrotransposons, two families of unusual non-coding elements are identified in the genome of this species. Further phylogenetic analysis of RT domain indicates that the elements of B.mori show high diversity and can form different clades in each group. An analysis of sequence variation from different families reveals distinct patterns of variation for the elements belonging to three groups. The analysis of the domesticated silkworm LTR retrotransposons should assist in our understanding of the roles of retroelement in lepidopteron insect genome evolution.     

Genome-wide analysis of cytochrome P450 monooxygenase genes in the silkworm, Bombyx mori

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed. A total of 84 CYP-related sequences were identified and could be classified into 26 families and 47 subfamilies according to standard nomenclature. Seventy eight of the eighty four genes appear to be functional and six are probable pseudogenes. The distribution of Bombyx mori P450s in the genome shows that most of them are tandem arranged on chromosomes, only 34 genes are present as singletons, with 8 clusters including 3 or more than 3 genes. Sequence alignments were used to reconstruct phylogenetic trees and to analyze the intron-exon organizations of the functional genes. The conserved intron positioning agrees perfectly with their common grouping on the tree. The presence of three extremely ancient introns which are conserved across different clans indicates that a few introns are still highly conserved after they have undergone extensive evolutionary changes of B. mori P450 duplication and divergence. Comparison of the P450s from B. mori to the P450s from Drosophila melanogaster shows that the expansion is not uniform across the gene families. Remarkably, two mitochondrial families, the B. mori CYP333 and D. melanogaster Cyp12, formed two orthologous groups in the phylogenetic tree. All CYP333s can be proposed to be related to xenobiotic metabolism in accordance with the D. melanogaster Cyp12s. The characterization and evolutionary analysis of P450s from B. mori in the current study provide useful information for understanding the characteristics and diversity of P450s from B. mori and the baseline for functional analyses of individual P450s in this model Lepidopteran insect.     

The search of miRNA genes in Bombyx mori nuclear polyhedrosis virus genomes regions complementary to the latest genes

The search of miRNA genes in Bombyx mori nuclear polyhedrosis virus genome region complementary to very late genes has been carried out. The search miRNA algorithm in silico was developed by us. It was shown that NPV B. mori genome region containing orf4 gene complementary to ph gene encodes the potential miRNA. NPV B. mori genome region containing p74 gene complementary to p10 gene encodes mature miRNA and potential miRNA. The genome region containing orf1629 encodes two small non-coding RNAs complementary to orf 5'-end of polyhedrin miRNA. From obtained results it is proposed that two small noncoding RNAs complementary to regions of polyhedrin miRNA are included in polyhedra.     

Comparative genetic diversity and genetic structure of three Chinese silkworm species Bombyx mori L. (Lepidoptera: Bombycidae), Antheraea pernyi Guérin-Meneville and Samia cynthia ricini Donovan (Lepidoptera: Saturniidae)

The genetic diversity and genetic structure of three Chinese silkworm species Bombyx mori L., Antheraea pernyi Guérin-Meneville and Samia cynthia ricini Donovan were comparatively assessed based on RAPD markers. At the species level, A. pernyi and B. mori showed high levels of genetic diversity, whereas S. cynthia ricini showed low level of genetic diversity. However, at the strain level, A. pernyi had relatively highest genetic diversity and B. mori had lowest genetic diversity. Analysis of molecular variance (AMOVA) suggested that 60% and 72% of genetic variation resided within strains in A. pernyi and S. cynthia ricini, respectively, whereas only 16% of genetic variation occurred within strains in B. mori. In UPGMA dendrogram, individuals of A. pernyi and B. mori formed the strain-specific genetic clades, whereas those of S. cynthia ricini were distributed in a mixed way. The implications of these results for the conservation and utilization in breeding programs of three silkworm species are discussed.     

家蚕微孢子虫中小RNAs的全基因组鉴定与分析

【目的】家蚕Bombyx mori微粒子病是蚕业生产上的毁灭性病害,家蚕微孢子虫Nosema bombycis是该病的病原,可经卵垂直传播和经口水平传播。为了探索家蚕微孢子虫中对重复元件的抵御以及对基因转录调控的潜在方式,本研究拟在基因组水平上对该物种的小RNAs进行全面系统的分析,鉴定与转座子相关的小RNAs和潜在的miRNAs。【方法】从感染家蚕微孢子虫的家蚕中肠中提取总RNA,分离小片段RNA并反转录后,进行Solexa高通量测序。通过生物信息学方法对小RNAs进行分类及功能注释,鉴定起源于家蚕微孢子虫不同类型转座子的小RNAs,并对潜在的miRNA进行预测分析。【结果】家蚕微孢子虫小RNAs的长度主要是24和25 nt,其中大部分序列表现出5′末端的尿嘧啶偏好性。家蚕微孢子虫中存在丰富的与转座子相关联的小RNAs,并且与转座子标准序列匹配的反义小RNAs明显多于正义小RNAs。同时,鉴定获得了31个候选miRNAs,部分为Nosema属的其他孢子虫中所共有,暗示其在微孢子虫基因组进化上具有保守性。【结论】首次鉴定到家蚕微孢子虫的转座子相关性小RNAs,暗示小RNAs在家蚕微孢子虫基因组对转座子防御过程中起到作用,31个潜在的miRNAs为家蚕微孢子虫miRNAs的功能验证提供了后续靶标。     

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.     

Synteny and chromosome evolution in the lepidoptera: evidence from mapping in Heliconius melpomene

The extent of conservation of synteny and gene order in the Lepidoptera has been investigated previously only by comparing a small subset of linkage groups between the moth Bombyx mori and the butterfly Heliconius melpomene. Here we report the mapping of 64 additional conserved genes in H. melpomene, which contributed 47 markers to a comparative framework of 72 orthologous loci spanning all 21 H. melpomene chromosomes and 27 of the 28 B. mori chromosomes. Comparison of the maps revealed conserved synteny across all chromosomes for the 72 loci, as well as evidence for six cases of chromosome fusion in the Heliconius lineage that contributed to the derived 21-chromosome karyotype. Comparisons of gene order on these fused chromosomes revealed two instances of colinearity between H. melpomene and B. mori, but also one instance of likely chromosomal rearrangement. B. mori is the first lepidopteran species to have its genome sequenced, and the finding that there is conserved synteny and gene order among Lepidoptera indicates that the genomic tools developed in B. mori will be broadly useful in other species.     

Mitochondrial genome nucleotide substitution pattern between domesticated silkmoth, Bombyx mori, and its wild ancestors, Chinese Bombyx mandarina and Japanese Bombyx mandarina

Bombyx mori and Bombyx mandarina are morphologically and physiologically similar. In this study, we compared the nucleotide variations in the complete mitochondrial (mt) genomes between the domesticated silkmoth, B. mori, and its wild ancestors, Chinese B. mandarina (ChBm) and Japanese B. mandarina (JaBm). The sequence divergence and transition mutation ratio between B. mori and ChBm are significantly smaller than those observed between B. mori and JaBm. The preference of transition by DNA strands between B. mori and ChBm is consistent with that between B. mori and JaBm, however, the regional variation in nucleotide substitution rate shows a different feature. These results suggest that the ChBm mt genome is not undergoing the same evolutionary process as JaBm, providing evidence for selection on mtDNA. Moreover, investigation of the nucleotide sequence divergence in the A+T-rich region of Bombyx mt genomes also provides evidence for the assumption that the A+T-rich region might not be the fastest evolving region of the mtDNA of insects.     

Molecular phylogeny of silkmoths reveals the origin of domesticated silkmoth, Bombyx mori from Chinese Bombyx mandarina and paternal inheritance of Antheraea proylei mitochondrial DNA

Molecular phylogeny of some of the economically important silkmoths was derived using three mitochondrial genes, 12S rRNA, 16S rRNA, and COI, and the control region (CR). Maximum likelihood (ML) analyses showed two distinct clades, one consisting of moths from Bombycidae family and the other from Saturniidae family. The mitochondrial CR showed length polymorphisms with indels. The ML analyses for complete mitochondrial genome sequences of Bombyx mori (strains Aojuku, C108, Backokjam, and Xiafang), Japanese and Chinese strains of B. mandarina (Japanese mandarina and Chinese mandarina) and, Antheraea pernyi revealed two distinct clades, one comprising of B. mori strains and the other with B. mandarina, and A. pernyi forming an outgroup. Pairwise distances revealed that all of the strains of B. mori studied are closer to Chinese than to Japanese mandarina. Phylogenetic analyses based on whole mitochondrial genome sequences, the finding of a tandem triplication of a 126bp repeat element only in Japanese mandarina, and chromosome number variation in B. mandarina suggest that B. mori must have shared its recent common ancestor with Chinese mandarina. Another wild species of the Bombycidae family, Theophila religiosa, whose phylogenetic status was not clear, clustered together with the other bombycid moths in the study. Analysis of the interspecific hybrid, A. proylei gave evidence for paternal inheritance of mitochondrial DNA.     

Isolation and partial characterization of U1-U6 small RNAs from Bombyx mori

We have used a variety of techniques to characterize the U-series small nuclear RNAs from the posterior silk gland of Bombyx mori. Six molecular species have been identified which correspond to the vertebrate U1-U6 RNAs by the following criteria: (a) presence of the RNAs in ribonucleoprotein particles which can be immunoprecipitated by lupus Sm antisera; (b) presence of a 2,2,7-trimethylguanosine cap, as assayed by immunoprecipitation with anti-2,2,7-trimethylguanosine IgG; (c) size, as assayed by acrylamide/urea gel electrophoresis using HeLa cell U-RNA markers; and (d) primary nucleotide sequence, as determined by chemical/enzymatic cleavage of end-labeled molecules. The high conservation of primary sequence (66-81% homology based on partial sequences) relative to the corresponding vertebrate U-RNAs has permitted unambiguous identification of each molecule. With the exception of two subspecies of U3 RNA, the U-snRNAs of Bombyx exhibit a striking conservation of secondary structure relative to the proposed structures of the U-RNAs of vertebrates. This conservation is best exemplified by several compensatory base alterations that result in the maintenance of hairpin structures. These are particularly evident in U1 and U5 RNAs. Bombyx U3 is interesting in that two subspecies (of a total of four that were sequenced) diverge considerably in sequence (and presumably in structure) relative to the U3 RNA of vertebrates. The most abundant U-RNAs in the posterior silk gland appear to be U1 and U2, while U3-U6 are present in relatively small amounts.     

A covariant change of the two highly conserved bases in the GTPase-associated center of 28 S rRNA in silkworms and other moths

The GTPase-associated center in 23/28 S rRNA is one of the most conserved functional domains throughout all organisms. We detected a unique sequence of this domain in Bombyx mori species in which the bases at positions 1094 and 1098 (numbering from Escherichia coli 23 S rRNA) are C and G instead of the otherwise universally conserved bases U and A, respectively. These changes were also observed in four other species of moths, but not in organisms other than the moths. Characteristics of the B. mori rRNA domain were investigated by native polyacrylamide gel electrophoresis using RNA fragments containing residues 1030-1128. Although two bands of protein-free RNA appeared on gel, they shifted to a single band when bound to Bombyx ribosomal proteins Bm-L12 and Bm-P complex, equivalent to E. coli L11 and L8, respectively. Bombyx RNA showed lower binding capacity than rat RNA for the ribosomal proteins and anti-28 S autoantibody, specific for a folded structure of the eukaryotic GTPase-associated domain. When the C(1094)/G(1098) bases in Bombyx RNA were replaced by the conserved U/A bases, the protein-free RNA migrated as a single band, and the complex formation with Bm-L12, Bm-P complex, and anti-28 S autoantibody was comparable to that of rat RNA. The results suggest that the GTPase-associated domain of moth-type insects has a labile structural feature that is caused by an unusual covariant change of the U(1094)/A(1098) bases to C/G.     

A genome-wide analysis of genes and gene families involved in innate immunity of Bombyx mori

A genome-wide analysis of innate immunity-related genes and gene families was conducted using the silkworm, Bombyx mori. We identified orthologs for a large number of genes involved in insect immunity that have been reported from Drosophila melanogaster (Diptera), Anopheles gambiae (Diptera), Apis mellifera (Hymenoptera) and Tribolium castaneum (Coleoptera). B. mori has a unique recognition gene and antimicrobial peptide genes that are not present in the Drosophila, Anopheles, Apis and Tribolium genomes, suggesting a lineage-specific gene evolution for lepidopteran insects. The comparative analysis of the insect immune repertoires indicated a dynamic and flexible gene expansion in recognition, modulation and effector mechanisms due to different selection pressures. Differential gene regulation by different bacterial species was found in PGRP and Serpin genes, suggesting that Bombyx has a highly selective gene regulation system depending on bacterial species.     

Characterization of unique small RNA populations from rice grain

Small RNAs (approximately 20 to 24 nucleotides) function as naturally occurring molecules critical in developmental pathways in plants and animals. Here we analyze small RNA populations from mature rice grain and seedlings by pyrosequencing. Using a clustering algorithm to locate regions producing small RNAs, we classified hotspots of small RNA generation within the genome. Hotspots here are defined as 1 kb regions within which small RNAs are significantly overproduced relative to the rest of the genome. Hotspots were identified to facilitate characterization of different categories of small RNA regulatory elements. Included in the hotspots, we found known members of 23 miRNA families representing 92 genes, one trans acting siRNA (ta-siRNA) gene, novel siRNA-generating coding genes and phased siRNA generating genes. Interestingly, over 20% of the small RNA population in grain came from a single foldback structure, which generated eight phased 21-nt siRNAs. This is reminiscent of a newly arising miRNA derived from duplication of progenitor genes. Our results provide data identifying distinct populations of small RNAs, including phased small RNAs, in mature grain to facilitate characterization of small regulatory RNA expression in monocot species.     

Coronavirus multiplication strategy. II. Mapping the avian infectious bronchitis virus intracellular RNA species to the genome.

Avian infectious bronchitis virus, a coronavirus, directed the synthesis of six major single-stranded polyadenylated RNA species in infected chicken embryo kidney cells. These RNAs include the intracellular form of the genome (RNA F) and five smaller RNA species (RNAs A, B, C, D, and E). Species A, B, C, and D are subgenomic RNAs and together with the genome form a nested sequence set, with the sequences of each RNA contained within every larger RNA species (D. F. Stern and S. I. T. Kennedy, J. Virol 34:665-674, 1980). In the present paper we show by RNase T1 oligonucleotide fingerprinting that RNA E is also a member of the nested set. Partial alkaline fragmentation of the genome followed by sucrose fractionation, oligodeoxythymidylate-cellulose chromatography, and RNase T1 fingerprinting gave a partial 3'-to-5' oligonucleotide spot order. A comparison of the oligonucleotides of each of the five subgenomic RNAs with this spot order established that all of the RNAs are comprised of nucleotide sequences inward from the 3' end of the genome. This result is discussed in relation to the multiplication strategy both of coronaviruses and of other RNA-containing viruses.     

Endogenous siRNAs and piRNAs derived from transposable elements and genes in the malaria vector mosquito Anopheles gambiae

Background

The siRNA and piRNA pathways have been shown in insects to be essential for regulation of gene expression and defence against exogenous and endogenous genetic elements (viruses and transposable elements). The vast majority of endogenous small RNAs produced by the siRNA and piRNA pathways originate from repetitive or transposable elements (TE). In D. melanogaster, TE-derived endogenous siRNAs and piRNAs are involved in genome surveillance and maintenance of genome integrity. In the medically relevant malaria mosquito Anopheles gambiae TEs constitute 12-16% of the genome size. Genetic variations induced by TE activities are known to shape the genome landscape and to alter the fitness in An. gambiae.

Results

Here, using bioinformatics approaches we analyzed the small RNA data sets from 6 libraries formally reported in a previous study and examined the expression of the mixed germline/somatic siRNAs and piRNAs produced in adult An. gambiae females. We characterized a large population of TE-derived endogenous siRNAs and piRNAs, which constitutes 56-60% of the total siRNA and piRNA reads in the analysed libraries. Moreover, we identified a number of protein coding genes producing gene-specific siRNAs and piRNAs that were generally expressed at much lower levels than the TE-associated small RNAs. Detailed sequence analysis revealed that An. gambiae piRNAs were produced by both “ping-pong” dependent (TE-associated piRNAs) and independent mechanisms (genic piRNAs). Similarly to D. melanogaster, more than 90% of the detected piRNAs were produced from TE-associated clusters in An. gambiae. We also found that biotic stress as blood feeding and infection with Plasmodium parasite, the etiological agent of malaria, modulated the expression levels of the endogenous siRNAs and piRNAs in An. gambiae.

Conclusions

We identified a large and diverse set of the endogenously derived siRNAs and piRNAs that share common and distinct aspects of small RNA expression across insect species, and inferred their impact on TE and gene activity in An. gambiae. The TE-specific small RNAs produced by both the siRNA and piRNA pathways represent an important aspect of genome stability and genetic variation, which might have a strong impact on the evolution of the genome and vector competence in the malaria mosquitoes.

Electronic supplementary material

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