Signal-exon trap: a novel method for the identification of signal sequences from genomic DNA

We describe a genomic DNA-based signal sequence trap method, signal-exon trap (SET), for the identification of genes encoding secreted and membrane-bound proteins. SET is based on the coupling of an exon trap to the translation of captured exons, which allows screening of the exon-encoded polypeptides for signal peptide function. Since most signal sequences are expected to be located in the 5′-terminal exons of genes, we first demonstrate that trapping of these exons is feasible. To test the applicability of SET for the screening of complex genomic DNA, we evaluated two critical features of the method. Specificity was assessed by the analysis of random genomic DNA and efficiency was demonstrated by screening a 425 kb YAC known to contain the genes of four secretory or membrane-bound proteins. All trapped clones contained a translation initiation signal followed by a hydrophobic stretch of amino acids representing either a known signal peptide, transmembrane domain or novel sequence. Our results suggest that SET is a potentially useful method for the isolation of signal sequence-containing genes and may find application in the discovery of novel members of known secretory gene clusters, as well as in other positional cloning approaches.     

Analysis of progressive deletions of the transmembrane and cytoplasmic domains of influenza hemagglutinin

Site-directed oligonucleotide mutagenesis has been used to introduce chain termination codons into the cloned DNA sequences encoding the carboxy-terminal transmembrane (27 amino acids) and cytoplasmic (10 amino acids) domains of influenza virus hemagglutinin (HA). Four mutant genes were constructed which express truncated forms of HA that lack the cytoplasmic domain and terminate at amino acids 9, 14, 17, or 27 of the wild-type hydrophobic domain. Analysis of the biosynthesis and intracellular transport of these mutants shows that the cytoplasmic tail is not needed for the efficient transport of HA to the cell surface; the stop-transfer sequences are located in the hydrophobic domain; 17 hydrophobic amino acids are sufficient to anchor HA stably in the membrane; and mutant proteins with truncated hydrophobic domains show drastic alterations in transport, membrane association, and stability.     

Selection of Arabidopsis genes encoding secreted and plasma membrane proteins

Secreted and plasma membrane proteins play crucial roles in a variety of physiological and developmental processes of multicellular organisms. Systematic cloning of the genes encoding these proteins is therefore of general interest. An effective method of trapping signal sequences was first described by Tashiro et al. (1993), and a similar yet more efficient method was reported by Klein et al. (1996) and Jacobs et al. (1997). In this study, we carried out the latter yeast-based signal sequence trap to clone genes from Arabidopsis thaliana encoding secreted and plasma membrane proteins. Of 144 sequenced cDNA clones, 18% are identical to previously cloned Arabidopsis thaliana genes, 12% are homologous to genes identified from various organisms, and 46% are novel. All of the isolated genes identical or homologous to previously reported genes are either secreted or plasma membrane proteins, and the remaining novel genes appear to contain functional signal sequences based on computer-aided sequence analysis. The full-length cDNA clones of one homologous gene and another novel gene were isolated and sequenced. The deduced amino acid sequences suggest that the former encodes a secreted protein, and the latter encodes a type 1 membrane protein. These results indicate that the signal sequence trap method is effective and useful for the isolation of plant genes encoding secreted and plasma membrane proteins.     

On filtering false positive transmembrane protein predictions

While helical transmembrane (TM) region prediction tools achieve high (>90%) success rates for real integral membrane proteins, they produce a considerable number of false positive hits in sequences of known nontransmembrane queries. We propose a modification of the dense alignment surface (DAS) method that achieves a substantial decrease in the false positive error rate. Essentially, a sequence that includes possible transmembrane regions is compared in a second step with TM segments in a sequence library of documented transmembrane proteins. If the performance of the query sequence against the library of documented TM segment-containing sequences in this test is lower than an empirical threshold, it is classified as a non-transmembrane protein. The probability of false positive prediction for trusted TM region hits is expressed in terms of E-values. The modified DAS method, the DAS-TMfilter algorithm, has an unchanged high sensitivity for TM segments ( approximately 95% detected in a learning set of 128 documented transmembrane proteins). At the same time, the selectivity measured over a non-redundant set of 526 soluble proteins with known 3D structure is approximately 99%, mainly because a large number of falsely predicted single membrane-pass proteins are eliminated by the DAS-TMfilter algorithm.     

Identification of putative parasitism genes expressed in the esophageal gland cells of the soybean cyst nematode Heterodera glycines

Cloning parasitism genes encoding secretory proteins expressed in the esophageal gland cells is the key to understanding the molecular basis of nematode parasitism of plants. Suppression subtractive hybridization (SSH) with the microaspirated contents from Heterodera glycines esophageal gland cells and intestinal region was used to isolate genes expressed preferentially in the gland cells of parasitic stages. Twenty-three unique cDNA sequences from a SSH cDNA library were identified and hybridized to the genomic DNA of H. glycines in Southern blots. Full-length cDNAs of 21 clones were obtained by screening a gland-cell long-distance polymerase chain reaction cDNA library. Deduced proteins of 10 clones were preceded by a signal peptide for secretion, and PSORT II computer analysis predicted eight proteins as extracellular, one as nuclear, and one as plasmalemma localized. In situ hybridization showed that four of the predicted extracellular clones were expressed specifically in the dorsal gland cell, one in the subventral gland cells, and three in the intestine in H. glycines. The predicted nuclear clone and the plasmalemma-localized clone were expressed in the subventral gland cells and the dorsal gland cell, respectively. SSH is an efficient method for cloning putative parasitism genes encoding esophageal gland cell secretory proteins that may have a role in H. glycines parasitism of soybean.     

Internal gene duplication in the evolution of prokaryotic transmembrane proteins

We investigated the evolution of transmembrane (TM) topology by detecting partial sequence repeats in TM protein sequences and analyzing them in detail. A total of 377 sequences that seem to have evolved by internal gene duplication events were found among 38,124 predicted TM protein sequences (except for single-spannings) from 87 prokaryotic genomes. Various types of internal duplication patterns were identified in these sequences. The majority of them are diploid-type (including quasi-diploid-type) duplication in which a primordial protein sequence was duplicated internally to become an extant TM protein with twice as many TM segments as the primordial one, and the remaining ones are partial duplications including triploid-type. The diploid-type repeats are recognized in many 8-tms, 10-tms and 12-tms TM protein sequences, suggesting the diploid-type duplication was a principle mechanism in the evolutionary development of these types of TM proteins. The "positive-inside" rule is satisfied in whole sequences of both 10-tms and 8-tms TM proteins and in both halves of 10-tms proteins while not necessarily in the second half of 8-tms proteins, providing fit examples of "internal divergent topology evolution" likely occurred after a diploid-type internal duplication event. From analyzing the partial duplication patterns, several evolutionary pathways were recognized for 6-tms TM proteins, i.e. from primordial 2-tms, 3-tms and 4-tms TM proteins to extant 6-tms proteins. Similarly, the duplication pattern analysis revealed plausible evolution scenarios that 7-tms TM proteins have arisen from 3-tms, 4-tms and 5-tms TM protein precursors via partial internal gene duplications.     

Effect of sequence hydrophobicity and bilayer width upon the minimum length required for the formation of transmembrane helices in membranes

The minimum hydrophobic length necessary to form a transmembrane (TM) helix in membranes was investigated using model membrane-inserted hydrophobic helices. The fluorescence of a Trp at the center of the sequence and its sensitivity to quenching were used to ascertain helix position within the membrane. Peptides with hydrophobic cores composed of poly(Leu) were compared to sequences containing a poly 1:1 Leu:Ala core (which have a hydrophobicity typical of natural TM helices). Studies varying bilayer width revealed that the poly(Leu) core peptides predominately formed a TM state when the bilayer width exceeded hydrophobic sequence length by (i.e. when negative mismatch was) up to ∼ 11-12 Å (e.g. the case of a 11-12 residue hydrophobic sequence in bilayers with a biologically relevant width, i.e. dioleoylphosphatidylcholine (DOPC) bilayers), while poly(LeuAla) core peptides formed predominantly TM state with negative mismatch of up to 9 Å (a 13 residue hydrophobic sequence in DOPC bilayers). This indicates that minimum length necessary to form a predominating amount of a TM state (minimum TM length) is only modestly hydrophobicity-dependent for the sequences studied here, and a formula that defines the minimum TM length as a function of hydrophobicity for moderately-to-highly hydrophobic sequences was derived. The minimum length able to form a stable TM helix for alternating LeuAla sequences, and that for sequences with a Leu block followed by an Ala block, was similar, suggesting that a hydrophobicity gradient along the sequence may not be an important factor in TM stability. TM stability was also similar for sequences flanked by different charged ionizable residues (Lys, His, Asp). However, ionizable flanking residues destabilized the TM configuration much more when charged than when uncharged. The ability of short hydrophobic sequences to form TM helices in membranes in the presence of substantial negative mismatch implies that lipid bilayers have a considerable ability to adjust to negative mismatch, and that short TM helices may be more common than generally believed. Factors that modulate the ability of bilayers to adjust to mismatch may strongly affect the configuration of short hydrophobic helices.     

The specificity of association of the IgD molecule with the accessory proteins BAP31/BAP29 lies in the IgD transmembrane sequence.

Mature B cells co-express on their cell surface two classes of antigen receptor, the IgM and IgD immunoglobulins. The structural and functional differences between the two receptor classes are poorly understood. Recently two proteins of 29 and 31 kDa (BAP29 and BAP31) have been described that are preferentially associated with membrane IgD but only weakly with membrane IgM. We describe here the cloning of full-length murine and human BAP31 cDNAs encoding proteins of 245 and 246 amino acids respectively. The two BAP31 proteins are 95% identical. The BAP31 gene is ubiquitously expressed in murine tissues and is located on the X chromosome in both mouse and man. The murine BAP31 protein has 43% sequence identity to murine BAP29. Both proteins have a hydrophobic N-terminus and an alpha-helical C-terminus which ends with a KKXX motif implicated in vesicular transport. By a mutational analysis we have identified amino acids in the transmembrane sequence of the delta m chain that are critical for binding to BAP31/BAP29. A structural model of the BAPs and their potential functions are discussed.