Incipient diploidization of the medicinal plant Perilla within 10,000 years.

Perilla is a young allotetraploid Lamiaceae species widely used in East Asia as herb and oil plant. Here, we report the high-quality, chromosome-scale genomes of the tetraploid (Perilla frutescens) and the AA diploid progenitor (Perilla citriodora). Comparative analyses suggest post Neolithic allotetraploidization within 10,000 years, and nucleotide mutation in tetraploid is 10% more than in diploid, both of which are dominated by G:C → A:T transitions. Incipient diploidization is characterized by balanced swaps of homeologous segments, and subsequent homeologous exchanges are enriched towards telomeres, with excess of replacements of AA genes by fractionated BB homeologs. Population analyses suggest that the crispa lines are close to the nascent tetraploid, and involvement of acyl-CoA: lysophosphatidylcholine acyltransferase gene for high α-linolenic acid content of seed oil is revealed by GWAS. These resources and findings provide insights into incipient diploidization and basis for breeding improvement of this medicinal plant.

Genetic diversity and population structure of Alternaria species from tomato and potato in North Carolina and Wisconsin.

Early blight (EB) caused by Alternaria linariae or Alternaria solani and leaf blight (LB) caused by A. alternata are economically important diseases of tomato and potato. Little is known about the genetic diversity and population structure of these pathogens in the United States. A total of 214 isolates of A. alternata (n = 61), A. linariae (n = 96), and A. solani (n = 57) were collected from tomato and potato in North Carolina and Wisconsin and grouped into populations based on geographic locations and tomato varieties. We exploited 220 single nucleotide polymorphisms derived from DNA sequences of 10 microsatellite loci to analyse the population genetic structure between species and between populations within species and infer the mode of reproduction. High genetic variation and genotypic diversity were observed in all the populations analysed. The null hypothesis of the clonality test based on the index of association [Formula: see text] was rejected, and equal frequencies of mating types under random mating were detected in some studied populations of Alternaria spp., suggesting that recombination can play an important role in the evolution of these pathogens. Most genetic differences were found between species, and the results showed three distinct genetic clusters corresponding to the three Alternaria spp. We found no evidence for clustering of geographic location populations or tomato variety populations. Analyses of molecular variance revealed high (> 85%) genetic variation within individuals in a population, confirming a lack of population subdivision within species. Alternaria linariae populations harboured more multilocus genotypes (MLGs) than A. alternata and A. solani populations and shared the same MLG between populations within a species, which was suggestive of gene flow and population expansion. Although both A. linariae and A. solani can cause EB on tomatoes and potatoes, these two species are genetically differentiated. Our results provide new insights into the evolution and structure of Alternaria spp. and can lead to new directions in optimizing management strategies to mitigate the impact of these pathogens on tomato and potato production in North Carolina and Wisconsin.

A single nucleotide substitution in the coding region of Ogura male sterile gene, orf138, determines effectiveness of a fertility restorer gene, Rfo, in radish.

Cytoplasmic male sterility (CMS) observed in many plants leads defect in the production of functional pollen, while the expression of CMS is suppressed by a fertility restorer gene in the nuclear genome. Ogura CMS of radish is induced by a mitochondrial orf138, and a fertility restorer gene, Rfo, encodes a P-type PPR protein, ORF687, acting at the translational level. But, the exact function of ORF687 is still unclear. We found a Japanese variety showing male sterility even in the presence of Rfo. We examined the pollen fertility, Rfo expression, and orf138 mRNA in progenies of this variety. The progeny with Type H orf138 and Rfo showed male sterility when their orf138 mRNA was unprocessed within the coding region. By contrast, all progeny with Type A orf138 were fertile though orf138 mRNA remained unprocessed in the coding region, demonstrating that ORF687 functions on Type A but not on Type H. In silico analysis suggested a specific binding site of ORF687 in the coding region, not the 5' untranslated region estimated previously, of Type A. A single nucleotide substitution in the putative binding site diminishes affinity of ORF687 in Type H and is most likely the cause of the ineffectiveness of ORF687. Furthermore, fertility restoration by RNA processing at a novel site in some progeny plants indicated a new and the third fertility restorer gene, Rfs, for orf138. This study clarified that direct ORF687 binding to the coding region of orf138 is essential for fertility restoration by Rfo.

Identification of miRNA-mRNA Regulatory Modules Involved in Lipid Metabolism and Seed Development in a Woody Oil Tree (Camellia oleifera).

Tea oil camellia (Camellia oleifera), an important woody oil tree, is a source of seed oil of high nutritional and medicinal value that is widely planted in southern China. However, there is no report on the identification of the miRNAs involved in lipid metabolism and seed development in the high- and low-oil cultivars of tea oil camellia. Thus, we explored the roles of miRNAs in the key periods of oil formation and accumulation in the seeds of tea oil camellia and identified miRNA-mRNA regulatory modules involved in lipid metabolism and seed development. Sixteen small RNA libraries for four development stages of seed oil biosynthesis in high- and low-oil cultivars were constructed. A total of 196 miRNAs, including 156 known miRNAs from 35 families, and 40 novel miRNAs were identified, and 55 significantly differentially expressed miRNAs were found, which included 34 upregulated miRNAs, and 21 downregulated miRNAs. An integrated analysis of the miRNA and mRNA transcriptome sequence data revealed that 10 miRNA-mRNA regulatory modules were related to lipid metabolism; for example, the regulatory modules of ath-miR858b-MYB82/MYB3/MYB44 repressed seed oil biosynthesis, and a regulation module of csi-miR166e-5p-S-ACP-DES6 was involved in the formation and accumulation of oleic acid. A total of 23 miRNA-mRNA regulatory modules were involved in the regulation of the seed size, such as the regulatory module of hpe-miR162a_L-2-ARF19, involved in early seed development. A total of 12 miRNA-mRNA regulatory modules regulating growth and development were identified, such as the regulatory modules of han-miR156a_L+1-SPL4/SBP2, promoting early seed development. The expression changes of six miRNAs and their target genes were validated using quantitative real-time PCR, and the targeting relationship of the cpa-miR393_R-1-AFB2 regulatory module was verified by luciferase assays. These data provide important theoretical values and a scientific basis for the genetic improvement of new cultivars of tea oil camellia in the future.

Reproductive phasiRNAs regulate reprogramming of gene expression and meiotic progression in rice.

Plant spermatogenesis is a complex process that directly affects crop breeding. A rapid change in gene abundance occurs at early meiosis prophase, when gene regulation is selective. However, how these genes are regulated remains unknown. Here, we show that rice reproductive phasiRNAs are essential for the elimination of a specific set of RNAs during meiotic prophase I. These phasiRNAs cleave target mRNAs in a regulatory manner such that one phasiRNA can target more than one gene, and/or a single gene can be targeted by more than one phasiRNA to efficiently silence target genes. Our investigation of phasiRNA-knockdown and PHAS-edited transgenic plants demonstrates that phasiRNAs and their nucleotide variations are required for meiosis progression and fertility. This study highlights the importance of reproductive phasiRNAs for the reprogramming of gene expression during meiotic progression and establishes a basis for future studies on the roles of phasiRNAs with a goal of crop improvement.

Identification and Molecular Mechanisms of Key Nucleotides Causing Attenuation in Pathogenicity of Dahlia Isolate of Potato Spindle Tuber Viroid.

While the potato spindle tuber viroid (PSTVd) variant, PSTVd-Dahlia (PSTVd-D or PSTVd-Dwt) induces very mild symptoms in tomato cultivar 'Rutgers', PSTVd-Intermediate (PSTVd-I or PSTVd-Iwt) induces severe symptoms. These two variants differ by nine nucleotides, of which six mutations are located in the terminal left (TL) to the pathogenicity (P) domains. To evaluate the importance of mutations located in the TL to the P domains, ten types of point mutants were created by swapping the nucleotides between the two viroid variants. Bioassay in tomato plants demonstrated that two mutants created on PSTVd-Iwt at positions 42 and 64 resulted in symptom attenuation. Phenotypic and RT-qPCR analysis revealed that mutation at position 42 of PSTVd-Iwt significantly reduced disease severity and accumulation of the viroid, whereas mutation at position 64 showed a significant reduction in stunting when compared to the PSTVd-Iwt infected plant. RT-qPCR analysis on pathogenesis-related protein 1b1 and chalcone synthase genes showed a direct correlation with symptom severity whereas the expansin genes were down-regulated irrespective of the symptom severity. These results indicate that the nucleotides at positions 42 and 64 are in concert with the ones at positions 43, 310, and 311/312, which determines the slower and stable accumulation of PSTVd-D without eliciting excessive host defense responses thus contributing in the attenuation of disease symptom.

Luteolin escape mutants of dengue virus map to prM and NS2B and reveal viral plasticity during maturation.

We previously showed that luteolin, a well-known plant-derived component found in the "heat clearing" class of Traditional Chinese Medicine (TCM) herbs, is an uncompetitive inhibitor (Ki 58.6 µM) of the host proprotein convertase furin, an endoprotease that is required for maturation of flaviviruses in the trans-Golgi compartment. Luteolin also weakly inhibited recombinant dengue virus NS2B/NS3 protease (Ki 140.36 µM) non-competitively. In order to further explore the mechanism of inhibition we isolated resistant mutants by continuous passaging of DENV2 in the presence of increasing concentrations of luteolin. Nucleotide sequence analysis of the luteolin-resistant escape mutants revealed nucleotide changes that lead to amino acid substitutions in the prM (T79R) and NS2B (I114M) genes. These mutations were introduced into a DENV2 infectious clone and tested for replication in Huh-7 cells. Interestingly we found that the replication kinetics of prM T19R-NS2B I114M double-mutant (DM) was similar to wild-type virus (WT). On the other hand the prM T79R single mutant (SM1) was attenuated and the NS2B I114M single mutant (SM2) showed enhanced replication. Time of drug addition assay with luteolin showed that the mutant viruses were able to produce more mature virions than WT in the order DM > SM2>SM1>WT. Exogenous addition of furin to purified immature WT or mutant viruses revealed that luteolin blocked the prM cleavage of WT and SM2 at a similar level. On the other hand the SM1 immature virus showed some cleavage while the DM immature virus revealed efficient furin cleavage of prM even in the presence of 50 µM luteolin. Our findings suggest that luteolin inhibition of furin may occur at host/pathogen interface that permits the virus to escape the suppression by mutating key residue that may lead to an altered interface.

An origin of the immunogenicity of in vitro transcribed RNA.

The emergence of RNA-based therapeutics demands robust and economical methods to produce RNA with few byproducts from aberrant activity. While in vitro transcription using the bacteriophage T7 RNA polymerase is one such popular method, its transcripts are known to display an immune-stimulatory activity that is often undesirable and uncontrollable. We here showed that the immune-stimulatory activity of T7 transcript is contributed by its aberrant activity to initiate transcription from a promoter-less DNA end. This activity results in the production of an antisense RNA that is fully complementary to the intended sense RNA product, and consequently a long double-stranded RNA (dsRNA) that can robustly stimulate a cytosolic pattern recognition receptor, MDA5. This promoter-independent transcriptional activity of the T7 RNA polymerase was observed for a wide range of DNA sequences and lengths, but can be suppressed by altering the transcription reaction with modified nucleotides or by reducing the Mg2+ concentration. The current work thus not only offers a previously unappreciated mechanism by which T7 transcripts stimulate the innate immune system, but also shows that the immune-stimulatory activity can be readily regulated.

Complete mitochondrial genome and evolutionary analysis of Turritopsis dohrnii, the "immortal" jellyfish with a reversible life-cycle.

Turritopsis dohrnii (Cnidaria, Hydrozoa, Hydroidolina, Anthoathecata) is the only known metazoan that is capable of reversing its life cycle via morph rejuvenation from the adult medusa stage to the juvenile polyp stage. Here, we present a complete mitochondrial (mt) genome sequence of T. dohrnii, which harbors genes for 13 proteins, two transfer RNAs, and two ribosomal RNAs. The T. dohrnii mt genome is characterized by typical features of species in the Hydroidolina subclass, such as a high A+T content (71.5%), reversed transcriptional orientation for the large rRNA subunit gene, and paucity of CGN codons. An incomplete complementary duplicate of the cox1 gene was found at the 5' end of the T. dohrnii mt chromosome, as were variable repeat regions flanking the chromosome. We identified species-specific variations (nad5, nad6, cob, and cox1 genes) and putative selective constraints (atp8, nad1, nad2, and nad5 genes) in the mt genes of T. dohrnii, and predicted alterations in tertiary structures of respiratory chain proteins (NADH4, NADH5, and COX1 proteins) of T. dohrnii. Based on comparative analyses of available hydrozoan mt genomes, we also determined the taxonomic relationships of T. dohrnii, recovering Filifera IV as a paraphyletic taxon, and assessed intraspecific diversity of various Hydrozoa species.

Evidence of a Strong Domestication Bottleneck in the Recently Cultivated New Zealand Endemic Root Crop, Arthropodium cirratum (Asparagaceae).

We use chloroplast DNA sequencing to examine aspects of the pre-European Maori cultivation of an endemic New Zealand root crop, Arthropodium cirratum (rengarenga). Researching the early stages of domestication is not possible for the majority of crops, because their cultivation began many thousands of years ago and/or they have been substantially altered by modern breeding methods. We found high levels of genetic variation and structuring characterised the natural distribution of A. cirratum, while the translocated populations only retained low levels of this diversity, indicating a strong bottleneck even at the early stages of this species' cultivation. The high structuring detected at four chloroplast loci within the natural A. cirratum range enabled the putative source(s) of the translocated populations to be identified as most likely located in the eastern Bay of Plenty/East Cape region. The high structuring within A. cirratum also has implications for the conservation of genetic diversity within this species, which has undergone recent declines in both its natural and translocated ranges.

Maximal dinucleotide and trinucleotide circular codes.

We determine here the number and the list of maximal dinucleotide and trinucleotide circular codes. We prove that there is no maximal dinucleotide circular code having strictly less than 6 elements (maximum size of dinucleotide circular codes). On the other hand, a computer calculus shows that there are maximal trinucleotide circular codes with less than 20 elements (maximum size of trinucleotide circular codes). More precisely, there are maximal trinucleotide circular codes with 14, 15, 16, 17, 18 and 19 elements and no maximal trinucleotide circular code having less than 14 elements. We give the same information for the maximal self-complementary dinucleotide and trinucleotide circular codes. The amino acid distribution of maximal trinucleotide circular codes is also determined.

Maximal dinucleotide comma-free codes.

The problem of retrieval and maintenance of the correct reading frame plays a significant role in RNA transcription. Circular codes, and especially comma-free codes, can help to understand the underlying mechanisms of error-detection in this process. In recent years much attention has been paid to the investigation of trinucleotide circular codes (see, for instance, Fimmel et al., 2014; Fimmel and Strüngmann, 2015a; Michel and Pirillo, 2012; Michel et al., 2012, 2008), while dinucleotide codes had been touched on only marginally, even though dinucleotides are associated to important biological functions. Recently, all maximal dinucleotide circular codes were classified (Fimmel et al., 2015; Michel and Pirillo, 2013). The present paper studies maximal dinucleotide comma-free codes and their close connection to maximal dinucleotide circular codes. We give a construction principle for such codes and provide a graphical representation that allows them to be visualized geometrically. Moreover, we compare the results for dinucleotide codes with the corresponding situation for trinucleotide maximal self-complementary C(3)-codes. Finally, the results obtained are discussed with respect to Crick׳s hypothesis about frame-shift-detecting codes without commas.

Loss of function of OsMADS3 via the insertion of a novel retrotransposon leads to recessive male sterility in rice (Oryza sativa).

Natural mutation is the source of natural variation, which is the fundamental basis for the genetic improvement of crops. During the process of developing a recombinant inbred line (RI), a spontaneous mutagenesis in RI127 led to the production of the recessive male-sterile line RI127S. Via a map-based cloning approach, the gene controlling the male sterility was identified as OsMADS3, which was previously reported to be associated with floral organ development and male sterility. Thermal asymmetric interlaced PCR isolated one 1633-bp insertion in OsMADS3 in RI127S, which damaged its function due to failed transcription. The 1633-bp insertion was derived from a fragment flanked by retrotransposon genes on chromosome 5. Seven haplotypes of OsMADS3 were observed among 529 cultivars and 107 wild rice accessions, and 98% of the investigated genotypes carried the same H2 haplotype, indicating that OsMADS3 is highly conserved. RI127S has the combined genome constitution of its parents, indica rice Teqing and japonica 02428, and carries the widely compatible S5 gene donated by 02428. RI127 exhibits good performance in regard to its agronomic traits and has a wide compatibility. Therefore, RI127S would be an elite mediator for recurrent breeding in cases requiring a tedious hand-crossing-based inter-crossing phase. RI127S can be crossed not only with indica rice but also with japonica rice, thus providing breeders with flexible arrangements in recurrent breeding programs.

Metabolic rescue in pluripotent cells from patients with mtDNA disease.

Mitochondria have a major role in energy production via oxidative phosphorylation, which is dependent on the expression of critical genes encoded by mitochondrial (mt)DNA. Mutations in mtDNA can cause fatal or severely debilitating disorders with limited treatment options. Clinical manifestations vary based on mutation type and heteroplasmy (that is, the relative levels of mutant and wild-type mtDNA within each cell). Here we generated genetically corrected pluripotent stem cells (PSCs) from patients with mtDNA disease. Multiple induced pluripotent stem (iPS) cell lines were derived from patients with common heteroplasmic mutations including 3243A>G, causing mitochondrial encephalomyopathy and stroke-like episodes (MELAS), and 8993T>G and 13513G>A, implicated in Leigh syndrome. Isogenic MELAS and Leigh syndrome iPS cell lines were generated containing exclusively wild-type or mutant mtDNA through spontaneous segregation of heteroplasmic mtDNA in proliferating fibroblasts. Furthermore, somatic cell nuclear transfer (SCNT) enabled replacement of mutant mtDNA from homoplasmic 8993T>G fibroblasts to generate corrected Leigh-NT1 PSCs. Although Leigh-NT1 PSCs contained donor oocyte wild-type mtDNA (human haplotype D4a) that differed from Leigh syndrome patient haplotype (F1a) at a total of 47 nucleotide sites, Leigh-NT1 cells displayed transcriptomic profiles similar to those in embryo-derived PSCs carrying wild-type mtDNA, indicative of normal nuclear-to-mitochondrial interactions. Moreover, genetically rescued patient PSCs displayed normal metabolic function compared to impaired oxygen consumption and ATP production observed in mutant cells. We conclude that both reprogramming approaches offer complementary strategies for derivation of PSCs containing exclusively wild-type mtDNA, through spontaneous segregation of heteroplasmic mtDNA in individual iPS cell lines or mitochondrial replacement by SCNT in homoplasmic mtDNA-based disease.

Cloning and polymorphism analysis of the 2-methyl-6-phytyl-1,4-benzoquinol methyltransferase gene (VTE3) in Arachis hypogaea, A. duranensis, and A. ipaënsis.

One of the critical enzymes involved in vitamin E biosynthesis in plants is 2-methyl-6-phytyl-1,4-benzoquinol methyltransferase (MPBQ MT). The full-length VTE3 cDNA (designated rVTE3-1 and -2) encoding MPBQ MT and the full-length DNA of VTE3 (designated gVTE3-1 and -2) were isolated from cultivated peanuts (Arachis hypogaea). The full-length DNA of VTE3 (designated gVTE3-A and -B) was isolated from the wild groundnut species A. duranensis (A-genome) and A. ipaënsis (B-genome), and polymorphism analysis of VTE3 was performed. The results demonstrated that rVTE3-1 and -2 both have a DNA sequence that is 1059 bp long and encodes 351 amino acids; the homology of the 2 amino acid sequences was 98.6%. The gVTE3-1 of cultivated Fenghua 2 peanut samples was 2710 bp long, with 3 introns located at 44-163, 772-1295, and 1603-2437 bp, and the Fenghua 2 gVTE3-2 was 2706 bp long, with 3 introns located at 44-169, 778-1291, and 1599-2433 bp. The homology for gVTE3-1 and -2 across 13 cultivated peanut samples was 99.9 and 100%, respectively. gVTE3-1 and -2 were from A- and B-genome, respectively, with 96.6% homology between the 2 sequences. The present study demonstrated that abundant polymorphisms were present in the VTE3 genes from different genomes. Additionally, polymorphisms were observed in the gVTE3-1 alleles of the 13 cultivars and wild A. duranensis species but not in the gVTE3-2 alleles of the 13 cultivars and wild A. ipaënsis species.

Effect of a novel mutation in a Δ9-stearoyl-ACP-desaturase on soybean seed oil composition.

Soybean [Glycine max (L.) Merr.] oil typically contains 2-4% stearic acid. Seed oil with 20% stearic acid would be useful for solid fat applications, both for its cooking properties and health benefits. Breeding lines with high stearic acid have been developed, but many suffer from agronomic problems. This study identifies a new source of high stearic acid, determines its relationship with another high stearic locus and presents molecular markers for it is use in breeding. TCJWB03-806-7-19, a 'Holladay' mutant with high stearic acid, was crossed to two FAM94-41-derived lines that contained a point mutation in a seed-specific isoform of a Δ9-stearoyl-acyl carrier protein-desaturase (SACPD-C). Fatty acid analysis was performed over two growing seasons with F(2)-derived lines and transgressive segregation for stearic acid content was observed. Sequencing of SACPD isoforms in TCJWB03-806-7-19 revealed the deletion of an 'A' nucleotide in exon 3 of SACPD-B, which results in a protein whose final 28 amino acids are predicted to differ from Williams 82 SACPD-B. Sorting intolerant from tolerant (SIFT) analysis revealed that this frameshift mutation may affect SACPD-B protein function. Allele-specific genotyping for the SACPD-C point mutation and SACPD-B nucleotide deletion was performed in both populations. Additive effects and R(2) for stearic acid were +3.3 and 0.55 for SACPD-C and +1.9 and 0.19 for SACPD-B. Average stearic acid in lines homozygous for both mutations was 14.6%. This SACPD-B mutation represents a novel high stearic allele.

[Geno- and phenotypic characteristic of Bacillus strains--components of endosporin].

Endosporin is used in veterinary for the prophylaxis and treatment of disbacteriosis, intestinal infections, festering wounds and postpartum pyoinflammatory complications in agricultural animals. The probiotic is based on two Bacillus strains which inhibit growth of a broad spectrum of pathogenic microorganisms and synthesise proteolytic enzymes and other biologically active secondary metabolites, particularly - polysaccharides. The activity of these two strains was supplementary. For the species identification of these strains, sequences of 16S rRNA genes and fatty acid content of cell walls were analysed. It was found that the both strains belong to B. velezensis. Limitations of application of 16S rRNA sequences for identification of closely related species are discussed in the paper. A method of 16S rRNA sequence profiling by polymorphic nucleotides was proposed. It was also shown that usefulness of Bacillus strains in probiotics is mostly based on their unique strain specific properties rather than on general species characteristics.

SCAR markers for discriminating species of two genera of medicinal plants, Liriope and Ophiopogon.

The development of DNA markers that can closely discriminate between Liriope and Ophiopogon species is vital for efficient and accurate identification of these species, and to ensure the quality, safety, and efficacy of medicines made from these plants. We developed species-specific molecular markers for these two genera. Forty RAPD primers were tested to detect polymorphism; species-specific RAPD bands were gel-purified, cloned, and sequenced. Primers for sequence-characterized amplified regions (SCARs) were then designed, based on nucleotide sequences of specific RAPD primers. SCAR markers SA06 and SB05, specific to Ophiopogon japonicus, amplified 460- and 553-bp DNA fragments, respectively. The marker SA12 amplified a 485-bp fragment specific to Liriope platyphylla. This is the first report of a species-specific SCAR marker for this group. These markers will be useful for rapid identification of closely related Liriope and Ophiopogon species.

[Study on optimization of SRAP-PCR reaction system for Pinellia ternata in Suzhou].

OBJECTIVE: To investigate the optimization system of SRAP-PCR molecular marker technology in the analysis on Pinellia ternata. METHOD: SRAP-PCR reaction system for P. ternata was optimized by L16 (5(4)) orthogonal design with five elements (dNTPs, Mg2+, the template DNA, primers, Taq enzyme) and four standards. RESULT: The most suitable forward primer for SRAP for Pinellia ternata was 5'-TGAGTCCAAACCGGAAG-3', while the reverse primer was 5'-GACTGCGTACGAATTACG-3'. The optimized reaction system contained 70 ng DNA template, 0.9 micromol x L(-1) primer, 0.20 mmol x L(-1) dNTP s, 1.5 - 2.0 mmol x L(-1) Mg2+, and 2 U Taq enzyme. CONCLUSION: SRAP-PCR system for P. ternata is established to lay a foundation for future construction of SRAP genetic map of P. ternata.

[Spontaneous nucleotide substitution in potato virus X genes].

Elucidation of the factors underlying nucleotide substitution types and rates in genomes is of significance for improving the prevalent models for molecular evolution, and is essential in studying the evolution process, distinguishing the coding and noncoding RNAs, identifying biologically significant motifs such as replication signals, control regions and so on. However, many questions concerning the factors that affect the nucleotide substitutions in plant virus genomes still remain unclear. In this paper a comparative computer analysis of nucleotide substitution in genomic sequences of 10 pairs of closely related potato virus X (PVX) strains having equal genome and gene length has been carried out. A significant frequency variation of nucleotide substitutions, depending on nucleotide and strain pairs, length and localization of genome and gene regions as well as substitution types and their codon positions, was found. The highest relations of substitution frequency are found in different nucleotide pairs (u-->c/c-->g, 1240), codon positions (third/second, 17.3) and substitution types (transitions/transversions, 5.8), the least ones - in forward and back substitutions in the same nucleotide pairs (1.03-1.6). The transversion frequencies are generally significantly lower but more variable than the transition one, and there is no direct relationship between gene (region) length and transversion frequency.