Differential gene expression and microRNA profile in corpora allata-corpora cardiaca of Aedes aegypti mosquitoes with weak juvenile hormone signalling.

The corpora allata-corpora cardiaca (CA-CC) is an endocrine gland complex that regulates mosquito development and reproduction through the synthesis of juvenile hormone (JH). Epoxidase (Epox) is a key enzyme in the production of JH. We recently utilized CRISPR/Cas9 to establish an epoxidase-deficient (epox-/-) Aedes aegypti line. The CA from epox-/- mutants do not synthesize epoxidated JH III but methyl farneosate (MF), a weak agonist of the JH receptor, and therefore have reduced JH signalling. Illumina sequencing was used to examine the differences in gene expression between the CA-CC from wild type (WT) and epox-/- adult female mosquitoes. From 18,034 identified genes, 317 were significantly differentially expressed. These genes are involved in many biological processes, including the regulation of cell proliferation and apoptosis, energy metabolism, and nutritional uptake. In addition, the same CA-CC samples were also used to examine the microRNA (miRNA) profiles of epox-/- and WT mosquitoes. A total of 197 miRNAs were detected, 24 of which were differentially regulated in epox-/- mutants. miRNA binding sites for these particular miRNAs were identified using an in silico approach; they target a total of 101 differentially expressed genes. Our results suggest that a lack of epoxidase, besides affecting JH synthesis, results in the diminishing of JH signalling that have significant effects on Ae. aegypti CA-CC transcriptome profiles, as well as its miRNA repertoire.

Aedes aegypti gut transcriptomes respond differently to microbiome transplants from field-caught or laboratory-reared mosquitoes.

The mosquito microbiome is critical for host development and plays a major role in many aspects of mosquito biology. While the microbiome is commonly dominated by a small number of genera, there is considerable variation in composition among mosquito species, life stages, and geography. How the host controls and is affected by this variation is unclear. Using microbiome transplant experiments, we asked whether there were differences in transcriptional responses when mosquitoes of different species were used as microbiome donors. We used microbiomes from four different donor species spanning the phylogenetic breadth of the Culicidae, collected either from the laboratory or the field. We found that when recipients received a microbiome from a donor reared in the laboratory, the response was remarkably similar regardless of donor species. However, when the donor had been collected from the field, many more genes were differentially expressed. We also found that while the transplant procedure did have some effect on the host transcriptome, this is likely to have had a limited effect on mosquito fitness. Overall, our results highlight the possibility that variation in mosquito microbiome communities is associated with variability in host-microbiome interactions and further demonstrate the utility of the microbiome transplantation technique for investigating host-microbe interactions in mosquitoes.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros).

BACKGROUND: An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation. RESULTS: High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles. CONCLUSIONS: The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

Transcriptional response of Wolbachia-transinfected Aedes aegypti mosquito cells to dengue virus at early stages of infection.

Mosquito-borne flaviviruses are responsible for viral infections and represent a considerable public health burden. Aedes aegypti is the principal vector of dengue virus (DENV), therefore understanding the intrinsic virus-host interactions is vital, particularly in the presence of the endosymbiont Wolbachia, which blocks virus replication in mosquitoes. Here, we examined the transcriptional response of Wolbachia-transinfected Ae. aegypti Aag2 cells to DENV infection. We identified differentially expressed immune genes that play a key role in the activation of anti-viral defence such as the Toll and immune deficiency pathways. Further, genes encoding cytosine and N6-adenosine methyltransferases and SUMOylation, involved in post-transcriptional modifications, an antioxidant enzyme, and heat-shock response were up-regulated at the early stages of DENV infection and are reported here for the first time. Additionally, several long non-coding RNAs were among the differentially regulated genes. Our results provide insight into Wolbachia-transinfected Ae. aegypti's initial virus recognition and transcriptional response to DENV infection.

Transcription Profile and Genomic Variations of Oryctes Rhinoceros Nudivirus in Coconut Rhinoceros Beetles.

Oryctes rhinoceros nudivirus (OrNV) is a double-stranded DNA (dsDNA) virus which has been used as a biocontrol agent to suppress the coconut rhinoceros beetle (Oryctes rhinoceros) in Southeast Asia and the Pacific Islands. A new wave of O. rhinoceros incursions in Oceania is thought to be related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. In this study, chronically infected beetles were collected from Philippines, Fiji, Papua New Guinea (PNG), and the Solomon Islands (SI). RNA sequencing (RNA-seq) was performed to investigate the global viral gene expression profiles and for comparative genomic analysis of structural variations. Maximum likelihood phylogenic analysis indicated that OrNV strains from the SI and Philippines are closely related, while OrNV strains from PNG and Fiji formed a distinct adjacent clade. We detected several polymorphic sites with a frequency higher than 35% in 892 positions of the viral genome. Nonsynonymous mutations were detected in several hypothetical proteins and 15 nudivirus core genes, such as gp034, lef-8, lef-4, and vp91 We found limited evidence of variation in viral gene expression among geographic populations. Only a few genes, such as gp01, gp022, and gp107, were differentially expressed among different strains. Additionally, small RNA sequencing from the SI population suggested that OrNV is targeted by the host RNA interference (RNAi) response with abundant 21-nucleotide small RNAs. Some of these genomic changes are specific to the geographic population and could be related to particular phenotypic characteristics of the strain, such as viral pathogenicity or transmissibility, and this requires further investigation.IMPORTANCE Oryctes rhinoceros nudivirus has been an effective biocontrol agent against the coconut rhinoceros beetle in Southeast Asia and the Pacific Islands for decades. The recent outbreak of these beetles in many South Pacific islands has had a significant impact on livelihoods in the region. It has been suggested that the resurgence and spread of the pest are related to the presence of low-virulence isolates of OrNV or virus-tolerant haplotypes of beetles. We examined viral genomic and transcriptional variations in chronically infected beetles from different geographical populations. A high number of polymorphic sites among several geographical strains of OrNV were identified, but potentially only a few of these variations in the genome are involved in functional changes and can potentially alter the typical function. These findings provide valuable resources for future studies to improve our understanding of the OrNV genetic variations in different geographic regions and their potential link to virus pathogenicity.

A new dicistro-like virus from soldier fly, Inopus flavus (Diptera: Stratiomyidae), a pest of sugarcane.

Native Australian soldier flies, Inopus spp. (Diptera: Stratiomyidae), are agricultural pests of economic importance to the sugarcane industry. A screen of the salivary gland transcriptome of Inopus flavus (James) revealed the presence of viral RNA belonging to a potentially novel member of the family Dicistroviridae. The complete genome sequence consists of 9793 nucleotides with two open reading frames. The genome includes two potential internal ribosomal entry sites (IRESs): one within the 5' UTR and the other in the intergenic region (IGR). Virus particles purified from infected larvae and visualised by electron microscopy were found to be icosahedral, non-enveloped, and 30 nm in diameter.

Understanding the role of microRNAs in the interaction of Aedes aegypti mosquitoes with an insect-specific flavivirus.

The Flavivirus genus contains some of the most prevalent vector-borne viruses, such as the dengue, Zika and yellow fever viruses that cause devastating diseases in humans. However, the insect-specific clade of flaviviruses is restricted to mosquito hosts, albeit they have retained the general features of the genus, such as genome structure and replication. The interactions between insect-specific flaviviruses (ISFs) and their mosquito hosts are largely unknown. Pathogenic flaviviruses are known to modulate host-derived microRNAs (miRNAs), a class of non-coding RNAs that are important in controlling gene expression. Alterations in miRNAs may represent changes in host gene expression and promote understanding of virus-host interactions. The role of miRNAs in ISF-mosquito interactions is largely unknown. A recently discovered Australian ISF, Palm Creek virus (PCV), has the ability to suppress medically relevant flaviviruses. Here, we investigated the potential involvement of miRNAs in PCV infection using the model mosquito Aedes aegypti. By combining small-RNA sequencing and bioinformatics analysis, differentially expressed miRNAs were determined. Our results indicated that PCV infection hardly affects host miRNAs. Out of 101 reported miRNAs of Ae. aegypti, only aae-miR-2940-5p had a significantly altered expression over the course of infection. However, further analysis of aae-miR-2940-5p revealed that this miRNA does not have any direct impact on PCV replication in vitro. Thus, overall the results suggest that PCV infection has a limited effect on the mosquito miRNA profile and therefore miRNAs may not play a significant role in the PCV-Ae. aegypti interaction.

Discovery of new orbiviruses and totivirus from Anopheles mosquitoes in Eastern Australia.

Three new viruses classifiable within the Totivirus and Orbivirus genera were detected from Anopheles mosquito species collected in Eastern Australia. The viruses could not be isolated in C6/36 mosquito cell cultures but were shown to replicate in their mosquito hosts by small RNA analysis. The viruses grouped phylogenetically with other viruses recently detected in insects. These discoveries contribute to a better understanding of commensal viruses in Australian mosquitoes and the evolution of these viruses.

Wolbachia suppresses cell fusing agent virus in mosquito cells.

The genus Flavivirus contains a large number of positive-sense ssRNA viruses. While some are transmitted by mosquitoes or other arthropods and are pathogenic to humans and animals (e.g. dengue and Zika viruses), some are insect-specific and do not replicate in vertebrate cells. These are known as insect-specific flaviviruses (ISFs). Cell fusing agent virus (CFAV) was the first described ISF, which was detected in an Aedes aegypti cell line, Aag2. Here, we investigated the effect of Wolbachia, a widespread endosymbiont of many insect species, that is known to block replication of several pathogenic flaviviruses, on CFAV. Our results demonstrated that, in mosquito cells, Wolbachia vastly suppresses replication of CFAV, with significantly less CFAV viral interfering small RNAs produced in the cells. However, removal of Wolbachia with tetracycline led to increased CFAV replication. These results suggest that Wolbachia is also able to suppress an ISF.

Transcriptome Analysis Reveals a Diverse Range of Novel Viruses in Australian Sugarcane Soldier Fly (Inopus flavus) Larvae.

In Australia, Soldier flies (Inopus spp.) are economically significant pests of sugarcane that currently lack a viable management strategy. Despite various research efforts, the mechanisms underlying the damage caused by soldier fly larvae remain poorly understood. Our study aims to explore whether this damage is associated with the transmission of plant viruses during larval feeding. We also explore the larval transcriptome to identify any entomopathogenic viruses with the potential to be used as biocontrol agents in future pest management programs. Seven novel virus sequences are identified and characterised using de novo assembly of RNA-Seq data obtained from salivary glands of larvae. The novel virus sequences belong to different virus families and are tentatively named SF-associated anphevirus (SFaAV), SF-associated orthomyxo-like virus (SFaOV), SF-associated narna-like virus (SFaNV), SF-associated partiti-like virus (SFaPV), SF-associated toti-like virus (SFaTV-1 and SFaTV-2) and SF-associated densovirus (SFaDV). These newly identified viruses are more likely insect-associated viruses, as phylogenetic analyses show that they cluster with other insect-specific viruses. Small RNA analysis indicates prominent peaks at both 21 nt and 26-29 nt, suggesting the activation of host siRNA and piwiRNA pathways. Our study helps to improve understanding of the virome of soldier flies and could identify insect viruses for deployment in novel pest management strategies.

Conserved microRNA miR-8 blocks activation of the Toll pathway by upregulating Serpin 27 transcripts.

microRNAs (miRNAs) play significant regulatory roles in gene expression at the post-transcriptional level. This includes modulating processes such as development, immunity, cancer, and host-pathogen interactions. It was recently shown that the phylogenetically deeply conserved miRNA, miR-8, plays a role in maintaining the homeostasis of immunity by suppressing the production of anti-microbial peptides. In this study, we show that miR-8 from the insect Plutella xylostella positively regulates the transcript levels of the serine protease inhibitor Serpin 27, which has been shown to regulate activation of the Toll pathway and prophenoloxidase involved in the melanization response in insects. Interestingly, miR-8 is downregulated following parasitization by Diadegma semiclausum leading to significant declines in Serpin 27 transcript levels. This allows upregulation of antimicrobial peptides, such as gloverin, that are controlled by the Toll pathway and activation of proteolytic cascades essential for humoral immune responses to foreign invasion.

Analysing inhibition of dengue virus in Wolbachia-infected mosquito cells following the removal of Wolbachia.

Wolbachia pipientis is known to block replication of positive sense RNA viruses. Previously, we created an Aedes aegypti Aag2 cell line (Aag2.wAlbB) transinfected with the wAlbB strain of Wolbachia and a matching tetracycline-cured Aag2.tet cell line. While dengue virus (DENV) was blocked in Aag2.wAlbB cells, we found significant inhibition of DENV in Aag2.tet cells. RNA-Seq analysis of the cells confirmed removal of Wolbachia and lack of expression of Wolbachia genes that could have been due to lateral gene transfer in Aag2.tet cells. However, we noticed a substantial increase in the abundance of phasi charoen-like virus (PCLV) in Aag2.tet cells. When RNAi was used to reduce the PCLV levels, DENV replication was significantly increased. Further, we found significant changes in the expression of antiviral and proviral genes in Aag2.tet cells. Overall, the results reveal an antagonistic interaction between DENV and PCLV and how PCLV-induced changes could contribute to DENV inhibition.

Sex dependent transcriptome responses of the diamondback moth, Plutella xylostella L. to cold stress.

Temperature has fundamental influences on the performance and distribution of insects. While considerable attention has been devoted to extreme conditions, particularly extreme cold conditions, few studies have investigated effects of mild cold conditions on insects. We examined the transcriptomic changes in mid-fourth instar larvae of both sexes reared at 10 °C and 25 °C to investigate sex-dependent responses of Plutella xylostella to mild cold stress. There were 624 differentially expressed genes (DEGs) in females, the majority of which (n = 386) were down-regulated. In males 3239 genes were differentially expressed and the majority (n = 2341) were up-regulated. Only 280 DEGs were common to both sexes. In females, there were no DEGs encoding heat shock or cold shock proteins, but six of these DEGs were found in males. These differences suggest that females and males might adopt some different strategies to cope with cold stress and/or that they were affected by rearing under cold conditions to different degrees and in different ways. In addition, DEGs encoding antimicrobial peptides, cytochrome P450 monooxygenases, fatty acid-related enzymes, cuticle proteins, myofilament, and hormone-related proteins were found in both sexes under cold stress. The transcriptome study reveals unexpected sex-dependent thermal responses and provides new information of how an insect that does not diapause copes with low temperatures.

Aedes aegypti gut transcriptomes respond differently to microbiome transplants from field-caught or laboratory-reared mosquitoes.

The mosquito microbiome is critical for host development and plays a major role in many aspects of mosquito biology. While the microbiome is commonly dominated by a small number of genera, there is considerable variation in composition among mosquito species, life stages, and geography. How the host controls and is affected by this variation is unclear. Using microbiome transplant experiments, we asked whether there were differences in transcriptional responses when mosquitoes of different species were used as microbiome donors. We used microbiomes from four different donor species spanning the phylogenetic breadth of the Culicidae, collected either from the laboratory or field. We found that when recipients received a microbiome from a donor reared in the laboratory, the response was remarkably similar regardless of donor species. However, when the donor had been collected from the field, far more genes were differentially expressed. We also found that while the transplant procedure did have some effect on the host transcriptome, this is likely to have had a limited effect on mosquito fitness. Overall, our results highlight the possibility that variation in mosquito microbiome communities are associated with variability in host-microbiome interactions and further demonstrate the utility of the microbiome transplantation technique.

N6-methyladenosine modification of the Aedes aegypti transcriptome and its alteration upon dengue virus infection in Aag2 cell line.

The N6-methyladenosine (m6A) modification of RNA has been reported to affect viral infections. Studies have confirmed the role of m6A in replication of several vector-borne flaviviruses, including dengue virus (DENV), in mammalian cells. Here, we explored the role of m6A in DENV replication in the mosquito Aedes aegypti Aag2 cell line. We first determined the presence of m6A on the RNAs from mosquito cells and using methylated RNA immunoprecipitation and sequencing (MeRIP-Seq) identified m6A modification of the mosquito transcriptome and those that changed upon DENV infection. Depletion of m6A methyltransferases and the m6A binding protein YTHDF3 RNAs decreased the replication of DENV. In particular, we found that the Ae. aegypti ubiquitin carrier protein 9 (Ubc9) is m6A modified and its expression increases after DENV infection. Silencing of the gene and ectopic expression of Ubc9 led to reduced and increased DENV replication, respectively. The abundance of Ubc9 mRNA and its stability were reduced with the inhibition of m6A modification, implying that m6A modification of Ubc9 might enhance expression of the gene. We also show that the genome of DENV is m6A modified at five sites in mosquito cells. Altogether, this work reveals the involvement of m6A modification in Ae. aegypti-DENV interaction.

Discovery of a Novel Jingmenvirus in Australian Sugarcane Soldier Fly (Inopus flavus) Larvae.

In Australia, soldier flies are major pests of sugarcane, and they can cause significant yield losses in some areas, possibly due to the virus' transmission to the plants. We sequenced fly larvae salivary glands and identified a novel jingmenvirus, putatively named Inopus flavus jingmenvirus 1 (IFJV1). Phylogenetic trees confirmed that IFJV1 groups with insect-associated jingmenviruses, newly identified flavivirus-like viruses with a segmented genome. After the design and the validation of molecular detection systems for IFJV1, larval homogenates were passaged on insect and vertebrate cells, but IFJV1 could only be detected in the first two passages in insect cells and not at all in vertebrate cells. Despite this lack of consistent replication in laboratory models, this virus does replicate in its host Inopus flavus, as sequenced, small RNA from the larvae matched the IFJV1 sequences. Moreover, they were found to be predominantly 21 nucleotides long and map to the whole sequences on both strands, which is typical of an actively replicating virus. This discovery confirms the worldwide presence of jingmenviruses which, until now, had only been detected on four continents. However, the study of IFJV1 tropism and the possible pathogenicity to its host or the sugarcane it parasitizes requires the development of a stable replication model.

Transcriptomics Reveal Several Novel Viruses from Canegrubs (Coleoptera: Scarabaeidae) in Central Queensland, Australia.

Canegrubs (Coleoptera: Scarabaeidae) are major pests of sugarcane crops in Australia, but despite long-term and intensive research, no commercially viable biological control agents have been identified. We used the RNA-Seq approach to explore the viriomes of three different species of canegrubs from central Queensland, Australia to identify potential candidates for biological control. We identified six novel RNA viruses, characterized their genomes, and inferred their evolutionary relationships with other closely related viruses. These novel viruses showed similarity to other known members from picornaviruses, benyviruses, sobemoviruses, totiviruses, and reoviruses. The abundance of viral reads varied in these libraries; for example, Dermolepida albohirtum picorna-like virus (9696 nt) was built from 83,894 assembled reads while only 1350 reads mapped to Lepidiota negatoria beny-like virus (6371 nt). Future studies are essential to determine their natural incidence in different life stages of the host, biodiversity, geographical distributions, and potential as biological control agents for these important pests of sugarcane.

Exploration of RNA-Seq data to identify a potential pathogen of the leaf-mining moth, Stomphastis thraustica (Meyrick, 1908) (Lepidoptera: Gracillariidae).

The leaf-mining moth, Stomphastis thraustica (Meyrick, 1908) was imported to Australia as a potential biological control agent of an exotic weed, bellyache bush (Jatropha gossypiifolia), from Peru. The insect colony has been maintained in the quarantine facility for over eight years but recently, significant mortality was observed in the culture. The larvae demonstrated swollen intersegments with a fragile integument. The infected larvae are cloudy muted green or yellowish whereas a healthy late instar larva is a vivid green. They slowly dehydrate and eventually die, at which point the larval body becomes rubbery and turns to black. We used next generation sequencing to identify the cause of mortality in the insects. Total RNA was extracted from 20 larvae in two cohorts, one with and one without apparent symptoms of disease, for deep sequencing on NovaSeq platform after eukaryote ribosomal RNA depletion. We identified several non-insect sequences belonging to viruses, bacteria, and fungi, but none of those showed significant abundance or enrichment in the infected dataset. The sequences related to a unicellular yeast, Saccharomyces cerevisiae, and they were among the highly expressed non-insect contigs; more than 5% of reads in both libraries mapped to the genome of this opportunistic microorganism.

Differential temperature responses between Plutella xylostella and its specialist endo-larval parasitoid Diadegma semiclausum-Implications for biological control.

Understanding the thermal dynamics of host-parasitoid interactions is crucial to predicting how biological control of pest insects by parasitoids might be affected by geographic location and climate change. We compared performance traits of Plutella xylostella (Lepidoptera: Plutellidae) and its solitary endo-larval parasitoid Diadegma semiclausum (Hymenoptera: Ichneumonidae), over a wide range of constant rearing temperatures (10-30°C). Parasitoids reared at 30°C experienced reductions in pupation rate, pupal mass, egg load, and adult life span when compared with those reared at lower temperatures. Our analyses of the fate of parasitoids and their hosts and intergenerational population growth at different rearing temperatures show that D. semiclausum and P. xylostella respond differently to temperature, leading to divergent outcomes under different temperature conditions. Some parasitoid larvae could not complete development at 30°C, the temperature at which the host biomass was least and the metabolic demands of the parasitoid could be high, suggesting that parasitoid development might be constrained by lack of host resources at higher temperatures. We discuss the potential mechanisms of parasitoid susceptibility to elevated temperatures, which likely explain the pronounced seasonal dynamics of D. semiclausum in subtropical regions and its failure to establish in lowland tropical regions, where P. xylostella is a serious pest. Similar interactions in other host-parasitoid associations would constrain the efficacy of parasitoids as biological control agents as global temperatures increase.

Corrigendum to ``Examination of population genetics of the Coconut Rhinoceros Beetle (Oryctes rhinoceros) and the incidence of its biocontrol agent (Oryctes rhinoceros nudivirus) in the South Pacific Islands Current'' [Research in Insect Science 1 (2021) 100015].

[This corrects the article DOI: 10.1016/j.cris.2021.100015.].