BiSearch: primer-design and search tool for PCR on bisulfite-treated genomes

Bisulfite genomic sequencing is the most widely used technique to analyze the 5-methylation of cytosines, the prevalent covalent DNA modification in mammals. The process is based on the selective transformation of unmethylated cytosines to uridines. Then, the investigated genomic regions are PCR amplified, subcloned and sequenced. During sequencing, the initially unmethylated cytosines are detected as thymines. The efficacy of bisulfite PCR is generally low; mispriming and non-specific amplification often occurs due to the T richness of the target sequences. In order to ameliorate the efficiency of PCR, we developed a new primer-design software called BiSearch, available on the World Wide Web. It has the unique property of analyzing the primer pairs for mispriming sites on the bisulfite-treated genome and determines potential non-specific amplification products with a new search algorithm. The options of primer-design and analysis for mispriming sites can be used sequentially or separately, both on bisulfite-treated and untreated sequences. In silico and in vitro tests of the software suggest that new PCR strategies may increase the efficiency of the amplification.     

Characterization of an <Emphasis Type="Italic">Eco</Emphasis>RI family of satellite DNA from two species

We have cloned and sequenced a 321bp band of repetitive DNA from Eptesicus fuscus and E. serotinus observed after gel electrophoresis of EcoRI digested genomic DNA in both species. Southern blot analysis of genomic DNA (from both species) digested with the same enzyme showed the existence of a ladder pattern indicating that the repetitive DNA is arrayed in tandem. The repetitive sequences have a monomer unit of 321bp which is composed of two subunits of 160bp, suggested by the existence of a 160bp band in the ladder of E. fuscus and by the presence of some direct repeats found in the analysis of the consensus sequence. Analysis of the methylation status demonstrated that cytosines in CCGG sequences in this satellite DNA are methylated in E. fuscus but not in the E. serotinus. Alignment of the sequenced clones showed that several nucleotide positions are diagnostic species-specific and consequently the phylogenetic analysis grouped the monomer units from both species in two clearly separated groups.     

Extraction and purification of microbial DNA from petroleum-contaminated soils and detection of low numbers of toluene, octane and pesticide degraders by multiplex polymerase chain reaction and Southern analysis

We investigated the use of multiplex polymerase chain reaction (FCR) techniques coupled with Southern analysis to detect xenobiotic-degrading organisms that had been added to three soils. Two soils highly contaminated with petroleum hydrocarbons and a less contaminated control soil were amended with tenfold dilutions of Pseudomonas putida mt-2 (pWWO), P. oleovorans (OCT), and Alcaligenes eutrophus JMP134 (pJP4), or, for controls, phosphate buffer alone. Total DNA was then isolated from the soils and purified using a sequential precipitation and dissolution purification procedure. This DNA was subjected to multiplex polymerase chain reaction (PCR) using primers that amplify regions of xylM (PCR product = 631 bp), alkB (546 bp) and tfdA (710 bp), which are found on pWWO, OCT and pJP4, respectively. The sizes of the amplified DNA fragments were designed to permit simultaneous amplification and detection of the target genes. Ethidium bromide-stained gels of the initial PCR reaction indicated detectable amplification of between 10* to 10* cells per gram soil, depending on the soil and the target gene. Southern analysis of the PCR amplified DNA improved detection limits to between 1 and 10 cells of each target species per gram of soil, and confirmed the identity of the PCR products. For some samples that were initially resistant to PCR, dilution of the environmental DNA resulted in positive PCR results. This treatment presumably overcame the inhibition of the PCR by diluting coextracted contaminants in the environmental DNA. A second PCR on an aliquot (1 μL) of the first reaction increased the ethidium bromide-based detection limits for one of the soils to six cells per gram of soil; it did not increase the detection limits for the other soils. Therefore, the DNA extraction procedure and multiplex PCR permitted the simultaneous detection of three types of biodegradarJve cells, at a lower detection limit of = > 10 cells per gram of highly contaminated, organic soil. However, due to kinetic limitations of multiplex PCR, the amplified signals did not follow a close dose response to the numbers of added target cells.     

A taxon-specific oligonucleotide probe for temperate zone soil isolates of Glomus mosseae

 The 5.8 S subunit and flanking internal transcribed spacer (ITS) regions in nuclear ribosomal DNA (rDNA) from spores of Glomus mosseae FL156 and UK118 were amplified by polymerase chain reaction (PCR) using ITS1 and ITS4 as primers. The amplification product from template DNA of UK118 was cloned and sequenced (569 bp); the amplified DNA from FL156 was sequenced directly (582 bp). There was a 95% sequence similarity between DNAs amplified from the two isolates; in contrast, major dissimilarities with partial sequences of seven other glomalean taxa were observed. Four oligonucleotide sequences unique to Glomus mosseae were identified as potential primers. Their specificity to Glomus mosseae was assessed by PCR amplification of genomic DNA from spores from 36 glomalean fungi: 13 isolates of Glomus mosseae, two Glomus monosporum, 10 other Glomus isolates, and 11 other glomalean taxa from each of four other genera. The Glomus mosseae isolates were from a broad range of temperate zone agricultural soils. Oligonucleotide pair GMOS1 : GMOS2 primed specific amplification of an oligonucleotide sequence (approximately 400 bp) present in all Glomus mosseae isolates and two isolates of the closely related Glomus monosporum. This primer pair did not prime PCR when the template consisted of DNA from any of the other glomalean fungi or any of the nonmycorrhizal controls. In addition, a 24-mer oligonucleotide, designated GMOS5, hybridized with Glomus mosseae and Glomus monosporum DNA amplified by PCR using primer pairs ITS1 : ITS4 and GMOS1 : GMOS2. Colony-blot assays showed that GMOS5 hybridized to 100% and 97% of E. coli pUC19 clones of amplification products from Glomus mosseae FL156 and UK118 DNA templates, respectively, indicating that nearly all clones contained an homologous sequence. GMOS5 was used successfully to detect specifically Glomus mosseae in DNA extracted from colonized sudan grass (Sorghum sudanense L.) roots and amplified by PCR using the primer pair GMOS1 : GMOS2. The results confirm several previous indications that Glomus mosseae and Glomus monosporum are indistinguishable taxonomic entities. Accepted: 14 February 1998     

Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis

We describe a large-scale random approach termed reduced representation bisulfite sequencing (RRBS) for analyzing and comparing genomic methylation patterns. BglII restriction fragments were size-selected to 500–600 bp, equipped with adapters, treated with bisulfite, PCR amplified, cloned and sequenced. We constructed RRBS libraries from murine ES cells and from ES cells lacking DNA methyltransferases Dnmt3a and 3b and with knocked-down (kd) levels of Dnmt1 (Dnmt[1^kd,3a^−/−,3b^−/−]). Sequencing of 960 RRBS clones from Dnmt[1^kd,3a^−/−,3b^−/−] cells generated 343 kb of non-redundant bisulfite sequence covering 66212 cytosines in the genome. All but 38 cytosines had been converted to uracil indicating a conversion rate of >99.9%. Of the remaining cytosines 35 were found in CpG and 3 in CpT dinucleotides. Non-CpG methylation was >250-fold reduced compared with wild-type ES cells, consistent with a role for Dnmt3a and/or Dnmt3b in CpA and CpT methylation. Closer inspection revealed neither a consensus sequence around the methylated sites nor evidence for clustering of residual methylation in the genome. Our findings indicate random loss rather than specific maintenance of methylation in Dnmt[1^kd,3a^−/−,3b^−/−] cells. Near-complete bisulfite conversion and largely unbiased representation of RRBS libraries suggest that random shotgun bisulfite sequencing can be scaled to a genome-wide approach.     

Sequence comparison of distal and proximal ribosomal DNA arrays in rice (Oryza sativa L.) chromosome 9S and analysis of their flanking regions

Rice (Oryza sativa ssp. japonica cv. Nipponbare) harbors a ribosomal RNA gene (rDNA) cluster in the nucleolar-organizing region at the telomeric end of the short arm of chromosome 9. We isolated and sequenced two genomic clones carrying rice rDNA fragments from this region. The rice rDNA repeat units could be classified into three types based on length, which ranged from 7,928 to 8,934 bp. This variation was due to polymorphism in the number of 254-bp subrepeats in the intergenic spacer (IGS). Polymerase chain reaction (PCR) analysis suggested that the rDNA units in rice vary widely in length and that the copy number of the subrepeats in the IGS ranges from 1 to 12 in the rice genome. PCR and Southern blot analyses showed that most rDNA units have three intact and one truncated copies of the subrepeats in the IGS, and distal (telomere-side) rDNA units have more subrepeats than do proximal (centromere-side) ones. Both genomic clones we studied contained rDNA-flanking DNA sequences of either telomeric repeats (5′-TTTAGGG-3′) or a chromosome-specific region, suggesting that they were derived from the distal or proximal end, respectively, of the rDNA cluster. A similarity search indicated that retrotransposons appeared more frequently in a 500-kb portion of the proximal rDNA-flanking region than in other subtelomeric regions or sequenced regions of the genome. This study reveals the repetitive nature of the telomeric end of the short arm of chromosome 9, which consists of telomeric repeats, an rDNA array, and a retrotransposon-rich chromosomal region.Sequence accession numbers in DDBJ assigned for OSJNOa063K24 and OSJNBb0013K10 are AP009051 and AP008245, respectively.     

Detection of section-specific random amplified polymorphic DNA (RAPD) markers in Lilium

Random amplified polymorphic DNA (RAPD) markers were utilized for the identification of Lilium species and inter-specific hybrids. The optimum annealing temperature of the polymerase chain reaction (PCR) for the RAPD assay in Lilium was 54 °C, which is relatively higher than the temperature used for other genera reported by previous researchers. Among 76 primers used to amplify genomic DNA by PCR, 18 primers (24%) generated polymorphic DNA fragments in Lilium species and hybrids. Cultivars were also identified by RAPD markers. Some amplified fragments were unique to species of each section and to hybrids derived from these species; that is, they were the section-specific DNA markers. Sections, Sinomartagon, Leucolirion b, Leucolirion a and Archelirion could be identified by 6 section-specific markers amplified with five primers. Seven inter-section hybrids showed the section-specific bands of both parental sections, indicating that these markers would be useful for identifying the parental sections of inter-section hybrids.     

PCR 扩增人型结核杆菌 DNA 及结核杆菌属 DNA 探针的制备

用多聚酶链反应(PCR)方法扩增人型、牛型结核杆菌基因组 DNA,获得特异的158 bpDNA 片段,而从另外十三种分枝杆菌未见到特异的扩增产物.回收158 bpDNA 片段作探针,它除与人型、牛型结核杆菌有特异的杂交信号外,与金黄色葡萄球菌、绿脓杆菌及一些分枝杆菌皆没有杂交反应.结果表明,PCR 可用于检测结核杆菌基因组 DNA,扩增产物158 bp DNA 片段可作为探针用于检测人型、牛型结核杆菌并鉴别结核杆菌与其它分枝杆菌.     

Construction of poplar (Populus tremula) chromosome 1-specific DNA library by using a microdissection technique

A method for single-chromosome microdissection and microcloning was established in forest plants using poplar (Populus tremula) as a model. By use of meristematic cell division in root tip and the wall degradation hypotonic method, well-spread poplar metaphase chromosome spreads showing low contamination were quickly prepared and fitted for chromosome microdissection. An individual chromosome 1 was microdissected from the metaphase spreads of poplar root-tip cells with a fine glass needle controlled by a micromanipulator. The dissected chromosome was amplified in vitro by theSau3A linker adaptor-mediated PCR technique, by which 200- to 3000-bp smear DNA fragments were obtained. Southern hybridization results showed that the PCR products from the single poplar chromosome were homogeneous with poplar genomic DNA, indicating that DNA from the single chromosome has been successfully amplified. Next, the second-round PCR products from the single chromosome 1 were cloned into T-easy vectors to generate a DNA library of the chromosome 1. About 3×10⁵ recombinant clones were obtained. Evaluation based on 160 randomly selected clones showed that the sizes of the cloned inserts varied from 230–2200 bp, with an average of 800 bp. Therefore, this research suggests that microdissection and microcloning of single small chromosomes in forest plants is feasible.     

Molecular identification of isolates of <Emphasis Type="Italic">Peronosclerospora sorghi</Emphasis> from maize using PCR-based SCAR marker

The fungus Peronosclerospora sorghi [Weston and Uppal (Shaw)] infects both sorghum and maize and incites downy mildew disease. Pathogenic and molecular variability among isolates of P. sorghi from sorghum and maize has been reported. In the present study we developed a DNA sequence characterized amplified region (SCAR) marker for identification of isolates of P. sorghi from maize by using polymerase chain reaction (PCR). The random amplified polymorphic DNA (RAPD) primer OPB15 consistently amplified a 1,000 base pairs (bp) product in PCR only from DNA of P. sorghi isolates from maize and not from isolates of sorghum. The PCR-amplified 1,000-bp product was cloned and sequenced. The sequence of the SCAR marker was used for designing specific primers for identification of maize isolates of P. sorghi. The SCAR primers amplified a 800 bp fragment only from genomic DNA of maize isolates of P. sorghi. The SCAR primers developed in this study are highly specific and reproducible, and proved to be powerful tool for identification of P. sorghi isolates from maize.     

Sequence subfamilies of satellite repeats related to rDNA intergenic spacer are differentially amplified on<Emphasis Type="Italic"> Vicia sativa</Emphasis> chromosomes

We cloned and sequenced the Vicia sativa 25S-18S rDNA intergenic spacer (IGS) and the satellite repeat S12, thought to be related to the spacer sequence. The spacer was shown to contain three types of subrepeats (A, B, and C) with monomers of 173 bp (A), 10 bp (B), and 66 bp (C), separated by unique or partially duplicated sequences. Two spacer variants were detected in V. sativa that differed in length (2990 and 3168 bp) owing to an extra copy of the subrepeat A. The A subrepeats were also shown to be highly homologous to the satellite repeat S12, which is located in large clusters on chromosomes 4, 5, and 6, and is not associated with the rDNA loci. Sequencing of additional S12 clones retrieved from a shotgun genomic library allowed definition of three subfamilies of this repeat based on minor differences in their nucleotide sequences. Two of these subfamilies could be discriminated from the rest of the S12 sequences as well as from the IGS A subrepeats using specific oligonucleotide primers that labeled only a subset of the S12 loci when used in the primed in situ DNA labeling (PRINS) reaction on mitotic chromosomes. These experiments showed that, in spite of the high overall similarity of the IGS A subrepeats and the S12 satellite repeats, there are S12 subfamilies that are divergent from the common consensus and are present at only some of the chromosomes containing the S12 loci. Thus, the subfamilies may have evolved at these loci following the spreading of the A subrepeats from the IGS to genomic regions outside the rDNA clusters.Electronic Supplementary Material Supplementary material is available in the online version of this article at Accession numbers: GenBank AY234364–AY234374. The monomer sequences and additional information about the family of IGS-like repeat S12 will also appear in the PlantSat database (Macas et al. 2002, ) under Accession name Vicia_sativa_IGS-like     

Noninvasive Method of DNA Isolation From Fecal Epithelial Tissue of Dairy Animals

A novel noninvasive genomic DNA isolation protocol from fecal tissue, by the proteinase K digestion and guanidine hydrochloride extraction method, was assessed for the genotyping of cattle and buffalo. The epithelial tissues present on the surface of the feces were used as source for isolation of genomic DNA. The DNA isolated from fecal tissue was found to be similar as those obtained from other body tissues such as skin, brain, liver, kidney, and muscle. The quality of DNA was checked by agarose gel electrophoresis and polymerase chain reaction (PCR). We successfully amplified a 320 bp MHC class II DRB gene and a 125 bp mt-DNA D-loop region from isolated genomic DNA of cattle. Thus, the DNA isolated using this method was suitable for common molecular biology methods, such as restriction enzyme digestion and genotyping of dairy animals through PCR.     

The Behaviour of 5-Hydroxymethylcytosine in Bisulfite Sequencing

Background

We recently showed that enzymes of the TET family convert 5-mC to 5-hydroxymethylcytosine (5-hmC) in DNA. 5-hmC is present at high levels in embryonic stem cells and Purkinje neurons. The methylation status of cytosines is typically assessed by reaction with sodium bisulfite followed by PCR amplification. Reaction with sodium bisulfite promotes cytosine deamination, whereas 5-methylcytosine (5-mC) reacts poorly with bisulfite and is resistant to deamination. Since 5-hmC reacts with bisulfite to yield cytosine 5-methylenesulfonate (CMS), we asked how DNA containing 5-hmC behaves in bisulfite sequencing.

Methodology/Principal Findings

We used synthetic oligonucleotides with different distributions of cytosine as templates for generation of DNAs containing C, 5-mC and 5-hmC. The resulting DNAs were subjected in parallel to bisulfite treatment, followed by exposure to conditions promoting cytosine deamination. The extent of conversion of 5-hmC to CMS was estimated to be 99.7%. Sequencing of PCR products showed that neither 5-mC nor 5-hmC undergo C-to-T transitions after bisulfite treatment, confirming that these two modified cytosine species are indistinguishable by the bisulfite technique. DNA in which CMS constituted a large fraction of all bases (28/201) was much less efficiently amplified than DNA in which those bases were 5-mC or uracil (the latter produced by cytosine deamination). Using a series of primer extension experiments, we traced the inefficient amplification of CMS-containing DNA to stalling of Taq polymerase at sites of CMS modification, especially when two CMS bases were either adjacent to one another or separated by 1–2 nucleotides.

Conclusions

We have confirmed that the widely used bisulfite sequencing technique does not distinguish between 5-mC and 5-hmC. Moreover, we show that CMS, the product of bisulfite conversion of 5-hmC, tends to stall DNA polymerases during PCR, suggesting that densely hydroxymethylated regions of DNA may be underrepresented in quantitative methylation analyses.     

利用甲基化特异性引物高通量检测DNA甲基化

建立一种基于甲基化特异性引物和SAGE技术的高通量DNA甲基化定量检测新方法（MSP-SAGE）,首先利用亚硫酸氢钠对基因组DNA进行处理,使未甲基化的C转变为U,而甲基化的CpG不变.将处理和未处理的DNA双链变性后用随机引物PNNNNCG对存在含有CG的单链进行延伸,而无甲基化CG的单链处则不能延伸;将差异延伸的单链序列和频次信息经过系列分子操作后,引入PCR扩增模板;对中间带有未知序列的PCR扩增产物进行串连克隆测序.将来自于未处理组和处理组的某一CpG位点的序列出现的次数定义为［Tags］A和［Tags］B,将标准系列的实际甲基化水平和［Tags］B/［Tags］A之间建立线性回归方程.根据每一CpG位点的［Tags］B/［Tags］A比值可反推该位点的甲基化水平.MSP-SAGE具有良好的线性,基于标准系列的［Tags］B/［Tags］A与其实际甲基化水平的标准曲线方程为y＝1.455x（R²＝0.984,P<0.01）.MSP-SAGE的回收率在95%到110%之间,精确度位于4.2％和10.5%,检测限在3％左右,单次检测通量可达24个CpG位点.MSP-SAGE是一种很有应用前途的高通量DNA甲基化定量检测方法.     

Evidence that two genomic species of Rhizobium are associated with Medicago truncatula

Seventy-three isolates of rhizobia sampled from root nodules of Medicago truncatula were analyzed by restriction fragment length polymorphism (RFLP) of DNA regions amplified by the polymerase chain reaction (PCR) targeting the symbiotic plasmid (nifD-K, nodD1, and nodD2 genes) and the chromosome (16S rDNA plus intergenic spacer). Two genotypic groups were found, regardless of the DNA region targeted. These two groups were given the status of genomic species based on results of DNA/DNA hybridization. Received: 1 August 1995 / Accepted: 13 October 1995     

Molecular detection of Colletotrichum falcatum causing red rot disease of sugarcane (Saccharum officinarum) using a SCAR marker

Red rot disease of sugarcane caused by Colletotrichum falcatum is one of the most destructive diseases of sugarcane (Saccharum officinarum) worldwide. The pathogen spreads primarily through infected sugarcane setts and hence the use of disease‐free planting materials is essential for preventing disease development in the field. In the present study a polymerase chain reaction (PCR) assay was developed for accurate and sensitive detection of C. falcatum in planting materials. Randomly amplified polymorphic DNA (RAPD) analysis identified a 566 bp PCR fragment that was specific to C. falcatum. The DNA sequence of this fragment was determined and used to design oligonucleotides amplifying a 442 bp sequence characterised amplified region (SCAR). The specificity of the SCAR primers was evaluated using purified DNA from C. falcatum and other Colletotrichum spp. as templates in PCR. The results indicated that the SCAR primers were highly specific to C. falcatum since the 442 bp fragment was amplified only from DNA of isolates and races of C. falcatum but not from any other Colletotrichum spp. tested. The detection sensitivity of C. falcatum was 0.1 ng for genomic DNA of C. falcatum and 5 ng for DNA extracted from infected sugarcane tissue. This new PCR‐based assay is a convenient tool for detection of this important pathogen in seed canes to ensure production of sugarcane.     

A polymerase chain reaction (PCR) application for free‐living nematodes (Rhabditida) discrimination

Synthetic primers were designed and constructed based on Caenorhabditis elegans genomic sequences, targeting HSP70 related sequences. Genomic DNA derived from several free‐living bacterial feeding nematode species was subjected to polymerase chain reaction (PCR), attempting discrimination among them. The amplified DNA fragments exhibited a distinct and reproducible pattern that characterizes different nematode species and populations.     

Analysis of the methylation status of imprinted genes based on methylation-specific polymerase chain reaction combined with denaturing high-performance liquid chromatography

A procedure for the analysis of the methylation status of imprinted genes is described. The method offers a rapid and reliable alternative to conventional methods such as Southern blots and methylation-specific polymerase chain reaction (PCR) (i.e., allele-specific methylation-specific PCR). The efficient resolution of the differentially methylated alleles is demonstrated for three human imprinted genes: SNRPN, LIT1 (alias KCNQ1OT1), and H19. Abnormal imprinting of SNRPN is associated with the Angelman/Prader-Willi syndromes, and that of LIT1 and H19 with the Beckwith-Wiedemann syndrome. The method is based on methylation-specific PCR followed by denaturing high-performance liquid chromatography (MSP/DHPLC). Briefly, genomic DNA is initially subjected to an in vitro bisulfite treatment, whereby unmethylated cytosines are deaminated. Subsequent PCR amplifications, using primers specific for modified DNA, are aimed at DNA segments that show parent-of-origin-specific methylation. PCR conditions are chosen that allow an efficient amplification of both alleles. The PCR products representing the two alleles are identical in size; they differ, however, at a number of positions within the amplified DNA segment. The DHPLC analysis allows very efficient resolution of the two populations of PCR products, providing qualitative and quantitative results.