Direct Haplotype‐Specific DNA Sequencing

Abstract

Determining haplotype‐specific DNA sequence information is very important in a wide range of research fields. However, no simple and robust approaches are currently available for determining haplotype‐specific sequence information. We have addressed this problem by developing a very simple and robust haplotype‐specific sequencing approach. We utilise the fact that DNA sequencing polymerases are sensitive to 3′end mismatches in the sequencing primer. By using two sequencing primers with 3′end corresponding to the two alleles in a given SNP locus, we are able to obtain allele‐specific DNA sequences from both alleles.

We evaluated this direct haplotype‐specific approach by determining haplotypes within the intron 2 sequence of the fructan‐6‐fructosyltransferase (6‐ft) gene in Lolium perenne L. We obtained reliable haplotype‐specific sequences for all primers and genotypes evaluated. We conclude that the haplotype‐specific sequencing is robust, and that the approach has a potentially very wide application range for any diploid organism.     

New perspectives in diet analysis based on DNA barcoding and parallel pyrosequencing: the trnL approach

The development of DNA barcoding (species identification using a standardized DNA sequence), and the availability of recent DNA sequencing techniques offer new possibilities in diet analysis. DNA fragments shorter than 100-150 bp remain in a much higher proportion in degraded DNA samples and can be recovered from faeces. As a consequence, by using universal primers that amplify a very short but informative DNA fragment, it is possible to reliably identify the plant taxon that has been eaten. According to our experience and using this identification system, about 50% of the taxa can be identified to species using the trnL approach, that is, using the P6 loop of the chloroplast trnL (UAA) intron. We demonstrated that this new method is fast, simple to implement, and very robust. It can be applied for diet analyses of a wide range of phytophagous species at large scales. We also demonstrated that our approach is efficient for mammals, birds, insects and molluscs. This method opens new perspectives in ecology, not only by allowing large-scale studies on diet, but also by enhancing studies on resource partitioning among competing species, and describing food webs in ecosystems.     

Shotgun haplotyping: a novel method for surveying allelic sequence variation

Haplotypic sequences contain significantly more information than genotypes of genetic markers and are critical for studying disease association and genome evolution. Current methods for obtaining haplotypic sequences require the physical separation of alleles before sequencing, are time consuming and are not scaleable for large surveys of genetic variation. We have developed a novel method for acquiring haplotypic sequences from long PCR products using simple, high-throughput techniques. This method applies modified shotgun sequencing protocols to sequence both alleles concurrently, with read-pair information allowing the two alleles to be separated during sequence assembly. Although the haplotypic sequences can be assembled manually from the resultant data using pre-existing sequence assembly software, we have devised a novel heuristic algorithm to automate assembly and remove human error. We validated the approach on two long PCR products amplified from the human genome and confirmed the accuracy of our sequences against full-length clones of the same alleles. This method presents a simple high-throughput means to obtain full haplotypic sequences potentially up to 20 kb in length and is suitable for surveying genetic variation even in poorly-characterized genomes as it requires no prior information on sequence variation.     

Specific-primer-directed DNA sequencing using automated fluorescence detection.

Automated fluorescence-based DNA sequence analysis offers the possibility to undertake very large scale sequencing projects. Directed strategies, such as the specific-primer-directed sequencing approach ('gene walking'), should prove useful in such projects. Described herein is a study involving the use of this approach in conjunction with automated fluorescence detection on a commercial instrument (ABI 370A DNA sequencer). This includes procedures for the rapid chemical synthesis and purification of labeled primers, the design of primer sequences that are compatible with the commercial analysis software, and automated DNA sequence analysis using such primers. A set of four fluorophore-labeled primers can be reliably synthesized in a twenty-four hour period, and greater than 300 nucleotides of analyzed new sequence obtained using this set in an additional twenty-four hours. Scale-up of these procedures to take advantage of the full capabilities of the sequencer is, at present, too slow and costly to be suitable for routine sequencing, and therefore the use of specific-primers is best suited to the closure of gaps in extended sequence produced using random cloning and sequencing strategies.     

Not all sequence tags are created equal: designing and validating sequence identification tags robust to indels

Ligating adapters with unique synthetic oligonucleotide sequences (sequence tags) onto individual DNA samples before massively parallel sequencing is a popular and efficient way to obtain sequence data from many individual samples. Tag sequences should be numerous and sufficiently different to ensure sequencing, replication, and oligonucleotide synthesis errors do not cause tags to be unrecoverable or confused. However, many design approaches only protect against substitution errors during sequencing and extant tag sets contain too few tag sequences. We developed an open-source software package to validate sequence tags for conformance to two distance metrics and design sequence tags robust to indel and substitution errors. We use this software package to evaluate several commercial and non-commercial sequence tag sets, design several large sets (max_count = 7,198) of edit metric sequence tags having different lengths and degrees of error correction, and integrate a subset of these edit metric tags to polymerase chain reaction (PCR) primers and sequencing adapters. We validate a subset of these edit metric tagged PCR primers and sequencing adapters by sequencing on several platforms and subsequent comparison to commercially available alternatives. We find that several commonly used sets of sequence tags or design methodologies used to produce sequence tags do not meet the minimum expectations of their underlying distance metric, and we find that PCR primers and sequencing adapters incorporating edit metric sequence tags designed by our software package perform as well as their commercial counterparts. We suggest that researchers evaluate sequence tags prior to use or evaluate tags that they have been using. The sequence tag sets we design improve on extant sets because they are large, valid across the set, and robust to the suite of substitution, insertion, and deletion errors affecting massively parallel sequencing workflows on all currently used platforms.     

Patterns of variation in MHC class II beta loci of the little greenbul (Andropadus virens) with comments on MHC evolution in birds

We have isolated major histocompatibility complex (MHC) class II beta loci from the little greenbul (Andropadus virens), an African songbird. We utilized preexisting information about conserved regions of the avian MHC to design primers to amplify a pool of sequences representing multiple loci. From this pool, a unique locus spanning 1109 bp that we designate as Anvi-DAB1 was cloned and sequenced. We designed locus-specific primers based on this sequence information and amplified six alleles from seven individuals. Compared to other A. virens MHC sequences obtained from genomic DNA or cDNA, the variability of sequences from Anvi-DAB1 was low and the ratio of nonsynonymous to synonymous substitution was much less than one, suggesting that Anvi-DAB1 may either be a pseudogene or a nonclassical MHC locus. Phylogenetic analysis revealed that the Anvi-DAB1 locus was highly divergent when compared with other passerine or A. virens genomic or transcribed MHC sequences. The use of conserved MHC primers followed by analysis of cloned sequences allows rapid isolation of MHC loci from exotic species and avoids laborious large-scale cloning and sequencing.     

Classification of Nucleotide Sequences Using Support Vector Machines

Species identification is one of the most important issues in biological studies. Due to recent increases in the amount of genomic information available and the development of DNA sequencing technologies, the applicability of using DNA sequences to identify species (commonly referred to as “DNA barcoding”) is being tested in many areas. Several methods have been suggested to identify species using DNA sequences, including similarity scores, analysis of phylogenetic and population genetic information, and detection of species-specific sequence patterns. Although these methods have demonstrated good performance under a range of circumstances, they also have limitations, as they are subject to loss of information, require intensive computation and are sensitive to model mis-specification, and can be difficult to evaluate in terms of the significance of identification. Here, we suggest a new DNA barcoding method in which support vector machine (SVM) procedures are adopted. Our new method is nonparametric and thus is expected to be robust for a wide range of evolutionary scenarios as well as multilocus analyses. Furthermore, we describe bootstrap procedures that can be used to test the significances of species identifications. We implemented a novel conversion technique for transforming sequence data to real-valued vectors, and therefore, bootstrap procedures can be easily combined with our SVM approach. In this study, we present the results of simulation studies and empirical data analyses to demonstrate the performance of our method and discuss its properties.     

Group-specific polymerase chain reaction for DNA-based analysis of species diversity and identity in dietary samples

Unique DNA sequences are present in all species and can be used as biomarkers for the detection of cells from that species. These DNA sequences can most easily be detected using the polymerase chain reaction (PCR), which allows very small quantities of target DNA sequence to be amplified even when the target is mixed with large amounts of nontarget DNA. PCR amplification of DNA markers that are present in a wide range of species has proven very useful for studies of species diversity in environmental samples. The taxonomic range of species to be identified from environmental samples may often need to be restricted to simplify downstream analyses and to ensure that less abundant sequences are amplified. Group-specific PCR primer sets are one means of specifying the range of taxa that produce an amplicon in a PCR. We have developed a range of group-specific PCR primers for studying the prey diversity found in predator stomach contents and scats. These primers, their design and their application to studying prey diversity and identity in predator diet are described.     

Clonal polymerase chain reaction-single-strand conformational polymorphism analysis: an effective approach for identifying cloned sequences

The amplification of target sequences from genomic DNA can result in more than one amplicon sequence being produced even when highly specific primers are used. Here we present a clonal polymerase chain reaction-single-strand conformational polymorphism (PCR-SSCP) approach for screening cloned amplicons and identifying particular clones prior to sequence determination. Comparison of the PCR-SSCP patterns of the cloned amplicons with the PCR-SSCP patterns observed for the DNA templates from which the clones were derived allows PCR artifacts, different alleles, and even different loci to be differentiated prior to sequencing. Using this approach, the number of clones required for reliable sequence determination is minimized, and complex “mixed” amplicons can be resolved easily, cost-effectively, and reliably.     

Development and characterization of simple sequence repeat (SSR) markers and their use in determining relationships among Lycopersicon esculentum cultivars

The simple sequence repeat (SSR) or microsatellite marker is currently the preferred molecular marker due to its highly desirable properties. The aim of this study was to develop and characterize more SSR markers because the number of SSR markers currently available in tomato is very limited. Five hundred DNA sequences of tomato were searched for SSRs and analyzed for the design of PCR primers. Of the 158 pairs of SSR primers screened against a set of 19 diverse tomato cultivars, 129 pairs produced the expected DNA fragments in their PCR products, and 65 of them were polymorphic with the polymorphism information content (PIC) ranging from 0.09 to 0.67. Among the polymorphic loci, 2-6 SSR alleles were detected for each locus with an average of 2.7 alleles per locus; 49.2% of these loci had two alleles and 33.8% had three alleles. The vast majority (93.8%) of the microsatellite loci contained di- or tri-nucleotide repeats and only 6.2% had tetra- and penta-nucleotide repeats. It was also found that TA/AT was the most frequent type of repeat, and the polymorphism information content (PIC) was positively correlated with the number of repeats. The set of 19 tomato cultivars were clustered based on the banding patterns generated by the 65 polymorphic SSR loci. Since the markers developed in this study are primarily from expressed sequences, they can be used not only for molecular mapping, cultivar identification and marker-assisted selection, but for identifying gene-trait relations in tomato.     

Solutions for PCR,cloning and sequencing errors in population genetic analysis

PCR and sequencing artefacts can seriously bias population genetic analyses, particularly of populations with low genetic variation such as endangered vertebrate populations. Here, we estimate the error rates, discuss their population genetics implications, and propose a simple detection method that helps to reduce the risk of accepting such errors. We study the major histocompatibility complex (MHC) class IIB of guppies, Poecilia reticulata and find that PCR base misincorporations inflate the apparent sequence diversity. When analysing neutral genes, such bias can inflate estimates of effective population size. Previously suggested protocols for identifying genuine alleles are unlikely to exclude all sequencing errors, or they ignore genuine sequence diversity. We present a novel and statistically robust method that reduces the likelihood of accepting PCR artefacts as genuine alleles, and which minimises the necessity of repeated genotyping. Our method identifies sequences that are unlikely to be a PCR artefact, and which need to be independently confirmed through additional PCR of the same template DNA. The proposed methods are recommended particularly for population genetic studies that involve multi-template DNA and in studies on genes with low genetic diversity.     

Characterizing homologues of crop domestication genes in poorly described wild relatives by high‐throughput sequencing of whole genomes

Wild crop relatives represent a source of novel alleles for crop genetic improvement. Screening biodiversity for useful or diverse gene homologues has often been based upon the amplification of targeted genes using available sequence information to design primers that amplify the target gene region across species. The crucial requirement of this approach is the presence of sequences with sufficient conservation across species to allow for the design of universal primers. This approach is often not successful with diverse organisms or highly variable genes. Massively parallel sequencing (MPS) can quickly produce large amounts of sequence data and provides a viable option for characterizing homologues of known genes in poorly described genomes. MPS of genomic DNA was used to obtain species‐specific sequence information for 18 rice genes related to domestication characteristics in a wild relative of rice, Microlaena stipoides. Species‐specific primers were available for 16 genes compared with 12 genes using the universal primer method. The use of species‐specific primers had the potential to cover 92% of the sequence of these genes, while traditional universal primers could only be designed to cover 80%. A total of 24 species‐specific primer pairs were used to amplify gene homologues, and 11 primer pairs were successful in capturing six gene homologues. The 23 million, 36‐base pair (bp) paired end reads, equated to an average of 2X genome coverage, facilitated the successful amplification and sequencing of six target gene homologues, illustrating an important approach to the discovery of useful genes in wild crop relatives.     

Genotyping HapSTR loci: phase determination from direct sequencing of PCR products

HapSTRs combine information from a microsatellite (or simple tandem repeat, STR) with one or more single nucleotide polymorphisms in the DNA sequence immediately flanking the STR. These loci may offer increased power for the estimation of demographic parameters, but also present some challenges for data collection and analysis. We describe a process for inferring HapSTR alleles, including the flanking haplotypes, STR alleles and their phase relative to each other, directly from DNA sequence electropherograms of PCR products from heterozygous individuals. Our approach eliminates the need for more costly and time-consuming processes, such as cloning or acrylamide gel electrophoresis to separate alleles prior to sequencing.     

A pragmatic approach to the analysis of diets of generalist predators: the use of next‐generation sequencing with no blocking probes

Predicting whether a predator is capable of affecting the dynamics of a prey species in the field implies the analysis of the complete diet of the predator, not simply rates of predation on a target taxon. Here, we employed the Ion Torrent next‐generation sequencing technology to investigate the diet of a generalist arthropod predator. A complete dietary analysis requires the use of general primers, but these will also amplify the predator unless suppressed using a blocking probe. However, blocking probes can potentially block other species, particularly if they are phylogenetically close. Here, we aimed to demonstrate that enough prey sequence could be obtained without blocking probes. In communities with many predators, this approach obviates the need to design and test numerous blocking primers, thus making analysis of complex community food webs a viable proposition. We applied this approach to the analysis of predation by the linyphiid spider Oedothorax fuscus in an arable field. We obtained over two million raw reads. After discarding the low‐quality and predator reads, the libraries still contained over 61 000 prey reads (3% of the raw reads; 6% of reads passing quality control). The libraries were rich in Collembola, Lepidoptera, Diptera and Nematoda. They also contained sequences derived from several spider species and from horticultural pests (aphids). Oedothorax fuscus is common in UK cereal fields, and the results showed that it is exploiting a wide range of prey. Next‐generation sequencing using general primers but without blocking probes provided ample sequences for analysis of the prey range of this spider and proved to be a simple and inexpensive approach.     

Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates.

A PCR-based approach to sequencing complete mitochondrial genomes is described along with a set of 86 primers designed primarily for avian mitochondrial DNA (mtDNA). This PCR-based approach allows an accurate determination of complete mtDNA sequences that is faster than sequencing cloned mtDNA. The primers are spaced at about 500-base intervals along both DNA strands. Many of the primers incorporate degenerate positions to accommodate variation in mtDNA sequence among avian taxa and to reduce the potential for preferential amplification of nuclear pseudogenes. Comparison with published vertebrate mtDNA sequences suggests that many of the primers will have broad taxonomic utility. In addition, these primers should make available a wider variety of mitochondrial genes for studies based on smaller data sets.     

SSCP is not so difficult: the application and utility of single-stranded conformation polymorphism in evolutionary biology and molecular ecology

All genetic markers are estimators of DNA nucleotide sequence variation. Rather than obtaining DNA sequence data, it is cheaper and faster to use techniques that estimate sequence variation, although this usually results in the loss of some information. SSCP (single-stranded conformation polymorphism) offers a sensitive but inexpensive, rapid, and convenient method for determining which DNA samples in a set differ in sequence, so that only an informative subset need be sequenced. In short, most DNA sequence variation can be detected with relatively little sequencing. SSCP has been widely applied in medical diagnosis, yet few studies have been published in population genetics. The utility and convenience of SSCP is far from fully appreciated by molecular population biologists. We hope to help redress this by illustrating the application of a single simple SSCP protocol to mitochondrial genes, nuclear introns, microsatellites, and anonymous nuclear sequences, in a range of vertebrates and invertebrates.     

Bidirectional sequencing of supercoiled plasmid DNA

In this paper we show that restriction DNA fragments can prime DNA synthesis of a homologous supercoiled plasmid DNA. Using the dideoxyribonucleotide chain terminator method, newly synthesized truncated chains can be detached from the primers by restriction enzyme digestion. Therefore, by choosing DNA fragments flanked by two different restriction enzymes sites, nucleotide sequence information can be simultaneously obtained on both regions of the DNA surrounding the restriction fragment. The advantage of this sequencing approach over current methods is that no prior knowledge of the primary sequence is needed to find the nucleotide sequence of a given DNA fragment. Thus, synthetic primers are not required and internal sequences of a given clone can be easily accessed without the need of fragmenting the original construct. The method has been used with rapid plasmid preparations, thus considerable time and effort can be saved in the gathering of nucleotide sequence information.     

大麦基因组中的微卫星标记及其应用

微卫星是以少数几个核苷酸为单位多次串联重复的DNA序列,是一种简单序列重复(simple sequence repeats,SSR),两侧一般是保守序列。由于它具有多态性高、共显性、容易用PCR检测和结果稳定可靠等特点,因此是一种十分理想的分子标记。大麦的微卫星DNA随机分布于基因组中,平均每一个微卫星基因座有3～18个等位基因,最高可达37个。SSR标记已广泛用于分子遗传图谱的构建、遗传多样性研究、种质鉴定、主要性状基因的定位及分子标记辅助选择育种等。大多数SSR标记集中在着丝粒附近区域,1HL、5HL和6HS明显缺乏SSR标记。大麦的SSR标记还有待进一步的开发。 Microsatellite Markers and Applications in the Barley Genome FENG Zong-yun1,2,3,ZHANG Yi-zheng1,LING Hong-qing3 1.College of Life Sciences,Sichuan University,Chengdu 610065,China; 2.College of Agronomy,Sichuan Agricultural University,Ya'an 625014,China; 3.The State Key Laboratory of Plant Cell & Chromosome Engineering,Institute of Genetics & Developmental Biology,Chinese Academy of Sciences,Beijing 100101,China Abstract:Microsatellites,also called simple sequence repeats (SSR),are simple,tandemly repeated DNA sequences with a repeat length of a few base pairs,and are very ideally used as molecular markers because of their abundance,high level of polymorphism,co-dominance and ease of assay with the polymerase chain reaction (PCR) by selecting primers as the conserved DNA sequences flanking the SSRs,as well as better stability.The experiments showed that SSRs are randomly distributed throughout the barley genome,and there are 3~18 alleles at a single SSR locus,up to 37 alleles/locus.SSR markers have being widely applied in the construction of molecular genetic map,the study of genetic diversity,the identification of germplasm,gene mapping for important traits and molecular marker-assisted selection.Meanwhile,most of markers are strongly clustered around the centromeric regions of all seven linkage groups.As a result of the clustering,genome coverage with SSRs remains incomplete with an obvious lack of markers on the long arms of chromosomes 1H and 5H and short arm of chromosome 6H.Therefore,it is very potential and necessary to further develop SSR markers in barley. Key words：barley;microsatellite marker;simple sequence repeats;genetic diversity;molecular mapping