RNA-mediated interference and reverse transcription control the persistence of RNA viruses in the insect model Drosophila.

How persistent viral infections are established and maintained is widely debated and remains poorly understood. We found here that the persistence of RNA viruses in Drosophila melanogaster was achieved through the combined action of cellular reverse-transcriptase activity and the RNA-mediated interference (RNAi) pathway. Fragments of diverse RNA viruses were reverse-transcribed early during infection, which resulted in DNA forms embedded in retrotransposon sequences. Those virus-retrotransposon DNA chimeras produced transcripts processed by the RNAi machinery, which in turn inhibited viral replication. Conversely, inhibition of reverse transcription hindered the appearance of chimeric DNA and prevented persistence. Our results identify a cooperative function for retrotransposons and antiviral RNAi in the control of lethal acute infection for the establishment of viral persistence.

Histone-derived piRNA biogenesis depends on the ping-pong partners Piwi5 and Ago3 in Aedes aegypti.

The piRNA pathway is of key importance in controlling transposable elements in most animal species. In the vector mosquito Aedes aegypti, the presence of eight PIWI proteins and the accumulation of viral piRNAs upon arbovirus infection suggest additional functions of the piRNA pathway beyond genome defense. To better understand the regulatory potential of this pathway, we analyzed in detail host-derived piRNAs in A. aegypti Aag2 cells. We show that a large repertoire of protein-coding genes and non-retroviral integrated RNA virus elements are processed into genic piRNAs by different combinations of PIWI proteins. Among these, we identify a class of genes that produces piRNAs from coding sequences in an Ago3- and Piwi5-dependent fashion. We demonstrate that the replication-dependent histone gene family is a genic source of ping-pong dependent piRNAs and that histone-derived piRNAs are dynamically expressed throughout the cell cycle, suggesting a role for the piRNA pathway in the regulation of histone gene expression. Moreover, our results establish the Aag2 cell line as an accessible experimental model to study gene-derived piRNAs.

piRNA pathway is not required for antiviral defense in Drosophila melanogaster.

Since its discovery, RNA interference has been identified as involved in many different cellular processes, and as a natural antiviral response in plants, nematodes, and insects. In insects, the small interfering RNA (siRNA) pathway is the major antiviral response. In recent years, the Piwi-interacting RNA (piRNA) pathway also has been implicated in antiviral defense in mosquitoes infected with arboviruses. Using Drosophila melanogaster and an array of viruses that infect the fruit fly acutely or persistently or are vertically transmitted through the germ line, we investigated in detail the extent to which the piRNA pathway contributes to antiviral defense in adult flies. Following virus infection, the survival and viral titers of Piwi, Aubergine, Argonaute-3, and Zucchini mutant flies were similar to those of wild type flies. Using next-generation sequencing of small RNAs from wild type and siRNA mutant flies, we showed that no viral-derived piRNAs were produced in fruit flies during different types of viral infection. Our study provides the first evidence, to our knowledge, that the piRNA pathway does not play a major role in antiviral defense in adult Drosophila and demonstrates that viral-derived piRNA production depends on the biology of the host-virus combination rather than being part of a general antiviral process in insects.

Uncovering the Repertoire of Endogenous Flaviviral Elements in Aedes Mosquito Genomes.

Endogenous viral elements derived from nonretroviral RNA viruses have been described in various animal genomes. Whether they have a biological function, such as host immune protection against related viruses, is a field of intense study. Here, we investigated the repertoire of endogenous flaviviral elements (EFVEs) in Aedes mosquitoes, the vectors of arboviruses such as dengue and chikungunya viruses. Previous studies identified three EFVEs from Aedes albopictus cell lines and one from Aedes aegypti cell lines. However, an in-depth characterization of EFVEs in wild-type mosquito populations and individual mosquitoes in vivo has not been performed. We detected the full-length DNA sequence of the previously described EFVEs and their respective transcripts in several A. albopictus and A. aegypti populations from geographically distinct areas. However, EFVE-derived proteins were not detected by mass spectrometry. Using deep sequencing, we detected the production of PIWI-interacting RNA-like small RNAs, in an antisense orientation, targeting the EFVEs and their flanking regions in vivo The EFVEs were integrated in repetitive regions of the mosquito genomes, and their flanking sequences varied among mosquito populations. We bioinformatically predicted several new EFVEs from a Vietnamese A. albopictus population and observed variation in the occurrence of those elements among mosquitoes. Phylogenetic analysis of an A. aegypti EFVE suggested that it integrated prior to the global expansion of the species and subsequently diverged among and within populations. The findings of this study together reveal the substantial structural and nucleotide diversity of flaviviral integrations in Aedes genomes. Unraveling this diversity will help to elucidate the potential biological function of these EFVEs.IMPORTANCE Endogenous viral elements (EVEs) are whole or partial viral sequences integrated in host genomes. Interestingly, some EVEs have important functions for host fitness and antiviral defense. Because mosquitoes also have EVEs in their genomes, characterizing these EVEs is a prerequisite for their potential use to manipulate the mosquito antiviral response. In the study described here, we focused on EVEs related to the Flavivirus genus, to which dengue and Zika viruses belong, in individual Aedes mosquitoes from geographically distinct areas. We show the existence in vivo of flaviviral EVEs previously identified in mosquito cell lines, and we detected new ones. We show that EVEs have evolved differently in each mosquito population. They produce transcripts and small RNAs but not proteins, suggesting a function at the RNA level. Our study uncovers the diverse repertoire of flaviviral EVEs in Aedes mosquito populations and contributes to an understanding of their role in the host antiviral system.

Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor.

BACKGROUND: Klebsiella pneumoniae strains are pathogenic to animals and humans, in which they are both a frequent cause of nosocomial infections and a re-emerging cause of severe community-acquired infections. K. pneumoniae isolates of the capsular serotype K2 are among the most virulent. In order to identify novel putative virulence factors that may account for the severity of K2 infections, the genome sequence of the K2 reference strain Kp52.145 was determined and compared to two K1 and K2 strains of low virulence and to the reference strains MGH 78578 and NTUH-K2044. RESULTS: In addition to diverse functions related to host colonization and virulence encoded in genomic regions common to the four strains, four genomic islands specific for Kp52.145 were identified. These regions encoded genes for the synthesis of colibactin toxin, a putative cytotoxin outer membrane protein, secretion systems, nucleases and eukaryotic-like proteins. In addition, an insertion within a type VI secretion system locus included sel1 domain containing proteins and a phospholipase D family protein (PLD1). The pld1 mutant was avirulent in a pneumonia model in mouse. The pld1 mRNA was expressed in vivo and the pld1 gene was associated with K. pneumoniae isolates from severe infections. Analysis of lipid composition of a defective E. coli strain complemented with pld1 suggests an involvement of PLD1 in cardiolipin metabolism. CONCLUSIONS: Determination of the complete genome of the K2 reference strain identified several genomic islands comprising putative elements of pathogenicity. The role of PLD1 in pathogenesis was demonstrated for the first time and suggests that lipid metabolism is a novel virulence mechanism of K. pneumoniae.

Complete nucleotide sequence of two multidrug-resistant IncR plasmids from Klebsiella pneumoniae.

We report here the complete nucleotide sequence of two IncR replicons encoding multidrug resistance determinants, including ß-lactam (blaDHA-1, blaSHV-12), aminoglycoside (aphA1, strA, strB), and fluoroquinolone (qnrB4, aac6'-1b-cr) resistance genes. The plasmids have backbones that are similar to each other, including the replication and stability systems, and contain a wide variety of transposable elements carrying known antibiotic resistance genes. This study confirms the increasing clinical importance of IncR replicons as resistance gene carriers.

Limitations in harnessing oral RNA interference as an antiviral strategy in Aedes aegypti.

Mosquitoes, particularly Aedes aegypti, are critical vectors for globally significant pathogenic viruses. This study examines the limitations of oral RNA interference (RNAi) as a strategy to disrupt viral transmission by Ae. aegypti. We hypothesized that double-stranded RNA (dsRNA) targeting the Zika virus (ZIKV) or chikungunya virus (CHIKV) genomes produced by engineered bacterial symbionts could trigger an antiviral response. Mosquitoes mono-colonized with Escherichia coli producing dsZIK or dsCHIK did not display reduced viral titers following exposure to virus-contaminated bloodmeals and failed to generate dsZIK- or dsCHIK-derived small interfering RNAs. To address potential limitations of bacterial dsRNA release, we explored dsRNA inoculation via feeding and injection. Although viral replication was impeded in mosquitoes injected with dsZIK or dsCHIK, no antiviral effect was observed in dsRNA-fed mosquitoes. These findings highlight complexities of implementing oral RNAi as an antiviral strategy in Ae. aegypti and warrant further exploration of local and systemic RNAi mechanisms.

Viral infection disrupts intestinal homeostasis via Sting-dependent NF-κB signaling in Drosophila.

Host-microbe interactions influence intestinal stem cell (ISC) activity to modulate epithelial turnover and composition. Here, we investigated the functional impacts of viral infection on intestinal homeostasis and the mechanisms by which viral infection alters ISC activity. We report that Drosophila A virus (DAV) infection disrupts intestinal homeostasis in Drosophila by inducing sustained ISC proliferation, resulting in intestinal dysplasia, loss of gut barrier function, and reduced lifespan. We found that additional viruses common in laboratory-reared Drosophila also promote ISC proliferation. The mechanism of DAV-induced ISC proliferation involves progenitor-autonomous epidermal growth factor receptor (EGFR) signaling, c-Jun N-terminal kinase (JNK) activity in enterocytes, and requires Sting-dependent nuclear factor κB (NF-κB) (Relish) activity. We further demonstrate that activating Sting-Relish signaling is sufficient to induce ISC proliferation, promote intestinal dysplasia, and reduce lifespan in the absence of infection. Our results reveal that viral infection can significantly disrupt intestinal physiology, highlight a novel role for Sting-Relish signaling, and support a role for viral infection in aging.

Complete nucleotide sequence of the first KPC-2- and SHV-12-encoding IncX plasmid, pKpS90, from Klebsiella pneumoniae.

We report the complete nucleotide sequence of the pKpS90 plasmid, carrying the bla(KPC-2) and bla(SHV-12) genes. This plasmid was isolated from a sequence type 258 (ST258) Klebsiella pneumoniae strain responsible for an outbreak in a French university hospital in 2009. pKpS90 is a 53,286-bp plasmid that belongs to the IncX incompatibility group. pKpS90 consists of a backbone from IncX plasmids, in which the KPC-2-encoding Tn4401 transposon and a bla(SHV-12)-encoding region have been inserted.

Complete nucleotide sequence of the large conjugative pTC2 multireplicon plasmid encoding the VIM-1 metallo-ß-lactamase.

OBJECTIVES: To determine the complete nucleotide sequence of the VIM-1-encoding plasmid pTC2, which was isolated from a Greek Providencia stuartii multiresistant strain. METHODS: The pTC2 plasmid was extracted and sequenced using shotgun and 3 kb paired-end DNA libraries and a 454 sequencing approach. Following assembly into a unique scaffold, gaps were closed by PCR followed by Sanger DNA sequencing. Gene predictions and annotations were performed using the CAAT-Box tool and the Conserved Domain Search service. Sequence comparisons were performed using the Artemis Comparison Tool. RESULTS: Plasmid pTC2 (180,184 bp) was found to be a multireplicon plasmid (IncA/C and IncR), with a large IncA/C backbone and a mosaic multidrug resistance (MDR) region, in which was inserted a 13 kb IncR fragment. Gene annotation allowed the identification of a complete IncA/C-type transfer system and of several putative maintenance modules, both on the IncA/C backbone and on the IncR fragment. The complex MDR region contained 9 insertion sequences (7 IS26, 1 IS1 and 1 IS6100), 10 resistance genes and a mercury resistance operon integrated into unit transposons, composite transposons or integrons. CONCLUSIONS: The pTC2 combines a broad host range, transfer and maintenance capacities, plasticity of the MDR region and a wide variety of resistance genes, properties that may contribute to the spreading of resistance determinants.

Ribosomal RNA (rRNA) sequences from 33 globally distributed mosquito species for improved metagenomics and species identification.

Total RNA sequencing (RNA-seq) is an important tool in the study of mosquitoes and the RNA viruses they vector as it allows assessment of both host and viral RNA in specimens. However, there are two main constraints. First, as with many other species, abundant mosquito ribosomal RNA (rRNA) serves as the predominant template from which sequences are generated, meaning that the desired host and viral templates are sequenced far less. Second, mosquito specimens captured in the field must be correctly identified, in some cases to the sub-species level. Here, we generate mosquito rRNA datasets which will substantially mitigate both of these problems. We describe a strategy to assemble novel rRNA sequences from mosquito specimens and produce an unprecedented dataset of 234 full-length 28S and 18S rRNA sequences of 33 medically important species from countries with known histories of mosquito-borne virus circulation (Cambodia, the Central African Republic, Madagascar, and French Guiana). These sequences will allow both physical and computational removal of rRNA from specimens during RNA-seq protocols. We also assess the utility of rRNA sequences for molecular taxonomy and compare phylogenies constructed using rRNA sequences versus those created using the gold standard for molecular species identification of specimens-the mitochondrial cytochrome c oxidase I (COI) gene. We find that rRNA- and COI-derived phylogenetic trees are incongruent and that 28S and concatenated 28S+18S rRNA phylogenies reflect evolutionary relationships that are more aligned with contemporary mosquito systematics. This significant expansion to the current rRNA reference library for mosquitoes will improve mosquito RNA-seq metagenomics by permitting the optimization of species-specific rRNA depletion protocols for a broader range of species and streamlining species identification by rRNA sequence and phylogenetics.

Multifaceted contributions of Dicer2 to arbovirus transmission by Aedes aegypti.

Arthropod-borne viruses (arboviruses) transmitted by Aedes aegypti mosquitoes are an increasing threat to global health. The small interfering RNA (siRNA) pathway is considered the main antiviral immune pathway of insects, but its effective impact on arbovirus transmission is surprisingly poorly understood. Here, we use CRISPR-Cas9-mediated gene editing in vivo to mutate Dicer2, a gene encoding the RNA sensor and key component of the siRNA pathway. The loss of Dicer2 enhances early viral replication and systemic viral dissemination of four medically significant arboviruses (chikungunya, Mayaro, dengue, and Zika viruses) representing two viral families. However, Dicer2 mutants and wild-type mosquitoes display overall similar levels of vector competence. In addition, Dicer2 mutants undergo significant virus-induced mortality during infection with chikungunya virus. Together, our results define a multifaceted role for Dicer2 in the transmission of arboviruses by Ae. aegypti mosquitoes and pave the way for further mechanistic investigations.

Comparative genomics of protoploid Saccharomycetaceae.

Our knowledge of yeast genomes remains largely dominated by the extensive studies on Saccharomyces cerevisiae and the consequences of its ancestral duplication, leaving the evolution of the entire class of hemiascomycetes only partly explored. We concentrate here on five species of Saccharomycetaceae, a large subdivision of hemiascomycetes, that we call "protoploid" because they diverged from the S. cerevisiae lineage prior to its genome duplication. We determined the complete genome sequences of three of these species: Kluyveromyces (Lachancea) thermotolerans and Saccharomyces (Lachancea) kluyveri (two members of the newly described Lachancea clade), and Zygosaccharomyces rouxii. We included in our comparisons the previously available sequences of Kluyveromyces lactis and Ashbya (Eremothecium) gossypii. Despite their broad evolutionary range and significant individual variations in each lineage, the five protoploid Saccharomycetaceae share a core repertoire of approximately 3300 protein families and a high degree of conserved synteny. Synteny blocks were used to define gene orthology and to infer ancestors. Far from representing minimal genomes without redundancy, the five protoploid yeasts contain numerous copies of paralogous genes, either dispersed or in tandem arrays, that, altogether, constitute a third of each genome. Ancient, conserved paralogs as well as novel, lineage-specific paralogs were identified.

Evaluation of high-throughput sequencing for identifying known and unknown viruses in biological samples.

High-throughput sequencing furnishes a large number of short sequence reads from uncloned DNA and has rapidly become a major tool for identifying viruses in biological samples, and in particular when the target sequence is undefined. In this study, we assessed the analytical sensitivity of a pipeline for detection of viruses in biological samples based on either the Roche-454 genome sequencer or Illumina genome analyzer platforms. We sequenced biological samples artificially spiked with a wide range of viruses with genomes composed of single or double-stranded DNA or RNA, including linear or circular single-stranded DNA. Viruses were added at a very low concentration most often corresponding to 3 or 0.8 times the validated level of detection of quantitative reverse transcriptase PCRs (RT-PCRs). For the viruses represented, or resembling those represented, in public nucleotide sequence databases, we show that the higher output of Illumina is associated with a much greater sensitivity, approaching that of optimized quantitative (RT-)PCRs. In this blind study, identification of viruses was achieved without incorrect identification. Nevertheless, at these low concentrations, the number of reads generated by the Illumina platform was too small to facilitate assembly of contigs without the use of a reference sequence, thus precluding detection of unknown viruses. When the virus load was sufficiently high, de novo assembly permitted the generation of long contigs corresponding to nearly full-length genomes and thus should facilitate the identification of novel viruses.

Interactions of the Insect-Specific Palm Creek Virus with Zika and Chikungunya Viruses in Aedes Mosquitoes.

Palm Creek virus (PCV) is an insect-specific flavivirus that can interfere with the replication of mosquito-borne flaviviruses in Culex mosquitoes, thereby potentially reducing disease transmission. We examined whether PCV could interfere with arbovirus replication in Aedes (Ae.) aegypti and Ae. albopictus mosquitoes, major vectors for many prominent mosquito-borne viral diseases. We infected laboratory colonies of Ae. aegypti and Ae. albopictus with PCV to evaluate infection dynamics. PCV infection was found to persist to at least 21 days post-infection and could be detected in the midguts and ovaries. We then assayed for PCV-arbovirus interference by orally challenging PCV-infected mosquitoes with Zika and chikungunya viruses. For both arboviruses, PCV infection had no effect on infection and transmission rates, indicating limited potential as a method of intervention for Aedes-transmitted arboviruses. We also explored the hypothesis that PCV-arbovirus interference is mediated by the small interfering RNA pathway in silico. Our findings indicate that RNA interference is unlikely to underlie the mechanism of arbovirus inhibition and emphasise the need for empirical examination of individual pairs of insect-specific viruses and arboviruses to fully understand their impact on arbovirus transmission.

Genome-wide computational prediction of tandem gene arrays: application in yeasts.

BACKGROUND: This paper describes an efficient in silico method for detecting tandem gene arrays (TGAs) in fully sequenced and compact genomes such as those of prokaryotes or unicellular eukaryotes. The originality of this method lies in the search of protein sequence similarities in the vicinity of each coding sequence, which allows the prediction of tandem duplicated gene copies independently of their functionality. RESULTS: Applied to nine hemiascomycete yeast genomes, this method predicts that 2% of the genes are involved in TGAs and gene relics are present in 11% of TGAs. The frequency of TGAs with degenerated gene copies means that a significant fraction of tandem duplicated genes follows the birth-and-death model of evolution. A comparison of sequence identity distributions between sets of homologous gene pairs shows that the different copies of tandem arrayed paralogs are less divergent than copies of dispersed paralogs in yeast genomes. It suggests that paralogs included in tandem structures are more recent or more subject to the gene conversion mechanism than other paralogs. CONCLUSION: The method reported here is a useful computational tool to provide a database of TGAs composed of functional or nonfunctional gene copies. Such a database has obvious applications in the fields of structural and comparative genomics. Notably, a detailed study of the TGA catalog will make it possible to tackle the fundamental questions of the origin and evolution of tandem gene clusters.

Genome evolution in yeasts.

Identifying the mechanisms of eukaryotic genome evolution by comparative genomics is often complicated by the multiplicity of events that have taken place throughout the history of individual lineages, leaving only distorted and superimposed traces in the genome of each living organism. The hemiascomycete yeasts, with their compact genomes, similar lifestyle and distinct sexual and physiological properties, provide a unique opportunity to explore such mechanisms. We present here the complete, assembled genome sequences of four yeast species, selected to represent a broad evolutionary range within a single eukaryotic phylum, that after analysis proved to be molecularly as diverse as the entire phylum of chordates. A total of approximately 24,200 novel genes were identified, the translation products of which were classified together with Saccharomyces cerevisiae proteins into about 4,700 families, forming the basis for interspecific comparisons. Analysis of chromosome maps and genome redundancies reveal that the different yeast lineages have evolved through a marked interplay between several distinct molecular mechanisms, including tandem gene repeat formation, segmental duplication, a massive genome duplication and extensive gene loss.

A comparison of random sequence reads versus 16S rDNA sequences for estimating the biodiversity of a metagenomic library.

The construction of metagenomic libraries has permitted the study of microorganisms resistant to isolation and the analysis of 16S rDNA sequences has been used for over two decades to examine bacterial biodiversity. Here, we show that the analysis of random sequence reads (RSRs) instead of 16S is a suitable shortcut to estimate the biodiversity of a bacterial community from metagenomic libraries. We generated 10,010 RSRs from a metagenomic library of microorganisms found in human faecal samples. Then searched them using the program BLASTN against a prokaryotic sequence database to assign a taxon to each RSR. The results were compared with those obtained by screening and analysing the clones containing 16S rDNA sequences in the whole library. We found that the biodiversity observed by RSR analysis is consistent with that obtained by 16S rDNA. We also show that RSRs are suitable to compare the biodiversity between different metagenomic libraries. RSRs can thus provide a good estimate of the biodiversity of a metagenomic library and, as an alternative to 16S, this approach is both faster and cheaper.

Differential Small RNA Responses against Co-Infecting Insect-Specific Viruses in Aedes albopictus Mosquitoes.

The mosquito antiviral response has mainly been studied in the context of arthropod-borne virus (arbovirus) infection in female mosquitoes. However, in nature, both female and male mosquitoes are frequently infected with insect-specific viruses (ISVs). ISVs are capable of infecting the reproductive organs of both sexes and are primarily maintained by vertical transmission. Since the RNA interference (RNAi)-mediated antiviral response plays an important antiviral role in mosquitoes, ISVs constitute a relevant model to study sex-dependent antiviral responses. Using a naturally generated viral stock containing three distinct ISVs, Aedes flavivirus (AEFV), Menghai rhabdovirus (MERV), and Shinobi tetra virus (SHTV), we infected adult Aedes albopictus females and males and generated small RNA libraries from ovaries, testes, and the remainder of the body. Overall, both female and male mosquitoes showed unique small RNA profiles to each co-infecting ISV regardless of the sex or tissue tested. While all three ISVs generated virus-derived siRNAs, only MERV generated virus-derived piRNAs. We also studied the expression of PIWI genes in reproductive tissues and carcasses. In contrast to Piwi5-9, Piwi1-4 were abundantly expressed in ovaries and testes, suggesting that Piwi5-9 are involved in exogenous viral piRNA production. Together, our results show that ISV-infected Aedes albopictus produce viral small RNAs in a virus-specific manner and that male mosquitoes mount a similar small RNA-mediated antiviral response to that of females.

Evidence For Long-Lasting Transgenerational Antiviral Immunity in Insects.

Transgenerational immune priming (TGIP) allows memory-like immune responses to be transmitted from parents to offspring in many invertebrates. Despite increasing evidence for TGIP in insects, the mechanisms involved in the transfer of information remain largely unknown. Here, we show that Drosophila melanogaster and Aedes aegypti transmit antiviral immunological memory to their progeny that lasts throughout generations. We observe that TGIP, which is virus and sequence specific but RNAi independent, is initiated by a single exposure to disparate RNA viruses and also by inoculation of a fragment of viral double-stranded RNA. The progeny, which inherit a viral DNA that is only a fragment of the viral RNA used to infect the parents, display enriched expression of genes related to chromatin and DNA binding. These findings represent a demonstration of TGIP for RNA viruses in invertebrates, broadly increasing our understanding of the immune response, host genome plasticity, and antiviral memory of the germline.