Computerized polymorphic marker identification: Experimental validation and a predicted human polymorphism catalog

A computational system for the prediction of polymorphic loci directly and efficiently from human genomic sequence was developed and verified. A suite of programs, collectively called pompous (polymorphic marker prediction of ubiquitous simple sequences) detects tandem repeats ranging from dinucleotides up to 250 mers, scores them according to predicted level of polymorphism, and designs appropriate flanking primers for PCR amplification. This approach was validated on an approximately 750-kilobase region of human chromosome 3p21.3, involved in lung and breast carcinoma homozygous deletions. Target DNA from 36 paired B lymphoblastoid and lung cancer lines was amplified and allelotyped for 33 loci predicted by pompous to be variable in repeat size. We found that among those 36 predominately Caucasian individuals 22 of the 33 (67%) predicted loci were polymorphic with an average heterozygosity of 0.42. Allele loss in this region was found in 27/36 (75%) of the tumor lines using these markers. pompous provides the genetic researcher with an additional tool for the rapid and efficient identification of polymorphic markers, and through a World Wide Web site, investigators can use pompous to identify polymorphic markers for their research. A catalog of 13,261 potential polymorphic markers and associated primer sets has been created from the analysis of 141,779,504 base pairs of human genomic sequence in GenBank. This data is available on our Web site (pompous.swmed.edu) and will be updated periodically as GenBank is expanded and algorithm accuracy is improved.

Identification of lumichrome as a Sinorhizobium enhancer of alfalfa root respiration and shoot growth

Sinorhizobium meliloti bacteria produce a signal molecule that enhances root respiration in alfalfa (Medicago sativa L.) and also triggers a compensatory increase in whole-plant net carbon assimilation. Nuclear magnetic resonance, mass spectrometry, and ultraviolet–visible absorption identify the enhancer as lumichrome, a common breakdown product of riboflavin. Treating alfalfa roots with 3 nM lumichrome increased root respiration 21% (P < 0.05) within 48 h. A closely linked increase in net carbon assimilation by the shoot compensated for the enhanced root respiration. For example, applying 5 nM lumichrome to young alfalfa roots increased plant growth by 8% (P < 0.05) after 12 days. Soaking alfalfa seeds in 5 nM lumichrome before germination increased growth by 18% (P < 0.01) over the same period. In both cases, significant growth enhancement (P < 0.05) was evident only in the shoot. S. meliloti requires exogenous CO2 for growth and may benefit directly from the enhanced root respiration that is triggered by lumichrome. Thus Sinorhizobium–alfalfa associations, which ultimately form symbiotic N2-reducing root nodules, may be favored at an early developmental stage by lumichrome, a previously unrecognized mutualistic signal. The rapid degradation of riboflavin to lumichrome under many physiological conditions and the prevalence of riboflavin release by rhizosphere bacteria suggest that events demonstrated here in the S. meliloti–alfalfa association may be widely important across many plant–microbe interactions.

Identification of genes required for Drosophila eye development using a phenotypic enhancer-trap

A novel method of P-element mutagenesis is described for the isolation of mutants affecting the development of the Drosophila compound eye. It exploits the interaction between the Bride of Sevenless (Boss) ligand and the Sevenless (Sev) receptor tyrosine kinase that triggers the formation of the UV-sensitive photoreceptor neuron, R7. Transposition of a boss cDNA transgene, in an otherwise boss mutant background, was used as a “phenotypic trap” in live flies to identify enhancers expressed during a narrow time window in eye development. Using a rapid behavioral screen, more than 400,000 flies were tested for restoration of R7. Some 1,800 R7-containing flies were identified. Among these, 21 independent insertions with expression of the boss reporter gene in the R8 cell were identified by a external eye morphology and staining with an antibody against Boss. Among 900 lines with expression of the boss reporter gene in multiple cells assessed for homozygous mutant phenotypes, insertions in the marbles, glass, gap1, and fasciclin II genes were isolated. This phenotypic enhancer-trap facilitates (i) the isolation of enhancer-traps with a specific expression pattern, and (ii) the recovery of mutants disrupting development of specific tissues. Because the temporal and tissue specificity of the phenotypic trap is dependent on the choice of the marker used, this approach can be extended to other tissues and developmental stages.

Rapid down-regulation of mammalian Period genes during behavioral resetting of the circadian clock

The pervasive role of circadian clocks in regulating physiology and behavior is widely recognized. Their adaptive value is their ability to be entrained by environmental cues such that the internal circadian phase is a reliable predictor of solar time. In mammals, both light and nonphotic behavioral cues can entrain the principal oscillator of the hypothalamic suprachiasmatic nuclei (SCN). However, although light can advance or delay the clock during circadian night, behavioral events trigger phase advances during the subjective day, when the clock is insensitive to light. The recent identification of Period (Per) genes in mammals, homologues of dperiod, which encodes a core element of the circadian clockwork in Drosophila, now provides the opportunity to explain circadian timing and entrainment at a molecular level. In mice, expression of mPer1 and mPer2 in the SCN is rhythmic and acutely up-regulated by light. Moreover, the temporal relations between mRNA and protein cycles are consistent with a clock based on a transcriptional/translational feedback loop. Here we describe circadian oscillations of Per1 and Per2 in the SCN of the Syrian hamster, showing that PER1 protein and mRNA cycles again behave in a manner consistent with a negative-feedback oscillator. Furthermore, we demonstrate that nonphotic resetting has the opposite effect to light: acutely down-regulating these genes. Their sensitivity to nonphotic resetting cues supports their proposed role as core elements of the circadian oscillator. Moreover, this study provides an explanation at the molecular level for the contrasting but convergent effects of photic and nonphotic cues on the clock.

Reciprocal subtraction differential RNA display: An efficient and rapid procedure for isolating differentially expressed gene sequences

A reciprocal subtraction differential RNA display (RSDD) approach has been developed that permits the rapid and efficient identification and cloning of both abundant and rare differentially expressed genes. RSDD comprises reciprocal subtraction of cDNA libraries followed by differential RNA display. The RSDD strategy was applied to analyze the gene expression alterations resulting during cancer progression as adenovirus-transformed rodent cells developed an aggressive transformed state, as documented by elevated anchorage-independence and enhanced in vivo oncogenesis in nude mice. This approach resulted in the identification and cloning of both known and a high proportion (>65%) of unknown sequences, including cDNAs displaying elevated expression as a function of progression (progression-elevated gene) and cDNAs displaying suppressed expression as a function of progression (progression-suppressed gene). Sixteen differentially expressed genes, including five unknown progression-elevated genes and six unknown progression-suppressed genes, have been characterized. The RSDD scheme should find wide application for the effective detection and isolation of differentially expressed genes.

Identification of fission yeast nuclear markers using random polypeptide fusions with green fluorescent protein

We describe a method for identifying genes encoding proteins with stereospecific intracellular localizations in the fission yeast Schizosaccharomyces pombe. Yeast are transformed with a gene library in which S. pombe genomic sequences are fused to the gene encoding the Aequorea victoria green fluorescent protein (GFP), and intracellular localizations are subsequently identified by rapid fluorescence screening in vivo. In a model application of these methods to the fission yeast nucleus, we have identified several novel genes whose products are found in specific nuclear regions, including chromatin, the nucleolus, and the mitotic spindle, and sequence similarities between some of these genes and previously identified genes encoding nuclear proteins have validated the approach. These methods will be useful in identifying additional components of the S. pombe nucleus, and further extensions of this approach should also be applicable to a more comprehensive identification of the elements of intracellular architecture in fission yeast.

A novel approach for the identification of protein–protein interaction with integral membrane proteins

Protein–protein interaction plays a major role in all biological processes. The currently available genetic methods such as the two-hybrid system and the protein recruitment system are relatively limited in their ability to identify interactions with integral membrane proteins. Here we describe the development of a reverse Ras recruitment system (reverse RRS), in which the bait used encodes a membrane protein. The bait is expressed in its natural environment, the membrane, whereas the protein partner (the prey) is fused to a cytoplasmic Ras mutant. Protein–protein interaction between the proteins encoded by the prey and the bait results in Ras membrane translocation and activation of a viability pathway in yeast. We devised the expression of the bait and prey proteins under the control of dual distinct inducible promoters, thus enabling a rapid selection of transformants in which growth is attributed solely to specific protein–protein interaction. The reverse RRS approach greatly extends the usefulness of the protein recruitment systems and the use of integral membrane proteins as baits. The system serves as an attractive approach to explore novel protein–protein interactions with high specificity and selectivity, where other methods fail.

Identification of interdomain sequences promoting the intronless evolution of a bacterial protein family.

In the evolution of eukaryotic genes, introns are believed to have played a major role in increasing the probability of favorable duplication events, chance recombinations, and exon shuffling resulting in functional hybrid proteins. As a rule, prokaryotic genes lack introns, and the examples of prokaryotic introns described do not seem to have contributed to gene evolution by exon shuffling. Still, certain protein families in modern bacteria evolve rapidly by recombination of genes, duplication of functional domains, and as shown for protein PAB of the anaerobic bacterial species Peptostreptococcus magnus, by the shuffling of an albumin-binding protein module from group C and G streptococci. Characterization of a protein PAB-related gene in a P. magnus strain with less albumin-binding activity revealed that the shuffled module was missing. Based on this fact and observations made when comparing gene sequences of this family of bacterial surface proteins interacting with albumin and/or immunoglobulin, a model is presented that can explain how this rapid intronless evolution takes place. A new kind of genetic element is introduced: the recer sequence promoting interdomain, in frame recombination and acting as a structure-less flexibility-promoting spacer in the corresponding protein. The data presented also suggest that antibiotics could represent the selective pressure behind the shuffling of protein modules in P. magnus, a member of the indigenous bacterial flora.

Complementation of an Escherichia coli adhE mutant by the Entamoeba histolytica EhADH2 gene provides a method for the identification of new antiamebic drugs.

The pathogenic protozoan parasite Entamoeba histolytica, the cause of amebic dysentery and amebic liver abscess, is an obligate anaerobe, and derives energy from the fermentation of glucose to ethanol with pyruvate and acetyl coenzyme A as intermediates. We have isolated EhADH2, a key enzyme in this pathway, that is a NAD+- and Fe2+-dependent bifunctional enzyme with acetaldehyde dehydrogenase and alcohol dehydrogenase activities. EhADH2 is the only known eukaryotic member of a newly defined family of prokaryotic multifunctional enzymes, which includes the Escherichia coli AdhE enzyme, an enzyme required for anaerobic growth of E. coli. Because of the critical role of EhADH2 in the amebic fermentation pathway and the lack of known eukaryotic homologues of the EhADH2 enzyme, EhADH2 represents a potential target for antiamebic chemotherapy. However, screening of compounds for antiamebic activity is hampered by the cost of large scale growth of Ent. histolytica, and difficulties in quantitating drug efficacy in vitro. To approach this problem, we expressed the EhADH2 gene in a mutant strain of E. coli carrying a deletion of the adhE gene. Expression of EhADH2 restored the ability of the mutant E. coli strain to grow under anaerobic conditions. By screening compounds for the ability to inhibit the anaerobic growth of the E. coli/EhADH2 strain, we have developed a rapid assay for identifying compounds with anti-EhADH2 activity. Using bacteria to bypass the need for parasite culture in the initial screening process for anti-parasitic agents could greatly simplify and reduce the cost of identifying new therapeutic agents effective against parasitic diseases.

Mutagenesis and gene identification in Dictyostelium by shotgun antisense.

We have developed a mutagenesis technique that uses antisense cDNA to identify genes required for development in Dictyostelium discoideum. We transformed Dictyostelium cells with a cDNA library made from the mRNA of vegetative and developing cells. The cDNA was cloned in an antisense orientation immediately downstream of a vegetative promoter, so that in transformed cells the promoter will drive the synthesis of an antisense RNA transcript. We find that individual transformants typically contain one or occasionally two antisense cDNAs. Using this mutagenesis technique, we have generated mutants that fail to aggregate, aggregate but fail to form fruiting bodies, or aggregate but form abnormal fruiting bodies. The individual cDNA molecules from the mutants were identified and cloned using PCR. Initial sequence analysis of the PCR products from 35 mutants has identified six novel Dictyostelium genes, each from a transformant with one antisense cDNA. When the PCR-isolated antisense cDNAs were ligated into the antisense vector and the resulting constructs transformed into cells, the phenotypes of the transformed cells matched those of the original mutants from which each cDNA was obtained. We made homologous recombinant gene disruption transformants for three of the novel genes, in each case generating mutants with phenotypes indistinguishable from those of the original antisense transformants. Shotgun antisense thus is a rapid way to identify genes in Dictyostelium and possibly other organisms.

Identification of the promoter of the mouse obese gene.

Primer extension and RACE (rapid amplification of cDNA ends) assays were used to identify and sequence the 5' terminus of mouse ob mRNA. This sequence was used to obtain a recombinant bacteriophage containing the first exon of the encoding gene. DNA sequence analysis of the region immediately upstream of the first exon of the mouse ob gene revealed DNA sequences corresponding to presumptive cis-regulatory elements. A canonical TATA box was observed 30-34 base pairs upstream from the start site of transcription and a putative binding site for members of the C/EBP family of transcription factors was identified immediately upstream from the TATA box. Nuclear extracts prepared from primary adipocytes contained a DNA binding activity capable of avid and specific interaction with the putative C/EBP response element; antibodies to C/EBP alpha neutralized the DNA binding activity present in adipocyte nuclear extracts. When linked to a firefly luciferase reporter and transfected into primary adipocytes, the presumptive promoter of the mouse ob gene facilitated luciferase expression. When transfected into HepG2 cells, which lack C/EBP alpha, the mouse ob promoter was only weakly active. Supplementation of C/EBP alpha by cotransfection with a C/EBP alpha expression vector markedly stimulated luciferase expression. Finally, an ob promoter variant mutated at the C/EBP response element was inactive in both primary adipocytes and HepG2 cells. These observations provide evidence for identification of a functional promoter capable of directing expression of the mouse ob gene.

Sequence scanning: A method for rapid sequence acquisition from large-fragment DNA clones.

A strategy of "sequence scanning" is proposed for rapid acquisition of sequence from clones such as bacteriophage P1 clones, cosmids, or yeast artificial chromosomes. The approach makes use of a special vector, called LambdaScan, that reliably yields subclones with inserts in the size range 8-12 kb. A number of subclones, typically 96 or 192, are chosen at random, and the ends of the inserts are sequenced using vector-specific primers. Then long-range spectrum PCR is used to order and orient the clones. This combination of shotgun and directed sequencing results in a high-resolution physical map suitable for the identification of coding regions or for comparison of sequence organization among genomes. Computer simulations indicate that, for a target clone of 100 kb, the scanning of 192 subclones with sequencing reads as short as 350 bp results in an approximate ratio of 1:2:1 of regions of double-stranded sequence, single-stranded sequence, and gaps. Longer sequencing reads tip the ratio strongly toward increased double-stranded sequence.

Analysis of RNA-binding proteins by in vitro genetic selection: identification of an amino acid residue important for locking U1A onto its RNA target.

An in vitro genetic system was developed as a rapid means for studying the specificity determinants of RNA-binding proteins. This system was used to investigate the origin of the RNA-binding specificity of the mammalian spliceosomal protein U1A. The U1A domain responsible for binding to U1 small nuclear RNA was locally mutagenized and displayed as a combinatorial library on filamentous bacteriophage. Affinity selection identified four U1A residues in the mutagenized region that are important for specific binding to U1 hairpin II. One of these residues (Leu-49) disproportionately affects the rates of binding and release and appears to play a critical role in locking the protein onto the RNA. Interestingly, a protein variant that binds more tightly than U1A emerged during the selection, showing that the affinity of U1A for U1 RNA has not been optimized during evolution.

Identification of human granzyme B promoter regulatory elements interacting with activated T-cell-specific proteins: implication of Ikaros and CBF binding sites in promoter activation.

Granzyme B serine protease is found in the granules of activated cytotoxic T cells and in natural and lymphokine-activated killer cells. This protease plays a critical role in the rapid induction of target cell DNA fragmentation. The DNA regulatory elements that are responsible for the specificity of granzyme B gene transcription in activated T-cells reside between nt -148 and +60 (relative to the transcription start point at +1) of the human granzyme B gene promoter. This region contains binding sites for the transcription factors Ikaros, CBF, Ets, and AP-1. Mutational analysis of the human granzyme B promoter reveals that the Ikaros binding site (-143 to -114) and the AP-1/CBF binding site (-103 to -77) are essential for the activation of transcription in phytohemagglutinin-activated peripheral blood lymphocytes, whereas mutation of the Ets binding site does not affect promoter activity in these cells.

Identification of human brain regions underlying responses to resistive inspiratory loading with functional magnetic resonance imaging.

Compensatory ventilatory responses to increased inspiratory loading are essential for adequate breathing regulation in a number of pulmonary diseases; however, the human brain sites mediating such responses are unknown. Midsagittal and axial images were acquired in 11 healthy volunteers during unloaded and loaded (30 cmH2O; 1 cmH2O = 98 Pa) inspiratory breathing, by using functional magnetic resonance imaging (fMRI) strategies (1.5-tesla MR; repetition time, 72 msec; echo time, 45 msec; flip angle, 30 degrees; field of view, 26 cm; slice thickness, 5 mm; number of excitations, 1; matrix, 128 x 256). Digital image subtractions and region of interest analyses revealed significantly increased fMRI signal intensity in discrete areas of the ventral and dorsal pons, interpeduncular nucleus, basal forebrain, putamen, and cerebellar regions. Upon load withdrawal, certain regions displayed a rapid fMRI signal off-transient, while in others, a slower fMRI signal decay emerged. Sustained loading elicited slow decreases in fMRI signal across activated regions, while second application of an identical load resulted in smaller signal increases compared to initial signal responses (P < 0.001). A moderate inspiratory load is associated with consistent regional activation of discrete brain locations; certain of these regions have been implicated in mediation of loaded breathing in animal models. We speculate that temporal changes in fMRI signal may indicate respiratory after-discharge and/or habituation phenomena.