Comparative DNA sequence and methylation analyses of orthologous genes in humans and non-human primates: Genetic and epigenetic footprints of evolution

Welche genetische Unterschiede machen uns verschieden von unseren nächsten Verwandten, den Schimpansen, und andererseits so ähnlich zu den Schimpansen? Was wir untersuchen und auch verstehen wollen, ist die komplexe Beziehung zwischen den multiplen genetischen und epigenetischen Unterschieden, deren Interaktion mit diversen Umwelt- und Kulturfaktoren in den beobachteten phänotypischen Unterschieden resultieren. Um aufzuklären, ob chromosomale Rearrangements zur Divergenz zwischen Mensch und Schimpanse beigetragen haben und welche selektiven Kräfte ihre Evolution geprägt haben, habe ich die kodierenden Sequenzen von 2 Mb umfassenden, die perizentrischen Inversionsbruchpunkte flankierenden Regionen auf den Chromosomen 1, 4, 5, 9, 12, 17 und 18 untersucht. Als Kontrolle dienten dabei 4 Mb umfassende kollineare Regionen auf den rearrangierten Chromosomen, welche mindestens 10 Mb von den Bruchpunktregionen entfernt lagen. Dabei konnte ich in den Bruchpunkten flankierenden Regionen im Vergleich zu den Kontrollregionen keine höhere Proteinevolutionsrate feststellen. Meine Ergebnisse unterstützen nicht die chromosomale Speziationshypothese für Mensch und Schimpanse, da der Anteil der positiv selektierten Gene (5,1% in den Bruchpunkten flankierenden Regionen und 7% in den Kontrollregionen) in beiden Regionen ähnlich war. Durch den Vergleich der Anzahl der positiv und negativ selektierten Gene per Chromosom konnte ich feststellen, dass Chromosom 9 die meisten und Chromosom 5 die wenigsten positiv selektierten Gene in den Bruchpunkt flankierenden Regionen und Kontrollregionen enthalten. Die Anzahl der negativ selektierten Gene (68) war dabei viel höher als die Anzahl der positiv selektierten Gene (17). Eine bioinformatische Analyse von publizierten Microarray-Expressionsdaten (Affymetrix Chip U95 und U133v2) ergab 31 Gene, die zwischen Mensch und Schimpanse differentiell exprimiert sind. Durch Untersuchung des dN/dS-Verhältnisses dieser 31 Gene konnte ich 7 Gene als negativ selektiert und nur 1 Gen als positiv selektiert identifizieren. Dieser Befund steht im Einklang mit dem Konzept, dass Genexpressionslevel unter stabilisierender Selektion evolvieren. Die meisten positiv selektierten Gene spielen überdies eine Rolle bei der Fortpflanzung. Viele dieser Speziesunterschiede resultieren eher aus Änderungen in der Genregulation als aus strukturellen Änderungen der Genprodukte. Man nimmt an, dass die meisten Unterschiede in der Genregulation sich auf transkriptioneller Ebene manifestieren. Im Rahmen dieser Arbeit wurden die Unterschiede in der DNA-Methylierung zwischen Mensch und Schimpanse untersucht. Dazu wurden die Methylierungsmuster der Promotor-CpG-Inseln von 12 Genen im Cortex von Menschen und Schimpansen mittels klassischer Bisulfit-Sequenzierung und Bisulfit-Pyrosequenzierung analysiert. Die Kandidatengene wurden wegen ihrer differentiellen Expressionsmuster zwischen Mensch und Schimpanse sowie wegen Ihrer Assoziation mit menschlichen Krankheiten oder dem genomischen Imprinting ausgewählt. Mit Ausnahme einiger individueller Positionen zeigte die Mehrzahl der analysierten Gene keine hohe intra- oder interspezifische Variation der DNA-Methylierung zwischen den beiden Spezies. Nur bei einem Gen, CCRK, waren deutliche intraspezifische und interspezifische Unterschiede im Grad der DNA-Methylierung festzustellen. Die differentiell methylierten CpG-Positionen lagen innerhalb eines repetitiven Alu-Sg1-Elements. Die Untersuchung des CCRK-Gens liefert eine umfassende Analyse der intra- und interspezifischen Variabilität der DNA-Methylierung einer Alu-Insertion in eine regulatorische Region. Die beobachteten Speziesunterschiede deuten darauf hin, dass die Methylierungsmuster des CCRK-Gens wahrscheinlich in Adaption an spezifische Anforderungen zur Feinabstimmung der CCRK-Regulation unter positiver Selektion evolvieren. Der Promotor des CCRK-Gens ist anfällig für epigenetische Modifikationen durch DNA-Methylierung, welche zu komplexen Transkriptionsmustern führen können. Durch ihre genomische Mobilität, ihren hohen CpG-Anteil und ihren Einfluss auf die Genexpression sind Alu-Insertionen exzellente Kandidaten für die Förderung von Veränderungen während der Entwicklungsregulation von Primatengenen. Der Vergleich der intra- und interspezifischen Methylierung von spezifischen Alu-Insertionen in anderen Genen und Geweben stellt eine erfolgversprechende Strategie dar.

No longitudinal mitochondrial DNA sequence changes in HIV-infected individuals with and without lipoatrophy

The potential for mitochondrial (mt) DNA mutation accumulation during antiretroviral therapy (ART), and preferential accumulation in patients with lipoatrophy compared with control participants, remains controversial. We sequenced the entire mitochondrial genome, both before ART and after ART exposure, in 29 human immunodeficiency virus (HIV)-infected Swiss HIV Cohort Study participants initiating a first-line thymidine analogue-containing ART regimen. No accumulation of mtDNA mutations or deletions was detected in 13 participants who developed lipoatrophy or in 16 control participants after significant and comparable ART exposure (median duration, 3.3 and 3.7 years, respectively). In HIV-infected persons, the development of lipoatrophy is unlikely to be associated with accumulation of mtDNA mutations detectable in peripheral blood.

Patentability and Scope of Protection for DNA Sequence

Since the mapping of the human genome and the technical innovations in the field of biotechnology, patent law has gone through great controversies. Protection is required for an investor to make an investment but how broad should the given protection be? Whether the invention is a mi- cro-organism capable of dissolving crude oil, or the gene of a soya plant, the genetic engineering required for their production entails vast amounts of capi- tal. The policy in that respect is tailored by legislative acts and judicial decisions, ensuring a fair balance be- tween the interests of patent right holders and third parties. However, the policy differs from jurisdiction to jurisdiction, thus creating inconsistencies with re- gards to the given protection to the same invention, and as a result this could deter innovation and pro- mote stagnation. The most active actors shaping the patent policy on an international level are the patent offices of the United States of America, Japan and the European Patent Organization. These three patent offices have set up a cooperation programme in order to promote and improve efficiency with regards to their patent policies on a global scale. However, recent judicial de- velopments have shown that the policy in respect to the field of biotechnology differs between the patent regimes of the United States of America and the two- layer system of the European Patent Organisation/ the European Union.

Models of DNA sequence evolution

Models of DNA sequence evolution and methods for estimating evolutionary distances are needed for studying the rate and pattern of molecular evolution and for inferring the evolutionary relationships of organisms or genes. In this dissertation, several new models and methods are developed.^ The rate variation among nucleotide sites: To obtain unbiased estimates of evolutionary distances, the rate heterogeneity among nucleotide sites of a gene should be considered. Commonly, it is assumed that the substitution rate varies among sites according to a gamma distribution (gamma model) or, more generally, an invariant+gamma model which includes some invariable sites. A maximum likelihood (ML) approach was developed for estimating the shape parameter of the gamma distribution $(\alpha)$ and/or the proportion of invariable sites $(\theta).$ Computer simulation showed that (1) under the gamma model, $\alpha$ can be well estimated from 3 or 4 sequences if the sequence length is long; and (2) the distance estimate is unbiased and robust against violations of the assumptions of the invariant+gamma model.^ However, this ML method requires a huge amount of computational time and is useful only for less than 6 sequences. Therefore, I developed a fast method for estimating $\alpha,$ which is easy to implement and requires no knowledge of tree. A computer program was developed for estimating $\alpha$ and evolutionary distances, which can handle the number of sequences as large as 30.^ Evolutionary distances under the stationary, time-reversible (SR) model: The SR model is a general model of nucleotide substitution, which assumes (i) stationary nucleotide frequencies and (ii) time-reversibility. It can be extended to SRV model which allows rate variation among sites. I developed a method for estimating the distance under the SR or SRV model, as well as the variance-covariance matrix of distances. Computer simulation showed that the SR method is better than a simpler method when the sequence length $L>1,000$ bp and is robust against deviations from time-reversibility. As expected, when the rate varies among sites, the SRV method is much better than the SR method.^ The evolutionary distances under nonstationary nucleotide frequencies: The statistical properties of the paralinear and LogDet distances under nonstationary nucleotide frequencies were studied. First, I developed formulas for correcting the estimation biases of the paralinear and LogDet distances. The performances of these formulas and the formulas for sampling variances were examined by computer simulation. Second, I developed a method for estimating the variance-covariance matrix of the paralinear distance, so that statistical tests of phylogenies can be conducted when the nucleotide frequencies are nonstationary. Third, a new method for testing the molecular clock hypothesis was developed in the nonstationary case. ^

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

The theoretical limits of DNA sequence discrimination by linked polyamides

Linked polyamides bind in the minor groove of double-stranded DNA in a partially sequence-specific manner. This report analyzes the theoretical limits of DNA sequence discrimination by linked polyamides composed of two to four different types of heterocyclic rings, determining (i) the optimal choice of base-binding specificity for each ring and (ii) the optimal design for a polyamide composed of these rings to target a given DNA sequence and designed to maximize the fraction of the total polyamide binding to the specified target sequence relative to all other sequences. The results show that, fortuitously, polyamides composed of pyrrole, a naturally occurring G-excluding element, and imidazole, a rationally designed G-favoring element, have features similar to the theoretical optimum design for polyamides composed of two different rings. The results also show that, in polyamides composed of two or three types of heterocyclic rings, choosing a nonspecific “placeholder” ring, which binds equally strongly to each of the four bases, along with one or two base-specific rings will often enhance sequence specificity over a polyamide composed entirely of base-specific rings.

The active digestion of uniparental chloroplast DNA in a single zygote of Chlamydomonas reinhardtii is revealed by using the optical tweezer

The non-Mendelian inheritance of organelle genes is a phenomenon common to almost all eukaryotes, and in the isogamous alga Chlamydomonas reinhardtii, chloroplast (cp) genes are transmitted from the mating type positive (mt+) parent. In this study, the preferential disappearance of the fluorescent cp nucleoids of the mating type negative (mt−) parent was observed in living young zygotes. To study the change in cpDNA molecules during the preferential disappearance, the cpDNA of mt+ or mt− origin was labeled separately with bacterial aadA gene sequences. Then, a single zygote with or without cp nucleoids was isolated under direct observation by using optical tweezers and investigated by nested PCR analysis of the aadA sequences. This demonstrated that cpDNA molecules are digested completely during the preferential disappearance of mt− cp nucleoids within 10 min, whereas mt+ cpDNA and mitochondrial DNA are protected from the digestion. These results indicate that the non-Mendelian transmission pattern of organelle genes is determined immediately after zygote formation.

Dual roles for DNA sequence identity and the mismatch repair system in the regulation of mitotic crossing-over in yeast

Sequence divergence acts as a potent barrier to homologous recombination; much of this barrier derives from an antirecombination activity exerted by mismatch repair proteins. An inverted repeat assay system with recombination substrates ranging in identity from 74% to 100% has been used to define the relationship between sequence divergence and the rate of mitotic crossing-over in yeast. To elucidate the role of the mismatch repair machinery in regulating recombination between mismatched substrates, we performed experiments in both wild-type and mismatch repair defective strains. We find that a single mismatch is sufficient to inhibit recombination between otherwise identical sequences, and that this inhibition is dependent on the mismatch repair system. Additional mismatches have a cumulative negative effect on the recombination rate. With sequence divergence of up to approximately 10%, the inhibitory effect of mismatches results mainly from antirecombination activity of the mismatch repair system. With greater levels of divergence, recombination is inefficient even in the absence of mismatch repair activity. In both wild-type and mismatch repair defective strains, an approximate log-linear relationship is observed between the recombination rate and the level of sequence divergence.

Chloroplast DNA footprints of postglacial recolonization by oaks

Recolonization of Europe by forest tree species after the last glaciation is well documented in the fossil pollen record. This spread may have been achieved at low densities by rare events of long-distance dispersal, rather than by a compact wave of advance, generating a patchy genetic structure through founder effects. In long-lived oak species, this structure could still be discernible by using maternally transmitted genetic markers. To test this hypothesis, a fine-scale study of chloroplast DNA (cpDNA) variability of two sympatric oak species was carried out in western France. The distributions of six cpDNA length variants were analyzed at 188 localities over a 200 × 300 km area. A cpDNA map was obtained by applying geostatistics methods to the complete data set. Patches of several hundred square kilometers exist which are virtually fixed for a single haplotype for both oak species. This local systematic interspecific sharing of the maternal genome strongly suggests that long-distance seed dispersal events followed by interspecific exchanges were involved at the time of colonization, about 10,000 years ago.

A Chloroplast DNA Helicase II from Pea That Prefers Fork-Like Replication Structures

A DNA helicase, called chloroplast DNA (ctDNA) helicase II, was purified to apparent homogeneity from pea (Pisum sativum). The enzyme contained intrinsic, single-stranded, DNA-dependent ATPase activity and an apparent molecular mass of 78 kD on sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The DNA helicase was markedly stimulated by DNA substrates with fork-like replication structures. A 5′-tailed fork was more active than the 3′-tailed fork, which itself was more active than substrates without a fork. The direction of unwinding was 3′ to 5′ along the bound strand, and it failed to unwind blunt-ended duplex DNA. DNA helicase activity required only ATP or dATP hydrolysis. The enzyme also required a divalent cation (Mg2+>Mn2+>Ca2+) for its unwinding activity and was inhibited at 200 mm KCl or NaCl. This enzyme could be involved in the replication of ctDNA. The DNA major groove-intercalating ligands nogalamycin and daunorubicin were inhibitory to unwinding (Ki approximately 0.85 μm and 2.2 μm, respectively) and ATPase (Ki approximately 1.3 μm and 3.0 μm, respectively) activities of pea ctDNA helicase II, whereas ellipticine, etoposide (VP-16), and camptothecin had no effect on the enzyme activity. These ligands may be useful in further studies of the mechanisms of chloroplast helicase activities.

DNA sequence insertion and evolutionary variation in gene regulation.

Current evidence on the long-term evolutionary effect of insertion of sequence elements into gene regions is reviewed, restricted to cases where a sequence derived from a past insertion participates in the regulation of expression of a useful gene. Ten such examples in eukaryotes demonstrate that segments of repetitive DNA or mobile elements have been inserted in the past in gene regions, have been preserved, sometimes modified by selection, and now affect control of transcription of the adjacent gene. Included are only examples in which transcription control was modified by the insert. Several cases in which merely transcription initiation occurred in the insert were set aside. Two of the examples involved the long terminal repeats of mammalian endogenous retroviruses. Another two examples were control of transcription by repeated sequence inserts in sea urchin genomes. There are now six published examples in which Alu sequences were inserted long ago into human gene regions, were modified, and now are central in control/enhancement of transcription. The number of published examples of Alu sequences affecting gene control has grown threefold in the last year and is likely to continue growing. Taken together, all of these examples show that the insertion of sequence elements in the genome has been a significant source of regulatory variation in evolution.

Binding of a hairpin polyamide in the minor groove of DNA: sequence-specific enthalpic discrimination.

Hairpin polyamides are synthetic ligands for sequence-specific recognition in the minor groove of double-helical DNA. A thermodynamic characterization of the DNA-binding properties exhibited by a six-ring hairpin polyamide, ImPyPy-gamma-PyPyPy-beta-Dp (where Im = imidazole, Py = pyrrole, gamma = gamma-aminobutyric acid, beta = beta-alanine, and Dp = dimethylaminopropylamide), reveals an approximately 1-2 kcal/mol greater affinity for the designated match site, 5'-TGTTA-3', relative to the single base pair mismatch sites, 5'-TGGTA-3' and 5'-TATTA-3'. The enthalpy and entropy data at 20 degrees C reveal this sequence specificity to be entirely enthalpic in origin. Correlations between the thermodynamic driving forces underlying the sequence specificity exhibited by ImPyPy-gamma-PyPyPy-beta-Dp and the structural properties of the heterodimeric complex of PyPyPy and ImPyPy bound to the minor groove of DNA provide insight into the molecular forces that govern the affinity and specificity of pyrrole-imidazole polyamides.

Detection of JC virus DNA sequence and expression of the viral oncoprotein, tumor antigen, in brain of immunocompetent patient with oligoastrocytoma.

We describe molecular and clinical findings in an immunocompetent patient with an oligoastrocytoma and the concomitant presence of the human papovavirus, JC virus (JCV), which is the etiologic agent of the subacute, debilitating demyelinating disease, progressive multifocal leukoencephalopathy. Histologic review revealed a glial neoplasm consisting primarily of a moderately cellular oligodendroglioma with distinct areas of a fibrillary astrocytoma. Immunohistochemical analysis revealed nuclear staining of tumor cells with antibodies against the viral oncoprotein [tumor antigen (T antigen)], the proliferation marker (Ki67), and the cellular proliferation regulator (p53). Using primers specific to the JCV control region, PCR yielded amplified DNA that was identical to the control region of the Mad-4 strain of the virus. PCR analysis demonstrated the presence of the genome for the viral oncoprotein, T antigen, and results from primer extension studies revealed synthesis of the viral early RNA for T antigen in the tumor tissues. The presence of viral T antigen in the tumor tissue was further demonstrated by immunoblot assay. To our knowledge, this is the first report of the presence of JCV DNA, RNA, and T antigen in tissue in which viral T antigen is localized to tumor cell nuclei and suggests the possible association of JCV with some glial neoplasms.

Chloroplast DNA evidence of colonization, adaptive radiation, and hybridization in the evolution of the Macaronesian flora.

Most evolutionary studies of oceanic islands have focused on the Pacific Ocean. There are very few examples from the Atlantic archipelagos, especially Macaronesia, despite their unusual combination of features, including a close proximity to the continent, a broad range of geological ages, and a biota linked to a source area that existed in the Mediterranean basin before the late Tertiary. A chloroplast DNA (cpDNA) restriction site analysis of Argyranthemum (Asteraceae: Anthemideae), the largest endemic genus of plants of any volcanic archipelago in the Atlantic Ocean, was performed to examine patterns of plant evolution in Macaronesia. cpDNA data indicated that Argyranthemum is a monophyletic group that has speciated recently. The cpDNA tree showed a weak correlation with the current sectional classification and insular distribution. Two major cpDNA lineages were identified. One was restricted to northern archipelagos--e.g., Madeira, Desertas, and Selvagens--and the second comprised taxa endemic to the southern archipelago--e.g., the Canary Islands. The two major radiations identified in the Canaries are correlated with distinct ecological habitats; one is restricted to ecological zones under the influence of the northeastern trade winds and the other to regions that are not affected by these winds. The patterns of phylogenetic relationships in Argyranthemum indicate that interisland colonization between similar ecological zones is the main mechanism for establishing founder populations. This phenomenon, combined with rapid radiation into distinct ecological zones and interspecific hybridization, is the primary explanation for species diversification.