Population structure and recent evolution of Plasmodium falciparum

Plasmodium falciparum is the agent of malignant malaria, one of mankind's most severe maladies. The parasite exhibits antigenic polymorphisms that have been postulated to be ancient. We have proposed that the extant world populations of P. falciparum have derived from one single parasite, a cenancestor, within the last 5,000–50,000 years. This inference derives from the virtual or complete absence of synonymous nucleotide polymorphisms at genes not involved in immune or drug responses. Seeking to conciliate this claim with extensive antigenic polymorphism, we first note that allele substitutions or polymorphisms can arise very rapidly, even in a single generation, in large populations subject to strong natural selection. Second, new alleles can arise not only by single-nucleotide mutations, but also by duplication/deletion of short simple-repeat DNA sequences, a process several orders of magnitude faster than single-nucleotide mutation. We analyze three antigenic genes known to be extremely polymorphic: Csp, Msp-1, and Msp-2. We identify regions consisting of tandem or proximally repetitive short DNA sequences, including some previously unnoticed. We conclude that the antigenic polymorphisms are consistent with the recent origin of the world populations of P. falciparum inferred from the analysis of nonantigenic genes.

Rates of nucleotide substitution in sexual and anciently asexual rotifers

The class Bdelloidea of the phylum Rotifera is the largest well studied eukaryotic taxon in which males and meiosis are unknown, and the only one for which these indications of ancient asexuality are supported by cytological and molecular genetic evidence. We estimated the rates of synonymous and nonsynonymous substitutions in the hsp82 heat shock gene in bdelloids and in facultatively sexual rotifers of the class Monogononta, employing distance based and maximum likelihood methods. Relative-rate tests, using acanthocephalan rotifers as an outgroup, showed slightly higher rates of nonsynonymous substitution and slightly lower rates of synonymous substitution in bdelloids as compared with monogononts. The opposite trend, however, was seen in intraclass pairwise comparisons. If, as it seems, bdelloids have evolved asexually, an equality of bdelloid and monogonont substitution rates would suggest that the maintenance of sexual reproduction in monogononts is not attributable to an effect of sexual reproduction in limiting the load of deleterious nucleotide substitutions.

Pattern of mutation in the genome of influenza A virus on adaptation to increased virulence in the mouse lung: Identification of functional themes

The genetic basis for virulence in influenza virus is largely unknown. To explore the mutational basis for increased virulence in the lung, the H3N2 prototype clinical isolate, A/HK/1/68, was adapted to the mouse. Genomic sequencing provided the first demonstration, to our knowledge, that a group of 11 mutations can convert an avirulent virus to a virulent variant that can kill at a minimal dose. Thirteen of the 14 amino acid substitutions (93%) detected among clonal isolates were likely instrumental in adaptation because of their positive selection, location in functional regions, and/or independent occurrence in other virulent influenza viruses. Mutations in virulent variants repeatedly involved nuclear localization signals and sites of protein and RNA interaction, implicating them as novel modulators of virulence. Mouse-adapted variants with the same hemagglutinin mutations possessed different pH optima of fusion, indicating that fusion activity of hemagglutinin can be modulated by other viral genes. Experimental adaptation resulted in the selection of three mutations that were in common with the virulent human H5N1 isolate A/HK/156/97 and that may be instrumental in its extreme virulence. Analysis of viral adaptation by serial passage appears to provide the identification of biologically relevant mutations.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as “bags of genes” and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

Structural basis for fibroblast growth factor receptor 2 activation in Apert syndrome

Apert syndrome (AS) is characterized by craniosynostosis (premature fusion of cranial sutures) and severe syndactyly of the hands and feet. Two activating mutations, Ser-252 → Trp and Pro-253 → Arg, in fibroblast growth factor receptor 2 (FGFR2) account for nearly all known cases of AS. To elucidate the mechanism by which these substitutions cause AS, we determined the crystal structures of these two FGFR2 mutants in complex with fibroblast growth factor 2 (FGF2) . These structures demonstrate that both mutations introduce additional interactions between FGFR2 and FGF2, thereby augmenting FGFR2–FGF2 affinity. Moreover, based on these structures and sequence alignment of the FGF family, we propose that the Pro-253 → Arg mutation will indiscriminately increase the affinity of FGFR2 toward any FGF. In contrast, the Ser-252 → Trp mutation will selectively enhance the affinity of FGFR2 toward a limited subset of FGFs. These predictions are consistent with previous biochemical data describing the effects of AS mutations on FGF binding. Alterations in FGFR2 ligand affinity and specificity may allow inappropriate autocrine or paracrine activation of FGFR2. Furthermore, the distinct gain-of-function interactions observed in each crystal structure provide a model to explain the phenotypic variability among AS patients.

Activation of the Maize Anthocyanin Gene a2 Is Mediated by an Element Conserved in Many Anthocyanin Promoters1

Two transcription factors, C1 (a Myb-domain protein) and B (a basic-helix-loop-helix protein), mediate transcriptional activation of the anthocyanin-biosynthetic genes of maize (Zea mays). To begin to assess the mechanism of activation, the sequences required for C1- and B-mediated induction have been determined for the a2 promoter, which encodes an anthocyanin-biosynthetic enzyme. Analysis of a series of 7- to 13-base-pair substitutions revealed two regions crucial for activation. One region, centered at −99, contained a C1-binding site that abolished C1 binding. The other crucial region was adjacent, centered at −91. C1 binding was not detected at this site, and mutation of this site did not prevent C1 binding at −99. An oligonucleotide dimer containing these two crucial elements was sufficient for C1 and B activation of a heterologous promoter. These data suggest that activation of the anthocyanin genes involves C1 and another factor binding at closely adjacent sites. Mutating a previously postulated anthocyanin consensus sequence within a2 did not significantly reduce activation by C1 and B. However, sequence comparisons of the crucial a2 regions with sequences important for C1- and B-mediated activation in two other anthocyanin promoters led to a revised consensus element shared by these promoters.

Deficiency of rds/peripherin causes photoreceptor death in mouse models of digenic and dominant retinitis pigmentosa

Retinitis pigmentosa (RP) is a group of inherited blinding diseases caused by mutations in multiple genes including RDS. RDS encodes rds/peripherin (rds), a 36-kDa glycoprotein in the rims of rod and cone outer-segment (OS) discs. Rom1 is related to rds with similar membrane topology and the identical distribution in OS. In contrast to RDS, no mutations in ROM1 alone have been associated with retinal disease. However, an unusual digenic form of RP has been described. Affected individuals in several families were doubly heterozygous for a mutation in RDS causing a leucine 185 to proline substitution in rds (L185P) and a null mutation in ROM1. Neither mutation alone caused clinical abnormalities. Here, we generated transgenic/knockout mice that duplicate the amino acid substitutions and predicted levels of rds and rom1 in patients with RDS-mediated digenic and dominant RP. Photoreceptor degeneration in the mouse model of digenic RP was faster than in the wild-type and monogenic controls by histological, electroretinographic, and biochemical analysis. We observed a positive correlation between the rate of photoreceptor loss and the extent of OS disorganization in mice of several genotypes. Photoreceptor degeneration in RDS-mediated RP appears to be caused by a simple deficiency of rds and rom1. The critical threshold for the combined abundance of rds and rom1 is ≈60% of wild type. Below this value, the extent of OS disorganization results in clinically significant photoreceptor degeneration.

Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius

Microbes whose genomes are encoded by DNA and for which adequate information is available display similar genomic mutation rates (average 0.0034 mutations per chromosome replication, range 0.0025 to 0.0046). However, this value currently is based on only a few well characterized microbes reproducing within a narrow range of environmental conditions. In particular, no genomic mutation rate has been determined either for a microbe whose natural growth conditions may extensively damage DNA or for any member of the archaea, a prokaryotic lineage deeply diverged from both bacteria and eukaryotes. Both of these conditions are met by the extreme thermoacidophile Sulfolobus acidocaldarius. We determined the genomic mutation rate for this species when growing at pH 3.5 and 75°C based on the rate of forward mutation at the pyrE gene and the nucleotide changes identified in 101 independent mutants. The observed value of about 0.0018 extends the range of DNA-based microbes with rates close to the standard rate simultaneously to an archaeon and to an extremophile whose cytoplasmic pH and normal growth temperature greatly accelerate the spontaneous decomposition of DNA. The mutations include base pair substitutions (BPSs) and additions and deletions of various sizes, but the S. acidocaldarius spectrum differs from those of other DNA-based organisms in being relatively poor in BPSs. The paucity of BPSs cannot yet be explained by known properties of DNA replication or repair enzymes of Sulfolobus spp. It suggests, however, that molecular evolution per genome replication may proceed more slowly in S. acidocaldarius than in other DNA-based organisms examined to date.

An AU-box motif upstream of the SD sequence of light-dependent psbA transcripts confers mRNA instability in darkness in cyanobacteria

The psbA2 gene of a unicellular cyanobacterium, Microcystis aeruginosa K-81, encodes a D1 protein homolog in the reaction center of photosynthetic Photosystem II. The expression of the psbA2 transcript has been shown to be light-dependent as assessed under light and dark (12/12 h) cycling conditions. We aligned the 5′-untranslated leader regions (UTRs) of psbAs from different photosynthetic organisms and identified a conserved sequence, UAAAUAAA or the ‘AU-box’, just upstream of the SD sequences. To clarify the role of 5′-upstream cis-elements containing the AU-box for light-dependent expression of psbA2, a series of deletion and point mutations in the region were introduced into the genome of heterologous cyanobacterium Synechococcus sp. strain PCC 7942, and psbA2 expression was examined. A clear pattern of light-dependent expression was observed in recombinant cyanobacteria carrying the K-81 psbA2 –38/+36 region (which includes the minimal promoter element and a light-dependent cis-element with the AU-box), +1 indicating the transcription start site. A constitutive pattern of expression, in which the transcripts remained almost stable under dark conditions, was obtained in cells harboring the –38/+14 region (the minimal element), indicating that the +14/+36 region with the AU-box is important for the observed light-dependent expression. Point mutations analyses within the AU-box also revealed that changes in number, direction and identity (as assayed by adenine/uridine nucleotide substitutions) influenced the light-dependent pattern of expression. The level of psbA2 transcripts increased markedly in CG- or deletion-box mutants in the dark, strongly indicating that the AU- (AT-) box acts as a negative cis-element. Furthermore, characterization of transcript accumulation in cells treated with rifampicin suggests that psbA2 5′-mRNA is unstable in the dark, supporting the view that the light-dependent expression is controlled at the post-transcriptional level. We discuss various mechanisms that may lead to altered mRNA stability such as the binding of factor(s) or ribosomes to the 5′-UTR and possible roles of the AU-box motif and the SD sequence.

Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

Human brain contains high levels of heteroplasmy in the noncoding regions of mitochondrial DNA.

We have analyzed the level of intraindividual sequence variability (heteroplasmy) of mtDNA in human brain by denaturing gradient gel electrophoresis and sequencing. Single base substitutions, as well as insertions or deletions of single bases, were numerous in the noncoding control region (D-loop), and 35-45% of the molecules from a single tissue showed sequence differences. By contrast, heteroplasmy in coding regions was not detected. The lower level of heteroplasmy in the coding regions is indicative of selection against deleterious mutations. Similar levels of heteroplasmy were found in two brain regions from the same individual, while no heteroplasmy was detected in blood. Thus, heteroplasmy seems to be more frequent in nonmitotic tissues. We observed a 7.7-fold increase in the frequency of deletions/insertions and a 2.2-fold increase in the overall frequency of heteroplasmic mutations in two individuals aged 96 and 99, relative to an individual aged 28. Our results show that intraindividual sequence variability occurs at a high frequency in the noncoding regions of normal human brain and indicate that small insertions and deletions might accumulate with age at a lower rate than large rearrangements.

Determinants of RNA polymerase alpha subunit for interaction with beta, beta', and sigma subunits: hydroxyl-radical protein footprinting.

Escherichia coli RNA polymerase (RNAP) alpha subunit serves as the initiator for RNAP assembly, which proceeds according to the pathway 2 alpha-->alpha 2-->alpha 2 beta-->alpha 2 beta beta'-->alpha 2 beta beta' sigma. In this work, we have used hydroxyl-radical protein footprinting to define determinants of alpha for interaction with beta, beta', and sigma. Our results indicate that amino acids 30-75 of alpha are protected from hydroxyl-radical-mediated proteolysis upon interaction with beta (i.e., in alpha 2 beta, alpha 2 beta beta', and alpha 2 beta beta' sigma), and amino acids 175-210 of alpha are protected from hydroxyl-radical-mediated proteolysis upon interaction with beta' (i.e., in alpha 2 beta beta' and alpha 2 beta beta' sigma). The protected regions are conserved in the alpha homologs of prokaryotic, eukaryotic, archaeal, and chloroplast RNAPs and contain sites of substitutions that affect RNAP assembly. We conclude that the protected regions define determinants of alpha for direct functional interaction with beta and beta'. The observed maximal magnitude of protection upon interaction with beta and the observed maximal magnitude of protection upon interaction with beta' both correspond to the expected value for complete protection of one of the two alpha protomers of RNAP (i.e., 50% protection). We propose that only one of the two alpha protomers of RNAP interacts with beta and that only one of the two alpha protomers of RNAP interacts with beta'.

Random mutagenesis of Thermus aquaticus DNA polymerase I: concordance of immutable sites in vivo with the crystal structure.

Expression of Thermus aquaticus (Taq) DNA polymerase I (pol I) in Escherichia, coli complements the growth defect caused by a temperature-sensitive mutation in the host pol I. We replaced the nucleotide sequence encoding amino acids 659-671 of the O-helix of Taq DNA pol I, corresponding to the substrate binding site, with an oligonucleotide containing random nucleotides. Functional Taq pol I mutants were selected based on colony formation at the nonpermissive temperature. By using a library with 9% random substitutions at each of 39 positions, we identified 61 active Taq pol I mutants, each of which contained from one to four amino acid substitutions. Some amino acids, such as alanine-661 and threonine-664, were tolerant of several or even many diverse replacements. In contrast, no replacements or only conservative replacements were identified at arginine-659, lysine-663, and tyrosine-671. By using a library with totally random nucleotides at five different codons (arginine-659, arginine-660, lysine-663, phenylalanine-667, and glycine-668), we confirmed that arginine-659 and lysine-663 were immutable, and observed that only tyrosine substituted for phenylalanine-667. The two immutable residues and the two residues that tolerate only highly conservative replacements lie on the side of O-helix facing the incoming deoxynucleoside triphosphate, as determined by x-ray analysis. Thus, we offer a new approach to assess concordance of the active conformation of an enzyme, as interpreted from the crystal structure, with the active conformation inferred from in vivo function.

Identification of ciliary neurotrophic factor (CNTF) residues essential for leukemia inhibitory factor receptor binding and generation of CNTF receptor antagonists.

Ciliary neurotrophic factor (CNTF) drives the sequential assembly of a receptor complex containing the ligand-specific alpha-receptor subunit (CNTFR alpha) and the signal transducers gp130 and leukemia inhibitory factor receptor-beta (LIFR). The D1 structural motif, located at the beginning of the D-helix of human CNTF, contains two amino acid residues, F152 and K155, which are conserved among all cytokines that signal through LIFR. The functional importance of these residues was assessed by alanine mutagenesis. Substitution of either F152 or K155 with alanine was found to specifically inhibit cytokine interaction with LIFR without affecting binding to CNTFR alpha or gp130. The resulting variants behaved as partial agonists with varying degrees of residual bioactivity in different cell-based assays. Simultaneous alanine substitution of both F152 and K155 totally abolished biological activity. Combining these mutations with amino acid substitutions in the D-helix, which enhance binding affinity for the CNTFR alpha, gave rise to a potent competitive CNTF receptor antagonist. This protein constitutes a new tool for studies of CNTF function in normal physiology and disease.

Probing the pore region of recombinant N-methyl-D-aspartate channels using external and internal magnesium block.

Mg2+ ions block N-methyl-D-aspartate (NMDA) channels by entering the pore from either the extracellular or the cytoplasmic side of the membrane in a voltage-dependent manner. We have used these two different block phenomena to probe the structure of the subunits forming NMDA channels. We have made several amino acid substitutions downstream of the Q/R/N site in the TMII region of both NR1 and NR2A subunits. Mutant NR1 subunits were coexpressed with wild-type NR2A subunits and vice versa in Xenopus oocytes. We found that individually mutating the first two amino acid residues downstream to the Q/R/N site affects mostly the block by external Mg2+. Mutations of residues five to seven positions downstream of the Q/R/N site do not influence the external Mg2+ block, but clearly influence the block by internal Mg2+. These data add support to the hypothesis that there are two separate binding sites for external and internal Mg2+ block. They also indicate that the C-terminal end of TMII contributes to the inner vestibule of the pore of NMDA channels and thus provide additional evidence that TMII forms a loop that reemerges toward the cytoplasmic side of the membrane.