Monitoring of representational difference analysis subtraction procedures by global microarrays

Various approaches to the study of differential gene expression are applied to compare cell lines and tissue samples in a wide range of biological contexts. The compromise between focusing on only the important genes in certain cellular processes and achieving a complete picture is critical for the selection of strategy. We demonstrate how global microarray technology can be used for the exploration of the differentially expressed genes extracted through representational difference analysis (RDA). The subtraction of ubiquitous gene fragments from the two samples was demonstrated using cDNA microarrays including more than 32 000 spotted, PCR-amplified human clones. Hybridizations indicated the expression of 9100 of the microarray elements in a macrophage/foam cell atherosclerosis model system, of which many were removed during the RDA process. The stepwise subtraction procedure was demonstrated to yield an efficient enrichment of gene fragments overrepresented in either sample (18% in the representations, 86% after the first subtraction, and 88% after the second subtraction), many of which were impossible to detect in the starting material. Interestingly, the method allowed for the observation of the differential expression of several members of the low-abundant nuclear receptor gene family. We also observed a certain background level in the difference products of nondifferentially expressed gene fragments, warranting a verification strategy for selected candidate genes. The differential expression of several genes was verified by real-time PCR.     

Detection of differentially expressed genes in primary tumor tissues using representational differences analysis coupled to microarray hybridization.

The identification of differential gene expressionbetween cells is a frequent goal in modern biological research. Here we demonstrate the coupling of representational difference analysis (RDA) of cDNA with microarray analysis of the output for high throughput screening. Two primary Ewing's sarcoma tissue samples with different biological behavior in vivo were compared by RDA: one which was metastatic and progressed rapidly; the other localized and successfully treated. A modified RDA protocol that minimizes the necessary starting material was employed. After a reduced number of subtractive rounds, the output of RDA was shotgun cloned into a plasmid vector. Inserts from individual colonies from the subtracted library were amplified with vector-specific primers and arrayed at high density on glass slides. The arrays were then hybridized with differentially fluorescently labeled starting amplicons from the two tissues and fluorescent signals were measured at each DNA spot. We show that the relative amounts of fluorescent signal correlate well with the abundance of fragments in the RDA amplicon and in the starting mRNA. In our system, we analyzed 192 products and 173 (90%) were appropriately detected as being >2-fold differentially expressed. Fifty unique, differentially expressed clones were identified. Therefore, the use of RDA essentially provides an enriched library of differentially expressed genes, while analysis of this library with microarrays allows rapid and reproducible screening of thousands of DNA molecules simultaneously. The coupling of these two techniques in this system resulted in a large pool of differentially expressed genes.     

Genomic subtraction to identify and characterize sequences of Shiga toxin-producing Escherichia coli O91:H21

To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.     

Direct cloning of cell differential expression genes with full-length by a new strategy based on the multiple rounds of 'long distance' polymerase chain reaction and magnetic beads mediated subtraction.

The paper described a new cDNA subtractive cloning strategy. This strategy was based on the 'cap-finder' method, 'long distance' polymerase chain reaction (PCR), streptavidin magnetic beads mediated subtraction, and spin column chromatography. When PCR products were resolved on agarose gel after three rounds of subtraction, the 'single gene difference' group displayed a predominantly enriched band, but the 'multiple gene difference' group did not display any apparent difference. Of 200 clones inserted with 0.7-2 kb fragments from the 'multiple gene difference' group, 50% were identified to be new sequences and 35% were known. Of 100 new sequences, 35% contained coding regions and 75% were confirmed by dot-blotting to be differentially expressed by genes of the target cell. The results suggested that this strategy might be very efficient for full-length cloning of differentially expressed genes of the cell.     

Genomic Subtraction To Identify and Characterize Sequences of Shiga Toxin-Producing Escherichia coli O91:H21

To identify Shiga toxin-producing Escherichia coli genes associated with severe human disease, a genomic subtraction technique was used with hemolytic-uremic syndrome-associated O91:H21 strain CH014 and O6:H10 bovine strains. The method was adapted to the Shiga toxin-producing E. coli genome: three rounds of subtraction were used to isolate DNA fragments specific to strain CH014. The fragments were characterized by genetic support analysis, sequencing, and hybridization to the genome of a collection of Shiga toxin-producing E. coli strains. A total of 42 fragments were found, 19 of which correspond to previously identified unique DNA sequences in the enterohemorrhagic E. coli EDL933 reference strain, including 7 fragments corresponding to prophage sequences and others encoding candidate virulence factors, such a SepA homolog protein and a fimbrial usher protein. In addition, the subtraction procedure yielded plasmid-related sequences from Shigella flexneri and enteropathogenic and Shiga toxin-producing E. coli virulence plasmids. We found that lateral gene transfer is extensive in strain CH014, and we discuss the role of genomic mobile elements, especially bacteriophages, in the evolution and possible transfer of virulence determinants.     

小菜蛾对溴氰菊酯抗性相关的DNA片段的克隆与分析

采用代表性差异分析（representational difference analysis, RDA ）方法,以小菜蛾Plutella xylostella (L.)的敏感品系作为对照组,溴氰菊酯抗性近等基因系（near isogenic line）作为测试组,通过三轮消减杂交得到大小为 150～300 bp的差异扩增带,经亚克隆测序并与GenBank等数据库同源比较,发现差异片段与已知抗性相关基因无同源性;以随机选取的差异片段作探针进行Southern blot分析,在测试组扩增子和三轮差异产物中都获得了阳性结果,而在对照组扩增子中的结果是阴性的。这些结果表明测序片段可能是与抗性相关的新基因序列或调控序列。     

A cDNA RDA protocol using solid-phase technology suited for analysis in small tissue samples

cDNA representational difference analysis (cDNA RDA) is a PCR-based subtractive enrichment procedure for the cloning of differentially expressed genes. In this study, we have further developed the procedure to take advantage of solid-phase technology, and to facilitate the use of RDA when starting material is limited. Several parameters of the PCR-based generation of cDNA representations were investigated, and a solid-phase based purification step was introduced to simplify removal of digested adapter-ends and uncleaved fragments. The use of magnetic particles increased the speed of the method, and also eliminated the risk of carry-over contamination between iterative steps of subtraction and PCR amplification. The modified protocol was evaluated in monitoring differences in gene expression in (i) a rat system consisting of livers with and without growth hormone treatment, and in (ii) a human system consisting of normal colon and colon cancer.     

ILG1 : a new integrase-like gene that is a marker of bacterial contamination by the laboratory Escherichia coli strain TOP10F'

BACKGROUND: Identification of differentially expressed genes between normal and diseased states is an area of intense current medical research that can lead to the discovery of new therapeutic targets. However, isolation of differentially expressed genes by subtraction often suffers from unreported contamination of the resulting subtraction library with clones containing DNA sequences not from the original RNA samples. MATERIALS AND METHODS: Subtraction using cDNA representational difference analysis (RDA) was performed on human B cells from normal or common variable immunodeficiency patients. The material remaining after the subtraction was cloned and individual clones were sequenced. The sequence of one clone with similarity to integrases (ILG1, integrase-like gene-1) was used to obtain the full length cDNA sequence and as a probe for the presence of this sequence in RNA or genomic DNA samples. RESULTS: After five rounds of cDNA RDA, 23.3% of the clones from the resulting subtraction library contained Escherichia coli DNA. In addition, three clones contained the sequence of a new integrase, ILG1. The full length cDNA sequence of ILG1 exhibits prokaryotic, but not eukaryotic, features. At the DNA level, ILG1 is not similar to any known gene. At the protein level, ILG1 has 58% similarity to integrases from the cryptic P4 bacteriophage family (S clade). The catalytic domain of ILG1 contains the conserved features found in site-specific recombinases. The critical residues that form the catalytic active site pocket are conserved, including the highly conserved R-H-R-Y hallmark of these recombinases. Interestingly, ILG1 was not present in the original B cell populations. By probing genomic DNA, ILG1 could only be detected in the E. coli TOP10F' strain used in our laboratory for molecular cloning, but not in any of its precursor strains, including TOP10. Furthermore, bacteria cultured from the mouth of the laboratory worker who performed cDNA RDA were also positive for ILG1. CONCLUSIONS: In the course of our studies using cDNA RDA, we have isolated and identified ILG1, a likely active site-specific recombinase and new member of the bacteriophage P4 family of integrases. This family of integrases is implicated in the horizontal DNA transfer of pathogenic genes between bacterial species, such as those found in pathogenic strains of E. coli, Shigella, Yersinia, and Vibrio cholera. Using ILG1 as a marker of our laboratory E. coli strain TOP10F', our evidence suggests that contaminating bacterial DNA in our subtraction experiment is due to this laboratory bacterial strain, which colonized exposed surfaces of the laboratory worker. Thus, identification of differentially expressed genes between normal and diseased states could be dramatically improved by using extra precaution to prevent bacterial contamination of samples.     

银额果蝇B染色体特异性基因的初步研究

采用代表性差异分析法（RDA）研究了银额果蝇两个单雌系AKM46（含B染色体）和AGZ2（不含B染色体）两基因组间的差异。用AKM46作检测（tester）扩增子,AGZ2作驱赶（driver）扩增子,通过三轮消减杂交后,获得了6个差异片段（100bp～300bp）。亚克隆后,对11个片段测序并与GenBank数据库进行同源性比较分析,获得了9个新的序列。选择clone22及clone42进行Southern杂交分析,这两个片段仅在检测扩增子及第一、第二、第三轮差异片段中检测到杂交信号,而在驱赶扩增子检测不到杂交信号。证实了这两个片段来自含有B染色体的单雌系AKM46,而且可能是B染色体上的特异基因片段。     

cDNA representational difference analysis of ileal Peyer's patches in lambs after oral inoculation with scrapie

cDNA representational difference analysis (RDA) was used to study gene expression profiles in the ileal Peyer's patch of a lamb 1 week after oral inoculation with the scrapie agent. Twenty-five differentially expressed cDNA fragments were identified and cloned. Sequence analysis indicated seven novel gene sequences. Other clones shared sequence homology with genes encoding ribosomal and mitochondrial proteins, the translation initiation factor EIF4GII and the bovine pancreatic thread protein. Reverse Northern was used to confirm the differential expression in another four lambs inoculated with scrapie and the tissue distribution of the novel genes was examined using Northern blot analysis.     

Targeted mapping of rye chromatin in wheat by representational difference analysis.

A targeted mapping strategy using representational difference analysis (RDA) was employed to isolate new restriction fragment length polymorphism probes for the long arm of chromosome 6 in rye (6RL), which carries a gene for resistance to Hessian fly larvae. Fragments from the 6RL arm were specifically isolated using a 'Chinese Spring' (CS) wheat - rye ditelosomic addition line (CSDT6RL) as tester, and CS and (or) CS4R as the driver for the genomic subtraction. Three RDA experiments were performed using BamHI amplicons, two of which were successful in producing low-copy clones. All low-copy clones were confirmed to have originated from 6RL, indicating substantial enrichment for target sequences. Two mapping populations, both of which are derived from a cross between two similar wheat-rye translocation lines, were used to map five RDA probes as well as five wheat probes. One of the populations was prescreened for recombinants by C-banding analysis. Fifteen loci, including seven new RDA markers, were placed on a map of the distal half of 6RL. The Hessian fly resistance gene was localized by mapping and C-banding analysis to approximately the terminal 1% of the arm. The utility of RDA as a method of targeted mapping in cereals and prospects for map-based cloning of the resistance gene are discussed.     

Construction and use of signal sequence selection vectors in Escherichia coli and Bacillus subtilis.

To study the diversity and efficiency of signal peptides for secreted proteins in gram-positive bacteria, two plasmid vectors were constructed which were used to probe for export signal-coding regions in Bacillus subtilis. The vectors contained genes coding for extracellular proteins (the alpha-amylase gene from Bacillus licheniformis and the beta-lactamase gene from Escherichia coli) which lacked a functional signal sequence. By shotgun cloning of restriction fragments from B. subtilis chromosomal DNA, a great variety of different export-coding regions were selected. These regions were functional both in B. subtilis and in E. coli. In a number of cases where protein export had been restored, intracellular precursor proteins of increased size could be detected, which upon translocation across the cellular membrane were processed to mature products. The high frequency with which export signal-coding regions were obtained suggests that, in addition to natural signal sequences, many randomly cloned sequences can function as export signal.     

Trypanosoma cruzi: subtractive hybridization as a molecular strategy to generate new targets to distinguish groups and hybrids

RAPD analysis and sequences of the mini-exon and ribosomal genes show that Trypanosoma cruzi can be clustered into two phylogenetic groups-T. cruzi I and II. Herein, the Representational Difference Analysis (RDA) method was used, providing new targets specific for each group. After three rounds of RDA hybridizing F strain (tester) with Y strain (driver) and vice-versa, an F-specific (F#30) and Y-specific (Y#22) clone were obtained specifically recognizing isolates from Amazonas (T. cruzi I) and Piauí (T. cruzi II). These segments corresponded to an unspecified protein (F#30) and a trans-sialidase (Y#22). Analysis of the F#30 sequence in T. cruzi I, T. cruzi II and zymodeme 3 samples displayed negligible specific differences that distinguished each group. In addition this F#30 gene has great potential as a hybrid marker.     

A genetic analysis of neural progenitor differentiation

Genetic mechanisms regulating CNS progenitor function and differentiation are not well understood. We have used microarrays derived from a representational difference analysis (RDA) subtraction in a heterogeneous stem cell culture system to systematically study the gene expression patterns of CNS progenitors. This analysis identified both known and novel genes enriched in progenitor cultures. In situ hybridization in a subset of clones demonstrated that many of these genes were expressed preferentially in germinal zones, some showing distinct ventricular or subventricular zone labeling. Several genes were also enriched in hematopoietic stem cells, suggesting an overlap of gene expression in neural and hematopoietic progenitors. This combination of methods demonstrates the power of using custom microarrays derived from RDA-subtracted libraries for both gene discovery and gene expression analysis in the central nervous system.     

Identifying differences in mRNA expression by representational difference analysis of cDNA.

Detection of differentially regulated genes has been severely hampered by technical limitations. In an effort to overcome these problems, the PCR-coupled subtractive process of representational difference analysis (RDA) [Lisitsyn, N. et al. (1993) Science 259, 946-951] has been adapted for use with cDNA. In a model system, RAG-1 and RAG-2, the genes responsible for activating V(D)J recombination, were identified in a genomic transfectant by cDNA RDA in a small fraction of the time taken by conventional means. The system was also modified to eliminate expected difference products to facilitate the identification of novel genes. Additional alterations to the conditions allowed isolation of differentially expressed fragments. Several caffeine up-regulated clones were obtained from the pre-B cell line 1-8, including IGF-1B, and a predicted homologue of the natural killer cell antigen, NKR-P1. The approach was found to be fast, extremely sensitive, reproducible, and predominantly lacked false positives. cDNA RDA has the capacity and adaptability to be applied to a wide range of biological problems, including the study of single gene disorders, characterization of mutant and complemented cell types, developmental or post-event expression time courses, and examination of pathogen-host interactions.     

Y chromosome and other heterochromatic sequences of the Drosophila melanogaster genome: how far can we go?

Whole genome shotgun assemblies have proven remarkably successful in reconstructing the bulk of euchromatic genes, with the only limit appearing to be determined by the sequencing depth. For genes imbedded in heterochromatin, however, the low cloning efficiency of repetitive sequences, combined with the computational challenges, demand that additional clues be used to annotate the sequences. One approach that has proven very successful in identifying protein coding genes in Y-linked heterochromatin of Drosophila melanogaster has been to make a BLASTable database of the small, unmapped contigs and fragments leftover at the end of a shotgun assembly, and to attempt to capture these by blasting with an appropriate query sequence. This approach often yields a staggered alignment of contigs from the unmapped set to the query sequence, as though the disjoint contigs represent small portions of the gene. Further inspection frequently shows that the contigs are broken by very large, heterochromatic introns. Methods of this sort are being expanded to make best use of all available clues to determine which unmapped contigs are associated with genes. These include use of EST libraries, and, in the case of the Y chromosome, testing of male specific genes and reduced shotgun depth of relevant contigs. It appears much more hopeful than anyone would have imagined that whole genome shotgun assemblies can recover the great bulk of even heterochromatic genes.