Evolving mtDNA populations within cells

Oxidative stress and mitochondrial dysfunction of granulosa cells in polycystic ovarian syndrome

Polycystic ovarian syndrome (PCOS) is one of the leading causes of anovulatory infertility in women, affecting 5%-15% of women of reproductive age worldwide. The clinical manifestations of patients include ovulation disorders, amenorrhea, hirsutism, and obesity. Life-threatening diseases, such as endometrial cancer, type 2 diabetes, hyperlipidaemia, hypertension, and cardiovascular disease, can be distant complications of PCOS. PCOS has diverse etiologies and oxidative stress (OS) plays an important role. Mitochondria, as the core organelles of energy production, are the main source of reactive oxygen species (ROS). The process of follicular growth and development is extremely complex, and the granulosa cells (GCs) are inextricably linked to follicular development. The abnormal function of GCs may directly affect follicular development and alter many symptoms of PCOS. Significantly higher levels of OS markers and abnormal mitochondrial function in GCs have been found in patients with PCOS compared to healthy subjects, suggesting that increased OS is associated with PCOS progression. Therefore, the aim of this review was to summarize and discuss the findings suggesting that OS and mitochondrial dysfunction in GCs impair ovarian function and induce PCOS.

Cellular and environmental dynamics influence species-specific extents of organelle gene retention

Mitochondria and plastids rely on many nuclear-encoded genes, but retain small subsets of the genes they need to function in their own organelle DNA (oDNA). Different species retain different numbers of oDNA genes, and the reasons for these differences are not completely understood. Here, we use a mathematical model to explore the hypothesis that the energetic demands imposed by an organism's changing environment influence how many oDNA genes it retains. The model couples the physical biology of cell processes of gene expression and transport to a supply-and-demand model for the environmental dynamics to which an organism is exposed. The trade-off between fulfilling metabolic and bioenergetic environmental demands, and retaining genetic integrity, is quantified for a generic gene encoded either in oDNA or in nuclear DNA. Species in environments with high-amplitude, intermediate-frequency oscillations are predicted to retain the most organelle genes, whereas those in less dynamic or noisy environments the fewest. We discuss support for, and insight from, these predictions with oDNA data across eukaryotic taxa, including high oDNA gene counts in sessile organisms exposed to day-night and intertidal oscillations (including plants and algae) and low counts in parasites and fungi.

Evolutionary inference across eukaryotes identifies universal features shaping organelle gene retention

Mitochondria and plastids power complex life. Why some genes and not others are retained in their organelle DNA (oDNA) genomes remains a debated question. Here, we attempt to identify the properties of genes and associated underlying mechanisms that determine oDNA retention. We harness over 15k oDNA sequences and over 300 whole genome sequences across eukaryotes with tools from structural biology, bioinformatics, machine learning, and Bayesian model selection. Previously hypothesized features, including the hydrophobicity of a protein product, and less well-known features, including binding energy centrality within a protein complex, predict oDNA retention across eukaryotes, with additional influences of nucleic acid and amino acid biochemistry. Notably, the same features predict retention in both organelles, and retention models learned from one organelle type quantitatively predict retention in the other, supporting the universality of these features-which also distinguish gene profiles in more recent, independent endosymbiotic relationships. A record of this paper's transparent peer review process is included in the supplemental information.

Organelle bottlenecks facilitate evolvability by traversing heteroplasmic fitness valleys

Bioenergetic organelles—mitochondria and plastids—retain their own genomes (mtDNA and ptDNA), and these organelle DNA (oDNA) molecules are vital for eukaryotic life. Like all genomes, oDNA must be able to evolve to suit new environmental challenges. However, mixed oDNA populations in cells can challenge cellular bioenergetics, providing a penalty to the appearance and adaptation of new mutations. Here we show that organelle “bottlenecks,” mechanisms increasing cell-to-cell oDNA variability during development, can overcome this mixture penalty and facilitate the adaptation of beneficial mutations. We show that oDNA heteroplasmy and bottlenecks naturally emerge in evolutionary simulations subjected to fluctuating environments, demonstrating that this evolvability is itself evolvable. Usually thought of as a mechanism to clear damaging mutations, organelle bottlenecks therefore also resolve the tension between intracellular selection for pure cellular oDNA populations and the “bet-hedging” need for evolvability and adaptation to new environments. This general theory suggests a reason for the maintenance of organelle heteroplasmy in cells, and may explain some of the observed diversity in organelle maintenance and inheritance across taxa.

Mitochondrial DNA variants segregate during human preimplantation development into genetically different cell lineages that are maintained postnatally

Humans present remarkable diversity in their mitochondrial DNA (mtDNA) in terms of variants across individuals as well as across tissues and even cells within one person. We have investigated the timing of the first appearance of this variant-driven mosaicism. For this, we deep-sequenced the mtDNA of 254 oocytes from 85 donors, 158 single blastomeres of 25 day-3 embryos, 17 inner cell mass and trophectoderm samples of 7 day-5 blastocysts, 142 bulk DNA and 68 single cells of different adult tissues. We found that day-3 embryos present blastomeres that carry variants only detected in that cell, showing that mtDNA mosaicism arises very early in human development. We classified the mtDNA variants based on their recurrence or uniqueness across different samples. Recurring variants had higher heteroplasmic loads and more frequently resulted in synonymous changes or were located in non-coding regions than variants unique to one oocyte or single embryonic cell. These differences were maintained through development, suggesting that the mtDNA mosaicism arising in the embryo is maintained into adulthood. We observed a decline in potentially pathogenic variants between day-3 and day-5 of development, suggesting early selection. We propose a model in which closely clustered mitochondria carrying specific mtDNA variants in the ooplasm are asymmetrically distributed throughout the cell divisions of the preimplantation embryo, resulting in the earliest form of mtDNA mosaicism in human development.

Single-cell mtDNA heteroplasmy in colorectal cancer

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

Cytoplasmic Transfer Methods for Studying the Segregation of Mitochondrial DNA in Mice

Heteroplasmic mice represent a valuable tool to study the segregation of different mtDNA haplotypes (mtDNAs with differing alleles) in vivo against a defined nuclear background. We describe two methods for the creation of such models, differing in the resulting initial heteroplasmy levels: (a) transfer of ooplasm and (b) fusion of two blastomeres. These methods result in typical heteroplasmy of 5% and 50% donor mtDNA , respectively. The choice of method depends on the aim of the study. By means of breeding even 100% donor mtDNA can be reached within a few generations.

Defining coalescent genes: Theory meets practice in organelle phylogenomics

The species tree paradigm that dominates current molecular systematic practice infers species trees from collections of sequences under assumptions of the multispecies coalescent (MSC), i.e., that there is free recombination between the sequences and no (or very low) recombination within them. These coalescent genes (c-genes) are thus defined in an historical rather than molecular sense, and can in theory be as large as an entire genome or as small as a single nucleotide. A debate about how to define c-genes centers on the contention that nuclear gene sequences used in many coalescent analyses undergo too much recombination, such that their introns comprise multiple c-genes, violating a key assumption of the MSC. Recently a similar argument has been made for the genes of plastid (e.g., chloroplast) and mitochondrial genomes, which for the last 30 or more years have been considered to represent a single c-gene for the purposes of phylogeny reconstruction because they are non-recombining in a historical sense. Consequently, it has been suggested that these genomes should be analyzed using coalescent methods that treat their genes-over 70 protein-coding genes in the case of most plastid genomes (plastomes)-as independent estimates of species phylogeny, in contrast to the usual practice of concatenation, which is appropriate for generating gene trees. However, although recombination certainly occurs in the plastome, as has been recognized since the 1970's, it is unlikely to be phylogenetically relevant. This is because such historically effective recombination can only occur when plastomes with incongruent histories are brought together in the same plastid. However, plastids sort rapidly into different cell lineages and rarely fuse. Thus, because of plastid biology, the plastome is a more canonical c-gene than is the average multi-intron mammalian nuclear gene. The plastome should thus continue to be treated as a single estimate of the underlying species phylogeny, as should the mitochondrial genome. The implications of this long-held insight of molecular systematics for studies in the phylogenomic era are explored.

Sexually antagonistic evolution of mitochondrial and nuclear linkage

Across eukaryotes, genes encoding bioenergetic machinery are located in both mitochondrial and nuclear DNA, and incompatibilities between the two genomes can be devastating. Mitochondria are often inherited maternally, and theory predicts sex-specific fitness effects of mitochondrial mutational diversity. Yet how evolution acts on linkage patterns between mitochondrial and nuclear genomes is poorly understood. Using novel mito-nuclear population-genetic models, we show that the interplay between nuclear and mitochondrial genes maintains mitochondrial haplotype diversity within populations, and selects both for sex-independent segregation of mitochondrion-interacting genes and for paternal leakage. These effects of genetic linkage evolution can eliminate male-harming fitness effects of mtDNA mutational diversity. With maternal mitochondrial inheritance, females maintain a tight mitochondrial-nuclear match, but males accumulate mismatch mutations because of the weak statistical associations between the two genomic components. Sex-independent segregation of mitochondria-interacting loci improves the mito-nuclear match. In a sexually antagonistic evolutionary process, male nuclear alleles evolve to increase the rate of recombination, whereas females evolve to suppress it. Paternal leakage of mitochondria can evolve as an alternative mechanism to improve the mito-nuclear linkage. Our modelling framework provides an evolutionary explanation for the observed paucity of mitochondrion-interacting genes on mammalian sex chromosomes and for paternal leakage in protists, plants, fungi and some animals.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Avoiding organelle mutational meltdown across eukaryotes with or without a germline bottleneck

Mitochondrial DNA (mtDNA) and plastid DNA (ptDNA) encode vital bioenergetic apparatus, and mutations in these organelle DNA (oDNA) molecules can be devastating. In the germline of several animals, a genetic “bottleneck” increases cell-to-cell variance in mtDNA heteroplasmy, allowing purifying selection to act to maintain low proportions of mutant mtDNA. However, most eukaryotes do not sequester a germline early in development, and even the animal bottleneck remains poorly understood. How then do eukaryotic organelles avoid Muller’s ratchet—the gradual buildup of deleterious oDNA mutations? Here, we construct a comprehensive and predictive genetic model, quantitatively describing how different mechanisms segregate and decrease oDNA damage across eukaryotes. We apply this comprehensive theory to characterise the animal bottleneck with recent single-cell observations in diverse mouse models. Further, we show that gene conversion is a particularly powerful mechanism to increase beneficial cell-to-cell variance without depleting oDNA copy number, explaining the benefit of observed oDNA recombination in diverse organisms which do not sequester animal-like germlines (for example, sponges, corals, fungi, and plants). Genomic, transcriptomic, and structural datasets across eukaryotes support this mechanism for generating beneficial variance without a germline bottleneck. This framework explains puzzling oDNA differences across taxa, suggesting how Muller’s ratchet is avoided in different eukaryotes.

A comprehensive model for mitochondrial and plasmid DNA segregation, supported by with genomic, transcriptomic, and single-cell data, shows how the attritional effects of Muller’s ratchet can be avoided in the organelles of diverse eukaryotes.

Beyond the "selfish mitochondrion" theory of uniparental inheritance: A unified theory based on mutational variance redistribution

"Selfish" gene theories have offered invaluable insight into eukaryotic genome evolution, but they can also be misleading. The "selfish mitochondrion" hypothesis, developed in the 90s explained uniparental organelle inheritance as a mechanism of conflict resolution, improving cooperation between genetically distinct compartments of the cell. But modern population genetic models provided a more general explanation for uniparental inheritance based on mutational variance redistribution, modulating the efficiency of both purifying and adaptive selection. Nevertheless, as reviewed here, "selfish" conflict theories still dominate the literature. While these hypotheses are rich in metaphor and highly intuitive, selective focus on only one type of mitochondrial mutation limits the generality of our understanding and hinders progress in mito-nuclear evolution theory. Recognizing that uniparental inheritance may have evolved-and is maintained across the eukaryotic tree of life-because of its influence on mutational variance and improved selection will only increase the generality of our evolutionary reasoning, retaining "selfish" conflict explanations as a special case of a much broader theory.

A Naturally Heteroplasmic Clam Provides Clues about the Effects of Genetic Bottleneck on Paternal mtDNA

Mitochondrial DNA (mtDNA) is present in multiple copies within an organism. Since these copies are not identical, a single individual carries a heterogeneous population of mtDNAs, a condition known as heteroplasmy. Several factors play a role in the dynamics of the within-organism mtDNA population: among them, genetic bottlenecks, selection, and strictly maternal inheritance are known to shape the levels of heteroplasmy across mtDNAs.

In Metazoa, the only evolutionarily stable exception to the strictly maternal inheritance of mitochondria is the doubly uniparental inheritance (DUI), reported in 100+ bivalve species. In DUI species, there are two highly divergent mtDNA lineages, one inherited through oocyte mitochondria (F-type) and the other through sperm mitochondria (M-type). Having both parents contributing to the mtDNA pool of the progeny makes DUI a unique system to study the dynamics of mtDNA populations. Since, in bivalves, the spermatozoon has few mitochondria (4–5), M-type mtDNA faces a tight bottleneck during embryo segregation, one of the narrowest mitochondrial bottlenecks investigated so far.

Here, we analyzed the F- and M-type mtDNA variability within individuals of the DUI species Ruditapes philippinarum and investigated for the first time the effects of such a narrow bottleneck affecting mtDNA populations. As a potential consequence of this narrow bottleneck, the M-type mtDNA shows a large variability in different tissues, a condition so pronounced that it leads to genotypes from different tissues of the same individual not to cluster together. We believe that such results may help understanding the effect of low population size on mtDNA bottleneck.

Alternative Applications of Genotyping Array Data Using Multivariant Methods

One of the forerunners that pioneered the revolution of high-throughput genomic technologies is the genotyping microarray technology, which can genotype millions of single-nucleotide variants simultaneously. Owing to apparent benefits, such as high speed, low cost, and high throughput, the genotyping array has gained lasting applications in genome-wide association studies (GWAS) and thus accumulated an enormous amount of data. Empowered by continuous manufactural upgrades and analytical innovation, unconventional applications of genotyping array data have emerged to address more diverse genetic problems, holding promise of boosting genetic research into human diseases through the re-mining of the rich accumulated data. Here, we review several unconventional genotyping array analysis techniques that have been built on the idea of large-scale multivariant analysis and provide empirical application examples. These unconventional outcomes of genotyping arrays include polygenic score, runs of homozygosity (ROH)/heterozygosity ratio, distant pedigree computation, and mitochondrial DNA (mtDNA) copy number inference.

MtDNA sequence features associated with 'selfish genomes' predict tissue-specific segregation and reversion

Mitochondrial DNA (mtDNA) encodes cellular machinery vital for cell and organism survival. Mutations, genetic manipulation, and gene therapies may produce cells where different types of mtDNA coexist in admixed populations. In these admixtures, one mtDNA type is often observed to proliferate over another, with different types dominating in different tissues. This ‘segregation bias’ is a long-standing biological mystery that may pose challenges to modern mtDNA disease therapies, leading to substantial recent attention in biological and medical circles. Here, we show how an mtDNA sequence’s balance between replication and transcription, corresponding to molecular ‘selfishness’, in conjunction with cellular selection, can potentially modulate segregation bias. We combine a new replication-transcription-selection (RTS) model with a meta-analysis of existing data to show that this simple theory predicts complex tissue-specific patterns of segregation in mouse experiments, and reversion in human stem cells. We propose the stability of G-quadruplexes in the mtDNA control region, influencing the balance between transcription and replication primer formation, as a potential molecular mechanism governing this balance. Linking mtDNA sequence features, through this molecular mechanism, to cellular population dynamics, we use sequence data to obtain and verify the sequence-specific predictions from this hypothesis on segregation behaviour in mouse and human mtDNA.

Mitochondrial dysfunction is involved in aristolochic acid I-induced apoptosis in renal proximal tubular epithelial cells

Aristolochic acid (AA) is a compound extracted from the Aristolochia species of herbs. AA exposure is associated with kidney injury known as aristolochic acid nephropathy (AAN). Proximal tubular epithelial cell (PTEC) is the primary target of AA and rich in mitochondria. Recently, increasing evidence suggests that mitochondrial dysfunction plays a critical role in the pathogenesis of kidney disease. However, the status of mitochondrial function in PTEC after exposure to AA remains largely unknown. The aim of this study was to explore the effect of aristolochic acid I (AAI) on cell apoptosis and mitochondrial function in PTEC. Normal rat kidney-52E (NRK-52E) cells were exposed to different concentrations of AAI for different time periods. Cell viability was detected by 3-(4,5-dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide, cell apoptosis was analyzed by flow cytometry, and the expression of cleaved caspase-3 by Western blotting. Mitochondrial function was evaluated by reactive oxygen species (ROS), mitochondrial membrane potential (MMP), mitochondrial DNA (mtDNA) copy number, and adenosine triphosphate (ATP). It was found that AAI reduced cell viability and increased cell apoptosis in a dose- and time-dependent manner. In parallel to increased apoptosis, NRK-52E cell manifested signs of mitochondrial dysfunction in response to AAI treatment. The data indicated that AAI could increase ROS level, lower MMP, decrease mtDNA copy number, and reduce ATP production. In addition, Szeto-Schiller 31, a mitochondria-targeted antioxidant peptide, attenuated AAI-induced mitochondrial dysfunction and apoptosis. Our study depicted significant aberrant of mitochondrial function in AAI-treated NRK-52E cell, which suggested that mitochondrial dysfunction may be involved in AAI-induced apoptosis in PTEC.

Varied Mechanisms and Models for the Varying Mitochondrial Bottleneck

Mitochondrial DNA (mtDNA) molecules exist in populations within cells, and may carry mutations. Different cells within an organism, and organisms within a family, may have different proportions of mutant mtDNA in these cellular populations. This diversity is often thought of as arising from a “genetic bottleneck.” This article surveys approaches to characterize and model the generation of this genetic diversity, aiming to provide an introduction to the range of concepts involved, and to highlight some recent advances in understanding. In particular, differences between the statistical “genetic bottleneck” (mutant proportion spread) and the physical mtDNA bottleneck and other cellular processes are highlighted. Particular attention is paid to the quantitative analysis of the “genetic bottleneck,” estimation of its magnitude from observed data, and inference of its underlying mechanisms. Evidence that the “genetic bottleneck” (mutant proportion spread) varies with age, between individuals and species, and across mtDNA sequences, is described. The interpretation issues that arise from sampling errors, selection, and different quantitative definitions are also discussed.