Katsiaryna Tarbashevich1
Erez Raz- 1Institute of Cell Biology, Center for Molecular Biology of Inflammation (ZMBE), Muenster, Germany
- 2Max Planck Institute for Molecular Biomedicine, Münster, Germany
To maintain a range of cellular functions and to ensure cell survival, cells must control their levels of reactive oxygen species (ROS). The main source of these molecules is the mitochondrial respiration machinery, and the first line of defense against these toxic substances is the mitochondrial enzyme superoxide dismutase 2 (Sod2). Thus, investigating early expression patterns and functions of this protein is critical for understanding how an organism develops ways to protect itself against ROS and enhance tissue fitness. Here, we report on expression pattern and function of zebrafish Sod2, focusing on the role of the protein in migration and maintenance of primordial germ cells during early embryonic development. We provide evidence that Sod2 is involved in purifying selection of vertebrate germ cells, which can contribute to the fitness of the organism in the following generations.
Introduction
The control over the level of reactive oxygen species (ROS) in the cell is important for a range of physiological and pathological processes (reviewed in (Bowling et al., 2019; Lawlor et al., 2020; Sie et al., 2020; Shields et al., 2021)). A major source of these harmful ROS compounds is oxidative phosphorylation in the mitochondria, but mitochondria also have a way to handle them: ROS are neutralized by detoxifying enzymes such as superoxide dismutase 2 (SOD2), a mitochondrial enzyme that converts superoxide radicals into the less harmful compound hydrogen peroxide. Consequently, the function of SOD2 is critical for cell survival in Drosophila (Celotto et al., 2012) and for the immune response in zebrafish (Peterman et al., 2015). Mice that have a reduced SOD2 function are more likely to develop epilepsy (Liang and Patel, 2004), cardiomyopathy defects (Van Remmen et al., 2001), and oncological and aging-related phenotypes ((Justilien et al., 2007; Velarde et al., 2012; Weyemi et al., 2012); reviewed in (Lee et al., 2009)).
While SOD function has been shown to be important in different cell types, the role of the protein during the early development of vertebrate embryos is not well understood, in particular regarding its role in selecting for fitter germline cells that ensure the robustness of the organism across generations. So far, the process of “purifying selection” based on cell competition, namely, the selection for cells with optimal energy supply or metabolic status, has been described for somatic cells during early mouse embryogenesis ((Döhla et al., 2022; Lima et al., 2021)). In the mouse, natural selection of male germ cells was reported to be based on differentiation efficiency but the impact of metabolic parameters on the process was not analyzed ((Nguyen and Laird, 2021; Nguyen et al., 2020)). Thus, not all aspects of competition among germ cells have been studied. To explore these issues, we employed the zebrafish model, benefiting from the extrauterine development of the embryo and the available genetic and reverse genetic tools. We examined the role of Sod2 in early embryogenesis with respect to cell maintenance as well as collective and single cell migration, focusing on its function in germ cells. We show that Sod2 is expressed in germ cells and is important for their development and provide evidence that Sod2 protein is involved in the selection for fitter germ cells, thereby contributing to the fitness of the organism across generations.
Materials and methods
Fish lines and husbandry
The following zebrafish (Danio rerio) lines were employed: wild-type of the AB background, Tg(kop:EGFP-F’nos3′UTR) (Blaser et al., 2005), Tg(kop:mcherry-F-nos3′UTR) (Tarbashevich et al., 2015), Tg(gsc:GFP) (Doitsidou et al., 2002), Tg (cldnB:lynEGFP)+/+ (Haas and Gilmour, 2006), mu13 (MZsod2STOP31) Tg(kop:EGFP-F’nos3′UTR), mu14 (MZsod2STOP19) Tg (cldnB:lynEGFP)+/+.
WT and manipulated embryos were collected, kept in 0.3× Danieau’s solution [17.4 mM NaCl, 0.21 mM KCl, 0.12 mM MgSO4·7H2O, 0.18 mM Ca(NO3)2, and 1.5 mM Hepes (pH 7.6)], and raised at 28°C. The general zebrafish maintenance was performed in compliance with the German, North-Rhine-Westphalia state law, following the regulations of the Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen and was supervised by the veterinarian office of the city of Muenster.
Microinjections into the zebrafish embryos
For all experiments except for the one presented in Figure 7, embryos were microinjected into the yolk with 2 nL of the injection mixture. The following capped sense mRNAs synthesized using the mMESSAGE mMACHINE (Thermo Fisher Scientific): GrpEL-EGFP (C131.GrpEL.EGFP.nos) (used as mitochondria marker for co-localization studies), Sod2-mCherry (C247.Sod2-mCherry-nos) (used to determine intracellular protein localization), Sod2 (C382.Sod2.nos, E046.Sod2ORF.sod2 3′UTR) (both used for rescue experiments), mCherry-H2B (B325. mCherry.H2B.globinUTR) (labels nuclei of all cells in the embryo, used for the correction for the gastrulation movement in the cell migration analysis), Mito-mCherry (C433.mCherryHA-Bcl2TM.nos3′UTR) (used for the fluorescence-based analysis of the mitochondria numbers in the PGCs).The mRNAs containing nos 3′UTR are cleared from somatic cells via miRNA430-mediated degradation and are, thus, predominantly expressed in germ cells (Kedde et al., 2007).
Translation-blocking morpholinos (MO, Gene tools) utilized for this work were: CntrMO (CCTCTTACCTCAGTTACAATTTATA), Sod2MO (CATGCTCTAGTCCGTCACACAGTGA). Both morpholinos were injected in 0.5 µM concentration.
For the experiment described in Figure 7, 8-cell-stage embryos were microinjected into one of the middle blastomeres with 1 nL injection mix containing Cas9 nuclease (0.5 μg/μL, Integrated DNA technologies (IDT), cat. Nr. 1081059), guide crRNA egfp (18 ng/μL, 5′-AAGGGCGAGGAGCTGTTCAC-3′), guide crRNA sod2 (18 ng/μL, 5′-ACAAATCTGTCACCCAATGGCGG-3′), tracrRNA (33.5 ng/μL, cat. Nr. 1072534) according to the manufacturer’s recommendations (IDT).
Transplantation experiments
For the posterior lateral line primordium (pLLP) transplantation experiments, wild-type and MZsod2STOP31 embryos were injected with 500 pg of dextran conjugates (Invitrogen™, Dextran, Alexa Fluor™ 568; 10,000 MW, Anionic, Fixable, Catalogue number: D22912 or Alexa Fluor™ 680; 10,000 MW, Anionic, Fixable, Catalogue number: D34680) of distinct colors and used as donors. Uninjected wild-type embryos carrying a cldnB:lynEGFP transgene served as hosts.
Equal numbers of both types of donor embryos were grown to 4 h post fertilization (hpf), before removing the yolks. The cell caps were combined in a tube and carefully dissociated by tapping, until a homogeneous cell mixture was obtained. This cell mixture was transplanted into shield stage host embryos. At 31 hpf, dextran-labeled donor cells were distributed throughout the host embryos. Transplanted embryos were screened for cases in which both wild-type and MZsod2STOP19 donor cells had integrated into the host pLLP. The primordium was then imaged with Z-steps of 20 µm.
For the PGC transplantations, wild-type (WT) (Tg(kop:EGFP-F’nos3′UTR)+/+) or sod2 knockout (KO) (MZsod2STOP31 Tg(kop:EGFP-F’nos3′UTR)+/+) germ cells from 4.7-5 hpf donor embryos expressing EGFP-F′ transgene and dominant marker Tg(cry:DsRed) (red eyes) were transplanted into 4 hpf hosts (Tg(kop:mcherry-F-nos3′UTR)) expressing mCherry-F′ transgene and dominant marker Tg(cmlc:EGFP) (green heart) or into “dark” MZsod2STOP31 embryos. Embryos with mosaic germ cell clusters were selected at 24 hpf and raised to adulthood.
Immunohistochemistry
Whole-mount in situ hybridization using the DIG-labeled probe for the sod2, egr2a, mab21l, en2a and foxa3d was performed as previously described (Weidinger et al., 2002).
RNAscope was performed as described previously (Gross-Thebing et al., 2014) using the following probes: sod2 (cat. Nr. 300031-C2/320269-C2), vasa (cat. Nr. 407271-C3), cxcl12a (cat. Nr. 406481), noto (cat. Nr. 4835511-C2), gata1(cat. Nr. 473371-C3), neurog1 (cat. Nr. 505081-C2), myoD (cat. Nr. 402461-C2), pax8 (cat. Nr. 494721-C2), pcdh8 (cat. Nr. 494741-C3).
Fluorescence whole-mount immunostaining of the telomere repeats was performed using TelC-Cy3 (cat. Nr. TP-007) and TelC-Cy5 (cat. Nr. TP-009) according to the manufacturer’s instructions (PNA Biotech).
Image acquisition and microscopy
For live imaging, embryos were dechorionated, and those older than 20 hpf were anesthetized using tricaine (concentration, A5040, Sigma-Aldrich) in 0.3× Danieau’s solution, mounted in agarose-coated ramps covered with Danieau’s solution and manually oriented. For experiments aimed at determining migration parameters, samples were incubated at 28°C on a heated stage (PeCon, TempController 2000-2).
Spinning disk confocal microscopy was performed using Carl Zeiss Axio Imager Z1 and M1 microscopes equipped with a Yokogawa CSUX1FW-06P-01 spinning disk units. Imaging was performed using a Hamamatsu ORCA-Flash4.0 LT C11440 camera and Visitron Systems acquisition software (VisiView). Imaging of RNAscope samples (×10 objective) was conducted by acquiring 400-μm z-stacks (80 Z planes, 5 μm apart). For the analysis of PGC speed, ×10 objective multi-stage time-lapses were acquired for 2 h with 2-min acquisition intervals. For the PGC blebbing analysis, ×63 objective high-magnification movies of individual cells were acquired with 5-s intervals for 5 min. For the neomycin treatment, dechorionated embryos were incubated in a 400 µM neomycin (Sigma, cat. Nr. N1142) solution in Danieau’s buffer from 3 hpf until the end of the experiment. For the control treatments, the same amount (in µl) of DMSO (Sigma, cat. Nr. D8418) was used.
For the analysis of the pLLP migration, time-lapse movies were acquired at a rate of one frame every 30 min between 24 hpf 31 hpf. The samples were imaged using a1×0 objective to obtain150-µm Z-stacks of 10 µm optical slices.
High-magnification images for determining cell positioning and volume of the pLLP (40x) were performed by acquiring 2-μm optical slices. Analysis of the pLLP volumes was performed using the surface function of the Imaris software (Bitplane, version 9.5.1.).
Confocal laser scanning image acquisition was performed using LSM710 (Zeiss) upright microscope and ZEN software (Zeiss). If not stated otherwise, imaging was performed using a ×63 water-dipping objective with 2-μm optical slices.
Determination of the embryo size
Determination of zebrafish embryo size at the 5-somite stage was performed using the Fiji distribution function of ImageJ (NIH). First, the nuclei signal (Hoechst) of RNAscope image data, was projected along the Z-axis (400 µm in total), performing a maximum intensity projection. Rolling ball background-subtraction (radius = 50 pixels) was subsequently applied to the projected image. Next, a binary image was created by thresholding the image using the Otsu algorithm and adjusting the threshold such that the entire embryo was included. The “Fill holes” function was used to fill the area between the nuclei that were previously not segmented. The area of the obtained embryo segmentations was measured (for wild-type and MZsodSTOP31 embryos).
Analysis of cell migration
Analysis of the PGC migration speed was performed by tracking the germ cells using Imaris software (Bitplane, version 9.5.1.), with the correction of the gastrulation movement of the surrounding somatic cells. Analysis of PGC blebbing activity was performed manually (bleb count) using Fiji software. To determine the portion of blebs generated by PGCs at the back of a cell relative to the total number of blebs, the cells were divided into front and back halves and the total number of blebs, as well as the number of blebs at the back was counted for 5-min time-lapse movies.
For the analysis of the PGC cluster length, the value obtained for a cluster on one side of an embryo was normalized to the length of the yolk extension of the given embryo. Analysis was performed using Fiji software.
Analysis of the pLLP migration speed and distance were performed using Fiji software. The image orientation was adjusted, such that the lateral line was perfectly horizontal, with the primordium traveling to the right. The anterior point of the yolk extension was used as the point from which the distance to the tip of the traveling primordium was measured. The obtained values were normalized to the average wild-type pLLP speed/distance of travel for the given experiment.
qPCR
qPCR was performed as described in Cawthon (2009). Telomere primers: 5′ CGGTTTGTTTGGGTTTGGGTTTGGGTTTGGGTTTGGGTT-3′ and 5′- GGCTTGCCTTACCCTTACCCTTACCCTTACCCTTACCCT-3’. The house keeping gene hmg (high mobility group) was used as a reference RNA amplified with primers 5′-GTGGAAGACCCCGAAAACA-3′ and 5′- TTCCTCCTCTTCTTCCTCCTCG-3’. Additional normalization was performed to the values of the WT group.
Western blotting analysis
Western blotting analysis was performed on protein extracts from wild-type and MZsod2STOP embryos of 24hpf. The following primary antibodies were used: anti-Sod2 antibody (Proteintech, #24127-1-AP, in 1:3000 dilution) and a pan-Actin antibody (Thermofisher, #MA1-744, in 1:3000 dilution). For detection the following secondary antibodies were employed: goat-anti-mouse (Li-Cor, 926-68071, IRDye 680RD) and donkey-anti-rabbit (Li-Cor, 926-32213, IRDye 800RD) in 1:5000 dilution.
ROS detection
For the detection of the level of the oxidative stress, 1hpf wild-type and MZsod2STOP embryos were manually dechorionated and incubated for 3 h in 10 µM CellRox (ThermoFisher, C10422) in Danieau’s at 28°C followed by removal of the chemical and confocal imaging.
Statistical analysis
Statistical analysis (one-way ANOVA or Student’s t-test) was performed using GraphPad Prism software (version 8). If not stated otherwise, all experiments were performed in three independent replicates. Where applicable, multiple comparisons were statistically analyzed using Dunn’s multiple comparison test.
Results
sod2 mRNA is enriched in migratory cell populations
Upon screening for RNA molecules enriched in zebrafish primordial germ cells (deep sequencing data published in (Hartwig et al., 2014; Paksa et al., 2016)), we identified superoxide dismutase type 2 (sod2) (ZFIN:ZDB-GENE-030131-7742). During early stages of embryogenesis (4 and 7 h post fertilization (hpf)), we found that sod2 mRNA was expressed globally, but enriched in primordial germ cells (PGCs) (Figures 1A,B, Supplementary Figure S1A). At later stages of embryogenesis, sod2 expression was detected in several structures such as the caudal somites as well as in specific regions within the developing brain (Figure 1C, Supplementary Figure S1A and Supplementary Figure S2). Intriguingly, in addition to being expressed in the PGCs (blue arrowhead in Figure 1C, Supplementary Figure S1, Supplementary Figure S2), sod2 RNA was strongly expressed within two other populations of migratory cells: caudal hematopoietic trunk cells (CHT) (yellow arrowhead in Figure 1C, Supplementary Figure S1A) and cells of the posterior lateral line primordium (pLLP) (purple arrowhead in Figure 1C, Supplementary Figure S1A and Supplementary Figure S2). Consistent with its function, the Sod2-mCherry fusion protein localized to the mitochondria of cells of early embryos (e.g., in germ cells expressing the mitochondrial protein GrpEL (Figure 1D-D’’)) (Srivastava et al., 2017).
FIGURE 1. sod2 RNA expression pattern and localization of the protein during early embryogenesis. (A, B) RNAscope-based localization of sod2 mRNA and the germ cell marker vasa mRNA at (A) 4 hpf and (B) 7 hpf. (C) Whole-mount in situ hybridization image probing for sod2 mRNA in 24 hpf embryos. Blue arrowhead points at the PGC cluster, yellow arrowhead at hematopoietic precursor cells, and purple arrowhead at the migrating pLLP. (D–D’’) Representative confocal images of Sod2-mCherry protein in germ cells relative to the mitochondrial marker (GrpEL-EGFP). Scale bar 10 µm.
Generation of sod2 maternal zygotic mutant fish lines
To study the function of Sod2 during early zebrafish embryogenesis, we generated two knockout (KO) alleles employing CRISPR/Cas9 technology. Both lines carry mutations at the beginning of the catalytic domain of the wild-type (WT) protein, resulting in frameshifts and stop codons shortly following the mutation sites (Figure 2A). Thus, we considered the two mutated alleles to result in full loss of function of Sod2 (see below).
FIGURE 2. Generation of sod2 maternal zygotic fish lines. (A) The schematic representation of the position of the two Cas9-induced mutations. The mutations resulted in stop codons at amino acid position 31 (sod2STOP31) and 19 (sod2STOP19) (marked by asterisks). (B) Western blotting analysis of the Sod2 protein expression in wild type and MZsod2 embryos.
Since sod2 transcripts (and presumably the protein) are maternally provided to the embryo, we aimed at generating adult mutant fish that were homozygous for both knockout alleles. In this way, one can obtain embryos that do not contain maternally-provided material that are also incapable of expressing functional RNA from their own genome.
Interestingly, despite the central role that the protein plays in controlling ROS levels in the cell, under lab conditions the fish lacking the function of Sod2 were viable and developed into adults that proved to be fertile. This allowed us to generate and analyze embryos lacking both maternal and zygotically transcribed sod2 RNAs (maternal-zygotic mutants (MZ), we designated as MZsod2STOP19 and MZsod2STOP31). The schemes of crosses and genotypes analyzed in this manuscript are presented in Figure 3A. In such embryos, presumably as a result of nonsense-mediated RNA decay, no sod2 RNA (whole-mount in situ hybridization, Supplementary Figure S1B) or protein (Western blotting analysis, Figure 2B) could be detected.
FIGURE 3. MZsod2 KO embryos exhibit reduction in the germ cell migration speed. (A) The scheme of crosses and genotypes analyzed in this manuscript. (B, C) PGC migration speed was determined by tracking cell migration for 70 min using Imaris software. (D) Normalized PGC cluster lengths at 24 hpf in WT and MZsod2 KO embryos treated with neomycin (enhanced ROS production condition). (E, F) Representative images of embryos analyzed for the PGC cluster length (quantification shown in (D)) ****p < 0.0001, *p < 0.05 was determined by Student’s t-test. n - number of PGCs analyzed in 38 WT and 36 MZsod2 KO embryos (B), 24 MZsod2 KO + control (contr) mRNA and 29 MZsod2 KO + sod2 mRNA embryos (C). N - number of PGC clusters (one cluster per embryo) analyzed.
Nevertheless, a phenotype we did observe in MZ mutant embryos at early somitogenesis stages (11 hpf) was that of a mild reduction in body size (Supplementary Figure S3A). However, as judged by RNA expression, patterning events progressed normally in all germ layers, i.e., ectoderm, mesoderm and endoderm (Supplementary Figure S3B, C), and the embryos appeared morphologically normal at 5 days of development. Relevant for this report, the expression of cxcl12a, namely, the chemokine guiding the migration of both pLLP and PGCs, was expressed normally in the mutant embryos (Supplementary Figure S3C).
Sod2 function is dispensable for the collective migration of the pLLP cell cluster
The posterior lateral line primordium (pLLP) is a cluster of 100–150 cells that migrate collectively towards the tail of the zebrafish embryo (Haas and Gilmour, 2006). Interestingly, sod2 mRNA was expressed within cells of the pLLP (Figure 1C, Supplementary Figure S1A), with elevated levels detected at the front of the cell cluster compared to the rear (Supplementary Figure S2). This expression pattern prompted us to examine the migration of the cell cluster in embryos lacking Sod2 function (MZsod2STOP19). As judged by the position of the pLLP at 24 hpf to 35 hpf (Supplementary Figure S4A), the cell clusters in mutant embryos migrated at the same speed as the non-manipulated wild-type clusters. Similarly, we did not observe an effect on the cluster size, suggesting a normal rate of cell division and death (S4B-B’’). Last, in mosaic pLLP clusters in which only some of the cells lacked Sod2 function and the rest were wild-type, we found no role of the protein in dictating the position of the cells within the primordium (Supplementary Figure S4C–E).
Thus, based on the parameters investigated in Supplementary Figure S4, Sod2 function was not essential for proper migration of the pLLP, nor for cell behavior within the cluster.
Elevation of ROS affects PGC migration
To investigate the possible function of Sod2 in another cell population that expresses the protein, we turned to primordial germ cells (PGCs). In different organisms PGCs are specified during early embryonic development and migrate as single cells from their site of specification to the region where the gonad develops to give rise to gametes, sperm and egg ((Richardson and Lehmann, 2010; Aalto et al., 2021; Grimaldi and Raz, 2020)). We examined the migration speed of PGCs in seven to nine hpf embryos and observed a mild reduction in the speed of sod2 mutant PGCs as compared with their wild-type counterparts (Figure 3B). In line with this finding, cells devoid of Sod2 function were less polar as manifested by an increase in the rate of bleb-type protrusion formation, particularly at the cell rear (Figures 4A–C, Supplementary Movie S1). Conversely, overexpression of Sod2 led to an increase in migration speed, suggesting that the Sod2 activity level constitutes a limiting factor for germ cell translocation (Figure 3C). Despite the defects in protrusion formation and the slower migration of PGCs in MZsod2STOP31 embryos, normal size of PGC cluster formed at the region where the gonad develops (Supplementary Figure S5A).
FIGURE 4. Defects in protrusion formation and germ cell polarity in sod2 KO and neomycin-treated embryos. Analysis of (A, D) the number of blebs in PGCs per minute and (B, E) the proportion of blebs at the cell back. In graphs A and B, MZsod2STOP31 embryos were injected with control RNA or with RNA directing Sod2 expression to the germ cells (see Methods); the same cells were analyzed in panels A and (B). In graphs D and E, wild-type embryos were treated with DMSO or with neomycin; the same cells were analyzed in panels D and (E). (C, F) Snapshots from the representative time-lapse movies (movies S1 and S2) analyzed in (A,B,D, E). Yellow asterisks point at blebs at the rear of cells. n - the number of PGCs analyzed. One PGC per embryo was imaged and analyzed. **p < 0.01, *p < 0.05, were calculated by Student’s t-test.
In addition to the finding presented above, we detected higher ROS levels in cells deficient for Sod2 function (Supplementary Figure S5B). To examine whether Sod2 confers robustness to the migration process under conditions of increased oxidative stress in the PGCs, we subjected MZsod2STOP31 and wild-type embryos to neomycin (Neo), which was shown to increase oxidative stress in cells ((Harris et al., 2003; Esterberg et al., 2016)). Increasing the oxidative stress in this way led to an increase in bleb formation and a reduction in cell polarity (Figures 4D–F, Movie S2) in wild-type embryos. We next checked whether the function of Sod2 could confer robustness to the arrival of germ cells at their target under these conditions. Indeed, increasing ROS level by subjecting embryos to Neomycin further enhances the oxidative stress in cells lacking SOD2 function (Supplementary Figure S5C). Interestingly, loss of Sod2 function coupled to neomycin treatment affected PGC cluster formation at the region where the gonad develops (Figures 3D–F).
Sod2 function is important for selection of fitter germ cells
To further examine the function of Sod2 in the germline, we monitored the mitochondria, the organelle within which the enzyme functions. Measuring the fluorescence intensity levels of a mitochondria-directed mCherry protein in 8 and 24 hpf embryos revealed a very strong reduction in the average fluorescence levels particularly in MZsod2STOP31 PGCs at 24 hpf as compared with their wild-type counterparts (Figure 5, Supplementary Figure S6).
FIGURE 5. MZsod2 KO germ cells exhibit reduced mitochondria-derived fluorescence signal. Average intensity of the fluorescently-labeled mitochondria in WT and MZsod2 KO germ cells at 8 and 24 hpf. n - number of PGCs analyzed in N–number of embryos. ****p < 0.0001, *p < 0.05 was determined by ANOVA test. Scale bar 10 µm.
It has previously been reported that, as quantified by the amount of mitochondrial DNA, the number of mitochondria decreases exponentially from the time of egg fertilization until late larval stages, when replication of these organelles reinitiates (Otten et al., 2016). In light of this, we hypothesized that the further decrease in mitochondria-derived fluorescence signal in MZsod2STOP germ cells could result in reduced fitness of these cells relative to wild-type cells during later stages of development (Figure 9).
To test this hypothesis, we transplanted germ cells from wild-type and MZsod2STOP31 lines into either wild-type or MZsod2STOP31 hosts, with the transplanted donor germ cells carrying a dominant marker (fluorescent eyes, Supplementary Figure S7). Chimeric embryos containing both host and donor germ cells were raised and mated with wild-type fish, and we scored for the embryos expressing the donor dominant marker (fluorescent eyes), which signifies the fitness of transplanted PGCs (Figure 6).
FIGURE 6. Analysis of zebrafish with mosaic gonads: PGC transplantation approach. Embryos with mosaic germ cell clusters (white cells are the transplanted donor PGCs, red cells are the host PGCs in the upper panels) were selected at 24 hpf and raised to adulthood and crossed to wild-type fish. If clutches contained embryos with the dominant marker of the donor line (red eyes), the outcrossed adult fish would be counted as “positive.” Lower graph shows the number of such positive fish relative to the total number of adults tested is presented in percentage for each background. Numbers within the bars represent the number of F0 adults with donor progeny divided by the total number of F0 adults raised for this line.
Following this experimental scheme, we detected F1 progeny derived from the transplanted germ cells in 7 of 13 (54%) adult fish when wild-type donor PGCs were transplanted into wild-type hosts. Similarly, 8 out of 16 (50%) adult fish had progeny expressing the donor domain marker when MZsod2STOP31 PGCs were transplanted into MZsod2STOP31 hosts. The result of the experiment changed when germ cells were transplanted into heterologous backgrounds. In these cases, F1 progeny expressing donor dominant makers were observed in 11 of 16 (69%) animals in which wild-type germ cells transplanted into the MZsod2STOP31 background. Conversely, the number of adult wild-type hosts giving rise to dominant-marker carrying progeny was only 4 of 11 fish (36%) in the case of MZsod2STOP31 cells transplanted into wild-type hosts.
Due to the complexity and low throughput of the experiment presented above, only a relatively small number of chimeric animals could be generated. Thus, as an additional test for the hypothesis, we generated fish with chimeric gonads containing wild-type and MZsod2STOP19 germ cells using a different approach (Figure 7A). To this end, we co-injected Cas9 protein and guide RNAs targeting egfp and sod2 sequences into 1 blastomere of 8-cell-stage wild-type and MZsod2STOP19 embryos. The injected embryos were homozygous for an EGFP transgene. The injected embryos were raised, and adult males were crossed with wild-type female fish, with the progeny analyzed for the percentage of embryos lacking a GFP signal. Thus, knocking out sod2 according to this scheme resulted in adult fish with a chimeric germline, where progeny lacking the EGFP signal representing the proportion of sod2 KO germ cells that contributed to the germline. Intriguingly, we found that the potency of cells lacking Sod2 to contribute to the germline and form gametes was very strongly reduced in wild-type animals as compared to their ability to do so when located within a gonad containing only Sod2 mutant cells (Figure 7B). We thus conclude that while Sod2 mutant cells can generate gametes, they are much less likely to do so when positioned within an environment containing wild-type cells. Thus, when located within an environment of wild-type cells, the Sod2-deficient germline cells are outcompeted and contribute less to the gametes (WT in Figure 7B). Importantly, the effect of germ cell competition was not observed in similarly generated control chimeric gonads of embryos injected only with guide RNA for the egfp sequence (no sod2 knockout) (Figure 7C). Together, Sod2 deficient cells contribute less to the germline when located among wild-type germ cells. In contrast, the negative controls where Sod2 activity was not altered (anti gfp guides, or anti sod2 guides in sod2 mutant cells) the treated cells showed no difference in their contribution to the germline as compared with non-treated cells. Consistent with the importance of reduced ROS levels in germline cells for proper gametogenesis, MZsod2 homozygous mutant fish are less fertile (Supplementary Figure S8).
FIGURE 7. Generation and analysis of zebrafish with mosaic gonads: Cas9-based approach. (A) Embryos of WT or sod2 KO backgrounds homozygous for Tg(cldnB:lynEGFP) (Tg(GFP) +/+) were injected at the 8-cell stage into one of the middle blastomeres with Cas9 protein and two guide RNAs targeting egfp and sod2 sequences (sgEGFP and sgSod2). This manipulation affected only a fraction of germ cells per embryo, allowing for a high-throughput generation of animals with mosaic gonads containing WT germ cells (expressing Sod2 and EGFP transgene) and KO germs cells (where both Sod2 and EGFP functions were eliminated). Such embryos were raised to adulthood. Adult males were outcrossed to WT females, and the percentage of GFP-negative embryos (lacking EGFP and Sod2 expression) within the clutch was calculated. Results are presented in graph (B). The results of the negative control experiment with animals injected only with sgGFP (manipulating GFP expression only) are presented in graph (C). n–number of embryos and N - number of adult fish analyzed. ****p < 0.0001, n.s.p > 0.05 was determined by Student’s t-test.
An interesting phenotype we observed, which may be relevant for the reduced potential of Sod2-depleted germ cells to generate gametes as compared to non-manipulated neighboring germ cells, is a reduction in the number of telomeric repeats in MZsod2 KO embryos. This analysis was performed using qPCR (Figure 8A), or fluorescence in situ hybridization using the telomere-specific probe (TelC-Cy3, Figures 8B–D, Supplementary Figure S9). A possible basis for the shortening of telomeres is DNA damage that telomere length reports on. The reduction in telomere length we observe in sod2 mutants could affect the fitness of cells that carry the genetic information to future generations.
FIGURE 8. The number of telomeric repeats is reduced in sod2 KO fish lines. (A) The relative length of telomeric repeats in wild type and MZsod2 mutants was determined by qPCR (graph). The qPCR values were normalized to the WT condition. (B) The relative number of telomeric repeats was determined by fluorescence in situ hybridization (FISH) using Cy3-conjugated TelC probes. The graph presents intensity measurements of the TelC signal in individual PGC nuclei. n–number of PGCs analyzed in 32 WT and 35 MZsod2 embryos (C, D) Representative confocal images of TelC-FISH staining. PGCs are identified by the expression of EGFP-F′ on their cell membranes (represented by yellow dotted lines). ****p < 0.0001was determined by Student’s t-test. Scale bar 10 µm.
Discussion
The Sod2 protein has been shown to be involved in a range of physiological and pathological processes. For example, the Drosophila sod2 mutant (SOD2Bewildered) exhibits aberrant brain morphology, abnormal axonal targeting, and neurodegeneration (Celotto et al., 2012). We observed strong expression of sod2 mRNA in the zebrafish brain as well, but, based on the molecular markers we examined and the apparent normal behavior of MZsod2 adult fish, the loss of Sod2 in zebrafish did not result in a pronounced phenotype in this tissue. Since the ROS level has been shown to affect lipid metabolism ((Quijano et al., 2016; Zhuang et al., 2021)), the main energy source of the early embryo, the reduction in body size of 12-hour-old MZsod2STOP mutant embryos could stem from defects in the availability of these molecules (reviewed in (Sant and Timme-Laragy, 2018; Quinlivan and Farber, 2017)).
The PGC migration phenotype we observed in MZsod2 mutants and the similar phenotype of neomycin-treated embryos could reflect defects in energy production as well, which in this case could stem from defects in the function of the mitochondria (Esterberg et al., 2016). We attribute the lack of migration defects in the pLLP to the more robust nature of collective cell migration in which mild locomotion aberrations in some cells can be compensated by other cells in the cluster ((Knutsdottir et al., 2017; Yamaguchi et al., 2022)).
In contrast to the very mild or lack of phenotypes we observed in PGCs during early stages of embryonic development (first 1–2 days post fertilization), the phenotypes we detected at later stages were much more pronounced. A possible explanation could be that the PGC competition phenotype we observe reflects a cumulative effect of lack of Sod2 function over the 3-month period between the establishment of the germline and the generation of gametes.
A frequently documented type of cellular heterogeneity that is subjected to selection has been described in somatic cells in the form of mitochondrial heteroplasmy (e.g., mitochondrial function in human cells (Wei et al., 2019)). In this case, cells with an elevated proportion of defective mitochondria are eliminated by more robust neighboring cells in a process termed “purifying selection” (Lima et al., 2021). Heterogeneities in cellular functions have also been identified in other cell populations including stem cells ((Ellis et al., 2019; Lawlor et al., 2020; Lima et al., 2021), reviewed in (Muller-Sieburg et al., 2012; Krieger and Simons, 2015)). In germ cells, such heterogeneities can be especially consequential, as they are transmitted to the next-generation (Zhang et al., 2022).
Since during the embryonic stages studied here there was no de novo production of mitochondria, the number of the maternally provided organelles per cell declined with each cell division, reaching a very low number (estimated to be 100 per cell (Otten et al., 2016)). Indeed, the time of PGC migration we examined here corresponds to developmental stages with very low numbers of mitochondria per cell (Otten et al., 2016). The small number of mitochondria is an obvious parameter relevant for mitochondrial heteroplasmy in PGCs and for the so called mitochondrial “bottleneck” effect (Otten et al., 2016). In the case of germ cells, this time point can serve as the stage when purifying selection takes place, contributing to the transmission of healthier mitochondria to the progeny. Consistent with the Sod2 level being important for mitochondria function and DNA integrity, in MZsod2STOP embryos the already low mitochondria number was further reduced (schematically presented in Figure 9), and telomeres were shorter. The low level of mitochondria observed in PGCs in MZsod2STOP embryos and the aberrations in migration motivated us to study the possible evolutionary role of Sod2 in PGC maintenance and natural selection. Indeed, the results of the germline transplantation experiments (Figure 6) and those of the chimeric germline analysis (Figure 7) show that germ cells lacking Sod2 were less likely to generate gametes as compared with wild-type PGCs when positioned together within the developing gonad. Since such phenomenon has not been observed among cells of similar genotypes, these results fit the definition of cell competition ((de la Cova et al., 2004; Moreno and Basler, 2004)), which in this case involves germline cells. In the context of this work, such competition could select for germline cells with lower ROS-generated damage, thereby controlling the quality of gametes and, thereby, the fitness of somatic cells in the next-generation. Thus, it is likely that the full manifestation of the MZsod2STOP phenotype may only be appreciated on the evolutionary scale when followed over several generations.
FIGURE 9. Schematic illustration of changes in mitochondria numbers during zebrafish development. Representation of mitochondria numbers per germ cell at 8 hpf and 24 hpf of zebrafish development in WT (black) and MZsod2 KO (red) animals. The dotted part of the graph is based on (33).
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by Landesamt für Natur, Umwelt und Verbraucherschutz Nordrhein-Westfalen. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
KT, LE, JW and JP performed the experiments and analyzed the data. KT and ER designed the study and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the German Research Counsil (DFG), CRC1348 and CRU 326 program and by the Medical Faculty of the University of Muenster.
Acknowledgments
We thank E.-M. Messerschmidt, I. Sandbote and U. Jordan for technical assistance and C. Brennecka for comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2023.1250643/full#supplementary-material
References
Aalto, A., Olguin-Olguin, A., and Raz, E. (2021). Zebrafish primordial germ cell migration. Front. Cell Dev. Biol. 9, 684460. doi:10.3389/fcell.2021.684460
Blaser, H., Eisenbeiss, S., Neumann, M., Reichman-Fried, M., Thisse, B., Thisse, C., et al. (2005). Transition from non-motile behaviour to directed migration during early PGC development in zebrafish. J. Cell Sci. 118 (17), 4027–4038. doi:10.1242/jcs.02522
Bowling, S., Lawlor, K., and Rodríguez, T. A. (2019). Cell competition: the winners and losers of fitness selection. Development 146 (13), dev167486. doi:10.1242/dev.167486
Cawthon, R. M. (2009). Telomere length measurement by a novel monochrome multiplex quantitative PCR method. Nucleic Acids Res. 37 (3), e21. doi:10.1093/nar/gkn1027
Celotto, A. M., Liu, Z., Vandemark, A. P., and Palladino, M. J. (2012). A novel Drosophila SOD2 mutant demonstrates a role for mitochondrial ROS in neurodevelopment and disease. Brain Behav. 2 (4), 424–434. doi:10.1002/brb3.73
de la Cova, C., Abril, M., Bellosta, P., Gallant, P., and Johnston, L. A. (2004). Drosophila myc regulates organ size by inducing cell competition. Cell 117 (1), 107–116. doi:10.1016/s0092-8674(04)00214-4
Döhla, J., Kuuluvainen, E., Gebert, N., Amaral, A., Englund, J. I., Gopalakrishnan, S., et al. (2022). Metabolic determination of cell fate through selective inheritance of mitochondria. Nat. Cell Biol. 24 (2), 148–154. doi:10.1038/s41556-021-00837-0
Doitsidou, M., Reichman-Fried, M., Stebler, J., Köprunner, M., Dörries, J., Meyer, D., et al. (2002). Guidance of primordial germ cell migration by the chemokine SDF-1. Cell 111 (5), 647–659. doi:10.1016/s0092-8674(02)01135-2
Ellis, S. J., Gomez, N. C., Levorse, J., Mertz, A. F., Ge, Y., and Fuchs, E. (2019). Distinct modes of cell competition shape mammalian tissue morphogenesis. Nature 569 (7757), 497–502. doi:10.1038/s41586-019-1199-y
Esterberg, R., Linbo, T., Pickett, S. B., Wu, P., Ou, H. C., Rubel, E. W., et al. (2016). Mitochondrial calcium uptake underlies ROS generation during aminoglycoside-induced hair cell death. J. Clin. Invest. 126 (9), 3556–3566. doi:10.1172/JCI84939
Grimaldi, C., and Raz, E. (2020). Germ cell migration—evolutionary issues and current understanding. Semin. Cell Dev. Biol. 100, 152–159. doi:10.1016/j.semcdb.2019.11.015
Gross-Thebing, T., Paksa, A., and Raz, E. (2014). Simultaneous high-resolution detection of multiple transcripts combined with localization of proteins in whole-mount embryos. BMC Biol. 12 (1), 55. doi:10.1186/s12915-014-0055-7
Haas, P., and Gilmour, D. (2006). Chemokine signaling mediates self-organizing tissue migration in the zebrafish lateral line. Dev. Cell 10 (5), 673–680. doi:10.1016/j.devcel.2006.02.019
Harris, J. A., Cheng, A. G., Cunningham, L. L., MacDonald, G., Raible, D. W., and Rubel, E. W. (2003). Neomycin-induced hair cell death and rapid regeneration in the lateral line of zebrafish (Danio rerio). J. Assoc. Res. Otolaryngol. JARO 4 (2), 219–234. doi:10.1007/s10162-002-3022-x
Hartwig, J., Tarbashevich, K., Seggewiß, J., Stehling, M., Bandemer, J., Grimaldi, C., et al. (2014). Temporal control over the initiation of cell motility by a regulator of G-protein signaling. Proc. Natl. Acad. Sci. 111 (31), 11389–11394. doi:10.1073/pnas.1400043111
Justilien, V., Pang, J. J., Renganathan, K., Zhan, X., Crabb, J. W., Kim, S. R., et al. (2007). SOD2 knockdown mouse model of early AMD. Invest. Ophthalmol. Vis. Sci. 48 (10), 4407–4420. doi:10.1167/iovs.07-0432
Kedde, M., Strasser, M. J., Boldajipour, B., Vrielink, JAFO, Slanchev, K., Le Sage, C., et al. (2007). RNA-binding protein Dnd1 inhibits MicroRNA access to target mRNA. Cell 131 (7), 1273–1286. doi:10.1016/j.cell.2007.11.034
Knutsdottir, H., Zmurchok, C., Bhaskar, D., Palsson, E., Dalle Nogare, D., Chitnis, A. B., et al. (2017). Polarization and migration in the zebrafish posterior lateral line system. PLOS Comput. Biol. 13 (4), 1005451. doi:10.1371/journal.pcbi.1005451
Krieger, T., and Simons, B. D. (2015). Dynamic stem cell heterogeneity. Development 142 (8), 1396–1406. doi:10.1242/dev.101063
Lawlor, K., Pérez-Montero, S., Lima, A., and Rodríguez, T. A. (2020). Transcriptional versus metabolic control of cell fitness during cell competition. Semin. Cancer Biol. 63, 36–43. doi:10.1016/j.semcancer.2019.05.010
Lee, Y. H., Lin, Q., Boelsterli, U. A., and Chung, M. C. M. (2009). The Sod2 mutant mouse as a model for oxidative stress: a functional proteomics perspective. Mass Spectrom. Rev. 29, 179–196. doi:10.1002/mas.20226
Liang, L. P., and Patel, M. (2004). Mitochondrial oxidative stress and increased seizure susceptibility in Sod2−/+ mice. Free Radic. Biol. Med. 36 (5), 542–554. doi:10.1016/j.freeradbiomed.2003.11.029
Lima, A., Lubatti, G., Burgstaller, J., Hu, D., Green, A. P., Di Gregorio, A., et al. (2021). Cell competition acts as a purifying selection to eliminate cells with mitochondrial defects during early mouse development. Nat. Metab. 3 (8), 1091–1108. doi:10.1038/s42255-021-00422-7
Moreno, E., and Basler, K. (2004). dMyc transforms cells into super-competitors. Cell 117 (1), 117–129. doi:10.1016/s0092-8674(04)00262-4
Muller-Sieburg, C. E., Sieburg, H. B., Bernitz, J. M., and Cattarossi, G. (2012). Stem cell heterogeneity: implications for aging and regenerative medicine. Blood 119 (17), 3900–3907. doi:10.1182/blood-2011-12-376749
Nguyen, D. H., and Laird, D. J. (2021). Natural selection at the cellular level: insights from male germ cell differentiation. Cell Death Differ. 28 (7), 2296–2299. doi:10.1038/s41418-021-00812-0
Nguyen, D. H., Soygur, B., Peng, S. P., Malki, S., Hu, G., and Laird, D. J. (2020). Apoptosis in the fetal testis eliminates developmentally defective germ cell clones. Nat. Cell Biol. 22 (12), 1423–1435. doi:10.1038/s41556-020-00603-8
Otten, A. B. C., Theunissen, T. E. J., Derhaag, J. G., Lambrichs, E. H., Boesten, I. B. W., Winandy, M., et al. (2016). Differences in strength and timing of the mtDNA bottleneck between zebrafish germline and non-germline cells. Cell Rep. 16 (3), 622–630. doi:10.1016/j.celrep.2016.06.023
Paksa, A., Bandemer, J., Hoeckendorf, B., Razin, N., Tarbashevich, K., Minina, S., et al. (2016). Repulsive cues combined with physical barriers and cell–cell adhesion determine progenitor cell positioning during organogenesis. Nat. Commun. 7 (1), 11288. doi:10.1038/ncomms11288
Peterman, E. M., Sullivan, C., Goody, M. F., Rodriguez-Nunez, I., Yoder, J. A., and Kim, C. H. (2015). Neutralization of mitochondrial superoxide by superoxide dismutase 2 promotes bacterial clearance and regulates phagocyte numbers in zebrafish. Infect. Immun. 83 (1), 430–440. doi:10.1128/IAI.02245-14
Quijano, C., Trujillo, M., Castro, L., and Trostchansky, A. (2016). Interplay between oxidant species and energy metabolism. Redox Biol. 8, 28–42. doi:10.1016/j.redox.2015.11.010
Quinlivan, V. H., and Farber, S. A. (2017). Lipid uptake, metabolism, and transport in the larval zebrafish. Front. Endocrinol. 8, 319. doi:10.3389/fendo.2017.00319
Richardson, B. E., and Lehmann, R. (2010). Mechanisms guiding primordial germ cell migration: strategies from different organisms. Nat. Rev. Mol. Cell Biol. 11 (1), 37–49. doi:10.1038/nrm2815
Sant, K. E., and Timme-Laragy, A. R. (2018). Zebrafish as a model for toxicological perturbation of yolk and nutrition in the early embryo. Curr. Environ. Health Rep. 5 (1), 125–133. doi:10.1007/s40572-018-0183-2
Shields, H. J., Traa, A., and Van Raamsdonk, J. M. (2021). Beneficial and detrimental effects of reactive oxygen species on lifespan: a comprehensive review of comparative and experimental studies. Front. Cell Dev. Biol. 9, 628157. doi:10.3389/fcell.2021.628157
Sies, H., and Jones, D. P. (2020). Reactive oxygen species (ROS) as pleiotropic physiological signalling agents. Nat. Rev. Mol. Cell Biol. 21 (7), 363–383. doi:10.1038/s41580-020-0230-3
Srivastava, S., Savanur, M. A., Sinha, D., Birje, A., Vigneshwaran, R., Saha, P. P., et al. (2017). Regulation of mitochondrial protein import by the nucleotide exchange factors GrpEL1 and GrpEL2 in human cells. J. Biol. Chem. 292 (44), 18075–18090. doi:10.1074/jbc.M117.788463
Tarbashevich, K., Reichman-Fried, M., Grimaldi, C., and Raz, E. (2015). Chemokine-dependent pH elevation at the cell front sustains polarity in directionally migrating zebrafish germ cells. Curr. Biol. 25 (8), 1096–1103. doi:10.1016/j.cub.2015.02.071
Van Remmen, H., Williams, M. D., Guo, Z., Estlack, L., Yang, H., Carlson, E. J., et al. (2001). Knockout mice heterozygous for Sod2 show alterations in cardiac mitochondrial function and apoptosis. Am. J. Physiol. Heart Circ. Physiol. 281 (3), H1422–H1432. doi:10.1152/ajpheart.2001.281.3.H1422
Velarde, M. C., Flynn, J. M., Day, N. U., Melov, S., and Campisi, J. (2012). Mitochondrial oxidative stress caused by Sod2 deficiency promotes cellular senescence and aging phenotypes in the skin. Aging 4 (1), 3–12. doi:10.18632/aging.100423
Wei, W., Tuna, S., Keogh, M. J., Smith, K. R., Aitman, T. J., Beales, P. L., et al. (2019). Germline selection shapes human mitochondrial DNA diversity. Science 364 (6442), 6520. doi:10.1126/science.aau6520
Weidinger, G., Wolke, U., Köprunner, M., Thisse, C., Thisse, B., and Raz, E. (2002). Regulation of zebrafish primordial germ cell migration by attraction towards an intermediate target. Development 129 (1), 25–36. doi:10.1242/dev.129.1.25
Weyemi, U., Parekh, P. R., Redon, C. E., and Bonner, W. M. (2012). SOD2 deficiency promotes aging phenotypes in mouse skin. Aging 4 (2), 116–118. doi:10.18632/aging.100433
Yamaguchi, N., Zhang, Z., Schneider, T., Wang, B., Panozzo, D., and Knaut, H. (2022). Rear traction forces drive adherent tissue migration in vivo. Nat. Cell Biol. 24 (2), 194–204. doi:10.1038/s41556-022-00844-9
Zhang, R., Tu, Y., Ye, D., Gu, Z., Chen, Z., and Sun, Y. (2022). A germline-specific regulator of mitochondrial fusion is required for maintenance and differentiation of germline stem and progenitor cells. Adv. Sci. 9 (36), 2203631. doi:10.1002/advs.202203631
Keywords: SOD2, mitochondria, germ cells, zebrafish, cell competition, cell migration, PLLp, posterior lateral line primordium
Citation: Tarbashevich K, Ermlich L, Wegner J, Pfeiffer J and Raz E (2023) The mitochondrial protein Sod2 is important for the migration, maintenance, and fitness of germ cells. Front. Cell Dev. Biol. 11:1250643. doi: 10.3389/fcell.2023.1250643
Received: 30 June 2023; Accepted: 16 October 2023;
Published: 26 October 2023.
Edited by:
Emilia Bramanti, National Research Council (CNR), ItalyReviewed by:
Wang Wang, University of Washington, United StatesMarta Marzullo, Sapienza University of Rome, Italy
Copyright © 2023 Tarbashevich, Ermlich, Wegner, Pfeiffer and Raz. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Erez Raz, erez.raz@uni-muenster.de