Search PubMed⌕ Search

Biomedical subjects

A von Haeseler

Publications and source records attributed to A von Haeseler.

At least 19 recordsLinked to original sources

Phylogenetic analysis of rubella virus including new genotype I isolates.

Infection during the first trimester of gestation with rubella virus (RV) is highly teratogenic. Embryopathy is a frequent outcome of the primary natural infection with wild type RV during pregnancy while accidental immunisation with life attenuated vaccine has apparently little or no adverse effect. Although the nucleic acid sequence of RV is exceptionally stable, differences between the vaccine and wild type viruses could play a role in the pathogenesis of intrauterine RV infection. Phylogenetic analysis of eight complete sequences (including two new isolates described in this paper, three vaccine strains, three cell culture adapted wild type viruses) confirms a striking low divergence of the RV genome. A variable region (amino acid residues 697-800) within the gene coding for the nonstructural protein NSP1 was defined. Phylogenetic analysis revealed a strong positive selection in this region. Multiple passages in vivo or in vitro did not account for this variability. As the function of the variable region has not yet been elucidated, reasons for and significance of positive selection are still speculative. It is conceivable that the variable region in NSP1 contributes to the molecular basis of RV embryopathy and other complications of postnatal RV infection.

Animals↗

DNA sequences from multiple amplifications reveal artifacts induced by cytosine deamination in ancient DNA.

We show that DNA molecules amplified by PCR from DNA extracted from animal bones and teeth that vary in age between 25 000 and over 50 000 years carry C-->T and G-->A substitutions. These substitutions can reach high proportions among the molecules amplified and are due to the occurrence of modified deoxycytidine residues in the template DNA. If the template DNA is treated with uracil N-glycosylase, these substitutions are dramatically reduced. They are thus likely to result from deamination of deoxycytidine residues. In addition, 'jumping PCR', i.e. the occurrence of template switching during PCR, may contribute to these substitutions. When DNA sequences are amplified from ancient DNA extracts where few template molecules initiate the PCR, precautions such as DNA sequence determination of multiple clones derived from more than one independent amplification are necessary in order to reduce the risk of determination of incorrect DNA sequences. When such precautionary measures are taken, errors induced by damage to the DNA template are unlikely to be more frequent than approximately 0.1% even under the unlikely scenario where each amplification starts from a single template molecule.

Animals↗

Assessing variability by joint sampling of alignments and mutation rates.

When two sequences are aligned with a single set of alignment parameters, or when mutation parameters are estimated on the basis of a single "optimal" sequence alignment, the variability of both the alignment and the estimated parameters can be seriously underestimated. To obtain a more realistic impression of the actual uncertainty, we propose sampling sequence alignments and mutation parameters simultaneously from their joint posterior distribution given the two original sequences. We illustrate our method with human and orangutan sequences from the hyper variable region I and with gene-pseudogene pairs.

Animals↗

Quartet-mapping, a generalization of the likelihood-mapping procedure.

Likelihood-mapping (LM) was suggested as a method of displaying the phylogenetic content of an alignment. However, statistical properties of the method have not been studied. Here we analyze the special case of a four-species tree generated under a range of evolution models and compare the results with those of a natural extension of the likelihood-mapping approach, geometry-mapping (GM), which is based on the method of statistical geometry in sequence space. The methods are compared in their abilities to indicate the correct topology. The performance of both methods in detecting the star topology is especially explored. Our results show that LM tends to reject a star tree more often than GM. When assumptions about the evolutionary model of the maximum-likelihood reconstruction are not matched by the true process of evolution, then LM shows a tendency to favor one tree, whereas GM correctly detects the star tree except for very short outer branch lengths with a statistical significance of >0.95 for all models. LM, on the other hand, reconstructs the correct bifurcating tree with a probability of >0.95 for most branch length combinations even under models with varying substitution rates. The parameter domain for which GM recovers the true tree is much smaller. When the exterior branch lengths are larger than a (analytically derived) threshold value depending on the tree shape (rather than the evolutionary model), GM reconstructs a star tree rather than the true tree. We suggest a combined approach of LM and GM for the evaluation of starlike trees. This approach offers the possibility of testing for significant positive interior branch lengths without extensive statistical and computational efforts.

Computer Simulation↗

A developmentally regulated aconitase related to iron-regulatory protein-1 is localized in the cytoplasm and in the mitochondrion of Trypanosoma brucei.

Mitochondrial energy metabolism and Krebs cycle activities are developmentally regulated in the life cycle of the protozoan parasite Trypanosoma brucei. Here we report cloning of a T. brucei aconitase gene that is closely related to mammalian iron-regulatory protein 1 (IRP-1) and plant aconitases. Kinetic analysis of purified recombinant TbACO expressed in Escherichia coli resulted in a K(m) (isocitrate) of 3 +/- 0.4 mM, similar to aconitases of other organisms. This was unexpected since an arginine conserved in the aconitase protein family and crucial for substrate positioning in the catalytic center and for activity of pig mitochondrial aconitase (Zheng, L., Kennedy, M. C., Beinert, H., and Zalkin, H. (1992) J. Biol. Chem. 267, 7895-7903) is substituted by leucine in the TbACO sequence. Expression of the 98-kDa TbACO was shown to be lowest in the slender bloodstream stage of the parasite, 8-fold elevated in the stumpy stage, and increased a further 4-fold in the procyclic stage. The differential expression of TbACO protein contrasted with only minor changes in TbACO mRNA, indicating translational or post-translational mechanisms of regulation. Whereas animal cells express two distinct compartmentalized aconitases, mitochondrial aconitase and cytoplasmic aconitase/IRP-1, TbACO accounts for total aconitase activity in trypanosomes. By cell fractionation and immunofluorescence microscopy, we show that native as well as a transfected epitope-tagged TbACO localizes in both the mitochondrion (30%) and in the cytoplasm (70%). Together with phylogenetic reconstructions of the aconitase family, this suggests that animal IRPs have evolved from a multicompartmentalized ancestral aconitase. The possible functions of a cytoplasmic aconitase in trypanosomes are discussed.

Aconitate Hydratase↗

A coalescent approach to study linkage disequilibrium between single-nucleotide polymorphisms.

We present the results of extensive simulations that emulate the development and distribution of linkage disequilibrium (LD) between single-nucleotide polymorphisms (SNPs) and a gene locus that is phenotypically stratified into two classes (disease phenotype and wild-type phenotype). Our approach, based on coalescence theory, allows an explicit modeling of the demographic history of the population without conditioning on the age of the mutation, and serves as an efficient tool to carry out simulations. More specifically, we compare the influence that a constant population size or an exponentially growing population has on the amount of LD. These results indicate that attempts to locate single disease genes are most likely successful in small and constant populations. On the other hand, if we consider an exponentially growing population that started to expand from an initially constant population of reasonable size, then our simulations indicate a lower success rate. The power to detect association is enhanced if haplotypes constructed from several SNPs are used as markers. The versatility of the coalescence approach also allows the analysis of other relevant factors that influence the chances that a disease gene will be located. We show that several alleles leading to the same disease have no substantial influence on the amount of LD, as long as the differences between the disease-causing alleles are confined to the same region of the gene locus and as long as each allele occurs in an appreciable frequency. Our simulations indicate that mapping of less-frequent diseases is more likely to be successful. Moreover, we show that successful attempts to map complex diseases depend crucially on the phenotype-genotype correlations of all alleles at the disease locus. An analysis of lipoprotein lipase data indicates that our simulations capture the major features of LD occurring in biological data.

Alleles↗

C-->U editing of apolipoprotein B mRNA in marsupials: identification and characterisation of APOBEC-1 from the American opossum Monodelphus domestica.

The C->U editing of RNA is widely found in plant and animal species. In mammals it is a discrete process confined to the editing of apolipoprotein B (apoB) mRNA in eutherians and the editing of the mitochondrial tRNA for glycine in marsupials. Here we have identified and characterised apoB mRNA editing in the American opossum Monodelphus domestica. The apoB mRNA editing site is highly conserved in the opossum and undergoes complete editing in the small intestine, but not in the liver or other tissues. Opossum APOBEC-1 cDNA was cloned, sequenced and expressed. The encoded protein is similar to APOBEC-1 of eutherians. Motifs previously identified as involved in zinc binding, RNA binding and catalysis, nuclear localisation and a C-terminal leucine-rich domain are all conserved. Opossum APOBEC-1 contains a seven amino acid C-terminal extension also found in humans and rabbits, but not present in rodents. The opossum APOBEC-1 gene has the same intron/exon organisation in the coding sequence as the eutherian gene. Northern blot and RT-PCR analyses and an editing assay indicate that no APOBEC-1 was expressed in the liver. Thus the far upstream promoter responsible for hepatic expression in rodents does not operate in the opossum. An APOBEC-1-like enzyme such as might be involved in C->U RNA editing of tRNA in marsupial mitochondria was not demonstrated. The activity of opossum APOBEC-1 in the presence of both chicken and rodent auxiliary editing proteins was comparable to that of other mammals. These studies extend the origins of APOBEC-1 back 170 000 000 years to marsupials and help bridge the gap in the origins of this RNA editing process between birds and eutherian mammals.

APOBEC-1 Deaminase↗

HvrBase: compilation of mtDNA control region sequences from primates.

HvrBase is a compilation of human and ape mtDNA control region sequences. Sequences and related information on individuals, such as from where the sequences were obtained, is stored in three ASCII files as described previously. Moreover, the collection is also available as Mac/PC database application with a graphical user interface. It can be accessed through the WWW at URL http://www.eva.mpg.de/hvrbase. The current collection comprises 5846 human sequences from hypervariable region I (HVRI) and 2302 human sequences from hypervariable region II (HVRII). From apes, 295 HVRI sequences and 13 HVRII sequences are available.

Animals↗

Distance measures in terms of substitution processes.

Phylogenetic reconstruction from DNA or amino acid sequences relies heavily on suitable distance measures. A number of new distance measures (asynchronous, LogDet, and paralinear distances) which possess the desired property of tree additivity under fairly general models of sequence evolution have been proposed recently, but they are not well understood from a mechanistic point of view. We review them here in a unifying framework, which is the substitution process in continuous time. The emerging interpretation will also clarify the relationship among these distance measures. We also tackle situations with site-to-site variation of substitution rates which is well known to cause non-additive distances and inconsistent branch lengths. For homogeneous, stationary, time-reversible models, this may be repaired provided that the distribution of rates is known. In contrast, we will show that, for non-stationary models, different tree topologies may produce identical joint distributions of letters in pairs of sequences, given the same distribution of rates. This precludes the existence of any tree-additive pairwise distance measure.

Genetic Heterogeneity↗

Toward assigning helical regions in alignments of ribosomal RNA and testing the appropriateness of evolutionary models.

We suggest a nucleotide substitution model that takes correlation between base-paired nucleotides into account. The model includes the estimation of the transition-transversion ratio and allows inference of the shape parameter of a discrete gamma distribution to include rate heterogeneity. A Cox-test statistic, applied to a diatom ribosomal RNA alignment, shows that the suggested correlation model explains evolution of the stem region better than usual independence models. Moreover, the Cox-test procedure is extended to shed some light upon the problem of assigning helical regions in a secondary structure based alignment. This approach provides an estimate of the percentage of stem positions that do not appear to be correlated.

Base Sequence↗

DNA sequence variation in a non-coding region of low recombination on the human X chromosome.

DNA sequence variation has become a major source of insight regarding the origin and history of our species as well as an important tool for the identification of allelic variants associated with disease. Comparative sequencing of DNA has to date focused mainly on mitochondrial (mt) DNA, which due to its apparent lack of recombination and high evolutionary rate lends itself well to the study of human evolution. These advantages also entail limitations. For example, the high mutation rate of mtDNA results in multiple substitutions that make phylogenetic analysis difficult and, because mtDNA is maternally inherited, it reflects only the history of females. For the history of males, the non-recombining part of the paternally inherited Y chromosome can be studied. The extent of variation on the Y chromosome is so low that variation at particular sites known to be polymorphic rather than entire sequences are typically determined. It is currently unclear how some forms of analysis (such as the coalescent) should be applied to such data. Furthermore, the lack of recombination means that selection at any locus affects all 59 Mb of DNA. To gauge the extent and pattern of point substitutional variation in non-coding parts of the human genome, we have sequenced 10 kb of non-coding DNA in a region of low recombination at Xq13.3. Analysis of this sequence in 69 individuals representing all major linguistic groups reveals the highest overall diversity in Africa, whereas deep divergences also exist in Asia. The time elapsed since the most recent common ancestor (MRCA) is 535,000+/-119,000 years. We expect this type of nuclear locus to provide more answers about the genetic origin and history of humans.

Africa↗

Pattern of nucleotide substitution and rate heterogeneity in the hypervariable regions I and II of human mtDNA.

This study provides a comprehensive survey of the complex pattern of nucleotide substitution in the control region of human mtDNA, which is of central importance to the studies of human evolution. A total of 1229 different hypervariable region I (HVRI) and 385 different hypervariable region II (HVRII) sequences were analyzed using a complex substitution model. Moreover, we suggest a new method to assign relative rates to each site in the sequence. Estimates are based on maximum-likelihood methods applied to randomly selected subsets of sequences. Our results indicate that the rate of substitution in HVRI is approximately twice as high as in HVRII and that this difference is mainly due to a higher frequency of pyrimidine transitions in HVRI. However, rate heterogeneity is more pronounced in HVRII.

DNA, Mitochondrial↗

Compilation of human mtDNA control region sequences.

This paper describes the organisation of a database for human mitochondrial control-region sequences. The data are divided into three ASCII files that contain aligned sequences from the hypervariable region I (HVRI), from the hypervariable region II (HVRII), and the available information about the individuals, from whom the sequences stem. The current collection comprises 4079 HVRI and 969 HVRII sequences. From 728 individuals sequences of both HVRI and HVRII are available. For easy access, the collection is made available to the scientific community via World Wide Web at URL http://www.zi.biologie.uni-muenchen.de/[symbol: see text]meyers/mtdna.html

Computer Communication Networks↗

Evolution of DNA or amino acid sequences with dependent sites.

A framework is outlined to study the evolution of DNA or amino acid sequences, if sequence sites do not evolve independently. The units of evolution are nonoverlapping subsequences of length l. Each subsequence evolves independently of the others, but within a subsequence the sequences show a Markov order one dependency. We describe an algorithm to mimic the evolution of such sequences. The influence of dependencies between sites on distance estimates and the reliability of tree reconstruction methods is investigated. We show that an inappropriate model of sequence evolution in the tree reconstruction process will lead to a nonempty Felsenstein zone. Finally, we describe a method to infer l from sequence data. Examples from the evolution of DNA sequences as well as from amino acids are given.

Algorithms↗

Inference of population history using a likelihood approach.

We introduce an approach to revealing the likelihood of different population histories that utilizes an explicit model of sequence evolution for the DNA segment under study. Based on a phylogenetic tree reconstruction method we show that a Tamura-Nei model with heterogeneous mutation rates is a fair description of the evolutionary process of the hypervariable region I of the mitochondrial DNA from humans. Assuming this complex model still allows the estimation of population history parameters, we suggest a likelihood approach to conducting statistical inference within a class of expansion models. More precisely, the likelihood of the data is based on the mean pairwise differences between DNA sequences and the number of variable sites in a sample. The use of likelihood ratios enables comparison of different hypotheses about population history, such as constant population size during the past or an increase or decrease of population size starting at some point back in time. This method was applied to show that the population of the Basques has expanded, whereas that of the Biaka pygmies is most likely decreasing. The Nuu-Chah-Nulth data are consistent with a model of constant population.

Base Composition↗

Mitochondrial DNA sequencing of shed hairs and saliva on robbery caps: sensitivity and matching probabilities.

Sequencing of mitochondrial DNA (mtDNA) has been used for human identification based on teeth and skeletal remains. Here, we describe an amplification system for the mtDNA control region (D-loop) suited for the analysis of shed hair, which constitutes the most common biological evidence material in forensic investigations. The success rate was over 90% when applied to evidence materials such as shed hair, saliva stains and saliva on stamps. The analysis of evidence materials collected from three similar robberies revealed the presence of mtDNA sequences identical to those of the suspects in the three crimes. The use of mtDNA control region sequences for individual identification was evaluated. The probability of identity by chance for the mtDNA types of the suspects in the robberies was found to vary between Pr = 0.017 - < 0.0017, depending on the reference population used, emphasizing the need for large population databases to obtain the appropriate estimate.

Base Sequence↗

A coalescent approach to the polymerase chain reaction.

A versatile algorithm is developed to model PCR on a computer. The method is based on a modification of the coalescent process and provides a general framework to analyse data from PCR. It allows for incorporation of the dynamics of the replication process as described in terms of the number of starting template molecules and cycle-dependent PCR efficiency. The simulation method generates, as a first step, the genealogy of a set of sequences sampled from a final PCR product. In a second step a mutation process is superimposed and the resulting data set is analysed. The efficiency of our algorithm enables us to get reliable approximations of various sample distributions. We demonstrate the relevance of our method with two applications: maximum likelihood estimation of the error rate in PCR and a test of homogeneity of the template.

Algorithms↗