Virome analyses of <Emphasis Type="Italic">Hevea brasiliensis</Emphasis> using small RNA deep sequencing and PCR techniques reveal the presence of a potential new virus

Background

Hevea brasiliensis is an important commercial crop due to the high quality of the latex it produces; however, little is known about viral infections in this plant. The only virus described to infect H. brasiliensis until now is a Carlavirus, which was described more than 30?years ago. Virus-derived small interfering RNA (vsiRNAs) are the product of the plant’s antiviral defense triggered by dsRNA viral intermediates generated, during the replication cycle. These vsiRNAs are complementar to viral genomes and have been widely used to identify and characterize viruses in plants.

Methods

In the present study, we investigated the virome of leaf and sapwood samples from native H. brasiliensis trees collected in two geographic areas in the Brazilian Amazon. Small RNA (sRNA) deep sequencing and bioinformatic tools were used to assembly, identify and characterize viral contigs. Subsequently, PCR amplification techniques were performed to experimentally verify the presence of the viral sequences. Finally, the phylogenetic relationship of the putative new virus with related viral genomes was analyzed.

Results

Our strategy allowed the identification of 32 contigs with high similarity to viral reference genomes, from which 23 exhibited homology to viruses of the Tymoviridae family. The reads showed a predominant size distribution at 21?nt derived from both strands, which was consistent with the vsiRNAs profile. The presence and genome position of the viral contigs were experimentally confirmed using droplet digital PCR amplifications. A 1913 aa long fragment was obtained and used to infer the phylogenetic relationship of the putative new virus, which indicated that it is taxonomically related to the Grapevine fleck virus, genus Maculavirus. The putative new virus was named Hevea brasiliensis virus (HBrV) in reference to its host.

Conclusion

The methodological strategy applied here proved to be efficient in detecting and confirming the presence of new viral sequences on a ‘very difficult to manage’ sample. This is the second time that viral sequences, that could be ascribed as a putative novel virus, associated to the rubber tree has been identified.

     

Mutant swarms of a totivirus‐like entities are present in the red macroalga Chondrus crispus and have been partially transferred to the nuclear genome

Chondrus crispus Stackhouse (Gigartinales) is a red seaweed found on North Atlantic rocky shores. Electrophoresis of RNA extracts showed a prominent band with a size of around 6,000 bp. Sequencing of the band revealed several sequences with similarity to totiviruses, double‐stranded RNA viruses that normally infect fungi. This virus‐like entity was named C. crispus virus (CcV). It should probably be regarded as an extreme viral quasispecies or a mutant swarm since low identity (<65%) was found between sequences. Totiviruses typically code for two genes: one capsid gene (gag) and one RNA‐dependent RNA polymerase gene (pol) with a pseudoknot structure between the genes. Both the genes and the intergenic structures were found in the CcV sequences. A nonidentical gag gene was also found in the nuclear genome of C. crispus, with associated expressed sequence tags (EST) and upstream regulatory features. The gene was presumably horizontally transferred from the virus to the alga. Similar dsRNA bands were seen in extracts from different life cycle stages of C. crispus and from all geographic locations tested. In addition, similar bands were also observed in RNA extractions from other red algae; however, the significance of this apparently widespread phenomenon is unknown. Neither phenotype caused by the infection nor any virus particles or capsid proteins were identified; thus, the presence of viral particles has not been validated. These findings increase the known host range of totiviruses to include marine red algae.     

Role of RNA structure in non-homologous recombination between genomic molecules of brome mosaic virus

Brome mosaic virus (BMV) is a tripartite genome, positive-sense RNA virus of plants. Previously it was demonstrated that local hybridization between BMV RNAs (RNA–RNA heteroduplex formation) efficiently promotes non-homologous RNA recombination. In addition, studies on the role of the BMV polymerase in RNA recombination suggested that the location of non-homologous crossovers depends mostly on RNA structure. As a result, a detailed analysis of a large number of non-homologous recombinants generated in the BMV-based system was undertaken. Recombination hot-spots as well as putative elements in RNA structure enhancing non-homologous crossovers and targeting them in a site-specific manner were identified. To verify these observations the recombinationally active sequence in BMV RNA3 derivative was modified. The results obtained with new RNA3 mutants suggest that the primary and secondary structure of the sequences involved in a heteroduplex formation rather than the length of heteroduplex plays the most important role in the recombination process. The presented data indicate that the sequences proximal to the heteroduplex may also affect template switching by BMV replicase. Moreover, it was shown that both short homologous sequences and a hairpin structure have to accompany a double-stranded region to target non-homologous crossovers in a site-specific manner.     

Characterisation of a novel Emaravirus identified in mosaic-diseased Eurasian aspen (Populus tremula)

Since Emaraviruses have been discovered in 2007 several new species were detected in a range of host plants. Five genome segments of a novel Emaravirus from mosaic-diseased Eurasian aspen (Populus tremula) have been completely determined. The monocistronic, segmented ssRNA genome of the virus shows a genome organisation typical for Emaraviruses encoding the viral RNA-dependent RNA polymerase (RdRP, 268.2 kDa) on RNA1 (7.1 kb), a glycoprotein precursor (GPP, 73.5 kDa) on RNA2 (2.3 kb), the viral nucleocapsid protein (N, 35.6 kDa) on RNA3 (1.6 kb), and a putative movement protein (MP, 41.0 kDa) on RNA4 (1.6 kb). The fifth identified genome segment (RNA5, 1.3 kb) encodes a protein of unknown function (P28, 28.1 kDa). We discovered that it is distantly related to proteins encoded by Emaraviruses, such as P4 of European mountain ash ringspot-associated virus. All proteins from this group contain a central hydrophobic region with a conserved secondary structure and a hydrophobic amino acid stretch, bordered by two highly conserved positions, thus clearly representing a new group of homologues of Emaraviruses. The virus identified in Eurasian aspen is closely associated with observed leaf symptoms, such as mottle, yellow blotching, variegation and chloroses along veins. All five viral RNAs were regularly detectable by RT-PCR in mosaic-diseased P. tremula in Norway, Finland and Sweden (Fennoscandia). Observed symptoms and testing of mosaic-diseased Eurasian aspen by virus-specific RT-PCR targeting RNA3 and RNA4 confirmed a wide geographic distribution of the virus in Fennoscandia. We could demonstrate that the mosaic-disease is graft-transmissible and confirmed that the virus is the causal agent by detection in symptomatic, graft-inoculated seedlings used as rootstocks as well as in the virus-infected scions used for graft-inoculation. Owing to these characteristics, the virus represents a novel species within the genus Emaravirus and was tentatively denominated aspen mosaic-associated virus.     

Identification of two essential aspartates for polymerase activity in parainfluenza virus L protein by a minireplicon system expressing secretory luciferase

Gene expression of nonsegmented negative‐strand RNA viruses (nsNSVs) such as parainfluenza viruses requires the RNA synthesis activity of their polymerase L protein; however, the detailed mechanism of this process is poorly understood. In this study, a parainfluenza minireplicon assay expressing secretory Gaussia luciferase (Gluc) was established to analyze large protein (L) activity. Measurement of Gluc expression in the culture medium of cells transfected with the minigenome and viral polymerase components enabled quick and concise calculation of L activity. By comparing the amino acid sequences in conserved region III (CRIII), a putative polymerase‐active domain of the L protein, two strictly conserved aspartates were identified in all families of nsNSV. A series of L mutants from human parainfluenza virus type 2 and parainfluenza virus type 5 showed that these aspartates are necessary for reporter gene expression. It was also confirmed that these aspartates are important for the production of viral mRNA and antigenome cRNA, but not for a polymerase‐complex formation. These findings suggest that these two aspartates are key players in the nucleotidyl transfer reaction using two metal ions.     

New DNA viruses identified in patients with acute viral infection syndrome

A sequence-independent PCR amplification method was used to identify viral nucleic acids in the plasma samples of 25 individuals presenting with symptoms of acute viral infection following high-risk behavior for human immunodeficiency virus type 1 transmission. GB virus C/hepatitis G virus was identified in three individuals and hepatitis B virus in one individual. Three previously undescribed DNA viruses were also detected, a parvovirus and two viruses related to TT virus (TTV). Nucleic acids in human plasma that were distantly related to bacterial sequences or with no detectable similarities to known sequences were also found. Nearly complete viral genome sequencing and phylogenetic analysis confirmed the presence of a new parvovirus distinct from known human and animal parvoviruses and of two related TTV-like viruses highly divergent from both the TTV and TTV-like minivirus groups. The detection of two previously undescribed viral species in a small group of individuals presenting acute viral syndrome with unknown etiology indicates that a rich yield of new human viruses may be readily identifiable using simple methods of sequence-independent nucleic acid amplification and limited sequencing.     

Transgenic cassava resistance to <Emphasis Type="Italic">African cassava mosaic virus</Emphasis> is enhanced by viral DNA-A bidirectional promoter-derived siRNAs

Expression of double-stranded RNA (dsRNA) homologous to virus sequences can effectively interfere with RNA virus infection in plant cells by triggering RNA silencing. Here we applied this approach against a DNA virus, African cassava mosaic virus (ACMV), in its natural host cassava. Transgenic cassava plants were developed to express small interfering RNAs (siRNA) from a CaMV 35S promoter-controlled, intron-containing dsRNA cognate to the common region-containing bidirectional promoter of ACMV DNA-A. In two of three independent transgenic lines, accelerated plant recovery from ACMV-NOg infection was observed, which correlates with the presence of transgene-derived siRNAs 21–24 nt in length. Overall, cassava mosaic disease symptoms were dramatically attenuated in these two lines and less viral DNA accumulation was detected in their leaves than in those of wild-type plants. In a transient replication assay using leaf disks from the two transgenic lines, strongly reduced accumulation of viral single-stranded DNA was observed. Our study suggests that a natural RNA silencing mechanism targeting DNA viruses through production of virus-derived siRNAs is turned on earlier and more efficiently in transgenic plants expressing dsRNA cognate to the viral promoter and common region.     

Efficient RNA virus control in Drosophila requires the RNA methyltransferase Dnmt2

Drosophila use small‐interfering RNA mechanisms to limit the amplification of viral genomes. However, it is unclear how small RNA interference components recognize and separate viral from cellular RNA. Dnmt2 enzymes are highly conserved RNA methyltransferases with substrate specificity towards cellular tRNAs. We report here that Dnmt2 is required for efficient innate immune responses in Drosophila. Dnmt2 mutant flies accumulate increasing levels of Drosophila C virus and show activated innate immune responses. Binding of Dnmt2 to DCV RNA suggests that Dnmt2 contributes to virus control directly, possibly by RNA methylation. These observations demonstrate a role for Dnmt2 in antiviral defence.     

Hybridization Selection and In Vitro Translation of Autographa californica Nuclear Polyhedrosis Virus mRNA

We isolated polyadenylated RNA from the cytoplasm of cells infected with Autographa californica nuclear polyhedrosis virus late after infection (21 h postinfection). At that time intracellular protein synthesis was directed almost exclusively toward infected cell-specific proteins. The polyadenylic acid-containing RNA sequences in the cytoplasm at 21 h postinfection were radiolabeled in vitro and hybridized to A. californica nuclear polyhedrosis virus DNA restriction fragments. The polyadenylic acid-containing RNA was derived from regions representing the entire viral genome. Translation in a reticulocyte cell-free protein-synthesizing system of cytoplasmic RNA selected by hybridization to viral DNA and polyadenylic acid-containing RNA produced almost identical polypeptide patterns, suggesting that late after infection almost all of the cytoplasmic polyadenylic acid-containing RNA present in infected cells was of viral origin. Polyhedrin protein (molecular weight, 33,000) and a number of virion structural proteins were among the translation products which were identified by immunoprecipitation and by comparing molecular weights. In addition, some tentative nonstructural infected cell-specific proteins were also detected. Using the hybridization selection technique, we determined that sequences complementary to the message coding for polyhedrin were located on EcoRI fragment I of A. californica nuclear polyhedrosis virus DNA, whereas sequences coding for a putative nonstructural protein (molecular weight, 39,000) were on EcoRI fragment J.     

Selenium-dependent glutathione peroxidase modules encoded by RNA viruses

Glutathione peroxidase (GPx) is the prototypical eukaryotic selenoprotein, with the rare amino acid selenocysteine (Sec) at the enzyme active site, encoded by the UGA codon in RNA. A DNA virus,Molluscum contagiosum, has now been shown to encode a functional selenium-dependent GPx enzyme. Using modifications of conventional sequence database searching techniques to locate potential viral GPx modules, combined with structurally guided comparative sequence analysis, we provide compelling evidence that Se-dependent GPx modules are encoded in a number of RNA viruses, including potentially serious human pathogens like HIV-1 and hepatitis C virus, coxsackievirus B3, HIV-2, and measles virus. Analysis of the sequences of multiple viral isolates reveals conservation of the putative GPx-related features, at least within viral subtypes or genotypes, supporting the hypothesis that these are functional GPx modules.     

The expression of viral and globin genes during differentiation of the friend cell

A complementary DNA probe has been prepared from the Friend murine erythroleukaemia virus complex (FV) released from Friend cells treated with dimethylsulphoxide (DMSO). The complementary DNA (cDNA) forms a hybrid specifically with the viral RNA genome. The availability of this viral probe together with mouse globin cDNA has made it possible to study the expression of both viral and globin genes in the Friend cell during differentiation using molecular hybridisation techniques. These specific probes have been used in an attempt to determine whether any connection exists between expression of Friend virus sequences and erythroid differentiation as measured by globin gene expression. A titration technique has been used to quantitate the levels of Friend viral- and globin-specific sequences in various Friend cell lines which differ in their ability to release Friend virus in response to DMSO although all produce haemoglobin under the same conditions. The results show: (a) that Friend cell lines unable to release virus nevertheless have a large pool of entire virus specific sequences in the polysomes; (b) an increase in virus release induced with DMSO is normally associated with a modest increase in viral sequence in the polysomes; (c) most cell lines show an early accumulation of viral and a later increase in globin mRNA sequences; (d) in an exceptional virus-negative, BUdR-resistant cell clone (B8/3), the accumulation of globin mRNA takes place very rapidly but there is no concomitant increase in viral RNA during differentiation.     

Emaravirus属新病毒:中国枣树花叶伴随病毒全基因组测序

[目的]2016年以来,新疆阿克苏等地区出现了一种新的枣树病害,严重威胁当地及周边红枣产业。本研究旨在鉴定引起此次病害的病原,探究病原体的传播方式,为生物防治策略的开发提供研究基础。[方法]对发病植株进行小RNA测序以鉴定病原体;对新鉴定的病毒,通过RNAseq和反转录PCR获取病毒全序列;体外表达重组的病毒结构蛋白并制备特异性抗体,通过Western斑点杂交法在发病植株中确证病毒蛋白;收集发病区域的媒介昆虫,通过反转录PCR在昆虫体内检测病毒的基因组,鉴定可能的传毒介体。[结果]本研究鉴定一种新的欧洲山梣环斑病毒属病毒为新疆新发枣树病害可能的病原体,命名为中国枣树花叶伴随病毒(Chinese date mosaic-associated virus,CDMaV)。CDMaV是一种多分段单链RNA病毒,基因组由5条负义RNA组成;RNA1-RNA5大小分别为7160、2224、1230、1493、971 nt,每条基因组RNA的互补链包含一个开放阅读框,共编码5个蛋白,依次为依赖RNA的RNA聚合酶、包膜糖蛋白、核衣壳蛋白和两个未知功能蛋白。在枣树寄生虫枣瘿螨体内扩增到病毒序列,表明该病毒可能以枣瘿螨为介体在枣树间进行传播。[结论]本研究为新疆新发枣树病害鉴定了相关病原体CDMaV,完成CDMaV全基因组测序,并鉴定枣瘿螨为可能的传毒介体。鉴定病原体和传播介体是建立病害防治方法的必要基础。     

Molecular and phylogenetic characterization of honey bee viruses,Nosema microsporidia,protozoan parasites,and parasitic mites in China

China has the largest number of managed honey bee colonies, which produce the highest quantity of honey and royal jelly in the world; however, the presence of honey bee pathogens and parasites has never been rigorously identified in Chinese apiaries. We thus conducted a molecular survey of honey bee RNA viruses, Nosema microsporidia, protozoan parasites, and tracheal mites associated with nonnative Apis mellifera ligustica and native Apis cerana cerana colonies in China. We found the presence of black queen cell virus (BQCV), chronic bee paralysis virus (CBPV), deformed wing virus (DWV), Israeli acute paralysis virus (IAPV), and sacbrood virus (SBV), but not that of acute bee paralysis virus (ABPV) or Kashmir bee virus (KBV). DWV was the most prevalent in the tested samples. Phylogenies of Chinese viral isolates demonstrated that genetically heterogeneous populations of BQCV, CBPV, DWV, and A. cerana‐infecting SBV, and relatively homogenous populations of IAPV and A. meliifera‐infecting new strain of SBV with single origins, are spread in Chinese apiaries. Similar to previous observations in many countries, Nosema ceranae, but not Nosema apis, was prevalent in the tested samples. Crithidia mellificae, but not Apicystis bombi was found in five samples, including one A. c. cerana colony, demonstrating that C. mellificae is capable of infecting multiple honey bee species. Based on kinetoplast‐encoded cytochrome b sequences, the C. mellificae isolate from A. c. cerana represents a novel haplotype with 19 nucleotide differences from the Chinese and Japanese isolates from A. m. ligustica. This suggests that A. c. cerana is the native host for this specific haplotype. The tracheal mite, Acarapis woodi, was detected in one A. m. ligustica colony. Our results demonstrate that honey bee RNA viruses, N. ceranae, C. mellificae, and tracheal mites are present in Chinese apiaries, and some might be originated from native Asian honey bees.