An improved method for subtractive cloning of differentially expressed genes in higher plants by protective exonuclease digestion and discriminating PCR amplification

An improved method for subtractive cloning with enhanced efficiency was developed by modifying the enzymatic degrading subtraction. The thionucleotide-modified tester cDNA fragments under control of one linker-primer were hybridized with excess driver cDNA fragments flanked by the other distinct linker-primer. After selective digestion of incompletely protected tester/driver and of unprotected driver/driver molecules with exonuclease III and VII, the protected tester/tester reassociates due to thionucleotides were exclusively amplified by PCR with the tester-cDNA-specific primer. The subtractively enriched target cDNA fragments, showing distinct bands in an agarose gel, were inserted into pUC19, and random colonies with inserts were screened by Northern hybridization to tester and driver RNA. Four distinct clones were confirmed to be up-regulated by the withdrawal of potassium from the nutrient solution of seedling barley growing hydroponically. The original protocol generated only smeared amplicons due to non-selective PCR amplification of the hybridized cDNA mixture including remains of undigested driver cDNA.Abbreviations EDS Enzymatic degrading subtraction - SET Subtractively enriched target     

The Isolation of differentially expressed mRNA sequences by selective amplification via biotin and restriction-mediated enrichment

Molecular analysis of development frequently implies the isolation and characterization of genes with specific spatial and temporal expression patterns. Several methods have been developed to identify such DNA sequences. The most comprehensive technique involves the genomewide probing of DNA sequence microarrays with mRNA sequences. However, at present this technology is limited to the few organisms for which the entire genome has been sequenced. Here, we describe a subtractive hybridization technique, called selective amplification via biotin and restriction-mediated enrichment (SABRE), which allows the selective amplification of cDNA fragments representing differentially expressed mRNA species. The method involves the competitive hybridization of an excess of driver cDNA fragments (D) to a trace of tester cDNA fragments (T), and the subsequent purification of tester homohybrids (in which both strands are contributed by the tester cDNA). After competitive hybridization, cDNA fragments that are more abundant in the tester than in the driver are enriched in the tester homohybrids. However, as the fraction of tester homohybrids is very small [T(2)/(D + T)(2)], their purification requires highly efficient procedures. In SABRE, the isolation of tester homohybrids is afforded by a combination of three successive steps: removal of biotinylated terminal sequences from most of the heterohybrids by S1 nuclease digestion, capture of biotinylated hybrids with streptavidin-coated paramagnetic beads, and specific release of homohybrids from the beads by restriction nuclease digestion. If several rounds of SABRE selection are performed in series, even relatively rare differentially expressed mRNA sequences may result in the production of predominant cDNA fragments in the final tester homohybrid population.     

Suppressive subtractive hybridization as a tool for identifying genetic diversity in an environmental metagenome: the rumen as a model

Molecular techniques previously used for genome comparisons of closely related bacterial species could prove extremely valuable for comparisons of complex microbial communities, or metagenomes. Our study aimed to determine the breadth and value of suppressive subtractive hybridization (SSH) in a pilot-scale analysis of metagenomic DNA from communities of microorganisms in the rumen. Suppressive subtractive hybridization was performed using total genomic DNA isolated from rumen fluid samples of two hay-fed steers, arbitrarily designated as tester or driver. Ninety-six subtraction DNA fragments from the tester metagenome were amplified, cloned and the DNA sequences were determined. Verification of the isolation of DNA fragments unique to the tester metagenome was accomplished through dot blot and Southern blot hybridizations. Tester-specific SSH fragments were found in 95 of 96 randomly selected clones. DNA sequences of subtraction fragments were analysed by computer assisted DNA and amino acid comparisons. Putative translations of 26 (32.1%) subtractive hybridization fragments exhibited significant similarity to Bacterial proteins, whereas 15 (18.5%) distinctive subtracted fragments had significant similarity to proteins from Archaea. The remainder of the subtractive hybridization fragments displayed no similarity to GenBank sequences. This metagenomic approach has exposed an unexpectedly large difference in Archaeal community structure between the rumen microbial populations of two steers fed identical diets and housed together. 16S rRNA dot blot hybridizations revealed similar proportions of Bacteria and Archaea in both rumen samples and suggest that the differences uncovered by SSH are the result of varying community structural composition. Our study demonstrates a novel approach to comparative analyses of environmental microbial communities through the use of SSH.     

Identification of repetitive, genome-specific probes in crucifer oilseed species.

Direct amplification of minisatellite DNA by PCR (DAMD PCR) was used to amplify and subsequently clone several fragments of DNA from crucifer species. The PCR-derived fragments of DNA were generated using known minisatellite core sequences as PCR primers. Southern hybridization of these putative minisatellite DNA fragments revealed that many were genome-specific; they hybridized with high affinity only to the genomic DNA of the species from which they were cloned. The DNA fragments were believed to be dispersed in the genome, based on smear-like hybridization signals on EcoRI-, BamHI-, and HindIII-digested genomic DNA. Genome-specific probes were specifically isolated from Brassica rapa (A genome), Brassica nigra (B genome), and Sinapis alba in addition to several other crucifer species. The sequence of a B. rapa specific probe (pBr17.1.3A) contained a minisatellite region that could be divided into three tandem repeats; each repeat contained between two and five subrepeats and each subrepeat shared a highly conserved core region of 29 bp. This minisatellite sequence also hybridized with high affinity to the A genome species B. napus and B. juncea. This research showed that dispersed, genome-specific probes can be isolated using DAMD PCR and that these probes could be used to detect and quantify alien DNA present in progeny from intergeneric or interspecific crosses.     

Detection and mapping of homologous, repeated and amplified DNA sequences by DNA renaturation in agarose gels.

A new molecular hybridization approach to the analysis of complex genomes has been developed. Tracer and driver DNAs were digested with the same restriction enzyme(s), and tracer DNA was labeled with 32P using T4 DNA polymerase. Tracer DNA was mixed with an excess amount of driver, and the mixture was electrophoresed in an agarose gel. Following electrophoresis, DNA was alkali-denatured in situ and allowed to reanneal in the gel, so that tracer DNA fragments could hybridize to the driver only when homologous driver DNA sequences were present at the same place in the gel, i.e. within a restriction fragment of the same size. After reannealing, unhybridized single-stranded DNA was digested in situ with S1 nuclease. The hybridized tracer DNA was detected by autoradiography. The general applicability of this technique was demonstrated in the following experiments. The common EcoRI restriction fragments were identified in the genomes of E. coli and four other species of bacteria. Two of these fragments are conserved in all Enterobacteriaceae. In other experiments, repeated EcoRI fragments of eukaryotic DNA were visualized as bands of various intensity after reassociation of a total genomic restriction digest in the gel. The situation of gene amplification was modeled by the addition of varying amounts of lambda phage DNA to eukaryotic DNA prior to restriction enzyme digestion. Restriction fragments of lambda DNA were detectable at a ratio of 15 copies per chicken genome and 30 copies per human genome. This approach was used to detect amplified DNA fragments in methotrexate (MTX)-resistant mouse cells and to identify commonly amplified fragments in two independently derived MTX-resistant lines.     

A simulation of subtractive hybridization.

Various strategies employed in genomic DNA cloning by subtractive hybridization have been examined by computer simulations, with the comparison between the predictions and the published results. The result shows that the efficiency of target sequence enrichment and the sensitivity to experimental conditions depend strongly on the enrichment strategy employed. The strategy selecting only tester/tester after hybridization can be very efficient to enrich targets. For successful target enrichment, however, the strategy requires a highly efficient subtraction method and proper hybridization conditions. The strategy also requires that the selected DNA be amplified by polymerase chain reaction (PCR) after each or each alternate subtraction. By contrast, the strategy selecting tester/tester plus single-stranded tester is less sensitive to various experimental factors, compared with the strategy selecting only tester/tester. However, it is not as efficient. With this strategy, the tester DNA selected may or may not be amplified by PCR before the next round. In the case of the strategy selecting single-stranded tester, the target DNA can be successfully enriched only when the selected DNA is directly used without PCR amplification in the next round. The strong features of existing methods can be combined to develop a protocol that is more efficient and more reliable.     

Isolation of pathogen-induced Chinese cabbage genes by subtractive hybridization employing selective adaptor ligation

We have developed a subtractive cloning method in which target sequences are effectively enriched by selective adaptor ligation and PCR after hybridization. In this method both tester and driver DNAs are digested with RsaI, ligated with the linker DNA containing a KpnI recognition site, and amplified by PCR. The tester DNA samples are divided into two aliquots, each digested with either RsaI or KpnI. The two DNA samples are then combined and hybridized with an excess of the driver DNA retaining the linker. After hybridization, the DNA mixture is ligated to a new adaptor compatible only with double-stranded tester/tester DNAs. Therefore, only the tester/tester is selectively amplified in subsequent PCR. This also leads to complete elimination of the tester DNA hybridized with driver DNA from the tester DNA population. Although our protocol employs enzymatic treatments, the efficiency of the enzymatic treatments does not affect the subtraction efficiency. This new subtractive enrichment method was applied to isolate Chinese cabbage defense-related genes induced by Pseudomonas syringae pv. tomato (Pst), which elicits a hypersensitive response in Chinese cabbage. After two or three rounds of subtractive hybridization, the sequences of enriched DNAs were determined and examined by BLAST analysis. Northern blot hybridization showed that 12 of the 19 genes analyzed were strongly induced by Pst treatment. Among the 12 Pst-induced genes five represent pathogenesis-related genes encoding PR1a, two chitinases, a thaumatin-like protein, and a PR4 protein. Other Pst-induced genes include two cytochrome P450 genes responsible for glucosinolate biosynthesis, a disease resistance gene homolog, and several genes encoding proteins with unknown functions.     

Isolation of DNA markers informative in purebred dog families by genomic representational difference analysis (gRDA)

Genomic Representational Difference Analysis (gRDA) is a subtractive DNA method to clone the differences between two related genomes, called tester and driver. We have evaluated this method to obtain polymorphic DNA markers for pedigree dogs. Amplified size-selected genomic restriction fragments (amplicons) of two dog littermates were repeatedly hybridized to each other in order to remove (subtract) those restriction fragments common to both sibs. Already after two rounds of subtractive hybridization, a clear enrichment of presumably tester-specific restriction fragments was observed, which was even more pronounced after the third round of subtraction. A plasmid library of 3000 recombinant clones was constructed of the second round and of the third round difference product. DNA sequence determination of randomly chosen clones of each difference product showed that approximately 1000 unique clones were obtained in the second-round difference product and approximately 500 in the third-round difference product. About half of the clones identified in the second-round difference product were also present in the third-round difference product. Of the second-round difference product, 39 different gRDA fragments could be identified, of which 21 were tester specific. In the third-round difference product, 22 different gRDA fragments were identified, of which 18 were tester specific. There were 13 fragments in common, resulting in a total of 48 different fragments. In order to establish the localization of these markers, we performed mapping using the dog radiation hybrid panel RHDF5000. Of 39 mapped clones, 29 were mapped to 20 existing RH groups, and 10 remained unlinked. It is concluded that gRDA is suitable to generate DNA markers to track disease genes within lines of pedigree dogs. Received: 26 April 2000 / Accepted: 11 May 2000     

A PCR-based amplification method retaining the quantitative difference between two complex genomes

With the increasing emergence of genome-wide analysis technologies (including comparative genomic hybridization (CGH), expression profiling on microarrays, differential display (DD), subtractive hybridization, and representational difference analysis (RDA)), there is frequently a need to amplify entire genomes or cDNAs by PCR to obtain enough material for comparisons among target and control samples. A major problem with PCR is that amplification occurs in a nonlinear manner and reproducibility is influenced by stray impurities. As a result, when two complex DNA populations are amplified separately, the quantitative relationship between two genes after amplification is generally not the same as their relation before amplification. Here we describe balanced PCR, a procedure that faithfully retains the difference among corresponding amplified genes by using a simple principle. Two distinct genomic DNA samples are tagged with oligonucleotides containing both a common and a unique DNA sequence. The genomic DNA samples are pooled and amplified in a single PCR tube using the common DNA tag. By mixing the two genomes, PCR loses the ability to discriminate among the different alleles and the influence of impurities is eliminated. The PCR-amplified pooled samples can be separated using the DNA tag unique to each individual genomic DNA sample. The principle of this method has been validated with synthetic DNA, genomic DNA, and cDNA applied on microarrays. By removing the bias of PCR, this method allows a balanced amplification of allelic fragments from two complex DNAs even after three sequential rounds of PCR. This balanced PCR approach should allow genetic analysis in minute laser-microdissected tissues, paraffin-embedded archived material, or single cells.     

The use of double fluorescence in situ hybridization to physically map the positions of 5S rDNA genes in relation to the chromosomal location of 18S-5.8S-26S rDNA and a C genome specific DNA sequence in the genus Avena.

A physical map of the locations of the 5S rDNA genes and their relative positions with respect to 18S-5.8S-26S rDNA genes and a C genome specific repetitive DNA sequence was produced for the chromosomes of diploid, tetraploid, and hexaploid oat species using in situ hybridization. The A genome diploid species showed two pairs of rDNA loci and two pairs of 5S loci located on both arms of one pair of satellited chromosomes. The C genome diploid species showed two major pairs and one minor pair of rDNA loci. One pair of subtelocentric chromosomes carried rDNA and 5S loci physically separated on the long arm. The tetraploid species (AACC genomes) arising from these diploid ancestors showed two pairs of rDNA loci and three pairs of 5S loci. Two pairs of rDNA loci and 2 pairs of 5S loci were arranged as in the A genome diploid species. The third pair of 5S loci was located on one pair of A-C translocated chromosomes using simultaneous in situ hybridization with 5S rDNA genes and a C genome specific repetitive DNA sequence. The hexaploid species (AACCDD genomes) showed three pairs of rDNA loci and six pairs of 5S loci. One pair of 5S loci was located on each of two pairs of C-A/D translocated chromosomes. Comparative studies of the physical arrangement of rDNA and 5S loci in polyploid oats and the putative A and C genome progenitor species suggests that A genome diploid species could be the donor of both A and D genomes of polyploid oats. Key words : oats, 5S rDNA genes, 18S-5.8S-26S rDNA genes, C genome specific repetitive DNA sequence, in situ hybridization, genome evolution.     

One-step PCR amplification of complete arthropod mitochondrial genomes

A new PCR primer set which enables one-step amplification of complete arthropod mitochondrial genomes was designed from two conserved 16S rDNA regions for the long PCR technique. For this purpose, partial 16S rDNAs amplified with universal primers 16SA and 16SB were newly sequenced from six representative arthropods: Armadillidium vulgare and Macrobrachium nipponense (Crustacea), Anopheles sinensis (Insecta), Lithobius forficatus and Megaphyllum sp. (Myriapoda), and Limulus polyphemus (Chelicerata). The genomic locations of two new primers, HPK16Saa and HPK16Sbb, correspond to positions 13314-13345 and 12951-12984, respectively, in the Drosophila yakuba mitochondrial genome. The usefulness of the primer set was experimentally examined and confirmed with five of the representative arthropods, except for A. vulgare, which has a linearized mitochondrial genome. With this set, therefore, we could easily and rapidly amplify complete mitochondrial genomes with small amounts of arthropod DNA. Although the primers suggested here were examined only with arthropod groups, a possibility of successful application to other invertebrates is very high, since the high degree of sequence conservation is shown on the primer sites in other invertebrates. Thus, this primer set can serve various research fields, such as molecular evolution, population genetics, and molecular phylogenetics based on DNA sequences, RFLP, and gene rearrangement of mitochondrial genomes in arthropods and other invertebrates.     

鼻咽癌上皮细胞株HNE1差异表达基因的分离与鉴定

为了分离鼻咽癌差异表达基因 ,应用抑制性扣除杂交技术 ,在正向抑制性扣除杂交中 ,以鼻咽癌上皮细胞株HNE1cDNA作为检测子 ,以人胚鼻咽上皮细胞cDNA作为驱赶子 ;在反向抑制性扣除杂交中 ,以人胚鼻咽上皮细胞cDNA作为检测子 ,以鼻咽癌上皮细胞株HNE1cDNA作为驱赶子 ,分别通过抑制性扣除杂交 ,构建了鼻咽癌上皮细胞株HNE1表达下调和表达上调的两个扣除cDNA文库 .从鼻咽癌相关的扣除cDNA文库中随机挑取 1 2 0 0个克隆 ,采用菌落PCR扩增其插入cDNA片段 ,自动点膜制备成cDNA微阵列膜 ,分别用鼻咽癌上皮细胞株HNE1、人胚鼻咽上皮mRNA经逆转录标记cDNA探针 ,分别与cDNA微阵列膜杂交 ,通过杂交信号的自动扫描分析 ,对杂交信号存在 5倍差异的克隆进行测序 ,获得了 1 0个鼻咽癌差异表达基因的cDNA片段 ,其中 3个为新基因序列 ,其GenBank登录号为 :AF5 1 0 1 88、AF5 1 0 1 89和AF5 1 0 1 90 ,7个代表已知基因序列 .采用RT PCR证实S1 0 0A8,CK1 9和RBP1基因在人胚鼻咽上皮中高表达而在鼻咽癌细胞株HNE1中低表达 .这些结果显示上述基因可能是鼻咽癌发生的重要因素     

Tailed pooled suppression subtractive hybridization (PSSH) adaptors do not alter efficiency

Suppression Subtractive Hybridization (SSH) and its derivative, Pooled Suppression Subtractive hybridization (PSSH), are powerful tools used to study variances larger than ~100 bp in prokaryotic genome structure. The initial steps involve ligating an oligonucleotide of known sequence (the “adaptor”) to a fragmented genome to facilitate amplification, subtraction and downstream sequencing. SSH results in the creation of a library of unique DNA fragments which have been traditionally analyzed via Sanger sequencing. Numerous next generation sequencing technologies have entered the market yet SSH is incompatible with these platforms. This is due to the high level of sequence conservation of the oligonucleotide used for SSH. This rigid adherence is partly because it has yet to be determined if alteration of this oligonucleotide will have a deleterious impact on subtraction efficiency. The subtraction occurs when non-unique fragments are inhibited by a secondary self-pairing structure which requires exact nucleotide sequence. We determine if appending custom sequence to the 5′ terminal ends of these oligonucleotides during the nested PCR stages of PSSH will reduce subtraction efficiency. We compare a pool of ten S. aureus clinical isolates with a standard PSSH and custom tailed-PSSH. We detected no statistically significant difference between their subtraction efficiencies. Our observations suggest that the adaptor’s terminal ends may be labeled during the nested PCR step. This produces libraries labeled with custom sequence. This does not lead to loss of subtraction efficiency and would be invaluable for groups wishing to combine SSH or PSSH with their own downstream applications, such as a high throughput sequencing platform.     

Identification of Genes Preferentially Expressed in Goat Hair Follicle Anagen-Catagen Transition Using Suppression Subtractive Hybridization

Suppressive subtraction hybridization (SSH) was used to identify differentially expressed genes in goat (Capra hircus) hair follicle anagen-catagen transition. The cDNA fragments, derived from SSH positive subtractive library (tester: anagen-catagen transition, driver: later anagen), were cloned into pEGM-T vector. Two hundred cDNA fragments screened from this library were subjected to identify forty-five unregulated isolates. Sequence analysis revealed that these fragments represented twenty-three genes. Blasting analysis with database in GenBank showed that twenty genes were previously clearly annotated, two were homologous to un-annotated expressed sequence tag (ESTs), and one might be novel. To identify characters of gene expression, seven genes in later anagen and anagen-catagen transition skin tissues were chosen for quantitative real-time PCR. Results indicated that expression of these seven genes varied much, reaching threefold among them, furthering indicating that expression of those genes was up-regulation in the anagen-catagen transition. We characterized expression levels of this potential novel gene and the goat ectodysplasin A during differential stages of hair cycle. These profiles suggested that these two genes might play a role in the goat secondary hair follicle cycle.     

Identification of genes preferentially expressed in goat hair follicle anagen-catagen transition using suppression subtractive hybridization

Suppressive subtraction hybridization (SSH) was used to identify differentially expressed genes in goat (Capra hircus) hair follicle anagen-catagen transition. The cDNA fragments, derived from SSH positive subtractive library (tester: anagen-catagen transition, driver: later anagen), were cloned into pEGM-T vector. Two hundred cDNA fragments screened from this library were subjected to identify forty-five unregulated isolates. Sequence analysis revealed that these fragments represented twenty-three genes. Blasting analysis with database in GenBank showed that twenty genes were previously clearly annotated, two were homologous to un-annotated expressed sequence tag (ESTs), and one might be novel. To identify characters of gene expression, seven genes in later anagen and anagen-catagen transition skin tissues were chosen for quantitative real-time PCR. Results indicated that expression of these seven genes varied much, reaching threefold among them, furthering indicating that expression of those genes was up-regulation in the anagen-catagen transition. We characterized expression levels of this potential novel gene and the goat ectodysplasin A during differential stages of hair cycle. These profiles suggested that these two genes might play a role in the goat secondary hair follicle cycle.     

A 9.6 kb intervening sequence in D. virilis rDNA, and sequence homology in rDNA interruptions of diverse species of Drosophila and other diptera.

A large proportion of the 28S ribosomal RNA genes in Drosophila virilis are interrupted by a DNA sequence 9.6 kilobase pairs long. As regards both its presence and its position in the 28S gene (about two thirds of the way in), the D. virilis rDNA intervening sequence is similar to that found in D. melanogaster rDNA, but lengths differ markedly between the two species. Degrees of nucleotide sequence homology have been detected bewteen rDNA interruptions of the two species. This homology extends to putative rDNA intervening sequences in diverse higher diptera (other Drosophila species, the house fly and the flesh fly), but hybridization of cloned D. melanogaster and D. virilis rDNA interruption segments to DNA of several lower diptera has been negative. As is the case with melanogaster rDNA interruptions, segments of the virilis rDNA intervening sequence hybridize with non-rDNA components of the virilis genome, and interspecific homology may involve these non-rDNA sequences as well as rDNA interruptions. There is, however, evidence from buoyant density fractionation of DNA that the distributions of interruption-related sequences are distinct in D. melanogaster and D. virilis genomes. Moreover, thermal denaturation studies have indicated differing extents of homology between hybridizable sequences in D. virilis DNA and different segments of the D. melanogaster rDNA intervening sequence. We infer from our studies that rDNA intervening sequences are prevalent among higher diptera; that in the course of the evolution of these organisms, elements of the intervening sequences have been moderately to highly conserved; and that this conservation extends in at least two distantly related species of Drosophila to similar sequences found elsewhere in the genomes.     

Comparison of mealybug (Planococcus lilacinus) and fruit fly genomes: isolation and analysis of conserved sequences and their utility in studying synteny in the mealybug

By using ligation-mediated PCR products from mealybug DNA as tester and biotinylated fly DNA as driver, we recovered a fraction of the tester that remains hybridized to driver following high-stringency washing conditions. This fraction is expected to contain mealybug sequences conserved in the fly (MCF). Reciprocal experiments enabled the isolation of fly sequences conserved in the mealybug (FCM). Coding sequences among MCF show amino acid identities >40% with fly proteins, allowing a reliable identification of orthologs. Three sequences from the fly cytogenetic positions 98-99 were hybridized onto mealybug chromosomes and the results identified differences in synteny between the two species. Taken together, our results present a method for direct isolation of sequences conserved between an 'orphan' (mealybug) genome and a 'reference' (fly) genome and showed that these sequences can be used to study chromosome synteny in the mealybug.     

Enhanced Methylation Analysis by Recovery of Unsequenceable Fragments

Bisulfite sequencing is a valuable tool for mapping the position of 5-methylcytosine in the genome at single base resolution. However, the associated chemical treatment causes strand scission, which depletes the number of sequenceable DNA fragments in a library and thus necessitates PCR amplification. The AT-rich nature of the library generated from bisulfite treatment adversely affects this amplification, resulting in the introduction of major biases that can confound methylation analysis. Here, we report a method that enables more accurate methylation analysis, by rebuilding bisulfite-damaged components of a DNA library. This recovery after bisulfite treatment (ReBuilT) approach enables PCR-free bisulfite sequencing from low nanogram quantities of genomic DNA. We apply the ReBuilT method for the first whole methylome analysis of the highly AT-rich genome of Plasmodium berghei. Side-by-side comparison to a commercial protocol involving amplification demonstrates a substantial improvement in uniformity of coverage and reduction of sequence context bias. Our method will be widely applicable for quantitative methylation analysis, even for technically challenging genomes, and where limited sample DNA is available.