Nonsynonymous substitution in abalone sperm fertilization genes exceeds substitution in introns and mitochondrial DNA

Strong positive Darwinian selection acts on two sperm fertilization proteins, lysin and 18-kDa protein, from abalone (Haliotis). To understand the phylogenetic context for this dramatic molecular evolution, we obtained sequences of mitochondrial cytochrome c oxidase subunit I (mtCOI), and genomic sequences of lysin, 18-kDa, and a G protein subunit. Based on mtDNA differentiation, four north Pacific abalone species diverged within the past 2 million years (Myr), and remaining north Pacific species diverged over a period of 4–20 Myr. Between-species nonsynonymous differences in lysin and 18-kDa exons exceed nucleotide differences in introns by 3.5- to 24-fold. Remarkably, in some comparisons nonsynonymous substitutions in lysin and 18-kDa genes exceed synonymous substitutions in mtCOI. Lysin and 18-kDa intron/exon segments were sequenced from multiple red abalone individuals collected over a 1,200-km range. Only two nucleotide changes and two sites of slippage variation were detected in a total of >29,000 nucleotides surveyed. However, polymorphism in mtCOI and a G protein intron was found in this species. This finding suggests that positive selection swept one lysin allele and one 18-kDa allele to fixation. Similarities between mtCOI and lysin gene trees indicate that rapid adaptive evolution of lysin has occurred consistently through the history of the group. Comparisons with mtCOI molecular clock calibrations suggest that nonsynonymous substitutions accumulate 2–50 times faster in lysin and 18-kDa genes than in rapidly evolving mammalian genes.

Modulation of a calcium-sensitive nonspecific cation channel by closely associated protein kinase and phosphatase activities

Regulation of nonspecific cation channels often underlies neuronal bursting and other prolonged changes in neuronal activity. In bag cell neurons of Aplysia, it recently has been suggested that an intracellular messenger-induced increase in the activity of a nonspecific cation channel may underlie the onset of a 30-min period of spontaneous action potentials referred to as the “afterdischarge.” In patch clamp studies of the channel, we show that the open probability of the channel can be increased by an average of 10.7-fold by application of ATP to the cytoplasmic side of patches. Duration histograms indicate that the increase is primarily a result of a reduction in the duration and percentage of channel closures described by the slowest time constant. The increase in open probability was not observed using 5′-adenylylimidodiphosphate, a nonhydrolyzable ATP analog, and was blocked in the presence of H7 or the more specific calcium/phospholipid-dependent protein kinase C (PKC) inhibitor peptide(19–36). Because the increase in activity observed in response to ATP occurred without application of protein kinase, our results indicate that a kinase endogenous to excised patches mediates the effect. The effect of ATP could be reversed by exogenously applied protein phosphatase 1 or by a microcystin-sensitive phosphatase also endogenous to excised patches. These results, together with work demonstrating the presence of a protein tyrosine phosphatase in these patches, suggest that the cation channel is part of a regulatory complex including at least three enzymes. This complex may act as a molecular switch to activate the cation channel and, thereby, trigger the afterdischarge.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

Inhibition of telomerase by 2′-O-(2-methoxyethyl) RNA oligomers: effect of length, phosphorothioate substitution and time inside cells

2′-O-(2-methoxyethyl) (2′-MOE) RNA possesses favorable pharmocokinetic properties that make it a promising option for the design of oligonucleotide drugs. Telomerase is a ribonucleoprotein that is up-regulated in many types of cancer, but its potential as a target for chemotherapy awaits the development of potent and selective inhibitors. Here we report inhibition of human telomerase by 2′-MOE RNA oligomers that are complementary to the RNA template region. Fully complementary oligomers inhibited telomerase in a cell extract with IC50 values of 5–10 nM at 37°C. IC50 values for mismatch-containing oligomers varied with length and phosphorothioate substitution. After introduction into DU 145 prostate cancer cells inhibition of telomerase activity persisted for up to 7 days, equivalent to six population doublings. Inside cells discrimination between complementary and mismatch-containing oligomers increased over time. Our results reveal two oligomers as especially promising candidates for initiation of in vivo preclinical trials and emphasize that conclusions regarding oligonucleotide efficacy and specificity in cell extracts do not necessarily offer accurate predictions of activity inside cells.

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.

Estimation of evolutionary distances under stationary and nonstationary models of nucleotide substitution

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Evidence for the role of proteoglycans in cation-mediated gene transfer.

We report evidence that gene complexes, consisting of polycations and plasmid DNA enter cells via binding to membrane-associated proteoglycans. Treatment of HeLa cells with sodium chlorate, a potent inhibitor of proteoglycan sulfation, reduced luciferase expression by 69%. Cellular treatment with heparinase and chondroitinase ABC inhibited expression by 78% and 20% with respect to control cells. Transfection was dramatically inhibited by heparin and heparan sulfate and to a smaller extent by chondroitan sulfate B. Transfection of mutant, proteoglycan deficient Chinese hamster ovary cells was 53 x lower than of wild-type cells. For each of these assays, the intracellular uptake of DNA at 37 degrees C and the binding of DNA to the cell membrane at 4 degrees C was impaired. Preliminary transfection experiments conducted in mutant and wild-type Chinese hamster ovary cells suggest that transfection by some cationic lipids is also proteoglycan dependent. The variable distribution of proteoglycans among tissues may explain why some cell types are more susceptible to transfection than others.

Cation-pi interactions in aromatics of biological and medicinal interest: electrostatic potential surfaces as a useful qualitative guide.

The cation-pi interaction is an important, general force for molecular recognition in biological receptors. Through the sidechains of aromatic amino acids, novel binding sites for cationic ligands such as acetylcholine can be constructed. We report here a number of calculations on prototypical cation-pi systems, emphasizing structures of relevance to biological receptors and prototypical heterocycles of the type often of importance in medicinal chemistry. Trends in the data can be rationalized using a relatively simple model that emphasizes the electrostatic component of the cation-pi interaction. In particular, plots of the electrostatic potential surfaces of the relevant aromatics provide useful guidelines for predicting cation-pi interactions in new systems.

The weaver mutation of GIRK2 results in a loss of inwardly rectifying K+ current in cerebellar granule cells.

The weaver mutation in mice results in a severe ataxia that is attributable to the degeneration of cerebellar granule cells and dopaminergic neurons in the substantia nigra. Recent genetic studies indicate that the GIRK2 gene is altered in weaver. This gene codes for a G-protein-activated, inwardly rectifying K+ channel protein (8). The mutation results in a single amino acid substitution (glycine-->serine) in the pore-forming H5 region of the channel. The functional consequences of this mutation appear to depend upon the co-expression of other GIRK subunits--leading to either a gain or loss of function. Here, we show that G-protein-activated inwardly rectifying K+ currents are significantly reduced in cerebellar granule cells from animals carrying the mutant allele. The reduction is most pronounced in homozygous neurons. These findings suggest that the death of neurons in weaver is attributable to the loss of GIRK2-mediated currents, not to the expression of a nonspecific cation current.

Up-regulation of pressure-activated Ca(2+)-permeable cation channel in intact vascular endothelium of hypertensive rats.

In endothelial cells, stretch-activated cation channels have been proposed to act as mechanosensors for changes in hemodynamic forces. We have identified a novel mechanosensitive pressure-activated channel in intact endothelium from rat aorta and mesenteric artery. The 18-pS cation channel responded with a multifold increase in channel activity when positive pressure was applied to the luminal cell surface with the patch pipette and inactivated at negative pipette pressure. Channel permeability ratio for K+, Na+, and Ca2+ ions was 1:0.98:0.23. Ca2+ influx through the channel was sufficient to activate a neighboring Ca2(+)-dependent K+ channel. Hemodynamic forces are chronically disturbed in arterial hypertension. Endothelial cell dysfunction has been implicated in the pathogenesis of arterial hypertension. In two comparative studies, density of the pressure-activated channel was found to be significantly higher in spontaneously hypertensive rats and renovascular hypertensive rats compared with their respective normotensive controls. Channel activity presumably leads to mechanosensitive Ca2+ influx and induces cell hyperpolarization by K+ channel activity. Both Ca2+ influx and hyperpolarization are known to induce a vasodilatory endothelial response by stimulating endothelial nitric oxide (NO) production. Up-regulation of channel density in hypertension could, therefore, represent a counterregulatory mechanism of vascular endothelium.

Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL.

A number of studies have noted that nucleotide substitution rates at the chloroplast-encoded rbcL locus violate the molecular clock principle. Substitution rate variation at this plastid gene is particularly pronounced between palms and grasses; for example, a previous study estimated that substitution rates in rbcL sequences are approximately 5-fold faster in grasses than in palms. To determine whether a proportionate change in substitution rates also occurs in plant nuclear genes, we characterized nucleotide substitution rates in palm and grass sequences for the nuclear gene Adh. In this article, we report that palm sequences evolve at a rate of 2.61 x 10(-9) substitution per synonymous site per year, a rate which is slower than most plant nuclear genes. Grass Adh sequences evolve approximately 2.5-fold faster than palms at synonymous sites. Thus, synonymous rates in nuclear Adh genes show a marked decrease in palms relative to grasses, paralleling the pattern found at the plastid rbcL locus. This shared pattern indicates that synonymous rates are correlated between a nuclear and a plastid gene. Remarkably, nonsynonymous rates do not show this correlation. Nonsynonymous rates vary between two duplicated grass Adh loci, and nonsynonymous rates at the palm Adh locus are not markedly reduced relative to grasses.

Substitution rate calibration of small subunit ribosomal RNA identifies chlorarachniophyte endosymbionts as remnants of green algae.

Chlorarachniophytes are amoeboid algae with chlorophyll a and b containing plastids that are surrounded by four membranes instead of two as in plants and green algae. These extra membranes form important support for the hypothesis that chlorarachniophytes have acquired their plastids by the ingestion of another eukaryotic plastid-containing alga. Chlorarachniophytes also contain a small nucleus-like structure called the nucleomorph situated between the two inner and the two outer membranes surrounding the plastid. This nucleomorph is a remnant of the endosymbiont's nucleus and encodes, among other molecules, small subunit ribosomal RNA. Previous phylogenetic analyses on the basis of this molecule provided unexpected and contradictory evidence for the origin of the chlorarachniophyte endosymbiont. We developed a new method for measuring the substitution rates of the individual nucleotides of small subunit ribosomal RNA. From the resulting substitution rate distribution, we derived an equation that gives a more realistic relationship between sequence dissimilarity and evolutionary distance than equations previously available. Phylogenetic trees constructed on the basis of evolutionary distances computed by this new method clearly situate the chlorarachniophyte nucleomorphs among the green algae. Moreover, this relationship is confirmed by transversion analysis of the Chlorarachnion plastid small subunit ribosomal RNA.

A new isoleucine substitution of Val-20 in transthyretin tetramers selectively impairs dimer-dimer contacts and causes systemic amyloidosis.

The most frequent form of inherited amyloidoses is associated with mutations in the transthyretin (TTR) gene coding for 127-amino acid residues of four identical, noncovalently linked subunits that form a pair of dimers in the plasma protein complex. Amyloid fibrils containing the variant and to a lesser extent the wild-type form of the TTR molecule are deposited in various organs, including peripheral nerves and the myocardium, with polyneuropathy and cardiomyopathy as major clinical manifestations. So far, more than 40 distinct amino acid substitutions distributed throughout the TTR sequence over 30 positions have been found to be correlated with an increased amyloidogenicity of TTR. Most of these amyloidogenic amino acid substitutions are suspected to alter the conformation and stability of the monomer. Here we identify and characterize by protein and DNA analysis a novel amyloidogenic Val-20 to Ile mutation in a German three-generation family. The index patient suffered from severe amyloid cardiomyopathy at the age of 60. Conformational stability and unfolding behavior of the Ile-20 monomer in urea gradients was found to be almost indistinguishable from that of wild-type TTR. In contrast, tetramer stability was significantly reduced in agreement with the expected change in the interactions between the two opposing dimers via the side chain of Ile-20. Our observations provide strong evidence for the view that amyloidogenic amino acid substitutions in TTR facilitate the conversion of tetrameric TTR complexes into those conformational intermediates of the TTR folding pathway that have an intrinsic amyloidogenic potential.

Novel human DNA alkyltransferases obtained by random substitution and genetic selection in bacteria.

DNA repair alkyltransferases protect organisms against the cytotoxic, mutagenic, and carcinogenic effects of alkylating agents by transferring alkyl adducts from DNA to an active cysteine on the protein, thereby restoring the native DNA structure. We used random sequence substitutions to gain structure-function information about the human O6-methylguanine-DNA methyltransferase (EC 2.1.1.63), as well as to create active mutants. Twelve codons surrounding but not including the active cysteine were replaced by a random nucleotide sequence, and the resulting random library was selected for the ability to provide alkyltransferase-deficient Escherichia coli with resistance to the methylating agent N-methyl-N'-nitro-N-nitrosoguanidine. Few amino acid changes were tolerated in this evolutionarily conserved region of the protein. One mutation, a valine to phenylalanine change at codon 139 (V139F), was found in 70% of the selected mutants; in fact, this mutant was selected much more frequently than the wild type. V139F provided alkyltransferase-deficient bacteria with greater protection than the wild-type protein against both the cytotoxic and mutagenic effects of N-methyl-N'-nitro-N-nitrosoguanidine, increasing the D37 over 4-fold and reducing the mutagenesis rate 2.7-5.5-fold. This mutant human alkyltransferase, or others similarly created and selected, could be used to protect bone marrow cells from the cytotoxic side effects of alkylation-based chemotherapeutic regimens.