Deleterious mutations destabilize ribosomal RNA in endosymbiotic bacteria

In populations that are small and asexual, mutations with slight negative effects on fitness will drift to fixation more often than in large or sexual populations in which they will be eliminated by selection. If such mutations occur in substantial numbers, the combined effects of long-term asexuality and small population size may result in substantial accumulation of mildly deleterious substitutions. Prokaryotic endosymbionts of animals that are transmitted maternally for very long periods are effectively asexual and experience smaller effective population size than their free-living relatives. The contrast between such endosymbionts and related free-living bacteria allows us to test whether a population structure imposing frequent bottlenecks and asexuality does lead to an accumulation of slightly deleterious substitutions. Here we show that several independently derived insect endosymbionts, each with a long history of maternal transmission, have accumulated destabilizing base substitutions in the highly conserved 16S rRNA. Stabilities of Domain I of this subunit are 15–25% lower in endosymbionts than in closely related free-living bacteria. By mapping destabilizing substitutions onto a reconstructed phylogeny, we show that decreased ribosomal stability has evolved separately in each endosymbiont lineage. Our phylogenetic approach allows us to demonstrate statistical significance for this pattern: becoming endosymbiotic predictably results in decreased stability of rRNA secondary structure.

Novel gene conversion between X-Y homologues located in the nonrecombining region of the Y chromosome in Felidae (Mammalia)

Genes located on the mammalian Y chromosome outside of the pseudoautosomal region do not recombine with those on the X and are predicted to either undergo selection for male function or gradually degenerate because of an accumulation of deleterious mutations. Here, phylogenetic analyses of X-Y homologues, Zfx and Zfy, among 26 felid species indicate two ancestral episodes of directed genetic exchange (ectopic gene conversion) from X to Y: once during the evolution of pallas cat and once in a common predecessor of ocelot lineage species. Replacement of the more rapidly evolving Y homologue with the evolutionarily constrained X copy may represent a mechanism for adaptive editing of functional genes on the nonrecombining region of the mammalian Y chromosome.

Male fitness increases when females are eliminated from gene pool: Implications for the Y chromosome

Because the two sexes share a common gene pool while performing many different biological functions, mutations benefiting one sex may not accumulate due to counter selection in the other sex. In these experiments 99% of a haploid genome of Drosophila melanogaster was constrained to segregate like a male-limited Y chromosome for 41 generations, thereby eliminating potential counter selection in females. The synthetic Y chromosomes rapidly accumulated genetic variation that increased male fitness and decreased female fitness. The survival and fertility of females declined when they were mated to males expressing the synthetic Y chromosomes. These results suggests that opposing selection between the sexes may substantially interfere with sex-specific adaptation. They also demonstrate how intersexual evolutionary conflict can lead to perpetual degeneration of the Y via genetic hitchhiking of deleterious mutations.

Metapopulation extinction caused by mutation accumulation

Theory suggests that the risk of extinction by mutation accumulation can be comparable to that by environmental stochasticity for an isolated population smaller than a few thousand individuals. Here we show that metapopulation structure, habitat loss or fragmentation, and environmental stochasticity can be expected to greatly accelerate the accumulation of mildly deleterious mutations, lowering the genetic effective size to such a degree that even large metapopulations may be at risk of extinction. Because of mutation accumulation, viable metapopulations may need to be far larger and better connected than would be required under just stochastic demography.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Reproductive systems and evolution in vascular plants

Differences in the frequency with which offspring are produced asexually, through self-fertilization and through sexual outcrossing, are a predominant influence on the genetic structure of plant populations. Selfers and asexuals have fewer genotypes within populations than outcrossers with similar allele frequencies, and more genetic diversity in selfers and asexuals is a result of differences among populations than in sexual outcrossers. As a result of reduced levels of diversity, selfers and asexuals may be less able to respond adaptively to changing environments, and because genotypes are not mixed across family lineages, their populations may accumulate deleterious mutations more rapidly. Such differences suggest that selfing and asexual lineages may be evolutionarily short-lived and could explain why they often seem to be of recent origin. Nonetheless, the origin and maintenance of different reproductive modes must be linked to individual-level properties of survival and reproduction. Sexual outcrossers suffer from a cost of outcrossing that arises because they do not contribute to selfed or asexual progeny, whereas selfers and asexuals may contribute to outcrossed progeny. Selfing and asexual reproduction also may allow reproduction when circumstances reduce opportunities for a union of gametes produced by different individuals, a phenomenon known as reproductive assurance. Both the cost of outcrossing and reproductive assurance lead to an over-representation of selfers and asexuals in newly formed progeny, and unless sexual outcrossers are more likely to survive and reproduce, they eventually will be displaced from populations in which a selfing or asexual variant arises.

Human brain contains high levels of heteroplasmy in the noncoding regions of mitochondrial DNA.

We have analyzed the level of intraindividual sequence variability (heteroplasmy) of mtDNA in human brain by denaturing gradient gel electrophoresis and sequencing. Single base substitutions, as well as insertions or deletions of single bases, were numerous in the noncoding control region (D-loop), and 35-45% of the molecules from a single tissue showed sequence differences. By contrast, heteroplasmy in coding regions was not detected. The lower level of heteroplasmy in the coding regions is indicative of selection against deleterious mutations. Similar levels of heteroplasmy were found in two brain regions from the same individual, while no heteroplasmy was detected in blood. Thus, heteroplasmy seems to be more frequent in nonmitotic tissues. We observed a 7.7-fold increase in the frequency of deletions/insertions and a 2.2-fold increase in the overall frequency of heteroplasmic mutations in two individuals aged 96 and 99, relative to an individual aged 28. Our results show that intraindividual sequence variability occurs at a high frequency in the noncoding regions of normal human brain and indicate that small insertions and deletions might accumulate with age at a lower rate than large rearrangements.

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.

Muller's ratchet decreases fitness of a DNA-based microbe.

Muller proposed that an asexual organism will inevitably accumulate deleterious mutations, resulting in an increase of the mutational load and an inexorable, ratchet-like, loss of the least mutated class [Muller, H.J. (1964) Mutat. Res. 1, 2-9]. The operation of Muller's ratchet on real populations has been experimentally demonstrated only in RNA viruses. However, these cases are exceptional in that the mutation rates of the RNA viruses are extremely high. We have examined whether Muller's ratchet operates in Salmonella typhimurium, a DNA-based organism with a more typical genomic mutation rate. Cells were grown asexually under conditions expected to result in high genetic drift, and the increase in mutational load was determined. S. typhimurium accumulated mutations under these conditions such that after 1700 generations, 1% of the 444 lineages tested had suffered an obvious loss of fitness, as determined by decreased growth rate. These results suggest that in the absence of sex and with high genetic drift, genetic mechanisms, such as back or compensatory mutations, cannot compensate for the accumulation of deleterious mutations. In addition, we measured the appearance of auxotrophs, which allowed us to calculate an average spontaneous mutation rate of approximately 0.3-1.5 x 10(-9) mutations per base pair per generation. This rate is measured for the largest genetic target studied so far, a collection of about 200 genes.

Mutations in the gene encoding the alpha subunit of the rod cGMP-gated channel in autosomal recessive retinitis pigmentosa.

Mutations in the genes encoding two proteins of the retinal rod phototransduction cascade, opsin and the beta subunit of rod cGMP phosphodiesterase, cause retinitis pigmentosa (RP) in some families. Here we report defects in a third member of this biochemical pathway in still other patients with this disease. We screened 94 unrelated patients with autosomal dominant RP and 173 unrelated patients with autosomal recessive RP for mutations in the gene encoding the alpha subunit of the rod cGMP-gated cation channel. Five mutant sequences cosegregated with disease among four unrelated families with autosomal recessive RP. Two of these were nonsense mutations early in the reading frame (Glu76End and Lys139End) and one was a deletion encompassing most if not all of the transcriptional unit; these three alleles would not be expected to encode a functional channel. The remaining two mutations were a missense mutation (Ser316Phe) and a frameshift [Arg654(1-bp del)] mutation truncating the last 32 aa in the C terminus. The latter two mutations were expressed in vitro and found to encode proteins that were predominantly retained inside the cell instead of being targeted to the plasma membrane. We conclude that the absence or paucity of functional cGMP-gated cation channels in the plasma membrane is deleterious to rod photoreceptors and is an uncommon cause of RP.

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.

Activating mutations in the extracellular domain of the fibroblast growth factor receptor 2 function by disruption of the disulfide bond in the third immunoglobulin-like domain

Multiple human skeletal and craniosynostosis disorders, including Crouzon, Pfeiffer, Jackson–Weiss, and Apert syndromes, result from numerous point mutations in the extracellular region of fibroblast growth factor receptor 2 (FGFR2). Many of these mutations create a free cysteine residue that potentially leads to abnormal disulfide bond formation and receptor activation; however, for noncysteine mutations, the mechanism of receptor activation remains unclear. We examined the effect of two of these mutations, W290G and T341P, on receptor dimerization and activation. These mutations resulted in cellular transformation when expressed as FGFR2/Neu chimeric receptors. Additionally, in full-length FGFR2, the mutations induced receptor dimerization and elevated levels of tyrosine kinase activity. Interestingly, transformation by the chimeric receptors, dimerization, and enhanced kinase activity were all abolished if either the W290G or the T341P mutation was expressed in conjunction with mutations that eliminate the disulfide bond in the third immunoglobulin-like domain (Ig-3). These results demonstrate a requirement for the Ig-3 cysteine residues in the activation of FGFR2 by noncysteine mutations. Molecular modeling also reveals that noncysteine mutations may activate FGFR2 by altering the conformation of the Ig-3 domain near the disulfide bond, preventing the formation of an intramolecular bond. This allows the unbonded cysteine residues to participate in intermolecular disulfide bonding, resulting in constitutive activation of the receptor.

Recessive and dominant mutations in the ethylene biosynthetic gene ACS5 of Arabidopsis confer cytokinin insensitivity and ethylene overproduction, respectively

We identified a set of cytokinin-insensitive mutants by using a screen based on the ethylene-mediated triple response observed after treatment with low levels of cytokinins. One group of these mutants disrupts ACS5, a member of the Arabidopsis gene family that encodes 1-aminocyclopropane-1-carboxylate synthase, the first enzyme in ethylene biosynthesis. The ACS5 isoform is mainly responsible for the sustained rise in ethylene biosynthesis observed in response to low levels of cytokinin and appears to be regulated primarily by a posttranscriptional mechanism. Furthermore, the dominant ethylene-overproducing mutant eto2 was found to be the result of an alteration of the carboxy terminus of ACS5, suggesting that this domain acts as a negative regulator of ACS5 function.

A phenotype-based screen for embryonic lethal mutations in the mouse

The genetic pathways that control development of the early mammalian embryo have remained poorly understood, in part because the systematic mutant screens that have been so successful in the identification of genes and pathways that direct embryonic development in Drosophila, Caenorhabditis elegans, and zebrafish have not been applied to mammalian embryogenesis. Here we demonstrate that chemical mutagenesis with ethylnitrosourea can be combined with the resources of mouse genomics to identify new genes that are essential for mammalian embryogenesis. A pilot screen for abnormal morphological phenotypes of midgestation embryos identified five mutant lines; the phenotypes of four of the lines are caused by recessive traits that map to single regions of the genome. Three mutant lines display defects in neural tube closure: one is caused by an allele of the open brain (opb) locus, one defines a previously unknown locus, and one has a complex genetic basis. Two mutations produce novel early phenotypes and map to regions of the genome not previously implicated in embryonic patterning.