Real-time detection of DNA hybridization and melting on oligonucleotide arrays by using optical wave guides.

The challenge of the Human Genome Project is to increase the rate of DNA sequence acquisition by two orders of magnitude to complete sequencing of the human genome by the year 2000. The present work describes a rapid detection method using a two-dimensional optical wave guide that allows measurement of real-time binding or melting of a light-scattering label on a DNA array. A particulate label on the target DNA acts as a light-scattering source when illuminated by the evanescent wave of the wave guide and only the label bound to the surface generates a signal. Imaging/visual examination of the scattered light permits interrogation of the entire array simultaneously. Hybridization specificity is equivalent to that obtained with a conventional system using autoradiography. Wave guide melting curves are consistent with those obtained in the liquid phase and single-base discrimination is facile. Dilution experiments showed an apparent lower limit of detection at 0.4 nM oligonucleotide. This performance is comparable to the best currently known fluorescence-based systems. In addition, wave guide detection allows manipulation of hybridization stringency during detection and thereby reduces DNA chip complexity. It is anticipated that this methodology will provide a powerful tool for diagnostic applications that require rapid cost-effective detection of variations from known sequences.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Subnuclear localization of the active variant surface glycoprotein gene expression site in Trypanosoma brucei

In Trypanosoma brucei, transcription by RNA polymerase II and 5′ capping of messenger RNA are uncoupled: a capped spliced leader is trans spliced to every RNA. This decoupling makes it possible to have protein-coding gene transcription driven by RNA polymerase I. Indeed, indirect evidence suggests that the genes for the major surface glycoproteins, variant surface glycoproteins (VSGs) in bloodstream-form trypanosomes, are transcribed by RNA polymerase I. In a single trypanosome, only one VSG expression site is maximally transcribed at any one time, and it has been speculated that transcription takes place at a unique site within the nucleus, perhaps in the nucleolus. We tested this by using fluorescence in situ hybridization. With probes that cover about 50 kb of the active 221 expression site, we detected nuclear transcripts of this site in a single fluorescent spot, which did not colocalize with the nucleolus. Analysis of marker gene-tagged active expression site DNA by fluorescent DNA in situ hybridization confirmed the absence of association with the nucleolus. Even an active expression site in which the promoter had been replaced by an rDNA promoter did not colocalize with the nulceolus. As expected, marker genes inserted in the rDNA array predominantly colocalize with the nucleolus, whereas the tubulin gene arrays do not. We conclude that transcription of the active VSG expression site does not take place in the nucleolus.

Subtraction hybridization identifies a transformation progression-associated gene PEG-3 with sequence homology to a growth arrest and DNA damage-inducible gene

Cancer is a progressive multigenic disorder characterized by defined changes in the transformed phenotype that culminates in metastatic disease. Determining the molecular basis of progression should lead to new opportunities for improved diagnostic and therapeutic modalities. Through the use of subtraction hybridization, a gene associated with transformation progression in virus- and oncogene-transformed rat embryo cells, progression elevated gene-3 (PEG-3), has been cloned. PEG-3 shares significant nucleotide and amino acid sequence homology with the hamster growth arrest and DNA damage-inducible gene gadd34 and a homologous murine gene, MyD116, that is induced during induction of terminal differentiation by interleukin-6 in murine myeloid leukemia cells. PEG-3 expression is elevated in rodent cells displaying a progressed-transformed phenotype and in rodent cells transformed by various oncogenes, including Ha-ras, v-src, mutant type 5 adenovirus (Ad5), and human papilloma virus type 18. The PEG-3 gene is transcriptionally activated in rodent cells, as is gadd34 and MyD116, after treatment with DNA damaging agents, including methyl methanesulfonate and γ-irradiation. In contrast, only PEG-3 is transcriptionally active in rodent cells displaying a progressed phenotype. Although transfection of PEG-3 into normal and Ad5-transformed cells only marginally suppresses colony formation, stable overexpression of PEG-3 in Ad5-transformed rat embryo cells elicits the progression phenotype. These results indicate that PEG-3 is a new member of the gadd and MyD gene family with similar yet distinct properties and this gene may directly contribute to the transformation progression phenotype. Moreover, these studies support the hypothesis that constitutive expression of a DNA damage response may mediate cancer progression.

Exploring drug-induced alterations in gene expression in Mycobacterium tuberculosis by microarray hybridization

Tuberculosis is a chronic infectious disease that is transmitted by cough-propelled droplets that carry the etiologic bacterium, Mycobacterium tuberculosis. Although currently available drugs kill most isolates of M. tuberculosis, strains resistant to each of these have emerged, and multiply resistant strains are increasingly widespread. The growing problem of drug resistance combined with a global incidence of seven million new cases per year underscore the urgent need for new antituberculosis therapies. The recent publication of the complete sequence of the M. tuberculosis genome has made possible, for the first time, a comprehensive genomic approach to the biology of this organism and to the drug discovery process. We used a DNA microarray containing 97% of the ORFs predicted from this sequence to monitor changes in M. tuberculosis gene expression in response to the antituberculous drug isoniazid. Here we show that isoniazid induced several genes that encode proteins physiologically relevant to the drug’s mode of action, including an operonic cluster of five genes encoding type II fatty acid synthase enzymes and fbpC, which encodes trehalose dimycolyl transferase. Other genes, not apparently within directly affected biosynthetic pathways, also were induced. These genes, efpA, fadE23, fadE24, and ahpC, likely mediate processes that are linked to the toxic consequences of the drug. Insights gained from this approach may define new drug targets and suggest new methods for identifying compounds that inhibit those targets.

Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

Keratinocyte growth factor (KGF) is a member of the fibroblast growth factor family. Portions of the gene encoding KGF were amplified during primate evolution and are present in multiple nonprocessed copies in the human genome. Nucleotide analysis of a representative sampling of these KGF-like sequences indicated that they were at least 95% identical to corresponding regions of the KGF gene. To localize these sequences to specific chromosomal sites in human and higher primates, we used fluorescence in situ hybridization. In human, using a cosmid probe encoding KGF exon 1, we assigned the location of the KGF gene to chromosome 15q15–21.1. In addition, copies of KGF-like sequences hybridizing only with a cosmid probe encoding exons 2 and 3 were localized to dispersed sites on chromosome 2q21, 9p11, 9q12–13, 18p11, 18q11, 21q11, and 21q21.1. The distribution of KGF-like sequences suggests a role for alphoid DNA in their amplification and dispersion. In chimpanzee, KGF-like sequences were observed at five chromosomal sites, which were each homologous to sites in human, while in gorilla, a subset of four of these homologous sites was identified; in orangutan two sites were identified, while gibbon exhibited only a single site. The chromosomal localization of KGF sequences in human and great ape genomes indicates that amplification and dispersion occurred in multiple discrete steps, with initial KGF gene duplication and dispersion taking place in gibbon and involving loci corresponding to human chromosomes 15 and 21. These findings support the concept of a closer evolutionary relationship of human and chimpanzee and a possible selective pressure for such dispersion during the evolution of higher primates.

Induction of inhibitory factor κBα mRNA in the central nervous system after peripheral lipopolysaccharide administration: An in situ hybridization histochemistry study in the rat

In this study we investigate the mRNA expression of inhibitory factor κBα (IκBα) in cells of the rat brain induced by an intraperitoneal (i.p.) injection of lipopolysaccharide (LPS). IκB controls the activity of nuclear factor κB, which regulates the transcription of many immune signal molecules. The detection of IκB induction, therefore, would reveal the extent and the cellular location of brain-derived immune molecules in response to peripheral immune challenges. Low levels of IκBα mRNA were found in the large blood vessels and in circumventricular organs (CVOs) of saline-injected control animals. After an i.p. LPS injection (2.5 mg/kg), dramatic induction of IκBα mRNA occurred in four spatio-temporal patterns. Induced signals were first detected at 0.5 hr in the lumen of large blood vessels and in blood vessels of the choroid plexus and CVOs. Second, at 1–2 hr, labeling dramatically increased in the CVOs and choroid plexus and spread to small vascular and glial cells throughout the entire brain; these responses peaked at 2 hr and declined thereafter. Third, cells of the meninges became activated at 2 hr and persisted until 12 hr after the LPS injection. Finally, only at 12 hr, induced signals were present in ventricular ependyma. Thus, IκBα mRNA is induced in brain after peripheral LPS injection, beginning in cells lining the blood side of the blood–brain barrier and progressing to cells inside brain. The spatiotemporal patterns suggest that cells of the blood–brain barrier synthesize immune signal molecules to activate cells inside the central nervous system in response to peripheral LPS. The cerebrospinal fluid appears to be a conduit for these signal molecules.

PCR-based subtractive hybridization and differences in gene content among strains of Helicobacter pylori

Genes that are characteristic of only certain strains of a bacterial species can be of great biologic interest. Here we describe a PCR-based subtractive hybridization method for efficiently detecting such DNAs and apply it to the gastric pathogen Helicobacter pylori. Eighteen DNAs specific to a monkey-colonizing strain (J166) were obtained by subtractive hybridization against an unrelated strain whose genome has been fully sequenced (26695). Seven J166-specific clones had no DNA sequence match to the 26695 genome, and 11 other clones were mixed, with adjacent patches that did and did not match any sequences in 26695. At the protein level, seven clones had homology to putative DNA restriction-modification enzymes, and two had homology to putative metabolic enzymes. Nine others had no database match with proteins of assigned function. PCR tests of 13 unrelated H. pylori strains by using primers specific for 12 subtracted clones and complementary Southern blot hybridizations indicated that these DNAs are highly polymorphic in the H. pylori population, with each strain yielding a different pattern of gene-specific PCR amplification. The search for polymorphic DNAs, as described here, should help identify previously unknown virulence genes in pathogens and provide new insights into microbial genetic diversity and evolution.

Yeast microarrays for genome wide parallel genetic and gene expression analysis

We have developed high-density DNA microarrays of yeast ORFs. These microarrays can monitor hybridization to ORFs for applications such as quantitative differential gene expression analysis and screening for sequence polymorphisms. Automated scripts retrieved sequence information from public databases to locate predicted ORFs and select appropriate primers for amplification. The primers were used to amplify yeast ORFs in 96-well plates, and the resulting products were arrayed using an automated micro arraying device. Arrays containing up to 2,479 yeast ORFs were printed on a single slide. The hybridization of fluorescently labeled samples to the array were detected and quantitated with a laser confocal scanning microscope. Applications of the microarrays are shown for genetic and gene expression analysis at the whole genome level.

Identification of transforming growth factor-β- regulated genes in Caenorhabditis elegans by differential hybridization of arrayed cDNAs

Members of the transforming growth factor-β family play critical roles in body patterning, in both vertebrates and invertebrates. One transforming growth factor-β-related gene, dbl-1, has been shown to regulate body length and male ray patterning in Caenorhabditis elegans. We screened arrayed cDNAs to identify downstream target genes for the DBL-1 signaling by using differential hybridization. C. elegans cDNAs representing 7,584 independent genes were arrayed on a nylon membrane at high density and hybridized with 33P-labeled DNA probes synthesized from the mRNAs of wild-type, dbl-1, sma-2, and lon-2 worms. Signals for all the spots representing hybridized DNA were quantified and compared among strains. The screening identified 22 and 2 clones, which were positively and negatively regulated, respectively, by the DBL-1 signal. Northern hybridization confirmed the expression profiles of most of the clones, indicating good reliability of the differential hybridization using arrayed cDNAs. In situ hybridization analysis revealed the spatial and temporal expression patterns of each clone and showed that at least four genes, including the gene for the type I receptor for DBL-1, sma-6, were transcriptionally regulated by the DBL-1 signal.

tRNA Nuclear Export in Saccharomyces cerevisiae: In Situ Hybridization Analysis

To understand the factors specifically affecting tRNA nuclear export, we adapted in situ hybridization procedures to locate endogenous levels of individual tRNA families in wild-type and mutant yeast cells. Our studies of tRNAs encoded by genes lacking introns show that nucleoporin Nup116p affects both poly(A) RNA and tRNA export, whereas Nup159p affects only poly(A) RNA export. Los1p is similar to exportin-t, which facilitates vertebrate tRNA export. A los1 deletion mutation affects tRNA but not poly(A) RNA export. The data support the notion that Los1p and exportin-t are functional homologues. Because LOS1 is nonessential, tRNA export in vertebrate and yeast cells likely involves factors in addition to exportin-t. Mutation of RNA1, which encodes RanGAP, causes nuclear accumulation of tRNAs and poly(A) RNA. Many yeast mutants, including those with the rna1-1 mutation, affect both pre-tRNA splicing and RNA export. Our studies of the location of intron-containing pre-tRNAs in the rna1-1 mutant rule out the possibility that this results from tRNA export occurring before splicing. Our results also argue against inappropriate subnuclear compartmentalization causing defects in pre-tRNA splicing. Rather, the data support “feedback” of nucleus/cytosol exchange to the pre-tRNA splicing machinery.

Peptide nucleic acid pre-gel hybridization: An alternative to Southern hybridization

We have found that it is possible to use labeled peptide nucleic acid (PNA)-oligomers as probes in pre-gel hybridization experiments, as an alternative for Southern hybridization. In this technique, the PNA probe is hybridized to a denatured DNA sample at low ionic strength and the mixture is loaded directly on to an electrophoresis system for size separation. Ensuing gel electrophoresis separates the single-stranded DNA fragments by length. The neutral backbone of PNA allows for hybridization at low ionic strength and imparts very low mobility to excess PNA. Detection of the bound PNA is possible by direct fluorescence detection with capillary electrophoresis, or the DNA/PNA hybrids can be blotted onto a membrane and detected with standard chemiluminescent techniques. Efficient single bp discrimination was achieved routinely using both capillary and slab-gel electrophoresis.

Intranuclear diffusion and hybridization state of oligonucleotides measured by fluorescence correlation spectroscopy in living cells

Fluorescein-labeled oligodeoxynucleotides (oligos) were introduced into cultured rat myoblasts, and their molecular movements inside the nucleus were studied by fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). FCS revealed that a large fraction of both intranuclear oligo(dT) (43%) and oligo(dA) (77%) moves rapidly with a diffusion coefficient of 4 × 10−7 cm2/s. Interestingly, this rate of intranuclear oligo movement is similar to their diffusion rates measured in aqueous solution. In addition, we detected a large fraction (45%) of the intranuclear oligo(dT), but not oligo(dA), diffusing at slower rates (≤1 × 10−7 cm2/s). The amount of this slower-moving oligo(dT) was greatly reduced if the oligo(dT) was prehybridized in solution with (unlabeled) oligo(dA) prior to introduction to cells, presumably because the oligo(dT) was then unavailable for subsequent hybridization to endogenous poly(A) RNA. The FCS-measured diffusion rate for much of the slower oligo(dT) population approximated the diffusion rate in aqueous solution of oligo(dT) hybridized to a large polyadenylated RNA (1.0 × 10−7 cm2/s). Moreover, this intranuclear movement rate falls within the range of calculated diffusion rates for an average-sized heterogeneous nuclear ribonucleoprotein particle in aqueous solution. A subfraction of oligo(dT) (15%) moved over 10-fold more slowly, suggesting it was bound to very large macromolecular complexes. Average diffusion coefficients obtained from FRAP experiments were in agreement with the FCS data. These results demonstrate that oligos can move about within the nucleus at rates comparable to those in aqueous solution and further suggest that this is true for large ribonucleoprotein complexes as well.

Hybridization of single-stranded DNA targets to immobilized complementary DNA probes: comparison of hairpin versus linear capture probes

A microtiter-based assay system is described in which DNA hairpin probes with dangling ends and single-stranded, linear DNA probes were immobilized and compared based on their ability to capture single-strand target DNA. Hairpin probes consisted of a 16 bp duplex stem, linked by a T2-biotin·dT-T2 loop. The third base was a biotinylated uracil (UB) necessary for coupling to avidin coated microtiter wells. The capture region of the hairpin was a 3′ dangling end composed of either 16 or 32 bases. Fundamental parameters of the system, such as probe density and avidin adsorption capacity of the plates were characterized. The target DNA consisted of 65 bases whose 3′ end was complementary to the dangling end of the hairpin or to the linear probe sequence. The assay system was employed to measure the time dependence and thermodynamic stability of target hybridization with hairpin and linear probes. Target molecules were labeled with either a 5′-FITC, or radiolabeled with [γ-33P]ATP and captured by either linear or hairpin probes affixed to the solid support. Over the range of target concentrations from 10 to 640 pmol hybridization rates increased with increasing target concentration, but varied for the different probes examined. Hairpin probes displayed higher rates of hybridization and larger equilibrium amounts of captured targets than linear probes. At 25 and 45°C, rates of hybridization were better than twice as great for the hairpin compared with the linear capture probes. Hairpin–target complexes were also more thermodynamically stable. Binding free energies were evaluated from the observed equilibrium constants for complex formation. Results showed the order of stability of the probes to be: hairpins with 32 base dangling ends > hairpin probes with l6 base dangling ends > 16 base linear probes > 32 base linear probes. The physical characteristics of hairpins could offer substantial advantages as nucleic acid capture moieties in solid support based hybridization systems.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.