Global-scale computational analysis of genomic sequences reveals the recombination pattern and coevolution dynamics of cereal-infecting geminiviruses

Molecular evolution and phylogeographic analysis of wheat dwarf virus

Wheat dwarf virus (WDV) has caused considerable economic loss in the global production of grain crops. Knowledge of the evolutionary biology and population history of the pathogen remain poorly understood. We performed molecular evolution and worldwide phylodynamic analyses of the virus based on the genes in the protein-coding region of the entire viral genome. Our results showed that host-driven and geography-driven adaptation are major factors that affects the evolution of WDV. Bayesian phylogenetic analysis estimates that the average WDV substitution rate was 4.240 × 10-4 substitutions/site/year (95% credibility interval, 2.828 × 10-4-5.723 × 10-4), and the evolutionary rates of genes encoding proteins with virion-sense transcripts and genes encoding proteins with complementary-sense transcripts were different. The positively selected sites were detected in only two genes encoding proteins with complementary-sense, and WDV-barley are subject to stronger purifying selection than WDV-wheat. The time since the most recent common WDV ancestor was 1746 (95% credibility interval, 1517-1893) CE. Further analyses identified that the WDV-barley population and WDV-wheat population experienced dramatic expansion-decline episodes, and the expansion time of the WDV-barley population was earlier than that of the WDV-wheat population. Our phylogeographic analysis showed that the WDV population originating in Iran was subsequently introduced to Europe, and then spread from Eastern Europe to China.

The Past, Present, and Future of Wheat Dwarf Virus Management-A Review

Wheat dwarf disease (WDD) is an important disease of monocotyledonous species, including economically important cereals. The causative pathogen, wheat dwarf virus (WDV), is persistently transmitted mainly by the leafhopper Psammotettix alienus and can lead to high yield losses. Due to climate change, the periods of vector activity increased, and the vectors have spread to new habitats, leading to an increased importance of WDV in large parts of Europe. In the light of integrated pest management, cultivation practices and the use of resistant/tolerant host plants are currently the only effective methods to control WDV. However, knowledge of the pathosystem and epidemiology of WDD is limited, and the few known sources of genetic tolerance indicate that further research is needed. Considering the economic importance of WDD and its likely increasing relevance in the coming decades, this study provides a comprehensive compilation of knowledge on the most important aspects with information on the causal virus, its vector, symptoms, host range, and control strategies. In addition, the current status of genetic and breeding efforts to control and manage this disease in wheat will be discussed, as this is crucial to effectively manage the disease under changing environmental conditions and minimize impending yield losses.

Comparative Genome Analysis of Old World and New World TYLCV Reveals a Biasness toward Highly Variable Amino Acids in Coat Protein

Begomoviruses, belonging to the family Geminiviridae and the genus Begomovirus, are DNA viruses that are transmitted by whitefly Bemisia tabaci (Gennadius) in a circulative persistent manner. They can easily adapt to new hosts and environments due to their wide host range and global distribution. However, the factors responsible for their adaptability and coevolutionary forces are yet to be explored. Among BGVs, TYLCV exhibits the broadest range of hosts. In this study, we have identified variable and coevolving amino acid sites in the proteins of Tomato yellow leaf curl virus (TYLCV) isolates from Old World (African, Indian, Japanese, and Oceania) and New World (Central and Southern America). We focused on mutations in the coat protein (CP), as it is highly variable and interacts with both vectors and host plants. Our observations indicate that some mutations were accumulating in Old World TYLCV isolates due to positive selection, with the S149N mutation being of particular interest. This mutation is associated with TYLCV isolates that have spread in Europe and Asia and is dominant in 78% of TYLCV isolates. On the other hand, the S149T mutation is restricted to isolates from Saudi Arabia. We further explored the implications of these amino acid changes through structural modeling. The results presented in this study suggest that certain hypervariable regions in the genome of TYLCV are conserved and may be important for adapting to different host environments. These regions could contribute to the mutational robustness of the virus, allowing it to persist in different host populations.

Development of Highly Efficient Resistance to Beet Curly Top Iran Virus (Becurtovirus) in Sugar Beet (B. vulgaris) via CRISPR/Cas9 System

Beet Curly Top Iran Virus (BCTIV, Becurtovirus) is a dominant and widespread pathogen responsible for great damage and yield reduction in sugar beet production in the Mediterranean and Middle East. CRISPR-based gene editing is a versatile tool that has been successfully used in plants to improve resistance against many viral pathogens. In this study, the efficiency of gRNA/Cas9 constructs targeting the expressed genes of BCTIV was assessed in sugar beet leaves by their transient expression. Almost all positive control sugar beets revealed systemic infection and severe disease symptoms (90%), with a great biomass reduction (68%) after BCTIV agroinoculation. On the other hand, sugar beets co-agronioculated with BCTIV and gRNA/Cas9 indicated much lower systemic infection (10–55%), disease symptoms and biomass reduction (13–45%). Viral inactivation was also verified by RCA and qPCR assays for gRNA/Cas9 treated sugar beets. PCR-RE digestion and sequencing assays confirmed the gRNA/Cas9-mediated INDEL mutations at the target sites of the BCTIV genome and represented high efficiencies (53–88%), especially for those targeting BCTIV’s movement gene and its overlapping region between capsid and ssDNA regulator genes. A multiplex CRISPR approach was also tested. The most effective four gRNAs targeting all the genes of BCTIV were cloned into a Cas9-containing vector and agroinoculated into virus-infected sugar beet leaves. The results of this multiplex CRISPR system revealed almost complete viral resistance with inhibition of systemic infection and mutant escape. This is the first report of CRSIPR-mediated broad-spectrum resistance against Becurtovirus in sugar beet.

Twenty years of evolution and diversification of digitaria streak virus in Digitaria setigera

Within the family Geminiviridae, the emergence of new species results from their high mutation and recombination rates. In this study, we report the variability and evolution of digitaria streak virus (DSV), a mastrevirus isolated in 1986 from the grass Digitaria setigera in an island of the Vanuatu archipelago. Viral DNA of DSV samples was amplified from D. setigera specimens, derived from the naturally infected original plant, which were propagated in different laboratories in France and Italy for more than 20 years. From the consensus sequences, the nucleotide substitution rate was estimated for the period between a sample and the original sequence published in 1987, as well as for the period between samples. In addition, the intra-host genetic complexity and diversity of 8 DSV populations with a total of 165 sequenced haplotypes was characterized. The evolutionary rate of DSV was estimated to be between 1.13 × 10⁻⁴ and 9.87 × 10⁻⁴ substitutions/site/year, within the ranges observed in other single-stranded DNA viruses and RNA viruses. Bioinformatic analyses revealed high variability and heterogeneity in DSV populations, which confirmed that mutant spectra are continuously generated and are organized as quasispecies. The analysis of polymorphisms revealed nucleotide substitution biases in viral genomes towards deamination and oxidation of single-stranded DNA. The differences in variability in each of the genomic regions reflected a dynamic and modular evolution in the mutant spectra that was not reflected in the consensus sequences. Strikingly, the most variable region of the DSV genome, encoding the movement protein, showed rapid fixation of the mutations in the consensus sequence and a concomitant dN/dS ratio of 6.130, which suggests strong positive selection in this region. Phylogenetic analyses revealed a possible divergence in three genetic lineages from the original Vanuatu DSV isolate.

Development of loop-mediated isothermal amplification assay for rapid detection of genetically different wheat dwarf virus isolates

Wheat dwarf virus (WDV) is considered as one of the most common viruses on cereal crops. Recently, severe outbreaks of WDV have been observed especially on winter wheat in southwestern part of Poland. Moreover, the presence of genetically different WDV-barley-specific and WDV-wheat-specific forms (WDV-B and WDV-W, respectively) was confirmed. In this study, a loop-mediated isothermal amplification assay (LAMP) was developed for the first time for efficient and rapid detection of WDV-B and WDV-W in infected plants. The reaction was performed using a set of three primer pairs: WDVF3/WDVB3, WDVFIB/WDVBIP and WDVLoopF/WDVLoopB specific for coat protein coding sequence. The amplified products were analyzed by direct staining of DNA, gel electrophoresis and real-time monitoring of the amplification curves. The sensitivity of optimized reaction was tenfold higher in comparison with conventional PCR. LAMP assay developed here is a useful and practical method for the rapid detection of different WDV isolates and can be implemented by phytosanitary services.

Monoclonal Antibody-Based Serological Detection Methods for Wheat Dwarf Virus

Wheat dwarf disease caused by wheat dwarf virus (WDV) is currently present in wheat growing regions in China and causes serious losses in wheat yield. To develop reliable and effective serological detection methods for WDV, the coat protein (CP) gene of WDV was cloned and expressed in Escherichia coli. The purified recombinant CP protein was immunized to BALB/c mice, and four hybridoma cell lines (i.e. 18G10, 9G4, 23F4 and 22A10) secreting anti-WDV monoclonal antibodies (MAbs) were obtained through the hybridoma technique. Using the prepared MAbs, an antigen-coated-plate enzyme-linked immunosorbent assay (ACP-ELISA) and a dot-ELISA were established for detecting WDV in wheat samples. The most sensitive ACP-ELISA based on MAb 23F4 or 22A10 was able to detect WDV in 1:163,840 (w/v, g/mL) diluted WDV-infected wheat plant crude extracts. The dot-ELISA based on MAb 23F4 was the most sensitive and able to detect the virus in 1:5,120 (w/v, g/mL) diluted wheat plant crude extracts. A total of 128 wheat samples were collected from wheat growing regions in the Shaanxi and Qinghai provinces, China, and were screened for the presence of WDV using two developed serological assays. Results from the survey showed that approximately 62% of the samples were infected with WDV. PCR followed by DNA sequencing and sequence alignment validated the results from the two serological assays. Therefore, we consider that these two serological detection methods can be significantly useful for the control of WDV in China.

Molecular Characterization and Distribution of Two Strains of Dasheen mosaic virus on Taro in Hawaii

Dasheen mosaic virus (DsMV) is one of the major viruses affecting taro (Colocasia esculenta) production worldwide. Whole genome sequences were determined for two DsMV strains, Hawaii Strain I (KY242358) and Hawaii Strain II (KY242359), from taro in Hawaii. They represent the first full-length coding sequences of DsMV reported from the United States. Hawaii Strains I and II were 77 and 85% identical, respectively, with other completely sequenced DsMV isolates. Hawaii Strain I was most closely related to vanilla mosaic virus (VanMV) (KX505964.1), a strain of DsMV infecting vanilla in the southern Pacific Islands. Hawaii Strain II was most closely related to a taro DsMV isolate CTCRI-II-14 (KT026108.1) from India. Phylogenetic analysis of all available DsMV isolates based on amino acid sequences of their coat protein showed some correlation between host plant and genetic diversity. Analyses of DsMV genome sequences detected three recombinants from China and India among the six isolates with known complete genome sequences. The DsMV strain NC003537.1 from China is a recombinant of KJ786965.1 from India and Hawaii Strain II. Another DsMV strain KT026108.1 is a recombinant of Hawaii Strain II and NC003537.1 from China. The third DsMV strain KJ786965.1 from India is a recombinant of Hawaii Strain II and NC003537.1 from China. To our knowledge, this is the first report of recombination events in DsMV. Both Hawaii Strains I and II of DsMV were found widespread throughout the Hawaiian islands.

Simultaneous detection of wheat dwarf virus, northern cereal mosaic virus, barley yellow striate mosaic virus and rice black-streaked dwarf virus in wheat by multiplex RT-PCR

Wheat dwarf virus (WDV), barley yellow striate mosaic virus (BYSMV), rice black-streaked dwarf virus (RBSDV) and northern cereal mosaic virus (NCMV) are four viruses infecting wheat and causing similar symptoms. In this paper, a multiplex reverse transcription polymerase chain reaction (m-RT-PCR) method has been developed for the simultaneous detection and discrimination of these viruses. The protocol uses specific primer set for each virus and produces four distinct fragments (273, 565, 783 and 1296bp), detecting the presence of RBSDV, BYSMV, WDV and NCMV, respectively. Annealing temperature, concentrations of dNTP, Taq polymerase and Mg²⁺ were optimized for the m-RT-PCR. The detection limit of the assay was up to 10^-2 dilution. The amplification specificity of these primers was tested against a range of field samples from different regions of China, where RBSDV, BYSMV, WDV have been detected. This study fulfills the need for a rapid and specific wheat virus detection that also has the potential for investigating the epidemiology of these new viral diseases.

The genomic and biological characterization of Citrullus lanatus cryptic virus infecting watermelon in China

A dsRNA virus was detected in the watermelon (Citrullus lanatus) samples collected from Kaifeng, Henan province, China through the use of next generation sequencing of small RNAs. The complete genome of this virus is comprised of dsRNA-1 (1603nt) and dsRNA-2 (1466nt), both of which are single open reading frames and potentially encode a 54.2kDa RNA-dependent RNA polymerase (RdRp) and a 45.9kDa coat protein (CP), respectively. The RdRp and CP share the highest amino acid identities 85.3% and 75.4% with a previously reported Israeli strain Citrullus lanatus cryptic virus (CiLCV), respectively. Genome comparisons indicate that this virus is the same species with CiLCV, whereas the reported sequences of the Israeli strain of CiLCV are partial, and our newly identified sequences can represent the complete genome of CiLCV. Futhermore, phylogenetic tree analyses based on the RdRp sequences suggest that CiLCV is one member in the genus Deltapartitivirus, family Partitiviridae. In addition, field investigation and seed-borne bioassays show that CiLCV commonly occurs in many varieties and is transmitted though seeds at a very high rate.

Phylogenetic analysis of Wheat dwarf virus isolates from Iran

Wheat dwarf virus (WDV) adversely affects cereal production in Asia, Europe, and North Africa. In this study, sequences of several WDV isolates from Iran which is located in the Fertile Crescent were analyzed. Analysis revealed a new geographic cluster for WDV-Wheat from Iran. Recombination analysis demonstrated the existence of several breakpoints in different regions of the viral genome. Data analysis demonstrated that WDV-Barley has an older history and lower diversity than WDV-Wheat. Sequence analysis identified a rare occasion of a co-infection of wheat with WDV-Wheat and WDV-Barley.