Sex-specific quantitative trait loci affecting longevity in Drosophila melanogaster

Senescence, the decline in survivorship and fertility with increasing age, is a near-universal property of organisms. Senescence and limited lifespan are thought to arise because weak natural selection late in life allows the accumulation of mutations with deleterious late-age effects that are either neutral (the mutation accumulation hypothesis) or beneficial (the antagonistic pleiotropy hypothesis) early in life. Analyses of Drosophila spontaneous mutations, patterns of segregating variation and covariation, and lines selected for late-age fertility have implicated both classes of mutation in the evolution of aging, but neither their relative contributions nor the properties of individual loci that cause aging in nature are known. To begin to dissect the multiple genetic causes of quantitative variation in lifespan, we have conducted a genome-wide screen for quantitative trait loci (QTLs) affecting lifespan that segregate among a panel of recombinant inbred lines using a dense molecular marker map. Five autosomal QTLs were mapped by composite interval mapping and by sequential multiple marker analysis. The QTLs had large sex-specific effects on lifespan and age-specific effects on survivorship and mortality and mapped to the same regions as candidate genes with fertility, cellular aging, stress resistance and male-specific effects. Late age-of-onset QTL effects are consistent with the mutation accumulation hypothesis for the evolution of senescence, and sex-specific QTL effects suggest a novel mechanism for maintaining genetic variation for lifespan.

Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics

The evolutionary relationships of 46 Shigella strains representing each of the serotypes belonging to the four traditional Shigella species (subgroups), Dysenteriae, Flexneri, Boydii, and Sonnei, were determined by sequencing of eight housekeeping genes in four regions of the chromosome. Analysis revealed a very similar evolutionary pattern for each region. Three clusters of strains were identified, each including strains from different subgroups. Cluster 1 contains the majority of Boydii and Dysenteriae strains (B1–4, B6, B8, B10, B14, and B18; and D3–7, D9, and D11–13) plus Flexneri 6 and 6A. Cluster 2 contains seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and Dysenteriae 2. Cluster 3 contains one Boydii strain (B12) and the Flexneri serotypes 1–5 strains. Sonnei and three Dysenteriae strains (D1, D8, and D10) are outside of the three main clusters but, nonetheless, are clearly within Escherichia coli. Boydii 13 was found to be distantly related to E. coli. Shigella strains, like the other pathogenic forms of E. coli, do not have a single evolutionary origin, indicating convergent evolution of Shigella phenotypic properties. We estimate the three main Shigella clusters to have evolved within the last 35,000 to 270,000 years, suggesting that shigellosis was one of the early infectious diseases of humans.

Identification of multiple quantitative trait loci linked to prion disease incubation period in mice

Polymorphisms in the prion protein gene are known to affect prion disease incubation times and susceptibility in humans and mice. However, studies with inbred lines of mice show that large differences in incubation times occur even with the same amino acid sequence of the prion protein, suggesting that other genes may contribute to the observed variation. To identify these loci we analyzed 1,009 animals from an F2 intercross between two strains of mice, CAST/Ei and NZW/OlaHSd, with significantly different incubation periods when challenged with RML scrapie prions. Interval mapping identified three highly significantly linked regions on chromosomes 2, 11, and 12; composite interval mapping suggests that each of these regions includes multiple linked quantitative trait loci. Suggestive evidence for linkage was obtained on chromosomes 6 and 7. The sequence conservation between the mouse and human genome suggests that identification of mouse prion susceptibility alleles may have direct relevance to understanding human susceptibility to bovine spongiform encephalopathy (BSE) infection, as well as identifying key factors in the molecular pathways of prion pathogenesis. However, the demonstration of other major genetic effects on incubation period suggests the need for extreme caution in interpreting estimates of variant Creutzfeldt–Jakob disease epidemic size utilizing existing epidemiological models.

Incipient speciation by sexual isolation in Drosophila: Concurrent evolution at multiple loci

Drosophila melanogaster from Zimbabwe and nearby regions shows strong but asymmetric sexual isolation from its cosmopolitan counterparts. By creating stable chromosome-substitution lines, earlier studies were able to show that the two major autosomes have very large effects on both male mating success and female mating preference. In this study, we genetically dissect this sexual isolation by recombination analysis between a whole-chromosome substitution line (which carries a Zimbabwe-derived third chromosome) and a strain with seven visible markers on that chromosome. Four loci are responsible for male mating success and three others are found to control female mating preference. Because male and female traits are not closely linked, their strong association among isofemale lines is most likely a reflection of sexual selection in nature. The results suggest that a large number of behavioral loci may evolve concurrently in the incipient stage of speciation before other aspects of reproductive isolation (such as hybrid sterility) have become evident. The results shed light on the population genetic processes underlying the formation of nascent species, as well as modes of speciation.

Integration of multiple repeats of geminiviral DNA into the nuclear genome of tobacco during evolution.

Integration of viral DNA into the host nuclear genome, although not unusual in bacterial and animal systems, has surprisingly not been reported for plants. We have discovered geminvirus-related DNA (GRD) sequences, in the form of distinct sets of multiple direct repeats comprising three related repeat classes, situated in a unique locus in the Nicotiana tabacum (tobacco) nuclear genome. The organization of these sequences is similar or identical in eight different tobacco cultivars we have examined. DNA sequence analysis reveals that each repeat has sequences most resembling those of the New World geminiviral DNA replication origin plus the adjacent AL1 gene, encoding the viral replication protein. We believe these GRD sequences originated quite recently in Nicotiana evolution through integration of geminiviral DNA by some combination of the processes of illegitimate recombination, amplification, deletions, and rearrangements. These events must have occurred in plant tissue that was subsequently able to contribute to meristematic tissue yielding gametes. GRD may have been retained in tobacco by selection or by random fixation in a small evolving population. Although we cannot detect transcription of these sequences, this does not exclude the possibility that they may originally have been expressed.

Mutually compensatory mutations during evolution of the tetramerization domain of tumor suppressor p53 lead to impaired hetero-oligomerization

We have measured the stability and stoichiometry of variants of the human p53 tetramerization domain to assess the effects of mutation on homo- and hetero-oligomerization. The residues chosen for mutation were those in the hydrophobic core that we had previously found to be critical for its stability but are not conserved in human p73 or p51 or in p53-related proteins from invertebrates or vertebrates. The mutations introduced were either single natural mutations or combinations of mutations present in p53-like proteins from different species. Most of the mutations were substantially destabilizing when introduced singly. The introduction of multiple mutations led to two opposite effects: some combinations of mutations that have occurred during the evolution of the hydrophobic core of the domain in p53-like proteins had additive destabilizing effects, whereas other naturally occurring combinations of mutations had little or no net effect on the stability, there being mutually compensating effects of up to 9.5 kcal/mol of tetramer. The triple mutant L332V/F341L/L344I, whose hydrophobic core represents that of the chicken p53 domain, was nearly as stable as the human domain but had impaired hetero-oligomerization with it. Thus, engineering of a functional p53 variant with a reduced capacity to hetero-oligomerize with wild-type human p53 can be achieved without any impairment in the stability and subunit affinity of the engineered homo-oligomer.

Evolution of microbial diversity during prolonged starvation

Models of evolutionary processes postulate that new alleles appear in populations through random spontaneous mutation. Alleles that confer a competitive advantage in particular environments are selected and populations can be taken over by individuals expressing these advantageous mutations. We have studied the evolutionary process by using Escherichia coli cultures incubated for prolonged periods of time in stationary phase. The populations of surviving cells were shown to be highly dynamic, even after many months of incubation. Evolution proceeded along different paths even when the initial conditions were identical. As cultures aged, the takeovers by fitter mutants were incomplete, resulting in the coexistence of multiple mutant forms and increased microbial diversity. Thus, the study of bacterial populations in stationary phase provides a model system for understanding the evolution of diversity in natural populations.

Fluorescence in situ hybridization analysis of keratinocyte growth factor gene amplification and dispersion in evolution of great apes and humans

Keratinocyte growth factor (KGF) is a member of the fibroblast growth factor family. Portions of the gene encoding KGF were amplified during primate evolution and are present in multiple nonprocessed copies in the human genome. Nucleotide analysis of a representative sampling of these KGF-like sequences indicated that they were at least 95% identical to corresponding regions of the KGF gene. To localize these sequences to specific chromosomal sites in human and higher primates, we used fluorescence in situ hybridization. In human, using a cosmid probe encoding KGF exon 1, we assigned the location of the KGF gene to chromosome 15q15–21.1. In addition, copies of KGF-like sequences hybridizing only with a cosmid probe encoding exons 2 and 3 were localized to dispersed sites on chromosome 2q21, 9p11, 9q12–13, 18p11, 18q11, 21q11, and 21q21.1. The distribution of KGF-like sequences suggests a role for alphoid DNA in their amplification and dispersion. In chimpanzee, KGF-like sequences were observed at five chromosomal sites, which were each homologous to sites in human, while in gorilla, a subset of four of these homologous sites was identified; in orangutan two sites were identified, while gibbon exhibited only a single site. The chromosomal localization of KGF sequences in human and great ape genomes indicates that amplification and dispersion occurred in multiple discrete steps, with initial KGF gene duplication and dispersion taking place in gibbon and involving loci corresponding to human chromosomes 15 and 21. These findings support the concept of a closer evolutionary relationship of human and chimpanzee and a possible selective pressure for such dispersion during the evolution of higher primates.

Origin and evolution of the slime molds (Mycetozoa)

The Mycetozoa include the cellular (dictyostelid), acellular (myxogastrid), and protostelid slime molds. However, available molecular data are in disagreement on both the monophyly and phylogenetic position of the group. Ribosomal RNA trees show the myxogastrid and dictyostelid slime molds as unrelated early branching lineages, but actin and β-tubulin trees place them together as a single coherent (monophyletic) group, closely related to the animal–fungal clade. We have sequenced the elongation factor-1α genes from one member of each division of the Mycetozoa, including Dictyostelium discoideum, for which cDNA sequences were previously available. Phylogenetic analyses of these sequences strongly support a monophyletic Mycetozoa, with the myxogastrid and dictyostelid slime molds most closely related to each other. All phylogenetic methods used also place this coherent Mycetozoan assemblage as emerging among the multicellular eukaryotes, tentatively supported as more closely related to animals + fungi than are green plants. With our data there are now three proteins that consistently support a monophyletic Mycetozoa and at least four that place these taxa within the “crown” of the eukaryote tree. We suggest that ribosomal RNA data should be more closely examined with regard to these questions, and we emphasize the importance of developing multiple sequence data sets.

Evolution of foraging behavior in Drosophila by density-dependent selection

One of the rare examples of a single major gene underlying a naturally occurring behavioral polymorphism is the foraging locus of Drosophila melanogaster. Larvae with the rover allele, forR, have significantly longer foraging path lengths on a yeast paste than do those homozygous for the sitter allele, fors. These variants do not differ in general activity in the absence of food. The evolutionary significance of this polymorphism is not as yet understood. Here we examine the effect of high and low animal rearing densities on the larval foraging path-length phenotype and show that density-dependent natural selection produces changes in this trait. In three unrelated base populations the long path (rover) phenotype was selected for under high-density rearing conditions, whereas the short path (sitter) phenotype was selected for under low-density conditions. Genetic crosses suggested that these changes resulted from alterations in the frequency of the fors allele in the low-density-selected lines. Further experiments showed that density-dependent selection during the larval stage rather than the adult stage of development was sufficient to explain these results. Density-dependent mechanisms may be sufficient to maintain variation in rover and sitter behavior in laboratory populations.

Intracellular gene transfer in action: Dual transcription and multiple silencings of nuclear and mitochondrial cox2 genes in legumes

The respiratory gene cox2, normally present in the mitochondrion, was previously shown to have been functionally transferred to the nucleus during flowering plant evolution, possibly during the diversification of legumes. To search for novel intermediate stages in the process of intracellular gene transfer and to assess the evolutionary timing and frequency of cox2 transfer, activation, and inactivation, we examined nuclear and mitochondrial (mt) cox2 presence and expression in over 25 legume genera and mt cox2 presence in 392 genera. Transfer and activation of cox2 appear to have occurred during recent legume evolution, more recently than previously inferred. Many intermediate stages of the gene transfer process are represented by cox2 genes in the studied legumes. Nine legumes contain intact copies of both nuclear and mt cox2, although transcripts could not be detected for some of these genes. Both cox2 genes are transcribed in seven legumes that are phylogenetically interspersed with species displaying only nuclear or mt cox2 expression. Inactivation of cox2 in each genome has taken place multiple times and in a variety of ways, including loss of detectable transcripts or transcript editing and partial to complete gene loss. Phylogenetic evidence shows about the same number (3–5) of separate inactivations of nuclear and mt cox2, suggesting that there is no selective advantage for a mt vs. nuclear location of cox2 in plants. The current distribution of cox2 presence and expression between the nucleus and mitochondrion in the studied legumes is probably the result of chance mutations silencing either cox2 gene.

Molecular evidence for multiple origins of woodiness and a New World biogeographic connection of the Macaronesian Island endemic Pericallis (Asteraceae: Senecioneae)

The prevalence of woody species in oceanic islands has attracted the attention of evolutionary biologists for more than a century. We used a phylogeny based on sequences of the internal-transcribed spacer region of nuclear ribosomal DNA to trace the evolution of woodiness in Pericallis (Asteraceae: Senecioneae), a genus endemic to the Macaronesian archipelagos of the Azores, Madeira, and Canaries. Our results show that woodiness in Pericallis originated independently at least twice in these islands, further weakening some previous hypotheses concerning the value of this character for tracing the continental ancestry of island endemics. The same data suggest that the origin of woodiness is correlated with ecological shifts from open to species-rich habitats and that the ancestor of Pericallis was an herbaceous species adapted to marginal habitats of the laurel forest. Our results also support Pericallis as closely related to New World genera of the tribe Senecioneae.

Stochastic processes strongly influence HIV-1 evolution during suboptimal protease-inhibitor therapy

It has long been assumed that HIV-1 evolution is best described by deterministic evolutionary models because of the large population size. Recently, however, it was suggested that the effective population size (Ne) may be rather small, thereby allowing chance to influence evolution, a situation best described by a stochastic evolutionary model. To gain experimental evidence supporting one of the evolutionary models, we investigated whether the development of resistance to the protease inhibitor ritonavir affected the evolution of the env gene. Sequential serum samples from five patients treated with ritonavir were used for analysis of the protease gene and the V3 domain of the env gene. Multiple reverse transcription–PCR products were cloned, sequenced, and used to construct phylogenetic trees and to calculate the genetic variation and Ne. Genotypic resistance to ritonavir developed in all five patients, but each patient displayed a unique combination of mutations, indicating a stochastic element in the development of ritonavir resistance. Furthermore, development of resistance induced clear bottleneck effects in the env gene. The mean intrasample genetic variation, which ranged from 1.2% to 5.7% before treatment, decreased significantly (P < 0.025) during treatment. In agreement with these findings, Ne was estimated to be very small (500–15,000) compared with the total HIV-1 RNA copy number. This study combines three independent observations, strong population bottlenecking, small Ne, and selection of different combinations of protease-resistance mutations, all of which indicate that HIV-1 evolution is best described by a stochastic evolutionary model.

The evolution of begging: Signaling and sibling competition

In many species, young solicit food from their parents, which respond by feeding them. Because of the difference in genetic make-up between parents and their offspring and the consequent conflict, this interaction is often studied as a paradigm for the evolution of communication. Existent theoretical models demonstrate that chick signaling and parent responding can be stable if solicitation is a costly signal. The marginal cost of producing stronger signals allows the system to converge to an equilibrium: young beg with intensity that reflects their need, and parents use this information to maximize their own inclusive fitness. However, we show that there is another equilibrium where chicks do not beg and parents’ provisioning effort is optimal with respect to the statistically probable distribution of chicks’ states. Expected fitness for parents and offspring at the nonsignaling equilibrium is higher than at the signaling equilibrium. Because nonsignaling is stable and it is likely to be the ancestral condition, we would like to know how natural systems evolved from nonsignaling to signaling. We suggest that begging may have evolved through direct sibling fighting before the establishment of a parental response, that is, that nonsignaling squabbling leads to signaling. In multiple-offspring broods, young following a condition-dependent strategy in the contest for resources provide information about their condition. Parents can use this information even though it is not an adaptation for communication, and evolution will lead the system to the signaling equilibrium. This interpretation implies that signaling evolved in multiple-offspring broods, but given that signaling is evolutionarily stable, it would also be favored in species which secondarily evolved single-chick broods.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.