The 5' coding region of Paramecium surface antigen genes controls mutually exclusive transcription.

Paramecium tetraurelia stock 51 can express at least 11 different types of surface antigens, yet only a single type is expressed on the surface of an individual cell at any one time. The differential expression of stock 51 type A and B surface antigen genes (51A and 51B) is regulated at the level of transcription. Previously, we reported that nucleotide sequences upstream of position -26 (relative to the start of translation) in the 51A and 51B surface antigen genes are necessary for transcriptional activity but are not sufficient to direct differential transcriptional control. In this report we demonstrate that at least some of the critical elements necessary for differential transcription of the 51A and 51B genes lie within the 5' coding region. A hybrid gene that contains 51B upstream sequences (-475 to +1) attached to the ATG start codon of 51A is not cotranscribed with the 51B gene. In contrast, further substitution with 51B sequences (-1647 to +885) allows the chimeric gene to be coexpressed with 51B. A different hybrid gene containing a substitution of 51B sequence from -26 to +885 in the 51A gene is also coexpressed with 51B, revealing that the critical elements within the coding region of 51B do not require 51B upstream sequences for their effect. Coinjection of the 51A gene with the chimeric gene that contains 51B up to +885 showed that the same sequences that allow coexpression with 51B prevent cotranscription with 51A. Together, these results demonstrate that a region downstream of the transcriptional start site between nucleotide positions +1 and +885 (relative to translational start) is necessary to control differential transcriptional activity.

Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis.

Competing hypotheses seek to explain the evolution of oxygenic and anoxygenic processes of photosynthesis. Since chlorophyll is less reduced and precedes bacteriochlorophyll on the modern biosynthetic pathway, it has been proposed that chlorophyll preceded bacteriochlorophyll in its evolution. However, recent analyses of nucleotide sequences that encode chlorophyll and bacteriochlorophyll biosynthetic enzymes appear to provide support for an alternative hypothesis. This is that the evolution of bacteriochlorophyll occurred earlier than the evolution of chlorophyll. Here we demonstrate that the presence of invariant sites in sequence datasets leads to inconsistency in tree building (including maximum-likelihood methods). Homologous sequences with different biological functions often share invariant sites at the same nucleotide positions. However, different constraints can also result in additional invariant sites unique to the genes, which have specific and different biological functions. Consequently, the distribution of these sites can be uneven between the different types of homologous genes. The presence of invariant sites, shared by related biosynthetic genes as well as those unique to only some of these genes, has misled the recent evolutionary analysis of oxygenic and anoxygenic photosynthetic pigments. We evaluate an alternative scheme for the evolution of chlorophyll and bacteriochlorophyll.

Equine rhinovirus 1 is more closely related to foot-and-mouth disease virus than to other picornaviruses.

Equine rhinovirus 1 (ERhV1) is a respiratory pathogen of horses which has an uncertain taxonomic status. We have determined the nucleotide sequence of the ERhV1 genome except for a small region at the 5' end. The predicted polyprotein was encoded by 6741 nucleotides and possessed a typical picornavirus proteolytic cleavage pattern, including a leader polypeptide. The genomic structure and predicted amino acid sequence of ERhV1 were more similar to those of foot-and-mouth disease viruses (FMDVs), the only members of the aphthovirus genus, than to those of other picornaviruses. Features which were most similar to FMDV included a 16-amino acid 2A protein which was 87.5% identical in sequence of FMDV 2A, a leader (L) protein similar in size to FMDV Lab and the possibility of a truncated L protein similar in size to FMDV Lb, and a 3C protease which recognizes different cleavage sites. However, unlike FMDV, ERhV1 had only one copy of the 3B (VPg) polypeptide. The phylogenetic relationships of the ERhV1 sequence and nucleotide sequences of representative species of the five genera of the family Picornaviridae were examined. Nucleotide sequences coding for the complete polyprotein, the RNA polymerase, and VP1 were analyzed separately. The phylogenetic trees confirmed that ERhV1 was more closely related to FMDV than to other picornaviruses and suggested that ERhV1 may be a member, albeit very distant, of the aphthovirus genus.

Molecular basis of human mitochondrial very-long-chain acyl-CoA dehydrogenase deficiency causing cardiomyopathy and sudden death in childhood.

beta-Oxidation of long-chain fatty acids provides the major source of energy in the heart. Defects in enzymes of the beta-oxidation pathway cause sudden, unexplained death in childhood, acute hepatic encephalopathy or liver failure, skeletal myopathy, and cardiomyopathy. Very-long-chain acyl-CoA dehydrogenase [VLCAD; very-long-chain-acyl-CoA:(acceptor) 2,3-oxidoreductase, EC 1.3.99.13] catalyzes the first step in beta-oxidation. We have isolated the human VLCAD cDNA and gene and determined the complete nucleotide sequences. Polymerase chain reaction amplification of VLCAD mRNA and genomic exons defined the molecular defects in two patients with VLCAD deficiency who presented with unexplained cardiac arrest and cardiomyopathy. In one, a homozygous mutation in the consensus dinucleotide of the donor splice site (g+1-->a) was associated with universal skipping of the prior exon (exon 11). The second patient was a compound heterozygote, with a missense mutation, C1837-->T, changing the arginine at residue 613 to tryptophan on one allele and a single base deletion at the intron-exon 6 boundary as the second mutation. This initial delineation of human mutations in VLCAD suggests that VLCAD deficiency reduces myocardial fatty acid beta-oxidation and energy production and is associated with cardiomyopathy and sudden death in childhood.

Electrophoretic variation in adenylate kinase of Neisseria meningitidis is due to inter- and intraspecies recombination.

In prokaryotic and eukaryotic organisms, the electrophoretic variation in housekeeping enzymes from natural populations is assumed to have arisen by the accumulation of stochastic predominantly neutral mutations. In the naturally transformable bacterium Neisseria meningitidis, we show that variation in the electrophoretic mobility of adenylate kinase is due to inter- and intraspecies recombination rather than mutation. The nucleotide sequences of the adenylate kinase gene (adk) from isolates that express the predominant slow electrophoretic variant were rather uniform, differing in sequence at an average of 1.1% of nucleotide sites. The adk sequences of rare isolates expressing the fast migrating variant were identical to each other but had a striking mosaic structure when compared to the adk genes from strains expressing the predominant variant. Thus the sequence from the fast variants was identical to those of typical slow variants in the first 158 bp of the gene but differed by 8.4% in the rest of the gene (nt 159-636). The fast electrophoretic variant appears to have arisen by the replacement of most of the meningococcal gene with the corresponding region from the adk gene of a closely related Neisseria species. The adk genes expressing the electrophoretic variant with intermediate mobility were perfect, or almost perfect, recombinants between the adk genes expressing the fast and slow variants. Recombination may, therefore, play a major role in the generation of electrophoretically detectable variation in housekeeping enzymes of some bacterial species.

Relationship between evolutionary rate and cellular location among the Inv/Spa invasion proteins of Salmonella enterica.

For 21 strains of Salmonella enterica, nucleotide sequences were obtained for three invasion genes, spaO, spaP, and spaQ, of the chromosomal inv/spa complex, the products of which form a protein export system required for entry of the bacteria into nonphagocytic host cells. These genes are present in all eight subspecies of the salmonellae, and homologues occur in a variety of other bacteria, including the enteric pathogens Shigella and Yersinia, in which they are plasmid borne. Evolutionary diversification of the invasion genes among the subspecies of S. enterica has been generally similar in pattern and average rate to that of housekeeping genes. However, the range of variation in evolutionary rate among the invasion genes is unusually large, and there is a relationship between the evolutionary rate and cellular location of the invasion proteins, possibly reflecting diversifying selection on exported proteins in adaptation to variable host factors in extracellular environments. The SpaO protein, which is hypervariable in S. enterica and exhibits only 24% sequence identity with its homologues in Shigella and Yersinia, is secreted. In contrast, the membrane-associated proteins SpaP, SpaQ, and InvA are weakly polymorphic and have > 60% sequence identity with the corresponding proteins of other enteric bacteria. Acquisition of the inv/spa genes may have been a key event in the evolution of the salmonellae as pathogens, following which the invention of flagellar phase shifting facilitated niche expansion to include warm-blooded vertebrates.

An oleate 12-hydroxylase from Ricinus communis L. is a fatty acyl desaturase homolog.

Recent spectroscopic evidence implicating a binuclear iron site at the reaction center of fatty acyl desaturases suggested to us that certain fatty acyl hydroxylases may share significant amino acid sequence similarity with desaturases. To test this theory, we prepared a cDNA library from developing endosperm of the castor-oil plant (Ricinus communis L.) and obtained partial nucleotide sequences for 468 anonymous clones that were not expressed at high levels in leaves, a tissue deficient in 12-hydroxyoleic acid. This resulted in the identification of several cDNA clones encoding a polypeptide of 387 amino acids with a predicted molecular weight of 44,407 and with approximately 67% sequence homology to microsomal oleate desaturase from Arabidopsis. Expression of a full-length clone under control of the cauliflower mosaic virus 35S promoter in transgenic tobacco resulted in the accumulation of low levels of 12-hydroxyoleic acid in seeds, indicating that the clone encodes the castor oleate hydroxylase. These results suggest that fatty acyl desaturases and hydroxylases share similar reaction mechanisms and provide an example of enzyme evolution.

Interaction with the recombination hot spot chi in vivo converts the RecBCD enzyme of Escherichia coli into a chi-independent recombinase by inactivation of the RecD subunit.

The RecBCD enzyme of Escherichia coli promotes recombination preferentially at chi nucleotide sequences and has in vivo helicase and strong duplex DNA exonuclease (exoV) activities. The enzyme without the RecD subunit, as in a recD null mutant, promotes recombination efficiently but independently of chi and has no nucleolytic activity. Employing phage lambda red gam crosses, phage T4 2- survival measurements, and exoV assays, it is shown that E. coli cells in which RecBCD has extensive opportunity to interact with linear chi-containing DNA (produced by rolling circle replication of a plasmid with chi or by bleomycin-induced fragmentation of the cellular chromosome) acquire the phenotype of a recD mutant and maintain this for approximately 2 h. It is concluded that RecBCD is converted into RecBC during interaction with chi by irreversible inactivation of RecD. After conversion, the enzyme is released and initiates recombination on other DNA molecules in a chi-independent fashion. Overexpression of recD+ (from a plasmid) prevented the phenotypic change and providing RecD after the change restored chi-stimulated recombination. The observed recA+ dependence of the downregulation of exoV could explain the previously noted "reckless" DNA degradation of recA mutants. It is proposed that chi sites are regulatory elements for the RecBCD to RecBC switch and thereby function as cis- and trans-acting stimulators of RecBC-dependent recombination.

Evolutionary origin of Plasmodium and other Apicomplexa based on rRNA genes.

We have explored the evolutionary history of the Apicomplexa and two related protistan phyla, Dinozoa and Ciliophora, by comparing the nucleotide sequences of small subunit ribosomal RNA genes. We conclude that the Plasmodium lineage, to which the malarial parasites belong, diverged from other apicomplexan lineages (piroplasmids and coccidians) several hundred million years ago, perhaps even before the Cambrian. The Plasmodium radiation, which gave rise to several species parasitic to humans, occurred approximately 129 million years ago; Plasmodium parasitism of humans has independently arisen several times. The origin of apicomplexans (Plasmodium), dinoflagellates, and ciliates may be > 1 billion years old, perhaps older than the three multicellular kingdoms of animals, plants, and fungi. Digenetic parasitism independently evolved several times in the Apicomplexa.

Accelerated evolution in the protein-coding regions is universal in crotalinae snake venom gland phospholipase A2 isozyme genes.

The nucleotide sequences of four genes encoding Trimeresurus gramineus (green habu snake, crotalinae) venom gland phospholipase A2 (PLA2; phosphatidylcholine 2-acylhydrolase, EC 3.1.1.4) isozymes were compared internally and externally with those of six genes encoding Trimeresurus flavoviridis (habu snake, crotalinae) venom gland PLA2 isozymes. The numbers of nucleotide substitutions per site (KN) for the noncoding regions including introns were one-third to one-eighth of the numbers of nucleotide substitutions per synonymous site (KS) for the protein-coding regions of exons, indicating that the noncoding regions are much more conserved than the protein-coding regions. The KN values for the introns were found to be nearly equivalent to those of introns of T. gramineus and T. flavoviridis TATA box-binding protein genes, which are assumed to be a general (nonvenomous) gene. Thus, it is evident that the introns of venom gland PLA2 isozyme genes have evolved at a similar rate to those of nonvenomous genes. The numbers of nucleotide substitutions per nonsynonymous site (KA) were close to or larger than the KS values for the protein-coding regions in venom gland PLA2 isozyme genes. All of the data combined reveal that Darwinian-type accelerated evolution has universally occurred only in the protein-coding regions of crotalinae snake venom PLA2 isozyme genes.

Preferential isolation of DNA fragments associated with CpG islands.

We describe a procedure for preferential isolation of DNA fragments with G+C-rich portions. Such fragments occur in known genes within or adjacent to CpG islands. Since about 56% of human genes are associated with CpG islands, isolation of these fragments permits detection and probing of many genes within much larger segments of DNA, such as cosmids or yeast artificial chromosomes, which have not been sequenced. Cloned DNA fragments digested with four restriction endonucleases were subjected to denaturing gradient gel electrophoresis. Long G+C-rich sections in fragments inhibit strand dissociation after the fragments reach retardation level in the gradient; such fragments are retained in the gel after most others disappear. Nucleotide sequences of the retained fragments show that about half of these fragments appear to be derived from CpG islands. Northern analysis indicated the presence of RNA complementary to most of the retained fragments. A heuristic approach to the relation between base sequence and the kinetics of strand dissociation of partly melted molecules appears to account for retention and nonretention. The expectation that CpG island fragments will be enriched among fragments retained in a denaturing gradient is supported by rate estimates based on melting theory applied to known sequences. This method, designated SPM for segregation of partly melted molecules, is expected to provide a means for convenient and efficient isolation of genes from unsequenced DNA.

Identification of two flavivirus-like genomes in the GB hepatitis agent.

A subtractive PCR methodology known as representational difference analysis was used to clone specific nucleotide sequences present in the infectious plasma from a tamarin infected with the GB hepatitis agent. Eleven unique clones were identified, seven of which were examined extensively. All seven clones appeared to be derived from sequences exogenous to the genomes of humans, tamarins, Saccharomyces cerevisiae, and Escherichia coli. In addition, sequences from these clones were not detected in plasma or liver tissue of tamarins prior to their inoculation with the GB agent. These sequences were detected by reverse transcription-PCR in acute-phase plasma of tamarins inoculated with the GB agent. Probes derived from two of the seven clones detected an RNA species of > or = 8.3 kb in the liver of a GB-agent-infected tamarin by Northern blot hybridization. Sequence analysis indicated that five of the seven clones encode polypeptides that possess limited amino acid identity with the nonstructural proteins of hepatitis C virus. Extension of the sequences found in the seven clones revealed that plasma from an infected tamarin contained two RNA molecules > 9 kb long. Limited sequence identity with various isolates of hepatitis C virus and the relative positions of putative RNA helicases and RNA-dependent RNA polymerases in the predicted protein products of these molecules suggested that the GB agent contains two unique flavivirus-like genomes.

Presence of activatable Shiga toxin genotype (stx2d) in Shiga toxigenic Escherichia coli from livestock sources

Stx2d is a recently described Shiga toxin whose cytotoxicity is activated 10- to 1,000-fold by the elastase present in mouse or human intestinal mucus. We examined Shiga toxigenic Escherichia coli (STEC) strains isolated from food and livestock sources for the presence of activatable stx(2d). The stx(2) operons of STEC were first analyzed by PCR-restriction fragment length polymorphism (RFLP) analysis and categorized as stx(2), stx(2c) (vha), stx(2c) (vhb), or stx(2d) (EH250). Subsequently, the stx(2c) (vha) and stx(2c) (vhb) operons were screened for the absence of a PstI site in the stx(2a) subunit gene, a restriction site polymorphism which is a predictive indicator for the stx(2d) (activatable) genotype. Twelve STEC isolates carrying putative stx(2d) operons were identified, and nucleotide sequencing was used to confirm the identification of these operons as stx(2d). The complete nucleotide sequences of seven representative stx(2d) operons were determined. Shiga toxin expression in stx(2d) isolates was confirmed by immunoblotting. stx(2d) isolates were induced for the production of bacteriophages carrying stx. Two isolates were able to produce bacteriophages phi1662a and phi1720a carrying the stx(2d) operons. RFLP analysis of bacteriophage genomic DNA revealed that phi1662a and phi1720a were highly related to each other; however, the DNA sequences of these two stx(2d) operons were distinct. The STEC strains carrying these operons were isolated from retail ground beef. Surveillance for STEC strains expressing activatable stx(2d) Shiga toxin among clinical cases may indicate the significance of this toxin subtype to human health.

HBV core sequence: definition of genotype-specific variability and correlation with geographical origin

There are eight genotypes and nine subtypes of HBV. Small differences in geographical origin are associated with sequence changes in the surface gene. Here, we compared core gene sequences from different genotypes and geographical regions. Specific combinations of 24 amino acid substitutions at nine residues allowed allocation of a sequence to a subtype. Six of these nine residues were located in different T cell epitopes depending on HBV geographical area and/or genotype. Thirty-seven nucleotide changes were associated uniquely with specific genotypes and subtypes. Unique amino acid and nucleotide variants were found in a majority of sequences from specific countries as well as within subtype ayw2 and adr. Specific nucleotide motifs were defined for Korean, Indian, Chinese, Italian and Pacific region isolates. Finally, we observed amino acid motifs that were common to either South-east Asian or Western populations, irrespective of subtype. We believe that HBV strains spread within constrained ethnic groups, result in selection pressures that define sequence variability within each subtype. It suggests that particular T cell epitopes are specific for geographical regions, and thus ethnic groups; this may affect the design of immunomodulatory therapies.

Nucleocapsid gene variability reveals two subgroups of Lettuce necrotic yellows virus

The complete nucleocapsid (N) genes of eight Australian isolates of Lettuce necrotic yellows virus (LNYV) were amplified by reverse transcription PCR, cloned and sequenced. Phylogenetic analyses of these sequences revealed two distinct subgroups of LNYV isolates. Nucleotide sequences within each subgroup were more than 96% identical but heterogeneity between groups was about 20% at the nucleotide sequence level. However, less than 4% heterogeneity was noted at the amino acid level, indicating mostly third nucleotide position changes and a strong conservation for N protein function. There was no obvious geographical or temporal separation of the subgroups in Australia.