Discovery and evolution of bunyavirids in arctic phantom midges and ancient bunyavirid-like sequences in insect genomes.

UNLABELLED: Bunyaviridae is a large family of RNA viruses chiefly comprised of vertebrate and plant pathogens. We discovered novel bunyavirids that are approximately equally divergent from each of the five known genera. We characterized novel genome sequences for two bunyavirids, namely, Kigluaik phantom virus (KIGV), from tundra-native phantom midges (Chaoborus), and Nome phantom virus (NOMV), from tundra-invading phantom midges, and demonstrated that these bunyavirid-like sequences belong to an infectious virus by passaging KIGV in mosquito cell culture, although the infection does not seem to be well sustained beyond a few passages. Virus and host gene sequences from individuals collected on opposite ends of North America, a region spanning 4,000 km, support a long-term, vertically transmitted infection of KIGV in Chaoborus trivittatus. KIGV-like sequences ranging from single genes to full genomes are present in transcriptomes and genomes of insects belonging to six taxonomic orders, suggesting an ancient association of this clade with insect hosts. In Drosophila, endogenous virus genes have been coopted, forming an orthologous tandem gene family that has been maintained by selection during the radiation of the host genus. Our findings indicate that bunyavirid-host interactions in nonbloodsucking arthropods have been much more extensive than previously thought. IMPORTANCE: Very little is known about the viral diversity in polar freshwater ponds, and perhaps less is known about the effects that climate-induced habitat changes in these regions will have on virus-host interactions in the coming years. Our results show that at the tundra-boreal boundary, a hidden viral landscape is being altered as infected boreal phantom midges colonize tundra ponds. Likewise, relatively little is known of the deeper evolutionary history of bunyavirids that has led to the stark lifestyle contrasts between some genera. The discovery of this novel bunyavirid group suggests that ancient and highly divergent bunyavirid lineages remain undetected in nature and may offer fresh insight into host reservoirs, potential sources of emerging disease, and major lifestyle shifts in the evolutionary history of viruses in the family Bunyaviridae.

A novel RNA binding protein affects rbcL gene expression and is specific to bundle sheath chloroplasts in C4 plants.

BACKGROUND: Plants that utilize the highly efficient C4 pathway of photosynthesis typically possess kranz-type leaf anatomy that consists of two morphologically and functionally distinct photosynthetic cell types, the bundle sheath (BS) and mesophyll (M) cells. These two cell types differentially express many genes that are required for C4 capability and function. In mature C4 leaves, the plastidic rbcL gene, encoding the large subunit of the primary CO2 fixation enzyme Rubisco, is expressed specifically within BS cells. Numerous studies have demonstrated that BS-specific rbcL gene expression is regulated predominantly at post-transcriptional levels, through the control of translation and mRNA stability. The identification of regulatory factors associated with C4 patterns of rbcL gene expression has been an elusive goal for many years. RESULTS: RLSB, encoded by the nuclear RLSB gene, is an S1-domain RNA binding protein purified from C4 chloroplasts based on its specific binding to plastid-encoded rbcL mRNA in vitro. Co-localized with LSU to chloroplasts, RLSB is highly conserved across many plant species. Most significantly, RLSB localizes specifically to leaf bundle sheath (BS) cells in C4 plants. Comparative analysis using maize (C4) and Arabidopsis (C3) reveals its tight association with rbcL gene expression in both plants. Reduced RLSB expression (through insertion mutation or RNA silencing, respectively) led to reductions in rbcL mRNA accumulation and LSU production. Additional developmental effects, such as virescent/yellow leaves, were likely associated with decreased photosynthetic function and disruption of associated signaling networks. CONCLUSIONS: Reductions in RLSB expression, due to insertion mutation or gene silencing, are strictly correlated with reductions in rbcL gene expression in both maize and Arabidopsis. In both plants, accumulation of rbcL mRNA as well as synthesis of LSU protein were affected. These findings suggest that specific accumulation and binding of the RLSB binding protein to rbcL mRNA within BS chloroplasts may be one determinant leading to the characteristic cell type-specific localization of Rubisco in C4 plants. Evolutionary modification of RLSB expression, from a C3 "default" state to BS cell-specificity, could represent one mechanism by which rbcL expression has become restricted to only one cell type in C4 plants.

Phylogeny, integration and expression of sigma virus-like genes in Drosophila.

The recent and surprising discovery of widespread NIRVs (non-retroviral integrated RNA viruses) has highlighted the importance of genomic interactions between non-retroviral RNA viruses and their eukaryotic hosts. Among the viruses with integrated representatives are the rhabdoviruses, a family of negative sense single-stranded RNA viruses. We identify sigma virus-like NIRVs of Drosophila spp. that represent unique cases where NIRVs are closely related to exogenous RNA viruses in a model host organism. We have used a combination of bioinformatics and laboratory methods to explore the evolution and expression of sigma virus-like NIRVs in Drosophila. Recent integrations in Drosophila provide a promising experimental system to study functionality of NIRVs. Moreover, the genomic architecture of recent NIRVs provides an unusual evolutionary window on the integration mechanism. For example, we found that a sigma virus-like polymerase associated protein (P) gene appears to have been integrated by template switching of the blastopia-like LTR retrotransposon. The sigma virus P-like NIRV is present in multiple retroelement fused open reading frames on the X and 3R chromosomes of Drosophila yakuba - the X-linked copy is transcribed to produce an RNA product in adult flies. We present the first account of sigma virus-like NIRVs and the first example of NIRV expression in a model animal system, and therefore provide a platform for further study of the possible functions of NIRVs in animal hosts.

Evolutionary maintenance of filovirus-like genes in bat genomes.

BACKGROUND: Little is known of the biological significance and evolutionary maintenance of integrated non-retroviral RNA virus genes in eukaryotic host genomes. Here, we isolated novel filovirus-like genes from bat genomes and tested for evolutionary maintenance. We also estimated the age of filovirus VP35-like gene integrations and tested the phylogenetic hypotheses that there is a eutherian mammal clade and a marsupial/ebolavirus/Marburgvirus dichotomy for filoviruses. RESULTS: We detected homologous copies of VP35-like and NP-like gene integrations in both Old World and New World species of Myotis (bats). We also detected previously unknown VP35-like genes in rodents that are positionally homologous. Comprehensive phylogenetic estimates for filovirus NP-like and VP35-like loci support two main clades with a marsupial and a rodent grouping within the ebolavirus/Lloviu virus/Marburgvirus clade. The concordance of VP35-like, NP-like and mitochondrial gene trees with the expected species tree supports the notion that the copies we examined are orthologs that predate the global spread and radiation of the genus Myotis. Parametric simulations were consistent with selective maintenance for the open reading frame (ORF) of VP35-like genes in Myotis. The ORF of the filovirus-like VP35 gene has been maintained in bat genomes for an estimated 13. 4 MY. ORFs were disrupted for the NP-like genes in Myotis. Likelihood ratio tests revealed that a model that accommodates positive selection is a significantly better fit to the data than a model that does not allow for positive selection for VP35-like sequences. Moreover, site-by-site analysis of selection using two methods indicated at least 25 sites in the VP35-like alignment are under positive selection in Myotis. CONCLUSIONS: Our results indicate that filovirus-like elements have significance beyond genomic imprints of prior infection. That is, there appears to be, or have been, functionally maintained copies of such genes in mammals. "Living fossils" of filoviruses appear to be selectively maintained in a diverse mammalian genus (Myotis).

Cellular production of a counterfeit viral protein confers immunity to infection by a related virus.

DNA copies of many non-retroviral RNA virus genes or portions thereof (NIRVs) are present in the nuclear genomes of many eukaryotes. These have often been preserved for millions of years of evolution, suggesting that they play an important cellular function. One possible function is resistance to infection by related viruses. In some cases, this appears to occur through the piRNA system, but in others by way of counterfeit viral proteins encoded by NIRVs. In the fungi, NIRVs may be as long as 1,400 uninterrupted codons. In one such case in the yeast Debaryomyces hansenii, one of these genes provides immunity to a related virus by virtue of expression of a counterfeit viral capsid protein, which interferes with assembly of viral capsids by negative complementation. The widespread occurrence of non-retroviral RNA virus genes in eukaryotes may reflect an underappreciated method of host resistance to infection. This work demonstrates for the first time that an endogenous host protein encoded by a gene that has been naturally acquired from a virus and fixed in a eukaryote can interfere with the replication of a related virus and do so by negative complementation.

A structural and primary sequence comparison of the viral RNA-dependent RNA polymerases.

A systematic bioinformatic approach to identifying the evolutionarily conserved regions of proteins has verified the universality of a newly described conserved motif in RNA-dependent RNA polymerases (motif F). In combination with structural comparisons, this approach has defined two regions that may be involved in unwinding double-stranded RNA (dsRNA) for transcription. One of these is the N-terminal portion of motif F and the second is a large insertion in motif F present in the RNA-dependent RNA polymerases of some dsRNA viruses.

Widespread mitovirus sequences in plant genomes.

The exploration of the evolution of RNA viruses has been aided recently by the discovery of copies of fragments or complete genomes of non-retroviral RNA viruses (Non-retroviral Endogenous RNA Viral Elements, or NERVEs) in many eukaryotic nuclear genomes. Among the most prominent NERVEs are partial copies of the RNA dependent RNA polymerase (RdRP) of the mitoviruses in plant mitochondrial genomes. Mitoviruses are in the family Narnaviridae, which are the simplest viruses, encoding only a single protein (the RdRP) in their unencapsidated viral plus strand. Narnaviruses are known only in fungi, and the origin of plant mitochondrial mitovirus NERVEs appears to be horizontal transfer from plant pathogenic fungi. At least one mitochondrial mitovirus NERVE, but not its nuclear copy, is expressed.

Evidence that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the Miocene.

An understanding of the timescale of evolution is critical for comparative virology but remains elusive for many RNA viruses. Age estimates based on mutation rates can severely underestimate divergences for ancient viral genes that are evolving under strong purifying selection. Paleoviral dating, however, can provide minimum age estimates for ancient divergence, but few orthologous paleoviruses are known within clades of extant viruses. For example, ebolaviruses and marburgviruses are well-studied mammalian pathogens, but their comparative biology is difficult to interpret because the existing estimates of divergence are controversial. Here we provide evidence that paleoviral elements of two genes (ebolavirus-like VP35 and NP) in cricetid rodent genomes originated after the divergence of ebolaviruses and cuevaviruses from marburgviruses. We provide evidence of orthology by identifying common paleoviral insertion sites among the rodent genomes. Our findings indicate that ebolaviruses and cuevaviruses have been diverging from marburgviruses since the early Miocene.

Selectively maintained paleoviruses in Holarctic water fleas reveal an ancient origin for phleboviruses.

The ecological model, Daphnia pulex (Cladocera: Daphniidae), is broadly distributed in Holarctic freshwater habitats and has been the subject of multidisciplinary study for over half a century, but never has a natural RNA virus infection been reported in daphnids. Here we report on a group of paleoviruses related to RNA dependent RNA polymerase in the genome of D. pulex. Phylogenetic analysis suggests that these paleoviruses are derived from a viral lineage within the genus Phlebovirus. Comparison of the genomic sequences flanking individual paleoviruses reveal that some are orthologous viral insertions having been present in the common ancestor of the D. pulex species complex, which is millions of years old. Still, we detected some sites that have the signature of purifying selection. In contrast, other paleoviruses in this group seem to be unique to specific host lineages and even contain undisrupted open reading frames, suggesting either more recent acquisition, or selective maintenance.

Virus-host co-evolution under a modified nuclear genetic code.

Among eukaryotes with modified nuclear genetic codes, viruses are unknown. However, here we provide evidence of an RNA virus that infects a fungal host (Scheffersomyces segobiensis) with a derived nuclear genetic code where CUG codes for serine. The genomic architecture and phylogeny are consistent with infection by a double-stranded RNA virus of the genus Totivirus. We provide evidence of past or present infection with totiviruses in five species of yeasts with modified genetic codes. All but one of the CUG codons in the viral genome have been eliminated, suggesting that avoidance of the modified codon was important to viral adaptation. Our mass spectroscopy analysis indicates that a congener of the host species has co-opted and expresses a capsid gene from totiviruses as a cellular protein. Viral avoidance of the host's modified codon and host co-option of a protein from totiviruses suggest that RNA viruses co-evolved with yeasts that underwent a major evolutionary transition from the standard genetic code.