Screening the V2R-type putative odorant receptor gene repertoire in bitterling Tanakia lanceolata

To study evolution of dinucleotide simple sequence repeats (diSSRs) we searched recently available mammalian genomes for UTR-localized diSSRs with conserved upstream flanking sequences (CFS). There were 252 reported Homo sapiens genes containing the repeats (AC)n, (GT)n, (AG)n or (CT)n in their UTRs including 22 (8.7%) with diSSR-upstream flanking sequences conserved comparing divergent mammalian lineages represented by Homo sapiens and the marsupial, Monodelphis domestica. Of these 22 genes, 19 had known functions including 18 (95%) that proved critical for mammalian nervous systems (Fishers exact test, P<0.0001). The remaining gene, Cd2ap, proved critical for development of kidney podocytes, cells that have multiple similarities to neurons. Gene functions included voltage and chloride channels, synapse-associated proteins, neurotransmitter receptors, axon and dendrite pathfinders, a NeuroD potentiator and other neuronal activities. Repeat length polymorphism was confirmed for 68% of CFS diSSRs even though these repeats were nestled among highly conserved sequences. This finding supports a hypothesis that SSR polymorphism has functional implications. A parallel study was performed on the self-complementary diSSRs (AT)n and (GC)n. When flanked by conserved sequences, the self-complementary diSSR (AT)n was also associated with genes expressed in the developing nervous system. Our findings implicate functional roles for diSSRs in nervous system development.     

Simple repeat evolution includes dramatic primary sequence changes that conserve folding potential

We previously demonstrated that many "weak-folding" simple repeats were replaced during evolution by alternative weak-folding repeats. This suggested repeat selection at the level of higher order structure potential. Here, we demonstrate similar phenomena for "strong-folding" simple repeats in non-coding DNA. The Rabgap1 gene's 3' UTR contained the self-complementary repeat (AT)n in Homo sapiens but, in Mus musculus, this site was occupied by the complementary repeats (GT)n and (AC)n. Similarly, primate Plag1 UTRs contained various (GT)n-(AC)n palindromes but in rodents, this site was occupied by (AT)n, preserving folding potential more than primary sequence. The Znf516, Senp1, Rock2, and other UTRs exhibited similar replacements. In the Bnc2 UTR, (AT)n was replaced by sequences that evolved with approximate symmetry about a central axis, a pattern difficult to explain without invoking selection to preserve secondary structure. These observations reflect a predictable evolutionary pattern for some common non-coding genomic sequences.     

Tandem repeats of a specific alternating purine-pyrimidine DNA sequence adjacent to protamine genes in the rainbow trout that can exist in the Z form

We have located an extensive (AC)n-rich but specific sequence downstream of three rainbow trout protamine genes. Although sharing considerable sequence homology, including a perfectly conserved 46 base pair repeat, the sequences exhibit a regular heterogeneity in the length of the (AC)n-rich tracts. Radioimmunoassay experiments, S1 nuclease sensitivity studies, two-dimensional electrophoretic analysis, and immunoelectron microscopy studies have been used to determine if the region could assume a Z DNA conformation. It was found that, in a supercoiled plasmid, the (AC)n-rich region has the ability to attain the Z DNA conformation under physiological conditions.     

Effects of sequence on repeat expansion during DNA replication

Small DNA repeat tracts are located throughout the human genome. The tracts are unstable, and expansions of certain repeat sequences cause neuromuscular disease. DNA expansions appear to be associated with lagging-strand DNA synthesis and DNA repair. At some sites of repeat expansion, e.g. the myotonic dystrophy type 2 (DM2) tetranucleotide repeat expansion site, more than one repeat tract with similar sequences lie side by side. Only one of the DM2 repeat tracts, however, is found to expand. Thus, DNA base sequence is a possible factor in repeat tract expansion. Here we determined the expansion potential, during DNA replication by human DNA polymerase β, of several tetranucleotide repeat tracts in which the repeat units varied by one or more bases. The results show that subtle changes, such as switching T for C in a tetranucleotide repeat, can have dramatic consequences on the ability of the nascent-strand repeat tract to expand during DNA replication. We also determined the relative stabilities of self-annealed 100mer repeats by melting-curve analysis. The relative stabilities did not correlate with the relative potentials of the analogous repeats for expansion during DNA replication, suggesting that hairpin formation is not required for expansion during DNA replication.     

Mung bean nuclease cleavage of a dA + dT-rich sequence or an inverted repeat sequence in supercoiled PM2 DNA depends on ionic environment

We have determined the nucleotide sequences around two alternative sites cleaved in supercoiled PM2 DNA by single-strand-specific mung bean nuclease in different ionic environments. In 10 mM Tris-HC1 (pH 7.0, 37 degrees C), the major site is a dA+dT-rich sequence which maps with a known early denaturation region at 0.75 map units. About 30 cleavages occurred in a 135 bp region. Cleavages were largely excluded at (dA)n . (dT)n (n = 3-7) sequences. Cleavage patterns of this type have not been previously observed in dA+dT-rich sequences. With the addition of 0.1 M NaC1 the major alternative site occurred in a hyphenated inverted repeat sequence 500 bp away (0.70 map units) and did not map to an early denaturation region. One major and 4 minor cleavages occurred in the region between the repeats, suggesting that a hairpin containing at most a 12 bp stem and 10 base loop is recognized. The basis for nuclease recognition of the dA+dT-rich sequence is not clear. The differences in the sequences and cleavage patterns at the alternative sites indicate that their secondary structures differ.     

The sequence specificity of homeodomain-DNA interaction

The Drosophila developmental gene, engrailed, encodes a sequence-specific DNA binding activity. Using deletion constructs expressed as fusion proteins in E. coli, we localized this activity to the conserved homeodomain (HD). The binding site consensus, TCAATTAAAT, is found in clusters in the engrailed regulatory region. Weak binding of the En HD to one copy of a synthetic consensus is enhanced by adjacent copies. The distantly related HD encoded by fushi tarazu binds to the same sites as the En HD, but differs in its preference for related sites. Both HDs bind a second type of sequence, a repeat of TAA. The similarity in sequence specificity of En and Ftz HDs suggests that, within families of DNA binding proteins, close relatives will exhibit similar specificities. Competition among related regulatory proteins might govern which protein occupies a given binding site and consequently determine the ultimate effect of cis-acting regulatory sites.     

Inference of splicing regulatory activities by sequence neighborhood analysis

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.     

Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

A novel transposon-like repeat interrupted by an LTR element occurs in a cluster of B1 repeats in the mouse c-mos locus.

The mouse genomic locus containing the oncogene c-mos was analyzed for repetitive DNA sequences. We found a single B1 repeat 10 kb upstream and three B1 repeats 0.6 kb, 2.7 kb, and 5.4 kb, respectively, downstream from c-mos. The B1 repeat closest to c-mos contains an internal 7-bp duplication and a 18-bp insertion. Localized between the last two B1 repeats is a copy of a novel mouse repeat. Sequence comparison of three copies of this novel repeat family shows that they a) contain a conserved BglII site, b) are approximately 420 bp long, c) possess internal 50-bp polypurine tracts, and d) have structural characteristics of transposable elements. They are present in about 1500 copies per haploid genome in the mouse, but are not detectable in DNA of other mammals. The BglII repeat downstream from c-mos is interrupted by a single 632-bp LTR element. We estimate that approximately 1200 copies of this element are present per haploid genome in BALB/c mice. It shares sequence homology in the R-U5 region with an LTR element found in 129/J mice.     

Identification and functional validation of a 5' upstream regulatory sequence in the human tyrosinase gene homologous to the locus control region of the mouse tyrosinase gene

Comparison analysis of the sequences of the mouse and human genomes has proven a powerful approach in identifying functional regulatory elements within the non-coding regions that are conserved through evolution between homologous mammalian loci. Here, we applied computational analysis to identify regions of homology in the 5' upstream sequences of the human tyrosinase gene, similar to the locus control region (LCR) of the mouse tyrosinase gene, located at -15 kb. We detected several stretches of homology within the first 30 kb 5' tyrosinase gene upstream sequences of both species that include the proximal promoter sequences, the genomic region surrounding the mouse LCR, and further upstream segments. We cloned and sequenced a 5' upstream regulatory sequence found between -8 and -10 kb of the human tyrosinase locus (termed h5'URS) homologous to the mouse LCR sequences, and confirmed the presence of putative binding sites at -9 kb, homologous to those described in the mouse tyrosinase LCR core. Finally, we functionally validated the presence of a tissue-specific enhancer in the h5'URS by transient transfection analysis in human and mouse cells, as compared with homologous DNA sequences from the mouse tyrosinase locus. Future experiments in cells and transgenic animals will help us to understand the in vivo relevance of this newly described h5'URS sequence as a potentially important regulatory element for the correct expression of the human tyrosinase gene.     

Compilation and analysis of Bacillus subtilis sigma A-dependent promoter sequences: evidence for extended contact between RNA polymerase and upstream promoter DNA.

Sequence analysis of 236 promoters recognized by the Bacillus subtilis sigma A-RNA polymerase reveals an extended promoter structure. The most highly conserved bases include the -35 and -10 hexanucleotide core elements and a TG dinucleotide at position -15, -14. In addition, several weakly conserved A and T residues are present upstream of the -35 region. Analysis of dinucleotide composition reveals A2- and T2-rich sequences in the upstream promoter region (-36 to -70) which are phased with the DNA helix: An tracts are common near -43, -54 and -65; Tn tracts predominate at the intervening positions. When compared with larger regions of the genome, upstream promoter regions have an excess of An and Tn sequences for n > 4. These data indicate that an RNA polymerase binding site affects DNA sequence as far upstream as -70. This sequence conservation is discussed in light of recent evidence that the alpha subunits of the polymerase core bind DNA and that the promoter may wrap around RNA polymerase.     

Evolutionarily conserved and divergent regulatory sequences in the fish rod opsin promoter

Fish have multiple types and subtypes of opsin genes that are expressed in a highly regulated manner in retinal photoreceptor cells. In the rod opsin proximal promoter region (RPPR) of zebrafish (Danio rerio), the BAT 1 regulatory region contains highly conserved OTX (GATTA) and OTX-like (TATTA) sequences that can be recognized by the mammalian cone–rod homeobox (CRX) protein. However, binding of zebrafish crx to the OTX sequence has remained elusive. In contrast to the BAT 1 region, the Ret 1 region, located approximately 20 bp upstream of the BAT 1 region in mammals, is not conserved in zebrafish. In the Ret 1 region, even the core OTX-like sequence (AATTA sequence in mammals) is destructed. We show in this study that a region between Ret 1 and BAT 1 (denoted IRB, Inter-Ret 1-BAT 1) is highly conserved among fish species. Using electrophoretic mobility shift assay (EMSA), we show that zebrafish crx binds to the conserved OTX sequence and that the fish-specific IRB region specifically binds elements present in both retinal and brain nuclear extracts of zebrafish. These results imply that the regulatory mechanisms of opsin gene expression consist not only of evolutionarily conserved but also of divergent machinery among different animal taxa.