Population Genetics of Polymorphism and Divergence

Frequencies of mutant sites are modeled as a Poisson random field in two species that share a sufficiently recent common ancestor. The selective effect of the new alleles can be favorable, neutral, or detrimental. The model is applied to the sample configurations of nucleotides in the alcohol dehydrogenase gene (Adh) in Drosophila simulans and Drosophila yakuba. Assuming a synonymous mutation rate of 1.5 x 10(-8) per site per year and 10 generations per year, we obtain estimates for the effective population size (N(e) = 6.5 x 10(6)), the species divergence time (tdiv = 3.74 million years), and an average selection coefficient (sigma = 1.53 x 10(-6) per generation for advantageous or mildly detrimental replacements), although it is conceivable that only two of the amino acid replacements were selected and the rest neutral. The analysis, which includes a sampling theory for the independent infinite sites model with selection, also suggests the estimate that the number of amino acids in the enzyme that are susceptible to favorable mutation is in the range 2-23 at any one time. The approach provides a theoretical basis for the use of a 2 x 2 contingency table to compare fixed differences and polymorphic sites with silent sites and amino acid replacements.     

Inferring Weak Selection from Patterns of Polymorphism and Divergence at ``silent' Sites in Drosophila DNA

Patterns of codon usage and ``silent'''' DNA divergence suggest that natural selection discriminates among synonymous codons in Drosophila. ``Preferred'''' codons are consistently found in higher frequencies within their synonymous families in Drosophila melanogaster genes. This suggests a simple model of silent DNA evolution where natural selection favors mutations from unpreferred to preferred codons (preferred changes). Changes in the opposite direction, from preferred to unpreferred synonymous codons (unpreferred changes), are selected against. Here, selection on synonymous DNA mutations is investigated by comparing the evolutionary dynamics of these two categories of silent DNA changes. Sequences from outgroups are used to determine the direction of synonymous DNA changes within and between D. melanogaster and Drosophila simulans for five genes. Population genetics theory shows that differences in the fitness effect of mutations can be inferred from the comparison of ratios of polymorphism to divergence. Unpreferred changes show a significantly higher ratio of polymorphism to divergence than preferred changes in the D. simulans lineage, confirming the action of selection at silent sites. An excess of unpreferred fixations in 28 genes suggests a relaxation of selection on synonymous mutations in D. melanogaster. Estimates of selection coefficients for synonymous mutations (3.6 <|N(e)s| < 1.3) in D. simulans are consistent with the reduced efficacy of natural selection (|N(e)s| < 1) in the three- to sixfold smaller effective population size of D. melanogaster. Synonymous DNA changes appear to be a prevalent class of weakly selected mutations in Drosophila.     

Sequence-Dependent Bending of DNA and Phasing of Nucleosomes

Abstract

Conformational analysis has revealed anisotropic flexibility of the B-DNA double helix: it bends most easily into the grooves, being the most rigid when bent in a perpendicular direction. This result implies that DNA in a nucleosome is curved by means of relatively sharp bends (“mini-kinks”) which are directed into the major and minor grooves alternatively and separated by 5–6 base pairs. The “mini-kink” model proved to be in keeping with the x-ray structure of the B-DNA dodecamer resolved later, which exhibits two “annealed kinks”, also directed into the grooves.

The anisotropy of B DNA is sequence-dependent: the pyrimidine-purine dimers (YR) favor bending into the minor groove, and the purine-pyrimidine dinucleotides (RY), into the minor one. The RR and YY dimers appear to be the most rigid dinucleotides. Thus, a DNA fragment consisting of the interchanging oligopurine and oligopyrmidine blocks 5–6 base pairs long should manifest a spectacular curvature in solution.

Similarly, a nucleotide sequence containing the RY and YR dimers separated by a half-pitch of the double helix is the most suitable for wrapping around the nucleosomal core. Analysis of the numerous examples demonstrating the specific alignment of nucleosomes on DNA confirms this concept. So, the sequence-dependent “mechanical” properties of the double helix influence the spatial arrangement of DNA in chromatin.     

Polymorphism and Divergence at a Drosophila Pseudogene Locus

The larval cuticle protein (Lcp) cluster in Drosophila melanogaster contains four functional genes and a closely related pseudogene. A 630-bp fragment including the larval cuticle pseudogene locus (Lcpψ) was nucleotide sequenced in 10 strains of D. melanogaster and a 458-bp Lcpψ fragment from D. simulans was also sequenced. We used these data to test the hypotheses that the rates of synonymous and nonsynonymous substitution are equal, that the absolute levels of variation are higher than in functional genes, and that intraspecific polymorphism is correlated with interspecific divergence. As predicted, synonymous and nonsynonymous substitution rates were equivalent, and overall nucleotide divergence between D. melanogaster and D. simulans (Jukes-Cantor distance = 0.149 +/- 0.150) was extremely high. However, within-species DNA sequence comparisons at Lcpψ revealed lower levels of polymorphism ( & = 0.001 +/- 0.001) than at many functional loci in D. melanogaster. Using the HUDSON, KREITMAN, and AGUADE (HKA) test, we show that the level of polymorphism in Lcpψ within D. melanogaster is lower than expected given the amount of divergence between D. melanogaster and D. simulans when the pseudogene data are compared to the Adh 5' flanking region. Because the Lcpψ lies in a region of relatively infrequent recombination, we suggest that the low level of within-species polymorphism is the result of background selection.     

Nucleosomes protect DNA from DNA methylation in vivo and in vitro

Positioned nucleosomes limit the access of proteins to DNA. However, the impact of nucleosomes on DNA methylation in vitro and in vivo is poorly understood. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the de novo methyltransferases. We show that compared to linker DNA, nucleosomal DNA is largely devoid of CpG methylation. ATP-dependent chromatin remodelling frees nucleosomal CpG dinucleotides and renders the remodelled nucleosome a 2-fold better substrate for Dnmt3a methyltransferase compared to free DNA. These results reflect the situation in vivo, as quantification of nucleosomal DNA methylation levels in HeLa cells shows a 2-fold decrease of nucleosomal DNA methylation levels compared to linker DNA. Our findings suggest that nucleosomal positions are stably maintained in vivo and nucleosomal occupancy is a major determinant of global DNA methylation patterns in vivo.     

Placeholder Nucleosomes Underlie Germline-to-Embryo DNA Methylation Reprogramming

1. Download high-res image (222KB)
2. Download full-size image

     

Initial stages of DNA Base Excision Repair in Nucleosomes

Molecular Biology - In mammalian cells, base excision repair (BER) is the main pathway responsible for the correction of a variety of chemically modified DNA bases. DNA packaging in chromatin...     

Repeatability of Habitat-Associated Divergence in Shell Shape of Turtles

Repeated patterns of phenotypic divergence between environments across disparate taxa provide strong evidence for the generation of adaptive phenotypes. Flow velocity is an important selective force in aquatic habitats; however, among vertebrates, the study of its effects on morphology has been limited almost exclusively to fully-aquatic bony fishes. We tested whether three confamilial species of semi-aquatic freshwater turtle (family Emydidae: Graptemys pseudogeographica, Graptemys nigrinoda, and Pseudemys concinna) displayed similar patterns of phenotypic divergence in carapace shape between fast- and slow-flowing aquatic environments. We used (1) geometric morphometrics to quantify shell shape, (2) multivariate analysis of variance to test the effects of species, sex, and flow, and (3) phenotypic trajectory analysis to compare patterns of divergence for six species-sex groups. We found significant effects on shell shape for all factors. In general, ecomorphs from fast-flowing habitats had flatter shells than those from slow-flowing habitats. Furthermore, results of trajectory analysis indicate that the degree to which, as well as the way in which, ecomorphs differed were concordant across all species. Our findings demonstrate that the effects of flow are not limited to fully-aquatic vertebrates, and provide evidence of the ability of flow to drive repeatable phenotypic divergence in tetrapods.     

Divergence and Polymorphism Under the Nearly Neutral Theory of Molecular Evolution

The nearly neutral theory attributes most nucleotide substitution and polymorphism to genetic drift acting on weakly selected mutants, and assumes that the selection coefficients for these mutants are drawn from a continuous distribution. This means that parameter estimation can require numerical integration, and this can be computationally costly and inaccurate. Furthermore, the leading parameter dependencies of important quantities can be unclear, making results difficult to understand. For some commonly used distributions of mutant effects, we show how these problems can be avoided by writing equations in terms of special functions. Series expansion then allows for their rapid calculation and, also, illuminates leading parameter dependencies. For example, we show that if mutants are gamma distributed, the neutrality index is largely independent of the effective population size. However, we also show that such results are not robust to misspecification of the functional form of distribution. Some implications of these findings are then discussed.     

Nucleosomes positioned by ORC facilitate the initiation of DNA replication

The packaging of eukaryotic DNA into nucleosomes is a critical regulator of nuclear events. To address the interplay between chromatin and replication initiation, we have assessed the determinants and function of the nucleosomal configuration of S. cerevisiae replication origins. Using in vitro and in vivo assays, we demonstrate that the yeast initiator, the origin recognition complex (ORC), is required to maintain the nucleosomal configuration adjacent to origins. Disruption of the ORC-directed nucleosomal arrangement at an origin interferes with initiation of replication, but does not alter the association of ORC with the origin. Instead, the nucleosomes positioned by ORC are important for prereplicative complex formation. These findings suggest that origin-proximal nucleosomes facilitate replication initiation, and that local chromatin structure affects origin function.     

Herpes Simplex Virus 1 DNA Is in Unstable Nucleosomes throughout the Lytic Infection Cycle,and the Instability of the Nucleosomes Is Independent of DNA Replication

Herpes simplex virus 1 (HSV-1) DNA is chromatinized during latency and consequently regularly digested by micrococcal nuclease (MCN) to nucleosome-size fragments. In contrast, MCN digests HSV-1 DNA in lytically infected cells to mostly heterogeneous sizes. Yet HSV-1 DNA coimmunoprecipitates with histones during lytic infections. We have shown that at 5 h postinfection, most nuclear HSV-1 DNA is in particularly unstable nucleoprotein complexes and consequently is more accessible to MCN than DNA in cellular chromatin. HSV-1 DNA was quantitatively recovered at this time in complexes with the biophysical properties of mono- to polynucleosomes following a modified MCN digestion developed to detect potential unstable intermediates. We proposed that most HSV-1 DNA is in unstable nucleosome-like complexes during lytic infections. Physiologically, nucleosome assembly typically associates with DNA replication, although DNA replication transiently disrupts nucleosomes. It therefore remained unclear whether the instability of the HSV-1 nucleoprotein complexes was related to the ongoing viral DNA replication. Here we tested whether HSV-1 DNA is in unstable nucleosome-like complexes before, during, or after the peak of viral DNA replication or when HSV-1 DNA replication is inhibited. HSV-1 DNA was quantitatively recovered in complexes fractionating as mono- to polynucleosomes from nuclei harvested at 2, 5, 7, or 9 h after infection, even if viral DNA replication was inhibited. Therefore, most HSV-1 DNA is in unstable nucleosome-like complexes throughout the lytic replication cycle, and the instability of these complexes is surprisingly independent of HSV-1 DNA replication. The specific accessibility of nuclear HSV-1 DNA, however, varied at different times after infection.     

Selective Anchoring of DNA Methyltransferases 3A and 3B to Nucleosomes Containing Methylated DNA

Proper DNA methylation patterns are essential for mammalian development and differentiation. DNA methyltransferases (DNMTs) primarily establish and maintain global DNA methylation patterns; however, the molecular mechanisms for the generation and inheritance of methylation patterns are still poorly understood. We used sucrose density gradients of nucleosomes prepared by partial and maximum micrococcal nuclease digestion, coupled with Western blot analysis to probe for the interactions between DNMTs and native nucleosomes. This method allows for analysis of the in vivo interactions between the chromatin modification enzymes and their actual nucleosomal substrates in the native state. We show that little free DNA methyltransferase 3A and 3B (DNMT3A/3B) exist in the nucleus and that almost all of the cellular contents of DNMT3A/3B, but not DNMT1, are strongly anchored to a subset of nucleosomes. This binding of DNMT3A/3B does not require the presence of other well-known chromatin-modifying enzymes or proteins, such as proliferating cell nuclear antigen, heterochromatin protein 1, methyl-CpG binding protein 2, Enhancer of Zeste homolog 2, histone deacetylase 1, and UHRF1, but it does require an intact nucleosomal structure. We also show that nucleosomes containing methylated SINE and LINE elements and CpG islands are the main sites of DNMT3A/3B binding. These data suggest that inheritance of DNA methylation requires cues from the chromatin component in addition to hemimethylation.Proper DNA methylation patterns are essential for mammalian development and differentiation. More than three decades ago, de novo cytosine DNA methylation and its maintenance were proposed to exist in eukaryotic cells (29, 54); however, the molecular mechanisms for the generation and inheritance of methylation patterns are still poorly understood. DNA methyltransferases (DNMTs) DNMT1, DNMT3A, and DNMT3B primarily establish and maintain global DNA methylation patterns (39, 48). DNMT1 preferentially methylates hemimethylated DNA in vitro (7) and is tethered to replication foci during S phase (38). In contrast, DNMT3A and DNMT3B (DNMT3A/3B) have no preference for hemimethylated DNA (49) and are required for de novo methylation of genomic DNA (48). It has been thought that DNMT1 acts mainly as a “maintenance methyltransferase” during DNA synthesis and that DNMT3A and DNMT3B act as “de novo” enzymes. However, more recent studies indicate that DNMT1 may also be required for de novo methylation of genomic DNA (17, 30) and that DNMT3A/3B are also required for maintenance functions (11, 40, 55). Furthermore, the different DNMTs cooperate in maintaining the methylation of some regions of the genome, particularly repetitive elements (40, 53).Recruitment of individual DNMT enzymes to different regions of chromatin in vivo, particularly to gene regulatory regions, may require interaction with auxiliary factors (28, 36). DNMT1, which is diffusely localized throughout nuclei in non-S-phase cells (38), is targeted to replication foci by interacting with proliferating cell nuclear antigen (PCNA) (15) and also physically interacts with UHRF1 (ubiquitinlike, containing PHD and RING finger domains 1) that binds to hemimethylated DNA (3, 4, 8, 27, 62). DNMT3 enzymes are usually found localized to heterochromatin regions in most transient-expression assays (5, 12). As genomic DNA in chromatin is packaged into nucleosomes which might limit the accessibility of target sites to the enzymes, the interaction of DNMTs with nucleosomes in a chromatin context is important for the regulation of genomic methylation.Genetic and biochemical studies have provided many insights into the distinct and cooperative functions of the DNMT enzymes; however, few of these studies have addressed how they interact with chromatin in vivo. Recombinant DNMT1 and DNMT3 enzymes can methylate the CpG sites on nucleosomes assembled in vitro (26, 50, 56, 65). Recently DNMT3L has been found to connect DNMT3A2 to nucleosomes in embryonic stem cells (52). However, DNMT3L is expressed only during gametogenesis and embryonic stages (1, 9), suggesting that other mechanisms might be necessary for directing the enzyme to specific chromatin regions in somatic cells.In the present study, we investigated how different DNMT enzymes interact with chromatin at the nucleosomal level in somatic cell lines. Micrococcal nuclease (MNase) treatment of nuclei in a low-ionic-strength buffer digests nucleosomal linker DNA regions, thereby minimizing the disruption of protein complexes on the nucleosomes. We prepared nucleosomes from partial or maximum MNase-digested nuclei and resolved them on sucrose density gradients to analyze their interactions with chromatin proteins. The results indicate that while DNMT1 interacts primarily with linker DNA, DNMT3A/3B enzymes interact strongly with nucleosomes containing methylated repetitive elements and also containing methylated CpG islands (CGIs) and may not require additional proteins for this strong binding. These data are particularly intriguing in that they provide insights into the mechanisms of the interaction of DNMTs with chromatin and maintenance of DNA methylation in somatic cells.     

Characterization of Dnmt1 Binding and DNA Methylation on Nucleosomes and Nucleosomal Arrays

The packaging of DNA into nucleosomes and the organisation into higher order structures of chromatin limits the access of sequence specific DNA binding factors to DNA. In cells, DNA methylation is preferentially occuring in the linker region of nucleosomes, suggesting a structural impact of chromatin on DNA methylation. These observations raise the question whether DNA methyltransferases are capable to recognize the nucleosomal substrates and to modify the packaged DNA. Here, we performed a detailed analysis of nucleosome binding and nucleosomal DNA methylation by the maintenance DNA methyltransferase Dnmt1. Our binding studies show that Dnmt1 has a DNA length sensing activity, binding cooperatively to DNA, and requiring a minimal DNA length of 20 bp. Dnmt1 needs linker DNA to bind to nucleosomes and most efficiently recognizes nucleosomes with symmetric DNA linkers. Footprinting experiments reveal that Dnmt1 binds to both DNA linkers exiting the nucleosome core. The binding pattern correlates with the efficient methylation of DNA linkers. However, the enzyme lacks the ability to methylate nucleosomal CpG sites on mononucleosomes and nucleosomal arrays, unless chromatin remodeling enzymes create a dynamic chromatin state. In addition, our results show that Dnmt1 functionally interacts with specific chromatin remodeling enzymes to enable complete methylation of hemi-methylated DNA in chromatin.     

Polymorphism and Divergence of Multilocus DNA Markers in Sibling Species Chironomus riparius Meigen and Chironomus piger Strenzke (Diptera,Chironomidae)

Intra- and interspecific variation and divergence of multilocus markers for genomic DNA of the sibling species from the thimmi group,Chironomus riparius and C. piger, were studied by PCR with arbitrary primers (RAPD). A high level of RAPD polymorphism was determined in both laboratory and natural populations of these species. The genetic distances were estimated between the C. riparius populations and between the sibling species C. riparius and C. piger. The genetic distance between C. riparius andC. piger was 4 to 5 times higher than that between the C. riparius populations. A comparison of the variation and divergence for the RAPD markers with those for other genomic markers—enzyme-coding genes and chromosomes (gene linkage groups)—showed that different components of the genome differed in their contribution to the genome divergence.     

Evolutionary Divergence and Convergence in Shape and Size Within African Antelope Proximal Phalanges

Morphological convergence amongst species inhabiting similar environments but having different evolutionary histories is a concept central to evolutionary biology. Cases of divergent evolution, where there is morphological divergence between closely related species exploiting different environments, are less well studied. Here we show divergent evolution in the morphology of the proximal phalanges of several closely related African antelope species inhabiting different environments. This morphological divergence was consistently observed in both a neutral morphospace and an externally ordinated morphospace. Divergence, but not convergence, was also observed when size and shape were considered independently. Finally, convergent evolution of the morphology of the proximal phalanges was observed, but only in the externally ordinated morphospace. Size shows less correlation with phylogeny than does shape. Therefore, we suggest that divergence in size will occur more readily when a species encounters new environmental conditions than divergence in shape. These findings are compatible with observations of rapid dwarfing on islands (Foster’s rule).     

Polymorphism of Natural DNA

X-ray fibre diffraction shows that natural very AT-rich DNA exists in at least four novel and distinct conformations depending on the sequence and salt content.     

Nucleosomes containing methylated DNA stabilize DNA methyltransferases 3A/3B and ensure faithful epigenetic inheritance

How epigenetic information is propagated during somatic cell divisions is still unclear but is absolutely critical for preserving gene expression patterns and cellular identity. Here we show an unanticipated mechanism for inheritance of DNA methylation patterns where the epigenetic mark not only recruits the catalyzing enzyme but also regulates the protein level, i.e. the enzymatic product (5-methylcytosine) determines the level of the methylase, thus forming a novel homeostatic inheritance system. Nucleosomes containing methylated DNA stabilize de novo DNA methyltransferases, DNMT3A/3B, allowing little free DNMT3A/3B enzymes to exist in the nucleus. Stabilization of DNMT3A/3B on nucleosomes in methylated regions further promotes propagation of DNA methylation. However, reduction of cellular DNA methylation levels creating more potential CpG substrates counter-intuitively results in a dramatic decrease of DNMT3A/3B proteins due to diminished nucleosome binding and subsequent degradation of the unstable free proteins. These data show an unexpected self-regulatory inheritance mechanism that not only ensures somatic propagation of methylated states by DNMT1 and DNMT3A/3B enzymes but also prevents aberrant de novo methylation by causing degradation of free DNMT3A/3B enzymes.