Genetic structure and biology of Xylella fastidiosa strains causing disease in citrus and coffee in Brazil

Xylella fastidiosa is a vector-borne, plant-pathogenic bacterium that causes disease in citrus (citrus variegated chlorosis [CVC]) and coffee (coffee leaf scorch [CLS]) plants in Brazil. CVC and CLS occur sympatrically and share leafhopper vectors; thus, determining whether X. fastidiosa isolates can be dispersed from one crop to another and cause disease is of epidemiological importance. We sought to clarify the genetic and biological relationships between CVC- and CLS-causing X. fastidiosa isolates. We used cross-inoculation bioassays and microsatellite and multilocus sequence typing (MLST) approaches to determine the host range and genetic structure of 26 CVC and 20 CLS isolates collected from different regions in Brazil. Our results show that citrus and coffee X. fastidiosa isolates are biologically distinct. Cross-inoculation tests showed that isolates causing CVC and CLS in the field were able to colonize citrus and coffee plants, respectively, but not the other host, indicating biological isolation between the strains. The microsatellite analysis separated most X. fastidiosa populations tested on the basis of the host plant from which they were isolated. However, recombination among isolates was detected and a lack of congruency among phylogenetic trees was observed for the loci used in the MLST scheme. Altogether, our study indicates that CVC and CLS are caused by two biologically distinct strains of X. fastidiosa that have diverged but are genetically homogenized by frequent recombination.     

A multigene phylogenetic study of clonal diversity and divergence in North American strains of the plant pathogen Xylella fastidiosa

Xylella fastidiosa is a pathogen that causes leaf scorch and related diseases in over 100 plant species, including Pierce's disease in grapevines (PD), phony peach disease (PP), plum leaf scald (PLS), and leaf scorch in almond (ALS), oak (OAK), and oleander (OLS). We used a high-resolution DNA sequence approach to investigate the evolutionary relationships, geographic variation, and divergence times among the X. fastidiosa isolates causing these diseases in North America. Using a large data set of 10 coding loci and 26 isolates, the phylogeny of X. fastidiosa defined three major clades. Two of these clades correspond to the recently identified X. fastidiosa subspecies piercei (PD and some ALS isolates) and X. fastidiosa subsp. multiplex (OAK, PP, PLS, and some ALS isolates). The third clade grouped all of the OLS isolates into a genetically distinct group, named X. fastidiosa subsp. sandyi. These well-differentiated clades indicate that, historically, X. fastidiosa has been a clonal organism. Based on their synonymous-site divergence ( approximately 3%), these three clades probably originated more than 15,000 years ago, long before the introduction of the nonnative plants that characterize most infections. The sister clades of X. fastidiosa subsp. sandyi and X. fastidiosa subsp. piercei have synonymous-site evolutionary rates 2.9 times faster than X. fastidiosa subsp. multiplex, possibly due to generation time differences. Within X. fastidiosa subsp. multiplex, a low level ( approximately 0.1%) of genetic differentiation indicates the recent divergence of ALS isolates from the PP, PLS, and OAK isolates due to host plant adaptation and/or allopatry. The low level of variation within the X. fastidiosa subsp. piercei and X. fastidiosa subsp. sandyi clades, despite their antiquity, suggests strong selection, possibly driven by host plant adaptation.     

Evolution of duplicate genes in a tetraploid animal, Xenopus laevis

To understand the evolution of duplicate genes, we compared rates of nucleotide substitution between 17 pairs of nonallelic duplicated genes in the tetraploid frog Xenopus laevis with rates between the orthologous loci of human and rodent. For all duplicated X. laevis genes, the number of synonymous substitutions per site (dS) was greater than the number of nonsynonymous substitutions per site (dN), indicating that these genes are subject to purifying selection. There was also a significant positive correlation (r = 0.915) between dN for the X. laevis genes and dN for the mammalian genes, suggesting that, at the amino acid level, the X. laevis genes and the mammalian genes are under similar constraints. Results of relative-rate tests showed nearly equal rates of nonsynonymous substitution in each copy of the X. laevis genes; apparently there are similar constraints on both copies. No correlation was found between dS for the X. laevis genes and dS for the mammalian genes. There was a significant positive correlation both between members of pairs of duplicated X. laevis genes (r = 0.951) and between human and rodent orthologues (r = 0.854) with respect to third- position G+C content but no such relationship between the X. laevis genes and either of their mammalian orthologues. The results indicate that both copies of a duplicate gene can be subject to purifying selection and thus support the hypothesis of selection against all genotypes containing a null allele at either of two duplicate loci.      

Multilocus sequence type system for the plant pathogen Xylella fastidiosa and relative contributions of recombination and point mutation to clonal diversity

Multilocus sequence typing (MLST) identifies and groups bacterial strains based on DNA sequence data from (typically) seven housekeeping genes. MLST has also been employed to estimate the relative contributions of recombination and point mutation to clonal divergence. We applied MLST to the plant pathogen Xylella fastidiosa using an initial set of sequences for 10 loci (9.3 kb) of 25 strains from five different host plants, grapevine (PD strains), oleander (OLS strains), oak (OAK strains), almond (ALS strains), and peach (PP strains). An eBURST analysis identified six clonal complexes using the grouping criterion that each member must be identical to at least one other member at 7 or more of the 10 loci. These clonal complexes corresponded to previously identified phylogenetic clades; clonal complex 1 (CC1) (all PD strains plus two ALS strains) and CC2 (OLS strains) defined the X. fastidiosa subsp. fastidiosa and X. fastidiosa subsp. sandyi clades, while CC3 (ALS strains), CC4 (OAK strains), and CC5 (PP strains) were subclades of X. fastidiosa subsp. multiplex. CC6 (ALS strains) identified an X. fastidiosa subsp. multiplex-like group characterized by a high frequency of intersubspecific recombination. Compared to the recombination rate in other bacterial species, the recombination rate in X. fastidiosa is relatively low. Recombination between different alleles was estimated to give rise to 76% of the nucleotide changes and 31% of the allelic changes observed. The housekeeping loci holC, nuoL, leuA, gltT, cysG, petC, and lacF were chosen to form the basis of a public database for typing X. fastidiosa (www.mlst.net). These loci identified the same six clonal complexes using the strain grouping criterion of identity at five or more loci with at least one other member.     

Molecular Evolution of the Genomic RNA of Apple Stem Grooving Capillovirus

The complete genome of the German isolate AC of Apple stem grooving virus (ASGV) was sequenced. It encodes two overlapping open reading frames (ORFs), similarly to previously described ASGV isolates. Two regions of high variability were detected between the ASGV isolates, variable region 1 (V1, from amino acids (aa) 532 to 570), and variable region 2 (V2, from aa 1,583 to 1,868). The phylogenetic analysis of the V1 and V2 regions suggested that the ASGV diversity was structured by host plant species rather than geographical origin. The dN/dS ratio between nonsynonymous and synonymous nucleotide substitution rates varied greatly along the ASGV genome. Most of ORF1 showed predominant negative selection except for the two regions V1 and V2. V1 showed an elevated dN and an average dS when compared to the ORF1 background but no significant positive selection was detected. The V2 region of ORF1 showed an elevated dN and a low dS when compared to the ORF1 background with an average dN/dS????3.0 indicative of positive selection. However, the V2 area includes overlapping ORFs, making the dN/dS estimate biased. Joint estimates of the selection intensity in the different ORFs by a recent method indicated that this region of ORF1 was in fact evolving close to neutrality. This was convergent with previous results showing that introduction of stop codons in this region of ORF1 did not impair plant infection. These data suggest that the elimination of a stop codon caused the overprinting of a novel coding region over the ancestral ORF.     

Mycobacterium avium subsp. paratuberculosis and M. avium subsp. avium are independently evolved pathogenic clones of a much broader group of M. avium organisms

Mycobacterium avium comprises organisms that share the same species designation despite considerable genomic and phenotypic variability. To determine the degree and nature of variability between subspecies and strains of M. avium, we used multilocus sequencing analysis, studying 56 genetically diverse strains of M. avium that included all described subspecies. In total, 8,064 bp of sequence from 10 gene loci were studied, with 205 (2.5%) representing variable positions. The majority (149/205) of these variations were found among M. avium subsp. hominissuis organisms. Recombination was also evident in this subspecies. In contrast, there was comparatively little variability and no evidence of recombination within the pathogenic subspecies, M. avium subsp. paratuberculosis, M. avium subsp. avium, and M. avium subsp. silvaticum. Phylogenetic analysis showed that M. avium subsp. avium and M. avium subsp. silvaticum strains clustered together on one branch, while a distinct branch defined M. avium subsp. paratuberculosis organisms. Despite the independent origin of these pathogenic subspecies, an analysis of their rates of nonsynonymous (dN) to synonymous (dS) substitutions showed increased dN/dS ratios for both: 0.67 for M. avium subsp. paratuberculosis and 0.50 for M. avium subsp. avium/M. avium subsp. silvaticum, while the value was 0.08 for M. avium subsp. hominissuis organisms. In conclusion, M. avium subsp. hominissuis represents a diverse group of organisms from which two pathogenic clones (M. avium subsp. paratuberculosis and M. avium subsp. avium/M. avium subsp. silvaticum) have evolved independently.     

Population genomic analysis of a bacterial plant pathogen: novel insight into the origin of Pierce's disease of grapevine in the U.S

Invasive diseases present an increasing problem worldwide; however, genomic techniques are now available to investigate the timing and geographical origin of such introductions. We employed genomic techniques to demonstrate that the bacterial pathogen causing Pierce's disease of grapevine (PD) is not native to the US as previously assumed, but descended from a single genotype introduced from Central America. PD has posed a serious threat to the US wine industry ever since its first outbreak in Anaheim, California in the 1880s and continues to inhibit grape cultivation in a large area of the country. It is caused by infection of xylem vessels by the bacterium Xylella fastidiosa subsp. fastidiosa, a genetically distinct subspecies at least 15,000 years old. We present five independent kinds of evidence that strongly support our invasion hypothesis: 1) a genome-wide lack of genetic variability in X. fastidiosa subsp. fastidiosa found in the US, consistent with a recent common ancestor; 2) evidence for historical allopatry of the North American subspecies X. fastidiosa subsp. multiplex and X. fastidiosa subsp. fastidiosa; 3) evidence that X. fastidiosa subsp. fastidiosa evolved in a more tropical climate than X. fastidiosa subsp. multiplex; 4) much greater genetic variability in the proposed source population in Central America, variation within which the US genotypes are phylogenetically nested; and 5) the circumstantial evidence of importation of known hosts (coffee plants) from Central America directly into southern California just prior to the first known outbreak of the disease. The lack of genetic variation in X. fastidiosa subsp. fastidiosa in the US suggests that preventing additional introductions is important since new genetic variation may undermine PD control measures, or may lead to infection of other crop plants through the creation of novel genotypes via inter-subspecific recombination. In general, geographically mixing of previously isolated subspecies should be avoided.     

Strain typing of Toxoplasma gondii: comparison of antigen-coding and housekeeping genes

Molecular characterization of Toxoplasma gondii isolates is central for understanding differences in disease transmission and manifestations. Only 3 subgroups (lineages) have been discerned with subtle within-lineage variation, permitting low-resolution classification of isolates. Because proteins, coding sequences, and especially antigen-coding genes have been used extensively in previous studies, we focused on sequence variation in introns of housekeeping genes, which may be more informative for phylogenetic analysis because they evolve under lower selection. We compared sequence variation in introns of 5 housekeeping genes with 2 antigen-coding genes. Introns of housekeeping genes were slightly more polymorphic than coding and noncoding regions of antigen-coding genes and only the former showed intralineage variation. Intragenic linkage disequilibrium was complete, but intergenic linkage, although highly significant, was incomplete, suggesting that genes are partially uncoupled. Six of 7 substitutions found within the region coding for the tachyzoite surface antigen, SAG2, were nonsynonymous, indicating that diversifying selection acts on this locus. Typing isolates on the basis of housekeeping and antigen-coding genes was consistent, but the phylogenetic relationships among the resulting groups was inconsistent. A cougar isolate typed as lineage II using a restriction fragment length polymorphism assay possessed multiple unique polymorphisms, suggesting that it represents a new lineage. We concluded that introns of housekeeping genes are preferred markers for phylogenetic study, and that multilocus genotyping is preferred for typing parasites, especially from feral or unstudied environments.     

Gene expression,chromosome heterogeneity and the fast-X effect in mammals

The higher rate of non-synonymous over synonymous substitutions (dN/dS) of the X chromosome compared with autosomes is often interpreted as a consequence of X hemizygosity. However, other factors, such as gene expression, are also known to vary between X and autosomes. Analysing 4800 orthologues in six mammals, we found that gene expression levels, associated with GC content, fully account for the variation in dN/dS between X and autosomes with no detectable effect of hemizygosity. We also report an extensive variance in dN/dS and gene expression between autosomes.     

Pseudomonas viridiflava, a multi host plant pathogen with significant genetic variation at the molecular level

The pectinolytic species Pseudomonas viridiflava has a wide host range among plants, causing foliar and stem necrotic lesions and basal stem and root rots. However, little is known about the molecular evolution of this species. In this study we investigated the intraspecies genetic variation of P. viridiflava amongst local (Cretan), as well as international isolates of the pathogen. The genetic and phenotypic variability were investigated by molecular fingerprinting (rep-PCR) and partial sequencing of three housekeeping genes (gyrB, rpoD and rpoB), and by biochemical and pathogenicity profiling. The biochemical tests and pathogenicity profiling did not reveal any variability among the isolates studied. However, the molecular fingerprinting patterns and housekeeping gene sequences clearly differentiated them. In a broader phylogenetic comparison of housekeeping gene sequences deposited in GenBank, significant genetic variability at the molecular level was found between isolates of P. viridiflava originated from different host species as well as among isolates from the same host. Our results provide a basis for more comprehensive understanding of the biology, sources and shifts in genetic diversity and evolution of P. viridiflava populations and should support the development of molecular identification tools and epidemiological studies in diseases caused by this species.     

Differential rates of evolution for the ZFY-related zinc finger genes,Zfy, Zfx,and Zfa in the mouse genus Mus

A comparative study of the last exon of the zinc finger genes Zfx, Zfy, and Zfa from species of mice in the genus Mus was conducted to assess the extent of gene-specific and chromosome-specific effects on the evolutionary patterns among related X-, Y-, and autosomal-linked genes. Phylogenetic analyses of 29 sequences from Zfx, Zfa, and Zfy from 10 taxa were performed to infer relatedness among the zinc finger loci, and codon-based maximum likelihood analyses were conducted to assess evolutionary pattern among genes. Five models of nucleotide sequence evolution were applied and compared using a likelihood ratio test. Estimates of nonsynonymous to synonymous changes (dN/dS) for these genes suggest that amino acid substitutions are occurring at a more rapid rate across the autosomal- and Y-specific lineages compared to the X-specific lineage, with the Y-specific lineage showing the highest rate under certain models. The data suggest the action of gene-specific effects on evolutionary pattern. In particular, Zfa and Zfy genes, both with presumed restricted expression, appear less functionally constrained relative to ubiquitously expressed Zfx. Slightly elevated dN/dS for Zfy genes in comparison to Zfa also suggest Y-specific effects.     

Preponderance of slightly deleterious polymorphism in mitochondrial DNA: nonsynonymous/synonymous rate ratio is much higher within species than between species

We estimated synonymous (dN) and nonsynonymous (dS) substitution rates for protein-coding genes of the mitochondrial genome from two individuals each of the species human, chimpanzee, and gorilla. The genes were analyzed both separately and in a combined data set. Pairwise sequence comparisons suggest that the dN/dS rate ratios are about 5-10 times higher in within-species comparisons than in between-species comparisons. This result is confirmed by a more rigorous likelihood ratio test, which rejected the null hypothesis that the dN/dS rate ratios are identical within and between species. The likelihood models account for the genetic code structure, transition/transversion rate ratio, and codon usage bias and are expected to produce more reliable results than the commonly used contingency test. Separate analyses of different genes show that the dN/dS rate ratios are higher within species than between species for all 13 mitochondrial genes, with the difference being statistically significant for all except three small or slowly evolving genes. Furthermore, in conserved genes, nonsynonymous rates within species tend to be higher than the between-species rates by a greater proportion than in fast-changing genes. Our findings confirm and extend earlier results obtained from smaller data sets and suggest the operation of slightly deleterious mutations throughout the mitochondrial genome in the hominoids. Implications of the results for evolutionary studies and, in particular, for studies of the origin of modern humans, are discussed.     

Context-dependent transmission of a generalist plant pathogen: host species and pathogen strain mediate insect vector competence

The specificity of pathogen–vector–host interactions is an important element of disease epidemiology. For generalist pathogens, different pathogen strains, vector species, or host species may all contribute to variability in disease incidence. One such pathogen is Xylella fastidiosa Wells et al., a xylem-limited bacterium that infects dozens of crop, ornamental, and native plants in the USA. This pathogen also has a diverse vector complex and multiple biologically distinct strains. We studied the implications of diversity in this pathogen–vector–host system, by quantifying variability in transmission efficiency of different X. fastidiosa strains (isolates from almond and grape genetic groups) for different host plants (grape, almond, and alfalfa) by two of the most important vectors in California: glassy-winged sharpshooter [ Homalodisca vitripennis (Germar)] and green sharpshooter ( Draeculacephala minerva Ball) (both Hemiptera: Cicadellidae). Transmission of isolates of the almond strain by H. vitripennis did not differ significantly, whereas transmission varied significantly among isolates from the grape strain (15–90%). Host plant species did not affect H. vitripennis transmission. Conversely, D. minerva efficiency was mediated by both host plant species and pathogen strain. No acquisition of an almond isolate occurred regardless of plant type (0/122), whereas acquisition of a grape isolate from alfalfa was 10-fold higher than from grape or almond plants. These results suggest that pathogen, vector, and host diversity impose contingencies on the transmission ecology of this complex disease system. Studies aimed at the development of management strategies for X. fastidiosa diseases should consider the complexity of these interactions as they relate to disease spread.     

Molecular phylogeny of the salmonellae: relationships among Salmonella species and subspecies determined from four housekeeping genes and evidence of lateral gene transfer events

The salmonellae are a diverse group of bacteria within the family Enterobacteriaceae that includes two species, Salmonella enterica and Salmonella bongori. In order to characterize the phylogenetic relationships of the species and subspecies of Salmonella, we analyzed four housekeeping genes, gapA, phoP, mdh and recA, comprising 3,459 bp of nucleotide sequence data for each isolate sequenced. Sixty-one isolates representing the most common serotypes of the seven subspecies of Salmonella enterica and six isolates of Salmonella bongori were included in this study. We present a robust phylogeny of the Salmonella species and subspecies that clearly defines the lineages comprising diphasic and monophasic subspecies. Evidence of intersubspecies lateral gene transfer of the housekeeping gene recA, which has not previously been reported, was obtained.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.     

Detecting genetic introgression: high levels of intersubspecific recombination found in Xylella fastidiosa in Brazil

Documenting the role of novel mutation versus homologous recombination in bacterial evolution, and especially in the invasion of new hosts, is central to understanding the long-term dynamics of pathogenic bacteria. We used multilocus sequence typing (MLST) to study this issue in Xylella fastidiosa subsp. pauca from Brazil, a bacterium causing citrus variegated chlorosis (CVC) and coffee leaf scorch (CLS). All 55 citrus isolates typed (plus one coffee isolate) defined three similar sequence types (STs) dominated by ST11 (85%), while the remaining 22 coffee isolates defined two STs, mainly ST16 (74%). This low level of variation masked unusually large allelic differences (>1% divergence with no intermediates) at five loci (leuA, petC, malF, cysG, and holC). We developed an introgression test to detect whether these large differences were due to introgression via homologous recombination from another X. fastidiosa subspecies. Using additional sequencing around these loci, we established that the seven randomly chosen MLST targets contained seven regions of introgression totaling 2,172 bp of 4,161 bp (52%), only 409 bp (10%) of which were detected by other recombination tests. This high level of introgression suggests the hypothesis that X. fastidiosa subsp. pauca became pathogenic on citrus and coffee (crops cultivated in Brazil for several hundred years) only recently after it gained genetic variation via intersubspecific recombination, facilitating a switch from native hosts. A candidate donor is the subspecies infecting plum in the region since 1935 (possibly X. fastidiosa subsp. multiplex). This hypothesis predicts that nonrecombinant native X. fastidiosa subsp. pauca (not yet isolated) does not cause disease in citrus or coffee.     

Evolutionary genetics of the isocitrate dehydrogenase gene (icd) in Escherichia coli and Salmonella enterica.

Sequences of the icd gene, encoding isocitrate dehydrogenase (IDH), were obtained for 33 strains representing the major phylogenetic lineages of Escherichia coli and Salmonella enterica. Evolutionary relationships of the strains based on variation in icd are generally similar to those previously obtained for several other housekeeping and for invasion genes, but the sequences of S. enterica subspecies V strains are unusual in being almost intermediate between those of the other S. enterica subspecies and E. coli. For S. enterica, the ratio of synonymous (silent) to nonsynonymous (replacement) nucleotide substitutions between pairs of strains was larger than comparable values for 12 other housekeeping and invasion genes, reflecting unusually strong purifying selection against amino acid replacement in the IDH enzyme. All amino acids involved in the catalytic activity and conformational changes of IDH are strictly conserved within and between species. In E. coli, the level of variation at the 3' end of the gene is elevated by the presence in some strains of a 165-bp replacement sequence supplied by the integration of either lambdoid phage 21 or defective prophage element e14. The 72 members of the E. coli Reference Collection (ECOR) and five additional E. coli strains were surveyed for the presence of phage 21 (as prophage) by PCR amplification of a phage 21-specific fragment in and adjacent to the host icd, and the sequence of the phage 21 segment extending from the 3' end of icd through the integrase gene (int) was determined in nine strains of E. coli. Phage 21 was found in 39% of E. coli strains, and its distribution among the ECOR strains is nonrandom. In two ECOR strains, the phage 21 int gene is interrupted by a 1,313-bp insertion element that has 99.3% nucleotide sequence identity with IS3411 of E. coli. The phylogenetic relationships of phage 21 strains derived from sequences of two different genomic regions were strongly incongruent, providing evidence of frequent recombination.     

Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

An excess of nonsynonymous substitutions over synonymous ones is an important indicator of positive selection at the molecular level. A lineage that underwent Darwinian selection may have a nonsynonymous/synonymous rate ratio (dN/dS) that is different from those of other lineages or greater than one. In this paper, several codon-based likelihood models that allow for variable dN/dS ratios among lineages were developed. They were then used to construct likelihood ratio tests to examine whether the dN/dS ratio is variable among evolutionary lineages, whether the ratio for a few lineages of interest is different from the background ratio for other lineages in the phylogeny, and whether the dN/dS ratio for the lineages of interest is greater than one. The tests were applied to the lysozyme genes of 24 primate species. The dN/dS ratios were found to differ significantly among lineages, indicating that the evolution of primate lysozymes is episodic, which is incompatible with the neutral theory. Maximum- likelihood estimates of parameters suggested that about nine nonsynonymous and zero synonymous nucleotide substitutions occurred in the lineage leading to hominoids, and the dN/dS ratio for that lineage is significantly greater than one. The corresponding estimates for the lineage ancestral to colobine monkeys were nine and one, and the dN/dS ratio for the lineage is not significantly greater than one, although it is significantly higher than the background ratio. The likelihood analysis thus confirmed most, but not all, conclusions Messier and Stewart reached using reconstructed ancestral sequences to estimate synonymous and nonsynonymous rates for different lineages.      

A Sliding Window-Based Method to Detect Selective Constraints in Protein-Coding Genes and Its Application to RNA Viruses

Here we present a new sliding window-based method specially designed to detect selective constraints in specific regions of a multiple protein-coding sequence alignment. In contrast to previous window-based procedures, our method is based on a nonarbitrary statistical approach to find the appropriate codon-window size to test deviations of synonymous (dS) and nonsynonymous (dN) nucleotide substitutions from the expectation. The probabilities of dN and dS are obtained from simulated data and used to detect significant deviations of dN and dS in a specific window region of the real sequence alignment. The nonsynonymous-to-synonymous rate ratio (w = dN/dS) was used to highlight selective constraints in any window wherein dS or dN was significantly different from the expectation. In these significant windows, w and its variance [V(w)] were calculated and used to test the neutral hypothesis. Computer simulations showed that the method is accurate even for highly divergent sequences. The main advantages of the new method are that it (i) uses a statistically appropriate window size to detect different selective patterns, (ii) is computationally less intensive than maximum likelihood methods, and (iii) detects saturation of synonymous sites, which can give deviations from neutrality. Hence, it allows the analysis of highly divergent sequences and the test of different alternative hypothesis as well. The application of the method to different human immunodeficiency virus type 1 and to foot-and-mouth disease virus genes confirms the action of positive selection on previously described regions as well as on new regions.