Restriction endonuclease mapping of Co1E2-P9 and Co1E3-CA38 plasmids

Using single and double restriction-endonuclease digestions, 16 and 17 cleavage sites have been mapped for the ColE2-P9 and ColE3-CA38 plasmids, respectively. One or more sites for AvaI, BglI, EcoRI, HincII, PvuI, PvuII, SmaI and XhoI endonucleases were found in both plasmids, two BglII sites were found only in ColE2-P9, and one KpnI site was unique to ColE3-CA38. ColE2-P9 was found to be slightly smaller than ColE3-CA38,4.4 Md compared to 4.6 Md. Eleven restriction sites are common to both plasmids in that they are identically placed relative to each other. These sites define a continuous DNA segment equal to over 60% of each plasmid. The remaining portions of the plasmids, which contain the non-homologous regions identified by Inselburg and Johns (1975) have no restriction sites in common, and differ in size by about 0.2 Md.     

Structure, polymorphism, and novel repeated DNA elements revealed by a complete sequence of the human alpha-fetoprotein gene

The human alpha-fetoprotein gene spans 19,489 base pairs from the putative "Cap" site to the polyadenylation site. It is composed of 15 exons separated by 14 introns, which are symmetrically placed within the three domains of alpha-fetoprotein. In the 5' region, a putative TATAAA box is at position -21, and a variant sequence, CCAAC, of the common CAT box is at -65. Enhancer core sequences GTGGTTTAAAG are found in introns 3 and 4, and several copies of glucocorticoid response sequences AGATACAGTA are found on the template strand of the gene. There are six polymorphic sites within 4690 base pairs of contiguous DNA derived from two allelic alpha-fetoprotein genes. This amounts to a measured polymorphic frequency of 0.13%, or 6.4 X 10(-4)/site, which is about 5-10 times lower than values estimated from studies on polymorphic restriction sites in other regions of the human genome. There are four types of repetitive sequence elements in the introns and flanking regions of the human alpha-fetoprotein gene. At least one of these is apparently a novel structure (designated Xba) and is found as a pair of direct repeats, with one copy in intron 7 and the other in intron 8. It is conceivable that within the last 2 million years the copy in intron 8 gave rise to the repeat in intron 7. Their present location on both sides of exon 8 gives these sequences a potential for disrupting the functional integrity of the gene in the event of an unequal crossover between them. There are three Alu elements, one of which is in intron 4; the others are located in the 3' flanking region. A solitary Kpn repeat is found in intron 3. The Xba and Kpn repeats were only detected by complete sequencing of the introns. Neither X, Xba, nor Kpn elements are present in the related human albumin gene, whereas Alu's are present in different positions. From phylogenetic evidence, it appears that Alu elements were inserted into the alpha-fetoprotein gene at some time postdating the mammalian radiation 85 million years ago.     

Conversion of the FokI endonuclease to a universal restriction enzyme: cleavage of phage M13mp7 DNA at predetermined sites

Endonuclease FokI belongs to class IIS of restriction enzymes, for which the DNA cut points lie outside the enzyme-recognition sites. This permitted conferring new cleavage specificities by combining FokI with tailored oligodeoxynucleotide adapters. Such adapters carry a single-stranded (ss) target-recognition domain, complementary to the selected ss target DNA, and a double-stranded (ds) enzyme-recognition site. Neither enzyme nor adapter alone has endonucleolytic activity toward phage M13mp7 ss DNA, whereas the enzymeadapter complex cleaves this ss target DNA at the particular sites foreordained by the sequence of the ss domain of the adapter. Two kinds of adapters (32 and 34 nucleotides long), with opposing orientations of the asymmetric FokI recognition site, were constructed and shown to direct specific cleavage under a variety of conditions. In addition to FokI, other class IIS enzymes, HphI, MboII and BbvI, which alone do not cleave ss DNA, are suitable for construction of tailored enzyme-adapter complexes with predictable cleavage specificities.

This report provides a preliminary experimental confirmation for the proposal of Szybalski [Gene 40 (1985) 169-173] for the design of adapter-enzyme complexes with novel and predictable specificities. Theoretically, using this approach any sequence could be precisely cleaved at a predetermined point.     

KpnI families of long, interspersed repetitive DNAs associated with the human β-globin gene cluster

KpnI families of long, interspersed repetitive DNAs are ubiquitous repetitive elements that occur in tens of thousands of copies in primate genomes. KpnI 1.2, 1.5 and two different KpnI 1.8-kb families were found within and flanking a 6.4-kb repeat beginning at 3 kb, 3' from the human β-globin gene. Thus, six different types of KpnI families have now been identified, and four of these are found next to each other in a specific 6.4-kb repeat. Clones of the distinct KpnI families were hybridized to clones of the 6.4-kb repeat and adjacent sequences encompassed within some 17.6 kb of DNA lying 3' to the β-globin gene cluster. The four KpnI families appear to make up the entire length of the 6.4-kb repeat. The linear order of the various cloned KpnI sequences in the repeat is 5'-pBK(1.8)26-pBK.(1.5)54-pBK(1.2)11-pBK(1.8)11-3'. KpnI 1.2-kb sequences were also detected downstream from the 6.4-kb repeat. As in the case of the KpnI 1.2 and 1.5-kb families, the two KpnI 1.8-kb sequence families described here each hybridized with about 15% of all plaques in two independently generated human genome libraries.     

Discrete size classes of monkey extrachromosomal circular DNA containing the L1 family of long interspersed nucleotide sequences are produced by a general non-sequence specific mechanism

The L1 family of long interspersed nucleotide sequences (LINES) has recently been identified and characterized in the small polydisperse circular DNA (spc-DNA) populations of monkey (1), human (2) and mouse (3) cells. In monkey spc-DNA, the L1 (also known as Kpn I) family is present in discrete size classes (ranging from 300 to 6000 base pairs (bp)) which appear to be generated by non homologous recombination events within chromosomal elements. In this communication it is shown that different regions of the consensus L1 family are present at different frequencies in monkey spc-DNA (as they are in chromosomal DNA), that all regions of the family are present in extrachromosomal DNA, and that each region appears to be represented in an identical discrete spc-DNA size distribution. This size distribution reflects a non-sequence specific mechanism that generates spc-DNA size classes by chromosomal DNA recombination events that are in some way constrained to occur between sites separated by relatively defined lengths.     

The construction of a versatile plasmid vector that allows direct selection of fragments cloned into six unique sites of the cI gene of coliphage 434

A new plasmid vector, pNSI, is described that allows positive selection for bacterial transformants carrying recombinant plasmids. It is a derivative of pBR327, and it includes a regulatory region from the lambdoid phage 434. The expression of the Tc^R gene of pNS1 is under the control of the o_Rp_R operator-promoter of phage 434, which is regulated by the represser gene c1. The cloning sites of pNSI (StuI, NdeI, HpaI, HindIII, AsuII and EcoRI) are situated within cI; hence insertion of foreign DNA into these sites causes derepressed expression of the Tc^R gene from p_R thus conferring the Tc^R phenotype on the harboring Escherichia coli strain. The use ofpNS1 is facilitated by the presence of another selectable marker, Ap^R its small size, and its known nucleotide sequence; no special host strain is required.     

Genomic representation of the Hind II 1.9 kb repeated DNA.

The genomic representation and organization of sequences homologous to a cloned Hind III 1.9 kb repeated DNA fragment were studied. Approximately 80% of homologous repeated DNA was contained in a genomic Hind III cleavage band of 1.9 kb. Double digestion studies indicated that the genomic family, in the majority, followed the arrangement of the sequenced clone, with minor restriction cleavage variations compatible with a few base changes. Common restriction sites external to the 1.9 kb sequence were mapped, and hybridization of segments of the cloned sequence indicated the 1.9 kb DNA was itself not tandemly repeated. Kpn I bands which were homologous to the sequence contained specific regions of the repeat, and the molecular weight of these larger fragments could be simply explained. Mapping of common external restriction sites indicated that in some but not all cases the repeat could be organized in larger defined blocks of greater than or equal to 5.5 kb. In some instances, flanking regions adjacent to the repeat may contain common DNA elements such as other repeated DNA sequences, or possibly rearranged segments of the 1.9 kb sequence. It is suggested that although the 1.9 kb sequence is not strictly contiguous, at least some of these repeated sequences in the human genome are arranged in clustered or intercalary arrays. A region of the 1.9 kb sequence hybridized to a mouse repeated DNA, indicating homology beyond the primates.     

Sequences near the 3'' secretion-specific polyadenylation site influence levels of secretion-specific and membrane-specific IgG2b mRNA in myeloma cells.

The expressed immunoglobulin gamma 2b (IgG2b) heavy-chain gene of 4T001 was cloned into the shuttle vector pSV2-gpt and transfected into myeloma J558L and lymphoma A20.2J. Northern blots indicated that the transfected gamma 2b gene was processed in a manner similar to the endogenous heavy chain in both lymphoma and myeloma cells. To identify sequences important for immunoglobulin mRNA processing, we constructed deletions around the secretion-specific polyadenylation site and introduced the deleted genes into J558L cells. The BAL deletion lacked 670 base pairs of intervening sequence between secreted and membrane regions; the Kpn deletion lacked 830 base pairs in this region. J558L cells transfected with either the entire gamma 2b gene or the delta BAL vector produced predominantly secretion-specific gamma 2b mRNA and protein. J558L cells transfected with the delta Kpn vector produced approximately equimolar amounts of secretion-specific and membrane-specific gamma 2b mRNA. Both 55,000-dalton secreted and 62,000-dalton putative surface IgG2b proteins were detected in the delta Kpn transfectants. We conclude that sequences absent in the Kpn deletion but present in the BAL deletion exert an important role in the production of secretion-specific mRNA. The Kpn deletion removes the normal site of cleavage and poly(A) addition, and it is possible that it is the absence of this site which changes the processing pattern. Alternatively, it is possible that sequences absent in the Kpn deletion but present in the BAL deletion function in regulating the production of predominantly secretion-specific mRNA in myeloma cells. The possible role of a highly conserved sequence found in this region is discussed.     

ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites.

We present a neural network based method (ChloroP) for identifying chloroplast transit peptides and their cleavage sites. Using cross-validation, 88% of the sequences in our homology reduced training set were correctly classified as transit peptides or nontransit peptides. This performance level is well above that of the publicly available chloroplast localization predictor PSORT. Cleavage sites are predicted using a scoring matrix derived by an automatic motif-finding algorithm. Approximately 60% of the known cleavage sites in our sequence collection were predicted to within +/-2 residues from the cleavage sites given in SWISS-PROT. An analysis of 715 Arabidopsis thaliana sequences from SWISS-PROT suggests that the ChloroP method should be useful for the identification of putative transit peptides in genome-wide sequence data. The ChloroP predictor is available as a web-server at http://www.cbs.dtu.dk/services/ChloroP/.     

Submicromolar free calcium modulates dexamethasone binding to the glucocorticoid receptor

The relative distribution of bound cis- and trans-(NH₃)₂PtCl₂ at specific sites in SV40 DNA is evaluated by monitoring the extent to which five restriction endonucleases, each of which cleave at a single, unique site, are inhibited as a result of the DNA modification. The order of cleavage inhibition is Bgl 1 ? Bam HI > Hpa II, Kpn I > Eco RI. Both isomers produce a comparable effect for any particular endonuclease. Inhibition correlates with the % (G+C) content within and about the recognition sequences. That modified sequences immediately adjacent to the recognition sequence influence cleavage is further supported by differential cleavage observed with the multicut Hind III endonuclease. The binding of cis-(NH₃)₂PtCl₂ at the hyper-reactive Bgl 1 site may well be directly responsible for inhibiting SV40 replication.     

Mitochondrial diversity within modern human populations

With the recent increase in the available number of high-quality, full-length mitochondrial sequences, it is now possible to construct and analyze a comprehensive human mitochondrial consensus sequence. Using a data set of 827 carefully selected sequences, it is shown that modern humans contain extremely low levels of divergence from the mitochondrial consensus sequence, differing by a mere 21.6nt sites on average. Fully 84.1% of the mitochondrial genome was found to be invariant and ‘private’ mutations accounted for 43.8% of the variable sites. Ninety eight percent of the variant sites had a primary nucleotide with an allele frequency of 0.90 or greater. Interestingly, the few truly ambiguous nucleotide sites could all be reliably assigned to either a purine or pyrimidine ancestral state. A comparison of this consensus sequence to several ancestral sequences derived from phylogenetic studies reveals a great deal of similarity, where, as expected, the most phylogenetically informative nucleotides in the ancestral studies tended to be the most variable nucleotides in the consensus. Allowing for this fact, the consensus approach provides variation data on the positions that do not contribute to phylogenetic reconstructions, and these data provide a baseline for measuring human mitochondrial variation in populations worldwide.     

LabCaS: Labeling calpain substrate cleavage sites from amino acid sequence using conditional random fields

The calpain family of Ca²⁺‐dependent cysteine proteases plays a vital role in many important biological processes which is closely related with a variety of pathological states. Activated calpains selectively cleave relevant substrates at specific cleavage sites, yielding multiple fragments that can have different functions from the intact substrate protein. Until now, our knowledge about the calpain functions and their substrate cleavage mechanisms are limited because the experimental determination and validation on calpain binding are usually laborious and expensive. In this work, we aim to develop a new computational approach (LabCaS) for accurate prediction of the calpain substrate cleavage sites from amino acid sequences. To overcome the imbalance of negative and positive samples in the machine‐learning training which have been suffered by most of the former approaches when splitting sequences into short peptides, we designed a conditional random field algorithm that can label the potential cleavage sites directly from the entire sequences. By integrating the multiple amino acid features and those derived from sequences, LabCaS achieves an accurate recognition of the cleave sites for most calpain proteins. In a jackknife test on a set of 129 benchmark proteins, LabCaS generates an AUC score 0.862. The LabCaS program is freely available at: http://www.csbio.sjtu.edu.cn/bioinf/LabCaS . Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Sequence and positional requirements for DNA sites in a mu transpososome.

Transposition of bacteriophage Mu uses two DNA cleavage sites and six transposase recognition sites, with each recognition site divided into two half-sites. The recognition sites can activate transposition of non-Mu DNA sequences if a complete set of Mu sequences is not available. We have analyzed 18 sequences from a non-Mu DNA molecule, selected in a functional assay for the ability to be transposed by MuA transposase. These sequences are remarkably diverse. Nonetheless, when viewed as a group they resemble a Mu DNA end, with a cleavage site and a single recognition site. Analysis of these "pseudo-Mu ends" indicates that most positions in the cleavage and recognition sites contribute sequence-specific information that helps drive transposition, though only the strongest contributors are apparent from mutagenesis data. The sequence analysis also suggests variability in the alignment of recognition half-sites. Transposition assays of specifically designed DNA substrates support the conclusion that the transposition machinery is flexible enough to permit variability in half-site spacing and also perhaps variability in the placement of the recognition site with respect to the cleavage site. This variability causes only local perturbations in the protein-DNA complex, as indicated by experiments in which altered and unaltered DNA substrates are paired.     

‘ATG vectors’ for regulated high-level expression of cloned genes in Escherichia coli

A plasmid cloning vector system has been constructed that allows for the production of large quantities of foreign proteins or fragments thereof, in an unfused state. These vectors provide strong regulated trp-lac fusion promoters and the lacZ ribosome-binding site (RBS) followed by an ATG translation initiation codon at an appropriate distance from the RBS. The ATG codon is located within a unique NcoI restriction site (CCATGG). Digestion with NcoI exposes the ATG for fusion. Gene fragments lacking a prokaryotic RBS and/or ATG start codons can be inserted in several ways. Expression experiments using a truncated cI gene of bacteriophage A or a large portion of the coding region of the Herpes simplex virus type l glycoprotein D gene have been performed. The results of these studies show that the vectors are useful for the high-level expression of prokaryotic and eukaryotic genes in Escherichia coli.     

A branch-and-bound algorithm for the inference of ancestral amino-acid sequences when the replacement rate varies among sites: Application to the evolution of five gene families

MOTIVATION: We developed an algorithm to reconstruct ancestral sequences, taking into account the rate variation among sites of the protein sequences. Our algorithm maximizes the joint probability of the ancestral sequences, assuming that the rate is gamma distributed among sites. Our algorithm probably finds the global maximum. The use of 'joint' reconstruction is motivated by studies that use the sequences at all the internal nodes in a phylogenetic tree, such as, for instance, the inference of patterns of amino-acid replacement, or tracing the biochemical changes that occurred during the evolution of a given protein family. RESULTS: We give an algorithm that guarantees finding the global maximum. The efficient search method makes our method applicable to datasets with large number sequences. We analyze ancestral sequences of five gene families, exploring the effect of the amount of among-site-rate-variation, and the degree of sequence divergence on the resulting ancestral states. AVAILABILITY AND SUPPLEMENTARY INFORMATION: http://evolu3.ism.ac.jp/~tal/ Contact: tal@ism.ac.jp     

Evolutionary conserved chromosomal segments in the human karyotype are bounded by unstable chromosome bands

In this paper an ancestral karyotype for primates, defining for the first time the ancestral chromosome morphology and the banding patterns, is proposed, and the ancestral syntenic chromosomal segments are identified in the human karyotype. The chromosomal bands that are boundaries of ancestral segments are identified. We have analyzed from data published in the literature 35 different primate species from 19 genera, using the order Scandentia, as well as other published mammalian species as out-groups, and propose an ancestral chromosome number of 2n = 54 for primates, which includes the following chromosomal forms: 1(a+c(1)), 1(b+c(2)), 2a, 2b, 3/21, 4, 5, 6, 7a, 7b, 8, 9, 10a, 10b, 11, 12a/22a, 12b/22b, 13, 14/15, 16a, 16b, 17, 18, 19a, 19b, 20 and X and Y. From this analysis, we have been able to point out the human chromosome bands more "prone" to breakage during the evolutionary pathways and/or pathology processes. We have observed that 89.09% of the human chromosome bands, which are boundaries for ancestral chromosome segments, contain common fragile sites and/or intrachromosomal telomeric-like sequences. A more in depth analysis of twelve different human chromosomes has allowed us to determine that 62.16% of the chromosomal bands implicated in inversions and 100% involved in fusions/fissions correspond to fragile sites, intrachromosomal telomeric-like sequences and/or bands significantly affected by X irradiation. In addition, 73% of the bands affected in pathological processes are co-localized in bands where fragile sites, intrachromosomal telomeric-like sequences, bands significantly affected by X irradiation and/or evolutionary chromosomal bands have been described. Our data also support the hypothesis that chromosomal breakages detected in pathological processes are not randomly distributed along the chromosomes, but rather concentrate in those important evolutionary chromosome bands which correspond to fragile sites and/or intrachromosomal telomeric-like sequences.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.