A new technique for selective identification and mapping of enhancers within long genomic sequences

We report a new experimental method of direct selection, identification, and mapping of potential enhancer sequences within extended stretches of genomic DNA. The method allows simultaneous cloning of a quantity of sequences instead of tedious screening of the separate ones, thus providing a robust and high-throughput approach to the mapping of enhancers. The selection procedure is based on the ability of such sequences to activate a minimal promoter that drives expression of a selective gene. To this end a mixture of short DNA fragments derived from the segment of interest was cloned in a retroviral vector containing the neomycin phosphotransferase II gene under control of a cytomegalovirus (CMV) minimal promoter. The pool of retroviruses obtained was used to infect HeLa cells and then to select neomycin-resistant colonies containing constructs with enhancer-like sequences. The pool of the genomic fragments was rescued by PCR and cloned, forming a library of the potential enhancers. Fifteen enhancer-like fragments were selected from 1-Mb human genome locus, and enhancer activity of 13 of them was verified in a transient transfection reporter gene assay. The sequences selected were found to be predominantly located near 5' regions of genes or within gene introns.     

Partial shotgun sequencing of the Boechera stricta genome reveals extensive microsynteny and promoter conservation with Arabidopsis

Comparative genomics provides insight into the evolutionary dynamics that shape discrete sequences as well as whole genomes. To advance comparative genomics within the Brassicaceae, we have end sequenced 23,136 medium-sized insert clones from Boechera stricta, a wild relative of Arabidopsis (Arabidopsis thaliana). A significant proportion of these sequences, 18,797, are nonredundant and display highly significant similarity (BLASTn e-value < or = 10(-30)) to low copy number Arabidopsis genomic regions, including more than 9,000 annotated coding sequences. We have used this dataset to identify orthologous gene pairs in the two species and to perform a global comparison of DNA regions 5' to annotated coding regions. On average, the 500 nucleotides upstream to coding sequences display 71.4% identity between the two species. In a similar analysis, 61.4% identity was observed between 5' noncoding sequences of Brassica oleracea and Arabidopsis, indicating that regulatory regions are not as diverged among these lineages as previously anticipated. By mapping the B. stricta end sequences onto the Arabidopsis genome, we have identified nearly 2,000 conserved blocks of microsynteny (bracketing 26% of the Arabidopsis genome). A comparison of fully sequenced B. stricta inserts to their homologous Arabidopsis genomic regions indicates that indel polymorphisms >5 kb contribute substantially to the genome size difference observed between the two species. Further, we demonstrate that microsynteny inferred from end-sequence data can be applied to the rapid identification and cloning of genomic regions of interest from nonmodel species. These results suggest that among diploid relatives of Arabidopsis, small- to medium-scale shotgun sequencing approaches can provide rapid and cost-effective benefits to evolutionary and/or functional comparative genomic frameworks.     

DNA sequence of an immediate-early gene (IEmRNA-5) of herpes simplex virus type I

We describe a 2560 base pair herpes simplex virus type 1 (HSV-1) DNA sequence containing the entire immediate-early mRNA-5 (IEmRNA-5) gene. The 3' and 5' termini of IEmRNA-5 were mapped within this DNA sequence by single-strand specific endonuclease protection experiments. The IEmRNA-5 gene contains DNA sequences from both the unique (Us) and reiterated (TRs/IRs) regions of the HSV-1 DNA short component and is interrupted by a single intron mapping in TRs/IRs. A search of the transcribed DNA sequence revealed no initiator codon within TRs/IRs. The first ATG was located 6 bases into Us sequences and this reading frame (316 codons) was also observed in the 3' transcribed region. The oligonucleotide sequences adjacent to the IEmRNA-5 termini are discussed in relation to those of the HSV-1 thymidine kinase gene and other genes transcribed by RNA polymerase II.     

Linkage arrangement of human placental lactogen and growth hormone genes

Human placental lactogen (hPL) and growth hormone (hGH) are two hormones thought to have evolved from a common ancestral gene (along with prolactin), yet they have quite different functions and specificities. The nucleic acid sequences of the respective cDNAs of the two genes share considerable homology, as well as the existence of multiple forms of each gene within the genome. In this study we report on the linkage arrangement of several genes from this group. Two hPL-like genes as well as an hGH gene are shown to be linked within a 38-kilobase pair region of DNA. Linkage between a variant hGH gene and an hPL gene is also shown. The orientation and structural organization of these genes was previously established using 5'- and 3'-specific probes from a placental lactogen cDNA clone and detailed restriction endonuclease mapping. Restriction fragments from the overlapping clones were verified by comparison to digests of high molecular weight genomic DNA. In addition, the location of a specific class of repetitive DNA sequences, the Alu family, was mapped on these clones using the recombinant clone BLUR 8. All members of this multigene family have Alu repeat sequences either immediately flanking their 3' or 5' untranslated regions or within their intervening sequences.     

Compact intergenic regions of the pufferfish genome facilitate isolation of gene promoters: characterization of Fugu 3'-phosphoadenosine 5'-phosphosulfate synthase 2 (fPapss2) gene promoter function in transgenic Xenopus

The highly compact nature of the pufferfish (Fugu rubripes) genome renders it a useful tool not only for annotating coding regions within vertebrate genomes, but also for the identification of sequences important to gene regulation. Indeed, owing to this compaction it will be feasible in many instances to initiate analyses using entire intergenic regions when mapping gene promoters; a strategy that is very rarely feasible with the expanded genomes of other species. Stemming from our interest in studying promoters expressed in chondrocytes, we selected for study the intergenic region upstream of Fugu 3'-phosphoadenosine 5'-phosphosulfate synthase 2, fPapss2, a gene required for the normal development of cartilage extracellular matrix. Functional characterization of the entire fPapss2 5' intergenic region was carried out by monitoring expression of the enhanced green fluorescent protein (EGFP) gene reporter in the developing cartilage of transgenic Xenopus laevis. By evaluating a series of 5' intergenic region deletions we defined a minimal fPapss2 sequence of approximately 300 bp that was essential for EGFP expression in tadpole cartilage. This functional analysis of an entire Fugu intergenic region, combined with the efficiency of Xenopus transgenesis, serves as a model for the rapid characterization of evolutionarily-conserved regulatory regions of other pufferfish genes.     

Paenibacillus larvae 16S-23S rDNA intergenic transcribed spacer (ITS) regions: DNA fingerprinting and characterization

Paenibacillus larvae is the causative agent of American foulbrood in honey bee (Apis mellifera) larvae. PCR amplification of the 16S-23S ribosomal DNA (rDNA) intergenic transcribed spacer (ITS) regions, and agarose gel electrophoresis of the amplified DNA, was performed using genomic DNA collected from 134 P. larvae strains isolated in Connecticut, six Northern Regional Research Laboratory stock strains, four strains isolated in Argentina, and one strain isolated in Chile. Following electrophoresis of amplified DNA, all isolates exhibited a common migratory profile (i.e., ITS-PCR fingerprint pattern) of six DNA bands. This profile represented a unique ITS-PCR DNA fingerprint that was useful as a fast, simple, and accurate procedure for identification of P. larvae. Digestion of ITS-PCR amplified DNA, using mung bean nuclease prior to electrophoresis, characterized only three of the six electrophoresis bands as homoduplex DNA and indicating three true ITS regions. These three ITS regions, DNA migratory band sizes of 915, 1010, and 1474 bp, signify a minimum of three types of rrn operons within P. larvae. DNA sequence analysis of ITS region DNA, using P. larvae NRRL B-3553, identified the 3' terminal nucleotides of the 16S rRNA gene, 5' terminal nucleotides of the 23S rRNA gene, and the complete DNA sequences of the 5S rRNA, tRNA(ala), and tRNA(ile) genes. Gene organization within the three rrn operon types was 16S-23S, 16S-tRNA(ala)-23S, and l6S-5S-tRNA(ile)-tRNA(ala)-23S and these operons were named rrnA, rrnF, and rrnG, respectively. The 23S rRNA gene was shown by I-CeuI digestion and pulsed-field gel electrophoresis of genomic DNA to be present as seven copies. This was suggestive of seven rrn operon copies within the P. larvae genome. Investigation of the 16S-23S rDNA regions of this bacterium has aided the development of a diagnostic procedure and has helped genomic mapping investigations via characterization of the ITS regions.     

Isolation and characterization of six different chicken actin genes.

Genes representing six different actin isoforms were isolated from a chicken genomic library. Cloned actin cDNAs as well as tissue-specific mRNAs enriched in different actin species were used as hybridization probes to group individual actin genomic clones by their relative thermal stability. Restriction maps showed that these actin genes were derived from separate and nonoverlapping regions of genomic DNA. Of the six isolated genes, five included sequences from both the 5' and 3' ends of the actin-coding area. Amino acid sequence analysis from both the NH2- and COOH-terminal regions provided for the unequivocal identification of these genes. The striated isoforms were represented by the isolated alpha-skeletal, alpha-cardiac, and alpha-smooth muscle actin genes. The nonmuscle isoforms included the beta-cytoplasmic actin gene and an actin gene fragment which lacked the 5' coding and flanking sequence; presumably, this region of DNA was removed from this gene during construction of the genomic library. Unexpectedly, a third nonmuscle chicken actin gene was found which resembled the amphibian type 5 actin isoform (J. Vandekerckhove, W. W. Franke, and K. Weber, J. Mol. Biol., 152:413-426). This nonmuscle actin type has not been previously detected in warm-blooded vertebrates. We showed that interspersed, repeated DNA sequences closely flanked the alpha-skeletal, alpha-cardiac, beta-, and type 5-like actin genes. The repeated DNA sequences which surround the alpha-skeletal actin-coding regions were not related to repetitious DNA located on the other actin genes. Analysis of genomic DNA blots showed that the chicken actin multigene family was represented by 8 to 10 separate coding loci. The six isolated actin genes corresponded to 7 of 11 genomic EcoRI fragments. Only the alpha-smooth muscle actin gene was shown to be split by an EcoRI site. Thus, in the chicken genome each actin isoform appeared to be encoded by a single gene.     

Genomic organization and the 5' upstream sequences associated with the specific spatio-temporal expression of HrEpiC, an epidermis-specific gene of the ascidian Halocynthia roretzi.

An epidermis-specific gene HrEpiC of the ascidian Halocynthia roretzi is activated in all presumptive blastomeres by the 64-cell stage. To explore the molecular mechanisms underlying the regulation of the timing of activation of HrEpiC, we studied the genomic organization and the 5' upstream sequences of HrEpiC associated with specific spatio-temporal expression of the gene. The restriction site mapping and sequencing of genomic clones showed that the H. roretzi genome contained two copies of HrEpiC gene, HrEpiC1 and HrEpiC2, aligned tandemly in about 8 kb of the genome. Analysis of various deletion constructs with the 5' flanking sequences of HrEpiC1 revealed that 103 bp of the 5' flanking region was sufficient for the minimal epidermis-specific expression of HrEpiC1 and that the region between -281 bp and -198 bp of the 5' flanking region was associated with the amplification of the minimal expression of the reporter gene in the epidermis. This module between -281 bp and -198 bp was also shown to be associated with the timing of the activation of HrEpiC1 by the 64-cell stage. We discussed how the spatio-temporal expression pattern of HrEpiC1 is regulated by the two modules.     

Structure and function of the non-coding regions of hepatitis C viral RNA

At the 5' and 3' end of genomic HCV RNA there are two highly conserved, untranslated regions, 5'UTR and 3'UTR. These regions are organized into spatially ordered structures and they play key functions in regulation of processes of the viral life cycle. Most nucleotides of the region located at the 5' side of the coding sequence serve as an internal ribosomal entry site, IRES, which directs cap-independent translation. The RNA fragment present at the 3' end of the genome is required for virus replication and probably contributes to translation of viral proteins. During virus replication its genomic strand is transcribed into a strand of minus polarity, the replicative strand. Its 3' terminus is responsible for initiation of synthesis of descendant genomic strands. This article summarizes our current knowledge on the structure and function of the non-coding regions of hepatitis C genomic RNA, 5'UTR and 3'UTR, and the complementary sequences of the replicative viral strand.     

Cloned embryonic DNA sequences flanking the mouse immunoglobulin C gamma 3 and C gamma 1 genes

To investigate the DNA surrounding genes for immunoglobulin heavy chain constant (CH) regions, we have isolated two clones bearing a C gamma 3 gene and two bearing a C gamma 1 gene from a library of mouse embryo DNA fragments. The C gamma 3 clones span 8.6 kilobase pairs (kb) on the 5' side of the gene and 6.7 kb on its 3' side, while the C gamma 1 clones together span 13 kb of 5' flanking sequence and 2.5 kb of 3' flanking sequence. Restriction mapping of the C gamma 3 gene indicates that intervening sequences divide the gene into segments of domain size, as in other CH genes. Hybridization of clone fragments to restriction digests of mouse DNA indicates that both the C gamma 1 and C gamma 3 genes probably occur as single copies in the genome. Moreover, the entire cloned sequences on the 5' side of both genes appear to be unique in the genome, indicating that no large common sequences flank CH genes. Restriction data suggest that the C gamma 3 gene is 37-40 kb 5' to the C gamma 1 gene.     

In silico identification and characterization of a putative phosphatidylinositol 4-phosphate 5-kinase (PIP5K) gene in Eimeria tenella

Phosphatidylinositol 4-phosphate 5-kinases (PIP5Ks) play diverse roles in the cellular biology of many organisms, including signal transduction, secretion and vesicular trafficking, and regulation of cytoskeleton assembly. Discovery of the PIP5K gene in Eimeria tenella may shed light on its role in the biology of this avian protozoan, and afford further understanding of the cell-host interaction, particularly during the invasion process. In this study, we report the identification of the PIP5K coding region in the genome sequence of Eimeria tenella using in silico gene prediction approaches. Prediction of the PIP5K coding sequence was confirmed by mapping the full-length cDNA sequence, generated via the Rapid Amplification of cDNA Ends (RACE) method, to the genomic sequence. The putative PIP5K gene of Eimeria tenella is located on the complementary strand of the E1080B12.b1 contig, and comprises 12 exons. Further analysis showed that the coding region spans from exon 1 to exon 7, with all exons obeying the adopted 'gt...ag' splicing rule of intronic sequences. Consensus of the hexameric 5' donor-splice site was deduced as GTRDBB... and the consensus for the 3' acceptor-splice sites as ...BHDYAG. The gene encodes a 252-amino acid residue protein. Domain search and protein fold recognition analyses provide compelling evidences that the deduced protein is a PIP5K.     

Novel rearrangements of herpes simplex virus DNA sequences resulting from duplication of a sequence within the unique region of the L component.

We constructed insertion mutants of herpes simplex virus type 1 that contained a duplication of DNA sequences from the BamHI-L fragment (map units 0.706 to 0.744), which is located in the unique region of the L component (UL) of the herpes simplex virus type 1 genome. The second copy of the BamHI-L sequence was inserted in inverted orientation into the viral thymidine kinase gene (map units 0.30 to 0.32), also located within UL. A significant fraction of the progeny produced by these insertion mutants had genomes with rearranged DNA sequences, presumably resulting from intramolecular or intermolecular recombination between the BamHI-L sequences at the two different genomic locations. The rearranged genomes either had an inversion of the DNA sequence flanked by the duplication or were recombinant molecules in which different regions of the genome had been duplicated and deleted. Genomic rearrangements similar to those described here have been reported previously but only for herpes simplex virus insertion mutants containing an extra copy of the repetitive a sequence. Such rearrangements have not been reported for insertion mutants that contain duplications of herpes simplex virus DNA sequences from largely unique regions of the genome. The implications of these results are discussed.