Amino acid sequence studies on the alpha chain of human fibrinogen. Covalent structure of the alpha-chain portion of fragment D.

The alpha-chain portion of fragment D has been purified from an exhaustive plasmic digest of human fibrinogen. The major polypeptide species has 91 amino acid residues, although a small amount of a 97-residue chain representing an earlier digestion stage remains. The amino acid sequence of the first 44 residues was determined by stepwise degradation with an automatic solid-phase sequencer. Another large stretch of sequence was revealed by the finding that the alpha chain of fragment D overlaps the cyanogen bromide fragments alphaCNIVA and alphaCNIII (Doolittle, R. F. Cassman, K. G., Cottrell, B. A., Friezner, S. J. Hucko, J. T., and Takagi, T. (1977), Biochemistry 16 (preceding paper in this issue)). The automatic sequencer results were confirmed and extended by the isolation and characterization of 18 of 19 expected tryptic peptides from the fragment D alpha chain. As a result, almost the entire sequence has been obtained. The overlap with key cyanogen bromide fragments has also allowed us to propose an order for the first 198 residues of the fibrinogen alpha chain. A striking homology with the gamma chain and beta chain is apparent which has interesting structural implications.     

A contig assembly program based on sensitive detection of fragment overlaps.

An effective computer program for assembling DNA fragments, the contig assembly program (CAP), has been developed. In the CAP program, a filter is used to eliminate quickly fragment pairs that could not possibly overlap, a dynamic programming algorithm is applied to compute the maximal-scoring overlapping alignment between each remaining pair of fragments, and a simple greedy approach is employed to assemble fragments in order of alignment scores. To identify the true fragment overlaps, the dynamic programming algorithm uses specially chosen sets of alignment parameters to tolerate sequencing errors and to penalize "mutational" changes between different copies of a repetitive sequence. The performance tests of the program on fragment data from genomic sequencing projects produced satisfactory results. The CAP program is efficient in computer time and memory; it took about 4 h to assemble a set of 1015 fragments into long contigs on a Sun workstation.     

Construction of the temperature-sensitive vectors pLUCH80 and pLUCH88 for delivery of Tn917::NotI/SmaI and use of these vectors to derive a circular map of Listeria monocytogenes Scott A, a serotype 4b isolate.

A physical map of Listeria monocytogenes Scott A was generated by the pulsed-field technique of contour-clamped-homogeneous-electric-field (CHEF) electrophoresis. The circular genome of this serotype 4b strain contains 12 AscI fragments (38 to 790 kb), 5 NotI fragments (55 to 1,400 kb), 3 SrfI fragments (110, 1,110, and 2,000 kb), and 2 SfiI fragments (1,320 and 1,920 kb). Summation of individually sized fragments derived by digestion of Scott A genomic DNA with each of these four enzymes provided an average estimated genome length of 3,210 +/- 60 kb. Efforts to assemble the macrorestriction map benefited greatly from the construction and use of pLUCH80 and pLUCH88, temperature-sensitive vectors for delivering transposon Tn917::NotI/SmaI to the chromosome of Scott A. As another component of this study, the positions of four known virulence genes (inlA, mpl, hly, and prf) and three L. monocytogenes-specific sequences (lisM44, lisM51, and lisM52) were localized on the physical map of Scott A by hybridization. Probes prepared from lisM44, lisM51, and the four virulence genes hybridized within a cluster on a 150-kb fragment of the Scott A genome that overlaps part of the NotI-B and AscI-D fragments. The lisM52 probe hybridized with the AscI-F2 (120-kb) fragment of Scott A, which is separated from the NotI-B-AscI-D region by about 300 kb. These results established the first physical and genetic map of a serotype 4b strain of L. monocytogenes and provided further insight on this important food-borne pathogen at the genome level.     

Coprolites as a source of information on the genome and diet of the cave hyena

We performed high-throughput sequencing of DNA from fossilized faeces to evaluate this material as a source of information on the genome and diet of Pleistocene carnivores. We analysed coprolites derived from the extinct cave hyena (Crocuta crocuta spelaea), and sequenced 90 million DNA fragments from two specimens. The DNA reads enabled a reconstruction of the cave hyena mitochondrial genome with up to a 158-fold coverage. This genome, and those sequenced from extant spotted (Crocuta crocuta) and striped (Hyaena hyaena) hyena specimens, allows for the establishment of a robust phylogeny that supports a close relationship between the cave and the spotted hyena. We also demonstrate that high-throughput sequencing yields data for cave hyena multi-copy and single-copy nuclear genes, and that about 50 per cent of the coprolite DNA can be ascribed to this species. Analysing the data for additional species to indicate the cave hyena diet, we retrieved abundant sequences for the red deer (Cervus elaphus), and characterized its mitochondrial genome with up to a 3.8-fold coverage. In conclusion, we have demonstrated the presence of abundant ancient DNA in the coprolites surveyed. Shotgun sequencing of this material yielded a wealth of DNA sequences for a Pleistocene carnivore and allowed unbiased identification of diet.     

Improving Phrap-based assembly of the rat using "reliable" overlaps

The assembly methods used for whole-genome shotgun (WGS) data have a major impact on the quality of resulting draft genomes. We present a novel algorithm to generate a set of "reliable" overlaps based on identifying repeat k-mers. To demonstrate the benefits of using reliable overlaps, we have created a version of the Phrap assembly program that uses only overlaps from a specific list. We call this version PhrapUMD. Integrating PhrapUMD and our "reliable-overlap" algorithm with the Baylor College of Medicine assembler, Atlas, we assemble the BACs from the Rattus norvegicus genome project. Starting with the same data as the Nov. 2002 Atlas assembly, we compare our results and the Atlas assembly to the 4.3 Mb of rat sequence in the 21 BACs that have been finished. Our version of the draft assembly of the 21 BACs increases the coverage of finished sequence from 93.4% to 96.3%, while simultaneously reducing the base error rate from 4.5 to 1.1 errors per 10,000 bases. There are a number of ways of assessing the relative merits of assemblies when the finished sequence is available. If one views the overall quality of an assembly as proportional to the inverse of the product of the error rate and sequence missed, then the assembly presented here is seven times better. The UMD Overlapper with options for reliable overlaps is available from the authors at http://www.genome.umd.edu. We also provide the changes to the Phrap source code enabling it to use only the reliable overlaps.     

Improved inverse PCR scheme for metagenome walking

Inverse PCR has been used for the recovery of genome regions flanking a known sequence, although its application to metagenome walking is limited due to inefficient amplification from low copy number fragments. Here we present an improved inverse PCR scheme that enables walking of rare fragments in environmental metagenomes. Our scheme includes the following steps: (i) inverse PCR in which one primer is connected to an affinity tag; (ii) affinity purification of PCR products for removing background metagenome; and (iii) nested PCR to recover target flanking regions (IAN-PCR). In a model experiment, flanking regions of a gene fragment in Ralstonia eutropha were recovered from mixtures of Ralstonia and Escherichia genomes by standard inverse PCR, inverse PCR coupled to nested PCR (IN-PCR), and IAN-PCR, showing that they were recovered when ratios of Ralstonia genome to the background Escherichia genome were greater than 10(-1), 10(-3), and 10(-5), respectively. The utility of IAN-PCR was also examined by recovering flanking regions of PCR-amplified putative chitinase gene fragments from a groundwater metagenome, showing that IAN-PCR obtained flanking regions for more diverse gene fragments than IN-PCR. Since rare sequences are a critical element of natural genetic diversity, IAN-PCR enables access to undiscovered diverse genes in the environment.     

Monoclonal antibody-based, selective isolation of DNA fragments containing an alkylated base to be quantified in defined gene sequences.

We have established a sensitive, monoclonal antibody (Mab)-based procedure permitting the selective enrichment of sequences containing the miscoding alkylation product O6-ethylguanine (O6-EtGua) from mammalian DNA. H5 rat hepatoma cells were reacted with the N-nitroso carcinogen N-ethyl-N-nitrosourea in vitro, to give overall levels of greater than or equal to 25 O6-EtGua residues per diploid genome (corresponding to O6-EtGua/guanine molar ratios of greater than or equal to 10(-8). For analysis, enzymatically restricted DNA from these cells is incubated with an antibody specific for O6-ethyl-2'-deoxyguanosine, the resulting Mab-DNA complexes are separated from (O6-EtGua)-free fragments by filtration through a nitrocellulose (NC) membrane, and the DNA is recovered from the filter-bound complexes quantitatively. The efficiency of Mab binding to DNA fragments containing O6-EtGua is constant over a range of O6-EtGua/guanine molar ratios between 10(-5) and 10(-8). (O6-EtGua)-containing restriction fragments encompassing known gene sequences (e.g., the immunoglobulin E heavy chain gene of H5 rat hepatoma cells used as a model in this study) are subsequently amplified by PCR and quantified by slot-blot hybridisation. The content and distribution of a specific carcinogen-DNA adduct in defined sequences of genomic DNA can thus be analyzed as well as the kinetics of intragenomic (toposelective) repair of any DNA lesion for which a suitable Mab is available.     

Molecular cloning and physical mapping of murine cytomegalovirus DNA.

Murine cytomegalovirus (MCMV) Smith strain DNA is cleaved by restriction endonuclease HindIII into 16 fragments, ranging in size from 0.64 to 22.25 megadaltons. Of the 16 HindIII fragments, 15 were cloned in plasmid pACYC177 in Escherichia coli HB101 (recA). The recombinant plasmid clones were characterized by cleavage with the enzymes XbaI and EcoRI. In addition, fragments generated by double digestion of cloned fragments with HindIII and XbaI were inserted into the plasmid vector pACYC184. The results obtained after hybridization of 32P-labeled cloned fragments to Southern blots of MCMV DNA cleaved with HindIII, XbaI, EcoRI, BamHI, ApaI, ClaI, EcoRV, or KpnI allowed us to construct complete physical maps of the viral DNA for the restriction endonucleases HindIII, XbaI, and EcoRI. On the basis of the cloning and mapping experiments, it was calculated that the MCMV genome spans about 235 kilobase pairs, corresponding to a molecular weight of 155,000,000. All fragments were found to be present in equimolar concentrations, and no cross-hybridization between any of the fragments was seen. We conclude that the MCMV DNA molecule consists of a long unique sequence without large terminal or internal repeat regions. Thus, the structural organization of the MCMV genome is fundamentally different from that of the human cytomegalovirus or herpes simplex virus genome.     

Nonspecific primer and PCR generated hybridization probes for physical ordering large restriction fragments in complex genome of S. aureus.

Pulsed field gel electrophoresis (PFGE) allows separation of large restriction fragments from bacterial genome. Restriction fragments obtained by digestion of Staphylococcus aureus DNA with rare cutting enzymes (Sma I, and Csp I) were separated by PFGE. To arrange the physical order of the fragments generated by digestion with one enzyme, probes were prepared by nonspecific priming and polymerase chain reaction (PCR), using individual fragments of the other enzymatic digest as a template. Probes were then used for Southern hybridization to the PFGE separated fragment distribution of the two infrequent cleaving enzymes (Sma I and Csp I). Using probes generated from four Sma I fragments and five Csp I fragments as individual templates, a partial physical order of Csp I fragments of the genome of S. aureus ISP8 has been determined in relation to a previously published Sma I map of S. aureus genome.     

DNABIT Compress - Genome compression algorithm

Data compression is concerned with how information is organized in data. Efficient storage means removal of redundancy from the data being stored in the DNA molecule. Data compression algorithms remove redundancy and are used to understand biologically important molecules. We present a compression algorithm, "DNABIT Compress" for DNA sequences based on a novel algorithm of assigning binary bits for smaller segments of DNA bases to compress both repetitive and non repetitive DNA sequence. Our proposed algorithm achieves the best compression ratio for DNA sequences for larger genome. Significantly better compression results show that "DNABIT Compress" algorithm is the best among the remaining compression algorithms. While achieving the best compression ratios for DNA sequences (Genomes),our new DNABIT Compress algorithm significantly improves the running time of all previous DNA compression programs. Assigning binary bits (Unique BIT CODE) for (Exact Repeats, Reverse Repeats) fragments of DNA sequence is also a unique concept introduced in this algorithm for the first time in DNA compression. This proposed new algorithm could achieve the best compression ratio as much as 1.58 bits/bases where the existing best methods could not achieve a ratio less than 1.72 bits/bases.     

A method for generating selective DNA probes for the analysis of C-negative regions in human chromosomes

Linker-adapter polymerase chain reaction (LA-PCR) is among the most efficient techniques for whole genome DNA amplification. The key stage in LA-PCR is the hydrolysis of a DNA sample with restriction endonucleases, and the choice of a restriction endonuclease (or several endonucleases) determines the composition of DNA probes generated in LA-PCR. Computer analysis of the localization of the restriction sites in human genome has allowed us to propose an efficient technique for generating DNA probes by LA-PCR using the restriction endonucleases HaeIII and RsaI. In silico hydrolysis of human genomic DNA with endonucleases HaeIII and RsaI demonstrate that 100- to 1,000-bp DNA fragments are more abundant in the gene-rich regions. Applying in situ hybridization to metaphase chromosomes, we demonstrated that the produced DNA probes predominantly hybridized to the C-negative chromosomal regions, whereas the FISH signal was almost absent in the C-positive regions. The described protocol for generating DNA probes may be successfully used in subsequent cytogenetic analysis of the C-negative chromosomal regions.     

Sorting by weighted reversals, transpositions, and inverted transpositions.

During evolution, genomes are subject to genome rearrangements that alter the ordering and orientation of genes on the chromosomes. If a genome consists of a single chromosome (like mitochondrial, chloroplast, or bacterial genomes), the biologically relevant genome rearrangements are (1) inversions--also called reversals--where a section of the genome is excised, reversed in orientation, and reinserted and (2) transpositions, where a section of the genome is excised and reinserted at a new position in the genome; if this also involves an inversion, one speaks of an inverted transposition. To reconstruct ancient events in the evolutionary history of organisms, one is interested in finding an optimal sequence of genome rearrangements that transforms a given genome into another genome. It is well known that this problem is equivalent to the problem of "sorting" a signed permutation into the identity permutation. In this paper, we provide a 1.5-approximation algorithm for sorting by weighted reversals, transpositions and inverted transpositions for biologically realistic weights.     

Synthesis of Penicillins and Cephalosporins by Penicillin Acylase Entrapped in Fibres

Tryptic peptides from two cyanogen bromide (CNBr) fragments CB II and CB III of the Ala chain of ricin D were sequenced by manual Edman degradation. Chymotryptic or peptic peptides from the two fragments were isolated by Dowex 1 x 2 column chromatography to obtain overlaps for the tryptic peptides, and the complete amino acid sequences of fragments CB II and III were established. The amino acid residues in fragments CB II and CB III accounted for 75 and 45 residues, respectively, of 260 residues in the Ala chain.

These sequences together with the sequence of fragment CBI described in the preceding paper established the complete sequence of the 260 amino acid residues in the Ala chain. Some structural characteristics of the protein are also discussed.     

Characterization of the guinea pig cytomegalovirus genome by molecular cloning and physical mapping.

Fragments of guinea pig cytomegalovirus (GPCMV) DNA produced by HindIII or EcoRI restriction endonuclease digestion were cloned into vectors pBR322 and pACYC184, and recombinant fragments representing ca. 97% of the genome were constructed. Hybridization of 32P-labeled cloned and gel-purified HindIII, EcoRI, and XbaI fragments to Southern blots of HindIII-, EcoRI-, and XbaI-cleaved GPCMV DNA verified the viral origin of cloned fragments and allowed construction of HindIII, EcoRI, and XbaI restriction maps. On the basis of the cloning and mapping experiments, the size of GPCMV DNA was calculated to include 239 kilobase pairs, corresponding to a molecular weight of 158 X 10(6). No cross-hybridization between any internal fragments was seen. We conclude that the GPCMV genome consists of a long unique sequence with terminal repeat sequences but without internal repeat regions. In addition, GPCMV DNA molecules exist in two forms. In the predominant form, the molecules demonstrate sequence homology between the terminal fragments; in the minor population, one terminal fragment is smaller by 0.7 X 10(6) daltons and is not homologous with the fragment at the other end of the physical map. The structural organization of GPCMV DNA is unique for a herpesvirus DNA, similar in its simplicity to the structure reported for murine cytomegalovirus DNA and quite dissimilar from that of human cytomegalovirus DNA.     

DNA Fragments Assembly Based on Nicking Enzyme System

A couple of DNA ligation-independent cloning (LIC) methods have been reported to meet various requirements in metabolic engineering and synthetic biology. The principle of LIC is the assembly of multiple overlapping DNA fragments by single-stranded (ss) DNA overlaps annealing. Here we present a method to generate single-stranded DNA overlaps based on Nicking Endonucleases (NEases) for LIC, the method was termed NE-LIC. Factors related to cloning efficiency were optimized in this study. This NE-LIC allows generating 3′-end or 5′-end ss DNA overlaps of various lengths for fragments assembly. We demonstrated that the 10 bp/15 bp overlaps had the highest DNA fragments assembling efficiency, while 5 bp/10 bp overlaps showed the highest efficiency when T4 DNA ligase was added. Its advantage over Sequence and Ligation Independent Cloning (SLIC) and Uracil-Specific Excision Reagent (USER) was obvious. The mechanism can be applied to many other LIC strategies. Finally, the NEases based LIC (NE-LIC) was successfully applied to assemble a pathway of six gene fragments responsible for synthesizing microbial poly-3-hydroxybutyrate (PHB).     

An in vitro strategy for the selective isolation of anomalous DNA from prokaryotic genomes

In sequenced genomes of prokaryotes, anomalous DNA (aDNA) can be recognized, among others, by atypical clustering of dinucleotides. We hypothesized that atypical clustering of hexameric endonuclease recognition sites in aDNA allows the specific isolation of anomalous sequences in vitro. Clustering of endonuclease recognition sites in aDNA regions of eight published prokaryotic genome sequences was demonstrated. In silico digestion of the Neisseria meningitidis MC58 genome, using four selected endonucleases, revealed that out of 27 of the small fragments predicted (<5 kb), 21 were located in known genomic islands. Of the 24 calculated fragments (>300 bp and <5 kb), 22 met our criteria for aDNA, i.e. a high dinucleotide dissimilarity and/or aberrant GC content. The four enzymes also allowed the identification of aDNA fragments from the related Z2491 strain. Similarly, the sequenced genomes of three strains of Escherichia coli assessed by in silico digestion using XbaI yielded strain-specific sets of fragments of anomalous composition. In vitro applicability of the method was demonstrated by using adaptor-linked PCR, yielding the predicted fragments from the N.meningitidis MC58 genome. In conclusion, this strategy allows the selective isolation of aDNA from prokaryotic genomes by a simple restriction digest–amplification–cloning–sequencing scheme.     

Molecular evolution of receptor-like kinase genes in hexaploid wheat. Independent evolution of orthologs after polyploidization and mechanisms of local rearrangements at paralogous loci

Hexaploid wheat is a young polyploid species and represents a good model to study mechanisms of gene evolution after polyploidization. Recent studies at the scale of the whole genome have suggested rapid genomic changes after polyploidization but so far the rearrangements that have occurred in terms of gene content and organization have not been analyzed at the microlevel in wheat. Here, we have isolated members of a receptor kinase (Lrk) gene family in hexaploid and diploid wheat, Aegilops tauschii, and barley (Hordeum vulgare). Phylogenetic analysis has allowed us to establish evolutionary relationships (orthology versus paralogy) between the different members of this gene family in wheat as well as with Lrk genes from barley. It also demonstrated that the sequences of the homoeologous Lrk genes evolved independently after polyploidization. In addition, we found evidence for gene loss during the evolution of wheat and barley. Analysis of large genomic fragments isolated from nonorthologous Lrk loci showed a high conservation of the gene content and gene organization at these loci on the homoeologous group 1 chromosomes of wheat and barley. Finally, sequence comparison of two paralogous fragments of chromosome 1B showed a large number of local events (sequence duplications, deletions, and insertions), which reveal rearrangements and mechanisms for genome enlargement at the microlevel.     

MIR: Family of repeats common to vertebrate genomes

We studied the occurrence of mammalian interspersed repeats (MIRs) in DNA and RNA of vertebrates, invertebrates, and bacteria using the data from GenBank. A special algorithm based on a weight position matrix with optimal alignment using dynamic programming was developed to search for the traces of MIR dissemination. This allowed us to search for highly divergent MIRs carrying deletions and insertions. MIRs were detected in genomes of various fishes, includingLatimeria. This suggests that the origin of MIRs dates back more than 400 million years. The method to search for similarity between highly divergent sequences may be used to find the genome fragments from various ancient repeat families and from various gene families.