Aberrant heteroplasmic transmission of mtDNA in cloned pigs arising from double nuclear transfer

Double nuclear transfer begins with the transfer of nuclear DNA from a donor cell into an enucleated recipient oocyte. This reconstructed oocyte is allowed to develop to the pronuclear stage, where the pronuclei are transferred into an enucleated zygote. This reconstructed zygote is then transferred to a surrogate sow. The genetic integrity of cloned offspring can be compromised by the transmission of mitochondrial DNA from the donor cell, the recipient oocyte and the recipient zygote. We have verified through the use of sequence analysis, restriction fragment length polymorphism analysis, allele specific PCR and primer extension polymorphism analysis that following double nuclear transfer the donor cell mtDNA is eliminated. However, it is likely that the recipient oocyte and zygote mitochondrial DNA are transmitted to the offspring, indicating bimaternal mitochondrial DNA transmission. This pattern of mtDNA inheritance is similar to that observed following cytoplasmic transfer and violates the strict unimaternal inheritance of mitochondrial DNA to offspring. This form of transmission raises concerns regarding the genetic integrity of cloned offspring and their uses in studies that require metabolic analysis or a stable genetic environment where only one variable is under analysis, such as in knockout technology.     

A subcloning strategy for DNA sequence analysis.

We describe here a new strategy of fragment preparation for sequencing procedures using endlabelled DNA fragments as substrates (2,3) which is directly applicable to DNA fragments cloned into the Pst I site of pBR322, or in modified form, to inserts into the BamH I or Sal I site of the same plasmid. Ordered sets of subclones of predetermined overlap are are generated. These can be sequenced directly without further strand- or fragment separation steps.     

Optical mapping of DNA: single-molecule-based methods for mapping genomes

The technologies associated with DNA sequencing are rapidly evolving. Indeed, single-molecule DNA sequencing strategies are cheaper and faster than ever before. Despite this progress, every sequencing platform to date relies on reading the genome in small, abstract fragments, typically of less than 1000 bases in length. The overarching aim of the optical map is to complement the information derived from DNA sequencing by providing long-range context on which these short sequence reads can be built. This is typically done using an enzyme to target and modify at short DNA sequences of, say, six bases in length throughout the genome. By accurately placing these short pieces of sequence on long genomic DNA fragments, up to several millions of bases in length, a scaffold for sequence assembly can be obtained. This review focuses on three enzymatic approaches to optical mapping. Optical mapping was first developed using restriction enzymes to sequence-specifically cleave DNA that is immobilized on a surface. More recently, nicking enzymes have found application in the sequence-specific fluorescent labeling of DNA for optical mapping. Such covalent modification allows the DNA to be imaged in solution, and this, in combination with developing nanofluidic technologies, is enabling new high-throughput approaches to mapping. And, finally, this review will discuss the recent development of mapping with subdiffraction-limit precision using methyltransferase enzymes to label the DNA with an ultrahigh density.     

PCR preparation of DNA inserts from lambda and plasmid vectors for RFLP mapping

A simplified method is described for preparing insert DNA for labelling reactions to be used in Southern hybridization. This method works with sequences cloned into both plasmid and lambda phage, and eliminates many of the steps leading to the labelling reaction. Small quantities of hostE. coli or lambda phage carrying a probe sequence are lysed and amplified via the polymerase chain reaction using standard sequencing primers. Unincorporated nucleotides are removed by ethanol precipitation or gel purification and insert DNA is ready for radio-labelling. This method reduces the time and expense associated with conventional insert preparation, and greatly simplifies the use of sequences cloned into lambda phage.     

Sequence alterations can mask each other's presence during screening with SSCP or heteroduplex analysis: BRCA genes as examples

For mutation detection, various screening techniques are widely used because DNA sequencing, the gold-standard method, is still considered to be expensive and laborious for high-throughput screening. Single-strand conformation polymorphism (SSCP) analysis, heteroduplex analysis (HA) and their variant techniques are popular and frequently used for this purpose. It is widely accepted that when searching for unknown sequence variations, any revealed distinct pattern should always be sequenced. We give examples here of the BRCA1 and BRCA2 genes where the SSCP/HA techniques can produce ambiguous predictions if used to detect known genetic variants compared to positive controls. Using direct DNA sequencing, we provide evidence that in such cases, mutations or polymorphisms can mask each other's presence. This phenomenon can often influence the results of any DNA testing because genetic variations such as single-nucleotide polymorphisms occur frequently in the human genome. We suggest that even in the case of known electrophoretic patterns of well-characterized genetic alterations, every sequence alteration should be confirmed by direct DNA sequencing, especially if genetic testing is carried out for diagnostic purposes.     

A parallel graph decomposition algorithm for DNA sequencing with nanopores

MOTIVATION: With the potential availability of nanopore devices that can sense the bases of translocating single-stranded DNA (ssDNA), it is likely that 'reads' of length approximately 10(5) will be available in large numbers and at high speed. We address the problem of complete DNA sequencing using such reads.We assume that approximately 10(2) copies of a DNA sequence are split into single strands that break into randomly sized pieces as they translocate the nanopore in arbitrary orientations. The nanopore senses and reports each individual base that passes through, but all information about orientation and complementarity of the ssDNA subsequences is lost. Random errors (both biological and transduction) in the reads create further complications. RESULTS: We have developed an algorithm that addresses these issues. It can be considered an extreme variation of the well-known Eulerian path approach. It searches over a space of de Bruijn graphs until it finds one in which (a) the impact of errors is eliminated and (b) both possible orientations of the two ssDNA sequences can be identified separately and unambiguously.Our algorithm is able to correctly reconstruct real DNA sequences of the order of 10(6) bases (e.g. the bacterium Mycoplasma pneumoniae) from simulated erroneous reads on a modest workstation in about 1 h. We describe, and give measured timings of, a parallel implementation of this algorithm on the Cray Multithreaded Architecture (MTA-2) supercomputer, whose architecture is ideally suited to this 'unstructured' problem. Our parallel implementation is crucial to the problem of rapidly sequencing long DNA sequences and also to the situation where multiple nanopores are used to obtain a high-bandwidth stream of reads.     

PCR-based accurate synthesis of long DNA sequences

Here we describe a simple and rapid method for assembly and PCR-based accurate synthesis (PAS) of long DNA sequences. The PAS protocol involves the following five steps: (i) design of the DNA sequence to be synthesized and of 60-bp overlapping oligonucleotides to cover the entire DNA sequence; (ii) purification of the oligonucleotides by PAGE; (iii) first PCR, to synthesize DNA fragments of 400-500 bp in length using 10 inner (template) and two outer (primer) oligonucleotides; (iv) second PCR, to assemble the products of the first PCR into the full-length DNA sequence; and (v) cloning and verification of the synthetic DNA by sequencing and, if needed, error correction using an overlap-extension PCR technique. This method, which takes approximately 1 wk, is suitable for synthesizing diverse types of long DNA molecule. We have successfully synthesized DNA fragments from 0.5 to 12.0 kb, with high G+C content, repetitive sequences or complex secondary structures. The PAS protocol therefore provides a simple, rapid, reliable and relatively inexpensive method for synthesizing long, accurate DNA sequences.     

The distribution of the frequency of occurrence of nucleotide subsequences, based on their overlap capability

DNA's genetic code can be represented as an alphabetic sequence composed of the four letters A, C, G, and T, which represent the four types of nucleotides--adenylic, cytidylic, guanylic, and thymidylic acid--of which DNA is composed. Now that these sequences have been identified for many genes and are available in computer-readable form, scientists can analyze these data and search for patterns in an attempt to learn more about the regulatory functions of the gene. One area of study is that of the frequency of occurrence of specific nucleotide subsequences (e.g., ACAC) within part or all of a nucleotide sequence. This paper derives the probability distribution of the frequency of occurrence of a subsequence within a nucleotide sequence, under the hypothesis that the four nucleotides occur at random and with equal probability. This distribution is nontrivial because different subsequences have different "overlap capability." For example, the subsequence AAAA can occur up to 17 times in a sequence of length 20 (which would happen if the sequence were composed solely of A's), but the subsequence ACGT cannot occur more than 5 times in a sequence of length 20. Thus, the frequency distributions are different for each type of overlap capability. It is of interest to assess and compare the degree of nonrandomness for different subsequences or among different portions of a sequence; the existence and degree of nonrandomness may be related to the type and degree of functionality of a nucleotide (sub)sequence. The frequency distributions provided here can be used to perform exact significance tests of the hypothesis of randomness. An approximate test is also described for use with long sequences; this can be used to test a more general null hypothesis of nucleotides occurring with unequal probabilities.     

Microevolution of the mitochondrial DNA control region in the Japanese brown bear (Ursus arctos) population.

We investigated nucleotide sequences of the mitochondrial DNA control region to describe natural genetic variations and to assess the relationships between subpopulations of the brown bear Ursus arctos on Hokkaido Island, Japan. Using the polymerase chain reaction product-direct sequencing technique, partial sequences (about 930 bases) of the control region were determined for 56 brown bears sampled throughout Hokkaido Island. A sequence alignment revealed that the brown bear control region included a variable sequence on the 5' side and a repetitive region on the 3' side. Phylogenetic trees reconstructed from the 5' variable region (696-702 bases) exhibited 17 haplotypes, which were clustered into three groups (Clusters A, B, and C). The distribution of each group did not overlap with those of the others, and the three different areas were located in separate mountainous forests of Hokkaido Island. Furthermore, most of the phylogenetically close haplotypes within each group were distributed geographically close to each other. In addition, the 3' repetitive region (arrays of 10 bases) exhibited a much faster mutation rate than the 5' variable region, resulting in heteroplasmy. Such mitochondrial DNA divergence in each group could have occurred after the brown bears migrated from the continent to Hokkaido and became fixed in the different areas.     

Microfluidic droplet enrichment for targeted sequencing

Targeted sequence enrichment enables better identification of genetic variation by providing increased sequencing coverage for genomic regions of interest. Here, we report the development of a new target enrichment technology that is highly differentiated from other approaches currently in use. Our method, MESA (Microfluidic droplet Enrichment for Sequence Analysis), isolates genomic DNA fragments in microfluidic droplets and performs TaqMan PCR reactions to identify droplets containing a desired target sequence. The TaqMan positive droplets are subsequently recovered via dielectrophoretic sorting, and the TaqMan amplicons are removed enzymatically prior to sequencing. We demonstrated the utility of this approach by generating an average 31.6-fold sequence enrichment across 250 kb of targeted genomic DNA from five unique genomic loci. Significantly, this enrichment enabled a more comprehensive identification of genetic polymorphisms within the targeted loci. MESA requires low amounts of input DNA, minimal prior locus sequence information and enriches the target region without PCR bias or artifacts. These features make it well suited for the study of genetic variation in a number of research and diagnostic applications.     

Pyrosequencing for Microbial Identification and Characterization

Pyrosequencing is a versatile technique that facilitates microbial genome sequencing that can be used to identify bacterial species, discriminate bacterial strains and detect genetic mutations that confer resistance to anti-microbial agents. The advantages of pyrosequencing for microbiology applications include rapid and reliable high-throughput screening and accurate identification of microbes and microbial genome mutations. Pyrosequencing involves sequencing of DNA by synthesizing the complementary strand a single base at a time, while determining the specific nucleotide being incorporated during the synthesis reaction. The reaction occurs on immobilized single stranded template DNA where the four deoxyribonucleotides (dNTP) are added sequentially and the unincorporated dNTPs are enzymatically degraded before addition of the next dNTP to the synthesis reaction. Detection of the specific base incorporated into the template is monitored by generation of chemiluminescent signals. The order of dNTPs that produce the chemiluminescent signals determines the DNA sequence of the template. The real-time sequencing capability of pyrosequencing technology enables rapid microbial identification in a single assay. In addition, the pyrosequencing instrument, can analyze the full genetic diversity of anti-microbial drug resistance, including typing of SNPs, point mutations, insertions, and deletions, as well as quantification of multiple gene copies that may occur in some anti-microbial resistance patterns.     

Sequencing by hybridization with the generic 6-mer oligonucleotide microarray: an advanced scheme for data processing

DNA sequencing by hybridization was carried out with a microarray of all 4(6) = 4,096 hexadeoxyribonucleotides (the generic microchip). The oligonucleotides immobilized in 100 x 100 x 20-microm polyacrylamide gel pads of the generic microchip were hybridized with fluorescently labeled ssDNA, providing perfect and mismatched duplexes. Melting curves were measured in parallel for all microchip duplexes with a fluorescence microscope equipped with CCD camera. This allowed us to discriminate the perfect duplexes formed by the oligonucleotides, which are complementary to the target DNA. The DNA sequence was reconstructed by overlapping the complementary oligonucleotide probes. We developed a data processing scheme to heighten the discrimination of perfect duplexes from mismatched ones. The procedure was united with a reconstruction of the DNA sequence. The scheme includes the proper definition of a discriminant signal, preprocessing, and the variational principle for the sequence indicator function. The effectiveness of the procedure was confirmed by sequencing, proofreading, and nucleotide polymorphism (mutation) analysis of 13 DNA fragments from 31 to 70 nucleotides long.     

Code domains in tandem repetitive DNA sequence structures

Traditionally, many people doing research in molecular biology attribute coding properties to a given DNA sequence if this sequence contains an open reading frame for translation into a sequence of amino acids. This protein coding capability of DNA was detected about 30 years ago. The underlying genetic code is highly conserved and present in every biological species studied so far. Today, it is obvious that DNA has a much larger coding potential for other important tasks. Apart from coding for specific RNA molecules such as rRNA, snRNA and tRNA molecules, specific structural and sequence patterns of the DNA chain itself express distinct codes for the regulation and expression of its genetic activity. A chromatin code has been defined for phasing of the histone-octamer protein complex in the nucleosome. A translation frame code has been shown to exist that determines correct triplet counting at the ribosome during protein synthesis. A loop code seems to organize the single stranded interaction of the nascent RNA chain with proteins during the splicing process, and a splicing code phases successive 5' and 3' splicing sites. Most of these DNA codes are not exclusively based on the primary DNA sequence itself, but also seem to include specific features of the corresponding higher order structures. Based on the view that these various DNA codes are genetically instructive for specific molecular interactions or processes, important in the nucleus during interphase and during cell division, the coding capability of tandem repetitive DNA sequences has recently been reconsidered.     

Optimal Sequencing Strategies for Surveying Molecular Genetic Diversity

Two commonly used measures of genetic diversity for intraspecies DNA sequence data are based, respectively, on the number of segregating sites, and on the average number of pairwise nucleotide differences. Expressions are derived for their variance in the presence of intragenic recombination for a panmictic population of fixed size that is at neutral equilibrium at the region sequenced. We show that, in contrast to the slow decrease in variance with increasing sample size, if the recombination rate is nonzero, the asymptotic rate of decrease of variance with increasing sequence length, for fixed sample size, is quite rapid. In particular, it is close to that which would be obtained by sequencing independent chromosome regions. The correlation between measures of diversity from linked regions is also examined. For a given total number of bases sequenced in a particular region, optimal sequencing strategies are derived. These typically involve sequencing relatively few (three to 10) long copies of the region. Under optimal strategies, the variances of the two measures are very similar for most parameter values considered. Results concerning optimal sequencing strategies will be sensitive to gross departures from the underlying assumptions, such as population bottlenecks, selective sweeps, and substantial population substructure.     

Advantages of genome sequencing by long-read sequencer using SMRT technology in medical area

PacBio RS II is the first commercialized third-generation DNA sequencer able to sequence a single molecule DNA in real-time without amplification. PacBio RS II’s sequencing technology is novel and unique, enabling the direct observation of DNA synthesis by DNA polymerase. PacBio RS II confers four major advantages compared to other sequencing technologies: long read lengths, high consensus accuracy, a low degree of bias, and simultaneous capability of epigenetic characterization. These advantages surmount the obstacle of sequencing genomic regions such as high/low G+C, tandem repeat, and interspersed repeat regions. Moreover, PacBio RS II is ideal for whole genome sequencing, targeted sequencing, complex population analysis, RNA sequencing, and epigenetics characterization. With PacBio RS II, we have sequenced and analyzed the genomes of many species, from viruses to humans. Herein, we summarize and review some of our key genome sequencing projects, including full-length viral sequencing, complete bacterial genome and almost-complete plant genome assemblies, and long amplicon sequencing of a disease-associated gene region. We believe that PacBio RS II is not only an effective tool for use in the basic biological sciences but also in the medical/clinical setting.     

Assembling millions of short DNA sequences using SSAKE

Novel DNA sequencing technologies with the potential for up to three orders magnitude more sequence throughput than conventional Sanger sequencing are emerging. The instrument now available from Solexa Ltd, produces millions of short DNA sequences of 25 nt each. Due to ubiquitous repeats in large genomes and the inability of short sequences to uniquely and unambiguously characterize them, the short read length limits applicability for de novo sequencing. However, given the sequencing depth and the throughput of this instrument, stringent assembly of highly identical sequences can be achieved. We describe SSAKE, a tool for aggressively assembling millions of short nucleotide sequences by progressively searching through a prefix tree for the longest possible overlap between any two sequences. SSAKE is designed to help leverage the information from short sequence reads by stringently assembling them into contiguous sequences that can be used to characterize novel sequencing targets. Availability: http://www.bcgsc.ca/bioinfo/software/ssake.     

System for DNA sequencing with resolution of up to 600 base pairs

A system capable of resolving about 500 bases is of interest for sequencing of longer DNA molecules. Studies on further optimization of resolution on DNA sequencing gels were carried out. The effect of physico-chemical properties of gels and buffers on resolution were tested, e.g. ionic strength and pH of buffers, different buffer systems, acrylamide concentration, crosslinker concentration, type of crosslinker, temperature of polymerization, denaturing conditions, gel length and thickness. Tested were as well different running conditions like electric field, gel temperature, dimension of sample slots. Gels 0.1-0.2 mm thick and up to 1.2 m long were cast and tested routinely. Gel lengths of 60-70 cm (for sequencing up to 350-400 bases) to about 100 cm (above 400 bases) are practicable. Little is gained in resolution by increasing the gel length from 1 to 1.2 m. Resolution was improved using 0.1 mm thick gels, at a higher pH value of 8.6-8.8, and molarity increased to 0.2 M. The sequencing pattern in the region of higher bases could be better resolved on a twice-magnified picture of that region on the autoradiogram. With the long gels (70-120 cm), it is advantageous to obtain the sequence overlap by running in parallel gels of different concentrations, without re-application of samples, all loaded at the same time. Buffer chamber for running of two of three gels and thermostating plates up to 1.2 m long were designed. In this way four to six thermostated gels can be run from a power supply with two inputs. Three 1 m long gels (concentrations: 4%, 6%, 12-16%) are loaded with several samples of DNA to be sequenced and run in parallel without re-application of the samples. With good samples, the sequence overlap from the gels could be counted up to 500 base pairs, with exceptionally good samples closer to 600 bases. At present this number seems to be near the limit of the resolving power of the polyacrylamide gels.     

Some possible codes for encrypting data in DNA

Three codes are reported for storing written information in DNA. We refer to these codes as the Huffman code, the comma code and the alternating code. The Huffman code was devised using Huffman's algorithm for constructing economical codes. The comma code uses a single base to punctuate the message, creating an automatic reading frame and DNA which is obviously artificial. The alternating code comprises an alternating sequence of purines and pyrimidines, again creating DNA that is clearly artificial. The Huffman code would be useful for routine, short-term storage purposes, supposing – not unrealistically – that very fast methods for assembling and sequencing large pieces of DNA can be developed. The other two codes would be better suited to archiving data over long periods of time (hundreds to thousands of years).     

HITAC-seq enables high-throughput cost-effective sequencing of plasmids and DNA fragments with identity

DNA sequencing is vital for many aspects of biological research and diagnostics. Despite the development of second and third generation sequencing technologies, Sanger sequencing has long been the only choice when required to precisely track each sequenced plasmids or DNA fragments. Here, we report a complete set of novel barcoding and assembling system, Highly-parallel Indexed Tagmentation-reads Assembled Consensus sequencing(HITAC-seq), that could massively sequence and track the identities of each individual sequencing sample. With the cost of much less than that of single read of Sanger sequencing,HITAC-seq can generate high-quality contiguous sequences of up to 10 kilobases or longer. The capability of HITAC-seq was confirmed through large-scale sequencing of thousands of plasmid clones and hundreds of amplicon fragments using approximately 100 pg of input DNAs. Due to its long synthetic length, HITACseq was effective in detecting relatively large structural variations, as demonstrated by the identification of a～1.3 kb Copia retrotransposon insertion in the upstream of a likely maize domestication gene. Besides being a practical alternative to traditional Sanger sequencing, HITAC-seq is suitable for many highthroughput sequencing and genotyping applications.