Evolution of the metazoan mitochondrial replicase

Comparative mitochondrial genomics in Nematoda reveal astonishing variation in compositional biases and substitution rates indicative of multi-level selection

BACKGROUND

Nematodes are the most abundant and diverse metazoans on Earth, and are known to significantly affect ecosystem functioning. A better understanding of their biology and ecology, including potential adaptations to diverse habitats and lifestyles, is key to understanding their response to global change scenarios. Mitochondrial genomes offer high species level characterization, low cost of sequencing, and an ease of data handling that can provide insights into nematode evolutionary pressures.

RESULTS

Generally, nematode mitochondrial genomes exhibited similar structural characteristics (e.g., gene size and GC content), but displayed remarkable variability around these general patterns. Compositional strand biases showed strong codon position specific G skews and relationships with nematode life traits (especially parasitic feeding habits) equal to or greater than with predicted phylogeny. On average, nematode mitochondrial genomes showed low non-synonymous substitution rates, but also high clade specific deviations from these means. Despite the presence of significant mutational saturation, non-synonymous (dN) and synonymous (dS) substitution rates could still be significantly explained by feeding habit and/or habitat. Low ratios of dN:dS rates, particularly associated with the parasitic lifestyles, suggested the presence of strong purifying selection.

CONCLUSIONS

Nematode mitochondrial genomes demonstrated a capacity to accumulate diversity in composition, structure, and content while still maintaining functional genes. Moreover, they demonstrated a capacity for rapid evolutionary change pointing to a potential interaction between multi-level selection pressures and rapid evolution. In conclusion, this study helps establish a background for our understanding of the potential evolutionary pressures shaping nematode mitochondrial genomes, while outlining likely routes of future inquiry.

A century of mitochondrial research, 1922-2022

Although recognized earlier as subcellular entities by microscopists, mitochondria have been the subject of functional studies since 1922, when their biochemical similarities with bacteria were first noted. In this overview I trace the history of research on mitochondria from that time up to the present day, focussing on the major milestones of the overlapping eras of mitochondrial biochemistry, genetics, pathology and cell biology, and its explosion into new areas in the past 25 years. Nowadays, mitochondria are considered to be fully integrated into cell physiology, rather than serving specific functions in isolation.

Slow crabs - fast genomes: Locomotory capacity predicts skew magnitude in crustacean mitogenomes

Base composition skews (G-C/G+C) of mitochondrial genomes are believed to be primarily driven by mutational pressure, which is positively correlated with metabolic rate. In marine animals, metabolic rate is also positively correlated with locomotory capacity. Given the central role of mitochondria in energy metabolism, we hypothesised that selection for locomotory capacity should be positively correlated with the strength of purifying selection (dN/dS), and thus be negatively correlated with the skew magnitude. Therefore, these two models assume diametrically opposite associations between the metabolic rate and skew magnitude: positive correlation in the prevailing paradigm, and negative in our working hypothesis. We examined correlations between the skew magnitude, metabolic rate, locomotory capacity, and several other variables previously associated with mitochondrial evolution on 287 crustacean mitogenomes. Weakly locomotory taxa had higher skew magnitude and ω (dN/dS) values, but not the gene order rearrangement rate. Skew and ω magnitudes were correlated. Multilevel regression analyses indicated that three competing variables, body size, gene order rearrangement rate, and effective population size, had negligible impacts on the skew magnitude. In most crustacean lineages selection for locomotory capacity appears to be the primary factor determining the skew magnitude. Contrary to the prevailing paradigm, this implies that adaptive selection outweighs nonadaptive selection (mutation pressure) in crustaceans. However, we found indications that effective population size (nonadaptive factor) may outweigh the impact of locomotory capacity in sessile crustaceans (Thecostraca). In conclusion, skew magnitude is a product of the interplay between adaptive and nonadaptive factors, the balance of which varies among lineages.

Evolutionary History of Inversions in Directional Mutational Pressures in Crustacean Mitochondrial Genomes: Implications for Evolutionary Studies

Inversions of the origin of replication (ORI) in mitochondrial genomes produce asymmetrical mutational pressures that can cause strong base composition skews. Due to skews often being overlooked, the total number of crustacean lineages that underwent ORI events remains unknown. We analysed skews, cumulative skew plots, conserved sequence motifs, and mitochondrial architecture of all 965 available crustacean mitogenomes (699 unique species). We found indications of an ORI in 159 (22.7%) species, and mapped these to 23 ORI events: 16 identified with confidence and 7 putative (13 newly proposed, and for 5 we improved the resolution). Two ORIs occurred at or above the order level: Isopoda and Copepoda. Shifts in skew plots are not a precise tool for identifying the replication mechanism. We discuss how ORIs can produce mutational bursts in mitogenomes and show how these can interfere with various types of evolutionary studies. Phylogenetic analyses were plagued by artefactual clustering, and ORI lineages exhibited longer branches, a higher number of synonymous substitutions, higher mutational saturation, and higher compositional heterogeneity. ORI events also affected codon usage and protein properties. We discuss how this may have caused erroneous interpretation of data in previous studies that did not account for skew patterns.

Roles of the mitochondrial replisome in mitochondrial DNA deletion formation

Mitochondrial DNA (mtDNA) deletions are a common cause of human mitochondrial diseases. Mutations in the genes encoding components of the mitochondrial replisome, such as DNA polymerase gamma (Pol γ) and the mtDNA helicase Twinkle, have been associated with the accumulation of such deletions and the development of pathological conditions in humans. Recently, we demonstrated that changes in the level of wild-type Twinkle promote mtDNA deletions, which implies that not only mutations in, but also dysregulation of the stoichiometry between the replisome components is potentially pathogenic. The mechanism(s) by which alterations to the replisome function generate mtDNA deletions is(are) currently under debate. It is commonly accepted that stalling of the replication fork at sites likely to form secondary structures precedes the deletion formation. The secondary structural elements can be bypassed by the replication-slippage mechanism. Otherwise, stalling of the replication fork can generate single- and double-strand breaks, which can be repaired through recombination leading to the elimination of segments between the recombination sites. Here, we discuss aberrances of the replisome in the context of the two debated outcomes, and suggest new mechanistic explanations based on replication restart and template switching that could account for all the deletion types reported for patients.

Amino and carboxy-terminal extensions of yeast mitochondrial DNA polymerase assemble both the polymerization and exonuclease active sites

Human and yeast mitochondrial DNA polymerases (DNAPs), POLG and Mip1, are related by evolution to bacteriophage DNAPs. However, mitochondrial DNAPs contain unique amino and carboxyl-terminal extensions that physically interact. Here we describe that N-terminal deletions in Mip1 polymerases abolish polymerization and decrease exonucleolytic degradation, whereas moderate C-terminal deletions reduce polymerization. Similarly, to the N-terminal deletions, an extended C-terminal deletion of 298 amino acids is deficient in nucleotide addition and exonucleolytic degradation of double and single-stranded DNA. The latter observation suggests that the physical interaction between the amino and carboxyl-terminal regions of Mip1 may be related to the spread of pathogenic POLG mutant along its primary sequence.

Characterization of the human homozygous R182W POLG2 mutation in mitochondrial DNA depletion syndrome

Mutations in mitochondrial DNA (mtDNA) have been linked to a variety of metabolic, neurological and muscular diseases which can present at any time throughout life. MtDNA is replicated by DNA polymerase gamma (Pol γ), twinkle helicase and mitochondrial single-stranded binding protein (mtSSB). The Pol γ holoenzyme is a heterotrimer consisting of the p140 catalytic subunit and a p55 homodimeric accessory subunit encoded by the nuclear genes POLG and POLG2, respectively. The accessory subunits enhance DNA binding and promote processive DNA synthesis of the holoenzyme. Mutations in either POLG or POLG2 are linked to disease and adversely affect maintenance of the mitochondrial genome, resulting in depletion, deletions and/or point mutations in mtDNA. A homozygous mutation located at Chr17: 62492543G>A in POLG2, resulting in R182W substitution in p55, was previously identified to cause mtDNA depletion and fatal hepatic liver failure. Here we characterize this homozygous R182W p55 mutation using in vivo cultured cell models and in vitro biochemical assessments. Compared to control fibroblasts, homozygous R182W p55 primary dermal fibroblasts exhibit a two-fold slower doubling time, reduced mtDNA copy number and reduced levels of POLG and POLG2 transcripts correlating with the reported disease state. Expression of R182W p55 in HEK293 cells impairs oxidative-phosphorylation. Biochemically, R182W p55 displays DNA binding and association with p140 similar to WT p55. R182W p55 mimics the ability of WT p55 to stimulate primer extension, support steady-state nucleotide incorporation, and suppress the exonuclease function of Pol γin vitro. However, R182W p55 has severe defects in protein stability as determined by differential scanning fluorimetry and in stimulating function as determined by thermal inactivation. These data demonstrate that the Chr17: 62492543G>A mutation in POLG2, R182W p55, severely impairs stability of the accessory subunit and is the likely cause of the disease phenotype.

Nuclear genes involved in mitochondrial diseases caused by instability of mitochondrial DNA

Mitochondrial diseases are defined by a respiratory chain dysfunction and in most of the cases manifest as multisystem disorders with predominant expression in muscles and nerves and may be caused by mutations in mitochondrial (mtDNA) or nuclear (nDNA) genomes. Most of the proteins involved in respiratory chain function are nuclear encoded, although 13 subunits of respiratory chain complexes (together with 2 rRNAs and 22 tRNAs necessary for their translation) encoded by mtDNA are essential for cell function. nDNA encodes not only respiratory chain subunits but also all the proteins responsible for mtDNA maintenance, especially those involved in replication, as well as other proteins necessary for the transcription and copy number control of this multicopy genome. Mutations in these genes can cause secondary instability of the mitochondrial genome in the form of depletion (decreased number of mtDNA molecules in the cell), vast multiple deletions or accumulation of point mutations which in turn leads to mitochondrial diseases inherited in a Mendelian fashion. The list of genes involved in mitochondrial DNA maintenance is long, and still incomplete.

Plant organellar DNA polymerases are replicative and translesion DNA synthesis polymerases

Genomes acquire lesions that can block the replication fork and some lesions must be bypassed to allow survival. The nuclear genome of flowering plants encodes two family-A DNA polymerases (DNAPs), the result of a duplication event, that are the sole DNAPs in plant organelles. These DNAPs, dubbed Plant Organellar Polymerases (POPs), resemble the Klenow fragment of bacterial DNAP I and are not related to metazoan and fungal mitochondrial DNAPs. Herein we report that replicative POPs from the plant model Arabidopsis thaliana (AtPolI) efficiently bypass one the most insidious DNA lesions, an apurinic/apyrimidinic (AP) site. AtPolIs accomplish lesion bypass with high catalytic efficiency during nucleotide insertion and extension. Lesion bypass depends on two unique polymerization domain insertions evolutionarily unrelated to the insertions responsible for lesion bypass by DNAP θ, an analogous lesion bypass polymerase. AtPolIs exhibit an insertion fidelity that ranks between the fidelity of replicative and lesion bypass DNAPs, moderate 3′-5′ exonuclease activity and strong strand-displacement. AtPolIs are the first known example of a family-A DNAP evolved to function in both DNA replication and lesion bypass. The lesion bypass capabilities of POPs may be required to prevent replication fork collapse in plant organelles.

Characterization of the complete mitochondrial genome of Ortleppascaris sinensis (Nematoda: Heterocheilidae) and comparative mitogenomic analysis of eighteen Ascaridida nematodes

Ortleppascaris sinensis (Nematoda: Ascaridida) is a dominant intestinal nematode of the captive Chinese alligator. However, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. In this study, the complete mitochondrial (mt) genome sequence of O. sinensis was first determined using a polymerase chain reaction (PCR)-based primer-walking strategy, and this is also the first sequencing of the complete mitochondrial genome of a member of the genus Ortleppascaris. The circular mitochondrial genome (13,828 bp) of O. sinensis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes, but lacked the ATP synthetase subunit 8 gene. Finally, phylogenetic analysis of mtDNAs indicated that the genus Ortleppascaris should be attributed to the family Heterocheilidae. It is necessary to sequence more mtNDAs of Ortleppascaris nematodes in the future to test and confirm our conclusion. The complete mitochondrial genome sequence of O. sinensis reported here should contribute to molecular diagnosis, epidemiological investigations and ecological studies of O. sinensis and other related Ascaridida nematodes.

Whole exome sequencing identifies a homozygous POLG2 missense variant in an infant with fulminant hepatic failure and mitochondrial DNA depletion

Mitochondrial DNA (mtDNA) depletion syndrome manifests as diverse early-onset diseases that affect skeletal muscle, brain and liver function. Mutations in several nuclear DNA-encoded genes cause mtDNA depletion. We report on a patient, a 3-month-old boy who presented with hepatic failure, and was found to have severe mtDNA depletion in liver and muscle. Whole-exome sequencing identified a homozygous missense variant (c.544C > T, p.R182W) in the accessory subunit of mitochondrial DNA polymerase gamma (POLG2), which is required for mitochondrial DNA replication. This variant is predicted to disrupt a critical region needed for homodimerization of the POLG2 protein and cause loss of processive DNA synthesis. Both parents were phenotypically normal and heterozygous for this variant. Heterozygous mutations in POLG2 were previously associated with progressive external ophthalmoplegia and mtDNA deletions. This is the first report of a patient with a homozygous mutation in POLG2 and with a clinical presentation of severe hepatic failure and mitochondrial depletion.

Animal Mitochondrial DNA Replication

Recent advances in the field of mitochondrial DNA (mtDNA) replication highlight the diversity of both the mechanisms utilized and the structural and functional organization of the proteins at mtDNA replication fork, despite the relative simplicity of the animal mtDNA genome. DNA polymerase γ, mtDNA helicase and mitochondrial single-stranded DNA-binding protein-the key replisome proteins, have evolved distinct structural features and biochemical properties. These appear to be correlated with mtDNA genomic features in different metazoan taxa and with their modes of DNA replication, although substantial integrative research is warranted to establish firmly these links. To date, several modes of mtDNA replication have been described for animals: rolling circle, theta, strand-displacement, and RITOLS/bootlace. Resolution of a continuing controversy relevant to mtDNA replication in mammals/vertebrates will have a direct impact on the mechanistic interpretation of mtDNA-related human diseases. Here we review these subjects, integrating earlier and recent data to provide a perspective on the major challenges for future research.

Mitochondrial DNA copy number and replication in reprogramming and differentiation

Until recently, it was thought that the role of the mitochondrial genome was confined to encoding key proteins that generate ATP through the process of oxidative phosphorylation in the electron transfer chain. However, with increasing new evidence, it is apparent that the mitochondrial genome has a major role to play in a number of diseases and phenotypes. For example, mitochondrial variants and copy number have been implicated in the processes of fertilisation outcome and development and the onset of tumorigenesis. On the other hand, mitochondrial DNA (mtDNA) haplotypes have been implicated in a variety of diseases and most likely account for the adaptation that our ancestors achieved in order that they were fit for their environments. The mechanisms, which enable the mitochondrial genome to either protect or promote the disease phenotype, require further elucidation. However, there appears to be significant 'crosstalk' between the chromosomal and mitochondrial genomes that enable this to take place. One such mechanism is the regulation of DNA methylation by mitochondrial DNA, which is often perturbed in reprogrammed cells that have undergone dedifferentiation and affects mitochondrial DNA copy number. Furthermore, it appears that the mitochondrial genome interacts with the chromosomal genome to regulate the transcription of key genes at certain stages during development. Additionally, the mitochondrial genome can accumulate a series of mtDNA variants, which can lead to diseases such as cancer. It is likely that a combination of certain mitochondrial variants and aberrant patterns of mtDNA copy number could indeed account for many diseases that have previously been unaccounted for. This review focuses on the role that the mitochondrial genome plays especially during early stages of development and in cancer.

Mitochondrial Single-stranded DNA-binding Proteins Stimulate the Activity of DNA Polymerase γ by Organization of the Template DNA

The activity of the mitochondrial replicase, DNA polymerase γ (Pol γ) is stimulated by another key component of the mitochondrial replisome, the mitochondrial single-stranded DNA-binding protein (mtSSB). We have performed a comparative analysis of the human and Drosophila Pols γ with their cognate mtSSBs, evaluating their functional relationships using a combined approach of biochemical assays and electron microscopy. We found that increasing concentrations of both mtSSBs led to the elimination of template secondary structure and gradual opening of the template DNA, through a series of visually similar template species. The stimulatory effect of mtSSB on Pol γ on these ssDNA templates is not species-specific. We observed that human mtSSB can be substituted by its Drosophila homologue, and vice versa, finding that a lower concentration of insect mtSSB promotes efficient stimulation of either Pol. Notably, distinct phases of the stimulation by both mtSSBs are distinguishable, and they are characterized by a similar organization of the template DNA for both Pols γ. We conclude that organization of the template DNA is the major factor contributing to the stimulation of Pol γ activity. Additionally, we observed that human Pol γ preferentially utilizes compacted templates, whereas the insect enzyme achieves its maximal activity on open templates, emphasizing the relative importance of template DNA organization in modulating Pol γ activity and the variation among systems.