Many illnesses, including influenza, hemorrhagic fever, and rabies, are caused by a group of viruses called negative-sense RNA viruses. The genetic information—or genome—of these viruses is encoded in strands of RNA that must be copied before they can be translated into the proteins needed to build new viruses. It is currently known that there are at least eight different families of these viruses, which have a wide range of shapes and sizes and arrange their RNA in different ways.

Insects, spiders, and other arthropods carry many different RNA viruses. Many of these viruses have not previously been studied, and those that have been studied so far are mainly those that cause diseases in humans and other vertebrates. Researchers therefore only know a limited amount about the diversity of the negative-sense RNA viruses that arthropods harbor and how these viruses evolved. Studying how viruses evolve helps scientists to understand what makes some viruses deadly and others harmless and can also help develop treatments or vaccines for the diseases caused by the viruses.

Li, Shi, Tian, Lin, Kang et al. collected 70 species of insects, spiders, centipedes, and other arthropods in China and sequenced all the negative-sense RNA viruses in the creatures. This revealed an enormous number of negative-sense RNA viruses, including 112 new viruses. Many of the newly discovered arthropod viruses appear to be the ancestors of disease-causing viruses, including influenza viruses and the filoviruses—the group that includes the Ebola virus. Indeed, it appears that arthropods host many—if not all—of the negative-sense RNA viruses that cause disease in vertebrates and plants.

While documenting the new RNA viruses and how they are related to each other, Li et al. found many different genome structures. Some genomes were segmented, which may play an important role in evolution as segments can be easily swapped to create new genetic combinations. Non-segmented and circular genomes were also found. This genetic diversity suggests that arthropods are likely to have played a key role in the evolution of new viruses by acting as a site where many different viruses can interact and exchange genetic information.

Negative-sense RNA viruses are important pathogens that cause a variety of diseases in humans including influenza, hemorrhagic fever, encephalitis, and rabies. Taxonomically, those negative-sense RNA viruses described to date comprise at least eight virus families and four unassigned genera or species (King et al., 2012). Although they share (i) a homologous RNA-dependent RNA polymerase (RdRp), (ii) inverted complementary genome ends, and (iii) an encapsidated negative-sense RNA genome, these viruses display substantial diversity in terms of virion morphology and genome organization (King et al., 2012). One key aspect of genome organization is the number of distinct segments, which is also central to virus classification. Among negative-sense RNA viruses, the number of segments varies from one (order Mononegavirales; unsegmented) to two (family Arenaviridae), three (Bunyaviridae), three-to-four (Ophioviridae), and six-to-eight (Orthomyxoviridae) and is further complicated by differences in the number, structure, and arrangement of the encoded genes.

Despite their diversity and importance in infectious disease, the origins and evolutionary history of the negative-sense RNA viruses are largely obscure. Arthropods harbor a diverse range of RNA viruses, which are often divergent from those that infect vertebrates (Marklewitz et al., 2011, 2013; Cook et al., 2013; Ballinger et al., 2014; Qin et al., 2014; Tokarz et al., 2014a, 2014b). However, those arthropod viruses sampled to date are generally those that have a relationship with vertebrates or are known to be agents of disease (Junglen and Drosten, 2013). To determine the extent of viral diversity harbored by arthropods, as well as their evolutionary history, we performed a systematic survey of negative-sense RNA viruses using RNA sequencing (RNA-seq) on a wide range of arthropods.

Our study suggests that arthropods are major reservoir hosts for many, if not all, of the negative-sense RNA viruses in vertebrates and plants, and hence have likely played a major role in their evolution. This is further supported by the high abundance of viral RNA in the arthropod transcriptome, as well as by the high frequencies of endogenous copies of these viruses in the arthropod genome, greatly expanding the known biodiversity of these genomic ‘fossils’ (Katzourakis and Gifford, 2010; Cui and Holmes, 2012). The often basal position of the arthropod viruses in our phylogenetic trees is also compatible with the idea that the negative-sense RNA viruses found in vertebrates and plants ultimately have their ancestry in arthropods, although this will only be confirmed with a far wider sample of virus biodiversity.

The rich genetic and phylogenetic diversity of arthropod RNA viruses may in part reflect the enormous species number and diversity of arthropods, and that they sometimes live in large and very dense populations that provide abundant hosts to fuel virus transmission. Furthermore, arthropods are involved in almost all ecological guilds and actively interact with other eukaryotes, including animals, plants, and fungi, such that it is possible that they serve as both sources and sinks for viruses present in the environment. In addition, not only were diverse viruses present, but they were often highly abundant. For example, in the pool containing 12 individuals (representing two species) from the Gerridae (Water striders) collected at the same site, we identified at least five negative-sense RNA viruses whose TPM values are well above 100, and where the viral RNA collectively made up more than 50% of the host total RNA (rRNA excluded). Determining why arthropods are able to carry such a large viral diversity and at such frequencies clearly merits further investigation.

The viruses discovered here also exhibited a huge variation in level of abundance. It is possible that this variation is in part due to the stage or severity of infection in individual viruses and may be significantly influenced by the process of pooling, since most of our libraries contain an uneven mixture of different host species or even genera. In addition, it is possible that some low abundance viruses may in fact be derived from other eukaryotic organisms present in the host sampled, such as undigested food or prey, gut micro flora, and parasites. Nevertheless, since the majority of the low abundance viruses appear in the same groups as the highly abundant ones in our phylogenetic analyses, these viruses are most likely associated with arthropods.

Viral infections in vertebrates and plants can be divided into two main categories: (i) arthropod-dependent infections, in which there is spill-over to non-arthropods but where continued virus transmission still requires arthropods, and (ii) arthropod-independent infections, in which the virus has shifted its host range to circulate among vertebrates exclusively (Figure 12). The first category of infections is often associated with major vector-borne diseases (Zhang et al., 2011, 2012). Given the biodiversity of arthropod viruses documented here, it seems likely that arthropod-independent viruses were ultimately derived from arthropod-dependent infections, with subsequent adaptation to vertebrate-only transmission (Figure 12).

Figure 12

Download asset Open asset

One of the most notable discoveries was that of a novel family, the Chuviridae. The identification of this diverse virus family provides a new perspective on the evolutionary origins of segmented and unsegmented viruses. In particular, the chuviruses occupy a phylogenetic position that is in some sense ‘intermediate’ between the segmented and unsegmented negative-sense RNA viruses and display genomic features of both. Indeed, our phylogenetic analysis reveals that genome segmentation has evolved multiple times within the diversity of chuviruses (Figure 7), such that this trait appears to be more flexible than previously anticipated. In addition, the majority of the chuviruses possess circular genomes. To date, the only known circular RNA virus is (hepatitis) deltavirus, although this potentially originated from the human genome (Salehi-Ashtiani et al., 2006) and requires hepatitis B virus for successful replication. As such, the chuviruses may represent the first report of autonomously replicating circular RNA viruses, which opens up an important line of future research.

Our results also provide insights into the evolution of genome segmentation. Within the bunya-arena-like viruses (Figures 3C, 4), the three-segment structure is the most common, with the viral polymerase, nucleoprotein, and surface glycoproteins present on different segments. Notably, our phylogenetic analysis seemingly revealed independent occurrences of both increasing (Tenuivirus and Emaravirus) and decreasing (Arenavirus) segment numbers from the three-segment form (Figure 4). Independent changes of genome segmentation numbers are also observed in the mononegavirales-like viruses (Figure 4) and, more frequently, in the chuviruses (Figure 7A). Consequently, the number of genome segments appears to be a relatively flexible trait at a broad evolutionary scale, although the functional relevance of these changes remains unclear. While the segmented viruses (bunya-arenaviruses, orthomyxoviruses, and ophioviruses) appear to be distinct from the largely unsegmented mononegavirales-like viruses in our phylogenetic analysis, this may be an artifact of under-sampling, especially given that only a tiny fraction of eukaryotes have been sampled to date. With a wider sample of eukaryotic viruses it will be possible to more accurately map changes in segment number onto phylogenetic trees and in so doing come to a more complete understanding of the patterns and determinants of the evolution of genome segmentation.

In sum, our results highlight the remarkable diversity of arthropod viruses. Because arthropods interact with a wide range of organisms including vertebrate animals and plants, they can be seen as the direct or indirect source of many clinically or economically important viruses. The viral genetic and phenotypic diversity documented in arthropods here therefore provides a new perspective on fundamental questions of virus origins, diversity, host range, genome evolution, and disease emergence.

1. 1. Ci-Xiu Li
   2. Mang Shi
   3. Jun-Hua Tian
   4. Xian-Dan Lin
   5. Yan-Jun Kang
   6. Liang-Jun Chen
   7. Xin-Cheng Qin
   8. Jianguo Xu
   9. Edward C Holmes
   10. Yong-Zhen Zhang

   (2014) Virus genome sequences

   ID KM817593-KM817764. Publicly available at NCBI GenBank (http://www.ncbi.nlm.nih.gov/genbank/).

   http://www.ncbi.nlm.nih.gov/genbank/
2. 1. Ci-Xiu Li
   2. Mang Shi
   3. Jun-Hua Tian
   4. Xian-Dan Lin
   5. Yan-Jun Kang
   6. Liang-Jun Chen
   7. Xin-Cheng Qin
   8. Jianguo Xu
   9. Edward C Holmes
   10. Yong-Zhen Zhang

   (2015) Sequence Reads

   ID SRP051790. Publicly available at NCBI Sequence Read Archive (http://www.ncbi.nlm.nih.gov/sra/).

   http://www.ncbi.nlm.nih.gov/sra/

1. 1. Ballinger MJ
   2. Bruenn JA
   3. Hay J
   4. Czechowski D
   5. Taylor DJ

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

   Journal of Virology 88:8783–8794.

   https://doi.org/10.1128/JVI.00531-14

   • Google Scholar
2. 1. Capella-Gutierrez S
   2. Silla-Martinez JM
   3. Gabaldon T

   (2009) trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses

   Bioinformatics 25:1972–1973.

   https://doi.org/10.1093/bioinformatics/btp348

   • Google Scholar
3. 1. Cook S
   2. Chung BY
   3. Bass D
   4. Moureau G
   5. Tang S
   6. McAlister E
   7. Culverwell CL
   8. Glucksman E
   9. Wang H
   10. Brown TD
   11. Gould EA
   12. Harbach RE
   13. de Lamballerie X
   14. Firth AE

   (2013) Novel virus discovery and genome reconstruction from field RNA samples reveals highly divergent viruses in dipteran hosts

   PLOS ONE 8:e80720.

   https://doi.org/10.1371/journal.pone.0080720

   • Google Scholar
4. 1. Cui J
   2. Holmes EC

   (2012) Endogenous RNA viruses of plants in insect genomes

   Virology 427:77–79.

   https://doi.org/10.1016/j.virol.2012.02.014

   • Google Scholar
5. 1. Falk BW
   2. Tsai JH

   (1998) Biology and molecular biology of viruses in the genus Tenuivirus

   Annual Review of Phytopathology 36:139–163.

   https://doi.org/10.1146/annurev.phyto.36.1.139

   • Google Scholar
6. 1. Folmer O
   2. Black M
   3. Hoeh W
   4. Lutz R
   5. Vrijenhoek R

   (1994)

   DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates

   Molecular Marine Biology and Biotechnology 3:294–299.

   • Google Scholar
7. 1. Grabherr MG
   2. Haas BJ
   3. Yassour M
   4. Levin JZ
   5. Thompson DA
   6. Amit I
   7. Adiconis X
   8. Fan L
   9. Raychowdhury R
   10. Zeng Q
   11. Chen Z
   12. Mauceli E
   13. Hacohen N
   14. Gnirke A
   15. Rhind N
   16. di Palma F
   17. Birren BW
   18. Nusbaum C
   19. Lindblad-Toh K
   20. Friedman N
   21. Regev A

   (2011) Full-length transcriptome assembly from RNA-Seq data without a reference genome

   Nature Biotechnology 29:644–652.

   https://doi.org/10.1038/nbt.1883

   • Google Scholar
8. 1. Guindon S
   2. Gascuel O

   (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood

   Systematic Biology 52:696–704.

   https://doi.org/10.1080/10635150390235520

   • Google Scholar
9. 1. Junglen S
   2. Drosten C

   (2013) Virus discovery and recent insights into virus diversity in arthropods

   Current Opinion in Microbiology 16:507–513.

   https://doi.org/10.1016/j.mib.2013.06.005

   • Google Scholar
10. 1. Katoh K
    2. Standley DM

    (2013) MAFFT multiple sequence alignment software version 7: improvements in performance and usability

    Molecular Biology and Evolution 30:772–780.

    https://doi.org/10.1093/molbev/mst010

    • Google Scholar
11. 1. Katzourakis A
    2. Gifford RJ

    (2010) Endogenous viral elements in animal genomes

    PLOS Genetics 6:e1001191.

    https://doi.org/10.1371/journal.pgen.1001191

    • Google Scholar
12. Book

    1. King AM
    2. Adams MJ
    3. Carstens EB
    4. Lefkowitz EJ

    (2012)

    Virus taxonomy: 9th report of the international committee on taxonomy of viruses

    Elsevier Academic Press.

    • Google Scholar
13. 1. Langmead B
    2. Salzberg SL

    (2012) Fast gapped-read alignment with Bowtie 2

    Nature Methods 9:357–359.

    https://doi.org/10.1038/nmeth.1923

    • Google Scholar
14. 1. Li B
    2. Ruotti V
    3. Stewart RM
    4. Thomson JA
    5. Dewey CN

    (2010) RNA-Seq gene expression estimation with read mapping uncertainty

    Bioinformatics 26:493–500.

    https://doi.org/10.1093/bioinformatics/btp692

    • Google Scholar
15. 1. Machida RJ
    2. Knowlton N

    (2012) PCR primers for metazoan nuclear 18S and 28S ribosomal DNA sequences

    PLOS ONE 7:e46180.

    https://doi.org/10.1371/journal.pone.0046180

    • Google Scholar
16. 1. Marklewitz M
    2. Handrick S
    3. Grasse W
    4. Kurth A
    5. Lukashev A
    6. Drosten C
    7. Ellerbrok H
    8. Leendertz FH
    9. Pauli G
    10. Junglen S

    (2011) Gouleako virus isolated from West African mosquitoes constitutes a proposed novel genus in the family Bunyaviridae

    Journal of Virology 85:9227–9234.

    https://doi.org/10.1128/JVI.00230-11

    • Google Scholar
17. 1. Marklewitz M
    2. Zirkel F
    3. Rwego IB
    4. Heidemann H
    5. Trippner P
    6. Kurth A
    7. Kallies R
    8. Briese T
    9. Lipkin WI
    10. Drosten C
    11. Gillespie TR
    12. Junglen S

    (2013) Discovery of a unique novel clade of mosquito-associated bunyaviruses

    Journal of Virology 87:12850–12865.

    https://doi.org/10.1128/JVI.01862-13

    • Google Scholar
18. 1. Qin XC
    2. Shi M
    3. Tian JH
    4. Lin XD
    5. Gao DY
    6. He JR
    7. Wang JB
    8. Li CX
    9. Kang YJ
    10. Yu B
    11. Zhou DJ
    12. Xu J
    13. Plyusnin A
    14. Holmes EC
    15. Zhang YZ

    (2014) A tick-borne segmented RNA virus contains genome segments derived from unsegmented viral ancestors

    Proceedings of the National Academy of Sciences of USA 111:6744–6749.

    https://doi.org/10.1073/pnas.1324194111

    • Google Scholar
19. 1. Quast C
    2. Pruesse E
    3. Yilmaz P
    4. Gerken J
    5. Schweer T
    6. Yarza P
    7. Peplies J
    8. Glockner FO

    (2013) The SILVA ribosomal RNA gene database project: improved data processing and web-based tools

    Nucleic Acids Research 41:D590–D596.

    https://doi.org/10.1093/nar/gks1219

    • Google Scholar
20. 1. Ronquist F
    2. Huelsenbeck JP

    (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models

    Bioinformatics 19:1572–1574.

    https://doi.org/10.1093/bioinformatics/btg180

    • Google Scholar
21. 1. Salehi-Ashtiani K
    2. Luptak A
    3. Litovchick A
    4. Szostak JW

    (2006) A genomewide search for ribozymes reveals an HDV-like sequence in the human CPEB3 gene

    Science 313:1788–1792.

    https://doi.org/10.1126/science.1129308

    • Google Scholar
22. 1. Tokarz R
    2. Sameroff S
    3. Leon MS
    4. Jain K
    5. Lipkin WI

    (2014a) Genome characterization of Long Island tick rhabdovirus, a new virus identified in Amblyomma americanum ticks

    Virology Journal 11:26.

    https://doi.org/10.1186/1743-422X-11-26

    • Google Scholar
23. 1. Tokarz R
    2. Williams SH
    3. Sameroff S
    4. Sanchez Leon M
    5. Jain K
    6. Lipkin WI

    (2014b) Virome analysis of Amblyomma americanum, Dermacentor variabilis, and Ixodes scapularis ticks reveals novel highly divergent vertebrate and invertebrate viruses

    Journal of Virology 88:11480–11492.

    https://doi.org/10.1128/JVI.01858-14

    • Google Scholar
24. 1. Zhang YZ
    2. Zhou DJ
    3. Qin XC
    4. Tian JH
    5. Xiong Y
    6. Wang JB
    7. Chen XP
    8. Gao DY
    9. He YW
    10. Jin D
    11. Sun Q
    12. Guo WP
    13. Wang W
    14. Yu B
    15. Li J
    16. Dai YA
    17. Li W
    18. Peng JS
    19. Zhang GB
    20. Zhang S
    21. Chen XM
    22. Wang Y
    23. Li MH
    24. Lu X
    25. Ye C
    26. de Jong MD
    27. Xu J

    (2012) The ecology, genetic diversity, and phylogeny of Huaiyangshan virus in China

    Journal of Virology 86:2864–2868.

    https://doi.org/10.1128/JVI.06192-11

    • Google Scholar
25. 1. Zhang YZ
    2. Zhou DJ
    3. Xiong Y
    4. Chen XP
    5. He YW
    6. Sun Q
    7. Yu B
    8. Li J
    9. Dai YA
    10. Tian JH
    11. Qin XC
    12. Jin D
    13. Cui Z
    14. Luo XL
    15. Li W
    16. Lu S
    17. Wang W
    18. Peng JS
    19. Guo WP
    20. Li MH
    21. Li ZJ
    22. Zhang S
    23. Chen C
    24. Wang Y
    25. de Jong MD
    26. Xu J

    (2011)

    Hemorrhagic fever caused by a novel tick-borne Bunyavirus in Huaiyangshan, China

    Zhonghua Liu Xing Bing Xue Za Zhi = Zhonghua Liuxingbingxue Zazhi 32:209–220.

    • Google Scholar

Author details

1. Ci-Xiu Li

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

   Contribution

   C-XL, Acquisition of data, Analysis and interpretation of data

   Contributed equally with

   Mang Shi, Jun-Hua Tian, Xian-Dan Lin and Yan-Jun Kang

   Competing interests

   The authors declare that no competing interests exist.

2. Mang Shi

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China
   3. Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Biological Sciences and Sydney Medical School, The University of Sydney, Sydney, Australia

   Contribution

   MS, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

   Contributed equally with

   Ci-Xiu Li, Jun-Hua Tian, Xian-Dan Lin and Yan-Jun Kang

   Competing interests

   The authors declare that no competing interests exist.

3. Jun-Hua Tian

   Wuhan Center for Disease Control and Prevention, Wuhan, China

   Contribution

   J-HT, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Contributed equally with

   Ci-Xiu Li, Mang Shi, Xian-Dan Lin and Yan-Jun Kang

   Competing interests

   The authors declare that no competing interests exist.

4. Xian-Dan Lin

   Wenzhou Center for Disease Control and Prevention, Wenzhou, China

   Contribution

   X-DL, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Contributed equally with

   Ci-Xiu Li, Mang Shi, Jun-Hua Tian and Yan-Jun Kang

   Competing interests

   The authors declare that no competing interests exist.

5. Yan-Jun Kang

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

   Contribution

   Y-JK, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Contributed equally with

   Ci-Xiu Li, Mang Shi, Jun-Hua Tian and Xian-Dan Lin

   Competing interests

   The authors declare that no competing interests exist.

6. Liang-Jun Chen

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

   Contribution

   L-JC, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Competing interests

   The authors declare that no competing interests exist.

7. Xin-Cheng Qin

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

   Contribution

   X-CQ, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Competing interests

   The authors declare that no competing interests exist.

8. Jianguo Xu

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

   Contribution

   JX, Acquisition of data, Analysis and interpretation of data, Contributed unpublished essential data or reagents

   Competing interests

   The authors declare that no competing interests exist.

9. Edward C Holmes

   1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
   2. Marie Bashir Institute for Infectious Diseases and Biosecurity, Charles Perkins Centre, School of Biological Sciences and Sydney Medical School, The University of Sydney, Sydney, Australia

   Contribution

   ECH, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article

   Competing interests

   The authors declare that no competing interests exist.

10. Yong-Zhen Zhang

    1. State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China
    2. Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China

    Contribution

    Y-ZZ, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article, Contributed unpublished essential data or reagents

    For correspondence

    zhangyongzhen@icdc.cn

    Competing interests

    The authors declare that no competing interests exist.

• 11,002

  views

• 2,638

  downloads

• 610

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.