Joaquín Dopazo

Genetic evolution and tropism of transmissible gastroenteritis coronaviruses

Abstract Transmissible gastroenteritis virus (TGEV) is an enteropathogenic coronavirus isolated for the first time in 1946. Nonenteropathogenic porcine respiratory coronaviruses (PRCVs) have been derived from TGEV. The genetic relationship among six European PRCVs and five coronaviruses of the TGEV antigenic cluster has been determined based on their RNA sequences. The S protein of six PRCVs have an identical deletion of 224 amino acids starting at position 21. The deleted area includes the antigenic sites C and B of TGEV S glycoprotein. Interestingly, two viruses (NEB72 and TOY56) with respiratory tropism have S proteins with a size similar to the enteric viruses. NEB72 and TOY56 viruses have in the S protein 2 and 15 specific amino acid differences with the enteric viruses. Four of the residues changed (aa 219 of NEB72 isolate and as 92, 94, and 218 of TOY56) are located within the deletion present in the PRCVs and may be involved in the receptor binding site (RBS) conferring enteric tropism to TGEVs. A second RBS used by the virus to infect ST cells might be located in a conserved area between sites A and D of the S glycoprotein, since monoclonal antibodies specific for these sites inhibit the binding of the virus to ST cells. An evolutionary tree relating 13 enteric and respiratory isolates has been proposed. According to this tree, a main virus lineage evolved from a recent progenitor virus which was circulating around 1941. From this, secondary lineages originated PUR46, NEB72, TOY56, MIL65, BR170, and the PRCVs, in this order. Least squares estimation of the origin of TGEV-related coronaviruses showed a significant constancy in the fixation of mutations with time, that is, the existence of a well-defined molecular clock. A mutation fixation rate of 7 ± 2 × 10−4 nucleotide substitutions per site and per year was calculated for TGEV-related viruses. This rate falls in the range reported for other RNA viruses. Point mutations and probably recombination events have occurred during TGEV evolution.

Phylogenetic relationships of European strains of porcine reproductive and respiratory syndrome virus (PRRSV) inferred from DNA sequences of putative ORF-5 and ORF-7 genes

The complete ORF-5 gene and a fragment of the ORF-7 gene from 14 different European porcine reproductive and respiratory syndrome virus (PRRSV) isolates were amplified by reverse transcriptase polymerase chain reaction (RT-PCR), and their nucleotide sequences were determined. The ORF-7 gene displayed nucleotide and amino acid identities of 94.1-99.6% and 95.3-100% among isolates from different countries. The ORF-5 gene showed higher nucleotide (87.1-99.2% identity) and amino acid (-88% identity) variability. The resulting sequences were aligned with other European and North American PRRSV strains and phylogenetic relationships among these strains were established by the maximum parsimony method. The phylogenetic trees inferred from both genes were in agreement and showed that European and North American PRRSV strains clearly represent two different genotypes. According to both trees, there is a perfect correlation between strains and the countries in which they were isolated. Additionally, the phylogenetic position of European and North American PRRSV strains within the recently proposed family Arteriviridae was also analyzed.

Estimating errors and confidence intervals for branch lengths in phylogenetic trees by a bootstrap approach

A method, based on the bootstrap procedure, is proposed for the estimation of branch-length errors and confidence intervals in a phylogenetic tree for which equal rates of substitution among lineages do not necessarily hold. The method can be used to test whether an estimated internodal distance is significantly greater than zero. In the application of the method, any estimator of genetic distances, as well as any tree reconstruction procedure (based on distance matrices), can be used. Also the method is not limited by the number of species involved in the phylogenetic tree. An example of the application of the method in the reconstruction of the phylogenetic tree for the four hominoid species—human, chimpanzee, gorilla, and orangutan—is shown.

Annotated draft genomic sequence from a Streptococcus pneumoniae type 19F clinical isolate.

The public availability of numerous microbial genomes is enabling the analysis of bacterial biology in great detail and with an unprecedented, organism-wide and taxon-wide, broad scope. Streptococcus pneumoniae is one of the most important bacterial pathogens throughout the world. We present here sequences and functional annotations for 2.1-Mbp of pneumococcal DNA, covering more than 90% of the total estimated size of the genome. The sequenced strain is a clinical isolate resistant to macrolides and tetracycline. It carries a type 19F capsular locus, but multilocus sequence typing for several conserved genetic loci suggests that the strain sequenced belongs to a pneumococcal lineage that most often expresses a serotype 15 capsular polysaccharide. A total of 2,046 putative open reading frames (ORFs) longer than 100 amino acids were identified (average of 1,009 bp per ORF), including all described two-component systems and aminoacyl tRNA synthetases. Comparisons to other complete, or nearly complete, bacterial genomes were made and are presented in a graphical form for all the predicted proteins.

Primer design for specific diagnosis by PCR of highly variable RNA viruses: Typing of foot-and-mouth disease virus

A PCR assay for the specific detection and identification of viral sequences that correlate with established serotypes of foot-and-mouth disease virus (FMDV) has been developed. A new analysis based on homology profiles among reported sequences was used for primer design. RNA replicase (3D) gene regions that showed high homology among FMDVs, and low homology to other picornaviruses, were used for PCR amplification. Specific and highly sensitive detection was achieved for RNA of FMDV types C, A, and O, either purified or extracted from vesicular fluids of infected animals, under reaction conditions permissive for the detection of variants present in the virus population. Similarly, serotype-specific primers were designed to amplify the carboxy-terminal end of VP1 gene of FMDV types either C, A, or O. The results of PCR amplification of 15 different FMDV RNAs using type-specific primers are in agreement with the serological typing of the corresponding viruses and show that the primer-selection procedure developed for FMDV constitutes a reliable method of viral diagnosis.

A RT-PCR assay for the differential diagnosis of vesicular viral diseases of swine

A RT-PCR assay based on specific amplification of RNA sequences from each of the etiological agents of three important vesicular diseases that affect swine, foot-and-mouth disease virus (FMDV), swine vesicular disease virus (SVDV), and vesicular stomatitis virus (VSV), was developed. Genotype-specific primers that amplified DNA fragments of differential size from SVDV 3D gene or VSV L gene were selected with the aid of a computer program. Experimental testing of the primers predicted as SVDV-specific identified a primer pair, SA2/SS4, that rendered a specific product from SVDV RNAs, but did not amplify RNA from either FMDV or coxsackie B5 virus (CV-B5), a highly related picornavirus. Primers SA2/SS4 were used in combination with primers 3D2/3D1, which amplify a product of different size on FMDV 3D gene (Rodriguez et al., 1992). This combined RT-PCR reaction allowed a sensitive and specific differential detection of FMDV and SVDV RNAs in a single tube, by means of the analysis of the amplified products in agarose gels. The results obtained were similar when RNA extracted from viral stocks or plastic wells coated with either viral supernatants or extracts from lesions of infected animals, were used as starting material in the reactions. Using a similar approach, VSV serotype-specific primers IA/IS and NA/NS were selected for the specific amplification of VSV-Indiana and VSV-New Jersey RNAs, respectively. The combined use of SVDV, FMDV and VSV specific primers in a single reaction resulted in a genotype-specific amplification of each of the viral RNAs. Thus, differential diagnosis of FMDV from SVDV and/or VSV can be carried out in a single RT-PCR reaction, using a rapid and simplified methodology.

Molecular epidemiology of foot-and-mouth disease virus type O

A phylogenetic tree based on the VPI sequences of type O foot-and-mouth disease virus (FMDV) has been derived. Direct sequencing of PCR products has been used to obtain the VP1 gene sequences of new isolates. The tree exhibits four main lineages that largely correlate with the geographical origin of isolates. The analysis supports a close relationship between European O1 field isolates and vaccine strains, with the exception of O Thalheim Aus/81 and O Wuppertal Ger/82 which were probably of non-European origin. Analysis of nucleotide substitutions indicates that synonymous mutations play a major role in FMDV evolution.

Evolutionary analysis of the picornavirus family

An exhaustive evolutionary analysis of the picornavirus family has been carried out using the amino acid sequences of several proteins of the viruses including: the capsid proteins (1D, 1B, and 1C) situated at the 5′ end of the genome and responsible for the serotype of the viruses, and the viral polymerase (3D), located at the 3′ end of the genome. The evolutionary relationships found among the viruses studied support the new classification, recently suggested, in contrast to the classical one, and the existence of a new genus for the picornavirus family. In the new taxonomic organization, five genera form the picornavirus family: (1) aphthoviruses, (2) cardioviruses, (3) hepatoviruses (previously classified as enteroviruses), (4) renteroviruses (which mainly constitute a combination of the previous genera rhinovirus and enterovirus), and (5) a new genus, with a new and unique representative: the echovirus 22. Our analysis also allowed us, for the first time, to propose the most probable sequence of speciation events to have given rise to the current picornavirus family.The bootstrap procedure was used to check the reliability of the phylogenetic trees obtained. The application of the method of the statistical geometry in distance space to internal branches of the tree revealed a high degree of evolutionary “noise,” which makes the resolution of some internal branching points difficult.

Supervised Neural Networks for Clustering Conditions in DNA Array Data After Reducing Noise by Clustering Gene Expression Profiles

In this paper we compare various applications of supervised and unsupervised neural networks to the analysis of the gene expression profiles produced using DNA microarrays. In particular we are interested in the classification of samples or conditions. We have found that if gene expression profiles are clustered at the optimal level, the classification of conditions obtained using the average gene expression profile of each cluster is better than that obtained directly using all the gene expression profiles. If a supervised method (a back propagation neural network) is used instead of an unsupervised method, the efficiency of the classification of conditions increases. We studied the relative efficiencies of different clustering methods for reducing the dimensionality of the gene expression profile data set and found that the Self-Organising Tree Algorithm (SOTA) is a good choice for this task.

A Method for Determining the Position and Size of Optimal Sequence Regions for Phylogenetic Analysis

The availability of fast and accurate sequencing procedures along with the use of PCR has led to a proliferation of studies of variability at the molecular level in populations. Nevertheless, it is often impractical to examine long genomic stretches and a large number of individuals at the same time. In order to optimize this kind of study, we suggest a heuristic procedure for detection of the shortest region whose informational content can be considered sufficient for significant phylogenetic reconstruction. The method is based on the comparison of the pairwise genetic distances obtained from a set of sequences of reference to those obtained for different windows of variable size and position by means of a simple index. We also present an approach for testing whether the informative content in the stretches selected in this way is significantly different from the corresponding content shown by the larger genomic regions used as reference. Application of this test to the analysis of the VP1 protein gene of foot-and-mouth-disease type C virus allowed us to define optimal stretches whose informative content is not significantly different from that displayed by the complete VP1 sequence. We showed that the predictions made for type C sequences are valid for type O sequences, indicating that the results of the procedure are consistent.