Associations of epithelial sodium channel genes with blood pressure: the GenSalt study.

In order to investigate the associations of SCNN1A, SCNN1G and SCNN1B genes with blood pressure (BP) in the Han Chinese population, we included 2880 participants did not use antihypertensive medication in the month prior to the baseline survey in the current analysis. Forty-four tag-single-nucleotide polymorphisms (SNPs) in epithelial sodium channel (ENaC) genes were selected and genotyped, and nine BP measurements were obtained during the 3-day examination. In the single-marker analyses, we identified significant associations of SCNN1A marker rs13306613 with diastolic BP (DBP) and SCNN1B marker rs12447134 with systolic BP (SBP) under codominant model after Bonferroni's correction (P=2.82 × 10(-5) and 4.63 × 10(-4), respectively). In addition, five SNPs in SCNN1G and four SNPs in SCNN1B achieved nominal significance for SBP, DBP or mean arterial pressure (MAP) under the additive model. For example, the minor C allele of rs5735 in SCNN1G gene was associated with decreased SBP, DBP and MAP (P=0.016, 5.41 × 10(-3) and 4.36 × 10(-3), respectively). Gene-based results showed significant associations of SCNN1G and SCNN1B with BP levels. This study suggested that ENaC genes have important roles in BP regulation in the Han Chinese population. Future studies are warranted to replicate these findings, and functional studies are needed to identify true causal variants in ENaC genes.

Contrasting levels of DNA polymorphism at the autosomal and X-linked visual color pigment loci in humans and squirrel monkeys.

The X-linked color pigment (opsin) locus is known to be highly polymorphic in the squirrel monkey and other New World monkeys. To see whether this is also the case for the autosomal (blue) opsin locus, we obtained 32 squirrel monkey and 30 human blue opsin gene sequences. No amino acid polymorphism was found in either the squirrel monkey sample or the human sample, contrary to the situation at the X-linked opsin locus. This sharp contrast in the level of polymorphism might be due to differences in gene expression between the autosomal and the X-linked loci. At the X-linked locus, heterozygote advantage can occur because, owing to X-inactivation, the two alleles in a heterozygote are expressed in different cone cells, producing two types of cone cell, whereas at the autosomal locus, heterozygote advantage cannot occur because the two alleles in a heterozygote are expressed in the same cone cells, producing only one type of cone cell (i.e., phenotypically a homozygote). From the sequence data, the levels of nucleotide diversity (pi, i.e., the number of nucleotide differences per site) are estimated: for the human sample, pi = 0.00% per nondegenerate site, 0.00% per twofold degenerate site, and 0.04% per fourfold degenerate site in the coding regions and 0.01% per site in intron 4; for the squirrel monkey sample, pi = 0.00% per nondegenerate site, 0.00% per twofold degenerate site, and 0.15% per fourfold degenerate site in the coding regions and 0.17% per site in intron 4. The blue opsin genes from the common and pygmy chimpanzees, the gorilla, the capuchin, and the howler monkey were also sequenced. Features critical to the function of the opsin are well conserved in all known mammalian sequences. However, the interhelical loops are, on average, actually more conservative than the transmembrane helical regions. In addition, these sequence data and those from some other genes indicate that the common and pygmy chimpanzees are not closely related, their divergence data being from one third to one half the date of the human-chimpanzee divergence.

Mutation pattern variation among regions of the primate genome.

We sequenced three argininosuccinate-synthetase-processed pseudogenes (PsiAS-A1, PsiAS-A3, PsiAS-3) and their noncoding flanking sequences in human, orangutan, baboon, and colobus. Our data showed that these pseudogenes were incorporated into the genome of the Old World monkeys after the divergence of the Old World and New World monkey lineages. These pseudogene flanking regions show variable mutation rates and patterns. The variation in the G/C to A/T mutation rate (u) can account for the unequal GC contents at equilibrium: 34.9, 36.9, and 41.7% in the pseudogene PsiAS-A1, PsiAS-A3, and PsiAS-3 flanking regions, respectively. The A/T to G/C mutation rate (v) seems stable and the u/v ratios equal 1.9, 1.7, and 1.4 in the flanking regions of PsiAS-A1, PsiAS-A3, and PsiAS-3, respectively. These "regional" variations of the mutation rate affect the evolution of the pseudogenes, too. The ratio u/v being greater than 1.0 in each case, the overall mutation rate in the GC-rich pseudogenes is, as expected, higher than in their GC-poor flanking regions. Moreover, a "sequence effect" has been found. In the three cases examined u and v are higher (at least 20%) in the pseudogene than in its flanking region-i.e., the pseudogene appears as mutation "hot" spots embedded in "cold" regions. This observation could be partly linked to the fact that the pseudogene flanking regions are long-standing unconstrained DNA sequences, whereas the pseudogenes were relieved of selection on their coding functions only around 30-40 million years ago. We suspect that relatively more mutable sites maintained unchanged during the evolution of the argininosuccinate gene are able to change in the pseudogenes, such sites being eliminated or rare in the flanking regions which have been void of strong selective constraints over a much longer period. Our results shed light on (1) the multiplicity of factors that tune the spontaneous mutation rate and (2) the impact of the genomic position of a sequence on its evolution.

Sequences and evolution of human and squirrel monkey blue opsin genes.

The sequences of the entire blue opsin gene in the squirrel monkey (Saimiri boliviensis) and the five introns of the human blue opsin gene were obtained. Intron 3 of these genes contains an Alu sequence and intron 4 contains a partial mer13 sequence. A comparison of the squirrel monkey opsin sequence with published mammalian opsin sequences shows that features believed to be functionally critical are all conserved. However, the blue opsin has evolved twice as fast as rhodopsin and is only as conservative as the beta globin, which has evolved at the average rate of mammalian proteins. Interestingly, the interhelical loops are, on average, actually more conservative than the transmembrane alpha helical regions. The introns of the blue opsin gene have evolved at the average rate of introns in primate genes.

Evolutionary divergence and salinity-mediated selection in halophilic archaea.

Halophilic (literally salt-loving) archaea are a highly evolved group of organisms that are uniquely able to survive in and exploit hypersaline environments. In this review, we examine the potential interplay between fluctuations in environmental salinity and the primary sequence and tertiary structure of halophilic proteins. The proteins of halophilic archaea are highly adapted and magnificently engineered to function in an intracellular milieu that is in ionic balance with an external environment containing between 2 and 5 M inorganic salt. To understand the nature of halophilic adaptation and to visualize this interplay, the sequences of genes encoding the L11, L1, L10, and L12 proteins of the large ribosome subunit and Mn/Fe superoxide dismutase proteins from three genera of halophilic archaea have been aligned and analyzed for the presence of synonymous and nonsynonymous nucleotide substitutions. Compared to homologous eubacterial genes, these halophilic genes exhibit an inordinately high proportion of nonsynonymous nucleotide substitutions that result in amino acid replacement in the encoded proteins. More than one-third of the replacements involve acidic amino acid residues. We suggest that fluctuations in environmental salinity provide the driving force for fixation of the excessive number of nonsynonymous substitutions. Tinkering with the number, location, and arrangement of acidic and other amino acid residues influences the fitness (i.e., hydrophobicity, surface hydration, and structural stability) of the halophilic protein. Tinkering is also evident at halophilic protein positions monomorphic or polymorphic for serine; more than one-third of these positions use both the TCN and the AGY serine codons, indicating that there have been multiple nonsynonymous substitutions at these positions. Our model suggests that fluctuating environmental salinity prevents optimization of fitness for many halophilic proteins and helps to explain the unusual evolutionary divergence of their encoding genes.

Conserved sequence elements involved in regulation of ribosomal protein gene expression in halophilic archaea.

A region of the Haloferax volcanii genome encoding ribosomal proteins L11e, L1e, L10e, and L12e was cloned and sequenced, and the transcripts derived from the cluster were characterized. Flanking and noncoding regions of the sequence were analyzed phylogenetically by comparison with the homologous sequences from two other halophilic archaea, i.e., Halobacterium cutirubrum and Haloarcula marismortui. Motifs, identified by high-level sequence conservation, include both transcriptional and translational regulatory elements and other elements of unknown function.

Contrasting rates of nucleotide substitution in the X-linked and Y-linked zinc finger genes.

We have sequenced the entire exon (approximately 1.180 bp) encoding the zinc finger domain of the X-linked and Y-linked zinc finger genes (ZFX and ZFY, respectively) in the orangutan, the baboon, the squirrel monkey, and the rat; a total of 9,442 bp were sequenced. The ratio of the rates of synonymous substitution in the ZFY and ZFX genes is estimated to be 2.1 in primates. This is close to the ratio of 2.3 estimated from primate ZFY and ZFX intron sequences and supports the view that the male-to-female ratio of mutation rate in humans in considerably higher than 1 but not extremely large. The ratio of synonymous substitution rates in ZFY and ZFX is estimated to be 1.3 in the rat lineage but 4.2 in the mouse lineage. The former is close to the estimate (1.4) from introns. The much higher ratio in the mouse lineage (not statistically significant) might have arisen from relaxation of selective constraints. The synonymous divergence between mouse and rat ZFX is considerably lower than that between mouse and rat autosomal genes, agreeing with previous observations and providing some evidence for stronger selective constraints on synonymous changes in X-linked genes than in autosomal genes. At the protein level ZFX has been highly conserved in all placental mammals studied while ZFY has been well conserved in primates and foxes but has evolved rapidly in mice and rats, possibly due to relaxation of functional constraints as a result of the development of X-inactivation of ZFX in rodents. The long persistence of the ZFY-ZFX gene pair in mammals provides some insight into the process of degeneration of Y-linked genes.

Structure and transcription of the L11-L1-L10-L12 ribosomal protein gene operon from the extreme thermophilic archaeon Sulfolobus acidocaldarius.

We have cloned and sequenced four ribosomal protein genes from the extreme thermophilic archaeon Sulfolobus acidocaldarius P1. These genes code for proteins equivalent to L11, L1, L10 and L12 from Escherichia coli. The genes for the Sulfolobus L11, L1, L10 and L12 proteins are arranged in the same order as the equivalent genes in E. coli, i.e. L11-L1-L10-L12, and are transcribed as a single unit. Sequences resembling the consensus sequence for archaeal promoters have been detected upstream of the transcription initiation site. Transcription ends at several sites following a pyrimidine-rich region. The genes for proteins L11, L10 and L1 start with unusual initiation codons: GUG in the case of the L1 and L10 genes; and UUG in the case of L11. There are overlapping stop/start codons between the L11 and L1 genes, and between the L1 and L10, suggesting that the translation of the four genes might be coupled as in the bacteria.

Weak male-driven molecular evolution in rodents.

In humans and rodents the male-to-female ratio of mutation rate (alpha m) has been suggested to be extremely large, so that the process of nucleotide substitution is almost completely male-driven. However, our sequence data from the last intron of the X chromosome-linked (Zfx) and Y chromosome-linked (Zfy) zinc finger protein genes suggest that alpha m is only approximately 2 in rodents with a 95% confidence interval from 1 to 3. Moreover, from published data on oogenesis and spermatogenesis we estimate the male-to-female ratio of the number of germ cell divisions per generation to be approximately 2 in rodents, confirming our estimate of alpha m and suggesting that errors in DNA replication are the primary source of mutation. As the estimated alpha m for rodents is only one-third of our previous estimate of approximately 6 for higher primates, there appear to be generation-time effects--i.e., alpha m decreases with decreasing generation time.

Potential problems in estimating the male-to-female mutation rate ratio from DNA sequence data.

It is commonly believed that the rate of mutation is much higher in males than in females because the number of germ-cell divisions per generation is much larger in males than in females. However, the precise magnitude of the male-to-female mutation rate ratio (alpha m) remains unknown. Recently there have been efforts to estimate alpha m by using DNA sequence data from different species. We have studied the potential problems in such an approach. We found that the rate of synonymous substitution varies about fivefold among X-linked genes, as large as the variation among autosomal genes. This large variation makes the assumption of selective neutrality of synonymous changes dubious, so one should be cautious in using the synonymous rates in X-linked and autosomal genes to estimate alpha m. A similar difficulty was also observed in using nonhomologous intron sequences to estimate alpha m. Contrary to the expectation that X-linked sequences should evolve more slowly than autosomal sequences, the Alu repeat in the last intron of the X-linked zinc finger gene has evolved faster than the four autosomal Alu repeats used in this study. It appears that the best way to estimate alpha m is to use homologous sequences. However, such sequences may be involved in gene conversion events. In fact, we found evidence that the Y-linked and X-linked zinc finger genes have been involved in multiple conversion events during primate evolution. Thus, the possibility of gene conversion should be considered when using homologous sequences to estimate alpha m.

Male-driven evolution of DNA sequences.

It is commonly believed that the mutation rate is much higher in the human male germ line than in the female germ line because the number of germ-cell divisions per generation is much larger in males than in females. But direct estimation of mutation rates is difficult, relying mainly on sex-linked genetic diseases, so the ratio (alpha m) of male to female mutation rates is not clear. It has been noted that if alpha m is very large, then the rate of synonymous substitution in X-linked genes should be only 2/3 of that in autosomal genes, and comparison of human and rodent genes supported this prediction. As the number of X-linked genes used in the study was small and the X-linked and autosomal sequences were non-homologous, and given that the synonymous rate varies among genes, we sequenced the last intron (approximately 1 kb) of the Y-linked and X-linked zinc-finger-protein genes (ZFY and ZFX) in humans, orang-utans, baboons and squirrel monkeys. The ratio Y/X of the substitution rate in the Y-linked intron to that in the X-linked intron is approximately 2.3, which is close to that estimated from synonymous rates in the ZFY and ZFX genes and implies alpha m approximately 6. This estimate of alpha m supports the view that the evolution of DNA sequences in higher primates is male-driven. It is, however, much lower than the previous estimate and therefore raises a number of issues.

Comparison of the structure of archaebacterial ribosomal proteins equivalent to proteins L11 and L1 from Escherichia coli ribosomes.

The sequences of two ribosomal proteins from two widely divergent species of archaebacteria, Halobacterium cutirubrum and Sulfolobus solfataricus, have been deduced from the structure of their respective genes. These two proteins were found to be equivalent to the L11 and L1 ribosomal proteins of the eubacterium Escherichia coli. Sequence comparison revealed that the archaebacterial L11e (equivalent to E. coli L11) proteins are longer than the eubacterial protein due to a C-terminal extension of about 30 residues. The archaebacterial L11e proteins, like the E. coli L11, are rich in proline residues; most of these are conserved. L11 is the most highly methylated protein in the E. coli ribosome. However, sites of methylation are generally not conserved in the archaebacterial L11e proteins. The region of highest sequence similarity between L11 and the archaebacterial L11e proteins is the N-terminal domain. This domain is believed to interact with release factor 1 during termination of translation. The amino acid sequences of the archaebacterial L1e proteins were compared to the eubacterial E. coli L1 and Bacillus stearothermophilus L1e sequences. The archaebacterial L1e proteins are slightly shorter at both their N- and C-termini. A region of high sequence similarity (7 of 14 residues) occurs near the center of the proteins.

A family of genes encode the multiple forms of the Saccharomyces cerevisiae ribosomal proteins equivalent to the Escherichia coli L12 protein and a single form of the L10-equivalent ribosomal protein.

The budding yeast Saccharomyces cerevisiae contains a family of genes that encodes four different but related small acidic ribosomal proteins designated L12eIA, L12eIB, L12eIIA, and L12eIIB and a single larger protein designated L10e. These proteins are equivalent (e) to the L12 and L10 proteins of Escherichia coli that assemble as a 4:1 complex onto the large ribosomal subunit. The five yeast genes (or their cDNAs) have been cloned and sequenced (M. Remacha, M. T. Saenz-Robles, M. D. Vilella, and J. P. G. Ballesta, J. Biol. Chem. 263:9044-9101, 1988; K. Mitsui and K. Tsurugi, Nucleic Acids Res. 16:3573, 3574, and 3575, 1988; this work). Here, the transcripts of these genes were characterized and quantitated and the proteins they encode were compared and aligned. Four of the genes, L12eIA, -IB, -IIA, and L10e, are uninterrupted, whereas the L12eIIB gene contains a 301-nucleotide-long intron between codons 38 and 39. The transcripts derived from each of these genes were analyzed by Northern (RNA) hybridization, primer extension, and S1 nuclease protection. All five genes are expressed, albeit at different levels. The transcript levels are coordinate and exhibit growth rate-dependent regulation in rich (glucose) and poor (ethanol) media. The five yeast proteins each contain a highly conserved acidic carboxy terminus of about 20 residues in length. This domain of unknown function is also present in archaebacterial but absent from eubacterial L10e and L12e proteins. Comparisons of the factor-binding domains in the yeast and other eucaryotic and archaebacterial L12e proteins indicate that the original duplication to produce the type I and II genes was a very ancient event. The evolutionary relationships between the eucaryotic, archaebacterial, and eubacterial L10e and L12e genes (and proteins) are discussed.

Sequence alignment and evolutionary comparison of the L10 equivalent and L12 equivalent ribosomal proteins from archaebacteria, eubacteria, and eucaryotes.

The genes corresponding to the L10 and L12 equivalent ribosomal proteins (L10e and L12e) of Escherichia coli have been cloned and sequenced from two widely divergent species of archaebacteria, Halobacterium cutirubrum and Sulfolobus solfataricus. The deduced amino acid sequences of the L10e and L12e proteins have been compared to each other and to available eubacterial and eucaryotic sequences. We have identified the human P0 protein as the eucaryotic L10e. The L10e proteins from the three kingdoms were found to be colinear. The eubacterial L10e protein is much shorter than the archaebacterial-eucaryotic proteins because of two large deletions, one internal and one at the carboxy terminus. The archaebacterial and eucaryotic L12e proteins were also colinear; the eubacterial protein is homologous to the archaebacterial and eucaryotic L12e proteins, but has suffered rearrangement through what appear to be gene fusion events. Intraspecies comparisons between L10e and L12e sequences indicate the archaebacterial and eucaryotic L10e proteins contain a partial copy of the L12e protein fused to their carboxy terminus. In the eubacteria most of this fusion has been removed by the carboxy terminal deletion. Within the L12e-derived region, a 26-amino acid-long internal modular sequence reiterated thrice in the archaebacterial L10e, twice in the eucaryotic L10e, and once in the eubacterial L10e was discovered. This modular sequence also appears to be present as a single copy in all L12e proteins and may play a role in L12e dimerization, L10e-L12e complex formation, and the function of L10e-L12e complex in translation.(ABSTRACT TRUNCATED AT 250 WORDS)

Characterization of the L11, L1, L10 and L12 equivalent ribosomal protein gene cluster of the halophilic archaebacterium Halobacterium cutirubrum.

We have cloned and characterized a 5.2 kb fragment of genomic Halobacterium cutirubrum DNA encoding two potential proteins of unknown function (ORF and NAB) and four proteins which are equivalent to the L11, L1, L10 and L12 ribosomal proteins of Escherichia coli (L11e, L1e, L10e and L12e). The ribosomal protein genes are clustered in the same order as that in E. coli although the transcription pattern differs. Transcripts characterized include (i) abundant monocistronic L11e and tricistronic L1e-L10e-L12e transcripts; (ii) less abundant bicistronic NAB-L11e and monocistronic NAB transcripts and (iii) a very rare ORF monocistronic transcript. The consensus sequence in the promoter region is TTCGA ... 4-10 nucleotides ... TTAA ... 25-26 nucleotides ... initiation site; termination generally occurs on poly(T) tracts following GC-rich regions. Poly(T) tracts in the sense strands within coding regions are notably absent; this is probably related to their participation in transcription termination and to the fact that these ribosomal protein genes are highly expressed and stoichiometrically balanced. In the third position of the codons G or C is utilized 87% of the time. The 74 nt long untranslated leader of the L1e-L10e-L12e transcript contains a region that has a sequence and structure almost identical to a region within the binding domain for the L1e protein in 23S rRNA and highly similar to the E. coli L11-L1 mRNA leader sequence that has been implicated in autogenous translational regulation. Other transcripts are initiated at or adjacent to the ATG translation initiation codon.

Organization of genes encoding the L11, L1, L10, and L12 equivalent ribosomal proteins in eubacteria, archaebacteria, and eucaryotes.

Archaebacterial and eucaryotic cytoplasmic ribosomes contain proteins equivalent to the L11, L1, L10, and L12 proteins of the eubacterium Escherichia coli. In E. coli the genes encoding these ribosomal proteins are clustered, cotranscribed, and autogenously regulated at the level of mRNA translation. Genomic restriction fragments encoding the L11e, L1e, L10e, and L12e (equivalent) proteins from two divergent archaebacteria. Halobacterium cutirubrum and Sulfolobus solfataricus, and the L10e and L12e proteins from the eucaryote Saccharomyces cerevisiae have been cloned, sequenced, and analyzed. In the archaebacteria, as in eubacteria, the four genes are clustered and the L11e, L1e, L10e, and L12e order is maintained. The transcription pattern of the H. cutirubrum cluster is different from the E. coli pattern and the flanking genes on either side of the tetragenic clusters in E. coli, H. cutirubrum, and Sulfolobus solfataricus are all unrelated to each other. In the eucaryote Saccharomyces cerevisiae there is a single L10e gene and four separate L12e genes that are designated L12eIA, L12eIB, L12eIIA, and L12eIIB. These five genes are not closely linked and each is transcribed as a monocistronic mRNA; the L10e, L12eIA, L12eIB, and the L12eIIA genes are contiguous and uninterrupted, whereas the L12eIIB gene is interrupted by a 301 nucleotide long intron located between codons 38 and 39.

Structure and evolution of the L11, L1, L10, and L12 equivalent ribosomal proteins in eubacteria, archaebacteria, and eucaryotes.

The genes corresponding to the L11, L1, L10, and L12 equivalent ribosomal proteins (L11e, L1e, L10e, and L12e) of Escherichia coli have been cloned and sequenced from two widely divergent species of archaebacteria, Halobacterium cutirubrum and Sulfolobus solfataricus, and the L10 and four different L12 genes have been cloned and sequenced from the eucaryote Saccharomyces cerevisiae. Alignments between the deduced amino acid sequences of these proteins and to other available homologous proteins of eubacteria and eucaryotes have been made. The data suggest that the archaebacteria are a distinct coherent phylogenetic group. Alignment of the proline-rich L11e proteins reveals that the N-terminal region, believed to be responsible for interaction with release factor 1, is the most highly conserved region and that there is specific conservation of most of the proline residues, which may be important in maintaining the highly elongated structure of the molecule. Although L11 is the most highly methylated protein in the E. coli ribosome, the sites of methylation are not conserved in the archaebacterial L11e proteins. The L1e proteins of eubacteria and archaebacteria show two regions of very high similarity near the center and the carboxy termini of the proteins. The L10e proteins of all kingdoms are colinear and contain approximately three fourths of an L12e protein fused to their carboxy terminus, although much of this fusion has been lost in the truncated eubacterial protein. The archaebacterial and eucaryotic L12e proteins are colinear, whereas the eubacterial protein has suffered a rearrangement through what appear to be gene fusion events. Within the L12e derived region of the L10e proteins there exists a repeated module of 26 amino acids, present in two copies in eucaryotes, three in archaebacteria, and one in eubacteria. This modular sequence is apparently also present in the L12e proteins of all kingdoms and may play a role in L12e dimerization, L10e-L12e complex formation, and the function of the L10e-L12e complex in translation.

Temperature-sensitive beta-lactam-tolerant mutants of Escherichia coli.

Seven temperature-sensitive penicillin-tolerant mutants of Escherichia coli strain LD5 (thi lysA dapD) were isolated and characterized. Treatment with beta-lactams caused lysis of the mutants at 30 degrees C. Although growth of the mutants at 42 degrees C was inhibited by beta-lactams, no lysis occurred. The mutants were also slightly tolerant to D-cycloserine at 42 degrees C but lysed readily when deprived of diaminopimelate or when treated with moenomycin. The minimum inhibitory concentrations of various antibiotics were the same for the mutants and their parent. The mutations conferring penicillin tolerance were phenotypically suppressed in the presence of a variety of compounds which may act as chaotropic or antichaotropic agents. No defects in penicillin-binding proteins and peptidoglycan hydrolases were detected. Temperature-resistant revertants of the mutants were no longer tolerant to penicillin-induced autolysis at 42 degrees C. The mutations in five isolates were localized to the 56 to 61 min region of the E. coli linkage map and to the 44 to 51 min region in the case of two other isolates.