Dynamics of the most common pathogenic mtDNA variant m.3243A > G demonstrate frequency-dependency in blood and positive selection in the germline.

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate 'on average', creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an 'arching' profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Workflow Optimization for Identification of Female Germline or Oogonial Stem Cells in Human Ovarian Cortex Using Single-Cell RNA Sequence Analysis.

In 2004, the identification of female germline or oogonial stem cells (OSCs) that can support post-natal oogenesis in ovaries of adult mice sparked a major paradigm shift in reproductive biology. Although these findings have been independently verified, and further extended to include identification of OSCs in adult ovaries of many species ranging from pigs and cows to non-human primates and humans, a recent study rooted in single-cell RNA sequence analysis (scRNA-seq) of adult human ovarian cortical tissue claimed that OSCs do not exist, and that other groups working with OSCs following isolation by magnetic-assisted or fluorescence-activated cell sorting have mistaken perivascular cells (PVCs) for germ cells. Here we report that rare germ lineage cells with a gene expression profile matched to OSCs but distinct from that of other cells, including oocytes and PVCs, can be identified in adult human ovarian cortical tissue by scRNA-seq after optimization of analytical workflow parameters. Deeper cell-by-cell expression profiling also uncovered evidence of germ cells undergoing meiosis-I in adult human ovaries. Lastly, we show that, if not properly controlled for, PVCs can be inadvertently isolated during flow cytometry protocols designed to sort OSCs because of inherently high cellular autofluorescence. However, human PVCs and human germ cells segregate into distinct clusters following scRNA-seq due to non-overlapping gene expression profiles, which would preclude the mistaken identification and use of PVCs as OSCs during functional characterization studies.

3GOLD: optimized Levenshtein distance for clustering third-generation sequencing data.

BACKGROUND: Third-generation sequencing offers some advantages over next-generation sequencing predecessors, but with the caveat of harboring a much higher error rate. Clustering-related sequences is an essential task in modern biology. To accurately cluster sequences rich in errors, error type and frequency need to be accounted for. Levenshtein distance is a well-established mathematical algorithm for measuring the edit distance between words and can specifically weight insertions, deletions and substitutions. However, there are drawbacks to using Levenshtein distance in a biological context and hence has rarely been used for this purpose. We present novel modifications to the Levenshtein distance algorithm to optimize it for clustering error-rich biological sequencing data. RESULTS: We successfully introduced a bidirectional frameshift allowance with end-user determined accommodation caps combined with weighted error discrimination. Furthermore, our modifications dramatically improved the computational speed of Levenstein distance. For simulated ONT MinION and PacBio Sequel datasets, the average clustering sensitivity for 3GOLD was 41.45% (S.D. 10.39) higher than Sequence-Levenstein distance, 52.14% (S.D. 9.43) higher than Levenshtein distance, 55.93% (S.D. 8.67) higher than Starcode, 42.68% (S.D. 8.09) higher than CD-HIT-EST and 61.49% (S.D. 7.81) higher than DNACLUST. For biological ONT MinION data, 3GOLD clustering sensitivity was 27.99% higher than Sequence-Levenstein distance, 52.76% higher than Levenshtein distance, 56.39% higher than Starcode, 48% higher than CD-HIT-EST and 70.4% higher than DNACLUST. CONCLUSION: Our modifications to Levenshtein distance have improved its speed and accuracy compared to the classic Levenshtein distance, Sequence-Levenshtein distance and other commonly used clustering approaches on simulated and biological third-generation sequenced datasets. Our clustering approach is appropriate for datasets of unknown cluster centroids, such as those generated with unique molecular identifiers as well as known centroids such as barcoded datasets. A strength of our approach is high accuracy in resolving small clusters and mitigating the number of singletons.

Dynamics of WNT signaling components in the human ovary from development to adulthood.

WNT signaling has been shown to play a pivotal role in mammalian gonad development and sex differentiation; however, its role in the developing human ovary has not been investigated. We analyzed a quantitative mass spectrometry dataset to determine the expression of WNT signaling components between 47 and 137 days of development and in adult ovarian cortex tissue. WNT signaling was identified within the top ten canonical pathways of proteins detected at every developmental stage examined. We further examined the specific localization of WNT signaling components glycogen synthase kinase 3 (GSK3B), frizzled 2 (FZD2), and ß-catenin (CTNNB1) within ovarian tissue. GSK3B was nearly ubiquitously expressed during fetal development, while FZD2 was specific to germ cell nests during early development. ß-catenin exhibited translocation from primarily membrane bound during early ovarian development to cytoplasmic and nuclear staining specifically in early primordial follicles in the fetal ovary. This cytoplasmic and nuclear ß-catenin persisted in primordial follicles in adult ovarian tissue, but returned to membrane-bound localization in secondary follicles. We conclude that WNT signaling components are expressed in the human ovary from early to mid-gestation and remain in the adult ovary, and observed evidence for canonical WNT signaling only in the oocytes of primordial follicles. Together, these data are indicative of a role for canonical WNT signaling via ß-catenin nuclear translocation during human follicle formation and follicle maintenance.

Impact of exercise on oocyte quality in the POLG mitochondrial DNA mutator mouse.

The mtDNA 'mutator' mouse, also called the 'POLG' mouse, is a well-characterized model frequently used for studies of progeroid aging. Harboring a mutation in the proofreading domain of the mitochondrial polymerase, polymerase-γ (Polg), POLG mice acquire mtDNA mutations at an accelerated rate. This results in premature mitochondrial dysfunction and a systemic aging phenotype. Previous work has demonstrated that the progeroid phenotype in POLG is attenuated following endurance exercise, the only reported intervention to extend health span and lifespan of these mice. Herein, oocyte quality was evaluated in sedentary and exercised POLG mice. In mice homozygous for the Polg mutation, litter size is dramatically reduced as compared to heterozygous Polg mice. Following ovarian hyper-stimulation, oocytes were retrieved until 9 months of age in exercised and sedentary groups, with no oocytes ovulated thereafter. Although ovulated oocyte numbers were not impacted by exercise, we did find a modest improvement in both the ovarian follicle reserve and in oocyte quality based on meiotic spindle assembly, chromosomal segregation and mitochondrial distribution at 7 months of age in exercised POLG mice as compared to sedentary counterparts. Of note, analysis of mtDNA mutational load revealed no differences between exercised and sedentary groups. Collectively, these data indicate that exercise differentially influences somatic tissues of the POLG mouse as compared to oocytes, highlighting important mechanistic differences between mitochondrial regulatory mechanisms in the soma and the germline.

A role for retinoids in human oocyte fertilization: regulation of connexin 43 by retinoic acid in cumulus granulosa cells.

Retinoids are essential for ovarian steroid production and oocyte maturation in mammals. Oocyte competency is known to positively correlate with efficient gap junction intercellular communication (GJIC) among granulosa cells in the cumulus-oocyte complex. Connexin 43 (C x 43) is the main subunit of gap junction channels in human cumulus granulosa cells (CGC) and is regulated by all-trans retinoic acid (ATRA) in other hormone responsive cell types. The objectives of this study were to quantify retinoid levels in human CGC obtained during IVF oocyte retrievals, to investigate the potential relationship between CGC ATRA levels and successful oocyte fertilization, and to determine the effects of ATRA on C x 43 protein expression in CGC. Results showed that CGC cultures actively metabolize retinol to produce ATRA. Grouped according to fertilization rate tertiles, mean ATRA levels were 2-fold higher in pooled CGC from women in the highest versus the lowest tertile (P < 0.05). ATRA induced a rapid dephosphorylation of C x 43 in CGC and granulosa cell line (KGN) cultures resulting in a >2-fold increase in the expression of the functional non-phosphorylated (P0) species (P < 0.02). Similar enhancement of P0 by ATRA was shown in CGC and KGN cultures co-treated with LH or hCG which, by themselves, enhanced the protein levels of C x 43 without altering its phosphorylation profile. Correspondingly, the combination of ATRA+hCG treatment of KGN caused a significant increase in GJIC compared with single agent treatments (P < 0.025) and a doubling of GJIC from that seen in untreated cells (P < 0.01). These findings indicate that CGC are a primary site of retinoid uptake and ATRA biosynthesis. Regulation of C x 43 by ATRA may serve an important role in folliculogenesis, development of oocyte competency, and successful fertilization by increasing GJIC in CGC.

Reanalysis of mtDNA mutations of human primordial germ cells (PGCs) reveals NUMT contamination and suggests that selection in PGCs may be positive.

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology 20: 144-51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between early-stage and late-stage PGCs. We re-analyzed Floros' et al. data and determined that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly, when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to bulk PGCs, Floros et al. reported the analysis of single-cell late-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, there were no NUMT-derived SNVs among single PGC mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. More specifically, nonsynonymous mutations show faster intracellular genetic drift towards higher mutant fraction than synonymous ones. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously andpotentially may be evolutionarily advantageous.

Broad proteomics analysis of seeding-induced aggregation of α-synuclein in M83 neurons reveals remodeling of proteostasis mechanisms that might contribute to Parkinson's disease pathogenesis.

Aggregation of misfolded α-synuclein (α-syn) is a key characteristic feature of Parkinson's disease (PD) and related synucleinopathies. The nature of these aggregates and their contribution to cellular dysfunction is still not clearly elucidated. We employed mass spectrometry-based total and phospho-proteomics to characterize the underlying molecular and biological changes due to α-syn aggregation using the M83 mouse primary neuronal model of PD. We identified gross changes in the proteome that coincided with the formation of large Lewy body-like α-syn aggregates in these neurons. We used protein-protein interaction (PPI)-based network analysis to identify key protein clusters modulating specific biological pathways that may be dysregulated and identified several mechanisms that regulate protein homeostasis (proteostasis). The observed changes in the proteome may include both homeostatic compensation and dysregulation due to α-syn aggregation and a greater understanding of both processes and their role in α-syn-related proteostasis may lead to improved therapeutic options for patients with PD and related disorders.

Assaying Mitochondrial Function by Multiparametric Flow Cytometry.

Flow cytometry has been a vital tool in cell biology for decades based on its versatile ability to detect and quantifiably measure both physical and chemical attributes of individual cells within a larger population. More recently, advances in flow cytometry have enabled nanoparticle detection. This is particularly applicable to mitochondria, which, as intracellular organelles have distinct subpopulations that can be evaluated based on differences in functional, physical, and chemical attributes, in a manner analogous to cells. This includes distinctions based on size, mitochondrial membrane potential (ΔΨm), chemical properties, and protein expression on the outer mitochondrial membrane in intact, functional organelles and internally in fixed samples. This method allows for multiparametric analysis of subpopulations of mitochondria, as well as collection for downstream analysis down to the level of a single organelle. The present protocol describes a framework for analysis and sorting mitochondria by flow cytometry, termed fluorescence activated mitochondrial sorting (FAMS), based on the separation of individual mitochondria belonging to subpopulations of interest using fluorescent dyes and antibody labeling.

Interpreting Sequence-Levenshtein distance for determining error type and frequency between two embedded sequences of equal length.

Levenshtein distance is a commonly used edit distance metric, typically applied in language processing, and to a lesser extent, in molecular biology analysis. Biological nucleic acid sequences are often embedded in longer sequences and are subject to insertion and deletion errors that introduce frameshift during sequencing. These frameshift errors are due to string context and should not be counted as true biological errors. Sequence-Levenshtein distance is a modification to Levenshtein distance that is permissive of frameshift error without additional penalty. However, in a biological context Levenshtein distance needs to accommodate both frameshift and weighted errors, which Sequence-Levenshtein distance cannot do. Errors are weighted when they are associated with a numerical cost that corresponds to their frequency of appearance. Here, we describe a modification that allows the use of Levenshtein distance and Sequence-Levenshtein distance to appropriately accommodate penalty-free frameshift between embedded sequences and correctly weight specific error types.

Lineage-Mismatched Mitochondrial Replacement in an Inducible Mitochondrial Depletion Model Effectively Restores the Original Proteomic Landscape of Recipient Cells.

In addition to critical roles in bioenergetics, mitochondria are key contributors to the regulation of many other functions in cells, ranging from steroidogenesis to apoptosis. Numerous studies further demonstrate that cell type-specific differences exist in mitochondria, with cells of a given lineage tailoring their endogenous mitochondrial population to suit specific functional needs. These findings, coupled with studies of the therapeutic potential of mitochondrial transplantation, provide a strong impetus to better understand how mitochondria can influence cell function or fate. Here an inducible mitochondrial depletion modelis used to study how cells lacking endogenous mitochondria respond, on a global protein expression level, to transplantation with lineage-mismatched (LM) mitochondria. It is shown that LM mitochondrial transplantation does not alter the proteomic profile in nonmitochondria-depleted recipient cells; however, enforced depletion of endogenous mitochondria results in dramatic changes in the proteomic landscape, which returns to the predepletion state following internalization of LM mitochondria. These data, derived from a cell system that can be rendered free of influence by endogenous mitochondria, indicate that transplantation of mitochondria-even from a source that differs significantly from the recipient cell population, effectively restores a normal proteomic landscape to cells lacking their own mitochondria.

Characterization of Oogonial Stem Cells in Adult Mouse Ovaries with Age and Comparison to In Silico Data on Human Ovarian Aging.

Many adult somatic stem cell lineages are comprised of subpopulations that differ in gene expression, mitotic activity, and differentiation status. In this study, we explored if cellular heterogeneity also exists within oogonial stem cells (OSCs), and how chronological aging impacts OSCs. In OSCs isolated from mouse ovaries by flow cytometry and established in culture, we identified subpopulations of OSCs that could be separated based on differential expression of stage-specific embryonic antigen 1 (SSEA1) and cluster of differentiation 61 (CD61). Levels of aldehyde dehydrogenase (ALDH) activity were inversely related to OSC differentiation, whereas commitment of OSCs to differentiation through transcriptional activation of stimulated by retinoic acid gene 8 was marked by a decline in ALDH activity and in SSEA1 expression. Analysis of OSCs freshly isolated from ovaries of mice between 3 and 20 months of age revealed that these subpopulations were present and persisted throughout adult life. However, expression of developmental pluripotency associated 3 (Dppa3), an epigenetic modifier that promotes OSC differentiation into oocytes, was lost as the mice transitioned from a time of reproductive compromise (10 months) to reproductive failure (15 months). Further analysis showed that OSCs from aged females could be established in culture, and that once established the cultured cells reactivated Dppa3 expression and the capacity for oogenesis. Analysis of single-nucleus RNA sequence data sets generated from ovaries of women in their 20s versus those in their late 40s to early 50s showed that the frequency of DPPA3-expressing cells decreased with advancing age, and this was paralleled by reduced expression of several key meiotic differentiation genes. These data support the existence of OSC subpopulations that differ in gene expression profiles and differentiation status. In addition, an age-related decrease in Dppa3/DPPA3 expression, which is conserved between mice and humans, may play a role in loss of the ability of OSCs to maintain oogenesis with age.

A transgenic zebrafish model of targeted oocyte ablation and de novo oogenesis.

We describe here a novel transgenic zebrafish, Tg(zpc:G4VP16/UAS:nfsB-mCherry) that effectively demonstrates the targeted oocyte ablation in the adult zebrafish ovary. This transgenic line expresses bacterial nitroreductase enzyme (nfsB) under the control of the oocyte-specific zona pellucida C (zpc) gene promoter. Adult transgenic females exposed to the prodrug metronidazole demonstrated near-complete ablation of growing oocytes, resulting in ovarian degeneration and complete cessation of reproductive function. Within 4 weeks of prodrug removal, treated fish demonstrated complete anatomical regeneration of the ovary and, within 7 weeks, ovarian function (fertility) was fully restored. Together, these results demonstrate functional renewal of the oocyte pool in the adult zebrafish ovary. Accordingly, this transgenic zebrafish model system provides a novel means to investigate ovarian growth dynamics in a genetically tractable vertebrate, and may be useful for evaluating signaling interactions that regulate gonadal development processes such as de novo oogenesis.

Non-neutral clonal selection and its potential role in mammalian germline stem cell dysfunction with advancing age.

The concept of natural selection, or "survival of the fittest", refers to an evolutionary process in nature whereby traits emerge in individuals of a population through random gene alterations that enable those individuals to better adapt to changing environmental conditions. This genetic variance allows certain members of the population to gain an advantage over others in the same population to survive and reproduce in greater numbers under new environmental pressures, with the perpetuation of those advantageous traits in future progeny. Here we present that the behavior of adult stem cells in a tissue over time can, in many respects, be viewed in the same manner as evolution, with each stem cell clone being representative of an individual within a population. As stem cells divide or are subjected to cumulative oxidative damage over the lifespan of the organism, random genetic alterations are introduced into each clone that create variance in the population. These changes may occur in parallel to, or in response to, aging-associated changes in microenvironmental cues perceived by the stem cell population. While many of these alterations will be neutral or silent in terms of affecting cell function, a small fraction of these changes will enable certain clones to respond differently to shifts in microenvironmental conditions that arise with advancing age. In some cases, the same advantageous genetic changes that support survival and expansion of certain clones over others in the population (viz. non-neutral competition) could be detrimental to the downstream function of the differentiated stem cell descendants. In the context of the germline, such a situation would be devastating to successful propagation of the species across generations. However, even within a single generation, the "evolution" of stem cell lineages in the body over time can manifest into aging-related organ dysfunction and failure, as well as lead to chronic inflammation, hyperplasia, and cancer. Increased research efforts to evaluate stem cells within a population as individual entities will improve our understanding of how organisms age and how certain diseases develop, which in turn may open new opportunities for clinical detection and management of diverse pathologies.

A FoxA2+ long-term stem cell population is necessary for growth plate cartilage regeneration after injury.

Longitudinal bone growth, achieved through endochondral ossification, is accomplished by a cartilaginous structure, the physis or growth plate, comprised of morphologically distinct zones related to chondrocyte function: resting, proliferating and hypertrophic zones. The resting zone is a stem cell-rich region that gives rise to the growth plate, and exhibits regenerative capabilities in response to injury. We discovered a FoxA2+group of long-term skeletal stem cells, situated at the top of resting zone, adjacent the secondary ossification center, distinct from the previously characterized PTHrP+ stem cells. Compared to PTHrP+ cells, FoxA2+ cells exhibit higher clonogenicity and longevity. FoxA2+ cells exhibit dual osteo-chondro-progenitor activity during early postnatal development (P0-P28) and chondrogenic potential beyond P28. When the growth plate is injured, FoxA2+ cells expand in response to trauma, and produce physeal cartilage for growth plate tissue regeneration.

TLR4 activates NF-κB in human ovarian granulosa tumor cells.

Previous studies have demonstrated expression of Toll-like receptors (TLRs) in the surface epithelium of normal ovaries (OSE) and in epithelial ovarian tumors. Most notably, OSE-derived cancers express TLR4, which activates the nuclear factor-kappa B (NF-κB) signaling cascade as a mediator of inflammatory response. Currently, there is considerable interest in elucidating the role of TLR-mediated signaling in cancers. Nevertheless, the expression of TLRs in granulosa cell tumors (GCTs) of the ovary, and the extent to which GCT expression of TLRs may influence cell-signaling pathways and/or modulate the efficacy of chemotherapeutics, has yet to be determined. In the present study, human GCT lines (COV434 and KGN) were utilized to evaluate expression of functional TLR4. TLR4 is expressed in GCT cell lines and ligation of TLR4 with bacterial lipopolysaccharide (LPS) led to IκB degradation and activation of NF-κB. NF-κB activation was confirmed by nuclear localization of NF-κB p65 following treatment with LPS and the naturally occurring ligand, HSP60. Notably, immunoneutralization of TLR4 blocked nuclear localization, and inhibition of NF-κB signaling attenuated LPS-induced TNFα plus increased doubling time in both cell lines. Contradictory to reports using human OSE cell lines, inhibition of NF-κB signaling failed to sensitize GCT lines to TRAIL or cisplatin. In summary, findings herein are the first to demonstrate a functional TLR-signaling pathway specifically in GCTs, and indicate that in contrast to OSE-derived cancers, inhibition of NF-κB does not sensitize GCTs to TRAIL or cisplatin.

Role of Granulosa Cells in the Aging Ovarian Landscape: A Focus on Mitochondrial and Metabolic Function.

Mitochondria are at the intersection of aging and fertility, with research efforts centered largely on the role that these specialized organelles play in the relatively rapid decline in oocyte quality that occurs as females approach reproductive senescence. In addition to various roles in oocyte maturation, fertilization, and embryogenesis, mitochondria are critical to granulosa cell function. Herein, we provide a review of the literature pertaining to the role of mitochondria in granulosa cell function, with emphasis on how mitochondrial aging in granulosa cells may impact reproduction in female mammals.

Biomechanical Strain Promotes the Differentiation of Murine Oogonial Stem Cells.

Cells within tissues are routinely subjected to physiological stress and strain, arising from direct interactions with neighboring cells as well as with extracellular matrix components. Accordingly, there is tremendous interest in deciphering how cells sense, and respond to, changes in biomechanical forces. In this study, we explored the effects of mechanostimulation on the differentiation of mouse female germline or oogonial stem cells (OSCs) as a model for adult stem cell function. We report that increasing levels, or repeated application of a subthreshold fixed level, of radial strain to OSCs in culture significantly increased rates of in vitro oocyte formation as a measure of stem cell differentiation. These responses involved changes in F-actin-mediated cytoskeletal tension as well as in activation of intracellular signaling by Rho-associated protein kinase (ROCK) and Yes-associated protein (YAP) phosphorylation. In addition, application of mechanical strain to OSCs enhanced association of YAP with muscle-specific cytidine-adenosine-thymidine (MCAT) response elements in the promoter stimulated by retinoic acid gene 8 (Stra8), the transcriptional activation of which is required for germline meiotic commitment. These data indicate that biomechanical strain directly promotes the differentiation of adult female germline stem cells through a signaling pathway involving F-actin, ROCK, YAP, and Stra8.