Selective pressure dominates the synonymous codon usage in parvoviridae.

Parvoviridae is a family of small non-enveloped viruses and divided into two subfamilies. The family members infect a wide range of organisms from insects to humans and some of the members (e.g., nonpathogenic adeno-associated viruses) are effective gene therapy delivery vectors. We detailed the synonymous codon usage pattern of Parvoviridae family from the available 58 sequenced genomes through multivariate statistical methods. Our results revealed that nine viruses showed some degree of strong codon bias, and the others possessed a general weak trend of codon bias. ENc-plot and neutrality plot results showed that selective pressure dominated over mutation in shapes coding sequence's composition. The overall GC content and GC content at the third synonymous codon position were the principal determinants behind the variations within the codon usage patterns, as they both significantly correlated with the first axis of correspondence analysis. In addition, gene length had no direct influence on the codon usage pattern. Densovirinae subfamily and Parvovirinae subfamily possessed nine identical preferred codons, though most of the two subfamilies codon usage frequencies were significantly different. The result of cluster analysis based on synonymous codon usage was discordant with that of taxonomic classification. Adeno-associated viruses formed a separated clade far from other Parvoviridae members in the dendrogram. Thus, we concluded that natural selection rather than mutation pressure accounts for the main factor that affects the codon bias in Parvoviridae family.

Characterization of a highly conserved Antheraea pernyi spermidine synthase gene.

In the present study, we isolated a spermidine synthase gene from Antheraea pernyi (ApSpds) using expressed sequence tag method. The obtained cDNA sequence of 1483 bp contains an open-reading frame of 864 bp encoding a polypeptide of 287 amino acids. Sequence analysis revealed that ApSpds belonged to class I of AdoMet-MTase family, and exhibited 30% identity to those from bacteria, 45-48% identity to fungi, 36-47% identity to plants, 52-54% identity to vertebrates and 53-80% identity to invertebrates. Phylogenetic analysis found that the used Spds protein sequences were well divided into five groups corresponding to bacteria, fungi, plants, invertebrates and vertebrates, respectively. These results further confirmed that Spds is highly conserved through evolution of life organisms. The ApSpds mRNA is expressed during all four developmental stages and is present in all examined tissues with the highest abundance in the muscle, in which the relative mRNA expression level was 1.6 times higher than in the fat body. Although not significant, the mRNA level decreased after high-temperature exposure suggesting that the Spds gene may not be involved in temperature stress tolerance in A. pernyi. Taken together, our results suggested that ApSpds play an important role in development of silkworm.

Transcriptomic analysis of the prothoracic gland from two lepidopteran insects, domesticated silkmoth Bombyx mori and wild silkmoth Antheraea pernyi.

The prothoracic gland (PG) is an important endocrine organ of synthesis and secretion of ecdysteroids that play critical roles in insects. Here, we used a comparative transcriptomic approach to characterize some common features of PGs from two lepidopteran species Bombyx mori and Antheraea pernyi. Functional and pathway annotations revealed an overall similarity in gene profile between the two PG transcriptomes. As expected, almost all steroid hormone biosynthesis genes and the prothoracicitropic hormone receptor gene (Torso) were well represented in the two PGs. Impressively, two ecdysone receptor genes, eleven juvenile hormone related genes, more than 10 chemosensory protein genes, and a set of genes involved in circadian clock were also presented in the two PGs. Quantitative real time -PCR (qRT-PCR) validated the expression of 8 juvenile hormone and 12 clock related genes in B. mori PG, and revealed a different expression pattern during development in whole fifth larval instar. This contribution to insect PG transcriptome data will extend our understanding of the function and regulation of this important organ.

Upregulation of a Trypsin-Like Serine Protease Gene in Antheraea pernyi (Lepidoptera: Saturniidae) Strains Exposed to Different Pathogens.

Antheraea pernyi Guérin-Méneville is used for silk production and as a food resource. Its infection by exogenous pathogens, including microsporidia, fungi, bacteria, and virus, can lead to silkworm diseases, causing major economic losses. A trypsin-like serine protease gene (TLS) was found in A. pernyi transcriptome data resulting from two different infection experiments. The cDNA sequence of ApTLS was 1,020 bp in length and contained an open reading frame of 774 bp encoding a 257-amino acid protein (GenBank KF779933). The present study investigated the expression patterns of ApTLS after exposure to different pathogens, and in four different A. pernyi strains. Semiquantitative RT-PCR indicated that ApTLS was expressed in all developmental stages and was most expressed in the midgut. Quantitative real-time PCR indicated ApTLS was upregulated in the midgut of A. pernyi exposed to nucleopolyhedrovirus (ApNPV), Nosema pernyi, Enterococcus pernyi, and Beauveria bassiana infections, and the highest gene expression level was found under ApNPV infection. The strain Shenhuang No. 2 presented the lowest infection rate and the highest ApTLS gene expression level when exposed to ApNPV. Thus, ApTLS seems to be involved in innate defense reactions in A. pernyi.

Codon usage in Alphabaculovirus and Betabaculovirus hosted by the same insect species is weak, selection dominated and exhibits no more similar patterns than expected.

Mutations shape synonymous codon usage bias in certain organism genomes, while selection shapes it in others. Lepidopteran-specific Alphabaculovirus and Betabaculovirus are two large genera in the family of Baculoviridae. In this study, we analyzed the codon usage patterns in 17 baculoviruses, including 10 alphabaculoviruses and 7 betabaculoviruses, which were isolated from seven insect species, and we characterized the codon usage patterns between Alphabaculovirus and Betabaculovirus. Our results show that all the baculoviruses possessed a general weak trend of codon bias. The differences of ENc (effective number of codons) values, nucleotide contents and the impacts of nucleotide content on ENc value within alpha-/betabaculovirus pairs were independent of whether the host species are the same or different. Furthermore, the majority of amino acid sequences adopted codons unequally in all viruses, but the numbers of common preferred codons between alpha- and betabaculoviruses hosted by the same insect species were not significantly different from the differences observed between alpha- and betabaculoviruses hosted by different insect species. In addition, the amino acids that adopt the same synonymous codon composition between alpha- and betabaculoviruses hosted by the same insect species were statistically as few as those between alpha- and betabaculoviruses hosted by different insect species. Correspondence analysis revealed that no major factors resulted in the codon bias in these baculoviruses, implying multiple minor influential factors exist. Neutrality plot analysis indicated that selection pressure dominated mutations in shaping the codon usage. However, the levels of selection pressure were not significantly different among viruses hosted by the same insect species. We expect that evolution would cause the alpha- and betabaculoviruses hosted by the same insect species to share more patterns, but this effect was not observed.

Transcriptome Analysis of the Midgut of the Chinese Oak Silkworm Antheraea pernyi Infected with Antheraea pernyi Nucleopolyhedrovirus.

The Antheraea pernyi nucleopolyhedrovirus (ApNPV) is an exclusive pathogen of A. pernyi. The intense interactions between ApNPV and A. pernyi cause a series of physiological and pathological changes to A. pernyi. However, no detailed report exists regarding the molecular mechanisms underlying the interactions between ApNPV and A. pernyi. In this study, four cDNA libraries of the A. pernyi midgut, including two ApNPV-infected groups and two control groups, were constructed for transcriptomic analysis to provide new clues regarding the molecular mechanisms that underlie these interactions. The transcriptome of the A. pernyi midgut was de novo assembled using the Trinity platform because of the lack of a genome resource for A. pernyi. Compared with the controls, a total of 5,172 differentially expressed genes (DEGs) were identified, including 2,183 up-regulated and 2,989 down-regulated candidates, of which 2,965 and 911 DEGs were classified into different GO categories and KEGG pathways, respectively. The DEGs involved in A. pernyi innate immunity were classified into several categories, including heat-shock proteins, apoptosis-related proteins, serpins, serine proteases and cytochrome P450s. Our results suggested that these genes were related to the immune response of the A. pernyi midgut to ApNPV infection via their essential roles in regulating a variety of physiological processes. Our results may serve as a basis for future research not only on the molecular mechanisms of ApNPV invasion but also on the anti-ApNPV mechanism of A. pernyi.