Evolution of chromosome bands: molecular ecology of noncoding DNA

Karyotypic comparison among Cebuella pygmaea, Callithrix jacchus and C. emiliae (Callitrichidae, Primates) and its taxonomic implications

The karyotype of Cebuella pygmaea (2n = 44) obtained by G-, C-banding, and NOR-staining is described. This species presents a heteromorphic C band in the intersticial region of the short arm of chromosome 2. The data obtained were compared with those previously described for the karyotypes of Callithrix jacchus and Callithrix emiliae. The three species differ in the amount and distribution of non-centromeric constitutive heterochromatin. The importance of the variation in constitutive heterochromatin for the phylogeny of the group is discussed. Comparison of the karyotypes in terms of G-banding patterns showed that C. pygmaea and C. emiliae differ from C. jacchus by a Robertsonian translocation and a paracentric inversion, whereas C. pygmaea and C. emiliae differ from each other by a reciprocal translocation between an acrocentric autosome and the short arm of the submetacentric chromosomes that distinguishes their karyotypes from that of C. jacchus. The possible evolutionary paths followed by the karyotypes of the three species are discussed.

Non-random distribution and co-localization of purine/pyrimidine-encoded information and transcriptional regulatory domains

In order to detect sequence-based information predictive for the location of eukaryotic transcriptional regulatory domains, the frequencies and distributions of the 36 possible purine/pyrimidine reverse complement hexamer pairs was determined for test sets of real and random sequences. The distribution of one of the hexamer pairs (RRYYRR/YYRRYY, referred to as M1) was further examined in a larger set of sequences (> 32 genes, 230 kb). Predominant clusters of M1 and the locations of eukaryotic transcriptional regulatory domains were found to be associated and non-randomly distributed along the DNA consistent with a periodicity of approximately 1.2 kb. In the context of higher ordered chromatin this would align promoters, enhancers and the predominant clusters of M1 longitudinally along one face of a 30 nm fiber. Using only information about the distribution of the M1 motif, 50-70% of a sequence could be eliminated as being unlikely to contain transcriptional regulatory domains with an 87% recovery of the regulatory domains present.

Evidence for selective evolution in codon usage in conserved amino acid segments of human alphaherpesvirus proteins

The genomes of human viruses herpes simplex 1 (HSV1) and varicella zoster (VZV), although similar in biology, largely concordant in gene order, and identical in many amino acid segments, differ widely in their genomic G + C (abbreviated S) content, which is high in HSV1 (68%) and low in VZV (46%). This paper analyzes several striking codon usage contrasts. The S difference in coding regions is dramatically large in codon site 3, S3, about 42%. The large difference in S3 is maintained at the same level in a subset of closely similar genes and even in corresponding identical amino acid blocks. A similar difference in S levels in silent site 1 (S1) is found in leucine and arginine. The difference in S3 levels occurs in every gene and in every multicodon amino acid form. The S difference also exists in amino acid usage, with HSV1 using significantly more codon types SSN, while VZV uses more codon types WWN (where W stands for A or T). The nonoverlapping and narrow histograms of S3 gene frequencies in both viruses suggest that the difference has arisen and been maintained by a process of selective rather than nonselective effects. This is in sharp contrast to the relatively large variance seen for highly similar genes in the human versus yeast analysis. Interpretations and hypotheses to explain the HSV1 vs VZV codon usage disparity relate to virus-host interactions, to the role of viral genes in DNA metabolism, to availability of molecular resources (molecular Gause exclusion principle), and to differences in genomic structure.

Evident diversity of codon usage patterns of human genes with respect to chromosome banding patterns and chromosome numbers; relation between nucleotide sequence data and cytogenetic data

The sequences of the human genome compiled in DNA databases are now about 10 megabase pairs (Mb), and thus the size of the sequences is several times the average size of chromosome bands at high resolution. By surveying this large quantity of data, it may be possible to clarify the global characteristics of the human genome, that is, correlation of gene sequence data (kb-level) to cytogenetic data (Mb-level). By extensively searching the GenBank database, we calculated codon usages in about 2000 human sequences. The highest G + C percentage at the third codon position was 97%, and that of about 250 sequences was 80% or more. The lowest G + C% was 27%, and that in about 150 sequences was 40% or less. A major portion of the GC-rich genes was found to be on special subsets of R-bands (T-bands and/or terminal R-bands). AT-rich genes, however, were mainly on G-bands or non-T-type internal R-bands. Average G + C% at the third position for individual chromosomes differed among chromosomes, and were related to T-band density, quinacrine dullness, and mitotic chiasmata density in the respective chromosomes.

Repair of UV-induced pyrimidine dimers in the individual genes Gart, Notch and white from Drosophila melanogaster cell lines

The excision repair of UV-induced pyrimidine dimers was investigated in three genes: Gart, Notch and white in a permanent Drosophila cell line Kc, derived from wild type Drosophila melanogaster embryonic cells. In this cell line Gart and Notch are actively transcribed, whereas white is not expressed. In all three genes UV-induced pyrimidine dimers were removed with the same rate and to the same extent: 60% removal within 16 hours, up to 80-100% in 24 hours after irradiation with 10 or 15 J/m2 UV. These kinetics are similar to the time course of dimer removal measured in the genome overall. No difference in repair of the inactive white locus compared to the active Gart and Notch genes was found. Similar results were obtained using a different wild type cell line, SL2, although repair appeared to be somewhat slower in this cell line. The results are discussed with respect to the data found for gene specific repair in other eukaryotic systems.

Co-amplification of L1 line elements with localised low copy repeats in Giemsa dark bands: implications for genome organisation

A repeat sequence island, located at the A3 Giemsa dark band on the mouse X chromosome and consisting of 50 copies of a localised long complex repeat unit (LCRU), features an unusually high concentration of L1 LINE repeat sequences juxtaposed and inserted within the LCRU. Sequence analysis of three independent genomic clones containing L1 LINE elements juxtaposed with the LCRU demonstrates a common junction sequence at the L1/LCRU boundary, suggesting that the high concentration of L1 LINE sequences in the repeat sequence island has arisen by association of an L1 element with an LCRU followed by amplification. The LCRU target site at this common junction sequence bears no resemblance to the target site of an L1 element inserted within one LCRU, indicating there is no specific preferential target site for L1 integration. We propose that co-amplification of L1 LINE elements with localised low copy repeat families throughout the genome could have a major effect on the chromosomal distribution of L1 LINE elements.

Systematic identification of LINE-1 repetitive DNA sequence differences having species specificity between Mus spretus and Mus domesticus

LINE-1 is a family of repetitive DNA sequences interspersed among mammalian genes. In the mouse haploid genome there are about 100,000 LINE-1 copies. We asked if the subspecies Mus spretus and Mus domesticus have developed species-specific LINE-1 subfamilies. Sequences from 14 M. spretus LINE-1 elements were obtained and compared to M. domesticus LINE-1 sequences. Using a molecular phylogenetic tree we identified several differences shared among a subset of young repeats in one or the other species as candidates for species-specific LINE-1 variants. Species specificity was tested using oligonucleotide probes complementary to each putative species-specific variant. When hybridized to genomic DNAs, single-variant probes detected an expanded number of elements in the expected mouse. In the other species these probes detected a smaller number of matches consistent with the average rate of random divergence among LINE-1 elements. It was further found that the combination of two species-specific sequence differences in the same probe reduced the detection background in the wrong species below our detection limit.

Deletion of chromosome 21 and normal intelligence: molecular definition of the lesion

Application of a method for the fine structure analysis of unbalanced chromosomal rearrangements using quantitative Southern blot analysis has established that an individual of normal intelligence and largely normal appearance has a significant interstitial deletion of chromosome 21. Using high resolution cytogenetic analysis and molecular analysis with five single copy DNA sequences unique to chromosome 21 and a probe for human SOD1 (CuZn, superoxide dismutase), we find that the deletion extends to the border of bands 21q11.1-11.2. and extends to the border of bands 21q21.2-q21.3. The latter border is established molecularly by the presence of two copies of SOD1, previously mapped to band 21q22.1, and of four single copy sequences known to be located distal to this region. The presence of SOD1 was confirmed by enzyme dosage analysis. These findings demonstrate that deletion of close to 20,000 kb of autosomal material is compatible with normal intelligence. Further, they suggest that chromosome 21 may include a large region of relative developmental neutrality whose molecular basis may now be investigated. Because of the limits of even high resolution cytogenetic analysis, fine structure molecular analyses of this type will be necessary to reliably detect and define similar small chromosomal deletions or insertions. The molecular definition of such aneuploidy provides the basis for increasing the resolution of the human physical genetic map.

Replication asynchrony between homologs 15q11.2: cytogenetic evidence for genomic imprinting

Replication kinetics of the Prader-Willi syndrome critical region (15q11.2) was investigated in seven normal healthy adult females using RBG replication bands. Replication asynchrony between homologs 15q11.2 was identified consistently in about 40% of cells in all individuals. It was limited to the stages in which Xp22, Xp11, Xq13 and Xq24/26 were visible in the late-replicating X chromosome. This asynchrony suggested that replication timing overlapped between 15q11.2 and the early replicating R-bands of the late X chromosome in some cells, and that the difference in replication timing between homologs was probably related to genomic imprinting; the latter has been suggested as a pathogenetic basis of Prader-Willi syndrome. As a result of an analysis of the proportions of asynchronous and synchronous cells in each replication stage, two types of cells were deduced providing 1:1 methylation mosaicism of genomic imprinting was assumed. The first type was composed of cells with normal replication in one homolog and delayed replication in the other. The second type was composed of cells with normal replication in both homologs. Our results provide cytogenetic evidence of methylation mosaicism for mammalian genomic imprinting.

Mammalian DNA replication: mutation biases and the mutation rate

Experimental studies have shown that the fidelity of DNA replication can be affected by the concentrations of free deoxyribonucleotides present in the cell. Replication of mammalian chromosomes is achieved using pools of newly-synthesized deoxyribonucleotides which fluctuate during the cell cycle. Since regions of mammalian chromosomes are replicated sequentially, there is the potential for differences among mammalian loci in both the relative and absolute frequencies of the various transitional and transversional mutations which may occur. Where these mutations are effectively neutral, at silent sites in genes and in non-coding sequences, this may result in different rates of evolution and in different base compositions, as have been observed in data from mammalian genes. A simple model of the DNA replication process is developed to describe how the mutation rate could be affected by the G + C contents of the deoxyribonucleotide pools and of the replicating DNA. Mutation rates are predicted to vary from locus to locus; only in the particular case of identical G + C contents in the DNA locus and the deoxyribonucleotide pools, and no proofreading, will the mutation rate be uniform over all loci.

Chromosomal organization of simple repeated DNA sequences used for DNA fingerprinting

Stretches of short, simple DNA sequences are widespread in all eukaryote genomes studied so far. Simple sequences are thought to undergo frequent expansion and deletion due to intrinsic genomic mechanisms. Some of the simple sequences were used successfully to detect hypervariable loci in various genomes. Hybridization experiments using synthetic probes not only revealed the informative simple repeats suitable for DNA fingerprinting in a particular species, but also reflected the wide range of distribution of the simple sequences among eukaryotes. The organization of these simple repetitive sequences at the chromosomal loci was investigated using in situ hybridization with chemically synthesized, pure oligonucleotide probes. Both biotin- and digoxigenin-attached probes detected specific chromosomal sites that are enriched in the respective simple-repeat blocks. Depending on the organism and probe used, accumulation of simple DNA sequences at individual or multiple sites on the chromosomes of different vertebrates could be demonstrated. The simple repetitive DNA sequences are located in different chromosomal regions (e.g., heterochromatin on the sex chromosomes, nucleolus organizer regions, and R-band sites), which are constrained considerably during evolution.

Co-polymer tracts in eukaryotic, prokaryotic, and organellar DNA

Large variations in DNA base composition and noticeable strand asymmetries are known to occur between different organisms and within different regions of the genomes of single organisms. Apparently such composition and sequence biases occur to fulfill structural rather than informational requirements. Here we report the wide occurrence of a more subtle biasing of DNA sequence that can have structural consequences: an increase or a suppression of the number of long tracts of two-base co-polymers. Strong biases were observed when the DNA sequences of the longest eukaryotic, prokaryotic, and organellar entries in the GenBank data base (totaling 773 kilobases) were analyzed for the number of occurrences of tracts of the two-base co-polymers (A,T)n, (G,C)n, and (A,C)n as a function of tract length. (The expression (A,T)n is used here to denote an uninterrupted tract, n nucleotides in length, of A and T bases in any proportion or order, terminated at each end by a G or C residue.) Characteristic differences are also observed in tract biases of eukaryotic vs. prokaryotic organisms.

A view of interphase chromosomes

Metaphase chromosomes are dynamically modified in interphase. This review focuses on how these structures can be modified, and explores the functional mechanisms and significance of these changes. Current analyses of genes often focus on relatively short stretches of DNA and consider chromatin conformations that incorporate only a few kilobases of DNA. In interphase nuclei, however, orderly transcription and replication can involve highly folded chromosomal domains containing hundreds of kilobases of DNA. Specific "junk" DNA sequences within selected chromosome domains may participate in more complex levels of chromosome folding, and may index different genetic compartments for orderly transcription and replication. Three-dimensional chromosome positions within the nucleus may also contribute to phenotypic expression. Entire chromosomes are maintained as discrete, reasonably compact entities in the nucleus, and heterochromatic coiled domains of several thousand kilobases can acquire unique three-dimensional positions in differentiated cell types. Some aspects of neoplasia may relate to alterations in chromosome structure at several higher levels of organization.

Molecular cytotaxonomy of primates by chromosomal in situ suppression hybridization

A new strategy for analyzing chromosomal evolution in primates is presented using chromosomal in situ suppression (CISS) hybridization. Biotin-labeled DNA libraries from flow-sorted human chromosomes are hybridized to chromosome preparations of catarrhines, platyrrhines, and prosimians. By this approach rearrangements of chromosomes that occurred during hominoid evolution are visualized directly at the level of DNA sequences, even in primate species with pronounced chromosomal shuffles.

Giant G+C% mosaic structures of the human genome found by arrangement of GenBank human DNA sequences according to genetic positions

To determine the overall variation in the G+C% distribution over long ranges of the human genome, DNA sequences of human genes, which were closely linked genetically or physically, were surveyed from the GenBank Data Bank. A total of 72 sequences longer than 2 kb, which were mutually linked within 500 kb, were identified. The sequences belonged to 17 linkage groups and were ordered in each group according to their genetic positions. Analyses of the G+C% distribution along the ordered sequences showed that sequences within each group almost always had similar G+C% levels, but those belonging to different groups often had different levels. Similar analyses of more distantly linked sequences (e.g., greater than 10 Mb) showed mosaic structures of G+C% distribution. These findings are consistent with predictions made from the "isochore" structures found by CsCl equilibrium centrifugation, in that the structures having homogeneous base compositions stretched over at least several hundred kilobases. A possible boundary of the giant G+C% mosaic structures was identified between X-linked G6PD and F8C.

Contrasts in codon usage of latent versus productive genes of Epstein-Barr virus: data and hypotheses

Epstein-Barr virus (EBV) has two different modes of existence: latent and productive. There are eight known genes expressed during latency (and hardly at all during the productive phase) and about 70 other ("productive") genes. It is shown that the EBV genes known to be expressed during latency display codon usage strikingly different from that of genes that are expressed during lytic growth. In particular, the percentage of S3 (G or C in codon site 3) is persistently lower (about 20%) in all latent genes than in nonlatent genes. Moreover, S3 is lower in each multicodon amino acid form. Also, the percentage of S in silent codon sites 1 of leucine and arginine is lower in latent than in nonlatent genes. The largest absolute differences in amino acid usage between latent and nonlatent genes emphasize codon types SSN and WWN (W means nucleotide A or T and N is any nucleotide). Two principal explanations to account for the EBV latent versus productive gene codon disparity are proposed. Latent genes have codon usage substantially different from that of host cell genes to minimize the deleterious consequences to the host of viral gene expression during latency. (Productive genes are not so constrained.) It is also proposed that the latency genes of EBV were acquired recently by the viral genome. Evidence and arguments for these proposals are presented.

A population genetic model of selection that maintains specific trinucleotides at a specific location

Periodic appearances of specific trinucleotides along the DNA sequence have been reported in the chicken core DNA, and the phenomenon has been suggested to be related to the supercoiling of DNA around nucleosomes. A population genetic model is constructed in which selection is operating to maintain specific trinucleotides at a specific location on the DNA sequence. Assuming low mutation rates, equilibrium probabilities of the appearances of respective trinucleotides were computed. Vague patterns appeared if the product of the effective size and the selection coefficient was 0.1-2.0. The genetic load and substitution rates in the equilibrium state were also computed. When the model was applied to the chicken DNA data, the product of the effective size and the selection coefficient was estimated to be 0.1-0.2. With this intensity of selection, the substitution rate was hardly different from that in the case without selection. However, the genetic load became fairly large. Considering the large number of times that DNA coils about nucleosomes, the number of trinucleotide sites must be very large, and thus the total load might be too large. Epistasis among these sites to reduce the total load is suggested to exist if selection is responsible for this periodic pattern observed in the chicken core DNA.

Chloroplast genes transferred to the nuclear plant genome have adjusted to nuclear base composition and codon usage

During plant evolution, some plastid genes have been moved to the nuclear genome. These transferred genes are now correctly expressed in the nucleus, their products being transported into the chloroplast. We compared the base compositions, the distributions of some dinucleotides and codon usages of transferred, nuclear and chloroplast genes in two dicots and two monocots plant species. Our results indicate that transferred genes have adjusted to nuclear base composition and codon usage, being now more similar to the nuclear genes than to the chloroplast ones in every species analyzed.

A unified model of eukaryotic chromosomes

A revised model of DNA packaging into chromosomes is presented. Its features are consistent with observed structural dimensions and the molecular periodicities related to transcription, replication and matrix attachment domains. The transitions between euchromatic, heterochromatic and metaphase states are explained simply. Molecular and physical properties of chromosomal bands, and their correlation with specific DNA sequence motifs are discussed.

Maintenance of function without selection: Alu sequences as "cheap genes"

Continued insertion into the genome of functional Alu sequences is expected to compensate for the functional eclipse of older sequences attributable to structural adulteration and can be presumed to establish a renewable store of functional sequences at a relatively elevated numerical level. This store of functional sequences could be maintained at almost no selective cost. A strategy of maintaining function in multiple sequence copies with selection limited to a very few master (source) sequences may be resorted to also by other types of DNA sequences that are generated repeatedly during evolution and that are spread over many sectors of the genome.