Studies on the transmission and tissue distribution of Antheraea pernyi iflavirus in the Chinese oak silkmoth Antheraea pernyi

Expanding known viral diversity in plants: virome of 161 species alongside an ancient canal

BACKGROUND

Since viral metagenomic approach was applied to discover plant viruses for the first time in 2006, many plant viruses had been identified from cultivated and non-cultivated plants. These previous researches exposed that the viral communities (virome) of plants have still largely uncharacterized. Here, we investigated the virome in 161 species belonging to 38 plant orders found in a riverside ecosystem.

RESULTS

We identified 245 distinct plant-associated virus genomes (88 DNA and 157 RNA viruses) belonging to 27 known viral families, orders, or unclassified virus groups. Some viral genomes were sufficiently divergent to comprise new species, genera, families, or even orders. Some groups of viruses were detected that currently are only known to infect organisms other than plants. It indicates a wider host range for members of these clades than previously recognized theoretically. We cannot rule out that some viruses could be from plant contaminating organisms, although some methods were taken to get rid of them as much as possible. The same viral species could be found in different plants and co-infections were common.

CONCLUSIONS

Our data describe a complex viral community within a single plant ecosystem and expand our understanding of plant-associated viral diversity and their possible host ranges.

Larval crowding effects during early development in the Chinese oak silkmoth Antheraea pernyi (Lepidoptera: Saturniidae)

Chinese sericulture relies in part on the rearing of the Chinese oak silkmoth Antheraea pernyi, an insect with key cultural and ecological roles. While feeding primarily on oak, Antheraea species are known to accept alternative hosts such as birch Betula sp with little to no apparent negative fitness consequences. This opens up the range of hostplants that could be used for large‐scale rearing of A. pernyi for silk production and food, or used by this species in possible invasions. To date, however, the natural history and ecology of A. pernyi remain subject of investigation. For instance, we still do not know how individuals respond to crowding developmental environments, which is an important factor to consider for the ecology of the species as well as for commercial rearing. Here, I describe the implications of larval crowding to the survival and growth of A. pernyi larvae during early development. I show that higher crowding is associated with stronger negative effects on growth and survival, corroborating findings from other holometabolous insects. I then discuss the implications of this findings for our understanding of optimum larval crowding. Overall, the findings reveal important ecological information for an insect species key for provisioning and cultural ecosystem services.

Infection of cotton bollworm by Helicoverpa armigera iflavirus decreases larval fitness

Previously, we reported a novel iflavirus in Helicoverpa armigera (helicoverpa armigera iflavirus, HaIV) and here we report the effects of HaIV on its host. In a laboratory bioassay, HaIV-positive larvae and pupae developed more slowly and had higher mortality than HaIV-negative larvae, suggesting that the virus is pathogenic. The relative fitness of H. armigera decreased with HaIV infection by a ratio of 0.65. Transcriptional analysis indicated that infection significantly changed the expression levels of host genes, with more genes affected at 72 h after inoculation than at 48 h (138 up- and 229 downregulated at 48 h; 185 up- and 299 downregulated at 72 h). Interestingly, pathways related to digestion and absorption were significantly enriched, e.g., protein digestion and absorption, suggesting developmental regulation of the host by HaIV via these pathways. HaIV-infected H. armigera showed significantly downregulated expression of genes encoding cuticular proteins (CPs), essential for structural and protective functions, at 48 h and 72 h, suggesting that HaIV increased larval mortality by downregulating CP gene expression.

Identification and RNAi Profile of a Novel Iflavirus Infecting Senegalese Aedes vexans arabiensis Mosquitoes

The inland floodwater mosquito Aedes vexans (Meigen, 1830) is a competent vector of numerous arthropod-borne viruses such as Rift Valley fever virus (Phenuiviridae) and Zika virus (Flaviviridae). Aedes vexans spp. have widespread Afrotropical distribution and are common European cosmopolitan mosquitoes. We examined the virome of Ae. vexans arabiensis samples from Barkédji village, Senegal, with small RNA sequencing, bioinformatic analysis, and RT-PCR screening. We identified a novel 9494 nt iflavirus (Picornaviridae) designated here as Aedes vexans iflavirus (AvIFV). Annotation of the AvIFV genome reveals a 2782 amino acid polyprotein with iflavirus protein domain architecture and typical iflavirus 5' internal ribosomal entry site and 3' poly-A tail. Aedes vexans iflavirus is most closely related to a partial virus sequence from Venturia canescens (a parasitoid wasp) with 56.77% pairwise amino acid identity. Analysis of AvIFV-derived small RNAs suggests that AvIFV is targeted by the exogenous RNA interference pathway but not the PIWI-interacting RNA response, as ~60% of AvIFV reads corresponded to 21 nt Dicer-2 virus-derived small RNAs and the 24-29 nt AvIFV read population did not exhibit a "ping-pong" signature. The RT-PCR screens of archival and current (circa 2011-2020) Ae. vexans arabiensis laboratory samples and wild-caught mosquitoes from Barkédji suggest that AvIFV is ubiquitous in these mosquitoes. Further, we screened wild-caught European Ae. vexans samples from Germany, the United Kingdom, Italy, and Sweden, all of which tested negative for AvIFV RNA. This report provides insight into the diversity of commensal Aedes viruses and the host RNAi response towards iflaviruses.

Biological characterization of Euscelidius variegatus iflavirus 1

Virus-based biocontrol technologies represent sustainable alternatives to pesticides and insecticides. Phytoplasmas are prokaryotic plant pathogens causing severe losses to crops worldwide. Novel approaches are needed since insecticides against their insect vectors and rogueing of infected plants are the only available strategies to counteract phytoplasma diseases. A new iflavirus, named EVV-1, has been described in the leafhopper phytoplasma vector Euscelidius variegatus, raising the potential to use virus-based application strategies against phytoplasma disease. Here transmission routes of EVV-1 are characterized, and localization within the host reveals the mechanism of insect tolerance to virus infection. Both vertical and horizontal transmission of EVV-1 occur and vertical transmission was more efficient. The virus is systemic and occurs in all life-stages, with the highest loads measured in ovaries and first to third instar nymphs. The basic knowledge gained here on the biology of the virus is crucial for possible future application of iflaviruses as biocontrol agents.

Transcriptome analysis reveals gene expression changes of the fat body of silkworm (Bombyx mori L.) in response to selenium treatment

A comparative transcriptome analysis was conducted to investigate the gene expression changes in the fat body of silkworm after treatment with different concentrations (50 μM and 200 μM) of selenium (Se). 912 differential expression genes (DEGs) (371 up-regulated and 541 down-regulated) and 1420 DEGs (1078 up-regulated and 342 down-regulated) were identified in silkworm fat body treated with 50 μM and 200 μM of Se, respectively. In case of 50 μM group, DEGs were mainly enriched in the peroxisome pathway and fatty acid metabolism pathway, and later were associated with antioxidant defense and nutrition regulation. After 200 μM Se-treatment, DEGs were mainly located in the glycerolipid metabolism and arachidonic acid metabolism pathways, which further encoded detoxification related genes. Furthermore, 32 candidate DEGs from these pathways had been selected to confirm the RNA-seq data. Among these DEGs, 14 genes were up-regulated in the 50 μM Se-treated group (only three genes in the 200 μM Se-treated group) which were involved in lipid metabolism and antioxidant defense, and 13 up-regulated genes (only two genes were up-regulated in the 50 μM Se-treated group) were involved in detoxification of the 200 μM Se-treated group. These changes showed that lower concentration of Se could regulate the nutrition and promote antioxidation pathways; whereas, high levels of Se promoted the detoxification of silkworm. These findings can be helpful to understand the possible mechanisms of Se action and detoxification in silkworm and other insects.

Codon Usage in the Iflaviridae Family Is Not Diverse Though the Family Members Are Isolated from Diverse Host Taxa

All iflavirus members belong to the unique genus, Iflavirus, of the family, Iflaviridae. The host taxa and sequence identities of these viruses are diverse. A codon usage bias, maintained by a balance between selection, mutation, and genetic drift, exists in a wide variety of organisms. We characterized the codon usage patterns of 44 iflavirus genomes that were isolated from the classes, Insecta, Arachnida, Mammalia, and Malacostraca. Iflaviruses lack a strong codon usage bias when they are evaluated using an effective number of codons. The odds ratios of the majority of dinucleotides are within the normal range. However, the dinucleotides at the 1st-2nd codon positions are more biased than those at the 2nd-3rd codon positions. Plots of effective numbers of codons, relative neutrality analysis, and PR2 bias analysis all indicate that selection pressure dominates mutations in shaping codon usage patterns in the family, Iflaviridae. When these viruses were grouped into their host taxa, we found that the indices, including the nucleotide composition, effective number of codons, relative synonymous codon usage, and the influencing factors behind the codon usage patterns, all show that there are non-significant differences between the six host-taxa-groups. Our results disagree with our assumption that diverse viruses should possess diverse codon usage patterns, suggesting that the nucleotide composition and codon usage in the family, Iflaviridae, are not host taxa-specific signatures.

Transcriptome Analysis Reveals Potential Antioxidant Defense Mechanisms in Antheraea pernyi in Response to Zinc Stress

The growth and development of the Chinese oak silkworm, Antheraea pernyi, are strongly influenced by environmental conditions, including heavy metal pollution. An excess of heavy metals causes cellular damage through the production of free radical reactive oxygen species. In this study, transcriptome analysis was performed to investigate global gene expression when A. pernyi was exposed to zinc infection. With RNA sequencing (RNA-Seq), a total of 25 795 510 and 38 158 855 clean reads were obtained from zinc-treated and control fat body libraries, respectively. We identified 2399 differential expression genes (DEGs) (1845 upregulated and 544 downregulated genes) in the zinc-treated library. In addition, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were related to the peroxisome pathway that was associated with antioxidant defense. Our results suggest that fat bodies of A. pernyi constitute a strong antioxidant defense against heavy metal contamination.

Characterization of a novel RNA virus from Nesidiocoris tenuis related to members of the genus Iflavirus

The complete genome of a novel virus from Nesidiocoris tenuis was determined by RNA-seq and rapid amplification of cDNA ends. This virus has a single-stranded RNA genome of 10633 nucleotides (nt) in length, not including the poly(A) tail, and contains two putative open reading frames (ORFs). ORF1 encodes a polypeptide of 1320 amino acids (aa) with a predicted molecular mass of 147.92 kDa and theoretical isoelectric point (pI) of 6.96. ORF2 encodes a polypeptide of 1728 aa with a predicted molecular mass of 197.09 kDa and pI of 6.73. Phylogenetic analysis with the deduced aa sequences of the conserved RNA dependent RNA polymerase domain as well as whole genome nt sequences indicated that the virus clusters with viruses classified within the genus Iflavirus, with a high bootstrap value in the maximum-likelihood and neighbor-joining trees. However, this virus has a distinct genome structure with two ORFs, iflaviruses normally having one, suggesting the virus might be a prototype of a new genus. We named the virus isolate Nesidiocoris tenuis virus 1 (NtV-1). The prevalence of NtV-1 infection in wild samples of N. tenuis was at a low level (7.32%, 6 positive in 82 samples), suggesting a possible harmful effect to its host.