ggmsa: a visual exploration tool for multiple sequence alignment and associated data.

The identification of the conserved and variable regions in the multiple sequence alignment (MSA) is critical to accelerating the process of understanding the function of genes. MSA visualizations allow us to transform sequence features into understandable visual representations. As the sequence-structure-function relationship gains increasing attention in molecular biology studies, the simple display of nucleotide or protein sequence alignment is not satisfied. A more scalable visualization is required to broaden the scope of sequence investigation. Here we present ggmsa, an R package for mining comprehensive sequence features and integrating the associated data of MSA by a variety of display methods. To uncover sequence conservation patterns, variations and recombination at the site level, sequence bundles, sequence logos, stacked sequence alignment and comparative plots are implemented. ggmsa supports integrating the correlation of MSA sequences and their phenotypes, as well as other traits such as ancestral sequences, molecular structures, molecular functions and expression levels. We also design a new visualization method for genome alignments in multiple alignment format to explore the pattern of within and between species variation. Combining these visual representations with prime knowledge, ggmsa assists researchers in discovering MSA and making decisions. The ggmsa package is open-source software released under the Artistic-2.0 license, and it is freely available on Bioconductor (https://bioconductor.org/packages/ggmsa) and Github (https://github.com/YuLab-SMU/ggmsa).

Genomic analysis of the brassica pathogen turnip mosaic potyvirus reveals its spread along the former trade routes of the Silk Road.

Plant pathogens have agricultural impacts on a global scale and resolving the timing and route of their spread can aid crop protection and inform control strategies. However, the evolutionary and phylogeographic history of plant pathogens in Eurasia remains largely unknown because of the difficulties in sampling across such a large landmass. Here, we show that turnip mosaic potyvirus (TuMV), a significant pathogen of brassica crops, spread from west to east across Eurasia from about the 17th century CE. We used a Bayesian phylogenetic approach to analyze 579 whole genome sequences and up to 713 partial sequences of TuMV, including 122 previously unknown genome sequences from isolates that we collected over the past five decades. Our phylogeographic and molecular clock analyses showed that TuMV isolates of the Asian-Brassica/Raphanus (BR) and basal-BR groups and world-Brassica3 (B3) subgroup spread from the center of emergence to the rest of Eurasia in relation to the host plants grown in each country. The migration pathways of TuMV have retraced some of the major historical trade arteries in Eurasia, a network that formed the Silk Road, and the regional variation of the virus is partly characterized by different type patterns of recombinants. Our study presents a complex and detailed picture of the timescale and major transmission routes of an important plant pathogen.

Rice Stripe Tenuivirus Has a Greater Tendency To Use the Prime-and-Realign Mechanism in Transcription of Genomic than in Transcription of Antigenomic Template RNAs.

Most segmented negative-sense RNA viruses employ a process termed cap snatching, during which they snatch capped RNA leaders from host cellular mRNAs and use the snatched leaders as primers for transcription, leading to the synthesis of viral mRNAs with 5' heterogeneous sequences (HSs). With traditional methods, only a few HSs can be determined, and identification of their donors is difficult. Here, the mRNA 5' ends of Rice stripe tenuivirus (RSV) and Rice grassy stunt tenuivirus (RGSV) and those of their host rice were determined by high-throughput sequencing. Millions of tenuiviral HSs were obtained, and a large number of them mapped to the 5' ends of corresponding host cellular mRNAs. Repeats of the dinucleotide AC, which are complementary to the U1G2 of the tenuiviral template 3'-U1G2U3G4UUUCG, were found to be prevalent at the 3' termini of tenuiviral HSs. Most of these ACs did not match host cellular mRNAs, supporting the idea that tenuiviruses use the prime-and-realign mechanism during cap snatching. We previously reported a greater tendency of RSV than RGSV to use the prime-and-realign mechanism in transcription with leaders cap snatched from a coinfecting reovirus. Besides confirming this observation in natural tenuiviral infections, the data here additionally reveal that RSV has a greater tendency to use this mechanism in transcribing genomic than in transcribing antigenomic templates. The data also suggest that tenuiviruses cap snatch host cellular mRNAs from translation- and photosynthesis-related genes, and capped RNA leaders snatched by tenuiviruses base pair with U1/U3 or G2/G4 of viral templates. These results provide unprecedented insights into the cap-snatching process of tenuiviruses.IMPORTANCE Many segmented negative-sense RNA viruses (segmented NSVs) are medically or agriculturally important pathogens. The cap-snatching process is a promising target for the development of antiviral strategies against this group of viruses. However, many details of this process remain poorly characterized. Tenuiviruses constitute a genus of agriculturally important segmented NSVs, several members of which are major viral pathogens of rice. Here, we for the first time adopted a high-throughput sequencing strategy to determine the 5' heterogeneous sequences (HSs) of tenuiviruses and mapped them to host cellular mRNAs. Besides providing deep insights into the cap snatching of tenuiviruses, the data obtained provide clear evidence to support several previously proposed models regarding cap snatching. Curiously and importantly, the data here reveal that not only different tenuiviruses but also the same tenuivirus synthesizing different mRNAs use the prime-and-realign mechanism with different tendencies during their cap snatching.

Molecular characterization of a new potyvirus infecting passion fruit.

A potyvirus (isolate PFV-FJ) infecting passion fruit in China was identified by small-RNA sequencing. The complete genome sequence of PFV-FJ was determined to be 9974 nucleotides, excluding the poly(A) tail. PFV-FJ shares 70-72% nucleotide and 69-74% amino acid sequence identity at the polyprotein level with seven reported potyviruses, but 89% nucleotide and 91% amino acid sequence identity with an unreported potyvirus, tentatively named "passionfruit Vietnam potyvirus" (PVNV-DakNong). This suggests that PFV-FJ and PVNV-DakNong should belong to the same potyvirus species and that PFV-FJ is a new member of the genus Potyvirus. This new potyvirus was tentatively named "passion fruit severe mottle-associated virus".

The first complete genome sequence of narcissus latent virus from Narcissus.

The complete sequence of a narcissus virus isolated from the Netherlands (Narv-NL) was determined to be 8172 nucleotides in length with an open reading frame encoding for 2624 amino acids. Phylogenetic analysis indicates that Narv-NL is clustered with high confidence among representative members from the genus Macluravirus, including artichoke latent virus (ArLV) and Chinese yam necrotic mosaic virus (CYNMV). Sequence analyses indicated Narv-NL shares 67%-69% nucleotide and 51%-68% amino acid sequence identity with ArLV and CYNMV either in the complete ORF or the coat protein (CP) gene, whereas it had 81%-99 % nucleotide and 80%-99 % amino acid sequence identity with the corresponding CP sequences of narcissus latent virus (NLV) isolates, suggesting that Narv-NL is a member of NLV. To our knowledge, this is the first report of the complete sequence of a NLV isolate.

Metabolites Produced by an Endophytic Phomopsis sp. and Their Anti-TMV Activity.

The fermentation and isolation of metabolites produced by an endophytic fungus, which was identified as Phomopsis sp. FJBR-11, based on phylogenetic analysis, led to the identification of six compounds, including dothiorelones A-C, and H, and cytosporones C and U. Among these compounds, cytosporone U exhibited potent inhibitory activity against Tobacco mosaic virus (TMV). Moreover, the crude and a purified exopolysaccharide were proved to possess strong inhibitory effects against the virus infection.

A multigene phylogentic analysis of Potato virus Y and its application in strain identification.

The objective of this study is to develop a rapid and accurate multigene phylogenetic analysis to identify Potato virus Y (PVY) strains. The phylogenetic relationships of strains within the PVY species were evaluated with isolate-strain association using five datasets of concatenated sequences from the P1, HC-pro, VPg and CP genes to determine the best dataset for PVY strain identification. Results from phylogenetic analyses and Bayesian tip-association significance (BaTS) tests indicated that the major PVY strains could be distinguished using the P1, VPg and CP concatenated sequences datasets but not the remaining concatenated sequence datasets. Phylogenetic trees reconstructed from the concatenated sequences of P1, VPg and CP genes revealed that the ML and NJ trees had broadly similar topologies and that both were better than the maximum clade credibility tree (MCC). Additionally, the full genome of HLJ26, one isolate randomly selected for the multigene phylogenetic analysis, was clustered with high confidence among members of the PVYNTN-NW (SYR-Ⅱ) strain, which includes isolates of SYR-Ⅱ-2-8, SYR-Ⅱ-Be1 and SYR-Ⅱ-DrH. This suggests that it was a PVYNTN-NW (SYR-Ⅱ) isolate. Recombination analysis of this isolate identified four putative recombination joints in the P1, HC-pro/P3, VPg and the 5'-terminus of CP. This pattern is similar to that observed in the genomic structure of PVYNTN-NW (SYR-I), supporting the classification of this isolate as the PVYNTN-NW strain (SYR-Ⅱ). Simultaneously, two expected fragments of approximately 1 000 and 400 bp in size were also amplified from the isolate by a multiplex RT-PCR, consistent with the expected band pattern of the PVYNTN-NW (SYR-Ⅱ) strain. This further supports the utility of the multigene phylogenetic method in identifying PVY strains. We propose that the major PVY strains could be distinguished accurately using multigene phylogenetic analysis based on the concatenated sequences from the P1, VPg and CP genes.

Complete genome analysis of jasmine virus T from Jasminum sambac in China.

The genome of a potyvirus (isolate JaVT_FZ) recovered from jasmine (Jasminum sambac L.) showing yellow ringspot symptoms in Fuzhou, China, was sequenced. JaVT_FZ is closely related to seven other potyviruses with completely sequenced genomes, with which it shares 66-70 % nucleotide and 52-56 % amino acid sequence identity. However, the coat protein (CP) gene shares 82-92 % nucleotide and 90-97 % amino acid sequence identity with those of two partially sequenced potyviruses, named jasmine potyvirus T (JaVT-jasmine) and jasmine yellow mosaic potyvirus (JaYMV-India), respectively. This suggests that JaVT_FZ, JaVT-jasmine and JaYMV-India should be regarded as members of a single potyvirus species, for which the name "Jasmine virus T" has priority.

Genetic diversity and molecular evolution of arabis mosaic virus based on the CP gene sequence.

Arabis mosaic virus (ArMV) is a virus with a wide host range. In this study, the genetic diversity of ArMV and the molecular mechanisms underlying its evolution were investigated using the coat protein (CP) sequence. Of the 33 ArMV isolates studied, three were found to be recombinants. The other 30 recombination-free ArMV isolates could be separated into two major lineages with a significant F ST value (0.384) and tended to cluster according to their geographical origin. Different evolutionary constraints were detected for the two linages, pointing to a role of natural selection in the differentiation of ArMV.

Geographically driven adaptation of chilli veinal mottle virus revealed by genetic diversity analysis of the coat protein gene.

Chilli veinal mottle virus (ChiVMV) is an important plant pathogen with a wide host range. The genetic structure of ChiVMV was investigated by analyzing the coat protein (CP) genes of 87 ChiVMV isolates from seven Asian regions. Pairwise F ST values between ChiVMV populations ranged from 0.108 to 0.681, indicating a significant spatial structure for this pathogen. In phylogeny-trait association analysis, the viral isolates from the same region tended to group together, showing a distinct geographic feature. These results suggest that geographic driven adaptation may be an important determinant of the genetic diversity of ChiVMV.

Molecular characterization and detection of a recombinant isolate of bamboo mosaic virus from China.

The complete genome sequences of three isolates of bamboo mosaic virus (BaMV) from mainland China were determined and compared to those of BaMV isolates from Taiwan. Sequence analysis showed that isolate BaMV-JXYBZ1 from Fuzhou shares 98 % nucleotide sequence identity with BaMV-YTHSL14 from nucleotides 2586 to 6306, and more than 94 % nucleotide sequence identity with BaMV-MUZHUBZ2 in other regions. Recombination and phylogenetic analyses indicate that BaMV-JXYBZ1 is a recombinant with one recombination breakpoint. To our knowledge, this is the first report of a BaMV recombinant worldwide.

Anti-TMV activity of malformin A1, a cyclic penta-peptide produced by an endophytic fungus Aspergillus tubingensis FJBJ11.

Plant-associated microorganisms are known to produce a variety of metabolites with novel structures and interesting biological activities. An endophytic fungus FJBJ11, isolated from the plant tissue of Brucea javanica (L.) Merr. (Simaroubaceae), was proven to be significantly effective in producing metabolites with anti-Tobacco mosaic virus (TMV) activities. The isolate was identified as Aspergillus tubingensis FJBJ11 based on morphological characteristics and ITS sequence. Bioassay-guided isolation led to the identification of a cycli penta-peptide, malformin A1, along with two cyclic dipeptides, cyclo (Gly-L-Pro) and cyclo (Ala-Leu). Malformin A1 showed potent inhibitory effect against the infection and replication of TMV with IC50 values of 19.7 and 45.4 µg·mL⁻¹, as tested using local lesion assay and leaf-disc method, respectively. The results indicated the potential use of malformin A1 as a leading compound or a promising candidate of new viricide.

[Comparative analysis of population genetic structure of Potato virus Y from different hosts].

Nucleotide sequences of P3 and pipo genes of Potato virus Y (PVY) from potato and tobacco were compared to investigate the effect of hosts on the population genetic structure. Meanwhile, mutation, natural selection and gene flow were evaluated to determine evolutionary forces responsible for the population genetic dynamics. The fixation indices of population differentiation (FST) of PVY from tobacco and potato were 0.116 and 0.120, respectively with significant difference, suggesting a moderate genetic differentiation between the two populations. Genetic variation analysis showed that nucleotide identities in P3 and pipo genes among the viral isolates from tobacco were respectively in the range of 85.2%-100% and 76.5%-100% while that from potato were respectively in the range of 95.7%-100% and 93.0%-100%, indicating higher genetic variation in PVY from tobacco than that from potato. Moreover, purifying selection was detected on the majority of polymorphic sites within P3 gene, suggesting that most of mutations in the gene were harmful and consequently being eliminated by natural selection. Conversely, positive selection was detected on two polymorphic sites, suggesting that these two mutations were beneficial to PVY. Neither purifying nor positive selection was detected in pipo gene, indicating neutral evolution of the gene. The values of gene flow (Nm) between PVY populations from tobacco and potato in P3 and pipo genes were 1.91 and 1.83, respectively, suggesting strong gene flow also contributes significantly to the population genetic dynamics of PVY population. In summary, this study indicates there was a significant genetic variation in PVY hosted by tobacco and potato, and mutation, natural selection and gene flow all contribute to the genetic diversity and population dynamic of the virus.

Complete genome analysis of a novel recombinant isolate of potato virus Y from China.

The complete sequence of GF_YL20, a potato virus Y (PVY) isolate from China, encodes a polyprotein of 3,061 amino acids. Sequence analysis indicates that GF_YL20 has a genomic structure different from previously reported PVY strains. It shares 99 % nucleotide sequence identity with PB209 (PVY(N:O)) except in VPg, but more than 97 % nucleotide sequence identity with the VPg of Mont (PVY(N)), PB312 (PVY(NTN)) and HN2 (SYR-I). Phylogenetic analysis indicates that GF_YL20 is a novel N:O recombinant with three recombination breakpoints.

Sequence variation of P1 gene in Potato virus Y isolated from Fujian Province.

To understand the sequence variation and the putative protein structure of P1 gene in Potato virus Y (PVY) and to identify the sources of the variation, P1 gene in PVY isolated from Fujian Province was amplified by reverse-transcription polymerase chain reaction (RT-PCR) using a pair of degenerate primers designed from the conserved regions of published sequences. Sequence variation and putative protein structure were analyzed, and phylogenetic tree was reconstructed using Bayesian inference method. Expected fragments of 915 bp in size were amplified from 12 samples collected from Fujian Province by RT-PCR. The 12 sequences shared 73%-99% nucleotide identity with the reference sequences from GenBank. A strong recombination signal was identified at position 309 in sequences of isolates QK44, XT02, XT08 and LH05. Among the 12 sequences, 85 amino acid variants were detected, indicating high sequence variation in the P1 protein. However, positions 41-275 in the protein were highly conserved, especially in three active sites (H192, D201 and V235). Phylogenetic analysis grouped the sequences into three clades, each with different Coiled-coil domains and 3D-structures, suggesting divergent phylogenetic relationship among the groups. The above results show P1 gene in PVY is highly variable but contains 3 conserved active sites (H192, D201, V235) and the high genetic variation in the gene is primarily due to mutation and recombination.

[Sequence variation and protein structure of pipo gene in Potato virus Y].

The objectives of this study were to understand the sequence variation and the putative protein structure of pipo gene in the Potato virus Y (PVY) collected from Solanum tuberosum. The pipo gene in PVY was cloned using a pair of degenerate primers designed from its conserved region and its sequences were used to re-construct phylogenetic tree in Potyvirus genera by a Bayesian inference method. An expected fragment of 235 bp was amplified in all 20 samples by RT-PCR and the pipo genes in the 20 samples assayed shared more than 92% nucleotide sequence similarity with the published sequences of PVY strains. Among the 20 pipo gene sequences, 13 polymorphic sites were detected, including 4 parsimony informative sites and 9 singleton variable sites. These results indicate that PVY pipo gene is highly conserved but some sequence variations exist. Further analyses suggest that the pipo gene encodes a hydrophilic protein without signal peptide and transmembrane region. The protein has theoretical isoelectric points (pI) ranging from 11.26 to 11.62 and contains three highly conserved regions, especially between aa 10 and 59. The protein is likely located in the mitochondria and has a-helix secondary structure. Bayesian inference of phylogenetic trees reveals that PVY isolates are clustered in the same branch with high posterior probability, while Sunflower chlorotic mottle virus (SoCMoV) and Pepper severe mosaic virus (PepSMV) are closely related, consisting with the classification of Potyvirus genera using other approaches. Our analyses suggest that the pipo gene can be a new marker for phylogenetic analysis of the genera. The results reported in this paper provide useful insights in the genetic variation and the evolution of PVY and can stimulate further research on structure and function of the PIPO protein.

Using PhyloSuite for molecular phylogeny and tree-based analyses.

Phylogenetic analysis has entered the genomics (multilocus) era. For less experienced researchers, conquering the large number of software programs required for a multilocus-based phylogenetic reconstruction can be somewhat daunting and time-consuming. PhyloSuite, a software with a user-friendly GUI, was designed to make this process more accessible by integrating multiple software programs needed for multilocus and single-gene phylogenies and further streamlining the whole process. In this protocol, we aim to explain how to conduct each step of the phylogenetic pipeline and tree-based analyses in PhyloSuite. We also present a new version of PhyloSuite (v1.2.3), wherein we fixed some bugs, made some optimizations, and introduced some new functions, including a number of tree-based analyses, such as signal-to-noise calculation, saturation analysis, spurious species identification, and etc. The step-by-step protocol includes background information (i.e., what the step does), reasons (i.e., why do the step), and operations (i.e., how to do it). This protocol will help researchers quick-start their way through the multilocus phylogenetic analysis, especially those interested in conducting organelle-based analyses.

Complete genome sequence of narcissus late season yellows virus infecting Chinese narcissus in China.

The complete genome sequence of a Chinese narcissus isolate of narcissus late season yellows virus from Zhangzhou, China (NLSYV-ZZ), was determined to be 9,651 nucleotides in length, excluding the 3'-terminal poly (A) tail, by amplification and sequencing of virus RNA. The viral genome contains a single long open reading frame of 9,315 nucleotides encoding a polyprotein of 3,105 amino acids. The polyprotein was predicted to be cleaved into ten mature proteins by three viral proteases. Complete genome sequence comparison and phylogenetic analysis indicated that NLSYV-ZZ was most closely related to narcissus yellow stripe virus (NYSV), which was also isolated from narcissus. These viruses shared 69.9 % identity in their complete nucleotide sequences and 77.0 % identity in their polyprotein amino acid sequences.

Correction: Histone H3 gene is not a suitable marker to distinguish Alternaria tenuissima from A. alternata affecting potato.

[This corrects the article DOI: 10.1371/journal.pone.0231961.].

Bayesian Phylogeographic Inference Suggests Japan as the Center for the Origin and Dissemination of Rice Stripe Virus.

Rice stripe virus (RSV) is one of the most important viral pathogens of rice in East Asia. The origin and dispersal of RSV remain poorly understood, but an emerging hypothesis suggests that: (i) RSV originates from Yunnan, a southwest province of China; and (ii) some places of eastern China have acted as a center for the international dissemination of RSV. This hypothesis, however, has never been tested rigorously. Using a data set comprising more than 200 time-stamped coat protein gene sequences of RSV from Japan, China and South Korea, we reconstructed the phylogeographic history of RSV with Bayesian phylogeographic inference. Unexpectedly, the results did not support the abovementioned hypothesis. Instead, they suggested that RSV originates from Japan and Japan has been the major center for the dissemination of RSV in the past decades. Based on these data and the temporal dynamics of RSV reported recently by another group, we proposed a new hypothesis to explain the origin and dispersal of RSV. This new hypothesis may be valuable for further studies aiming to clarify the epidemiology of RSV. It may also be useful in designing management strategies against this devastating virus.