Analysis of genetic variation and diversity of Rice stripe virus populations through high-throughput sequencing

Genetic variability of rice stripe virus after its pandemic in Jiangsu

BACKGROUND

Rice stripe virus (RSV) caused a serious disease pandemic in rice in East China between 2001 and 2010. The continuous integrated managements reduced virus epidemic year by year until it was non-epidemic. As an RNA virus, its genetic variability after undergoing a long-term non-epidemic period was meaningful to study. While in 2019, the sudden occurrence of RSV in Jiangsu provided an opportunity for the study.

METHODS AND RESULTS

The complete genome of JY2019, an RSV isolate from Jiangyan, was determined. A genotype profile of 22 isolates from China, Japan and Korea indicated that the isolates from Yunnan formed the subtype II, and other isolates clustered the subtype I. RNA 1-3 of JY2019 isolate well-clustered in the subtype I clade, and RNA 4 was also in subtype I, but it had a slight separation from other intra-group isolates. After phylogenetic analyses, it was considered NSvc4 gene contributed to the tendency, because it exhibited an obvious trend towards the subtype II (Yunnan) group. High sequence identity (100%) of NSvc4 between JY2019 and barnyardgrass isolate from different regions demonstrated genetic variation of NSvc4 was consistent in RSV natural populations in Jiangsu in the non-epidemic period. In the phylogenetic tree of all 74 NSvc4 genes, JY2019 belonged to a minor subtype Ib, suggesting the subtype Ib isolates might have existed in natural populations before the non-epidemic period, but not a dominant population.

CONCLUSIONS

Our results suggested that NSvc4 gene was susceptible to selection pressure, and the subtype Ib might be more adaptable for the interaction between RSV and hosts in the non-epidemic ecological conditions.

Comparison of High Throughput Sequencing to Standard Protocols for Virus Detection in Berry Crops

We completed a comprehensive study comparing virus detection between high throughput sequencing (HTS) and standard protocols in 30 berry selections (12 Fragaria, 10 Vaccinium, and eight Rubus) with known virus profiles. The study examined temporal detection of viruses at four sampling times encompassing two growing seasons. Within the standard protocols, reverse transcription (RT) PCR proved better than biological indexing. Detection of known viruses by HTS and RT-PCR nearly mirrored each other. HTS provided superior detection compared with RT-PCR on a wide spectrum of variants and discovery of novel viruses. More importantly, in most cases in which the two protocols showed parallel virus detection, 11 viruses in 16 selections were not consistently detected by both methods at all sampling points. Based on these data, we propose a testing requirement of four sampling times over two growing seasons for berry and potentially other crops, to ensure that no virus remains undetected independent of titer, distribution, or other virus-virus or virus-host interactions.

Rice stripe virus: Exploring Molecular Weapons in the Arsenal of a Negative-Sense RNA Virus

Rice stripe disease caused by Rice stripe virus (RSV) is one of the most devastating plant viruses of rice and causes enormous losses in production. RSV is transmitted from plant to plant by the small brown planthopper (Laodelphax striatellus) in a circulative-propagative manner. The recent reemergence of this pathogen in East Asia since 2000 has made RSV one of the most studied plant viruses over the past two decades. Extensive studies of RSV have resulted in substantial advances regarding fundamental aspects of the virus infection. Here, we compile and analyze recent information on RSV with a special emphasis on the strategies that RSV has adopted to establish infections. These advances include RSV replication and movement in host plants and the small brown planthopper vector, innate immunity defenses against RSV infection, epidemiology, and recent advances in the management of rice stripe disease. Understanding these issues will facilitate the design of novel antiviral therapies for management and contribute to a more detailed understanding of negative-sense virus-host interactions at the molecular level.

The population structure of Rose rosette virus in the USA

Rose rosette virus (RRV) (genus Emaravirus) is the causal agent of the homonymous disease, the most destructive malady of roses in the USA. Although the importance of the disease is recognized, little sequence information and no full genomes are available for RRV, a multi-segmented RNA virus. To better understand the population structure of the virus we implemented a Hi-Plex PCR amplicon high-throughput sequencing approach to sequence all 7 segments and to quantify polymorphisms in 91 RRV isolates collected from 16 states in the USA. Analysis revealed insertion/deletion (indel) polymorphisms primarily in the 5' and 3' non-coding, but also within coding regions, including some resulting in changes of protein length. Phylogenetic analysis showed little geographical structuring, suggesting that topography does not have a strong influence on virus evolution. Overall, the virus populations were homogeneous, possibly because of regular movement of plants, the recent emergence of RRV and/or because the virus is under strong purification selection to preserve its integrity and biological functions.

High Throughput Sequencing For Plant Virus Detection and Discovery

Over the last decade, virologists have discovered an unprecedented number of viruses using high throughput sequencing (HTS), which led to the advancement of our knowledge on the diversity of viruses in nature, particularly unraveling the virome of many agricultural crops. However, these new virus discoveries have often widened the gaps in our understanding of virus biology; the forefront of which is the actual role of a new virus in disease, if any. Yet, when used critically in etiological studies, HTS is a powerful tool to establish disease causality between the virus and its host. Conversely, with globalization, movement of plant material is increasingly more common and often a point of dispute between countries. HTS could potentially resolve these issues given its capacity to detect and discover. Although many pipelines are available for plant virus discovery, all share a common backbone. A description of the process of plant virus detection and discovery from HTS data are presented, providing a summary of the different pipelines available for scientists' utility in their research.

Immune responses induced by different genotypes of the disease-specific protein of Rice stripe virus in the vector insect

Persistent plant viruses circulate between host plants and vector insects, possibly leading to the genetic divergence in viral populations. We analyzed the single nucleotide polymorphisms (SNPs) of Rice stripe virus (RSV) when it incubated in the small brown planthopper and rice. Two SNPs, which lead to nonsynonymous substitutions in the disease-specific protein (SP) of RSV, produced three genotypes, i.e., GG, AA and GA. The GG type mainly existed in the early infection period of RSV in the planthoppers and was gradually substituted by the other two genotypes during viral transmission. The two SNPs did not affect the interactions of SP with rice PsbP or with RSV coat protein. The GG genotype of SP induced stronger immune responses than those of the other two genotypes in the pattern recognition molecule and immune-responsive effector pathways. These findings demonstrated the population variations of RSV during the circulation between the vector insect and host plant.

Genome-Wide Variation in Potyviruses

Potyviruses (family Potyviridae, genus Potyvirus) are the result of an initial radiation event that occurred 6,600 years ago. The genus currently consists of 167 species that infect monocots or dicots, including domesticated and wild plants. Potyviruses are transmitted in a non-persistent way by more than 200 species of aphids. As indicated by their wide host range, worldwide distribution, and diversity of their vectors, potyviruses have an outstanding capacity to adapt to new hosts and environments. However, factors that confer adaptability are poorly understood. Viral RNA-dependent RNA polymerases introduce nucleotide substitutions that generate genetic diversity. We hypothesized that selection imposed by hosts and vectors creates a footprint in areas of the genome involved in host adaptation. Here, we profiled genomic and polyprotein variation in all species in the genus Potyvirus. Results showed that the potyviral genome is under strong negative selection. Accordingly, the genome and polyprotein sequence are remarkably stable. However, nucleotide and amino acid substitutions across the potyviral genome are not randomly distributed and are not determined by codon usage. Instead, substitutions preferentially accumulate in hypervariable areas at homologous locations across potyviruses. At a frequency that is higher than that of the rest of the genome, hypervariable areas accumulate non-synonymous nucleotide substitutions and sites under positive selection. Our results show, for the first time, that there is correlation between host range and the frequency of sites under positive selection. Hypervariable areas map to the N terminal part of protein P1, N and C terminal parts of helper component proteinase (HC-Pro), the C terminal part of protein P3, VPg, the C terminal part of NIb (RNA-dependent RNA polymerase), and the N terminal part of the coat protein (CP). Additionally, a hypervariable area at the NIb-CP junction showed that there is variability in the sequence of the NIa protease cleavage sites. Structural alignment showed that the hypervariable area in the CP maps to the N terminal flexible loop and includes the motif required for aphid transmission. Collectively, results described here show that potyviruses contain fixed hypervariable areas in key parts of the genome which provide mutational robustness and are potentially involved in host adaptation.

A Multi-Level Strategy Based on Metabolic and Molecular Genetic Approaches for the Characterization of Different Coptis Medicines Using HPLC-UV and RAD-seq Techniques

Coptis plants (Ranunculaceae) to have played an important role in the prevention and treatment human diseases in Chinese history. In this study, a multi-level strategy based on metabolic and molecular genetic methods was performed for the characterization of four Coptis herbs (C. chinensis, C. deltoidea, C. omeiensis and C. teeta) using high performance liquid chromatography-ultraviolet (HPLC-UV) and restriction site-associated DNA sequencing (RAD-seq) techniques. Protoberberine alkaloids including berberine, palmatine, coptisine, epiberberine, columbamine, jatrorrhizine, magnoflorine and groenlandicine in rhizomes were identified and determined based on the HPLC-UV method. Among them, berberine was demonstrated as the most abundant compound in these plants. RAD-seq was applied to discover single nucleotide polymorphisms (SNPs) data. A total of 44,747,016 reads were generated and 2,443,407 SNPs were identified in regarding to four plants. Additionally, with respect to complicated metabolic and SNP data, multivariable statistical methods of principal component analysis (PCA) and hierarchical cluster analysis (HCA) were successively applied to interpret the structure characteristics. The metabolic variation and genetic relationship among different Coptis plants were successfully illustrated based on data visualization. Summarily, this comprehensive strategy has been proven as a reliable and effective approach to characterize Coptis plants, which can provide additional information for their quality assessment.

Whole genome deep sequencing revealed host impact on population structure, variation and evolution of Rice stripe virus

High-throughput deep sequencing and variant detection showed that variations of Rice stripe virus (RSV) populations obtained from small brown planthopper-transmitted rice plants and sap-inoculated N. benthamiana plants were single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels). The SNPs were more uniform across RSV genome, but InDels occurred mainly in the intergenic regions (IRs) and in the 5' or 3' noncoding regions. There were no clear patterns of InDels, although the inserted sequences were all from virus itself. Six, one, and one non-synonymous substitutions were respectively observed in the RdRP ORF, IR and the movement protein ORF. These non-synonymous substitutions were found to be stable, resulting in new consensus sequences in the NBL11 RSV population. Furthermore, the numbers of SNPs and InDels in RSV genome from N. benthamiana plants were much higher than that from O. sativa plants. These differences are likely caused by selection pressures generated by different host plants.

SNP discovery of Korean short day onion inbred lines using double digest restriction site-associated DNA sequencing

Onion (Allium cepa L.) is an economically important vegetable crop around the world. Genetic and genomic research into various onion accessions will provide insights into the onion genome to enhance breeding strategies and improve crops. However, the onion’s large genome size means that studies of molecular markers are limited in onion. This study aimed to discover high quality single nucleotide polymorphisms (SNPs) from 192 onion inbred lines relating to short-day cultivation in Korea. Paired-end (PE) double digested restriction site-associated DNA sequencing (ddRAD-seq) was used to discover SNPs in onion. A total of 538,973,706 reads (25.9 GB), with an average of 2,658,491 high-quality reads, were generated using ddRAD-seq. With stringent filtering, 1904 SNPs were discovered based on onion reference scaffolds. Further, population structure and genetic relationship studies suggested that two well-differentiated sub-populations exist in onion lines. SNP-associated flanking sequences were also compared with a public non-redundant database for gene ontology and pathway analysis. To our knowledge, this is the first report to identify high-quality SNPs in onion based on reference sequences using the ddRAD-seq platform. The SNP markers identified will be useful for breeders and the research community to deepen their understanding, enhance breeding programs, and support the management of onion genomic resources.

Genomic variations in the 3'-termini of Rice stripe virus in the rotation between vector insect and host plant

A large number of plant RNA viruses circulate between plants and insects. For RNA viruses, host alternations may impose a differential selective pressure on viral populations and induce variations in viral genomes. Here, we report the variations in the 3'-terminal regions of the multiple-segment RNA virus Rice stripe virus (RSV) that were discovered through de novo assembly of the genome using RNA sequencing data from infected host plants and vector insects. The newly assembled RSV genome contained 16- and 15-nt extensions at the 3'-termini of two genome segments compared with the published reference RSV genome. Our study demonstrated that these extensional sequences were consistently observed in two RSV isolates belonging to distinct genetic subtypes in RSV-infected rice, wheat and tobacco. Moreover, the de novo assembled genome of Southern rice black-streaked dwarf virus also contained 3'-terminal extensions in five RNA segments compared with the reference genome. Time course experiments confirmed that the 3'-terminal extensions of RSV were enriched in the vector insects, were gradually eliminated in the host plant and potentially affected viral replication. These findings indicate that variations in the 3'-termini of viral genomes may be different adaptive strategies for plant RNA viruses in insects and plants.

Deep sequencing for discovery and evolutionary analysis of plant viruses

The advent of next generation sequencing (NGS), or deep sequencing, has allowed great advances to be made in discovery, diagnostics, and evolutionary studies in plant viruses. Various methods have been used for enrichment for virus-specific nucleic acids, each of which have some drawbacks. Many novel viruses have been discovered in plants by NGS technologies, and there is a good deal of promise for more comprehensive studies in virus evolution. However, each aspect of using NGS has its caveats that need to be considered, and there is still a need for better tools of analysis, as well as method for validation of sequence variation.

Viruliferous rate of small brown planthopper is a good indicator of rice stripe disease epidemics

Rice stripe virus (RSV), its vector insect (small brown planthopper, SBPH) and climatic conditions in Jiangsu, China were monitored between 2002 and 2012 to determine key biotic and abiotic factors driving epidemics of the disease. Average disease severity, disease incidence and viruliferous rate of SBPH peaked in 2004 and then gradually decreased. Disease severity of RSV was positively correlated with viruliferous rate of the vector but not with the population density of the insect, suggesting that the proportion of vectors infected by the virus rather than the absolute number of vectors plays an important role in RSV epidemics and could be used for disease forecasting. The finding of a positive correlation of disease severity and viruliferous rate among years suggests that local infection is likely the main source of primary inoculum of RSV. Of the two main climatic factors, temperature plays a more important role than rainfall in RSV epidemics.