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.

Comparative analysis of ribonuclease P RNA using gene sequences from natural microbial populations reveals tertiary structural elements.

PCR amplification of template DNAs extracted from mixed, naturally occurring microbial populations, using oligonucleotide primers complementary to highly conserved sequences, was used to obtain a large collection of diverse RNase P RNA-encoding genes. An alignment of these sequences was used in a comparative analysis of RNase P RNA secondary and tertiary structure. The new sequences confirm the secondary structure model based on sequences from cultivated organisms (with minor alterations in helices P12 and P18), providing additional support for nearly every base pair. Analysis of sequence covariation using the entire RNase P RNA data set reveals elements of tertiary structure in the RNA; the third nucleotides (underlined) of the GNRA tetraloops L14 and L18 are seen to interact with adjacent Watson-Crick base pairs in helix P8, forming A:G/C or G:A/U base triples. These experiments demonstrate one way in which the enormous diversity of natural microbial populations can be used to elucidate molecular structure through comparative analysis.

Human von Willebrand factor gene sequences target expression to a subpopulation of endothelial cells in transgenic mice.

The present study was undertaken to define the 5' and 3' regulatory sequences of human von Willebrand factor gene that confer tissue-specific expression in vivo. Transgenic mice were generated bearing a chimeric construct that included 487 bp of 5' flanking sequence and the first exon fused in-frame to the Escherichia coli lacZ gene. In situ histochemical analyses in independent lines demonstrated that the von Willebrand factor promoter targeted expression of LacZ to a subpopulation of endothelial cells in the yolk sac and adult brain. LacZ activity was absent in the vascular beds of the spleen, lung, liver, kidney, testes, heart, and aorta, as well as in megakaryocytes. In contrast, in mice containing the lacZ gene targeted to the thrombomodulin locus, the 5-bromo-4-chloro-3-indolyl beta-D-galactopyranoside reaction product was detected throughout the vascular tree. These data highlight the existence of regional differences in endothelial cell gene regulation and suggest that the 733-bp von Willebrand factor promoter may be useful as a molecular marker to investigate endothelial cell diversity.

The TIGR Gene Indices: analysis of gene transcript sequences in highly sampled eukaryotic species

While genome sequencing projects are advancing rapidly, EST sequencing and analysis remains a primary research tool for the identification and categorization of gene sequences in a wide variety of species and an important resource for annotation of genomic sequence. The TIGR Gene Indices (http://www.tigr.org/tdb/tgi.shtml) are a collection of species-specific databases that use a highly refined protocol to analyze EST sequences in an attempt to identify the genes represented by that data and to provide additional information regarding those genes. Gene Indices are constructed by first clustering, then assembling EST and annotated gene sequences from GenBank for the targeted species. This process produces a set of unique, high-fidelity virtual transcripts, or Tentative Consensus (TC) sequences. The TC sequences can be used to provide putative genes with functional annotation, to link the transcripts to mapping and genomic sequence data, to provide links between orthologous and paralogous genes and as a resource for comparative sequence analysis.

Nucleotide sequences provide evidence of genetic exchange among distantly related lineages of Trypanosoma cruzi

Simple phylogenetic tests were applied to a large data set of nucleotide sequences from two nuclear genes and a region of the mitochondrial genome of Trypanosoma cruzi, the agent of Chagas' disease. Incongruent gene genealogies manifest genetic exchange among distantly related lineages of T. cruzi. Two widely distributed isoenzyme types of T. cruzi are hybrids, their genetic composition being the likely result of genetic exchange between two distantly related lineages. The data show that the reference strain for the T. cruzi genome project (CL Brener) is a hybrid. Well-supported gene genealogies show that mitochondrial and nuclear gene sequences from T. cruzi cluster, respectively, in three or four distinct clades that do not fully correspond to the two previously defined major lineages of T. cruzi. There is clear genetic differentiation among the major groups of sequences, but genetic diversity within each major group is low. We estimate that the major extant lineages of T. cruzi have diverged during the Miocene or early Pliocene (3–16 million years ago).

The evolution of the vertebrate Dlx gene family.

The vertebrate Dlx gene family consists of homeobox-containing transcription factors distributed in pairs on the same chromosomes as the Hox genes. To investigate the evolutionary history of Dlx genes, we have cloned five new zebrafish family members and have provided additional sequence information for two mouse genes. Phylogenetic analyses of Dlx gene sequences considered in the context of their chromosomal arrangements suggest that an initial tandem duplication produced a linked pair of Dlx genes after the divergence of chordates and arthropods but prior to the divergence of tunicates and vertebrates. This pair of Dlx genes was then duplicated in the chromosomal events that led to the four clusters of Hox genes characteristic of bony fish and tetrapods. It is possible that a pair of Dlx genes linked to the Hoxc cluster has been lost from mammals. We were unable to distinguish between independent duplication and retention of the ancestral state of bony vertebrates to explain the presence of a greater number of Dlx genes in zebrafish than mammals. Determination of the linkage relationship of these additional zebrafish Dlx genes to Hox clusters should help resolve this issue.

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.

The genomic and physical organization of Ty1-copia-like sequences as a component of large genomes in Pinus elliottii var. elliottii and other gymnosperms.

A DNA sequence, TPE1, representing the internal domain of a Ty1-copia retroelement, was isolated from genomic DNA of Pinus elliottii Engelm. var. elliottii (slash pine). Genomic Southern analysis showed that this sequence, carrying partial reverse transcriptase and integrase gene sequences, is highly amplified within the genome of slash pine and part of a dispersed element >4.8 kbp. Fluorescent in situ hybridization to metaphase chromosomes shows that the element is relatively uniformly dispersed over all 12 chromosome pairs and is highly abundant in the genome. It is largely excluded from centromeric regions and intercalary chromosomal sites representing the 18S-5.8S-25S rRNA genes. Southern hybridization with specific DNA probes for the reverse transcriptase gene shows that TPE1 represents a large subgroup of heterogeneous Ty1-copia retrotransposons in Pinus species. Because no TPE1 transcription could be detected, it is most likely an inactive element--at least in needle tissue. Further evidence for inactivity was found in recombinant reverse transcriptase and integrase sequences. The distribution of TPE1 within different gymnosperms that contain Ty1-copia group retrotransposons, as shown by a PCR assay, was investigated by Southern hybridization. The TPE1 family is highly amplified and conserved in all Pinus species analyzed, showing a similar genomic organization in the three- and five-needle pine species investigated. It is also present in spruce, bald cypress (swamp cypress), and in gingko but in fewer copies and a different genomic organization.

An induction gene trap screen in embryonic stem cells: Identification of genes that respond to retinoic acid in vitro.

We have developed a novel induction gene trap approach that preselects in vitro for integrations into genes that lie downstream of receptor/ligand-mediated signaling pathways. Using this approach, we have identified 20 gene trap integrations in embryonic stem cells, 9 of which were induced and 11 of which were repressed after exposure to exogenous retinoic acid (RA). All but one of these integrations showed unique spatially restricted or tissue-specific patterns of expression between 8.5 and 11.5 days of embryogenesis. Interestingly, expression was observed in tissues that are affected by alterations in RA levels during embryogenesis. Sequence analysis of fusion transcripts from six integrations revealed five novel gene sequences and the previously identified protooncogene c-fyn. To date, germ-line transmission and breeding has uncovered one homozygous embryonic lethal and three homozygous viable insertions. These studies demonstrate the potential of this induction gene trap approach for identifying and mutating genes downstream of signal transduction pathways.

A nuclear gene of eubacterial origin in Euglena gracilis reflects cryptic endosymbioses during protist evolution.

Genes for glycolytic and Calvin-cycle glyceraldehyde-3-phosphate dehydrogenase (GAPDH) of higher eukaryotes derive from ancient gene duplications which occurred in eubacterial genomes; both were transferred to the nucleus during the course of endosymbiosis. We have cloned cDNAs encoding chloroplast and cytosolic GAPDH from the early-branching photosynthetic protist Euglena gracilis and have determined the structure of its nuclear gene for cytosolic GAPDH. The gene contains four introns which possess unusual secondary structures, do not obey the GT-AG rule, and are flanked by 2- to 3-bp direct repeats. A gene phylogeny for these sequences in the context of eubacterial homologues indicates that euglenozoa, like higher eukaryotes, have obtained their GAPDH genes from eubacteria via endosymbiotic (organelle-to-nucleus) gene transfer. The data further suggest that the early-branching protists Giardia lamblia and Entamoeba histolytica--which lack mitochondria--and portions of the trypanosome lineage have acquired GAPDH genes from eubacterial donors which did not ultimately give rise to contemporary membrane-bound organelles. Evidence that "cryptic" (possibly ephemeral) endosymbioses during evolution may have entailed successful gene transfer is preserved in protist nuclear gene sequences.

Generation of a high-titer retroviral vector capable of expressing high levels of the human beta-globin gene.

Retrovirus-mediated gene transfer into hematopoietic cells may provide a means of treating both inherited and acquired diseases involving hematopoietic cells. Implementation of this approach for disorders resulting from mutations affecting the beta-globin gene (e.g., beta-thalassemia and sickle cell anemia), however, has been hampered by the inability to generate recombinant viruses able to efficiently and faithfully transmit the necessary sequences for appropriate gene expression. We have addressed this problem by carefully examining the interactions between retroviral and beta-globin gene sequences which affect vector transmission, stability, and expression. First, we examined the transmission properties of a large number of different recombinant proviral genomes which vary both in the precise nature of vector, beta-globin structural gene, and locus control region (LCR) core sequences incorporated and in the placement and orientation of those sequences. Through this analysis, we identified one specific vector, termed M beta 6L, which carries both the human beta-globin gene and core elements HS2, HS3, and HS4 from the LCR and faithfully transmits recombinant proviral sequences to cells with titers greater than 10(6) per ml. Populations of murine erythroleukemia (MEL) cells transduced by this virus expressed levels of human beta-globin transcript which, on a per gene copy basis, were 78% of the levels detected in an MEL-derived cell line, Hu11, which carries human chromosome 11, the site of the beta-globin locus. Analysis of individual transduced MEL cell clones, however, indicated that, while expression was detected in every clone tested (n = 17), the levels of human beta-globin treatment varied between 4% and 146% of the levels in Hu11. This clonal variation in expression levels suggests that small beta-globin LCR sequences may not provide for as strict chromosomal position-independent expression of beta-globin as previously suspected, at least in the context of retrovirus-mediated gene transfer.

The active digestion of uniparental chloroplast DNA in a single zygote of Chlamydomonas reinhardtii is revealed by using the optical tweezer

The non-Mendelian inheritance of organelle genes is a phenomenon common to almost all eukaryotes, and in the isogamous alga Chlamydomonas reinhardtii, chloroplast (cp) genes are transmitted from the mating type positive (mt+) parent. In this study, the preferential disappearance of the fluorescent cp nucleoids of the mating type negative (mt−) parent was observed in living young zygotes. To study the change in cpDNA molecules during the preferential disappearance, the cpDNA of mt+ or mt− origin was labeled separately with bacterial aadA gene sequences. Then, a single zygote with or without cp nucleoids was isolated under direct observation by using optical tweezers and investigated by nested PCR analysis of the aadA sequences. This demonstrated that cpDNA molecules are digested completely during the preferential disappearance of mt− cp nucleoids within 10 min, whereas mt+ cpDNA and mitochondrial DNA are protected from the digestion. These results indicate that the non-Mendelian transmission pattern of organelle genes is determined immediately after zygote formation.

Evolutionary analyses of hedgehog and Hoxd-10 genes in fish species closely related to the zebrafish

The study of development has relied primarily on the isolation of mutations in genes with specific functions in development and on the comparison of their expression patterns in normal and mutant phenotypes. Comparative evolutionary analyses can complement these approaches. Phylogenetic analyses of Sonic hedgehog (Shh) and Hoxd-10 genes from 18 cyprinid fish species closely related to the zebrafish provide novel insights into the functional constraints acting on Shh. Our results confirm and extend those gained from expression and crystalline structure analyses of this gene. Unexpectedly, exon 1 of Shh is found to be almost invariant even in third codon positions among these morphologically divergent species suggesting that this exon encodes for a functionally important domain of the hedgehog protein. This is surprising because the main functional domain of Shh had been thought to be that encoded by exon 2. Comparisons of Shh and Hoxd-10 gene sequences and of resulting gene trees document higher evolutionary constraints on the former than on the latter. This might be indicative of more general evolutionary patterns in networks of developmental regulatory genes interacting in a hierarchical fashion. The presence of four members of the hedgehog gene family in cyprinid fishes was documented and their homologies to known hedgehog genes in other vertebrates were established.

Evaluating multiple alternative hypotheses for the origin of Bilateria: An analysis of 18S rRNA molecular evidence

Six alternative hypotheses for the phylogenetic origin of Bilateria are evaluated by using complete 18S rRNA gene sequences for 52 taxa. These data suggest that there is little support for three of these hypotheses. Bilateria is not likely to be the sister group of Radiata or Ctenophora, nor is it likely that Bilateria gave rise to Cnidaria or Ctenophora. Instead, these data reveal a close relationship between bilaterians, placozoans, and cnidarians. From this, several inferences can be drawn. Morphological features that previously have been identified as synapomorphies of Bilateria and Ctenophora, e.g., mesoderm, more likely evolved independently in each clade. The endomesodermal muscles of bilaterians may be homologous to the endodermal muscles of cnidarians, implying that the original bilaterian mesodermal muscles were myoepithelial. Placozoans should have a gastrulation stage during development. Of the three hypotheses that cannot be falsified with the 18S rRNA data, one is most strongly supported. This hypothesis states that Bilateria and Placozoa share a more recent common ancestor than either does to Cnidaria. If true, the simplicity of placozoan body architecture is secondarily derived from a more complex ancestor. This simplification may have occurred in association with a planula-type larva becoming reproductive before metamorphosis. If this simplification took place during the common history that placozoans share with bilaterians, then placozoan genes that contain a homeobox, such as Trox2, should be explored, for they may include the gene or genes most closely related to Hox genes of bilaterians.

rrndb: the Ribosomal RNA Operon Copy Number Database

The Ribosomal RNA Operon Copy Number Database (rrndb) is an Internet-accessible database containing annotated information on rRNA operon copy number among prokaryotes. Gene redundancy is uncommon in prokaryotic genomes, yet the rRNA genes can vary from one to as many as 15 copies. Despite the widespread use of 16S rRNA gene sequences for identification of prokaryotes, information on the number and sequence of individual rRNA genes in a genome is not readily accessible. In an attempt to understand the evolutionary implications of rRNA operon redundancy, we have created a phylogenetically arranged report on rRNA gene copy number for a diverse collection of prokaryotic microorganisms. Each entry (organism) in the rrndb contains detailed information linked directly to external websites including the Ribosomal Database Project, GenBank, PubMed and several culture collections. Data contained in the rrndb will be valuable to researchers investigating microbial ecology and evolution using 16S rRNA gene sequences. The rrndb web site is directly accessible on the WWW at http://rrndb.cme.msu.edu.