Differential cDNA cloning by enzymatic degrading subtraction (EDS)

Systematic study of single-cell isolation from musculoskeletal tissues for single-sell sequencing

Background

The single-cell platform provided revolutionary way to study cellular biology. Technologically, a sophistic protocol of isolating qualified single cells would be key to deliver to single-cell platform, which requires high cell viability, high cell yield and low content of cell aggregates or doublets. For musculoskeletal tissues, like bone, cartilage, nucleus pulposus and tendons, as well as their pathological state, which are tense and dense, it’s full of challenge to efficiently and rapidly prepare qualified single-cell suspension. Conventionally, enzymatic dissociation methods were wildly used but lack of quality control. In the present study, we designed the rapid cycling enzymatic processing method using tissue-specific enzyme cocktail to treat different human pathological musculoskeletal tissues, including degenerated nucleus pulposus (NP), ossifying posterior longitudinal ligament (OPLL) and knee articular cartilage (AC) with osteoarthritis aiming to rapidly and efficiently harvest qualified single-cell suspensions for single-cell RNA-sequencing (scRNA-seq).

Results

We harvested highly qualified single-cell suspensions from NP and OPLL with sufficient cell numbers and high cell viability using the rapid cycling enzymatic processing method, which significantly increased the cell viability compared with the conventional long-time continuous digestion group (P < 0.05). Bioanalyzer trace showed expected cDNA size distribution of the scRNA-seq library and a clear separation of cellular barcodes from background partitions were verified by the barcode-rank plot after sequencing. T-SNE visualization revealed highly heterogeneous cell subsets in NP and OPLL. Unfortunately, we failed to obtain eligible samples from articular cartilage due to low cell viability and excessive cell aggregates and doublets.

Conclusions

In conclusion, using the rapid cycling enzymatic processing method, we provided thorough protocols for preparing single-cell suspensions from human musculoskeletal tissues, which was timesaving, efficient and protective to cell viability. The strategy would greatly guarantee the cell heterogeneity, which is critical for scRNA-seq data analysis. The protocol to treat human OA articular cartilage should be further improved.

Isolated Bacillus subtilis strain 330-2 and its antagonistic genes identified by the removing PCR

Plant growth-promoting bacteria (PGPB) may trigger tolerance against biotic/abiotic stresses and growth enhancement in plants. In this study, an endophytic bacterial strain from rapeseed was isolated to assess its role in enhancing plant growth and tolerance to abiotic stresses, as well as banded leaf and sheath blight disease in maize. Based on 16S rDNA and BIOLOG test analysis, the 330-2 strain was identified as Bacillus subtilis. The strain produced indole-3-acetic acid, siderophores, lytic enzymes and solubilized different sources of organic/inorganic phosphates and zinc. Furthermore, the strain strongly suppressed the in vitro growth of Rhizoctonia solani AG1-IA, Botrytis cinerea, Fusarium oxysporum, Alternaria alternata, Cochliobolus heterostrophus, and Nigrospora oryzae. The strain also significantly increased the seedling growth (ranging 14–37%) of rice and maize. Removing PCR analysis indicated that 114 genes were differentially expressed, among which 10%, 32% and 10% were involved in antibiotic production (e.g., srfAA, bae, fen, mln, and dfnI), metabolism (e.g., gltA, pabA, and ggt) and transportation of nutrients (e.g., fhu, glpT, and gltT), respectively. In summary, these results clearly indicate the effectiveness and mechanisms of B. subtilis strain 330-2 in enhancing plant growth, as well as tolerance to biotic/abiotic stresses, which suggests that the strain has great potential for commercialization as a vital biological control agent.

Induction of wheat defense and stress-related genes in response toFusarium graminearum

Fusarium head blight (FHB), caused by species of the fungus Fusarium, is a worldwide disease of wheat (Triticum aestivum L.). The Chinese T. aestivum 'Ning7840' is one of few wheat cultivars with resistance to FHB. To identify differentially expressed genes corresponding to FHB resistance, a cDNA library was constructed using pooled mRNA isolated from glumes of 'Ning7840' harvested at 2, 6, 12, 24, 36, 72, and 96 h after inoculation (hai) with a conidia spore suspension of Fusarium graminearum. Suppressive subtractive hybridization (SSH) cDNA subtraction was carried out using pooled glume mRNAs from the tester and the control. The cDNA library was differentially screened using the forward subtracted cDNAs and the reverse subtracted cDNAs as probes. Twenty-four clones with significant matches to either plant (16 sequences) or fungal (8 sequences) genes were isolated based on their specific hybridization with forward subtracted cDNA and not reverse subtracted cDNA. Six putative defense-related genes were confirmed by real-time quantitative reverse-transcriptase PCR. Many-fold higher induction of three clones (A3F8, B10H1, and B11H3) in the resistant genotypes compared with susceptible genotypes indicates a putative role in the resistance response to Fusarium graminearum. Transcript accumulations of P450, chitinase (Chi1), and one unknown gene (clone B8Q9) in both resistant and susceptible genotypes suggest an involvement in a generalized resistance response to F. graminearum. Nucleotide sequence analysis showed that cDNA clone A4C6 encodes a cytochrome P450 gene (CYP709C3v2), including 14 N-terminal amino acids that have a membrane-associated helical motif. Other domains characteristic of eukaryotic P450 are also present in CYP709C3v2. The deduced polypeptide of cDNA clone B2H2 encodes an acidic isoform of class I chitinase containing a 960-bp coding region. Southern hybridization using aneuploid lines of T. aestivum 'Chinese Spring' indicated that CYP709C3v2 was located on the short arm of chromosomes 2B and 2D.Key words: Fusarium head blight (FHB), suppressive subtractive hybridization, defense response, real-time quantitative RT-PCR.

Gene expression in the brain and kidney of rainbow trout in response to handling stress

Background

Microarray technologies are rapidly becoming available for new species including teleost fishes. We constructed a rainbow trout cDNA microarray targeted at the identification of genes which are differentially expressed in response to environmental stressors. This platform included clones from normalized and subtracted libraries and genes selected through functional annotation. Present study focused on time-course comparisons of stress responses in the brain and kidney and the identification of a set of genes which are diagnostic for stress response.

Results

Fish were stressed with handling and samples were collected 1, 3 and 5 days after the first exposure. Gene expression profiles were analysed in terms of Gene Ontology categories. Stress affected different functional groups of genes in the tissues studied. Mitochondria, extracellular matrix and endopeptidases (especially collagenases) were the major targets in kidney. Stress response in brain was characterized with dramatic temporal alterations. Metal ion binding proteins, glycolytic enzymes and motor proteins were induced transiently, whereas expression of genes involved in stress and immune response, cell proliferation and growth, signal transduction and apoptosis, protein biosynthesis and folding changed in a reciprocal fashion. Despite dramatic difference between tissues and time-points, we were able to identify a group of 48 genes that showed strong correlation of expression profiles (Pearson r > |0.65|) in 35 microarray experiments being regulated by stress. We evaluated performance of the clone sets used for preparation of microarray. Overall, the number of differentially expressed genes was markedly higher in EST than in genes selected through Gene Ontology annotations, however 63% of stress-responsive genes were from this group.

Conclusions

1. Stress responses in fish brain and kidney are different in function and time-course. 2. Identification of stress-regulated genes provides the possibility for measuring stress responses in various conditions and further search for the functionally related genes.

Suppression subtractive hybridization coupled with microarray analysis to examine differential expression of genes in virus infected cells

High throughput detection of differential expression of genes is an efficient means of identifying genes and pathways that may play a role in biological systems under certain experimental conditions. There exist a variety of approaches that could be used to identify groups of genes that change in expression in response to a particular stimulus or environment. We here describe the application of suppression subtractive hybridization (SSH) coupled with cDNA microarray analysis for isolation and identification of chicken transcripts that change in expression on infection of host cells with a paramyxovirus. SSH was used for initial isolation of differentially expressed transcripts, a large-scale validation of which was accomplished by microarray analysis. The data reveals a large group of regulated genes constituting many biochemical pathways that could serve as targets for future investigations to explore their role in paramyxovirus pathogenesis. The detailed methods described herein could be useful and adaptable to any biological system for studying changes in gene expression.

Transcriptional analysis of the response of poultry species to respiratory pathogens

Respiratory tract diseases are the single most important cause of economic loss due to infections among poultry populations worldwide. However, the molecular mechanisms of the host response to infections remain unknown. Here, we review the literature and describe the adoption of a conceptually simple approach to understand the genetic and biochemical responses of host cells during infection with respiratory pathogens, such as avian pneumovirus (APV). The strategy that we have adopted integrates the powerful techniques of cDNA subtraction hybridization and microarray analysis for global transcriptional profiling. The results of our investigations identify the specific transcriptional alterations in host-cell gene expression that result from an attempt by the host to combat and limit the spread of the pathogen or by the pathogen to enhance its own survival and ability to reproduce. Our studies suggest that a molecular description of host-pathogen interactions in terms of differential gene expression will provide key insights on the molecular basis of disease pathogenesis, pathogen virulence, and host immunity. In addition, the results suggest that the identification of genes and pathways with a role in host response to infection has considerable practical implications for the future design and development of effective immunomodulators and vaccines.

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.

A modification of Representational Difference Analysis, with application to the cloning of a candidate in the Reelin signalling pathway

BACKGROUND

cDNA-RDA is one of the subtractive cloning techniques used to isolate differentially expressed genes between two complex cDNA populations. In the present study we present a modification of the protocol described by Hubank and Schatz.

RESULTS

In the post-hybridization mix, the 5'-ends of homoduplexes of interest (tester-tester) are filled-in with alpha-thio-deoxynucleotides. Unprotected duplexes, as well as the single-stranded DNA fragments, are degraded using ExoIII and Mung Bean Nuclease, prior to PCR subtraction, resulting in less complex difference products. We illustrate this modification by the cloning of a new gene which is differentially expressed in normal, reelin and Dab1 mutant mice and is a candidate member of the Reelin signalling pathway involved in brain development.

CONCLUSION

We propose a modification of cDNA-RDA that may reduce the complexity of the post-hybridization mix and thus facilitate the amplification of differentially expressed products.

Digital analysis of cDNA abundance; expression profiling by means of restriction fragment fingerprinting

BACKGROUND

Gene expression profiling among different tissues is of paramount interest in various areas of biomedical research. We have developed a novel method (DADA, Digital Analysis of cDNA Abundance), that calculates the relative abundance of genes in cDNA libraries.

RESULTS

DADA is based upon multiple restriction fragment length analysis of pools of clones from cDNA libraries and the identification of gene-specific restriction fingerprints in the resulting complex fragment mixtures. A specific cDNA cloning vector had to be constructed that governed missing or incomplete cDNA inserts which would generate misleading fingerprints in standard cloning vectors. Double stranded cDNA was synthesized using an anchored oligo dT primer, uni-directionally inserted into the DADA vector and cDNA libraries were constructed in E. coli. The cDNA fingerprints were generated in a PCR-free procedure that allows for parallel plasmid preparation, labeling, restriction digest and fragment separation of pools of 96 colonies each. This multiplexing significantly enhanced the throughput in comparison to sequence-based methods (e.g. EST approach). The data of the fragment mixtures were integrated into a relational database system and queried with fingerprints experimentally produced by analyzing single colonies. Due to limited predictability of the position of DNA fragments on the polyacrylamid gels of a given size, fingerprints derived solely from cDNA sequences were not accurate enough to be used for the analysis. We applied DADA to the analysis of gene expression profiles in a model for impaired wound healing (treatment of mice with dexamethasone).

CONCLUSIONS

The method proved to be capable of identifying pharmacologically relevant target genes that had not been identified by other standard methods routinely used to find differentially expressed genes. Due to the above mentioned limited predictability of the fingerprints, the method was yet tested only with a limited number of experimentally determined fingerprints and was able to detect differences in gene expression of transcripts representing 0.05% of the total mRNA population (e.g. medium abundant gene transcripts).

Subtractive hybridization of biotinylated DNA in phenol emulsion

A method of subtractive hybridization using biotinylated DNA and phenol emulsion reassociation technique (PERT) has been proposed. A possibility of combining these techniques has been shown for the first time. The effect of biotinylation degree on the formation of water-insoluble Bio-DNA aggregates was studied. The conditions when Bio-DNA aggregation is actually absent were revealed. A possible use of the above method in hybridization experiments in a wide range of DNA concentrations has been shown. The time of hybridization was 0.5-1 h. The method was approved on a model system, and its possible application for the enrichment of rare mRNA was shown. No less than 300-fold enrichment is achieved for a rare transcript (IL-2) in three cycles of subtractive hybridization.

Cloning and analysis of mold-specific genes in the dimorphic fungus Histoplasma capsulatum

A critical feature in the pathogenesis of the respiratory pathogen Histoplasma capsulatum is the conversion from the mold form (found in soil) to the yeast form in the lungs of the host. Little is known about the molecular biology of Histoplasma dimorphism. In particular, the possible roles of genes which are transcriptionally silent in yeast (i.e. mold-specific) have not been studied. We have produced a cDNA library highly enriched for mold-upregulated clones by fragmenting cDNA and removing yeast-specific and common sequences with a highly efficient enzyme degrading subtraction method. Screening of randomly selected clones identified cDNA fragments representing 16 different mold-upregulated genes. Because multiple cDNA fragments can be treated as alleles in a genetic screen, we were able to apply probability analysis to estimate the total number of mold-upregulated genes. We estimate that there are 27 upregulated genes; cDNA fragments of 16 have been isolated. Here we report the first isolation and analysis of cDNA from two mold-specific genes, MS8 (GenBank AF292398) and MS88 (GenBank AF357882). The MS8 transcript was very strongly expressed in mold but not detected on Northern blots with yeast RNA. The putative MS8 protein was predicted to be 21.3 kDa (203 aa), very rich in glutamine and glycine and had a calculated pI of 6.76. The MS88 transcript was weakly expressed in mold and not detected in yeast. The putative MS88 protein was predicted to be 22.5 kDa (219 aa) with a pI of 4.46. GenBank similarity searches revealed that the putative MS8 protein was similar to a glutamine-rich protein, of unknown function, from the fungus Colletotrichum gloeosporioides (GenBank U94186). No significant matches were found for the putative MS88 protein.

A critical evaluation of differential display as a tool to identify genes involved in legume nodulation: looking back and looking forward

Screening for differentially expressed genes is a straightforward approach to study the molecular basis of a biological system. In the last 10 years, differential screening technology has evolved rapidly and currently high-throughput tools for genome-wide transcript profiling, such as expressed sequence tags and microarray analysis, are becoming widely available. Here, an overview of this (r)evolution is given with emphasis on the differential display method, which for many years has been the preferred technique of scientists in diverse fields of research. Differential display has also been the method of choice for the identification of genes involved in the symbiotic interaction between Azorhizobium caulinodans and Sesbania rostrata. The advantages with respect to tissue specificity of this particular model system for legume nodulation and the results of a screening for early nodulation-related genes have been considered in the context of transcriptome analyses in other rhizobium-legume interactions.

Lambda exonuclease-based subtractive hybridization approach to isolate differentially expressed genes from leaf cultures of Paulownia kawakamii

Genes that are preferentially expressed in a particular developmental pathway can be isolated by subtractive hybridization (SH). We developed a PCR-based approach coupled with lambda exonuclease digestion that allows for generating single-stranded tester and driver nucleic acids suitable for SH starting from cDNA libraries. An efficient subtraction strategy was developed to overcome some of the problems in the previously described SH protocols, such as the need for large amounts of experimental tissue, RNase contamination during solution hybridization, and postsubtraction recovery of nucleic acids. We used this method to obtain cDNA corresponding to genes expressed during adventitious shoot regeneration from excised leaf cultures of the fast-growing tree Paulownia kawakamii. Over 36 cDNA clones were isolated and 1 of the differentially expressed clones codes for a leucine zipper transcription factor. This clone showed about sixfold higher level of expression in the shoot-forming tissues (tester) compared to that in the callus-forming tissues (driver) of Paulownia, suggesting that differentially expressed genes can be efficiently isolated using this simple lambda exonuclease-based subtractive hybridization method.

Detection of differentially expressed genes in human colon carcinoma cells treated with a selective COX-2 inhibitor

Numerous reports suggest that use of nonsteroidal anti-inflammatory drugs (NSAIDs) decrease mortality from colorectal cancer. To better understand all of the mechanisms underlying this effect, the global pattern of gene expression in colon carcinoma cells following treatment with NS-398, a selective cyclo-oxygenase-2 inhibitor was evaluated. We utilized suppression subtractive hybridization combined with differential screening to identify genes whose expression was affected following treatment. Among the subtracted cDNA fragments confirmed as differentially expressed, there were two which are known to be involved in the regulation of cell adhesion (human FAT and proto-cadherin-7). We identified two other genes whose levels were decreased and these are known to be involved in the regulation of cell proliferation (cyclin K and p-100). We identified additional genes which are involved in different signaling pathways which regulate programmed cell death (Dynamin 2, Pdcd4 and LIP.1). These results provide evidence that some of the effects of NS-398 on carcinoma cells may be due to modulation of genes which regulate programmed cell death, cell proliferation and cell-cell communication. Additional studies are underway to determine the biological function of the novel genes that were identified.

Target gene search for the metal-responsive transcription factor MTF-1

Activation of genes by heavy metals, notably zinc, cadmium and copper, depends on MTF-1, a unique zinc finger transcription factor conserved from insects to human. Knockout of MTF-1 in the mouse results in embryonic lethality due to liver decay, while knockout of its best characterized target genes, the stress-inducible metallothionein genes I and II, is viable, suggesting additional target genes of MTF-1. Here we report on a multi-pronged search for potential target genes of MTF-1, including microarray screening, SABRE selective amplification, a computer search for MREs (DNA-binding sites of MTF-1) and transfection of reporter genes driven by candidate gene promoters. Some new candidate target genes emerged, including those encoding alpha-fetoprotein, the liver-enriched transcription factor C/EBPalpha and tear lipocalin/von Ebner's gland protein, all of which have a role in toxicity/the cell stress response. In contrast, expression of other cell stress-associated genes, such as those for superoxide dismutases, thioredoxin and heat shock proteins, do not appear to be affected by loss of MTF-1. Our experiments have also exposed some problems with target gene searches. First, finding the optimal time window for detecting MTF-1 target genes in a lethal phenotype of rapid liver decay proved problematical: 12.5-day-old mouse embryos (stage E12.5) yielded hardly any differentially expressed genes, whereas at stage 13.0 reduced expression of secretory liver proteins probably reflected the onset of liver decay, i.e. a secondary effect. Likewise, up-regulation of some proliferation-associated genes may also just reflect responses to the concomitant loss of hepatocytes. Another sobering finding concerns gamma-glutamylcysteine synthetase(hc) (gamma-GCS(hc)), which controls synthesis of the antioxidant glutathione and which was previously suggested to be a target gene contributing to the lethal phenotype in MTF-1 knockout mice. gamma-GCS(hc) mRNA is reduced at the onset of liver decay but MTF-1 null mutant embryos manage to maintain a very high glutathione level until shortly before that stage, perhaps in an attempt to compensate for low expression of metallothioneins, which also have a role as antioxidants.

Variability in gene expression patterns of Ewing tumor cell lines differing in EWS-FLI1 fusion type

Type 1 and type 2 EWS-FLI1 fusion products result from variation in breakpoint locations arising from the t(11;22)(q24;q12) recurrent chromosomal translocation in Ewing's sarcoma family tumors (EFT). Previously, studies from our institution (updated in the present communication at a median follow-up of more than 6 years) and others suggested a prognostic difference for EFT patients with localized disease depending on the type of EWS-FLI1 fusion present in the tumor. It has been suggested that the observed clinical discrepancies result from different transactivation potentials of the various EWS-FLI1 fusion proteins. In an attempt to identify genes whose expression levels are differentially modulated by structurally different EWS-FLI1 transcription factors, we have used two related PCR-based subtractive approaches, cDNA representational difference analysis (cDNA-RDA) and linker-capture subtraction (LCS) to compare transcript representations in cDNA pools of type 1 versus type 2 EFT cell lines. About 800 clones obtained by the two approaches were analyzed by dot blot hybridization to cDNA pools. Eighty-six clones showing the highest variability in signal intensities on the dot blots were further hybridized to individual EFT cell line RNAs on Northern blots, and four of them were additionally studied by real-time quantitative PCR (RTQ-PCR). Although interindividual variations in gene expression patterns in the range of one- to several-fold were observed, no correlation to specific EWS-FLI1 fusion types could be identified. Among the genes differentially expressed in individual EFT cell lines are several previously implicated in tumor growth, invasion, and metastasis. Although our data may have revealed candidate genes whose composite expression pattern may be relevant for the biology of individual EFT, they do not support a role of distinct EWS-FLI1 fusion types for EFT prognosis based on different transactivation potentials.

Isolation of binary species-specific PCR-based markers and their value for diagnostic applications

Representational difference analysis (RDA), a technique for the isolation of differences between highly similar complex genomes, was employed for isolation of species-specific markers. These markers can be easily adapted for a high throughput PCR-based assay in which multiple specimens can be simultaneously identified based on the presence/absence of amplification products. One of the important features of RDA performed on genomes of different species (interspecific RDA) is its ability to preferentially isolate families of repetitive sequences that are unique to one of the compared genomes. Such families of repetitive DNA are homoplasy-free characters that can be used for cost-efficient, mass identification of specimens in a variety of situations ranging from mark-recapture studies to screenings of egg or larval stages.

Genome changes and gene expression in human solid tumors

Genome-wide analysis techniques such as chromosome painting, comparative genomic hybridization, representational difference analysis, restriction landmark genome scanning and high-throughput analysis of LOH are now accelerating high-resolution genome aberration localization in human tumors. These techniques are complemented by procedures for detection of differentially expressed genes such as differential display, nucleic acid subtraction, serial analysis of gene expression and expression microarray analysis. These efforts are enabled by work from the human genome program in physical map development, cDNA library production/sequencing and in genome sequencing. This review covers several commonly used large-scale genome and gene expression analysis techniques, outlines genomic approaches to gene discovery and summarizes information that has come from large-scale analyses of human solid tumors.

Gene expression in hepatomas

Gene expression changes in accordance with cell growth, differentiation and carcinogenesis. To elucidate the molecular mechanisms for hepatocarcinogenesis as well as maintenance of normal hepatocytes, it is important to identify the genes that have altered expression with carcinogenesis. We established a new and efficient cDNA subtraction method via two cDNA populations. By using this method along with rat hepatomas made by the Soh-Farber protocol, we identified a number of genes, some of which are activated in hepatocellular carcinoma (HCC). These genes include ones which code for a transcription factor and a metabolic enzyme. One particular gene can be used as a tumour marker. Our method is beneficial for the isolation of a wide range of HCC-related genes in rats which, in turn, enables easy identification of their human counterparts. In this review, we describe details of our method and the isolated genes. We also briefly describe transcription factors in the liver.

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 differentially expressed genes upon immunoglobulin class switching by a subtractive hybridization method using uracil DNA glycosylase

Immunoglobulin class switch recombination enables B lymphocytes to sequentially express antibodies that have identical specificities but that differ in class and effector function. Although several cis elements required for class switch recombination have been identified, few trans -acting factors which are directly related to class switching have been found. Previously we have developed an efficient in vitro class switching system using a cell line, CH12F3-2. To clarify the molecular mechanism of class switching, we intended to isolate genes induced with interleukin (IL)-4, transforming growth factor (TGF)-beta and CD40L using the in vitro class switching system. For that purpose, an improved method for making subtracted cDNA libraries, using uracil DNA glycosylase, has been developed. This method can overcome a general problem of subtraction, that rare cDNAs are easily lost. This new subtraction method was applied to the CH12F3-2 switching system to isolate genes induced by stimulations with IL-4, TGF-beta and CD40L, and cDNAs encoding deiodinase 1 and SS32, an alternatively spliced form of the muscle LIM protein, were obtained. Their expression patterns in response to various combinations of stimuli and the time courses of the induction supported the usefulness of this method.

A systematic molecular genetic approach to study mammalian germline development

It is difficult to study gene expression in mammalian embryonic germ cells as PGCs constitute only a minor proportion of the mouse embryo. We have overcome this problem by using a novel combination of established molecular and transgenic approaches. A line of mice has been generated in which the cells of the germ lineage express the beta-galactosidase reporter gene during embryogenesis. Using this line, germ cells have been purified to near homogeneity from embryos at discrete stages during germline development by use of a stain for beta-gal activity and a fluorescence activated cell sorter. Subsequently, cDNA libraries have been constructed from each germ cell population using a modified lone-linker PCR strategy. These combined cDNA libraries represent genes expressed in PGCs during mammalian germline development. To facilitate a molecular genetic approach to studying mammalian germline development, these cDNA libraries will be pooled to form an arrayed, addressed reference embryonic germ cell cDNA library. In parallel with large-scale cDNA sequencing efforts; genes that are differentially expressed in germ cells will be identified by screening the reference library with probes generated by subtractive hybridization. Complementary DNAs identified using this approach will be analyzed by sequencing, database comparison, genomic mapping and in situ hybridization to ascertain the potential functional importance of each gene to germline development. In addition to providing a wealth of novel information regarding patterns of gene expression during mammalian germline development, these results will form the basis for future experiments to determine the function of these genes in this process.

A high throughput screening for rarely transcribed differentially expressed genes

A novel method combining elements of suppression subtractive hybridization with high throughput differential screening permits the efficient and rapid cloning of rarely transcribed differentially expressed genes. The experimental strategy virtually excludes the possibility of isolating false positive clones. The potential of the method is demonstrated by the isolation of 625 differentially expressed cDNAs from the metastatic adenocarcinoma cell line Bsp73-ASML when subtracted from its non-metastatic counterpart Bsp73-1AS. Northern analysis of 72 randomly selected clones demonstrated that 68 were differentially expressed with respect to Bsp73-ASML, indicating a true positive rate of 94%. Additionally, a large proportion of these clones represented rare transcripts as determined by the exposure time required to detect a signal. Sequence data indicated that of the 625 clones obtained, 92 clones scored perfect or near perfect matches with already known genes. Two hundred and eighty one clones scored between 60 and 95% homology to known human and mouse genes, whereas 252 clones scored no match with any sequences in the public databases. The method we describe is ideally suited whenever subtle changes in gene expression profiles need to be determined.

Selective amplification via biotin- and restriction-mediated enrichment (SABRE), a novel selective amplification procedure for detection of differentially expressed mRNAs

We present a novel subtractive enrichment protocol for the identification of differentially expressed mRNA species. This procedure, called SABRE (selective amplification via biotin- and restriction-mediated enrichment), uses selective streptavidin-biotin affinity and restriction enzyme site reconstitution to enrich for cDNA species more abundant in one population than in another. Analysis of liver cDNA from a mouse strain expressing the neomycin resistance gene demonstrated that this procedure is capable of identifying species present in one population but absent from another. Furthermore, experiments to identify genes with circadian expression patterns in mouse liver demonstrated that SABRE is capable of detecting even modest 2- to 10-fold differences in accumulation of moderately rare mRNA species, representing as little as 0.03% of total mRNA. These experiments identified the gene encoding coumarin 7-hydroxylase as displaying circadian expression in mouse liver.

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.

Subtractive cloning: past, present, and future

Subtractive cloning is a powerful technique for isolating genes expressed or present in one cell population but not in another. This method and a related one termed positive selection have their origins in nucleic acid reassociation techniques. We discuss the history of subtractive techniques, and fundamental information about the nucleic acid composition of cells that came out of reassociation analyses. We then explore current techniques for subtractive cloning and positive selection, discussing the merits of each. These techniques include cDNA library-based techniques and PCR-based techniques. Finally, we briefly discuss the future of subtractive cloning and new approaches that may augment or supersede current methods.