Microsatellite variation in populations of Drosophila pseudoobscura and Drosophila persimilis

We have isolated, characterized and mapped 33 dinucleotide, three trinucleotide and one tetranucleotide repeat loci from the four major chromosomes of Drosophila pseudoobscura. Average inferred repeat unit length of the dinucleotide repeats is 12 repeat units, similar to D. melanogaster. Assays of D. pseudoobscura and populations of its sibling species, D. persimilis, using 10 of these loci show extremely high levels of variation compared with similar studies of dinucleotide repeat variation in D. melanogaster populations. The high levels of variation are consistent with an average mutation rate of approximately 10(-6) per locus per generation and an effective population size of D. pseudoobscura approximately four times larger than that of D. melanogaster. Consistent with allozymes and nucleotide sequence polymorphism, the dinucleotide repeat loci reveal minimal structure across four populations of D. pseudoobscura. Finally, our preliminary recombinational mapping of 24 of these microsatellites suggests that the total recombinational genome size may be larger than previously inferred using morphological mutant markers.     

The distribution and frequency of microsatellite loci in Drosophila melanogaster

We report the results of a comprehensive search of Drosophila melanogaster DNA sequences in GenBank for di-, tri-, and tetranucleotide repeats of more than four repeat units, and a DNA library screen for dinucleotide repeats. Dinucleotide repeats are more abundant (66%) than tri- (30%) or tetranucleotide (4%) repeats. We estimate that 1917 dinucleotide repeats with 10 or more repeat units are present in the euchromatic D. melanogaster genome and, on average, they occur once every 60 kb. Relative to many other animals, dinucleotide repeats in D. melanogaster are short. Tri- and tetranucleotide repeats have even fewer repeat units on average than dinucleotide repeats. Our WorldWide Web site (http://www.bio.cornell.edu/genetics/aquadro/aquadro.html) posts the complete list of 1298 microsatellites (≥ five repeat units) identified from the GenBank search. We also summarize assay conditions for 70 D. melanogaster microsatellites characterized in previous studies and an additional 56 newly characterized markers.     

High density of long dinucleotide microsatellites in Drosophila subobscura

We isolated 96 dinucleotide repeats with five or more tandemly repeated units from a subgenomic Drosophila subobscura library. The mean repeat unit length of microsatellite clones in D. subobscura is 15, higher than that observed in other Drosophila species. Population variation was assayed in 32-40 chromosomes from Barcelona, Spain, using 18 randomly chosen microsatellite loci. Positive correlation between measures of variation and perfect repeat length measures (mean size, most common, and longest allele) is consistent with a higher mutation rate in loci with longer repeat units. Levels of microsatellite variation measured as variance in repeat number and heterozygosity in D. subobscura were similar to those of Drosophila pseudoobscura and higher than those of Drosophila melanogaster and Drosophila simulans. Our data suggest that higher levels of microsatellite variation, and possibly density, in D. subobscura compared with D. melanogaster are due to both a higher average effective population and a higher intrinsic slippage rate in the former species.     

Stable binding of Drosophila heat shock factor to head-to-head and tail-to-tail repeats of a conserved 5 bp recognition unit

The minimal DNA sequence required for the formation of a stable complex with Drosophila heat shock factor (HSF) in vitro is an inverted repeat of a 5 bp recognition unit, -GAA-. Surprisingly, both permutations of this 5 bp unit, head-to-head and tail-to-tail, bind to HSF with similar affinity and with striking 2-fold symmetry. HSF also binds to longer arrays of inverted 5 bp units, and the size of the HSF footprint increases with the addition of each 5 bp unit to these arrays. However, the electrophoretic mobility of the HSF-DNA complexes decreases most distinctly with the addition of every three 5 bp units. Cross-linking of purified HSF in the absence of DNA generates complexes with the sizes expected of HSF trimers. We propose that trimers of HSF bind to DNA and that the number of HSF subunits in direct contact with DNA is determined by the number of correctly positioned 5 bp recognition units.     

Sequences of three closely related variants of a complex satellite DNA diverge at specific domains

Major differences among the sequences of the repeat units of a very complex satellite DNA are located in domains which are sensitive to S1 nuclease under torsional stress, indicating that the domains assume unusual secondary or tertiary structures. Repeat units of the satellite, which accounts for 3% of the DNA of a land crab, have been inserted into pBR322 and the primary sequences of three cloned variants determined. The variants selected for sequencing include 1) RU (2089 base pairs (bp) ), representative of the average size of repeat units of cellular satellite; 2) TRU (1674 bp), truncated at an extra EcoRI site; and 3) EXT (2639 bp), extended by a 5-fold amplification of a 142-bp segment, one copy of which is present in RU and TRU (Bonnewell, V., Fowler, R.F., and Skinner, D.M. (1983) Science 221, 862-865). It appears that every copy of the satellite may be different and that the variants do not arise from cloning accidents. Extensive domains, as long as approximately 560 bp, are greater than 95% homologous among RU, TRU, and EXT; these conserved domains are composed of DNA whose base composition and sequences do not have remarkable features. By contrast, the sequences that comprise the divergent domains are unusually rich in 1) tracts of (dG X dC) 13-23 and arrangements of similar but not identical repetitive oligonucleotides or 2) alternating purines and pyrimidines (pu/py).     

Concerted evolution of primate alpha satellite DNA. Evidence for an ancestral sequence shared by gorilla and human X chromosome alpha satellite

To understand evolutionary events in the formation of higher-order repeat units in alpha satellite DNA, we have examined gorilla sequences homologous to human X chromosome alpha satellite. In humans, alpha satellite on the X chromosome is organized as a tandemly repeated, 2.0 x 10(3) base-pairs (bp) higher-order repeat unit, operationally defined by the restriction enzyme BamHI. Each higher-order repeat unit is composed of 12 tandem approximately 171 base-pair monomer units that have been classified into five distinct sequence homology groups. BamHI-digested gorilla genomic DNA hybridized with the cloned human 2 x 10(3) bp X alpha satellite repeat reveals three bands of sizes approximately 3.2 x 10(3), 2.7 x 10(3) and 2 x 10(3) bp. Multiple copies of all three repeat lengths have been isolated and mapped to the centromeric region of the gorilla X chromosome by fluorescence in situ hybridization. Long-range restriction mapping using pulsed-field gel electrophoresis shows that the 2.7 x 10(3) and 3.2 x 10(3) bp repeat arrays exist as separate but likely neighboring arrays on the gorilla X, each ranging in size from approximately 200 x 10(3) to 500 x 10(3) bp, considerably smaller than the approximately 2000 x 10(3) to 4000 x 10(3) bp array found on human X chromosomes. Nucleotide sequence analysis has revealed that monomers within all three gorilla repeat units can be classified into the same five sequence homology groups as monomers located within the higher-order repeat unit on the human X chromosome, suggesting that the formation of the five distinct monomer types predates the divergence of the lineages of contemporary humans and gorillas. The order of 12 monomers within the 2 x 10(3) and 2.7 x 10(3) bp repeat units from the gorilla X chromosome is identical with that of the 2 x 10(3) bp repeat unit from the human X chromosome, suggesting an ancestral linear arrangement and supporting hypotheses about events largely restricted to single chromosome types in the formation of alpha satellite higher-order repeat units.     

Intrahelical pseudoknots and interhelical associations mediated by mispaired human minisatellite DNA sequences in vitro.

The human minisatellite arrays, 33.6 and 33.15, consist of tandem reiterations of a 37-nucleotide (nt) and a 16-nt repeat unit sequence, respectively, both of which contain a majority of purine bases on one strand. Knot-like tertiary structures, which mapped to the cloned arrays, were observed by electron microscopy (EM) in homoduplex molecules produced by denaturation and reannealing in vitro. They result from a primary hybridization between misaligned repeat units of the array, forming a slipped-strand structure with staggered single-stranded DNA loops, followed by a secondary hybridization between repeat units in the two loops. Depending on the relative alignment of the loops when they hybridize, a particular form of intrahelical pseudoknot is produced. Theta-shaped, figure-of-eight, and bow-shaped structures were the most common conformational isomers observed in homoduplexes flattened into two dimensions during EM preparation. At the site of a bow-shaped structure, a conformation-dependent bend of approximately 60 degrees between the flanking DNA segments is induced; the other conformations generally do not deflect the line of the main DNA axis. Paired loops, similar to the bow-shaped structure, were apically situated in some supercoiled plasmids containing the 33.6 array. Both plasmids formed intermolecular associations, consisting of two (or more) homoduplex molecules held together at or immediately adjacent to a nexus which mapped to the minisatellite sequences. These associations might arise either by interhelical hybridization between arrays or by knot-like structures interfering with branch migration of chi-form Holliday junctions.     

DNA Sequence Evolution of the Amylase Multigene Family in Drosophila Pseudoobscura

The alpha-Amylase locus in Drosophila pseudoobscura is a multigene family of one, two or three copies on the third chromosome. The nucleotide sequences of the three Amylase genes from a single chromosome of D. pseudoobscura are presented. The three Amylase genes differ at about 0.5% of their nucleotides. Each gene has a putative intron of 71 (Amy1) or 81 (Amy2 and Amy3) bp. In contrast, Drosophila melanogaster Amylase genes do not have an intron. The functional Amy1 gene of D. pseudoobscura differs from the Amy-p1 gene of D. melanogaster at an estimated 13.3% of the 1482 nucleotides in the coding region. The estimated rate of synonymous substitutions is 0.398 +/- 0.043, and the estimated rate of nonsynonymous substitutions is 0.068 +/- 0.008. From the sequence data we infer that Amy2 and Amy3 are more closely related to each other than either is to Amy1. From the pattern of nucleotide substitutions we reason that there is selection against synonymous substitutions within the Amy1 sequence; that there is selection against nonsynonymous substitutions within the Amy2 sequence, or that Amy2 has recently undergone a gene conversion with Amy1; and that Amy3 is nonfunctional and subject to random genetic drift.     

Evolution of Repeated Sequence Arrays in the D-Loop Region of Bat Mitochondrial DNA

Analysis of mitochondrial DNA control region sequences from 41 species of bats representing 11 families revealed that repeated sequence arrays near the tRNA-Pro gene are present in all vespertilionine bats. Across 18 species tandem repeats varied in size from 78 to 85 bp and contained two to nine repeats. Heteroplasmy ranged from 15% to 63%. Fewer repeats among heteroplasmic than homoplasmic individuals in a species with up to nine repeats indicates selection may act against long arrays. A lower limit of two repeats and more repeats among heteroplasmic than homoplasmic individuals in two species with few repeats suggests length mutations are biased. Significant regressions of heteroplasmy, θ and π, on repeat number further suggest that repeat duplication rate increases with repeat number. Comparison of vespertilionine bat consensus repeats to mammal control region sequences revealed that tandem repeats of similar size, sequence and number also occur in shrews, cats and bighorn sheep. The presence of two conserved protein-binding sequences in all repeat units indicates that convergent evolution has occurred by duplication of functional units. We speculate that D-loop region tandem repeats may provide signal redundancy and a primitive repair mechanism in the event of somatic mutations to these binding sites.     

Variation in Activities of Amylase Allozymes Associated with Chromosome Inversions in DROSOPHILA PSEUDOOBSCURA, D. PERSIMILIS and D. MIRANDA

Different electrophoretic alleles of amylase show associations with particular chromosome 3 inversions in D. pseudoobscura and D. persimilis. Relative adult amylase activities were compared in 37, 37 and 10 strains of D. pseudoobscura, D. persimilis and D. miranda, respectively. Strains carrying the same electrophoretic allele were compared by crossing these lines individually to a reference strain carrying a different electrophoretic mobility allele. This procedure allows comparisons among species, inversions, electromorphs and strains for genetic variation in amylase activity. F2 analysis established that the activity variation co-segregates with the structural amylase locus. This type of variation could be due to either structural gene differences or differences in closely linked, cis-acting regulatory regions. Variation has been detected among and within electrophoretic mobility classes. Moreover, this variation is clearly nonrandom and reveals more of the genetic structure associated with the chromosomal inversion phylogeny of D. pseudoobscura and D. persimilis. ----Some of the findings are: (1) Similar electromorphs in D. pseudoobscura and D. persimilis usually show different activities. These species show nearly complete differentiation of amylase alleles, based on activities. (2) D. persimilis has the broadest range of variation in amylase activity, about four-fold between the highest and lowest alleles. D. pseudoobscura and D. miranda are also polymorphic for activity, but have more constrained ranges of variation. D. miranda alleles show on the average about four times the activity of D. pseudoobscura alleles. (3) Some association of electrophoretic mobility and activity has been found. Alleles 1.09 of D. persimilis, as well as 1.43 and 1.55 of D. miranda, have relatively high activity. It may be that these high activity alleles are part of an adaptation to cooler habitats. (4) Within electrophoretic classes, associations of activities with inversions have been found. These are especially strong in D. persimilis. The 1.00 alleles in the ST, KL, MD and WT inversions, the 0.92 allele in the ST and MD inversions and the 1.09 allele in the WT and KL inversions have levels of activities that depend upon the arrangement in which they are located. These results demonstrate that suppression of recombination in inversion heterokaryotypes can result in extensive genic divergence between inversions.     

Centromere 3 specific tandem repeat from Chironomus pallidivittatus

A 155-bp tandem repeat was previously reported to be present in all centromeric regions of the dipteran Chironomus pallidivittatus. We have now isolated a second centromere specific tandem repeat, 375 bp long. Two blocks were found of the new unit, differing in size, probably representing allelic forms. The repeat is present only in chromosome 3, bordering 155-bp repeat arrays. There are about 100 repeats per genome, compared to 1300 units for the 155-bp repeat. The two units contain an identical 9-bp sequence which can form target-site duplications flanking a short mobile element, Cp1. An inversion within the tandem array was isolated, the breakpoint of which is within the 9-bp target sequence. Another short shared motif, 10-bp long, is also present at the insertion site for a mobile element. The two repeat units are similar in having long regions with more than 80% AT and an overall high AT content. Received: 16 February 1998; in revised form: 4 August 1998 / Accepted: 5 August 1998     

The complete cDNA and polypeptide sequences of human erythroid alpha-spectrin.

Overlapping human erythroid alpha-spectrin cDNA clones were isolated from lambda gt11 libraries constructed from cDNAs of human fetal liver and erythroid bone marrow. The composite 8001-base pair (bp) cDNA nucleotide sequence contains 187-bp 5'- and 528-bp 3'-untranslated regions and has a single long open reading frame of 7287 bp that encodes a polypeptide of 2429 residues. As previously described (Speicher, D. W., and Marchesi, V. T. (1984) Nature 311, 177-180), spectrin is composed largely of homologous 106-amino acid repeat units. From the amino acid sequence deduced from the cDNA, alpha-spectrin can be divided into 22 segments. Segments 1-9 and 12-19 are homologous and can therefore be considered repeats; the average number of identical residues in pairwise comparisons of these repeats is 22 out of 106, or 21%. Of these 17 repeats, 11 are exactly 106 amino acids in length, whereas five others differ from this length by a single residue. Segments 11, 20, and 21, although less homologous, appear to be related to the more highly conserved repeat units. The very N-terminal 22 residues, segment 10, which is atypical both in length and sequence, and the C-terminal 150 residues in segment 22 appear to be unrelated to the conserved repeat units. The sequence of the erythroid alpha-spectrin polypeptide chain is compared to that of human alpha-fodrin and chicken alpha-actinin to which it is related. alpha-Spectrin is more distantly related to dystrophin.     

Interspecific comparison of the period gene of Drosophila reveals large blocks of non-conserved coding DNA.

We have cloned and sequenced the coding region of the period (per) gene from Drosophila pseudoobscura and D. virilis. A comparison with that of D. melanogaster reveals that the conceptual translation products consist of interspersed blocks of conserved and non-conserved amino acid sequence. The non-conserved portion, comprising approximately 33% of the protein sequence, includes the perfect Thr-Gly repeat of D. melanogaster, which is absent from the D. pseudoobscura and D. virilis proteins. Based on these observations and cross-species transformation experiments, we suggest that the interspecific variability in the per primary amino acid sequence contributes to the control of species-specific behaviors.