Accelerated evolution and Muller's rachet in endosymbiotic bacteria.

Many bacteria live only within animal cells and infect hosts through cytoplasmic inheritance. These endosymbiotic lineages show distinctive population structure, with small population size and effectively no recombination. As a result, endosymbionts are expected to accumulate mildly deleterious mutations. If these constitute a substantial proportion of new mutations, endosymbionts will show (i) faster sequence evolution and (ii) a possible shift in base composition reflecting mutational bias. Analyses of 16S rDNA of five independently derived endosymbiont clades show, in every case, faster evolution in endosymbionts than in free-living relatives. For aphid endosymbionts (genus Buchnera), coding genes exhibit accelerated evolution and unusually low ratios of synonymous to nonsynonymous substitutions compared to ratios for the same genes for enterics. This concentration of the rate increase in nonsynonymous substitutions is expected under the hypothesis of increased fixation of deleterious mutations. Polypeptides for all Buchnera genes analyzed have accumulated amino acids with codon families rich in A+T, supporting the hypothesis that substitutions are deleterious in terms of polypeptide function. These observations are best explained as the result of Muller's ratchet within small asexual populations, combined with mutational bias. In light of this explanation, two observations reported earlier for Buchnera, the apparent loss of a repair gene and the overproduction of a chaperonin, may reflect compensatory evolution. An alternative hypothesis, involving selection on genomic base composition, is contradicted by the observation that the speedup is concentrated at nonsynonymous sites.

A pre-Columbian Y chromosome-specific transition and its implications for human evolutionary history.

A polymorphic C-->T transition located on the human Y chromosome was found by the systematic comparative sequencing of Y-specific sequence-tagged sites by denaturing high-performance liquid chromatography. The results of genotyping representative global indigenous populations indicate that the locus is polymorphic exclusively within the Western Hemisphere. The pre-Columbian T allele occurs at > 90% frequency within the native South and Central American populations examined, while its occurrence in North America is approximately 50%. Concomitant genotyping at the polymorphic tetranucleotide microsatellite DYS19 locus revealed that the C-->T mutation displayed significant linkage disequilibrium with the 186-bp allele. The data suggest a single origin of linguistically diverse native Americans with subsequent haplotype differentiation within radiating indigenous populations as well as post-Columbian European and African gene flow. The mutation may have originated either in North America at a very early time during the expansion or before it, in the ancestral population(s) from which all Americans may have originated. The analysis of linkage of the DYS199 and the DYS19 tetranucleotide loci suggests that the C-->T mutation may have occurred around 30,000 years ago. We estimate the nucleotide diversity over 4.2 kb of the nonrecombining portion of the Y chromosome to be 0.00014. compared to autosomes, the majority of variation is due to the smaller effective population size of the Y chromosome rather than selective sweeps. There begins to emerge a pattern of pronounced geographical localization of Y-specific nucleotide substitution polymorphisms.

Genetic absolute dating based on microsatellites and the origin of modern humans.

We introduce a new genetic distance for microsatellite loci, incorporating features of the stepwise mutation model, and test its performance on microsatellite polymorphisms in humans, chimpanzees, and gorillas. We find that it performs well in determining the relations among the primates, but less well than other distance measures (not based on the stepwise mutation model) in determining the relations among closely related human populations. However, the deepest split in the human phylogeny seems to be accurately reconstructed by the new distance and separates African and non-African populations. The new distance is independent of population size and therefore allows direct estimation of divergence times if the mutation rate is known. Based on 30 microsatellite polymorphisms and a recently reported average mutation rate of 5.6 x 10(-4) at 15 dinucleotide microsatellites, we estimate that the deepest split in the human phylogeny occurred about 156,000 years ago. Unlike most previous estimates, ours requires no external calibration of the rate of molecular evolution. We can use such calibrations, however, to test our estimate.

A genetic record of population isolation in pocket gophers during Holocene climatic change

A long-standing question in Quaternary paleontology is whether climate-induced, population-level phenotypic change is a result of large-scale migration or evolution in isolation. To directly measure genetic variation through time, ancient DNA and morphologic variation was measured over 2,400 years in a Holocene sequence of pocket gophers (Thomomys talpoides) from Lamar Cave, Yellowstone National Park, Wyoming. Ancient specimens and modern samples collected near Lamar Cave share mitochondrial cytochrome b sequences that are absent from adjacent localities, suggesting that the population was isolated for the entire period. In contrast, diastemal length, a morphologic character correlated with body size and nutritional level, changed predictably in response to climatic change. Our results demonstrate that small mammal populations can experience the long-term isolation assumed by many theoretical models of microevolutionary change.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

Sex ratio adjustment in relation to paternal attractiveness in a wild bird population.

When the relative fitness of sons and daughters differs, sex-allocation theory predicts that it would be adaptive for individuals to adjust their investment in different sexes of offspring. Sex ratio adjustment by females in response to the sexual attractiveness of their mate would be an example of this. In vertebrates the existence of this form of sex ratio adjustment is controversial and may be confounded with sex-biased mortality, particularly in sexually size-dimorphic species. Here we use PCR amplification of a conserved W-chromosome-linked gene to show that the sex ratio within broods of a natural population of sexually size-monomorphic collared flycatchers Ficedula albicollis is related to the size of their father's forehead patch, a heritable secondary sexual character implicated in female choice. Experimental manipulations of paternal investment, which influence the size of his character in future breeding attempts, result in corresponding changes in the sex ratio of offspring born to males in those breeding attempts. In contrast, manipulations of maternal investment have no effect on future sex ratios, and there is no relationship between variables predicting the reproductive value of the brood and nestling sex ratio. Analysis of recruitment of offspring reveals similar patterns of sex ratio bias. The results suggest that female collared flycatchers be able to adjust the sex ratio of eggs ovulated in response to the phenotype of their mate. This finding is most consistent with "genetic quality" models of sexual selection.

Correlates of the desired family size among Indian communities.

The People of India database of the Anthropological Survey of India documents 631 cultural, ecological, and economic traits of the 4635 communities to which the entire Indian population is assigned. Focusing on 1342 communities of South India, we looked for correlates of low (1 or 2 children) and high (4 or more children) desired family size (DFS) reported as the norm for any given community by key informants. We found 10 cultural and 18 economic traits to be significantly correlated to high DFS and 21 cultural and 9 economic traits to low DFS. The economic traits so identified are compatible with high family size being desired by parents who have little capability of investing in quality of offspring, but whose children contribute economically from an early age. In contrast, communities desiring low family size are part of the modern intensive agriculture/organized industry/services sector and invest heavily in educating their children. A composite index based on 27 economic traits (CEI) has a high predictive value with respect to the DFS for the entire set of 4635 Indian communities. The 31 cultural traits highly correlated to high or low DFS constitute 5 clusters that can be identified as characterizing scheduled tribes, scheduled castes, rural and landless lower castes, urban upper castes, and Moslems. Whereas economic traits have similar influence on DFS within each of these ethnic categories, Moslems demonstrate a significantly higher DFS for lower values of CEI.

Population dynamics of flaviviruses revealed by molecular phylogenies.

The phylogeny of 123 complete envelope gene sequences was reconstructed in order to understand the evolution of tick- and mosquito-borne flaviviruses. An analysis of phylogenetic tree structure reveals a continual and asymmetric branching process in the tick-borne flaviviruses, compared with an explosive radiation in the last 200 years in viruses transmitted by mosquitoes. The distinction between these two viral groups probably reflects differences in modes of dispersal, propagation, and changes in the size of host populations. The most serious implication of this work is that growing human populations are being exposed to an expanding range of increasingly diverse viral strains.

Neoplastic development: paradoxical relation between impaired cell growth at low population density and excessive growth at high density.

The role of heritable, population-wide cell damage in neoplastic development was studied in the 28 L subline of NIH 3T3 cells. These cells differ from the 17(3c) subline used previously for such studies in their lower frequency of "spontaneous" transformation at high population density and their greater capacity to produce large, dense transformed foci. Three cultures of the 28 L subline of NIH 3T3 cells were held under the constraint of confluence for 5 wk (5 wk 1 degree assay) and then assayed twice in succession (2 degrees and 3 degrees assays) for transformed foci and saturation density. After the 2 degrees assay, the cells were also passaged at low density to determine their exponential growth rates and cloned to determine the size and morphological features of the colonies. Concurrent measurements were made in each case with control cells that had been kept only in frequent low-density passages and cells that had been kept at confluence for only 2 wk (2 wk 1 degree). Two of the three cultures transferred from the 2 degrees assay of the 5 wk 1 degree cultures produced light transformed foci, and the third produced dense foci. The light focus-forming cultures grew to twice the control saturation density in their 2 degrees assay and 6-8 times the control density in the 3 degrees assay; saturation densities for the dense focus formers were about 10 times the control values in both assays. All three of the cultures transferred from the 2 degrees assay of the 5 wk 1 degree cultures multiplied at lower rates than controls at low densities, but the dense focus formers multiplied faster than the light focus formers. The reduced rates of multiplication of the light focus formers persisted for > 50 generations of exponential multiplication at low densities. Isolated colonies formed from single cells of the light focus formers were of a lower population density than controls; colonies formed by the dense focus formers were slightly denser than the controls but occupied only half the area. A much higher proportion of the colonies from the 5 wk 1 degree cultures than the controls consisted of giant cells or mixtures of giant and normal-appearing cells. The results reinforce the previous conclusion that the early increases in saturation density and light focus formation are associated with, and perhaps caused by, heritable, population-wide damage to cells that is essentially epigenetic in nature. The more advanced transformation characterized by large increases in saturation density and dense focus formation could have originated from rare genetic changes, such as chromosome rearrangements, known to occur at an elevated frequency in cells destabilized by antecedent cellular damage.

Heritable, population-wide damage to cells as the driving force of neoplastic transformation.

Prolonged incubation of NIH 3T3 cells under the growth constraint of confluence results in the death of some cells in a manner suggestive of apoptosis. Successive rounds of prolonged incubation at confluence of the surviving cells produce increasing neoplastic transformation in the form of increments in saturation density and transformed focus formation. Cells from the postconfluent cultures are given a recovery period of various lengths to remove the direct inhibitory effect of confluence before their growth properties are studied. It is found that with each round of confluence the exponential growth rate of the cells at low densities gets lower and the size of isolated colonies of the same cells shows a similar progressive reduction. The decreased growth rate of cells from the third round of confluence persists for > 60 generations of growth at low density. The proportion of colonies containing giant cells is much higher after a 2-day recovery from confluence than after a 7-day recovery. Retardation of growth at low density and increased saturation density appear to be two sides of the same coin: both occur in the entire population of cells and precede the formation of transformed foci. We propose that the slowdown in growth and the formation of giant cells result from heritable damage to the cells, which in turn drives their transformation. Similar results have been reported for the survivors of x-irradiation and of treatment with chemical carcinogens and are associated with the aging process in animals. We suggest that these changes result from free radical damage to membrane lipids with particular damage to lysosomes. Proteases and nucleases would then be released to progressively modify the growth behavior and genetic stability of the cells toward autonomous proliferation.

Quantitative analysis of senile plaques in Alzheimer disease: observation of log-normal size distribution and molecular epidemiology of differences associated with apolipoprotein E genotype and trisomy 21 (Down syndrome).

The discovery that the epsilon 4 allele of the apolipoprotein E (apoE) gene is a putative risk factor for Alzheimer disease (AD) in the general population has highlighted the role of genetic influences in this extremely common and disabling illness. It has long been recognized that another genetic abnormality, trisomy 21 (Down syndrome), is associated with early and severe development of AD neuropathological lesions. It remains a challenge, however, to understand how these facts relate to the pathological changes in the brains of AD patients. We used computerized image analysis to examine the size distribution of one of the characteristic neuropathological lesions in AD, deposits of A beta peptide in senile plaques (SPs). Surprisingly, we find that a log-normal distribution fits the SP size distribution quite well, motivating a porous model of SP morphogenesis. We then analyzed SP size distribution curves in genotypically defined subgroups of AD patients. The data demonstrate that both apoE epsilon 4/AD and trisomy 21/AD lead to increased amyloid deposition, but by apparently different mechanisms. The size distribution curve is shifted toward larger plaques in trisomy 21/AD, probably reflecting increased A beta production. In apoE epsilon 4/AD, the size distribution is unchanged but the number of SP is increased compared to apoE epsilon 3, suggesting increased probability of SP initiation. These results demonstrate that subgroups of AD patients defined on the basis of molecular characteristics have quantitatively different neuropathological phenotypes.