Selection pressure and founder effects constrain genetic variation in differentiated populations of soilborne bymovirus Wheat yellow mosaic virus (Potyviridae) in China

Wheat adaptation to environmental stresses under climate change: Molecular basis and genetic improvement

Wheat (Triticum aestivum) is a staple food for about 40% of the world's population. As the global population has grown and living standards improved, high yield and improved nutritional quality have become the main targets for wheat breeding. However, wheat production has been compromised by global warming through the more frequent occurrence of extreme temperature events, which have increased water scarcity, aggravated soil salinization, caused plants to be more vulnerable to diseases, and directly reduced plant fertility and suppressed yield. One promising option to address these challenges is the genetic improvement of wheat for enhanced resistance to environmental stress. Several decades of progress in genomics and genetic engineering has tremendously advanced our understanding of the molecular and genetic mechanisms underlying abiotic and biotic stress responses in wheat. These advances have heralded what might be considered a "golden age" of functional genomics for the genetic improvement of wheat. Here, we summarize the current knowledge on the molecular and genetic basis of wheat resistance to abiotic and biotic stresses, including the QTLs/genes involved, their functional and regulatory mechanisms, and strategies for genetic modification of wheat for improved stress resistance. In addition, we also provide perspectives on some key challenges that need to be addressed.

Resistance of QYm.nau-2D to wheat yellow mosaic virus was derived from an alien introgression into common wheat

The QYm.nau-2D locus conferring wheat yellow mosaic virus resistance is an exotic introgression and we developed 11 diagnostic markers tightly linked to QYm.nau-2D. Wheat yellow mosaic virus (WYMV) is a serious disease of winter wheat in China. Breeding resistant varieties is the most effective strategy for WYMV control. A WYMV resistant locus QYm.nau-2D on the chromosome arm 2DL has been repeatedly reported but the mapped region is large. In the present study, we screened recombinants using a biparental population and mapped QYm.nau-2D into an 18.8 Mb physical interval. By genome-wide association studies of 372 wheat varieties for WYMV resistance in four environments, we narrowed down QYm.nau-2D into a 16.4 Mb interval. Haplotype analysis indicated QYm.nau-2D were present as six different states due to recombination during hybridization breeding. QYm.nau-2D was finally mapped into a linkage block of 11.2 Mb. Chromosome painting using 2D specific probes and collinearity analysis among the published sequences corresponding to QYm.nau-2D region indicated the block was an exotic introgression. The Illumina-sequenced reads of four diploid Aegilops species were mapped to the sequence of Fielder, a variety having the introgression. The mapping reads were significantly increased at the putative introgression regions of Fielder. Ae. uniaristata (NN) had the highest mapping reads, suggesting that QYm.nau-2D was possibly an introgression from genome N. We investigated the agronomic performances of different haplotypes and observed no linkage drag of the alien introgression for the 15 tested traits. For marker-assisted selection of QYm.nau-2D, we developed 11 diagnostic markers tightly linked to the locus. This research provided a case study of an exotic introgression, which has been utilized in wheat improvement for WYMV resistance.

Changes in soil fungal communities after onset of wheat yellow mosaic virus disease

Rhizosphere-associated microbes have important implications for plant health, but knowledge of the association between the pathological conditions of soil-borne virus-infected wheat and soil microbial communities, especially changes in fungal communities, remains limited. We investigated the succession of fungal communities from bulk soil to wheat rhizosphere soil in both infected and healthy plants using amplicon sequencing methods, and assessed their potential role in plant health. The results showed that the diversity of fungi in wheat rhizosphere and bulk soils significantly differed post wheat yellow mosaic virus disease onset. The structure differences in fungal community at the two wheat health states or two compartment niches were evident, soil physicochemical properties (i.e., NH4 +) contribute to differences in fungal community structure and alpha diversity. Comparison analysis showed Mortierellomycetes and Dothideomycetes as dominant communities in healthy wheat soils at class level. The genus Pyronemataceae and Solicoccozyma were significantly are significantly enriched in rhizosphere soil of diseased plant, the genus Cystofilobasidium, Cladosporium, Mortierella, and Stephanonectria are significantly enriched in bulk soil of healthy plant. Co-occurrence network analysis showed that the fungi in healthy wheat soil has higher mutual benefit and connectivity compared with diseased wheat. The results of this study demonstrated that the occurrence of wheat yellow mosaic virus diseases altered both fungal community diversity and composition, and that NH4 + is the most important soil physicochemical factor influencing fungal diversity and community composition.

Determinants of Virus Variation, Evolution, and Host Adaptation

Virus evolution is the change in the genetic structure of a viral population over time and results in the emergence of new viral variants, strains, and species with novel biological properties, including adaptation to new hosts. There are host, vector, environmental, and viral factors that contribute to virus evolution. To achieve or fine tune compatibility and successfully establish infection, viruses adapt to a particular host species or to a group of species. However, some viruses are better able to adapt to diverse hosts, vectors, and environments. Viruses generate genetic diversity through mutation, reassortment, and recombination. Plant viruses are exposed to genetic drift and selection pressures by host and vector factors, and random variants or those with a competitive advantage are fixed in the population and mediate the emergence of new viral strains or species with novel biological properties. This process creates a footprint in the virus genome evident as the preferential accumulation of substitutions, insertions, or deletions in areas of the genome that function as determinants of host adaptation. Here, with respect to plant viruses, we review the current understanding of the sources of variation, the effect of selection, and its role in virus evolution and host adaptation.

Banana bunchy top virus genetic diversity in Pakistan and association of diversity with recombination in its genomes

Banana Bunchy top virus (BBTV) is a multipartite circular single strand DNA virus that belongs to genus Babuvirus and family Nanoviridae. It causes significant crop losses worldwide and also in Pakistan. BBTV is present in Pakistan since 1988 however, till now only few (about twenty only) sequence of genomic components have been reported from the country. To have insights into current genetic diversity in Pakistan fifty-seven genomic components including five complete genomes (comprises of DNA-R, -U3, -S, -M, -C and -N components) were sequenced in this study. The genetic diversity analysis of populations from Pakistan showed that DNA-R is highly conserved followed by DNA-N, whereas DNA-U3 is highly diverse with the most diverse Common Region Stem-loop (CR-SL) in BBTV genome, a functional region, which previously been reported to have undergone recombination in Pakistani population. A Maximum Likelihood (ML) phylogenetic analysis of entire genomes of isolates by using sequence of all the components concatenated together with the reported genomes around the world revealed deeper insights about the origin of the disease in Pakistan. A comparison of the genetic diversity of Pakistani and entire BBTV populations around the world indicates that there exists a correlation between genetic diversity and recombination. Population genetics analysis indicated that the degree of selection pressure differs depending on the area and genomic component. A detailed analysis of recombination across various components and functional regions suggested that recombination is closely associated with the functional parts of BBTV genome showing high genetic diversity. Both genetic diversity and recombination analyses suggest that the CR-SL is a recombination hotspot in all BBTV genomes and among the six components DNA-U3 is the only recombined component that has extensively undergone inter and intragenomic recombination. Diversity analysis of recombinant regions results on average one and half fold increase and, in some cases up to four-fold increase due to recombination. These results suggest that recombination is significantly contributing to the genetic diversity of BBTV populations around the world.

Full Genome Evolutionary Studies of Wheat Streak Mosaic-Associated Viruses Using High-Throughput Sequencing

Wheat streak mosaic (WSM), a viral disease affecting cereals and grasses, causes substantial losses in crop yields. Wheat streak mosaic virus (WSMV) is the main causal agent of the complex, but mixed infections with Triticum mosaic virus (TriMV) and High plains wheat mosaic emaravirus (HPWMoV) were reported as well. Although resistant varieties are effective for the disease control, a WSMV resistance-breaking isolate and several potential resistance-breaking isolates have been reported, suggesting that viral populations are genetically diverse. Previous phylogenetic studies of WSMV were conducted by focusing only on the virus coat protein (CP) sequence, while there is no such study for either TriMV or HPWMoV. Here, we studied the genetic variation and evolutionary mechanisms of natural populations of WSM-associated viruses mainly in Kansas fields and fields in some other parts of the Great Plains using high-throughput RNA sequencing. In total, 28 historic and field samples were used for total RNA sequencing to obtain full genome sequences of WSM-associated viruses. Field survey results showed WSMV as the predominant virus followed by mixed infections of WSMV + TriMV. Phylogenetic analyses of the full genome sequences demonstrated that WSMV Kansas isolates are widely distributed in sub-clades. In contrast, phylogenetic analyses for TriMV isolates showed no significant diversity. Recombination was identified as the major evolutionary force of WSMV and TriMV variation in KS fields, and positive selection was detected in some encoding genomic regions in the genome of both viruses. Furthermore, the full genome sequence of a second Kansas HPWMoV isolate was reported. Here, we also identified previously unknown WSMV isolates in the Great Plains sharing clades and high nucleotide sequence similarities with Central Europe isolates. The findings of this study will provide more insights into the genetic structure of WSM-associated viruses and, in turn, help in improving strategies for disease management.

Molecular characteristics of segment 5, a unique fragment encoding two partially overlapping ORFs in the genome of rice black-streaked dwarf virus

Rice black-streaked dwarf virus (RBSDV), a ds-RNA virus in Fijivirus genus with family Reoviridae, which is transmitted by the small brown planthopper, is responsible for incidence of maize rough dwarf disease (MRDD) and rice black-streaked dwarf disease (RBSDD). To understand the variation and evolution of S5, a unique fragment in the genome of RBSDV which encodes two partially overlapping ORFs (ORF5-1 and ORF5-2), we analyzed 127 sequences from maize and rice exhibiting symptoms of dwarfism. The nucleotide diversity of both ORF5-1 (π = 0.039) and ORF5-2 (π = 0.027) was higher than that of the overlapping region (π = 0.011) (P < 0.05). ORF5-2 was under the greatest selection pressure based on codon bias analysis, and its activation was possibly influenced by the overlapping region. The recombinant fragments of three recombinant events (14NM23, 14BM20, and 14NM17) cross the overlapping region. Based on neighbor-joining tree analysis, the overlapping region could represent the evolutionary basis of the full-length S5, which was classified into three main groups. RBSDV populations were expanding and haplotype diversity resulted mainly from the overlapping region. The genetic differentiation of combinations (T127-B35, T127-J34, A58-B35, A58-J34, and B35-J34) reached significant or extremely significant levels. Gene flow was most frequent between subpopulations A58 and B35, with the smallest |Fst| (0.02930). We investigated interactions between 13 RBSDV proteins by two-hybrid screening assays and identified interactions between P5-1/P6, P6/P9-1, and P3/P6. We also observed self-interactive effects of P3, P6, P7-1, and P10. In short, we have proven that RBSDV populations were expanding and the overlapping region plays an important role in the genetic variation and evolution of RBSDV S5. Our results enable ongoing research into the evolutionary history of RBSDV-S5 with two partly overlapping ORFs.

Endoplasmic reticulum remodeling induced by Wheat yellow mosaic virus infection studied by transmission electron microscopy

Plant virus was a kind of organism lived depending on infecting viable host cell and propagated their posterity by replicating its hereditary nucleotide, transcripting into protein, assembling protein and nucleotide into virion (Ortín and Parra, 2006; Sanfaçon, 2005). Viral infection usually induces remodeling of host cell, especially endoplasmic reticulum (ER) for generating membrane packed viral factory. During the infection of Bymovirus, a kind of membranous body (MB) was generated in host cells, which is thought as an ER aggregate. In present study we performed a study on Wheat yellow mosaic virus (WYMV) induced MB by several transmission electron microscopy (TEM) based methods, including cytological observation, component analysis by immuno-gold labeling and structural analysis by electron tomography (ET). WYMV infection induced at least two morphologies of MB, including the lamella dominated morphology (lamella-MB) looked like sprawling cirrus, and the tubule dominated morphology (tubule-MB) looked like latticed network. MB was verified composing of ER as revealed by immuno-gold labeling by antibody against endoplasmic reticulum (ER) retention signal as well as by detailed observation of MB construction modules as double layer membrane. By immuno-gold labeling, both two MB morphologies (lamella-MB and tubule-MB) had same components in viral derived protein and membrane origination (from ER). Structural analysis by ET reconstruction revealed the organization of ER in MB. Lamella-MB was composed of cesER like structures arranged irregularly whereas tubule-MB was composed of tubER like structures arranged regularly. This study provided insights into the structural details in how Bymovirus utilizing host membrane system.

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.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Genome-wide association mapping of virulence gene in rice blast fungus Magnaporthe oryzae using a genotyping by sequencing approach

Magnaporthe oryzae is a fungal pathogen causing blast disease in many plant species. In this study, seventy three isolates of M. oryzae collected from rice (Oryza sativa) in 1996-2014 were genotyped using a genotyping-by-sequencing approach to detect genetic variation. An association study was performed to identify single nucleotide polymorphisms (SNPs) associated with virulence genes using 831 selected SNP and infection phenotypes on local and improved rice varieties. Population structure analysis revealed eight subpopulations. The division into eight groups was not related to the degree of virulence. Association mapping showed five SNPs associated with fungal virulence on chromosome 1, 2, 3, 4 and 7. The SNP on chromosome 1 was associated with virulence against RD6-Pi7 and IRBL7-M which might be linked to the previously reported AvrPi7.

Ribotypes of Polymyxa graminis in Wheat Samples Infected with Soilborne Wheat Viruses in China

Polymyxa graminis is an obligate parasite and important vector of more than 14 soilborne plant viruses that pose a significant threat to cereal crops in Europe, North America, and Asia. Different ribotypes or formae speciales of P. graminis have been recognized and these may be associated with different cereal hosts or with transmission of different viruses. Two soilborne viruses infecting winter wheat in China have been reported and well studied (Wheat yellow mosaic virus [WYMV, genus Bymovirus] and Chinese wheat mosaic virus [CWMV, genus Furovirus]) but there has been no reported characterization of P. graminis isolates associated with them. In this study, the ribosomal DNA internal transcribed spacer (ITS) regions of P. graminis were examined from 63 wheat samples with apparent virus symptoms obtained from 16 sites within six Chinese provinces. Their associations with soilborne viruses were investigated. Ribotype I (P. graminis f. sp. temperata) and ribotype II (P. graminis f. sp. tepida) were confirmed in winter wheat regions of China for the first time. All 63 wheat root samples were infected with ribotype I of P. graminis and 11 were also infected with ribotype II. There was no obvious association between the ribotypes and infection by either WYMV or CWMV (or double infection). Phylogenetic analysis of the P. graminis ITS1-5.8S-ITS2 sequences revealed that ribotype I in China belongs to previously reported subgroup Ib, whereas ribotype II belongs to IIa. There was considerable sequence variation (pairwise distances from 0.0219 to 0.0319) between Chinese ribotype I isolates of different regions and previously reported ribotype I isolate Ken5 (accession number HE860055.1).

Molecular Genetic Analysis and Evolution of Segment 7 in Rice Black-Streaked Dwarf Virus in China

Rice black-streaked dwarf virus (RBSDV) causes maize rough dwarf disease or rice black-streaked dwarf disease and can lead to severe yield losses in maize and rice. To analyse RBSDV evolution, codon usage bias and genetic structure were investigated in 111 maize and rice RBSDV isolates from eight geographic locations in 2013 and 2014. The linear dsRNA S7 is A+U rich, with overall codon usage biased toward codons ending with A (A3s, S7-1: 32.64%, S7-2: 29.95%) or U (U3s, S7-1: 44.18%, S7-2: 46.06%). Effective number of codons (Nc) values of 45.63 in S7-1 (the first open reading frame of S7) and 39.96 in S7-2 (the second open reading frame of S7) indicate low degrees of RBSDV-S7 codon usage bias, likely driven by mutational bias regardless of year, host, or geographical origin. Twelve optimal codons were detected in S7. The nucleotide diversity (π) of S7 sequences in 2013 isolates (0.0307) was significantly higher than in 2014 isolates (0.0244, P = 0.0226). The nucleotide diversity (π) of S7 sequences in isolates from Jinan (0.0391) was higher than that from the other seven locations (P < 0.01). Only one S7 recombinant was detected in Baoding. RBSDV isolates could be phylogenetically classified into two groups according to S7 sequences, and further classified into two subgroups. S7-1 and S7-2 were under negative and purifying selection, with respective Ka/Ks ratios of 0.0179 and 0.0537. These RBSDV populations were expanding (P < 0.01) as indicated by negative values for Tajima's D, Fu and Li's D, and Fu and Li's F. Genetic differentiation was detected in six RBSDV subpopulations (P < 0.05). Absolute Fst (0.0790) and Nm (65.12) between 2013 and 2014, absolute Fst (0.1720) and Nm (38.49) between maize and rice, and absolute Fst values of 0.0085-0.3069 and Nm values of 0.56-29.61 among these eight geographic locations revealed frequent gene flow between subpopulations. Gene flow between 2013 and 2014 was the most frequent.

An evolutionary analysis of the Secoviridae family of viruses

The plant-infecting Secoviridae family of viruses forms part of the Picornavirales order, an important group of non-enveloped viruses that infect vertebrates, arthropods, plants and algae. The impact of the secovirids on cultivated crops is significant, infecting a wide range of plants from grapevine to rice. The overwhelming majority are transmitted by ecdysozoan vectors such as nematodes, beetles and aphids. In this study, we have applied a variety of computational methods to examine the evolutionary traits of these viruses. Strong purifying selection pressures were calculated for the coat protein (CP) sequences of nine species, although for two species evidence of both codon specific and episodic diversifying selection were found. By using Bayesian phylogenetic reconstruction methods CP nucleotide substitution rates for four species were estimated to range from between 9.29×10⁻³ to 2.74×10⁻³ (subs/site/year), values which are comparable with the short-term estimates of other related plant- and animal-infecting virus species. From these data, we were able to construct a time-measured phylogeny of the subfamily Comovirinae that estimated divergence of ninety-four extant sequences occurred less than 1,000 years ago with present virus species diversifying between 50 and 250 years ago; a period coinciding with the intensification of agricultural practices in industrial societies. Although recombination (modularity) was limited to closely related taxa, significant and often unique similarities in the protein domains between secovirid and animal infecting picorna-like viruses, especially for the protease and coat protein, suggested a shared ancestry. We discuss our results in a wider context and find tentative evidence to indicate that some members of the Secoviridae might have their origins in insects, possibly colonizing plants in a number of founding events that have led to speciation. Such a scenario; virus infection between species of different taxonomic kingdoms, has significant implications for virus emergence.

Analysis of small RNAs derived from Chinese wheat mosaic virus

The virus-derived small interfering RNAs (vsiRNAs) of Chinese wheat mosaic virus (CWMV), a member of the genus Furovirus, were characterised from wheat plants by deep sequencing. CWMV vsiRNAs of 21-22 nt in length predominated, suggesting that there might be a conserved mechanism of DCL2 and DCL4 involvement in the biogenesis of vsiRNAs, as well as a common RNA silencing pathway in CWMV-infected wheat plants. The 5'-terminal base of vsiRNAs was biased towards A/U, suggesting that CWMV vsiRNAs might be loaded into diverse AGO-containing RISCs to disturb the gene expression of host plants. Possible targets for some of the vsiRNAs were predicted.

Durable field resistance to wheat yellow mosaic virus in transgenic wheat containing the antisense virus polymerase gene

Wheat yellow mosaic virus (WYMV) has spread rapidly and causes serious yield losses in the major wheat-growing areas in China. Because it is vectored by the fungus-like organism Polymyxa graminis that survives for long periods in soil, it is difficult to eliminate by conventional crop management or fungicides. There is also only limited resistance in commercial cultivars. In this research, fourteen independent transgenic events were obtained by co-transformation with the antisense NIb8 gene (the NIb replicase of WYMV) and a selectable gene bar. Four original transgenic lines (N12, N13, N14 and N15) and an offspring line (N12-1) showed high and durable resistance to WYMV in the field. Four resistant lines were shown to have segregated and only contain NIb8 (without bar) by PCR and herbicide resistance testing in the later generations. Line N12-1 showed broad-spectrum resistance to WYMV isolates from different sites in China. After growing in the infested soil, WYMV could not be detected by tissue printing and Western blot assays of transgenic wheat. The grain yield of transgenic wheat was about 10% greater than the wild-type susceptible control. Northern blot and small RNA deep sequencing analyses showed that there was no accumulation of small interfering RNAs targeting the NIb8 gene in transgenic wheat plants, suggesting that transgene RNA silencing, a common mechanism of virus-derived disease resistance, is not involved in the process of WYMV resistance. This durable and broad-spectrum resistance to WYMV in transgenic wheat will be useful for alleviating the damage caused by WYMV.

Nucleo-cytoplasmic shuttling of VPg encoded by Wheat yellow mosaic virus requires association with the coat protein

VPg (virus protein, genome-linked) is a multifunctional protein that plays important roles in viral multiplication in the cytoplasm. However, a number of VPgs encoded by plant viruses target the nucleus and this appears to be biologically significant. These VPgs may therefore be translocated between nuclear and cytoplasmic compartments during virus infection, but such nucleo-cytoplasmic transport has not been demonstrated. We report that VPg encoded by Wheat yellow mosaic virus (WYMV, genus Bymovirus, family Potyviridae) accumulated in both the nucleus and cytoplasm of infected cells, but localized exclusively in the nucleus when expressed alone in plants. Computational analyses predicted the presence of a nuclear localization signal (NLS) and a nuclear export signal (NES) in WYMV VPg. Mutational analyses showed that both the N-terminal and the NLS domains of VPg contribute to the efficiency of nuclear targeting. In vitro and in planta assays indicated that VPg interacts with WYMV coat protein (CP) and proteinase 1 (P1) proteins. Observation of VPg fused to a fluorescent protein and subcellular fractionation experiments showed that VPg was translocated to the cytoplasm when co-expressed with CP, but not with P1. Mutations in the NES domain or treatment with leptomycin B prevented VPg translocation to the cytoplasm when co-expressed with CP. Our results suggest that association with CP facilitates the nuclear export of VPg during WYMV infection.