Subfamily relationships and clustering of rabbit C repeats

C repeats constitute the predominant family of short interspersed repeats (SINEs) in the rabbit genome. Determination of the nucleotide sequence 5' to rabbit zeta-globin genes reveals clusters of C repeats, and analysis of these and other sequenced regions of rabbit chromosomes shows that the C repeats have a strong tendency to insert within or in close proximity to other C repeats. An alignment of 44 members of the C repeat family shows that they are composites of different sequences, including a tRNA-like sequence, a conserved central core, a stretch of repeating CT dinucleotides, and an A-rich tract. Cladograms generated by both parsimony and cluster analysis subdivide the C repeats into at least three distinct subfamilies. Nucleotides at sites diagnostic for subfamilies appear to have changed in a punctuated and progressive manner during evolution, indicating that a limited number of progenitors have given rise to new repeats in waves of dispersion. C repeats that insert into preexisting C repeats belong to subfamilies that are proposed to have been propagated more recently; hence, these data support the model of dispersion in successive waves. The divergence among the oldest group of C repeats is greater than that observed for the analogous Alu repeats in humans, indicating that rabbit C repeats have been propagating longer than human Alu repeats. The improved consensus sequence for these repeats is similar to that of the predominant artiodactyl SINE in both the tRNA-like region and a central region. Because members of different subfamilies cross-hybridize very poorly, hybridization data with representatives of each subfamily provide a new minimal estimate, 234,000, for the copy number of C repeats in the rabbit haploid genome, although it is likely that the actual value is closer to 1 million.     

A Database of Experimental Results on Globin Gene Expression

Information on gene expression and regulation is expanding rapidly, and the complexity of the experimental design and data makes unique demands on databases to store the results. We describe a prototype database containing experimental results on the expression of mammalian β-like globin genes, along with several query methods for accessing the information. The database has tables for DNA transfer experiments, protein–DNA binding results, and positions of DNase hypersensitive sites, which make extensive use of nested data structures. Comparison of data from various mammals is accomplished by providing a common coordinate system via a simultaneous alignment of matching DNA sequences. Interactive access to the database is available at a site called theGlobin Gene Serveron the World Wide Web (http://globin.cse. psu.edu). This software should be useful for any genetic system in which DNA sequence data are available.     

Comparative analysis of the locus control region of the rabbit beta-like gene cluster: HS3 increases transient expression of an embryonic epsilon-globin gene.

The rabbit homolog to the locus control region (LCR) of the human beta-like globin gene cluster was isolated, and long segments containing the DNase I hypersensitive sites (HS) were sequenced. The order and spacing of HS4, HS3, HS2 and HS1 are conserved between rabbit and human. Alignment of these sequences with their homologs from human, goat, and mouse shows that very long segments of DNA match between species, for over a thousand base pairs on either side of the previously identified functional cores, indicating that some important functions are found outside the cores. The activity of rabbit HS2 and HS3 was tested by attaching each to a novel reporter gene constructed by inserting the luciferase coding region into the rabbit epsilon-globin gene. In contrast to previous reports showing no effect of human or mouse HS3 on transient expression, both the rabbit HS2 and HS3 DNA fragments separately increased transient expression from the epsilon-luciferase hybrid gene and expression from stably integrated constructs in K562 erythroleukemia cells.     

Nonlinear relationships among evolutionary rates identify regions of functional divergence in heat-shock protein 70 genes

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP- dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily C (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.      

Locating well-conserved regions within a pairwise alignment

Within a single alignment of two DNA sequences or two proteinsequences, some regions may be much better conserved than others.Such strong conservation may reveal a region that possessesan important function. When alignments are so long that it isinfeasible, or at least undesirable, to inspect them in completedetail, it is helpful to have an automatic process that computesinformation about the varying degree of conservation along thealignment and displays the information in a graphical representationthat is readily assimilated. This paper presents methods forcomputing several such ‘robustness measures’ ateach position of a given alignment. These methods are all veryspace-efficient; they use only space proportional to the sumof the two sequence lengths. To illustrate their effectiveness,one of the methods is used to locate particularly well-conservedregions in the ß-globin gene locus control regionand in the 5' flank of the -globin gene.     

Complete nucleotide sequence of the rabbit beta-like globin gene cluster. Analysis of intergenic sequences and comparison with the human beta-like globin gene cluster

The nucleotide sequence of the entire beta-like globin gene cluster of rabbits has been determined. This sequence of a continuous stretch of 44.5 x 10(3) base-pairs (bp) starts about 6 x 10(3) bp upstream from epsilon (the 5'-most gene) and ends about 12 x 10(3) bp downstream from beta (the 3'-most gene). Analysis of the sequence reveals that: (1) the sequence is relatively A + T rich (about 60%); (2) regions with high G + C content are associated with OcC repeats, a short interspersed repeated DNA in rabbits; (3) the distribution of polypurines, polypyrimidines and alternating purine/pyrimidine tracts is not random within the cluster; (4) most open reading frames are associated with known globin coding regions, OcC repeats or long interspersed repeats (L1 repeats); (5) the most prominent open reading frames are found in the L1 repeats; (6) different strand asymmetries in base composition are associated with embyronic and adult genes as well as the tandem L1 repeats at the 3' end of the cluster; and (7) essentially all the repeats appear to have been inserted by a transposon mechanism. A comparison of the sequence with itself by a dot-plot analysis has revealed nine new members of the OcC family of repeats in addition to the six previously reported. The OcC repeats tend to be clustered, particularly in the epsilon-gamma and gamma-psi delta intergenic regions. Dot-plot comparisons between the rabbit and the human clusters have revealed extensive sequence matches. Homology starts about 6 x 10(3) bp 5' to epsilon or as far upstream as the rabbit sequence is available. It continues throughout the entire cluster and stops about 0.7 x 10(3) bp 3' to beta, at which point several repeats have inserted in both rabbits and humans. Throughout the gene cluster, the homology is interrupted mainly by insertions or deletions in either the rabbit or the human genome. Almost all of the insertions are of known short or long repeated DNAs. The positions of the insertions are different in the two gene clusters, which indicates that both short and long repeats have been transposing throughout the genome for the time since the mammalian radiation. An alignment of rabbit and human sequences allows the calculation of the substitution rate around epsilon. Sequences far removed from the gene are evolving at a rate equivalent to the pseudogene rate, although some short regions show an apparently higher rate.(ABSTRACT TRUNCATED AT 400 WORDS)     

Directional and balancing selection in human beta-defensins

Background  

In primates, infection is an important force driving gene evolution, and this is reflected in the importance of infectious disease in human morbidity today. The beta-defensins are key components of the innate immune system, with antimicrobial and cell signalling roles, but also reproductive functions. Here we examine evolution of beta-defensins in catarrhine primates and variation within different human populations.