Genetic evidence for the involvement of DNA ligase IV in the DNA-PK-dependent pathway of non-homologous end joining in mammalian cells

Cells of vertebrates remove DNA double-strand breaks (DSBs) from their genome predominantly utilizing a fast, DNA-PKcs-dependent form of non-homologous end joining (D-NHEJ). Mutants with inactive DNA-PKcs remove the majority of DNA DSBs utilizing a slow, DNA-PKcs-independent pathway that does not utilize genes of the RAD52 epistasis group, is error-prone and can therefore be classified as a form of NHEJ (termed basic or B-NHEJ). We studied the role of DNA ligase IV in these pathways of NHEJ. Although biochemical studies show physical and functional interactions between the DNA-PKcs/Ku and the DNA ligase IV/Xrcc4 complexes suggesting operation within the same pathway, genetic evidence to support this notion is lacking in mammalian cells. Primary human fibroblasts (180BR) with an inactivating mutation in DNA ligase IV, rejoined DNA DSBs predominantly with slow kinetics similar to those observed in cells deficient in DNA-PKcs, or in wild-type cells treated with wortmannin to inactivate DNA-PK. Treatment of 180BR cells with wortmannin had only a small effect on DNA DSB rejoining and no effect on cell radiosensitivity to killing although it sensitized control cells to 180BR levels. This is consistent with DNA ligase IV functioning as a component of the D-NHEJ, and demonstrates the unperturbed operation of the DNA-PKcs-independent pathway (B-NHEJ) at significantly reduced levels of DNA ligase IV. In vitro, extracts of 180BR cells supported end joining of restriction endonuclease-digested plasmid to the same degree as extracts of control cells when tested at 10 mM Mg2+. At 0.5 mM Mg2+, where only DNA ligase IV is expected to retain activity, low levels of end joining (∼10% of 10 mM) were seen in the control but there was no detectable activity in 180BR cells. Antibodies raised against DNA ligase IV did not measurably inhibit end joining at 10 mM Mg2+ in either cell line. Thus, in contrast to the situation in vivo, end joining in vitro is dominated by pathways with properties similar to B-NHEJ that do not display a strong dependence on DNA ligase IV, with D-NHEJ retaining only a limited contribution. The implications of these observations to studies of NHEJ in vivo and in vitro are discussed.

Rapid, local adaptation of zooplankton behavior to changes in predation pressure in the absence of neutral genetic changes

Organisms producing resting stages provide unique opportunities for reconstructing the genetic history of natural populations. Diapausing seeds and eggs often are preserved in large numbers, representing entire populations captured in an evolutionary inert state for decades and even centuries. Starting from a natural resting egg bank of the waterflea Daphnia, we compare the evolutionary rates of change in an adaptive quantitative trait with those in selectively neutral DNA markers, thus effectively testing whether the observed genetic changes in the quantitative trait are driven by natural selection. The population studied experienced variable and well documented levels of fish predation over the past 30 years and shows correlated genetic changes in phototactic behavior, a predator-avoidance trait that is related to diel vertical migration. The changes mainly involve an increased plasticity response upon exposure to predator kairomone, the direction of the changes being in agreement with the hypothesis of adaptive evolution. Genetic differentiation through time was an order of magnitude higher for the studied behavioral trait than for neutral markers (DNA microsatellites), providing strong evidence that natural selection was the driving force behind the observed, rapid, evolutionary changes.

Evolution of codon usage bias in Drosophila

We first review what is known about patterns of codon usage bias in Drosophila and make the following points: (i) Drosophila genes are as biased or more biased than those in microorganisms. (ii) The level of bias of genes and even the particular pattern of codon bias can remain phylogenetically invariant for very long periods of evolution. (iii) However, some genes, even very tightly linked genes, can change very greatly in codon bias across species. (iv) Generally G and especially C are favored at synonymous sites in biased genes. (v) With the exception of aspartic acid, all amino acids contribute significantly and about equally to the codon usage bias of a gene. (vi) While most individual amino acids that can use G or C at synonymous sites display a preference for C, there are exceptions: valine and leucine, which prefer G. (vii) Finally, smaller genes tend to be more biased than longer genes. We then examine possible causes of these patterns and discount mutation bias on three bases: there is little evidence of regional mutation bias in Drosophila, mutation bias is likely toward A+T (the opposite of codon usage bias), and not all amino acids display the preference for the same nucleotide in the wobble position. Two lines of evidence support a selection hypothesis based on tRNA pools: highly biased genes tend to be highly and/or rapidly expressed, and the preferred codons in highly biased genes optimally bind the most abundant isoaccepting tRNAs. Finally, we examine the effect of bias on DNA evolution and confirm that genes with high codon usage bias have lower rates of synonymous substitution between species than do genes with low codon usage bias. Surprisingly, we find that genes with higher codon usage bias display higher levels of intraspecific synonymous polymorphism. This may be due to opposing effects of recombination.

Dynamic evolution of plant mitochondrial genomes: Mobile genes and introns and highly variable mutation rates

We summarize our recent studies showing that angiosperm mitochondrial (mt) genomes have experienced remarkably high rates of gene loss and concomitant transfer to the nucleus and of intron acquisition by horizontal transfer. Moreover, we find substantial lineage-specific variation in rates of these structural mutations and also point mutations. These findings mostly arise from a Southern blot survey of gene and intron distribution in 281 diverse angiosperms. These blots reveal numerous losses of mt ribosomal protein genes but, with one exception, only rare loss of respiratory genes. Some lineages of angiosperms have kept all of their mt ribosomal protein genes whereas others have lost most of them. These many losses appear to reflect remarkably high (and variable) rates of functional transfer of mt ribosomal protein genes to the nucleus in angiosperms. The recent transfer of cox2 to the nucleus in legumes provides both an example of interorganellar gene transfer in action and a starting point for discussion of the roles of mechanistic and selective forces in determining the distribution of genetic labor between organellar and nuclear genomes. Plant mt genomes also acquire sequences by horizontal transfer. A striking example of this is a homing group I intron in the mt cox1 gene. This extraordinarily invasive mobile element has probably been acquired over 1,000 times separately during angiosperm evolution via a recent wave of cross-species horizontal transfers. Finally, whereas all previously examined angiosperm mtDNAs have low rates of synonymous substitutions, mtDNAs of two distantly related angiosperms have highly accelerated substitution rates.

The role of genetic and genomic attributes in the success of polyploids

In 1950, G. Ledyard Stebbins devoted two chapters of his book Variation and Evolution in Plants (Columbia Univ. Press, New York) to polyploidy, one on occurrence and nature and one on distribution and significance. Fifty years later, many of the questions Stebbins posed have not been answered, and many new questions have arisen. In this paper, we review some of the genetic attributes of polyploids that have been suggested to account for the tremendous success of polyploid plants. Based on a limited number of studies, we conclude: (i) Polyploids, both individuals and populations, generally maintain higher levels of heterozygosity than do their diploid progenitors. (ii) Polyploids exhibit less inbreeding depression than do their diploid parents and can therefore tolerate higher levels of selfing; polyploid ferns indeed have higher levels of selfing than do their diploid parents, but polyploid angiosperms do not differ in outcrossing rates from their diploid parents. (iii) Most polyploid species are polyphyletic, having formed recurrently from genetically different diploid parents. This mode of formation incorporates genetic diversity from multiple progenitor populations into the polyploid “species”; thus, genetic diversity in polyploid species is much higher than expected by models of polyploid formation involving a single origin. (iv) Genome rearrangement may be a common attribute of polyploids, based on evidence from genome in situ hybridization (GISH), restriction fragment length polymorphism (RFLP) analysis, and chromosome mapping. (v) Several groups of plants may be ancient polyploids, with large regions of homologous DNA. These duplicated genes and genomes can undergo divergent evolution and evolve new functions. These genetic and genomic attributes of polyploids may have both biochemical and ecological benefits that contribute to the success of polyploids in nature.

Rates of nucleotide substitution in sexual and anciently asexual rotifers

The class Bdelloidea of the phylum Rotifera is the largest well studied eukaryotic taxon in which males and meiosis are unknown, and the only one for which these indications of ancient asexuality are supported by cytological and molecular genetic evidence. We estimated the rates of synonymous and nonsynonymous substitutions in the hsp82 heat shock gene in bdelloids and in facultatively sexual rotifers of the class Monogononta, employing distance based and maximum likelihood methods. Relative-rate tests, using acanthocephalan rotifers as an outgroup, showed slightly higher rates of nonsynonymous substitution and slightly lower rates of synonymous substitution in bdelloids as compared with monogononts. The opposite trend, however, was seen in intraclass pairwise comparisons. If, as it seems, bdelloids have evolved asexually, an equality of bdelloid and monogonont substitution rates would suggest that the maintenance of sexual reproduction in monogononts is not attributable to an effect of sexual reproduction in limiting the load of deleterious nucleotide substitutions.

RNA virus error catastrophe: Direct molecular test by using ribavirin

RNA viruses evolve rapidly. One source of this ability to rapidly change is the apparently high mutation frequency in RNA virus populations. A high mutation frequency is a central tenet of the quasispecies theory. A corollary of the quasispecies theory postulates that, given their high mutation frequency, animal RNA viruses may be susceptible to error catastrophe, where they undergo a sharp drop in viability after a modest increase in mutation frequency. We recently showed that the important broad-spectrum antiviral drug ribavirin (currently used to treat hepatitis C virus infections, among others) is an RNA virus mutagen, and we proposed that ribavirin's antiviral effect is by forcing RNA viruses into error catastrophe. However, a direct demonstration of error catastrophe has not been made for ribavirin or any RNA virus mutagen. Here we describe a direct demonstration of error catastrophe by using ribavirin as the mutagen and poliovirus as a model RNA virus. We demonstrate that ribavirin's antiviral activity is exerted directly through lethal mutagenesis of the viral genetic material. A 99.3% loss in viral genome infectivity is observed after a single round of virus infection in ribavirin concentrations sufficient to cause a 9.7-fold increase in mutagenesis. Compiling data on both the mutation levels and the specific infectivities of poliovirus genomes produced in the presence of ribavirin, we have constructed a graph of error catastrophe showing that normal poliovirus indeed exists at the edge of viability. These data suggest that RNA virus mutagens may represent a promising new class of antiviral drugs.

Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius

Microbes whose genomes are encoded by DNA and for which adequate information is available display similar genomic mutation rates (average 0.0034 mutations per chromosome replication, range 0.0025 to 0.0046). However, this value currently is based on only a few well characterized microbes reproducing within a narrow range of environmental conditions. In particular, no genomic mutation rate has been determined either for a microbe whose natural growth conditions may extensively damage DNA or for any member of the archaea, a prokaryotic lineage deeply diverged from both bacteria and eukaryotes. Both of these conditions are met by the extreme thermoacidophile Sulfolobus acidocaldarius. We determined the genomic mutation rate for this species when growing at pH 3.5 and 75°C based on the rate of forward mutation at the pyrE gene and the nucleotide changes identified in 101 independent mutants. The observed value of about 0.0018 extends the range of DNA-based microbes with rates close to the standard rate simultaneously to an archaeon and to an extremophile whose cytoplasmic pH and normal growth temperature greatly accelerate the spontaneous decomposition of DNA. The mutations include base pair substitutions (BPSs) and additions and deletions of various sizes, but the S. acidocaldarius spectrum differs from those of other DNA-based organisms in being relatively poor in BPSs. The paucity of BPSs cannot yet be explained by known properties of DNA replication or repair enzymes of Sulfolobus spp. It suggests, however, that molecular evolution per genome replication may proceed more slowly in S. acidocaldarius than in other DNA-based organisms examined to date.

Expression of error-prone polymerases in BL2 cells activated for Ig somatic hypermutation

High affinity antibodies are generated in mice and humans by means of somatic hypermutation (SHM) of variable (V) regions of Ig genes. Mutations with rates of 10−5–10−3 per base pair per generation, about 106-fold above normal, are targeted primarily at V-region hot spots by unknown mechanisms. We have measured mRNA expression of DNA polymerases ι, η, and ζ by using cultured Burkitt's lymphoma (BL)2 cells. These cells exhibit 5–10-fold increases in heavy-chain V-region mutations targeted only predominantly to RGYW (R = A or G, Y = C or T, W = T or A) hot spots if costimulated with T cells and IgM crosslinking, the presumed in vivo requirements for SHM. An ∼4-fold increase pol ι mRNA occurs within 12 h when cocultured with T cells and surface IgM crosslinking. Induction of pols η and ζ occur with T cells, IgM crosslinking, or both stimuli. The fidelity of pol ι was measured at RGYW hot- and non-hot-spot sequences situated at nicks, gaps, and double-strand breaks. Pol ι formed T⋅G mispairs at a frequency of 10−2, consistent with SHM-generated C to T transitions, with a 3-fold increased error rate in hot- vs. non-hot-spot sequences for the single-nucleotide overhang. The T cell and IgM crosslinking-dependent induction of pol ι at 12 h may indicate an SHM “triggering” event has occurred. However, pols ι, η, and ζ are present under all conditions, suggesting that their presence is not sufficient to generate mutations because both T cell and IgM stimuli are required for SHM induction.

Common endocrine and genetic mechanisms of behavioral development in male and worker honey bees and the evolution of division of labor.

Temporal polyethism is a highly derived form of behavioral development displayed by social insects. Hormonal and genetic mechanisms regulating temporal polyethism in worker honey bees have been identified, but the evolution of these mechanisms is not well understood. We performed three experiments with male honey bees (drones) to investigate how mechanisms regulating temporal polyethism may have evolved because, relative to workers, drones display an intriguing combination of similarities and differences in behavioral development. We report that behavioral development in drones is regulated by mechanisms common to workers. In experiment 1, drones treated with the juvenile hormone (JH) analog methoprene started flying at significantly younger ages than did control drones, as is the case for workers. In experiment 2, there was an age-related increase in JH associated with the onset of drone flight, as in workers. In experiment 3, drones derived from workers with fast rates of behavioral development themselves started flying at younger ages than drones derived from workers with slower rates of behavioral development. These results suggest that endocrine and genetic mechanisms associated with temporal polyethism did not evolve strictly within the context of worker social behavior.

Age-specific inbreeding depression and components of genetic variance in relation to the evolution of senescence.

Two major theories of the evolution of senescence (mutation accumulation and antagonistic pleiotropy) make different predictions about the relationships between age, inbreeding effects, and the magnitude of genetic variance components of life-history components. We show that, under mutation accumulation, inbreeding decline and three major components of genetic variance are expected to increase with age in randomly mating populations. Under the simplest version of the antagonistic pleiotropy model, no changes in the severity of inbreeding decline, dominance variance, or the genetic variance of chromosomal homozygotes are expected, but additive genetic variance may increase with age. Age-specific survival rates and mating success were measured on virgin males, using lines extracted from a population of Drosophila melanogaster. For both traits, inbreeding decline and several components of genetic variance increase with age. The results are consistent with the mutation accumulation model, but can only be explained by antagonistic pleiotropy if there is a general tendency for an increase with age in the size of allelic effects on these life-history traits.

Positive selection and rates of evolution in immunodeficiency viruses from humans and chimpanzees.

Evolutionary theory predicts the recent spread of primate immunodeficiency viruses (PIVs) to new human populations to be accompanied by positive selection in response to new host environments and/or by random genetic drift. I assess evidence for positive selection in human and chimpanzee PIVs type I (PIV1s), using ratios of synonymous to nonsynonymous nucleotide change based on branch lengths and outgroup rooting. Ratios are smaller for PIV1s from humans than for PIV1 from a chimpanzee for the pol, gag, and env glycoprotein 120 (gp120) regions, indicating greater effects of positive selection in PIV1s from humans. Parsimony-based relative rate tests for amino acid changes showed significant differences between PIV1s from humans and chimpanzees in 18 of 48 pairwise comparisons, with all 18 showing faster rates of change in PIV1s from humans. This study indicates that in some instances, the recent evolution of human PIV1s follows a speciational pattern, in which increased diversification of taxa is correlated with greater amounts of character change appearing and being maintained through time. This extends the generality of the speciational pattern to a group of organisms (viruses) having the fastest known rates of anagenetic change for nucleotide characters and indicates that comprehensive understanding of PIV1 evolution requires consideration of both anagenetic change within viral lineages and the relative historical success of different viral clades. Phylogenetic analyses show that neither PIV1s infecting humans nor those infecting chimpanzees represent monophyletic groups and suggest multiple host-species shifts for PIV1s.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Analysis of RNA-binding proteins by in vitro genetic selection: identification of an amino acid residue important for locking U1A onto its RNA target.

An in vitro genetic system was developed as a rapid means for studying the specificity determinants of RNA-binding proteins. This system was used to investigate the origin of the RNA-binding specificity of the mammalian spliceosomal protein U1A. The U1A domain responsible for binding to U1 small nuclear RNA was locally mutagenized and displayed as a combinatorial library on filamentous bacteriophage. Affinity selection identified four U1A residues in the mutagenized region that are important for specific binding to U1 hairpin II. One of these residues (Leu-49) disproportionately affects the rates of binding and release and appears to play a critical role in locking the protein onto the RNA. Interestingly, a protein variant that binds more tightly than U1A emerged during the selection, showing that the affinity of U1A for U1 RNA has not been optimized during evolution.

Genetic and comparative analyses reveal an alternative secondary structure in the region of nt 912 of Escherichia coli 16S rRNA.

Mutations at position 912 of Escherichia coli 16S rRNA result in two notable phenotypes. The C-->U transition confers resistance to streptomycin, a translational-error-inducing antibiotic, while a C-->G transversion causes marked retardation of cell growth rate. Starting with the slow-growing G912 mutant, random mutagenesis was used to isolate a second site mutation that restored growth nearly to the wild-type rate. The second site mutation was identified as a G-->C transversion at position 885 in 16S rRNA. Cells containing the G912 mutation had an increased doubling time, abnormal sucrose gradient ribosome/subunit profile, increased sensitivity to spectinomycin, dependence upon streptomycin for growth in the presence of spectinomycin, and slower translation rate, whereas cells with the G912/C885 double mutation were similar to wild type in these assays. Comparative analysis showed there was significant covariation between positions 912 and 885. Thus the second-site suppressor analysis, the functional assays, and the comparative data suggest that the interaction between nt 912 and nt 885 is conserved and necessary for normal ribosome function. Furthermore, the comparative data suggest that the interaction extends to include G885-G886-G887 pairing with C912-U911-C910. An alternative secondary structure element for the central domain of 16S rRNA is proposed.