Comparative studies on mammalian Hoxc8 early enhancer sequence reveal a baleen whale-specific deletion of a cis-acting element

Variations in regulatory regions of developmental control genes have been implicated in the divergence of axial morphologies. To find potentially significant changes in cis-regulatory regions, we compared nucleotide sequences and activities of mammalian Hoxc8 early enhancers. The nucleotide sequence of the early enhancer region is extremely conserved among mammalian clades, with five previously described cis-acting elements, A–E, being invariant. However, a 4-bp deletion within element C of the Hoxc8 early enhancer sequence is observed in baleen whales. When assayed in transgenic mouse embryos, a baleen whale enhancer (unlike other mammalian enhancers) directs expression of the reporter gene to more posterior regions of the neural tube but fails to direct expression to posterior mesoderm. We suggest that regulation of Hoxc8 in baleen whales differs from other mammalian species and may be associated with variation in axial morphology.

Interactions of the chaperone Hsp104 with yeast Sup35 and mammalian PrP

[PSI+] is a genetic element in yeast for which a heritable change in phenotype appears to be caused by a heritable change in the conformational state of the Sup35 protein. The inheritance of [PSI+] and the physical state of Sup35 in vivo depend on the protein chaperone Hsp104 (heat shock protein 104). Although these observations provide a strong genetic argument in support of the “protein-only” or “prion” hypothesis for [PSI+], there is, as yet, no direct evidence of an interaction between the two proteins. We report that when purified Sup35 and Hsp104 are mixed, the circular dichroism (CD) spectrum differs from that predicted by the addition of the proteins’ individual spectra, and the ATPase activity of Hsp104 is inhibited. Similar results are obtained with two other amyloidogenic substrates, mammalian PrP and β-amyloid 1-42 peptide, but not with several control proteins. With a group of peptides that span the PrP protein sequence, those that produced the largest changes in CD spectra also caused the strongest inhibition of ATPase activity in Hsp104. Our observations suggest that (i) previously described genetic interactions between Hsp104 and [PSI+] are caused by direct interaction between Hsp104 and Sup35; (ii) Sup35 and PrP, the determinants of the yeast and mammalian prions, respectively, share structural features that lead to a specific interaction with Hsp104; and (iii) these interactions couple a change in structure to the ATPase activity of Hsp104.

β subunits influence the biophysical and pharmacological differences between P- and Q-type calcium currents expressed in a mammalian cell line

Human epithelial kidney cells (HEK) were prepared to coexpress α1A, α2δ with different β calcium channel subunits and green fluorescence protein. To compare the calcium currents observed in these cells with the native neuronal currents, electrophysiological and pharmacological tools were used conjointly. Whole-cell current recordings of human epithelial kidney α1A-transfected cells showed small inactivating currents in 80 mM Ba2+ that were relatively insensitive to calcium blockers. Coexpression of α1A, βIb, and α2δ produced a robust inactivating current detected in 10 mM Ba2+, reversibly blockable with low concentration of ω-agatoxin IVA (ω-Aga IVA) or synthetic funnel-web spider toxin (sFTX). Barium currents were also supported by α1A, β2a, α2δ subunits, which demonstrated the slowest inactivation and were relatively insensitive to ω-Aga IVA and sFTX. Coexpression of β3 with the same combination as above produced inactivating currents also insensitive to low concentration of ω-Aga IVA and sFTX. These data indicate that the combination α1A, βIb, α2δ best resembles P-type channels given the rate of inactivation and the high sensitivity to ω-Aga IVA and sFTX. More importantly, the specificity of the channel blocker is highly influenced by the β subunit associated with the α1A subunit.

The E2F6 transcription factor is a component of the mammalian Bmi1-containing polycomb complex

The E2F transcription factors play a key role in the regulation of cellular proliferation and terminal differentiation. E2F6 is the most recently identified and the least well understood member of the E2F family. It is only distantly related to the other E2Fs and lacks the sequences responsible for both transactivation and binding to the retinoblastoma protein. Consistent with this finding, E2F6 can behave as a dominant negative inhibitor of the other E2F family members. In this study, we continue to investigate the possible role(s) of E2F6 in vivo. We report the isolation of RYBP, a recently identified member of the mammalian polycomb complex, as an E2F6-interacting protein. Mapping studies indicate that RYBP binds within the known “repression domain” of E2F6. Moreover, we demonstrate that endogenous E2F6 and polycomb group proteins, including RYBP, Ring1, MEL-18, mph1, and the oncoprotein Bmi1, associate with one another. These findings suggest that the biological properties of E2F6 are mediated through its ability to recruit the polycomb transcriptional repressor complex.

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.

Identification of sample-specific sequences in mammalian cDNA and genomic DNA by the novel ligation-mediated subtraction (Limes)

The representational difference analysis (RDA) and other subtraction techniques are used to enrich sample-specific sequences by elimination of ubiquitous sequences existing in both the sample of interest (tester) and the subtraction partner (driver). While applying the RDA to genomic DNA of cutaneous lymphoma cells in order to identify tumor relevant alterations, we predominantly isolated repetitive sequences and artificial repeat-mediated fusion products of otherwise independent PCR fragments (PCR hybrids). Since these products severely interfered with the isolation of tester-specific fragments, we developed a considerably more robust and efficient approach, termed ligation-mediated subtraction (Limes). In first applications of Limes, genomic sequences and/or transcripts of genes involved in the regulation of transcription, such as transforming growth factor β stimulated clone 22 related gene (TSC-22R), cell death and cytokine production (caspase-1) or antigen presentation (HLA class II sequences), were found to be completely absent in a cutaneous lymphoma line. On the assumption that mutations in tumor-relevant genes can affect their transcription pattern, a protocol was developed and successfully applied that allows the identification of such sequences. Due to these results, Limes may substitute/supplement other subtraction/comparison techniques such as RDA or DNA microarray techniques in a variety of different research fields.

Breaksite batch mapping, a rapid method for assay and identification of DNA breaksites in mammalian cells

DNA breaks occur during many processes in mammalian cells, including recombination, repair, mutagenesis and apoptosis. Here we report a simple and rapid method for assaying DNA breaks and identifying DNA breaksites. Breaksites are first tagged and amplified by ligation-mediated PCR (LM-PCR), using nested PCR primers to increase the specificity and sensitivity of amplification. Breaksites are then mapped by batch sequencing LM-PCR products. This allows easy identification of multiple breaksites per reaction without tedious fractionation of PCR products by gel electrophoresis or cloning. Breaksite batch mapping requires little starting material and can be used to identify either single- or double-strand breaks.

SpliceDB: database of canonical and non-canonical mammalian splice sites

A database (SpliceDB) of known mammalian splice site sequences has been developed. We extracted 43 337 splice pairs from mammalian divisions of the gene-centered Infogene database, including sites from incomplete or alternatively spliced genes. Known EST sequences supported 22 815 of them. After discarding sequences with putative errors and ambiguous location of splice junctions the verified dataset includes 22 489 entries. Of these, 98.71% contain canonical GT–AG junctions (22 199 entries) and 0.56% have non-canonical GC–AG splice site pairs. The remainder (0.73%) occurs in a lot of small groups (with a maximum size of 0.05%). We especially studied non-canonical splice sites, which comprise 3.73% of GenBank annotated splice pairs. EST alignments allowed us to verify only the exonic part of splice sites. To check the conservative dinucleotides we compared sequences of human non-canonical splice sites with sequences from the high throughput genome sequencing project (HTG). Out of 171 human non-canonical and EST-supported splice pairs, 156 (91.23%) had a clear match in the human HTG. They can be classified after sequence analysis as: 79 GC–AG pairs (of which one was an error that corrected to GC–AG), 61 errors corrected to GT–AG canonical pairs, six AT–AC pairs (of which two were errors corrected to AT–AC), one case was produced from a non-existent intron, seven cases were found in HTG that were deposited to GenBank and finally there were only two other cases left of supported non-canonical splice pairs. The information about verified splice site sequences for canonical and non-canonical sites is presented in SpliceDB with the supporting evidence. We also built weight matrices for the major splice groups, which can be incorporated into gene prediction programs. SpliceDB is available at the computational genomic Web server of the Sanger Centre: http://genomic.sanger.ac.uk/spldb/SpliceDB.html and at http://www.softberry.com/spldb/SpliceDB.html.

Distinct roles for the N- and C-terminal regions in the cytotoxicity of pierisin-1, a putative ADP-ribosylating toxin from cabbage butterfly, against mammalian cells

Pierisin-1 is an 850-aa cytotoxic protein found in the cabbage butterfly, Pieris rapae, and has been suggested to consist of an N-terminal region with ADP-ribosyltransferase domain and of a C-terminal region that might have a receptor-binding domain. To elucidate the role of each region, we investigated the functions of various fragments of pierisin-1. In vitro expressed polypeptide consisting of amino acid residues 1–233 or 234–850 of pierisin-1 alone did not show cytotoxicity against human cervical carcinoma HeLa cells. However, the presence of both polypeptides in the culture medium showed some of the original cytotoxic activity. Introduction of the N-terminal polypeptide alone by electroporation also induced cell death in HeLa cells, and even in the mouse melanoma MEB4 cells insensitive to pierisin-1. Thus, the N-terminal region has a principal role in the cytotoxicity of pierisin-1 inside mammalian cells. Analyses of incorporated pierisin-1 indicated that the entire protein, regardless of whether it consisted of a single polypeptide or two separate N- and C-terminal polypeptides, was incorporated into HeLa cells. However, neither of the terminal polypeptides was incorporated when each polypeptide was present separately. These findings indicate that the C-terminal region is important for the incorporation of pierisin-1. Moreover, presence of receptor for pierisin-1 in the lipid fraction of cell membrane was suggested. The cytotoxic effects of pierisin-1 were enhanced by previous treatment with trypsin, producing “nicked” pierisin-1. Generation of the N-terminal fragment in HeLa cells was detected after application of intact entire molecule of pierisin-1. From the above observations, it is suggested that after incorporation of pierisin-1 into the cell by interaction of its C-terminal region with the receptor in the cell membrane, the entire protein is cleaved into the N- and C-terminal fragments with intracellular protease, and the N-terminal fragment then exhibits cytotoxicity.

Mammalian inositol polyphosphate multikinase synthesizes inositol 1,4,5-trisphosphate and an inositol pyrophosphate

Using a consensus sequence in inositol phosphate kinase, we have identified and cloned a 44-kDa mammalian inositol phosphate kinase with broader catalytic capacities than any other member of the family and which we designate mammalian inositol phosphate multikinase (mIPMK). By phosphorylating inositol 4,5-bisphosphate, mIPMK provides an alternative biosynthesis for inositol 1,4,5-trisphosphate [Ins(1,4,5)P3]. mIPMK also can form the pyrophosphate disphosphoinositol tetrakisphosphate (PP-InsP4) from InsP5. Additionally, mIPMK forms InsP4 from Ins(1,4,5)P3 and InsP5 from Ins(1,3,4,5)P4.

Analysis of endoplasmic reticulum trafficking signals by combinatorial screening in mammalian cells

To improve the accuracy of predicting membrane protein sorting signals, we developed a general methodology for defining trafficking signal consensus sequences in the environment of the living cell. Our approach uses retroviral gene transfer to create combinatorial expression libraries of trafficking signal variants in mammalian cells, flow cytometry to sort cells based on trafficking phenotype, and quantitative trafficking assays to measure the efficacy of individual signals. Using this strategy to analyze arginine- and lysine-based endoplasmic reticulum localization signals, we demonstrate that small changes in the local sequence context dramatically alter signal strength, generating a broad spectrum of trafficking phenotypes. Finally, using sequences from our screen, we found that the potency of di-lysine, but not di-arginine, mediated endoplasmic reticulum localization was correlated with the strength of interaction with α-COP.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Enhanced activity of adenine-DNA glycosylase (Myh) by apurinic/apyrimidinic endonuclease (Ape1) in mammalian base excision repair of an A/GO mismatch

Adenine-DNA glycosylase MutY of Escherichia coli catalyzes the cleavage of adenine when mismatched with 7,8-dihydro-8-oxoguanine (GO), an oxidatively damaged base. The biological outcome is the prevention of C/G→A/T transversions. The molecular mechanism of base excision repair (BER) of A/GO in mammals is not well understood. In this study we report stimulation of mammalian adenine-DNA glycosylase activity by apurinic/apyrimidinic (AP) endonuclease using murine homolog of MutY (Myh) and human AP endonuclease (Ape1), which shares 94% amino acid identity with its murine homolog Apex. After removal of adenine by the Myh glycosylase activity, intact AP DNA remains due to lack of an efficient Myh AP lyase activity. The study of wild-type Ape1 and its catalytic mutant H309N demonstrates that Ape1 catalytic activity is required for formation of cleaved AP DNA. It also appears that Ape1 stimulates Myh glycosylase activity by increasing formation of the Myh–DNA complex. This stimulation is independent of the catalytic activity of Ape1. Consequently, Ape1 preserves the Myh preference for A/GO over A/G and improves overall glycosylase efficiency. Our study suggests that protein–protein interactions may occur in vivo to achieve efficient BER of A/GO.

RAD51 supports spontaneous non-homologous recombination in mammalian cells, but not the corresponding process induced by topoisomerase inhibitors

The RAD51 protein has been shown to participate in homologous recombination by promoting ATP-dependent homologous pairing and strand transfer reactions. In the present study, we have investigated the possible involvement of RAD51 in non-homologous recombination. We demonstrate that overexpression of CgRAD51 enhances the frequency of spontaneous non-homologous recombination in the hprt gene of Chinese hamster cells. However, the rate of non-homologous recombination induced by the topoisomerase inhibitors campothecin and etoposide was not altered by overexpression of RAD51. These results indicate that the RAD51 protein may perform a function in connection with spontaneous non-homologous recombination that is not essential to or not rate-limiting for non-homologous recombination induced by camptothecin or etoposide. We discuss the possibility that the role played by RAD51 in non-homologous recombination observed here may not be linked to non-homologous end-joining.