Neil Howell

Reduced-Median-Network Analysis of Complete Mitochondrial DNA Coding-Region Sequences for the Major African, Asian, and European Haplogroups

The evolution of the human mitochondrial genome is characterized by the emergence of ethnically distinct lineages or haplogroups. Nine European, seven Asian (including Native American), and three African mitochondrial DNA (mtDNA) haplogroups have been identified previously on the basis of the presence or absence of a relatively small number of restriction-enzyme recognition sites or on the basis of nucleotide sequences of the D-loop region. We have used reduced-median-network approaches to analyze 560 complete European, Asian, and African mtDNA coding-region sequences from unrelated individuals to develop a more complete understanding of sequence diversity both within and between haplogroups. A total of 497 haplogroup-associated polymorphisms were identified, 323 (65%) of which were associated with one haplogroup and 174 (35%) of which were associated with two or more haplogroups. Approximately one-half of these polymorphisms are reported for the first time here. Our results confirm and substantially extend the phylogenetic relationships among mitochondrial genomes described elsewhere from the major human ethnic groups. Another important result is that there were numerous instances both of parallel mutations at the same site and of reversion (i.e., homoplasy). It is likely that homoplasy in the coding region will confound evolutionary analysis of small sequence sets. By a linkage-disequilibrium approach, additional evidence for the absence of human mtDNA recombination is presented here.

Mammalian mitochondrial genetics: heredity, heteroplasmy and disease

Mammalian mitochondrial DNA (mtDNA) is present at high copy number (10(3)-10(4) copies) in virtually all cells of the body. The mitochondrial genome shows strict maternal inheritance and the vast majority of copies are identical at birth (homoplasmy). Occasionally, a subpopulation of mtDNA molecules carry a pathogenic mutation. When this heteroplasmic mtDNA is present during embryogenesis, it can lead to a variety of clinical symptoms predominantly affecting muscle and nerve, but also affecting other tissues. While the importance of mitochodrial heteroplasmy in human disease is unquestioned, we remain largely ignorant of many fundamental aspects of mitochondrial genetics. How do mutations arise and can they be repaired, what influences the segregation and fixation of heteroplasmic mtDNA, do levels of heteroplasmy fluctuate during life, is it possible to modulate these levels by external intervention and, finally, can we predict the segregation and transmission of a mutant genome? The aim of this article is to summarize and discuss recent observations that have addressed several of these fundamental issues and to reiterate how much we still have to learn about mitochondrial genetics.

The Epidemiology of Leber Hereditary Optic Neuropathy in the North East of England

We performed the first population-based clinical and molecular genetic study of Leber hereditary optic neuropathy (LHON) in a population of 2,173,800 individuals in the North East of England. We identified 16 genealogically unrelated families who harbor one of the three primary mitochondrial DNA (mtDNA) mutations that cause LHON. Two of these families were found to be linked genetically to a common maternal founder. A de novo mtDNA mutation (G3460A) was identified in one family. The minimum point prevalence of visual failure due to LHON within this population was 3.22 per 100,000 (95% CI 2.47-3.97 per 100,000), and the minimum point prevalence for mtDNA LHON mutations was 11.82 per 100,000 (95% CI 10.38-13.27 per 100,000). These results indicate that LHON is not rare but has a population prevalence similar to autosomally inherited neurological disorders. The majority of individuals harbored only mutant mtDNA (homoplasmy), but heteroplasmy was detected in approximately 12% of individuals. Overall, however, approximately 33% of families with LHON had at least one heteroplasmic individual. The high incidence of heteroplasmy in pedigrees with LHON raises the possibility that a closely related maternal relative of an index case may not harbor the mtDNA mutation, highlighting the importance of molecular genetic testing for each maternal family member seeking advice about their risks of visual failure.

Clinical expression of Leber hereditary optic neuropathy is affected by the mitochondrial DNA-haplogroup background.

Leber hereditary optic neuropathy (LHON) is due primarily to one of three common point mutations of mitochondrial DNA (mtDNA), but the incomplete penetrance implicates additional genetic or environmental factors in the pathophysiology of the disorder. Both the 11778G-->A and 14484T-->C LHON mutations are preferentially found on a specific mtDNA genetic background, but 3460G-->A is not. However, there is no clear evidence that any background influences clinical penetrance in any of these mutations. By studying 3,613 subjects from 159 LHON-affected pedigrees, we show that the risk of visual failure is greater when the 11778G-->A or 14484T-->C mutations are present in specific subgroups of haplogroup J (J2 for 11778G-->A and J1 for 14484T-->C) and when the 3460G-->A mutation is present in haplogroup K. By contrast, the risk of visual failure is significantly less when 11778G-->A occurs in haplogroup H. Substitutions on MTCYB provide an explanation for these findings, which demonstrate that common genetic variants have a marked effect on the expression of an ostensibly monogenic mtDNA disorder.

The Pedigree Rate of Sequence Divergence in the Human Mitochondrial Genome. There Is a Difference between Phylogenetic and Pedigree Rates

We have extended our previous analysis of the pedigree rate of control-region divergence in the human mitochondrial genome. One new germline mutation in the mitochondrial DNA (mtDNA) control region was detected among 185 transmission events (generations) from five Leber hereditary optic neuropathy (LHON) pedigrees. Pooling the LHON pedigree analyses yields a control-region divergence rate of 1.0 mutation/bp/10(6) years (Myr). When the results from eight published studies that used a similar approach were pooled with the LHON pedigree studies, totaling >2,600 transmission events, a pedigree divergence rate of 0.95 mutations/bp/Myr for the control region was obtained with a 99.5% confidence interval of 0.53-1.57. Taken together, the cumulative results support the original conclusion that the pedigree divergence rate for the control region is approximately 10-fold higher than that obtained with phylogenetic analyses. There is no evidence that any one factor explains this discrepancy, and the possible roles of mutational hotspots (rate heterogeneity), selection, and random genetic drift and the limitations of phylogenetic approaches to deal with high levels of homoplasy are discussed. In addition, we have extended our pedigree analysis of divergence in the mtDNA coding region. Finally, divergence of complete mtDNA sequences was analyzed in two tissues, white blood cells and skeletal muscle, from each of 17 individuals. In three of these individuals, there were four instances in which an mtDNA mutation was found in one tissue but not in the other. These results are discussed in terms of the occurrence of somatic mtDNA mutations.

The inheritance of mitochondrial DNA heteroplasmy: random drift, selection or both?

The mammalian mitochondrial genome (mtDNA) is a small double-stranded DNA molecule that is exclusively transmitted down the maternal line. Pathogenic mtDNA mutations are usually heteroplasmic, with a mixture of mutant and wild-type mtDNA within the same organism. A woman harbouring one of these mutations transmits a variable amount of mutant mtDNA to each offspring. This can result in a healthy child or an infant with a devastating and fatal neurological disorder. Understanding the biological basis of this uncertainty is one of the principal challenges facing scientists and clinicians in the field of mitochondrial genetics.

Comparative Genomics and the Evolution of Human Mitochondrial DNA: Assessing the Effects of Selection

This article provides evidence that selection has been a significant force during the evolution of the human mitochondrial genome. Both gene-by-gene and whole-genome approaches were used here to assess selection in the 560 mitochondrial DNA (mtDNA) coding-region sequences that were used previously for reduced-median-network analysis. The results of the present analyses were complex, in that the action of selection was not indicated by all tests, but this is not surprising, in view of the characteristics and limitations of the different analytical methods. Despite these limitations, there is evidence for both gene-specific and lineage-specific variation in selection. Whole-genome sliding-window approaches indicated a lack of selection in large-scale segments of the coding region. In other tests, we analyzed the ratio of nonsynonymous-to-synonymous substitutions in the 13 protein-encoding mtDNA genes. The most straightforward interpretation of those results is that negative selection has acted on the mtDNA during evolution. Single-gene analyses indicated significant departures from neutrality in the CO1, ND4, and ND6 genes, although the data also suggested the possible operation of positive selection on the AT6 gene. Finally, our results and those of other investigators do not support a simple model in which climatic adaptation has been a major force during human mtDNA evolution.

Mitochondrial DNA and survival after sepsis: a prospective study

BACKGROUND Human genome evolution has been shaped by infectious disease. Although most genetic studies have focused on the immune system, recovery after sepsis is directly related to physiological reserve that is critically dependent on mitochondrial function. We investigated whether haplogroup H, the most common type of mitochondrial DNA (mtDNA) in Europe, contributes to the subtle genetic variation in survival after sepsis. METHODS In a prospective study, we included 150 individuals who were sequentially admitted to the intensive care unit in a hospital in Newcastle upon Tyne, UK. After clinical data were obtained, patients underwent mtDNA haplotyping by analysis with PCR and restriction fragment length polymorphism. As endpoints, we used death during the 6-month period or survival at 6 months. FINDINGS Follow-up was complete for all study participants, although the haplotype of two patients could not be reliably determined. On admission to the intensive care unit, the frequency of mtDNA haplogroup H in study patients did not differ between study patients admitted with severe sepsis and 542 age-matched controls from the northeast of England. MtDNA haplogroup H was a strong independent predictor of outcome during severe sepsis, conferring a 2.12-fold (95% CI 1.02-4.43) increased chance of survival at 180 days compared with individuals without the haplogroup H. INTERPRETATION Although haplogroup H is the most recent addition to the group of European mtDNA, paradoxically it is also the most common. Increased survival after sepsis provides one explanation for this observation. MtDNA haplotyping offers a new means of risk stratification of patients with severe infections, which suggests new avenues for therapeutic intervention.

Evolutionary conservation of protein regions in the protonmotive cytochromeb and their possible roles in redox catalysis

SummaryThe amino acid sequences of the protonmotive cytochromeb from seven representative and phylogenetically diverse species have been compared to identify protein regions or segments that are conserved during evolution. The sequences analyzed included both prokaryotic and eukaryotic examples as well as mitochondrial cytochromeb and chloroplastb6 proteins. The principal conclusion from these analyses is that there are five protein regions-each comprising about 20 amino acid residues—that are consistently conserved during evolution. These domains are evident despite the low density of invariant residues. The two most highly conserved regions, spanning approximately consensus residues 130–150 and 270–290, are located in extramembrane loops and are hypothesized to constitute part of the Qo reaction center. The intramembrane, hydrophobic protein regions containing the heme-ligating histidines are also conserved during evolution. It was found, however, that the conservation of the protein segments extramembrane to the histidine residues ligating the low potential b566 heme group showed a higher degree of sequence conservation. The location of these conserved regions suggests that these extramembrane segments are also involved in forming the Qo reaction center. A protein segment putatively constituting a portion of the Qi reaction center, located approximately in the region spanned by consensus residues 20–40, is conserved in species as divergent as mouse andRhodobacter. This region of the protein shows substantially less sequence conservation in the chloroplast cytochromeb6. The catalytic role of these conserved regions is strongly supported by locations of residues that are altered in mutants resistant to inhibitors of cytochromeb electron transport.

Mutational analysis of the mouse mitochondrial cytochrome b gene

The protonmotive cytochrome b protein of the mitochondrial bc1 respiratory chain complex contains two reactions centers, designated Qo and Qi, which can be distinguished by the effects of different inhibitors. The nucleotide sequences have been determined of the mitochondrial cytochrome b genes from a series of mouse cell mutants selected for increased inhibitor resistance. Each mutant contains a single nucleotide change which results in an amino acid substitution. When the proximity of the altered amino acid residues to the histidines involved in heme ligation is considered, the results support a model for cytochrome b folding in which there are eight transmembrane domains rather than the nine of the Widger-Saraste model. Replacement of the Gly38 residue by valine results in resistance to the Qi inhibitors antimycin A and funiculosin but not 2-n-heptyl-hydroxyquinoline-N-oxide. Based upon sequence comparisons of mitochondrial and bacterial cytochrome b and chloroplast b6 proteins, the region of the molecule involved in antimycin binding is as highly conserved as those domains involved in heme ligation. It is suggested that the antimycin binding domain of cytochrome b is involved in forming the Qi reaction center. Alterations of the Gly142 and Thr147 residues result in resistance to myxothiazol and stimatellin, respectively. While both inhibitors block the Qo reaction center, the two mutations do not confer cross-resistance to each other. This region of cytochrome b is the most highly conserved during evolution and these inhibitor binding sites probably occur within the protein domain constituting the Qo reaction center. In addition, there is a less conserved region of the protein, defined by the Leu294 residue, which may function in binding the hydrophobic portions of Qo inhibitors.