Genotyping for Species Identification and Diversity Assessment Using Double-Digest Restriction Site-Associated DNA Sequencing (ddRAD-Seq).

Genotyping-by-sequencing (GBS) is a powerful approach for studying the genetic diversity of legume species. By using restriction enzymes or other methods to generate a reduced representation of the genome for sequencing, GBS can provide genome-wide single nucleotide polymorphisms (SNP) for diversity analysis at high throughput and low cost. Here we describe a novel double-digest restriction site-associated DNA sequencing (ddRAD-seq) approach. We also describe the downstream bioinformatic analysis of the sequencing data, including alignment to a reference genome, de novo assembly, SNP calling, phylogenetic analysis, and structure analysis.

Development of a universal double-digest RAD sequencing approach for a group of nonmodel, ecologically and economically important insect and fish taxa.

The generation of genome-scale data is critical for a wide range of questions in basic biology using model organisms, but also in questions of applied biology in nonmodel organisms (agriculture, natural resources, conservation and public health biology). Using a genome-scale approach on a diverse group of nonmodel organisms and with the goal of lowering costs of the method, we modified a multiplexed, high-throughput genomic scan technique utilizing two restriction enzymes. We analysed several pairs of restriction enzymes and completed double-digestion RAD sequencing libraries for nine different species and five genera of insects and fish. We found one particular enzyme pair produced consistently higher number of sequence-able fragments across all nine species. Building libraries off this enzyme pair, we found a range of usable SNPs between 4000 and 37 000 SNPS per species and we found a greater number of usable SNPs using reference genomes than de novo pipelines in STACKS. We also found fewer reads in the Read 2 fragments from the paired-end Illumina Hiseq run. Overall, the results of this study provide empirical evidence of the utility of this method for producing consistent data for diverse nonmodel species and suggest specific considerations for sequencing analysis strategies.

Application of Created Restriction Site PCR-RFLP to Identify POT1 Gene Polymorphism.

Protection of telomeres protein 1 (POT1) plays pivotal roles in protection of chromosome ends and regulation of telomere length with other telomere binding proteins; its genetic polymorphisms are associated with many diseases. In this study, we explored a novel PCR-RFLP method for typing the single nucleotide polymorphism (SNP) rs1034794 of the human POT1 gene. A new restriction enzyme site was introduced into a POT1 gene amplification product by created restriction site PCR (CRS-PCR). One primer was designed based on changed sequence; after PCR amplification, a new restriction enzyme site for AluI was introduced into the PCR products. One hundred and seventy eight samples from Han Chinese individuals were tested to evaluate this new method. The 3'-end of the forward primer was next to the polymorphic site, and the third base from the 3'-end was the mismatched base A. The final PCR product contained the AGCT sequence (AluI recognition site) when the ancestral POT1 alleles were amplified. The data obtained with the new method perfectly matched those obtained with the sequencing method. Thus, CRS-PCR is a new low-cost and high-efficiency alternative for rs1034794 typing.

Enhanced reduced representation bisulfite sequencing for assessment of DNA methylation at base pair resolution.

DNA methylation pattern mapping is heavily studied in normal and diseased tissues. A variety of methods have been established to interrogate the cytosine methylation patterns in cells. Reduced representation of whole genome bisulfite sequencing was developed to detect quantitative base pair resolution cytosine methylation patterns at GC-rich genomic loci. This is accomplished by combining the use of a restriction enzyme followed by bisulfite conversion. Enhanced Reduced Representation Bisulfite Sequencing (ERRBS) increases the biologically relevant genomic loci covered and has been used to profile cytosine methylation in DNA from human, mouse and other organisms. ERRBS initiates with restriction enzyme digestion of DNA to generate low molecular weight fragments for use in library preparation. These fragments are subjected to standard library construction for next generation sequencing. Bisulfite conversion of unmethylated cytosines prior to the final amplification step allows for quantitative base resolution of cytosine methylation levels in covered genomic loci. The protocol can be completed within four days. Despite low complexity in the first three bases sequenced, ERRBS libraries yield high quality data when using a designated sequencing control lane. Mapping and bioinformatics analysis is then performed and yields data that can be easily integrated with a variety of genome-wide platforms. ERRBS can utilize small input material quantities making it feasible to process human clinical samples and applicable in a range of research applications. The video produced demonstrates critical steps of the ERRBS protocol.

A simple strain typing assay for Trypanosoma cruzi: discrimination of major evolutionary lineages from a single amplification product.

BACKGROUND: Trypanosoma cruzi is the causative agent of Chagas' Disease. The parasite has a complex population structure, with six major evolutionary lineages, some of which have apparently resulted from ancestral hybridization events. Because there are important biological differences between these lineages, strain typing methods are essential to study the T. cruzi species. Currently, there are a number of typing methods available for T. cruzi, each with its own advantages and disadvantages. However, most of these methods are based on the amplification of a variable number of loci. METHODOLOGY/PRINCIPAL FINDINGS: We present a simple typing assay for T. cruzi, based on the amplification of a single polymorphic locus: the TcSC5D gene. When analyzing sequences from this gene (a putative lathosterol/episterol oxidase) we observed a number of interesting polymorphic sites, including 1 tetra-allelic, and a number of informative tri- and bi-allelic SNPs. Furthermore, some of these SNPs were located within the recognition sequences of two commercially available restriction enzymes. A double digestion with these enzymes generates a unique restriction pattern that allows a simple classification of strains in six major groups, corresponding to DTUs TcI-TcIV, the recently proposed Tcbat lineage, and TcV/TcVI (as a group). Direct sequencing of the amplicon allows the classification of strains into seven groups, including the six currently recognized evolutionary lineages, by analyzing only a few discriminant polymorphic sites. CONCLUSIONS/SIGNIFICANCE: Based on these findings we propose a simple typing assay for T. cruzi that requires a single PCR amplification followed either by restriction fragment length polymorphism analysis, or direct sequencing. In the panel of strains tested, the sequencing-based method displays equivalent inter-lineage resolution to recent multi- locus sequence typing assays. Due to their simplicity and low cost, the proposed assays represent a good alternative to rapidly screen strain collections, providing the cornerstone for the development of robust typing strategies.

In silico fingerprinting (ISIF): a user-friendly in silico AFLP program.

The Amplified fragment Length Polymorphism (AFLP) is one of the cost-effective and useful fingerprinting techniques to study non-model species. One crucial AFLP step in the AFLP procedure is the choice of restriction enzymes and selective bases providing good-quality AFLP profiles. Here, we present a user-friendly program (ISIF) that allows carrying out in silico AFLPs on species for which whole genome sequences are available. Carrying out in silico analyses as preliminary tests can help to optimize the experimental work by allowing a rapid screening of candidate restriction enzymes and the combinations of selective bases to be used. Furthermore, using in silico AFLPs is of great interest to limit homoplasy and amplification of repetitive elements to target genomic regions of interest or to optimize complex and costly high-throughput genomic experiments.

Quantitatively evaluating mistaken clone assignments by RFLP analysis of 16S rRNA genes: a case study.

We quantitatively evaluated the errors of clone assignment based on the restriction fragment length polymorphism (RFLP) pattern of 16S rRNA genes. Eighty clones were randomly selected from a 16S rRNA gene library and were categorized into 35 operational taxonomic units (OTU) based on their indistinguishable enzyme restriction patterns of 3 tetrameric restriction enzymes RsaI, BsuRI, and HinfI. All of these clones were then sequenced and were reassigned into 36-53 OTUs using the DOTUR program when sequence similarities of 95%-100% were used. The number of the identically assigned clones ranged from 53 to 61 and the percentage varied from 66.3% to 76.3%. The Shannon-Weaver index for the bacterial community observed by RFLP analysis was 2.75, equal to that estimated by DOTUR at a 97% sequence similarity. Compared with clones assigned with the DOTUR program at a 97% sequence similarity, only 61 clones (76.3%) were correctly assigned by RFLP analysis. Six clones (7.5%) were assigned mistakenly at the phylum level, and the positions of 13 clones (16.2%) were phylogenetically different at a lower taxonomic rank.

Effect of restriction endonucleases on assessment of biodiversity of cultivable polar marine planktonic bacteria by amplified ribosomal DNA restriction analysis.

To choose a suitable restriction endonuclease for quick assessment of bacterial diversity in polar environments by ARDRA, we investigated the effect of restriction enzymes on ARDRA patterns of cultivable marine planktonic bacteria isolated from polar region. Thirty-three isolates were analyzed by ARDRA using five enzymes (HinfI, HaeIII, AluI, and the mix AfaI/MspI), respectively, resulting in different groups, each group corresponding to a particular genotype. A comparison of the ARDRA patterns was carried out, and phylogenetic position of all thirty-three bacteria was obtained by 16S rDNA sequencing. Consistent with phylogenetic analysis, ARDRA pattern comparison revealed that AluI, being sensitive and reliable enough to generate species-specific patterns, was a suitable restriction enzyme used for evaluating bacterial diversity, suggesting a combination of ARDRA with AluI and 16S rDNA sequencing can provide a simple, fast and reliable means for bacterial identification and diversity assessment in polar environments.

Improved assessment of denitrifying, N2-fixing, and total-community bacteria by terminal restriction fragment length polymorphism analysis using multiple restriction enzymes.

A database of terminal restriction fragments (tRFs) of the 16S rRNA gene was set up utilizing 13 restriction enzymes and 17,327 GenBank sequences. A computer program, termed TReFID, was developed to allow identification of any of these 17,327 sequences by means of polygons generated from the specific tRFs of each bacterium. The TReFID program complements and exceeds in its data content the Web-based phylogenetic assignment tool recently described by A. D. Kent, D. J. Smith, B. J. Benson, and E. W. Triplett (Appl. Environ. Microb. 69:6768-6766, 2003). The method to identify bacteria is different, as is the region of the 16S rRNA gene employed in the present program. For the present communication the software of the tRF profiles has also been extended to allow screening for genes coding for N2 fixation (nifH) and denitrification (nosZ) in any bacterium or environmental sample. A number of controls were performed to test the reliability of the TReFID program. Furthermore, the TReFID program has been shown to permit the analysis of the bacterial population structure of bacteria by means of their 16S rRNA, nifH, and nosZ gene content in an environmental habitat, as exemplified for a sample from a forest soil. The use of the TReFID program reveals that noncultured denitrifying and dinitrogen-fixing bacteria might play a more dominant role in soils than believed hitherto.

BACFinder: genomic localisation of large insert genomic clones based on restriction fingerprinting.

We have developed software that allows the prediction of the genomic location of a bacterial artificial chromosome (BAC) clone, or other large genomic clone, based on a simple restriction digest of the BAC. The mapping is performed by comparing the experimentally derived restriction digest of the BAC DNA with a virtual restriction digest of the whole genome sequence. Our trials indicate that this program identified the genomic regions represented by BAC clones with a degree of accuracy comparable to that of end-sequencing, but at considerably less cost. Although the program has been developed principally for use with Arabidopsis BACs, it should align large insert genomic clones to any fully sequenced genome.

[Diagnostic and differential diagnostic potential of mitochondrial DNA assessment in patients with Leber's hereditary optic neuropathy].

OBJECTIVE: To study the primary mutations of mitochondrial DNA (mtDNA) associated with Leber's hereditary optic neuropathy (LHON) in patients with optic neuropathy. METHODS: Seventy-nine patients with a variety of bilateral optic neuropathy were examined. Mutations at np 3,460, np 11,778 and np 14,484 of mtDNA were tested by PCR-restriction fragment length polymorphism technique to detect DNA in peripheral blood. The samples were taken from 16 cases of clinically diagnosed LHON, 44 cases of suspected LHON, two cases of alcohol amblyopia, four cases of multiple sclerosis, five cases of autosomal dominant hereditary optic atrophy, 4 cases with primary open-angle glaucoma, three cases of spinocerebellar degeneration, and one case of ethambutol-induced optic neuropathy. RESULTS: The mutation at np 11,778 was identified in 31 cases (39.2%), consisting of all the 16 clinically diagnosed LHON cases, thirteen cases (29.5%) of the suspected LHON, and the two cases of alcohol amblyopia. The remaining 48 cases were negative for mtDNA mutations at np 3,460, np 11,778, or np 14,484. CONCLUSION: Assessment of mtDNA provides a useful diagnostic aid in confirming and excluding the diagnosis of LHON, particularly useful in cases without a family hereditary history and cases with cause unknown bilateral optic neuritis.

Extended stability of restriction enzymes at ambient temperatures.

The stability of restriction enzymes as supplied by manufacturers without any modification has been examined. No reduction in activity was observed for three enzymes (HindIII, EcoRI and Tsp509I) held at ambient temperature or 4 degrees C for the period of study (12 months), while activity was observed for up to 12 weeks after storage at 37 degrees C, which was considerably better than following desiccation with trehalose, a recognized preservation technique. A larger trial of 23 different restriction enzymes held at room temperature for one week showed that all enzymes retained significant activity. As a practical demonstration of the usefulness of this finding, enzymes were posted to Africa by conventional mail (cost $1 US) and shown to retain activity upon arrival after three weeks in transit (compared to a cost of $1000 US by cold-chain transportation). Supplying enzymes to third-world markets should now be possible by removing the necessity for cold-chain transport. After arrival, enzymes can simply be stored in a standard domestic refrigerator.

A variation of the amplified-fragment length polymorphism (AFLP) technique using three restriction endonucleases, and assessment of the enzyme combination BglII-MfeI for AFLP analysis of Salmonella enterica subsp. enterica isolates.

We have performed amplified-fragment length polymorphism (AFLP) fingerprinting on a collection of Salmonella enterica subsp. enterica serovar typhimurium strains with a restriction endonuclease combination (BglII and MfeI) that has previously been used successfully for typing Campylobacter jejuni isolates with high resolution. Additionally, a variation of the AFLP assay in which two rare cutting restriction enzymes (XbaI and BsrGI) in combination with the frequent cutter (HinP1I) was examined. The BglII and MfeI enzyme combination offered low resolution for genotyping Salmonella typhimurium isolates and is not recommended for this common serovar. The three-enzyme combination gave a higher discrimination, and is thus a new alternate way of performing AFLP fingerprinting of S. typhimurium.

Optimization of restriction fragment DNA mapping.

Consider a mapping project in which overlap of clonal segments is inferred from complete multiple restriction digests. The fragment sizes of the clones are measured with some error, potentially leading to a map with erroneous links. The number of errors in the map depends on the number and types of enzymes used to characterize the clones. The most critical parameter is the decision rule k, or the criterion for declaring clone overlap. Small changes in k may cause an order of magnitude change in the amount of work it takes to build a map of given completion. We observe that the cost of an optimal mapping strategy is approximately proportional to the target size. While this finding is encouraging, considerable effort is nonetheless required: for large-scale sequencing projects with up-front mapping, mapping will be a non-negligible fraction of the total sequencing cost.

Assessment by Southern blot analysis of UV-induced damage and repair in human immunoglobulin genes.

Irradiation of DNA with UV light induces pyrimidine dimers and (6-4) photoproducts. The presence of one of these photolesions in the restriction site of a given endonuclease inhibits DNA cleavage and induces the formation of fragments by incomplete DNA digestion which appear as additional, facultative bands in Southern hybridization autoradiograms. The number and size of these fragments show a positive correlation with the UV dose. The response to UV light of immunoglobulin light-chain constant kappa and heavy-chain constant mu genes was analyzed with 2 specific probes. Constant kappa and mu genes when irradiated as part of the chromatin of living lymphocytes showed a UV sensitivity similar to that of naked DNA. The same genes from granulocytes had 50-60 times lower UV sensitivity. When cells were allowed to repair photolesions for 24 h the facultative bands from granulocytes disappeared indicating that these cells were able to remove photolesions from constant kappa and mu genes. Facultative bands from lymphocytes showed a smaller decrease of density after 24 h repair. This suggests that lymphocytes are less efficient than granulocytes in removing UV damage from constant kappa and mu genes.

An economic "power supply" using a diode for agarose and polyacrylamide gel electrophoresis.

Gel electrophoresis using agarose and polyacrylamide is a technique widely used for the separation of DNA and protein of various molecular sizes. This type of electrophoresis involves use of a rather expensive dc power supply composed of a complicated circuit. A simple circuit consisting of a single diode was designed, which can be used as a substitute for the conventional power supply. There are no appreciable differences in the so-obtained electrophoretic patterns. This simple, economic system is expected to be applied for use in laboratory, industrial, educational, and clinical facilities.