Sequence dependent hypermutation of the immunoglobulin heavy chain in cultured B cells

The variable (V) regions of immunoglobulin heavy and light chains undergo high rates of somatic mutation during the immune response. Although point mutations accumulate throughout the V regions and their immediate flanking sequences, analysis of large numbers of mutations that have arisen in vivo reveal that the triplet AGC appears to be most susceptible to mutation. We have stably transfected B cell lines with γ2a heavy chain constructs containing TAG nonsense codons in their V regions that are part of either a putative (T)AGC hot spot or a (T)AGA non-hot spot motif. Using an ELISA spot assay to detect revertants and fluctuation analysis to determine rates of mutation, the rate of reversion of the TAG nonsense codon has been determined for different motifs in different parts of the V region. In the NSO plasma cell line, the (T)AGC hot spot motif mutates at rates of ≈6 × 10−4/bp per generation and ≈3 × 10−5/bp per generation at residues 38 and 94 in the V region. At each of these locations, the (T)AGC hot spot motif is 20–30 times more likely to undergo mutation than the (T)AGA non-hot spot motif. Moreover, the AGA non-hot spot motif mutates at as high a rate as the hot spot motif when it is located adjacent to hot spot motifs, suggesting that more extended sequences influence susceptibility to mutation.

Identification of a cytochrome b-type NAD(P)H oxidoreductase ubiquitously expressed in human cells

Cytochrome b-type NAD(P)H oxidoreductases are involved in many physiological processes, including iron uptake in yeast, the respiratory burst, and perhaps oxygen sensing in mammals. We have identified a cytosolic cytochrome b-type NAD(P)H oxidoreductase in mammals, a flavohemoprotein (b5+b5R) containing cytochrome b5 (b5) and b5 reductase (b5R) domains. A genetic approach, using blast searches against dbest for FAD-, NAD(P)H-binding sequences followed by reverse transcription–PCR, was used to clone the complete cDNA sequence of human b5+b5R from the hepatoma cell line Hep 3B. Compared with the classical single-domain b5 and b5R proteins localized on endoplasmic reticulum membrane, b5+b5R also has binding motifs for heme, FAD, and NAD(P)H prosthetic groups but no membrane anchor. The human b5+b5R transcript was expressed at similar levels in all tissues and cell lines that were tested. The two functional domains b5* and b5R* are linked by an approximately 100-aa-long hinge bearing no sequence homology to any known proteins. When human b5+b5R was expressed as c-myc adduct in COS-7 cells, confocal microscopy revealed a cytosolic localization at the perinuclear space. The recombinant b5+b5R protein can be reduced by NAD(P)H, generating spectrum typical of reduced cytochrome b with alpha, beta, and Soret peaks at 557, 527, and 425 nm, respectively. Human b5+b5R flavohemoprotein is a NAD(P)H oxidoreductase, demonstrated by superoxide production in the presence of air and excess NAD(P)H and by cytochrome c reduction in vitro. The properties of this protein make it a plausible candidate oxygen sensor.

Generation of longer cDNA fragments from serial analysis of gene expression tags for gene identification

We have developed a technique called the generation of longer cDNA fragments from serial analysis of gene expression (SAGE) tags for gene identification (GLGI), to convert SAGE tags of 10 bases into their corresponding 3′ cDNA fragments covering hundred bases. A primer containing the 10-base SAGE tag is used as the sense primer, and a single base anchored oligo(dT) primer is used as an antisense primer in PCR, together with Pfu DNA polymerase. By using this approach, a cDNA fragment extending from the SAGE tag toward the 3′ end of the corresponding sequence can be generated. Application of the GLGI technique can solve two critical issues in applying the SAGE technique: one is that a longer fragment corresponding to a SAGE tag, which has no match in databases, can be generated for further studies; the other is that the specific fragment corresponding to a SAGE tag can be identified from multiple sequences that match the same SAGE tag. The development of the GLGI method provides several potential applications. First, it provides a strategy for even wider application of the SAGE technique for quantitative analysis of global gene expression. Second, a combined application of SAGE/GLGI can be used to complete the catalogue of the expressed genes in human and in other eukaryotic species. Third, it can be used to identify the 3′ cDNA sequence from any exon within a gene. It can also be used to confirm the reality of exons predicted by bioinformatic tools in genomic sequences. Fourth, a combined application of SAGE/GLGI can be applied to define the 3′ boundary of expressed genes in the genomic sequences in human and in other eukaryotic genomes.

An ethylene-induced cDNA encoding a lipase expressed at the onset of senescence

A cDNA clone encoding a lipase (lipolytic acyl hydrolase) expressed at the onset of petal senescence has been isolated by screening a cDNA expression library prepared from carnation flowers (Dianthus caryophyllus). The cDNA contains the lipase consensus sequence, ITFAGHSLGA, and encodes a 447-amino acid polypeptide with a calculated molecular mass of 50.2 kDa that appears to be a cytosolic protein. Over-expression of the clone in Escherichia coli yielded a protein of the expected molecular weight that proved capable of deesterifying fatty acids from p-nitrophenylpalmitate, tri-linolein, soybean phospholipid, and Tween in both in vitro and in situ assays of enzyme activity. The abundance of the lipase mRNA increases just as carnation flowers begin to senesce, and expression of the gene is also induced by treatment with ethylene. Southern blot analyses of carnation genomic DNA have indicated that the lipase is a single copy gene. The lipase gene is also expressed in carnation leaves and is up-regulated when the leaves are treated with ethylene. Deesterification of membrane lipids and ensuing loss of membrane structural integrity are well established early events of plant senescence, and the expression pattern of this lipase gene together with the lipolytic activity of its cognate protein indicate that it plays a fundamentally central role in mediating the onset of senescence.

NEEDLY, a Pinus radiata ortholog of FLORICAULA/LEAFY genes, expressed in both reproductive and vegetative meristems

The LEAFY/FLORICAULA genes from Arabidopsis and Antirrhinum are necessary for normal flower development and play a key role in diverse angiosperm species. A homologue of these flower meristem-identity genes, NEEDLY (NLY), has been identified in Pinus radiata. Although the NLY protein shares extensive sequence similarity with its angiosperm counterparts, it is lacking the proline-rich and acidic motifs thought to function as transcriptional activation domains. NLY already is expressed during vegetative development at least 5 years before the transition to the reproductive phase. Expression of NLY in transgenic Arabidopsis promotes floral fate, demonstrating that, despite its sequence divergence, NLY encodes a functional ortholog of the FLORICAULA/LEAFY genes of angiosperms. Expression of the LFY∷NLY transgene can largely complement the defects in flower development caused by a severe lfy allele.

Diversity and phylogeny of gephyrin: Tissue-specific splice variants, gene structure, and sequence similarities to molybdenum cofactor-synthesizing and cytoskeleton-associated proteins

Gephyrin is essential for both the postsynaptic localization of inhibitory neurotransmitter receptors in the central nervous system and the biosynthesis of the molybdenum cofactor (Moco) in different peripheral organs. Several alternatively spliced gephyrin transcripts have been identified in rat brain that differ in their 5′ coding regions. Here, we describe gephyrin splice variants that are differentially expressed in non-neuronal tissues and different regions of the adult mouse brain. Analysis of the murine gephyrin gene indicates a highly mosaic organization, with eight of its 29 exons corresponding to the alternatively spliced regions identified by cDNA sequencing. The N- and C-terminal domains of gephyrin encoded by exons 3–7 and 16–29, respectively, display sequence similarities to bacterial, invertebrate, and plant proteins involved in Moco biosynthesis, whereas the central exons 8, 13, and 14 encode motifs that may mediate oligomerization and tubulin binding. Our data are consistent with gephyrin having evolved from a Moco biosynthetic protein by insertion of protein interaction sequences.

Identification of two novel transmembrane γ-carboxyglutamic acid proteins expressed broadly in fetal and adult tissues

The proline-rich γ-carboxyglutamic acid (Gla) proteins (PRGPs) 1 and 2 are the founding members of a family of vitamin K-dependent single-pass integral membrane proteins characterized by an extracellular amino terminal domain of approximately 45 amino acids that is rich in Gla. The intracellular carboxyl terminal region of these two proteins contains one or two copies of the sequence PPXY, a motif present in a variety of proteins involved in such diverse cellular functions as signal transduction, cell cycle progression, and protein turnover. In this report, we describe the cloning of the cDNAs for two additional human transmembrane Gla proteins (TMG) of 20–24 kDa named TMG3 and TMG4. These two proteins possess extracellular Gla domains with 13 or 9 potential Gla residues, respectively, followed by membrane-spanning hydrophobic regions and cytoplasmic carboxyl terminal regions that contain PPXY motifs. This emerging family of integral membrane Gla proteins includes proline-rich Gla protein (PRGP) 1, PRGP2, TMG3, and TMG4, all of which are characterized by broad and variable distribution in both fetal and adult tissues. Members of this family can be grouped into two subclasses on the basis of their gene organization and amino acid sequence. These observations suggest novel physiological functions for vitamin K beyond its known role in the biosynthesis of proteins involved in blood coagulation and bone development. The identification and characterization of these proteins may allow a more complete understanding of the teratogenic consequences of exposure in utero to vitamin K antagonists, such as warfarin-based anticoagulants.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

A shark antibody heavy chain encoded by a nonsomatically rearranged VDJ is preferentially expressed in early development and is convergent with mammalian IgG

In most vertebrate embryos and neonates studied to date unique antigen receptors (antibodies and T cell receptors) are expressed that possess a limited immune repertoire. We have isolated a subclass of IgM, IgM1gj, from the nurse shark Ginglymostoma cirratum that is preferentially expressed in neonates. The variable (V) region gene encoding the heavy (H) chain underwent V-D-J rearrangement in germ cells (“germline-joined”). Such H chain V genes were discovered over 10 years ago in sharks but until now were not shown to be expressed at appreciable levels; we find expression of H1gj in primary and secondary lymphoid tissues early in life, but in adults only in primary lymphoid tissue, which is identified in this work as the epigonal organ. H1gj chain associates covalently with light (L) chains and is most similar in sequence to IgM H chains, but like mammalian IgG has three rather than the four IgM constant domains; deletion of the ancestral IgM C2 domain thus defines both IgG and IgM1gj. Because sharks are the members of the oldest vertebrate class known to possess antibodies, unique or specialized antibodies expressed early in ontogeny in sharks and other vertebrates were likely present at the inception of the adaptive immune system.

Characterization of Two Divergent Endo-β-1,4-Glucanase cDNA Clones Highly Expressed in the Nonclimacteric Strawberry Fruit

Two cDNAs clones (Cel1 and Cel2) encoding divergent endo-β-1,4-glucanases (EGases) have been isolated from a cDNA library obtained from ripe strawberry (Fragaria x ananassa Duch) fruit. The analysis of the amino acid sequence suggests that Cel1 and Cel2 EGases have different secondary and tertiary structures and that they differ in the presence of potential N-glycosylation sites. By in vitro translation we show that Cel1 and Cel2 bear a functional signal peptide, the cleavage of which yields mature proteins of 52 and 60 kD, respectively. Phylogenetic analysis revealed that the Cel2 EGase diverged early in evolution from other plant EGases. Northern analysis showed that both EGases are highly expressed in fruit and that they have different temporal patterns of accumulation. The Cel2 EGase was expressed in green fruit, accumulating as the fruit turned from green to white and remaining at an elevated, constant level throughout fruit ripening. In contrast, the Cel1 transcript was not detected in green fruit and only a low level of expression was observed in white fruit. The level of Cel1 mRNA increased gradually during ripening, reaching a maximum in fully ripe fruit. The high levels of Cel1 and Cel2 mRNA in ripe fruit and their overlapping patterns of expression suggest that these EGases play an important role in softening during ripening. In addition, the early expression of Cel2 in green fruit, well before significant softening begins, suggests that the product of this gene may also be involved in processes other than fruit softening, e.g. cell wall expansion.

Purification of the Trehalase GMTRE1 from Soybean Nodules and Cloning of Its cDNA. GMTRE1 Is Expressed at a Low Level in Multiple Tissues1

Trehalose (α-d-glucopyranosyl-1,1-α-d-glucopyranoside), a disaccharide widespread among microbes and lower invertebrates, is generally believed to be nonexistent in higher plants. However, the recent discovery of Arabidopsis genes whose products are involved in trehalose synthesis has renewed interest in the possibility of a function of trehalose in higher plants. We previously showed that trehalase, the enzyme that degrades trehalose, is present in nodules of soybean (Glycine max [L.] Merr.), and we characterized the enzyme as an apoplastic glycoprotein. Here we describe the purification of this trehalase to homogeneity and the cloning of a full-length cDNA encoding this enzyme, named GMTRE1 (G. max trehalase 1). The amino acid sequence derived from the open reading frame of GMTRE1 shows strong homology to known trehalases from bacteria, fungi, and animals. GMTRE1 is a single-copy gene and is expressed at a low but constant level in many tissues.

Sequence-Specific Interaction between the Disintegrin Domain of Mouse ADAM 3 and Murine Eggs: Role of β1 Integrin-associated Proteins CD9, CD81, and CD98

ADAM 3 is a sperm surface glycoprotein that has been implicated in sperm-egg adhesion. Because little is known about the adhesive activity of ADAMs, we investigated the interaction of ADAM 3 disintegrin domains, made in bacteria and in insect cells, with murine eggs. Both recombinant proteins inhibited sperm-egg binding and fusion with potencies similar to that which we recently reported for the ADAM 2 disintegrin domain. Alanine scanning mutagenesis revealed a critical importance for the glutamine at position 7 of the disintegrin loop. Fluorescent beads coated with the ADAM 3 disintegrin domain bound to the egg surface. Bead binding was inhibited by an authentic, but not by a scrambled, peptide analog of the disintegrin loop. Bead binding was also inhibited by the function-blocking anti-α6 monoclonal antibody (mAb) GoH3, but not by a nonfunction blocking anti-α6 mAb, or by mAbs against either the αv or β3 integrin subunits. We also present evidence that in addition to the tetraspanin CD9, two other β1-integrin-associated proteins, the tetraspanin CD81 as well as the single pass transmembrane protein CD98 are expressed on murine eggs. Antibodies to CD9 and CD98 inhibited in vitro fertilization and binding of the ADAM 3 disintegrin domain. Our findings are discussed in terms of the involvement of multiple sperm ADAMs and multiple egg β1 integrin-associated proteins in sperm-egg binding and fusion. We propose that an egg surface “tetraspan web” facilitates fertilization and that it may do so by fostering ADAM–integrin interactions.

Excitation–contraction coupling is unaffected by drastic alteration of the sequence surrounding residues L720–L764 of the α1S II-III loop

The II-III loop of the skeletal muscle dihydropyridine receptor (DHPR) α1S subunit is responsible for bidirectional-signaling interactions with the ryanodine receptor (RyR1): transmitting an orthograde, excitation–contraction (EC) coupling signal to RyR1 and receiving a retrograde, current-enhancing signal from RyR1. Previously, several reports argued for the importance of two distinct regions of the skeletal II-III loop (residues R681–L690 and residues L720–Q765, respectively), claiming for each a key function in DHPR–RyR1 communication. To address whether residues 720–765 of the II-III loop are sufficient to enable skeletal-type (Ca2+ entry-independent) EC coupling and retrograde interaction with RyR1, we constructed a green fluorescent protein (GFP)-tagged chimera (GFP-SkLM) having rabbit skeletal (Sk) DHPR sequence except for a II-III loop (L) from the DHPR of the house fly, Musca domestica (M). The Musca II-III loop (75% dissimilarity to α1S) has no similarity to α1S in the regions R681–L690 and L720–Q765. GFP-SkLM expressed in dysgenic myotubes (which lack endogenous α1S subunits) was unable to restore EC coupling and displayed strongly reduced Ca2+ current densities despite normal surface expression levels and correct triad targeting (colocalization with RyR1). Introducing rabbit α1S residues L720–L764 into the Musca II-III loop of GFP-SkLM (substitution for Musca DHPR residues E724–T755) completely restored bidirectional coupling, indicating its dependence on α1S loop residues 720–764 but its independence from other regions of the loop. Thus, 45 α1S-residues embedded in a very dissimilar background are sufficient to restore bidirectional coupling, indicating that these residues may be a site of a protein–protein interaction required for bidirectional coupling.

A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

The Two Genes Encoding Starch-Branching Enzymes IIa and IIb Are Differentially Expressed in Barley1

The sbeIIa and sbeIIb genes, encoding starch-branching enzyme (SBE) IIa and SBEIIb in barley (Hordeum vulgare L.), have been isolated. The 5′ portions of the two genes are strongly divergent, primarily due to the 2064-nucleotide-long intron 2 in sbeIIb. The sequence of this intron shows that it contains a retro-transposon-like element. Expression of sbeIIb but not sbeIIa was found to be endosperm specific. The temporal expression patterns for sbeIIa and sbeIIb were similar and peaked around 12 d after pollination. DNA gel-blot analysis demonstrated that sbeIIa and sbeIIb are both single-copy genes in the barley genome. By fluorescence in situ hybridization, the sbeIIa and sbeIIb genes were mapped to chromosomes 2 and 5, respectively. The cDNA clones for SBEIIa and SBEIIb were isolated and sequenced. The amino acid sequences of SBEIIa and SBEIIb were almost 80% identical. The major structural difference between the two enzymes was the presence of a 94-amino acid N-terminal extension in the SBEIIb precursor. The (β/α)8-barrel topology of the α-amylase superfamily and the catalytic residues implicated in branching enzymes are conserved in both barley enzymes.