Extreme rates and heterogeneity in insect DNA evolution

Summary DNA-DNA hybridization studies of insects, more specificallyDrosophila and cave crickets, have revealed interesting patterns of genome evolution that contrast markedly with what has been seen in other taxa, especially mammals and birds. Insect genomes are composed of sections of single-copy DNA with extreme variation in rates of evolutionary change. This variation is more extreme than between introns and exons; introns fall into the relatively conserved fraction of the genome. Attempts to calculate absolute rates of change inDrosophila DNA have all led to estimates some 5–10 times faster than those found in most vertebrates; this is true even for the more conservative part of the nuclear genome. Finally we point out that morphological similarity, chromosomal similarity, and/or ability to form interspecific hybrids is often associated with quite high levels of single-copy DNA divergence in insects as compared to mammals and birds.     

Intraspecific DNA divergence in Drosophila: a study on parthenogenetic D. mercatorum

Drosophila mercatorum is a species that can give rise to totally homozygous parthenogenetic strains. Using the technique of DNA-DNA hybridization, we have assessed the overall single-copy DNA differences among three independently derived strains that represent three independent genomes. Among strains, the average difference between homoduplex and heteroduplex median melting temperatures is 1.3 degrees C. This represents greater than or equal to 1.3% base-pair mismatch. Normalized percent of reassociation indicates further genetic differences, probably reflecting insertion/deletion differences and/or regions of the genome that are highly variable. This overall intraspecific genetic variation is higher than generally is thought to exist but is consistent with growing evidence of extensive DNA diversity within species of invertebrates. High intraspecific DNA variation may be correlated with rapid phyletic rates of evolution. Because of this high level of variation, the technique of DNA-DNA hybridization may be used to study intraspecific variation in invertebrates but is limited in its usefulness for higher systematic studies.      

Rates of DNA Evolution in Drosophila Depend on Function and Developmental Stage of Expression

DNA-sequence divergence of genes expressed in the embryonic stage was compared with the divergence of genes expressed in adults for 13 species of Drosophila representing various degrees of relatedness. DNA-DNA hybridization experiments were conducted using as tracers complementary DNA (cDNA) reversed transcribed from poly(A)(+) mRNA isolated from different developmental stages. The results indicate: (1) cDNA is less diverged than total single-copy DNA; (2) cDNA sequences are not in the rapidly evolving fraction of the single-copy genome of Drosophila; (3) early in evolutionary divergence embryonic messages are about half as diverged as adult messages; sequence data from some of the species compared indicate this is likely due to differences in rates of silent substitutions in genes expressed at different stages of development; and (4) at greater evolutionary distance, the differences in embryonic and adult messages disappear; this could be due to lineage-specific shifts in codon usage.     

DNA Sequence Comparison among Closely Related Drosophila Species in the MULLERI Complex

DNA hybridization was used to establish DNA sequence relationships among seven Drosophila species. Single-copy DNA was isolated from four species within the Drosophila mulleri complex, D. mojavensis, D. arizonensis, D. ritae and D. starmeri. These single-copy DNAs were used as tracers to be hybridized with each other and one additional member of the mulleri complex, D. aldrichi, a member of a closely related complex, D. hydei, and a distantly related species, D. melanogaster. Two methods have been used to determine the relatedness between these species: (1) the extent of duplex formed as measured by binding to hydroxyapatite and (2) the thermal stability of the duplexed DNA. Moderately repetitive DNA was purified from these species and used similarly to determine the divergence of this family of sequences. The rate of nucleotide substitution was estimated to be 0.2 +/-, 0.1% base pair change per million years for both single-copy and middle-repetitive DNAs. The size of the D. arizonensis genome, a representative of the mulleri complex, was calculated to be 2.2 X 10(8) base pairs from its kinetic complexity similar to that of D. hydei. The relative amounts (18%) and average reiteration frequency (100 copies) of the middle-repetitive DNA are similar for all Drosophila species studied. Finally, the data are presented in a phylogenetic tree.     

Rates and Patterns of Scndna and Mtdna Divergence within the Drosophila Melanogaster Subgroup

Levels of DNA divergence among the eight species of the Drosophila melanogaster subgroup and D. takahashii have been determined using the technique of DNA-DNA hybridization. Two types of DNA were used: single-copy nuclear DNA (scnDNA) and mitochondrial DNA (mtDNA). The major findings are: (1) A phylogeny has been derived for the group based on scnDNA which is congruent with chromosomal data, morphology, and behavior. The three homosequential species, simulans, sechellia, and mauritiana, are very closely related; the scnDNA divergence indicate the two island species are a monophyletic group. (2) The rates of change of scnDNA and mtDNA are not greatly different; if anything scnDNA evolves faster than mtDNA. (3) The rates of scnDNA evolution are not closely correlated to chromosomal (inversion) evolution. (4) The Drosophila genome appears to consist of two distinct classes of scnDNA with respect to rate of evolutionary change, a very rapidly evolving fraction and a relatively conservative fraction. (5) The absolute rate of change was estimated to be at least 1.7% nucleotide substitution per one million years. (6) DNA distance estimates based on restriction site variation are correlated with distances based on DNA-DNA hybridization, although the correlation is not very strong.     

Relationships of the chromosomal species in the eurasian mole rats of theSpalax ehrenbergi group as determined by DNA-DNA hybridization,and an estimate of the spalacid-murid divergence time

Summary DNA-DNA hybridization was used to measure the average genomic divergence among the four chromosomal species of the Eurasian mole rats belonging to theSpalax ehrenbergi complex (Rodentia: Spalacidae). The percent nucleotide substitutions in the single-copy nuclear DNA among the species ranged from 0 to 5%, suggesting that speciation has occurred with minor genomic changes in these animals. The youngest chromosomal species appear to differ by 0.2–0.6% base pair mismatch, which is only between one and three base differences in a 500-bp fragment. The interspecific values of percent nucleotide differences permit the recognition of two well-separated speciation events in theS. ehrenbergi complex, the older (of Lower Pleistocene age) having isolated the chromosomal species 2n=54 before the divergence of the three other species.DNA-DNA hybridization was also used to compare the Spalacinae (Eurasian mole rats), Murinae (Old World rats and mice), and Arvicolinae (voles and lemmings). These data enabled us to estimate the time of divergence of the spalacids at ca. 19 million years ago. The dates of divergence among the other rodent lineages, as predicted by DNA hybridization results, agree well with paleontological data. These dates of divergence are obtained by the relation between geological time and single-copy nuclear DNA change, a relation that was calibrated by Catzeflis et al. (1987) through the use of fossil Arvicolinae and Murinae data.     

Divergence of unique DNA sequences of bivalve mollusks from Mytilinae subfamily (Bivalvia mytilidae)

The DNA-DNA hybridization method was used to determine the divergence degree of unique sequences of five bivalvia species belonging to the Mytilinae subfamily. The matrix of delta Tm values for heteroduplexes of unique sequences was found which made it possible to define three phylogenetic branches within the subfamily. Under non-stringent hybridization conditions (55 degrees C, 0.5M PB) the divergence between species of any two branches was about 14% of nucleotide substitutions. The hybridization of [3H] unique sequences of Mytilinae with DNA fragments of Modiolus modiolus, a representative of the closest relative Modiolinae subfamily, showed that the divergence rate of unique sequences in two phylogenetic lines of the Mytilus genusis higher than in the line of the Crenomytilus genus I0.35-0.23% and 0.1% of nucleotide substitutions per one million years, respectively). According to the matrix of delta Tm values for five species of Mytilinae and Modiolus modiolus, a phylogenetic tree was built reflecting the differences between the divergence rates in different branches.     

Global population genetic structure and male-mediated gene flow in the green sea turtle (Chelonia mydas): analysis of microsatellite loci

We assessed the degree of population subdivision among global populations of green sea turtles, Chelonia mydas, using four microsatellite loci. Previously, a single-copy nuclear DNA study indicated significant male-mediated gene flow among populations alternately fixed for different mitochondrial DNA haplotypes and that genetic divergence between populations in the Atlantic and Pacific Oceans was more common than subdivisions among populations within ocean basins. Even so, overall levels of variation at single-copy loci were low and inferences were limited. Here, the markedly more variable microsatellite loci confirm the presence of male-mediated gene flow among populations within ocean basins. This analysis generally confirms the genetic divergence between the Atlantic and Pacific. As with the previous study, phylogenetic analyses of genetic distances based on the microsatellite loci indicate a close genetic relationship among eastern Atlantic and Indian Ocean populations. Unlike the single-copy study, however, the results here cannot be attributed to an artifact of general low variability and likely represent recent or ongoing migration between ocean basins. Sequence analyses of regions flanking the microsatellite repeat reveal considerable amounts of cryptic variation and homoplasy and significantly aid in our understanding of population connectivity. Assessment of the allele frequency distributions indicates that at least some of the loci may not be evolving by the stepwise mutation model.     

DNA DIVERGENCE AMONG HOMINOIDS

We have determined the degree of single-copy DNA divergence among the extant members of the Hominoidea employing the technique of DNA-DNA hybridization. The species studied include humans, two species of chimpanzees, gorillas, two subspecies of orangutans, and two species of gibbons; as an outgroup we have used a member of the Old World monkeys (Cercopithecidae), the baboon. Our methods are different from those previously used and allow us to control for two factors other than base-pair mismatch that can affect the thermal stability of DNA duplexes: the base composition and duplex length. In addition, we have studied more than one individual for most species and thus are able to assess the effect of intraspecific variation on phylogenetic conclusions. The results indicate that the closest extant relatives of humans are the chimpanzees. Gorillas are the next closest, followed by orangutans and gibbons. This result is strongly supported statistically, as there is virtually no overlap in measurements between different taxa. Our conclusions are in agreement with a growing amount of molecular evidence supporting this pattern of relatedness. The data behave as a reasonably good molecular clock, and we do not see an indication of slowdown in molecular evolution in the clade containing humans and African apes, contrary to what has been documented for protein-coding regions. Because of the clocklike nature of the results, we have estimated that the divergence of humans and chimpanzees occurred about 6–8 million years ago. Results from orangutans indicate that the Borneo and Sumatra populations are genetically distinct, about as different as the named species of chimpanzees.     

Testing the utility of internally transcribed spacer sequences in coral phylogenetics

Reef-building corals often possess high levels of intraindividual and intraspecific ribosomal DNA (rDNA) variation that is largely polyphyletic between closely related species. Polyphyletic rDNA phylogenies coupled with high intraindividual rDNA variation have been taken as evidence of introgressive hybridization in corals. Interpreting the data is problematic because the rDNA cluster evolves in a complex fashion and polyphyletic lineages can be generated by a variety of processes--such as incomplete lineage sorting and slow concerted evolution--in addition to hybridization. Using the genetically characterized Caribbean Acropora hybridization system, we evaluate how well rDNA data perform in revealing patterns of recent introgressive hybridization in contrast to genetic data from four single-copy loci. While the rDNA data are broadly consistent with the unidirectional introgression seen in other loci, we show that the phylogenetic signature of recent introgressive hybridization is obscured in the Caribbean Acropora by ancient shared rDNA lineages that predate the divergence of the species.     

Evolution of muntjac DNA

The extent of nuclear single-copy DNA divergence between Muntiacus reevesi and Muntiacus muntjak vaginalis (Cervidae), a species pair showing extreme karyotype differences but striking morphological similarity, is 2%, as judged from the thermal stability of interspecific DNA-DNA hybrids. A comparison of the total nuclear DNA reassociation kinetics of the two species indicates a reduction of lowly repetitive sequences in M. m. vaginalis.     

Rare and fleeting: an example of interspecific recombination in animal mitochondrial DNA

Recombination is thought to occur only rarely in animal mitochondrial DNA (mtDNA). However, detection of mtDNA recombination requires that cells become heteroplasmic through mutation, intramolecular recombination or 'leakage' of paternal mtDNA. Interspecific hybridization increases the probability of detecting mtDNA recombinants due to higher levels of sequence divergence and potentially higher levels of paternal leakage. During a study of historical variation in Atlantic salmon (Salmo salar) mtDNA, an individual with a recombinant haplotype containing sequence from both Atlantic salmon and brown trout (Salmo trutta) was detected. The individual was not an F1 hybrid but it did have an unusual nuclear genotype which suggested that it was a later-generation backcross. No other similar recombinant haplotype was found from the same population or three neighbouring Atlantic salmon populations in 717 individuals collected during 1948-2002. Interspecific recombination may increase mtDNA variability within species and can have implications for phylogenetic studies.     

Genotypic Diversity among Strains of Bradyrhizobium japonicum Belonging to Serogroup 110

Thirty-three strains of Bradyrhizobium japonicum within serogroup 110 were examined for genotypic diversity by using DNA-DNA hybridization analyses. The analysis of the DNA from 15 hydrogen-uptake-negative strains with the bradyrhizobial uptake hydrogenase probe pHU52 showed variation in degree of homology and restriction fragment length polymorphism of EcoRI-restricted DNA. Clustering analysis of the 33 strains on the basis of DNA-DNA hybridization analysis with four restriction enzymes and with the bradyrhizobial nodulation locus, pRJUT10, as probe indicated the existence of four groups of strains, which were less than 70% similar. Restriction digestion of genomic DNA with BamHI and DNA-DNA hybridization with pRJUT10 permitted classification of each of the strains according to a specific fingerprint pattern.     

DNA evolutionary rates in nine-primaried passerine birds

Differences in single-copy nuclear-DNA sequences among 13 species of passerine birds were measured using DNA-DNA hybridization. A matrix of pairwise dissimilarity values (delta mode distances) was constructed from analysis of fitted thermal dissociation curves. A least-squares method of phylogenetic estimation was used to construct two topologies from the distance matrix, one constraining branch lengths of sister taxa to be equal and the other permitting such lengths to vary. These topologies were identical in the pattern of branching of taxa, and the difference in their sums of squares was not statistically significant, suggesting that rates of DNA evolution in sister groups of nine- primaried oscines are equal. A nonparametric test for nonrandom variation in distances of sister groups to outgroup taxa revealed no statistically significant deviation from random variation that would be expected as a result of measurement error. However, the level of measurement error was such that rates of DNA evolution in sister taxa could vary by as much as 10% without being detected with the statistical methods used here.      

Identifying a single-copy DNA sequence associated with the expression of a heterochromatic gene, the light locus of Drosophila melanogaster

Although little is known about the molecular organization of most genes within heterochromatin, the unusual properties of these chromosomal regions suggest that genes therein may be organized and expressed very differently from those in euchromatin. We report here the cloning, by P transposon tagging, of sequences associated with the expression of the light locus, an essential gene found in the heterochromatin of chromosome 2 of Drosophila melanogaster. We conclude that this DNA is either a segment of the light locus, or a closely linked, heterochromatic sequence affecting its expression. While other functional DNA sequences previously described in heterochromatin have been repetitive, light gene function may be associated, at least in part, with single-copy DNA. This conclusion is based upon analysis of DNA from mutations and reversions induced by P transposable elements. The cloned region is unusual in that this single-copy DNA is embedded within middle-repetitive sequences. The in situ hybridization experiments also show that, unlike most other sequences in heterochromatin, this light-associated DNA evidently replicates in polytene chromosomes, but its diffuse hybridization signal may suggest an unusual chromosomal organization.     

Homology between the invertible deoxyribonucleic acid sequence that controls flagellar-phase variation in Salmonella sp. and deoxyribonucleic acid sequences in other organisms.

The invertible deoxyribonucleic acid (DNA) segment cloned from Salmonella sp. was radioactively labeled and used as a probe to search for homologous sequences by Southern hybridization. Only one copy of the invertible segment could be found on the Salmonella sp. genome. Partial sequence homology with the invertible region was detected in bacteriophage Mu and P1 DNA by low-stringency hybridization. Under these conditions, no homology was detected with Escherichia coli DNA. A strain of Salmonella sp. defective in phase variation carrying the vH2- allele was also analyzed by DNA-DNA hybridization. The results show that there is sequence divergence between diphasic and vH2- strains within the invertible sequence.     

Evolution of the RNA polymerase B' subunit gene (rpoB') in Halobacteriales: a complementary molecular marker to the SSU rRNA gene

Many prokaryotes have multiple ribosomal RNA operons. Generally, sequence differences between small subunit (SSU) rRNA genes are minor (<1%) and cause little concern for phylogenetic inference or environmental diversity studies. For Halobacteriales, an order of extremely halophilic, aerobic Archaea, within-genome SSU rRNA sequence divergence can exceed 5%, rendering phylogenetic assignment problematic. The RNA polymerase B' subunit gene (rpoB') is a single-copy conserved gene that may be an appropriate alternative phylogenetic marker for Halobacteriales. We sequenced a fragment of the rpoB' gene from 21 species, encompassing 15 genera of Halobacteriales. To examine the utility of rpoB' as a phylogenetic marker in Halobacteriales, we investigated three properties of rpoB' trees: the variation in resolution between trees inferred from the rpoB' DNA and RpoB' protein alignment, the degree of mutational saturation between taxa, and congruence with the SSU rRNA tree. The rpoB' DNA and protein trees were for the most part congruent and consistently recovered two well-supported monophyletic groups, the clade I and clade II haloarchaea, within a collection of less well resolved Halobacteriales lineages. A comparison of observed versus inferred numbers of substitution revealed mutational saturation in the rpoB' DNA data set, particularly between more distant species. Thus, the RpoB' protein sequence may be more reliable than the rpoB' DNA sequence for inferring Halobacteriales phylogeny. AU tests of tree selection indicated the trees inferred from rpoB' DNA and protein alignments were significantly incongruent with the SSU rRNA tree. We discuss possible explanations for this incongruence, including tree reconstruction artifact, differential paralog sampling, and lateral gene transfer. This is the first study of Halobacteriales evolution based on a marker other than the SSU rRNA gene. In addition, we present a valuable phylogenetic framework encompassing a broad diversity of Halobacteriales, in which novel sequences can be inserted for evolutionary, ecological, or taxonomic investigations.     

Evolutionary divergence in Dolichopoda cave crickets: A comparison of single copy DNA hybridization data with allozymes and morphometric distances

In this paper we attempt to investigate relationships between the amount of genetic divergence in nuclear genes and the degree of morphological differentiation for different sets of characters in Dolichopoda cave crickets. Six populations representing five Dolichopoda species from Central and Southern Italy have been studied. The overall genetic divergence at nuclear genes was estimated both by single copy DNA-DNA hybridization and allozyme frequencies at 26 loci. Euclidean distances for two multivariate sets of morphometric variables: one describing body and appendage morphology, the other male epiphallus shape. Results showed a close agreement between the branching patterns of ΔTm values from DNA hybridization and Nei's allozyme distance values. On the other hand, patterns of morphological divergence revealed independent trends, although the branching pattern based on epiphallus morphology matched to some extent the phylogenies inferred from molecular data. The relative value of molecular and morphological characters as reliable phylogenetic tracers was evaluated in relation to their dependence on evolutionary factors. Implications of these findings on the calibration of molecular clocks are also discussed. The absolute rate of molecular change based on scDNA was estimated to be at least 0.98% divergence/my/lineage. This result is in agreement with calibrations attempted on other insects. Estimates of time of divergence based on allozymes (Nei's D) were highly consistent with the estimate from geological data.     

The probability of a gene tree topology within a phylogenetic network with applications to hybridization detection

Gene tree topologies have proven a powerful data source for various tasks, including species tree inference and species delimitation. Consequently, methods for computing probabilities of gene trees within species trees have been developed and widely used in probabilistic inference frameworks. All these methods assume an underlying multispecies coalescent model. However, when reticulate evolutionary events such as hybridization occur, these methods are inadequate, as they do not account for such events. Methods that account for both hybridization and deep coalescence in computing the probability of a gene tree topology currently exist for very limited cases. However, no such methods exist for general cases, owing primarily to the fact that it is currently unknown how to compute the probability of a gene tree topology within the branches of a phylogenetic network. Here we present a novel method for computing the probability of gene tree topologies on phylogenetic networks and demonstrate its application to the inference of hybridization in the presence of incomplete lineage sorting. We reanalyze a Saccharomyces species data set for which multiple analyses had converged on a species tree candidate. Using our method, though, we show that an evolutionary hypothesis involving hybridization in this group has better support than one of strict divergence. A similar reanalysis on a group of three Drosophila species shows that the data is consistent with hybridization. Further, using extensive simulation studies, we demonstrate the power of gene tree topologies at obtaining accurate estimates of branch lengths and hybridization probabilities of a given phylogenetic network. Finally, we discuss identifiability issues with detecting hybridization, particularly in cases that involve extinction or incomplete sampling of taxa.     

PHYLOGENY, REPRODUCTIVE ISOLATION AND KIN RECOGNITION IN THE SOCIAL AMOEBA DICTYOSTELIUM PURPUREUM

Little is known about the population structure of social microorganisms, yet such studies are particularly interesting for the ways that genetic variation impacts their social evolution. Dictyostelium , a eukaryotic microbe widely used as a developmental model, has a social fruiting stage in which some formerly independent individuals die to help others. To assess genetic variation within the social amoeba Dictyostelium purpureum , we sequenced ∼4000 base pairs of ribosomal DNA (rDNA) from 37 isolates collected in Texas, Virginia, and Japan. Our analysis showed extensive genetic variation between populations and clear evidence of phylogenetic structure. We identified three major phylogenetic groups that were more different than other accepted species pairs. Tests using pairs of clones showed that both sexual macrocyst and asexual fruiting body formation were influenced by genetic divergence. Macrocysts were less likely to form between pairs of clones from different groups than from the same group. There was also a correlation between the genetic divergence of a pair of clones and their degree of mixing within fruiting bodies. These observations suggest that cryptic species might occur within D. purpureum and, more importantly, reveal how genetic variation impacts social interactions.