Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

Detection of mitochondrial single nucleotide polymorphisms using a primer elongation reaction on oligonucleotide microarrays

We have developed a novel allele-specific primer elongation protocol using a DNA polymerase on oligonucleotide chips. Oligonucleotide primers carrying polymorphic sites at their free 3́end were covalently bound to glass slides. The generation of single-stranded targets of genomic DNA containing single nuclotide polymorphisms (SNPs) to be typed was achieved by an asymmetric PCR reaction or exonuclease treatment of phosphothioate (PTO)-modified PCR products. In the presence of DNA polymerase and all four dNTPs, with Cy3-dUTP replacing dTTP, allele-specific extension of the immobilized primers took place along a stretch of target DNA sequence. The yield of elongated products was increased by repeated reaction cycles. We performed multiplexed assays with many small DNA targets, or used single targets of up to 4.4 kb mitochondrial DNA (mtDNA) sequence to detect multiple SNPs in one reaction. The latter approach greatly simplifies preamplification of SNP-containing regions, thereby providing a framework for typing hundreds of mtDNA polymorphisms.

T cell activation up-regulates cyclic nucleotide phosphodiesterases 8A1 and 7A3

Agents that increase intracellular cAMP inhibit the activation and function of T cells and can lead to cell death. Recently, it has been postulated that cAMP inhibits T cell function in large part by acting as a brake on the T cell receptor and costimulatory receptor pathways. Therefore, for full activation of the T cell to occur, this inhibitory influence must be removed. One likely mechanism for accomplishing this is by up-regulation and/or activation of specific cyclic nucleotide phosphodiesterases (PDEs), and such a mechanism for one phosphodiesterase, PDE7A1, has been reported. In this paper, we extend this mechanism to another isozyme variant of the same PDE family, PDE7A3. We also report the full-length sequence of human PDE8A1 and show that it also is induced in response to a combination of T cell receptor and costimulatory receptor pathway activation. However, the time course for induction of PDE8A1 is slower than that of PDE7A1. The basal level measured and, therefore, the apparent fold induction of PDE7A1 mRNA and protein depend in large part on the method of isolation of the T cells. On the other hand, regardless of the isolation method, the basal levels of PDE7A3 and PDE8A1 are very low and fold activation is much higher. Constitutively expressed PDE8A1 and PDE7A3 also have been isolated from a human T cell line, Hut78.

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.

The complete genome of the crenarchaeon Sulfolobus solfataricus P2

The genome of the crenarchaeon Sulfolobus solfataricus P2 contains 2,992,245 bp on a single chromosome and encodes 2,977 proteins and many RNAs. One-third of the encoded proteins have no detectable homologs in other sequenced genomes. Moreover, 40% appear to be archaeal-specific, and only 12% and 2.3% are shared exclusively with bacteria and eukarya, respectively. The genome shows a high level of plasticity with 200 diverse insertion sequence elements, many putative nonautonomous mobile elements, and evidence of integrase-mediated insertion events. There are also long clusters of regularly spaced tandem repeats. Different transfer systems are used for the uptake of inorganic and organic solutes, and a wealth of intracellular and extracellular proteases, sugar, and sulfur metabolizing enzymes are encoded, as well as enzymes of the central metabolic pathways and motility proteins. The major metabolic electron carrier is not NADH as in bacteria and eukarya but probably ferredoxin. The essential components required for DNA replication, DNA repair and recombination, the cell cycle, transcriptional initiation and translation, but not DNA folding, show a strong eukaryal character with many archaeal-specific features. The results illustrate major differences between crenarchaea and euryarchaea, especially for their DNA replication mechanism and cell cycle processes and their translational apparatus.

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.

Isolation of a superfamily of candidate disease-resistance genes in soybean based on a conserved nucleotide-binding site.

The tobacco N and Arabidopsis RPS2 genes, among several recently cloned disease-resistance genes, share highly conserved structure, a nucleotide-binding site (NBS). Using degenerate oligonucleotide primers for the NBS region of N and RPS2, we have amplified and cloned the NBS sequences from soybean. Each of these PCR-derived NBS clones detected low-or moderate-copy soybean DNA sequences and belongs to 1 of 11 different classes. Sequence analysis showed that all PCR clones encode three motifs (P-loop, kinase-2, and kinase-3a) of NBS nearly identical to those in N and RPS2. The intervening region between P-loop and kinase-3a of the 11 classes has high (26% average) amino acid sequence similarity to the N gene although not as high (19% average) to RPS2. These 11 classes represent a superfamily of NBS-containing soybean genes that are homologous to N and RPS2. Each class or subfamily was assessed for its positional association with known soybean disease-resistance genes through near-isogenic line assays, followed by linkage analysis in F2 populations using restriction fragment length polymorphisms. Five of the 11 subfamilies have thus far been mapped to the vicinity of known soybean genes for resistance to potyviruses (Rsv1 and Rpv), Phytophthora root rot (Rps1, Rps2, and Rps3), and powdery mildew (rmd). The conserved N- or RPS2-homologous NBS sequences and their positional associations with mapped soybean-resistance genes suggest that a number of the soybean disease-resistance genes may belong to this superfamily. The candidate subfamilies of NBS-containing genes identified by genetic mapping should greatly facilitate the molecular cloning of disease-resistance genes.

The primary structures of the Archaeon Halobacterium salinarium blue light receptor sensory rhodopsin II and its transducer, a methyl-accepting protein.

Recently, a large family of transducer proteins in the Archaeon Halobacterium salinarium was identified. On the basis of the comparison of the predicted structural domains of these transducers, three distinct subfamilies of transducers were proposed. Here we report isolation, complete gene sequences, and analysis of the encoded primary structures of transducer gene htrII, a member of family B, and its blue light receptor gene (sopII) of sensory rhodopsin II (SRII). The start codon ATG of the 714-bp sopII gene is one nucleotide beyond the termination codon TGA of the 2298-bp htrII gene. The deduced protein sequence of HtrII predicts a eubacterial chemotaxis transducer type with two hydrophobic membrane-spanning segments connecting sizable domains in the periplasm and cytoplasm. HtrII has a common feature with HtrI, the sensory rhodopsin I transducer; like HtrI, HtrII possesses a hydrophilic loop structure just after the second transmembrane segment. The C-terminal 299 residues (765 amino acid residues total) of HtrII show strong homology to the signaling and methylation domain of eubacterial transducer Tsr. The hydropathy plot of the primary structure of SRII indicates seven membrane-spanning alpha-helical segments, a characteristic feature of retinylidene proteins ("rhodopsins") from a widespread family of photoactive pigments. SRII shows high identity with SRI (42%), bacteriorhodopsin (BR) (32%), and halorhodopsin (24%). The crucial positions for retinal binding sites in these proteins are nearly identical, with the exception of Met-118 (numbering according to the mature BR sequence), which is replaced by Val in SRII. In BR, residues Asp-85 and Asp-96 are crucial in proton pumping. In SRII, the position corresponding to Asp-85 in BR is conserved, but the corresponding position of Asp-96 is replaced by an aromatic Tyr. Coexpression of the htrII and sopII genes restores SRII phototaxis to a mutant (Pho81) that contains a deletion in the htrI/sopI and insertion in htrII/sopII regions. This paper describes the first example that both HtrI and HtrII exist in the same halobacterial cell, confirming that different sensory rhodopsins SRI and SRII in the same organism have their own distinct transducers.

Interaction of the 82-kDa subunit of the vaccinia virus early transcription factor heterodimer with the promoter core sequence directs downstream DNA binding of the 70-kDa subunit.

The vaccinia virus early transcription factor (VETF), a heterodimeric protein composed of 82- and 70-kDa subunits, interacts with viral early promoters at both a sequence-specific core region upstream and a sequence-independent region downstream of the RNA start site. To determine the VETF subunit-promoter interactions, 32P-labeled DNA targets were chemically synthesized with uniquely positioned phosphorothioates to which azidophenacyl bromide moieties were coupled. After incubating the derivatized promoter with VETF and exposing the complex to 302-nm light, the protein was denatured and the individual subunits with or without covalently bound DNA were isolated with specific antiserum and analyzed by SDS/polyacrylamide gel electrophoresis. Using a set of 26 duplex probes, with uniquely positioned aryl azide moieties on the coding or template strands, we found that the 82-kDa subunit interacted primarily with the core region of the promoter, whereas the 70-kDa subunit interacted with the downstream region. Nucleotide substitutions in the core region that downregulate transcription affected the binding of both subunits: the 82-kDa subunit no longer exhibited specificity for upstream regions of the promoter but also bound to downstream regions, whereas the binding of the 70-kDa subunit was abolished even though the mutations were far upstream of its binding site. These results suggested mechanisms by which the interaction of the 82-kDa subunit with the core sequence directs binding of the 70-kDa subunit to DNA downstream.

On signal sequence polymorphisms and diseases of distribution.

We report a previously unappreciated property of the signals that target organelle-specific proteins to their subcellular sites of action. Such targeting sequences are shown to be polymorphic. We discovered this polymorphism when we cloned the mitochondrial manganese-containing superoxide dismutase from cell lines of normal individuals and patients with genetic diseases of premature aging and compared their sequences to each other and to those previously reported. The polymorphism consists of a single nucleotide change in the region of the DNA that encodes the signal sequence such that either an alanine or valine is present. Subsequently, eight cell lines were analyzed and all three possible combinations of the two signal sequences were observed. Such signal sequence polymorphisms could result in diseases of distribution, where essential proteins are not properly targeted, thereby leading to absolute or relative deficiencies of critical enzymes within specific cellular compartments. Progeria and related syndromes may be diseases of distribution.

Unexpected beta2-microglobulin sequence diversity in individual rainbow trout.

For mammals beta2-microglobulin (beta2m), the light chain of major histocompatibility complex (MHC) class I molecules, is invariant (or highly conserved) and is encoded by a single gene unlinked to the MHC. We find that beta2m of a salmonid fish, the rainbow trout (Oncorhynchus mykiss), does not conform to the mammalian paradigm. Ten of 12 randomly selected beta2m cDNA clones from an individual fish have different nucleotide sequences. A complex restriction fragment length polymorphism pattern is observed with rainbow trout, suggesting multiple beta2m genes in the genome, in excess of the two genes expected from the ancestral salmonid tetraploidy. Additional duplication and diversification of the beta2m genes might have occurred subsequently. Variation in the beta2m cDNA sequences is mainly at sites that do not perturb the structure of the mature beta2m protein, showing that the observed diversity of the trout beta2m genes is not primarily a result of pathogen selection.

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.

Identification, sequence, and expression of caveolin-2 defines a caveolin gene family.

Caveolin, a 21- to 24-kDa integral membrane protein, is a principal component of caveolae membranes. Caveolin interacts directly with heterotrimeric guanine nucleotide binding proteins (G proteins) and can functionally regulate their activity. Here, an approximately 20-kDa caveolin-related protein, caveolin-2, was identified through microsequencing of adipocyte-derived caveolin-enriched membranes; caveolin was retermed caveolin-1. Caveolins 1 and 2 are similar in most respects. mRNAs for both caveolin-1 and caveolin-2 are most abundantly expressed in white adipose tissue and are induced during adipocyte differentiation. Caveolin-2 colocalizes with caveolin-1, indicating that caveolin-2 also localizes to caveolae. However, caveolin-1 and caveolin-2 differ in their functional interactions with heterotrimeric G proteins, possibly explaining why caveolin-1 and -2 are coexpressed within a single cell.

Proposed active site domain in estrogen sulfotransferase as determined by mutational analysis.

Point mutations were selectively introduced into a cDNA for guinea pig estrogen sulfotransferase (gpEST); each construct was then expressed in Chinese hamster ovary K1 cells. The molecular site chosen for study is a conserved GXXGXXK sequence that resembles the P-loop-type nucleotide-binding motif for ATP- and GTP-binding proteins and is located near the C terminus of all steroid and phenol(aryl) sulfotransferases for which the primary structures are known. Preliminary experiments demonstrated that the GXXGXXK motif is essential for binding the activated sulfonate donor 3'-phosphoadenosine 5'-phosphosulfate (PAPS). The present study was undertaken to ascertain the relative importance of each individual residue of the motif. While the mutation of a single motif residue had little effect on the interaction between gpEST and PAPS as determined by kinetic analysis and photoaffinity labeling, the mutation of any two residues in concert resulted in an approximate 10-fold increase in the Km for PAPS and reduced photoaffinity labeling. The mutation of all three motif residues resulted in an inactive enzyme and complete loss of photoaffinity labeling. Interestingly, several mutants also displayed a striking effect on the Km for the steroid substrate; double mutants, again, demonstrated greater perturbations (8- to 28-fold increase) than did single mutants. Unexpectedly, whereas the mutation of nonmotif residues had a negligible effect on the Km for PAPS, a marked increase in the Km for the estrogen substrate ( > 30-fold) was noted. On the basis of these findings, it is concluded that the sequence GISGDWKN within the C-terminal domain of gpEST represents a critical component of the active site.

Subunit interactions in coordination of Ni2+ in cyclic nucleotide-gated channels.

Cyclic nucleotide-gated (CNG) channels present a unique model for studying the molecular mechanisms of channel gating. We have studied the mechanism of potentiation of expressed rod CNG channels by Ni2+ as a first step toward understanding the channel gating process. Here we report that coordination of Ni2+ between histidine residues (H420) on adjacent channel subunits occurs when the channels are open. Mutation of H420 to lysine completely eliminated the potentiation by Ni2+ but did not markedly alter the apparent cGMP affinity of the channel, indicating that the introduction of positive charge at the Ni(2+)-binding site was not sufficient to produce potentiation. Deletion or mutation of most of the other histidines present in the channel did not diminish potentiation by Ni2+. We studied the role of subunit interactions in Ni2+ potentiation by generating heteromultimeric channels using tandem dimers of the rod CNG channel sequence. Injection of single heterodimers in which one subunit contained H420 and the other did not (wt/H420Q or H420Q/wt) resulted in channels that were not potentiated by Ni2+. However, coinjection of both heterodimers into Xenopus oocytes resulted in channels that exhibited potentiation. The H420 residues probably occurred predominantly in nonadjacent subunits when each heterodimer was injected individually, but, when the two heterodimers were coinjected, the H420 residues could occur in adjacent subunits as well. These results suggest that the mechanism of Ni2+ potentiation involves intersubunit coordination of Ni2+ by H420. Based on the preferential binding of Ni2+ to open channels, we suggest that alignment of H420 residues of neighboring subunits into the Ni(2+)-coordinating position may be associated with channel opening.