Application of partial restriction procedure in both shotgun and non-random strategies for nucleotide sequencing

Partial digestion of a target DNA fragment with 4-bp-recognition restriction enzymes followed by a forced ligation to an M13 vector was employed for the construction of a subfragment library. The library can be used for either shotgun or non-random nucleotide sequencing. Application of the partial digests generated with the 4-bp recognition restriction enzymes instead of DNase I in the improved non-random strategy for nucleotide sequencing (Li and Wu, 1987) made the procedure as easy as that of the random strategy. The library can also be used in shotgun nucleotide sequencing directly, and few self-ligated subfragments were found. The usefulness of this procedure was demonstrated by the sequencing of a goat 6.5-kb EcoRI fragment, which is located 5' to the globin gene.     

New strategy for mapping the human genome based on a novel procedure for construction of jumping libraries

A novel procedure for construction of jumping libraries is described. The essential features of this procedure are as follows: (1) two diphasmid vectors (lambda SK17 and lambda SK22) are simultaneously used in the library construction to improve representativity, (2) a partial filling-in reaction is used to eliminate cloning of artifactual jumping clones and to obviate the need for a selectable marker. The procedure has been used to construct a representative human NotI jumping library (220,000 independent recombinant clones) from the lymphoblastoid cell line CBMI-Ral-STO, which features a low level of methylation of its resident EBV genomes. A human chromosome 3-specific NotI jumping library (500,000 independent recombinant clones) from the human chromosome 3 x mouse hybrid cell line MCH 903.1 has also been constructed. Of these recombinant clones 50-80% represent jumps to the neighboring cleavable NotI site. With our previously published method for construction of linking libraries this procedure makes a new genome mapping strategy feasible. This strategy includes the determination of tagging sequences adjacent to NotI sites in random linking and jumping clones. Special features of the lambda SK17 and lambda SK22 vectors facilitate such sequencing. The STS (sequence tagged site) information obtained can be assembled by computer into a map representing the linear order of the NotI sites for a chromosome or for the entire genome. The computerized mapping data can be used to retrieve clones near a region of interest. The corresponding clones can be obtained from the panel of original clones, or necessary probes can be made from genomic DNA by PCR.     

Optimization of the genotyping‐by‐sequencing strategy for population genomic analysis in conifers

Flexibility and low cost make genotyping‐by‐sequencing (GBS) an ideal tool for population genomic studies of nonmodel species. However, to utilize the potential of the method fully, many parameters affecting library quality and single nucleotide polymorphism (SNP) discovery require optimization, especially for conifer genomes with a high repetitive DNA content. In this study, we explored strategies for effective GBS analysis in pine species. We constructed GBS libraries using HpaII, PstI and EcoRI‐MseI digestions with different multiplexing levels and examined the effect of restriction enzymes on library complexity and the impact of sequencing depth and size selection of restriction fragments on sequence coverage bias. We tested and compared UNEAK, Stacks and GATK pipelines for the GBS data, and then developed a reference‐free SNP calling strategy for haploid pine genomes. Our GBS procedure proved to be effective in SNP discovery, producing 7000–11 000 and 14 751 SNPs within and among three pine species, respectively, from a PstI library. This investigation provides guidance for the design and analysis of GBS experiments, particularly for organisms for which genomic information is lacking.     

Subunit interactions in myosin subfragment 1. Subunit scrambling is independent of catalytic center involvement

Evidence is presented that, under conditions of 4.7 M NH4Cl and 10 mM Mg-ATP where no subunit dissociation can be detected by transport methods, a dynamic equilibrium exists in subfragment 1 between the associated and dissociated subunits. This is readily discerned by the formation of hybrid subfragment 1 species when a subfragment 1 isozyme is incubated with excess free light chains of the alternate isozyme. A similar process occurs with p-N,N'-phenylenedimaleimide (pPDM)-modified subfragment 1 containing [14C]Mg-ADP, but in this case, although extensive amounts of hybrid are formed, no loss of the trapped nucleotide is observed. Subunit scrambling without loss of the trapped nucleotide is apparent from incubating pPDM-SF1(A2)-[14C]Mg-ADP with unmodified SF1(A1) under similar conditions since the mixture subsequently contains SF1(A1), SF1(A2)h, pPDM-SF1(A1)h-[14C]Mg-ADP and pPDM-SF1(A2)-[14C]Mg-ADP. These data show that the nucleotide trapped in the presumptive active site does not escape during the dissociation-reassociation cycle, and suggest that the ATPase site resides solely on the heavy chain.     

A new recombinant DNA strategy for the molecular cloning of rare membrane proteins.

We have constructed a cDNA library in the plasmid expression vector pUEX enriched in sequences encoding membrane proteins. The procedure involved positive selection of sequences common to two different rat tissues (thus excluding tissue-specific mRNA) followed by positive selection between this material and RNA extracted from membrane bound polysomes (thus excluding cytoplasmic proteins). The resultant library prepared from rat kidney cDNA hybridized with rat liver poly(A)+ RNA, contained 30,000 clones and was shown to be enriched in cDNAs encoding membrane proteins. Seventeen clones selected because they encode large fusion proteins were shown to be single copy in the library, and not present in nucleotide data banks. Thus the strategy is particularly suitable for cloning low abundance cDNAs encoding membrane proteins.     

Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing

ABSTRACT: BACKGROUND: Genetic mapping and QTL detection are powerful methodologies in plant improvement and breeding. Construction of a high-density and high-quality genetic map would be of great benefit in the production of superior grapes to meet human demand. High throughput and low cost of the recently developed next generation sequencing (NGS) technology have resulted in its wide application in genome research. Sequencing restriction-site associated DNA (RAD) might be an efficient strategy to simplify genotyping. Combining NGS with RAD has proven to be powerful for single nucleotide polymorphism (SNP) marker development. RESULTS: An F1 population of 100 individual plants was developed. In-silico digestion-site prediction was used to select an appropriate restriction enzyme for construction of a RAD sequencing library. Next generation RAD sequencing was applied to genotype the F1 population and its parents. Applying a cluster strategy for SNP modulation, a total of 1,814 high-quality SNP markers were developed: 1,121 of these were mapped to the female genetic map, 759 to the male map, and 1,646 to the integrated map. A comparison of the genetic maps to the published Vitis vinifera genome revealed both conservation and variations. CONCLUSIONS: The applicability of next generation RAD sequencing for genotyping a grape F1 population was demonstrated, leading to the successful development of a genetic map with high density and quality using our designed SNP markers. Detailed analysis revealed that this newly developed genetic map can be used for a variety of genome investigations, such as QTL detection, sequence assembly and genome comparison. KEYWORDS: Grape; Genetic map; Next generation sequencing (NGS); Restriction-site associated DNA (RAD).     

Genome Wide Sampling Sequencing for SNP Genotyping: Methods,Challenges and Future Development

Genetic polymorphisms, particularly single nucleotide polymorphisms (SNPs), have been widely used to advance quantitative, functional and evolutionary genomics. Ideally, all genetic variants among individuals should be discovered when next generation sequencing (NGS) technologies and platforms are used for whole genome sequencing or resequencing. In order to improve the cost-effectiveness of the process, however, the research community has mainly focused on developing genome-wide sampling sequencing (GWSS) methods, a collection of reduced genome complexity sequencing, reduced genome representation sequencing and selective genome target sequencing. Here we review the major steps involved in library preparation, the types of adapters used for ligation and the primers designed for amplification of ligated products for sequencing. Unfortunately, currently available GWSS methods have their drawbacks, such as inconsistency in the number of reads per sample library, the number of sites/targets per individual, and the number of reads per site/target, all of which result in missing data. Suggestions are proposed here to improve library construction, genotype calling accuracy, genome-wide marker density and read mapping rate. In brief, optimized GWSS library preparation should generate a unique set of target sites with dense distribution along chromosomes and even coverage per site across all individuals.     

An efficient strategy for large-scale high-throughput transposon-mediated sequencing of cDNA clones

We describe an efficient high-throughput method for accurate DNA sequencing of entire cDNA clones. Developed as part of our involvement in the Mammalian Gene Collection full-length cDNA sequencing initiative, the method has been used and refined in our laboratory since September 2000. Amenable to large scale projects, we have used the method to generate >7 Mb of accurate sequence from 3695 candidate full-length cDNAs. Sequencing is accomplished through the insertion of Mu transposon into cDNAs, followed by sequencing reactions primed with Mu-specific sequencing primers. Transposon insertion reactions are not performed with individual cDNAs but rather on pools of up to 96 clones. This pooling strategy reduces the number of transposon insertion sequencing libraries that would otherwise be required, reducing the costs and enhancing the efficiency of the transposon library construction procedure. Sequences generated using transposon-specific sequencing primers are assembled to yield the full-length cDNA sequence, with sequence editing and other sequence finishing activities performed as required to resolve sequence ambiguities. Although analysis of the many thousands (22 785) of sequenced Mu transposon insertion events revealed a weak sequence preference for Mu insertion, we observed insertion of the Mu transposon into 1015 of the possible 1024 5mer candidate insertion sites.     

Genotyping by Genome Reducing and Sequencing for Outbred Animals

Next-generation sequencing (NGS) approaches are widely used in genome-wide genetic marker discovery and genotyping. However, current NGS approaches are not easy to apply to general outbred populations (human and some major farm animals) for SNP identification because of the high level of heterogeneity and phase ambiguity in the haplotype. Here, we reported a new method for SNP genotyping, called genotyping by genome reducing and sequencing (GGRS) to genotype outbred species. Through an improved procedure for library preparation and a marker discovery and genotyping pipeline, the GGRS approach can genotype outbred species cost-effectively and high-reproducibly. We also evaluated the efficiency and accuracy of our approach for high-density SNP discovery and genotyping in a large genome pig species (2.8 Gb), for which more than 70,000 single nucleotide polymorphisms (SNPs) can be identified for an expenditure of only $80 (USD)/sample.     

The covalent modification of myosin's proteolytic fragments by a purine disulfide analog of adenosine triphosphate. Reaction at a binding site other than the active site.

A purine disulfide analog of ATP, 6,6'-dithiobis(inosinyl imidodiphosphate), forms mixed disulfide bonds between the 6 thiol group on the purine ring and certain key cysteines on myosin, heavy meromyosin, and subfragment one. The EDTA ATPase activities of myosin and heavy meromyosin were completely inactivated when 4 mol of thiopurine nucleotide was bound. When similarly inactivated, subfragment one, depending on its method of preparation, incorporated either 1 or 2 mol of thiopurine nucleotide. Modification of a single cysteine on subfragment one resulted in an inhibition of both the Ca2+ and the EDTA ATPase activities, but the latter always to a greater extent. Modification of two cysteines per head of heavy meromyosin had the same effect suggesting that the active sites were not blocked by the thiopurine nucleotides. Direct evidence for this suggestion was provided by equilibrium dialysis experiments. Heavy meromyosin and subfragment one bound 1.9 and 0.8 mol of [8-3H]adenylyl imidodiphosphate per mol of enzyme, respectively, with an average dissociation constant of 5 X 10(-7) M. Heavy meromyosin with four thiopurine nucleotides bound or subfragment one with two thiopurine nucleotides bound retained 65-80% of these tight adenylyl imidodiphosphate binding sites confirming the above suggestion. Thus previous work assuming reaction of thiopurine nucleotide analogs at the active site of myosin must be reevaluated. Ultracentrifugation studies showed that heavy meromyosin which had incorporated four thiopurine nucleotides did not bind to F-actin while subfragment one with one thiopurine nucleotide bound interacted only very weakly with F-actin. Thus reaction of 6,6'-dithiobis(inosinyl imidodiphosphate) at nucleotide binding sites other than the active sites reduces the rate of ATP hydrolysis and inhibits actin binding. It is suggested that these second sites may function as regulatory sites on myosin.     

Phagemid VPCS vectors for priming, cloning and sequencing

A phagemid was adapted for use as the vector in the vector-primer-cloner-sequencer cloning system. The use of this new vector markedly expanded the utility of this technology for the construction of cDNA libraries. Technological advantages and new capabilities include: (1) a greater number of unique restriction sites within the polylinker region; (2) the ability to produce single-stranded templates for nucleotide sequencing, and (3) a convenient means to synthesize strand-specific hybridization probes. With the use of this cloning system, a rat liver cDNA library (8.56 x 10(5) recombinants from 1 microgram of poly(A)+ RNA) was rapidly (in two days) constructed.     

Determination of the sequence of an expressible cDNA clone encoding ERp60/calregulin by the use of a novel nested set method

An analysis of the N-terminal sequence of the luminal endoplasmic reticulum protein, ERp60, showed that it was identical to the well-characterized Ca2+-binding protein, calregulin. A full-length, expressible cDNA clone encoding this protein was isolated from a mouse fibroblast cDNA library. A novel nested set strategy for the production of overlapping fragments for DNA sequencing was used to determine the complete nucleotide (nt) sequence of both strands of the ERp60 clone. This method utilizes a series of nonspecific deletion primers in conjunction with a specific site primer to generate the nested set fragments. This procedure possesses several advantages over other nested set techniques, since it does not require (i) the re-cloning of the DNA insert into other vectors, (ii) any prior knowledge of the restriction sites of the nt sequence, or (iii) the transformation and analysis of bacterial subclones. ERp60 has a 17-amino acid (aa) signal sequence and the mature protein contains 399 aa with a calculated M(r) of 46,347.     

Rapid verification of lambda-cDNA clones using mixed-base oligonucleotide as screening probe and sequencing primer

A rapid procedure has been developed for the isolation and verification of cDNA clones isolated from a cDNA library based on lambda vectors. Using information of the partial amino acid sequence of a protein, synthetic mixed-base oligonucleotides are first employed as a screening probe using the plaque hybridization procedure. The cDNA inserts of the clones obtained are then directly amplified by polymerase chain reaction (PCR) using primers flanking the cloning site of the vector. Besides being used for cloning into a plasmid vector, the amplified DNA's are also subjected to nucleotide sequence analysis using the same mixed-base oligonucleotides as sequencing primers. This approach allows sequencing through the region of the known amino acid sequence for direct verification of the authenticity of the clones obtained. This procedure has successfully been used for cloning and partial characterization of the gene coding for a platelet aggregation inhibitor.     

Kinetic mechanism of 1-N6-etheno-2-aza-ATP hydrolysis by bovine ventricular myosin subfragment 1 and actomyosin subfragment 1. The nucleotide binding steps

The large change in fluorescence emission of 1-N6-etheno-2-aza-ATP (epsilon-aza-ATP) has been used to investigate the kinetic mechanism of etheno-aza nucleotide binding to bovine cardiac myosin subfragment 1 (myosin-S1) and actomyosin subfragment 1 (actomyosin-S1). The time course of nucleotide fluorescence enhancement observed during epsilon-aza-ATP hydrolysis is qualitatively similar to the time course of tryptophan fluorescence enhancement observed during ATP hydrolysis. In single turnover experiments, the nucleotide fluorescence rapidly increases to a maximum level, then decreases with a rate constant of 0.045 s-1 to a final level, which is about 30% of the maximal enhancement; a similar fluorescence enhancement is obtained by adding epsilon-aza-ADP to cardiac myosin-S1 or actomyosin-S1 under the same conditions (100 mM KCl, 10 mM 4-morpholinepropanesulfonic acid, 5 mM MgCl2, 0.1 mM dithiothreitol, pH 7.0, 15 degrees C). The kinetic data are consistent with a mechanism in which there are two sequential (acto)myosin-S1 nucleotide complexes with enhanced nucleotide fluorescence following epsilon-aza-ATP binding. The apparent second order rate constants of epsilon-aza-ATP binding to cardiac myosin subfragment 1 and actomyosin subfragment 1 are 2-12 times slower than those for ATP. Actin increases the rate of epsilon-aza-ADP dissociation from bovine cardiac myosin-S1 from 1.9 to 110 s-1 at 15 degrees C which can be compared to 0.3 and 65 s-1 for ADP dissociation under similar conditions. Although there are quantitative differences between the rate and equilibrium constants of epsilon-aza- and adenosine nucleotides to cardiac actomyosin-S1 and myosin-S1, the basic features of the nucleotide binding steps of the mechanism are unchanged.     

Alternative dideoxy sequencing of double-stranded DNA by cyclic reactions using Taq polymerase.

The polymerase chain reaction (PCR) is a technique to amplify a specific DNA sequence millions of times. The thermostable enzyme Taq polymerase allows this procedure to take place under conditions of high specificity and automatization. By combining the techniques of PCR and dideoxy sequencing, it is possible to perform DNA sequencing independently of their structures. The cyclic sequencing reaction is carried out in the presence of an excess amount of sequencing primer and a radioactive nucleotide ([alpha-35S]dATP) using a DNA thermal cycler. Different reaction conditions were investigated and optimized including nucleotide ratios in each termination mix, primer/template ratios, amount of a radioactive nucleotide, and the program of the reaction. This method allows the detection of single base substitutions in heterozygous alleles, and the detection of homozygous deletions. A new RFLP of the human porphobilinogen deaminase (PBGD) gene was identified using this technique. This RFLP is created by one base difference (cytosine or adenine) that changes the restriction site for Apa LI. The alternative sequencing method described in this study is a simple and time-saving procedure that can also be used for large sequencing projects.     

A rapid,efficient, and economical inverse polymerase chain reaction-based method for generating a site saturation mutant library

With the development of deep sequencing methodologies, it has become important to construct site saturation mutant (SSM) libraries in which every nucleotide/codon in a gene is individually randomized. We describe methodologies for the rapid, efficient, and economical construction of such libraries using inverse polymerase chain reaction (PCR). We show that if the degenerate codon is in the middle of the mutagenic primer, there is an inherent PCR bias due to the thermodynamic mismatch penalty, which decreases the proportion of unique mutants. Introducing a nucleotide bias in the primer can alleviate the problem. Alternatively, if the degenerate codon is placed at the 5′ end, there is no PCR bias, which results in a higher proportion of unique mutants. This also facilitates detection of deletion mutants resulting from errors during primer synthesis. This method can be used to rapidly generate SSM libraries for any gene or nucleotide sequence, which can subsequently be screened and analyzed by deep sequencing.     

Massively parallel tag sequencing reveals the complexity of anaerobic marine protistan communities

Background  

Recent advances in sequencing strategies make possible unprecedented depth and scale of sampling for molecular detection of microbial diversity. Two major paradigm-shifting discoveries include the detection of bacterial diversity that is one to two orders of magnitude greater than previous estimates, and the discovery of an exciting 'rare biosphere' of molecular signatures ('species') of poorly understood ecological significance. We applied a high-throughput parallel tag sequencing (454 sequencing) protocol adopted for eukaryotes to investigate protistan community complexity in two contrasting anoxic marine ecosystems (Framvaren Fjord, Norway; Cariaco deep-sea basin, Venezuela). Both sampling sites have previously been scrutinized for protistan diversity by traditional clone library construction and Sanger sequencing. By comparing these clone library data with 454 amplicon library data, we assess the efficiency of high-throughput tag sequencing strategies. We here present a novel, highly conservative bioinformatic analysis pipeline for the processing of large tag sequence data sets.     

Rapid and reliable screening of a tomato YAC library exclusively based on PCR

An improved procedure is presented to select clones from a tomato yeast artificial chromosome (YAC) library. The procedure is based exlcusively on the polymerase chain reaction (PCR). We combined DNA from approximately 36,000 YAC clones in pools containing 96-single YAC clones from one master plate and further in super pools representing 10 master plates. This pooling strategy allows the selection of single YAC clones homologous to a target sequence after three rounds of PCR using super pools, single pools, and single YAC clones as a template. Single YAC clones were spheroplasted prior to the third PCR round in order to omit the conventional radioactive colony hybridization step. To date, we applied this PCR-based selection strategy to isolate clones homologousto ten different sequence-tagged sites (STS) that are linked to genes targeted for map-based cloning. The selection of YAC clones can be readily accomplished within three days. The PCR-based screening strategy is easy to set up and contributes to a further acceleration of the construction of YAC contigs.