1. Introduction
Recent studies have shown that up to 10-15% of human couples in the world are currently infertile, and about 30%-40% of these cases are due to male infertility factors, including asthenospermia and oligospermia [
1]. Multiple genetic manipulation tools in
Drosophila melanogaster provide a significant opportunity to study conserved signaling pathways during spermatogenesis [
2]. The spermatogenesis of
Drosophila mainly includes germline stem cells renewal, primary spermatocyte meiosis and spermatocyte-to-sperm morphological transformation [
3]. In the spermatocyte, each germline stem cells (GSC) divides asymmetrically to create two daughter cells, one daughter cells attached to the central cells to maintain its stem cells characteristics, and the other daughter cells, called the spermatogonium (GB), enters the differentiation process, and one GB undergoes four mitosis to create 16 spermatogonia cells, which grow into the primary spermatocytes (16S). It then undergoes two meiosis to create 64 round sperm cells connected by intercytoplasmic bridges. The sperm cells in the cyst undergo lengthening and individualization, eventually becoming haploid mature sperm, after which they coil and enter the spermatocyst until fertilization occurs [
4].
The balance between animal germ cells proliferation and differentiation is critical for germline homeostasis. It is tightly regulated by stem cells niche and spermatogonial transit-amplification (TA) division [
5]. Germline differentiation must be coordinated with GSCs self-renewal, and spermatogonium must undergo meiosis at the right time [
6]. Previous studies have shown that transcriptional activation is the primary mechanism of mitosis to meiosis transition, and the expression of multiple genes is changed during the transition from S16 to primary spermatocyte [
7]. Activation of tMAC (testis-specific meiotic arrest complex) promotes spermatocyte meiosis [
8,
9]. The mutation of
aly gene, as one of the tMAC members, could result in a rich spermatocyte in the testes and abnormal spermatocyte differentiation, thus eventually male infertility [
10].
Recent studies have shown significant mitochondrial recombination during TA transformation in
Drosophila spermatogenesis [
11]. In the early stages of spermatogenesis in many insects, the number and appearance of mitochondria are similar to those of somatic cells, dispersed in the cytoplasm. In the primary spermatocyte before meiosis, the mitochondria cluster near the nucleus and present a spherical structure. After meiosis II is completed, hundreds of mitochondria fuse into two long mitochondrial derivatives, which intertwine to form a spherical aggregate near the nucleus, known as the accessory nucleus (nebenkern). During the elongation of sperm cells, the behavior of mitochondrial derivatives also changes, in which the contents of major derivatives are relatively dense, while the volume of minor derivatives decreases and there is no accumulation of contents [
12].
Drosophila fzo was first discovered as a mitochondrial dynamic gene, belongs to one member of the evolutionally-conserved mitochondrial fusion protein family (Mfns), encoding a sizeable transmembrane GTPase associated with dynamin that mediates mitochondrial fusion during spermatogenesis. Deletion of the
Fzo (fuzzy onions) gene causes the mitochondria to fail to fuse properly, resulting in the formation of malformed nebenkern and male sterility [
13].
The human gene
VARS2 (
Valyl-tRNA synthetase 2) is related to
COXPD20 (
combined oxidative phosphorylation deficiency 20) [
14]. Studies have shown that
VARS2 mutation can cause autosomal recessive mitochondrial encephalomyopathy [
15,
16].
VARS2 mutations can also lead to varying degrees of developmental delay, axial hypotonia, limb spasms, and even premature death [
17]. The homolog of
VARS in flies is
ValRS-m gene, located on mitochondria, and highly expressed in the testis. Notably, we have previously found that the knockdown of
ocnus which encoding histidine phosphatase in fly testis can lead to abnormal development of testis and male sterility [
18]. Through comparative proteomics detection, ValRS-m protein was significantly down-regulated in
ocn knock-down fly testes. However, its specific function in
Drosophila has not been reported so far.
In the present study, we explored the role of ValRS-m in Drosophila spermatogenesis. ValRS-m knockdown in testes caused male sterility and early spermatogenesis defects in files. ValRS-m-knockdown testes exhibited severe defects in spermatocyte differentiation. Immunofluorescence staining and electron microscopic examination revealed ValRS-m-knockdown disrupts the mitochondrial fusion and ATP synthesis, excessive apoptosis occurs in ValRS-m-knockdown testes. Together, these findings indicated that ValRS-m is required in Drosophila spermatogenesis by controlling mitochondrial dynamics and function in spermatocyte differentiation. Our results can not only further understand the regulatory mechanism of animal reproductive development, but also provide a valuable theoretical reference for the further study of human spermatogenesis and mitochondrial diseases.
3. Discussion
Spermatogenesis is a complex process of cell differentiation, and the transition of cell state is regulated by a strict cascade [
6]. The transition process from spermatogonia to primary spermatocyte is very essential in
Drosophila spermatogenesis, in which large numbers of genes are involved along with transcriptional and morphological changes, and about 1500 genes are estimated to be expressed only in spermatocytes by genome-wide microarray data [
31,
32]. In this study, we found that the disappearance of both sperm bundles and mature sperm in the spermatogenesis have been observed in one day old
ValRS-m RNAi testes, implying the role of
ValRS-m in regulating spermatogenesis and testes development. Especially, we characterized
ValRS-m as a crucial factor in regulating the transition from primary spermatocytes to spermatid in the
Drosophila testes, through immunofluorescence staining of Vasa and Spectrin (
Figure 2).
ValRS-m gene is localized in mitochondria, which was the powerhouses of eukaryotic cells due to their ability to produce ATP through oxidative phosphorylation [
17,
33]. ATP production is essential for sperm motility, capacitation and acrosomal reactions [
34]. In somatic cells, as cells enter interphase, mitochondria exhibit an elongated network pattern and aggregate in the nucleus and pericellular. In contrast, during most mitosis, mitochondria exhibit a fragmented network pattern scattered throughout the cytoplasm [
35,
36]. The development and meiosis of germ cells are much more complex than somatic cells. For example, sperm and eggs undergo meiosis and the complex dynamic distribution of the cytoskeleton. These processes require large amounts of ATP to maintain [
37,
38]. In this study, we found that ATP signaling was very powerful in control testes and highly expressed in both round and elongated spermatids. However, ATP5A signaling was almost undetectable in
ValRS-m knocked-down testes (
Figure 3). Therefore, these results suggest that
ValRS-m knockdown leads to the obstruction of ATP synthesis in testes.
Studies have shown that mitochondria can participate in a series of activities in sperm formation, including the fusion and aggregation of early active mitochondria to form a nebenkern [
39]. Dorogova et al (2013) showed that the primary spermatocytes of wild-type
Drosophila underwent a pre-meiosis growth phase, and the number of mitochondria increased significantly, the mitochondria of mature spermatocytes are compact and evenly distributed in the cytoplasm [
40]. Through transmission electron microscopy, we speculated that the spermatogenesis of
ValRS-m knock-down testes may be stagnant in the primary spermatocyte stage without subsequent meiosis, because the axial and fibrous complex is missing and the nucleus is in a diffuse state in the absence of
ValRS-m.
This prediction was further confirmed by subsequent RNA-seq and qRT-PCR results. The RNA-seq analysis and qRT-PCR revealed that many meiosis-related genes were significantly decreased in
ValRS-m RNAi testes compared to the control testes including
Cbc and
Tsen54 (
Figure 6B). The expression of
Cbc gene encodes a polynucleotide 5'-hydroxyl-kinase required at the transition to meiosis in spermatogenesis, as well as the expression of
Tsen54 is conducive to the transition to meiosis [
30]. At the same time, the tMAC complex helps to activate meiosis of spermatocytes [
41]. QRT-PCR results also showed that these tMAC complex encoding genes, including
Aly,
Comr,
Topi, and
Wuc expression level were down-regulated significantly in the lack of
ValRS-m. Due to
Aly,
Comr and
Topi are only expressed in primary spermatocytes [
10,
32,
41], implying that the meiosis is suppressed on account of knockdown of
ValRS-m through inhibiting the transcription of these meiosis related genes in the primary spermatocyte, hence resulting in the stalling of spermatogenesis. Moreover, these spermatocytes that did not enter meiosis eventually underwent apoptosis, which is consistent with the increase of apoptotic signal tunnel observed (
Figure 4).
Furthermore, we found that
ValRS-m RNAi interrupts the expression of multiple mitochondrial genes (
Figure 6A). These mitochondrial genes play a critical role in mitochondrial function and spermatogenesis. For example, Cyt-c1, as a subunit of mitochondrial respiratory chain complex III, is involved in oxidative phosphorylation mediating electron transport from cytochrome b to cytochrome c [
42,
43].
Sprn is specifically located in the mitochondrial nebenskern, playing a role in sperm elongation [
44].
Mics1 encodes mitochondrial intima proteins and interacts with proteins encoded by
chchd2 to enhance oxidative phosphorylation [
45]. Studies have shown that Knon, an unusually large paralogic homologue of the ATP synthase subunit d, participates in the internal shaping of the nebenskern forming mitochondrial membrane during spermatogenesis in
Drosophila. Male flies lacking
knon are sterile [
46].
Cyt-c1L encodes a subunit of mitochondrial respiratory chain complex III, involved in mitochondrial ATP synthesis and coupled proton transport process. We previously found that the knockdown of
Cyt-c1L could result in abnormal axial filament, unable to form individual complex, and lead to abnormal spermatocyte apoptosis and immature sperm, resulting in male infertility [
29]. Taken together, these studies explore that the deficiency of
ValRS-m might inhibit spermatogenesis by affecting mitochondrial structure and function in flies.
The effect of
ValRS-m RNAi on
Drosophila spermatogenesis is extensive and significant. To further reveal its regulatory mechanism on male fly fertility, we have identified 7710 DEGs in
ValRS-m RNAi testes by RNA-Seq, including 4725 down-regulated genes and 2985 up-regulated genes. We found that the metabolic pathways associated with up-regulated DEGs were most involved in the ribosome pathway (104 genes), such as RPS3 (Ribosomal protein S3), RPL22 (Ribosomal protein L22), RPL36 (Ribosomal protein L36), etc. Studies have shown that most ribosomal proteins not only play a role in the biological reaction process of ribosomes and protein synthesis, but also affect the process of biological somatic activities by performing "extra-ribosomal" functions, which are involved in many biological processes, including growth and development, cell apoptosis and aging processes [
47]. The fruit fly ribosome contains 79 different proteins, among which RPS3 is not only involved in translation and ribosome maturation, but also can recognize DNA damage [
48]. Importantly, RpS3 plays a critical role in regulating spermatid elongation and individualization [
49]. In addition, the deletion of
RpL22 gene is also common in many types of cancer and cell lines [
50], and
RPL36 knock-down results in a reduced number of germ cells in the testes and few or no mature sperm in the seminal vesicles [
51].
Additionally, the down-regulated DEGs is most enriched in glycolysis/gluconeogenesis pathway and pyruvate metabolic pathway, and is also involved in citric acid cycle and pentose phosphate metabolic pathway, which are all carbohydrate metabolic pathways. Glycolysis is a vital component of semen activity because it is involved in capacitation, provides a second energy source for spermatocyte development. It is closely associated with sperm functional maturation and epididymal transport [
52]. Studies have shown that carbohydrate metabolism is closely related to spermatogenesis [
53,
54]. For example, the testicles promote carbohydrate metabolism in adjacent gut sites through JAK-STAT signaling in fruit flies [
10,
55]. The male intestine secretes citrate to the adjacent sperms, thus promoting sperm maturation [
53]. We found that
ValRS-m knockdown significantly inhibited the expression level of genes involved in carbohydrate synthesis, thus it may have an impact on the energy required for spermatogenesis. Therefore, it is speculated that the knockdown of
ValRS-m may affect the function of various ribosomal proteins and carbohydrate metabolism hence destroy the spermatogenesis.
4. Materials and Methods
4.1. Fly Stocks
Flies were reared on a standard cornmeal/yeast diet at 25°C and under non-crowded conditions (200 ± 10 eggs per 50-ml vial of media in the 150 ml conical flask). Transgenic ValRS-m-hp D. melanogaster line was obtained from the Tsinghua Fly Center (Beijing). The bamGal4 vp16 line was a kind gift from Professor Zhaohui Wang at the Institute of Genetics and Developmental Biology, Chinese Academy of Sciences. All flies were raised in the 150 mL conical flasks with standard corn/sugar medium at 25℃.
4.2. Fertility Test
Gene knockdown flies (bamGal4>ValRS-m-hp) were generated by crossing transgenic RNAi males with virgin bamGal4 females. Flies from the cross between w1118 males and bamGal4 females were used as the control (bamGal4>+). The 1-day-old gene knockdown males were crossed with 3-4-day-old w1118 females to assay male fertility. For each biological repeat, 15 males and 10 females were used. After 12 h of crossing, males were removed, and then eggs were collected and incubated at 25℃ for 24 h. Hatch rates were determined from the proportion of hatched eggs to total eggs. Three biological repeats per cross type were performed.
4.3. Immunofluorescence Staining
One-day-old fly testes were dissected in PBS, fixed for 30 min in 4% paraformaldehyde at room temperature, and washed in wash buffer [phosphate buffered saline (PBS)/0.1% Triton X-100] three times for 15 min each. The testis samples were blocked for 30 min in 5% normal goat serum, incubated overnight at 4°C with primary antibodies. Primary antibodies were used at the following dilutions: rabbit anti-Vasa (1:100, cat. no. AB760351, Developmental Studies Hybridoma Bank, Iowa, IA, USA), mouse anti-spectrin (1:100, cat.no. AB528473, Developmental Studies Hybridoma Bank), anti-ATP5A. Secondary antibodies were used at the following dilutions: phalloidin (1:200, cat. no. BMD00084, Abbkine), rabbit 594 and mouse 488 (1:200, cat. no. A23420 and A23210, Abbkine). All samples were mounted on a glass slide with 4’-6-diamidino-2-phenylindole (DAPI) (2 μg ml−1, cat. no. S2110, Solarbio, Beijing, China) solution [
56]. Fluorescence images were collected using a Leica SP8 laser confocal microscope (Leica, Germany).
4.4. Transmission Electron Microscopy
The testes were dissected from 1d bamGal4/ValRS-m-hp and bamGal4>+ males in PBS, and fixed in 2.5% glutaraldehyde (0.2 M phosphate buffer, pH 7.4) at 4°C overnight. Then these testes were washed in phosphate buffer and post-fixed in 1% OsO4 for one hour. After double fixation, the samples were dehydrated through an ascending series of ethanol (30%, 50%, 70%, 80%, 100%, ten minutes for each concentration and doubled for 100%, and then embedded in Araldite (EMbed 812, China). Ultrathin sections (80 nm) were obtained in Leica Ultracut 7. The sections were placed on copper grids, stained with 2% uranyl acetate for 15 min, rinsed twice with H2O for five minutes, then stained with lead citrate for 15 min. The stained samples were observed using H-8100 transmission electron microscope (Hitachi, Tokyo, Japan) operating at 100 kV.
4.5. TUNEL Staining
TUNEL staining (TdT-mediated dUTP nick end labeling) was performed as follows in this study. Thirty pairs of 1-day-old Drosophila testis were dissected and fixed in 4% paraformaldehyde solution at RT for 40 min, then washed four times (15 minutes each time) in PBST (phosphate buffer, PBS + 0.1% Triton X-100). The sample was then incubated with the TUNEL reaction mixture (5 μl TdT enzyme solution and 45 μl fluorescein-dUTP tag solution) in a dark environment at 37°C for three hours (Roche, Germany). After rinsing three times with PBST in the dark, the testes used an anti-fading medium (Solarbio, Beijing, China) containing 2 μg/ml DAPI to fix on the glass slide. The Leica SP5 laser confocal microscope was used to observe and take photos.
4.6. RNA Extraction, Library Preparation and Sequencing
To screen the genes related to ValRS-m during spermatogenesis, we collected 80-100 testes samples from bamGal4/ValRS-m-hp and bamGal4>+ males. These testes were used for transcriptomic analysis. Six cDNA libraries were constructed and sequenced on the Illumina Novaseq 6000 platform (San Diego, USA) following the standard protocols set by Majorbio Bio-Pharm Technology Co. Ltd. (Shanghai, China). Subsequently, clean reads were mapped to the genome of Drosophila melanogaster (Flybase, dmelr6.15 genome) using TopHat (http:// tophat. cbcb.umd. edu/) software. After assembly of the mapped reads, the unigenes were run against a non-redundant database.
To identify DEGs (differential expression genes) between two different groups, the expression level of each gene was calculated according to the transcripts per million reads (TPM) method. RSEM (
http://deweylab.biostat.wisc.edu/rsem/) was used to quantify gene abundances [
57]. DEGs were determined with a cutoff of |log2 fold change (RNAi/control) |≥ 1 and p-value ≤ 0.05, as the criteria used in the previous study [
58]. In addition, functional-enrichment analysis including GO (Gene Ontology,
http://www.geneontology.org) and KEGG (Kyoto Encyclopedia of Genes and Genomes,
http://www.genome.jp/kegg/) were performed to identify which DEGs were significantly enriched in GO terms and metabolic pathways at P value≤0.05 compared with the whole-transcriptome background. GO functional enrichment and KEGG pathway analysis were carried out by Goatools and KOBAS [
59].
4.7. Quantitative Real-Time PCR
The assay was performed as previously described [
18]. Total RNA was extracted using Trizol (cat. no. 15596026, Invitrogen, Waltham, MA, USA) from testes and ovaries of 1-day- old adult w1118, as well as testes from bamGal4/ValRS-m-hp and bamGal4>+ males. First-strand cDNA was synthesized from 2 μg of total RNA using EasyScript first-strand cDNA synthesis SuperMix Kit (cat. no. AT321-01, TransGen Biotech, Beijing, China). qPCR was performed with specific primers (Table S2) using a Miniopticon system (Bio-Rad, Hercules, CA, USA) with a Platinum SYBR Green qPCR SuperMix (cat. no. Q711-02, Vazyme, Shanghai, China) as described previously (Zheng et al., 2018). The qPCR cycling programmed was as follows: 95°C for 3 min, followed by 40 cycles of 95°C for 10 s, 55–60°C (depending on different primers) for 30 s and 65°C for 5 s. Then, a melting curve was constructed from 55°C to 98°C. The relative expression of the gene was calibrated against the reference gene
rp49 using the 2−ΔΔCT calculation method: ΔΔCt=(Ct
Target−Ct
rp49)
ValRS-m RNAi−(Ct
Target−Ct
rp49)control [
60].
4.8. Statistical Analysis
Statistical analysis was performed using GraphPad Prism software (GraphPad Inc., La Jolla, CA, USA). Student’s test was used for two-group comparisons. A value of P< 0.05 was considered significantly different. The statistical details were indicated in the corresponding figure legends.
Figure 1.
Knockdown of ValRS-m in testes caused Drosophila male infertility. (A) The ValRS-m expression levels in one day old adult testes was significantly increased than that in ovaries by qRT-PCR. (B) The ValRS-m expression level was significantly reduced in Bam Gal4>ValRS-m-hp (ValRS-m RNAi) testes. (C) The hatch rate of eggs from the cross of wild-type females and ValRS-m RNAi males was zero. *P<0.05; ** P<0.01; *** P<0.001. The immunofluorescent staining using DAPI for spermatogenesis in control testis (D-D’’) and ValRS-m-knockdown testis (E-E’’). D and E are the whole testes. D’ and E' are the areas where sperm clusters (indicated by arrowheads) are located. D’’ and E’’ are the seminal vesicles (SV). D’ and E’ are parts of the enlarged area of D and E. the yellow boxes represent sperm clusters. DNA was stained with DAPI (blue). Many mature acicular sperm nuclei are visible in the SV of the control group (D’’), but no sperm nuclei are visible in the SV of the knockdown testes (E’’). DNA was stained with DAPI (blue). Scale bars, 100 μm (D, E); 50 μm (D’, D’’, E’, E’’).
Figure 1.
Knockdown of ValRS-m in testes caused Drosophila male infertility. (A) The ValRS-m expression levels in one day old adult testes was significantly increased than that in ovaries by qRT-PCR. (B) The ValRS-m expression level was significantly reduced in Bam Gal4>ValRS-m-hp (ValRS-m RNAi) testes. (C) The hatch rate of eggs from the cross of wild-type females and ValRS-m RNAi males was zero. *P<0.05; ** P<0.01; *** P<0.001. The immunofluorescent staining using DAPI for spermatogenesis in control testis (D-D’’) and ValRS-m-knockdown testis (E-E’’). D and E are the whole testes. D’ and E' are the areas where sperm clusters (indicated by arrowheads) are located. D’’ and E’’ are the seminal vesicles (SV). D’ and E’ are parts of the enlarged area of D and E. the yellow boxes represent sperm clusters. DNA was stained with DAPI (blue). Many mature acicular sperm nuclei are visible in the SV of the control group (D’’), but no sperm nuclei are visible in the SV of the knockdown testes (E’’). DNA was stained with DAPI (blue). Scale bars, 100 μm (D, E); 50 μm (D’, D’’, E’, E’’).
Figure 2.
Reduction of ValRS-m restrained germline differentiation. (A-A’’): Immunostaining of Vasa (red) and DAPI (blue) in testes from control bamGal4>+ males and bamGal4>ValRS-m-hp males. Vasa is a pan-marker for germ cells (A’’, B’’). The yellow boxes represent sperm clusters (A’). The white boxes represent germ cell cysts (A’’, B’’). The antibody against α-spectrin labels the fusomes (A1, B1), which exhibited two shapes: the round fusomes (indicated by white arrowhead) and branched fusomes (yellow arrowheads). The number of punctate fusomes (C) or branched fusomes (D) in the control bamGal4>+ (n = 10) and ValRS-m RNAi (n = 8) testes was quantified. *P<0.05; ** P<0.01. Scale bars,100 μm.
Figure 2.
Reduction of ValRS-m restrained germline differentiation. (A-A’’): Immunostaining of Vasa (red) and DAPI (blue) in testes from control bamGal4>+ males and bamGal4>ValRS-m-hp males. Vasa is a pan-marker for germ cells (A’’, B’’). The yellow boxes represent sperm clusters (A’). The white boxes represent germ cell cysts (A’’, B’’). The antibody against α-spectrin labels the fusomes (A1, B1), which exhibited two shapes: the round fusomes (indicated by white arrowhead) and branched fusomes (yellow arrowheads). The number of punctate fusomes (C) or branched fusomes (D) in the control bamGal4>+ (n = 10) and ValRS-m RNAi (n = 8) testes was quantified. *P<0.05; ** P<0.01. Scale bars,100 μm.
Figure 3.
Knockdown of ValRS-m affects mitochondrial structure and function in spermatogenesis. (A-C): Control testis. (A'-C'): ValRS-m knockdown testis. (B, B'): White boxes indicate mitochondria (Mitochondrial markers). Yellow boxes indicate mitochondrial envelope morphology (VDAC1 markers). White arrows indicate ATP signaling (ATP5A markers). The mean mitochondrial staining (Mito) signals (n = 8) in testes (D). (F-I) Transmission electron microscopy (TEM) images showing mitochondrial morphogenesis during spermatogenesis. Pre-individualization spermatids in control testes (F, G) show synchronized development of their mitochondrial derivatives in the cysts, where paracrystalline material (blue arrow in G) accumulates in the major mitochondrial derivative as a contiguous dense material at the attachment site of the axoneme (red arrow in I). A few scattered mitochondrial derivatives existed (blue arrow in I) but any axoneme were not seen in the ValRS-m RNAi testes (H, I). G and I represent magnifications of F and H respectively. ** P<0.01. Scale bars, 100 μm (A-C, A’-C’); 1 μm (F, H); 500 nm (G, I).
Figure 3.
Knockdown of ValRS-m affects mitochondrial structure and function in spermatogenesis. (A-C): Control testis. (A'-C'): ValRS-m knockdown testis. (B, B'): White boxes indicate mitochondria (Mitochondrial markers). Yellow boxes indicate mitochondrial envelope morphology (VDAC1 markers). White arrows indicate ATP signaling (ATP5A markers). The mean mitochondrial staining (Mito) signals (n = 8) in testes (D). (F-I) Transmission electron microscopy (TEM) images showing mitochondrial morphogenesis during spermatogenesis. Pre-individualization spermatids in control testes (F, G) show synchronized development of their mitochondrial derivatives in the cysts, where paracrystalline material (blue arrow in G) accumulates in the major mitochondrial derivative as a contiguous dense material at the attachment site of the axoneme (red arrow in I). A few scattered mitochondrial derivatives existed (blue arrow in I) but any axoneme were not seen in the ValRS-m RNAi testes (H, I). G and I represent magnifications of F and H respectively. ** P<0.01. Scale bars, 100 μm (A-C, A’-C’); 1 μm (F, H); 500 nm (G, I).
Figure 4.
Knockdown of ValRS-m induces apoptosis in the testes. (A) The Tunel (green) signal at the apical region of the testes from bamGal4>+ males (white arrow in A’’’), is not only present at the apical, but more concentrated in the middle and rear segments in the ValRS-m knockdown testis (B, white box in B’’’). DAPI (blue) stained the nuclei. Vasa (red) stained the germ cells. The single channel images for DAPI staining were shown in (A’, B’). The single channel images for Vasa staining were shown in (A’’, B’’). The single channel images for Tunel staining were shown in (A’’’, B’’’). Scale bars: 50 µm.
Figure 4.
Knockdown of ValRS-m induces apoptosis in the testes. (A) The Tunel (green) signal at the apical region of the testes from bamGal4>+ males (white arrow in A’’’), is not only present at the apical, but more concentrated in the middle and rear segments in the ValRS-m knockdown testis (B, white box in B’’’). DAPI (blue) stained the nuclei. Vasa (red) stained the germ cells. The single channel images for DAPI staining were shown in (A’, B’). The single channel images for Vasa staining were shown in (A’’, B’’). The single channel images for Tunel staining were shown in (A’’’, B’’’). Scale bars: 50 µm.
Figure 5.
Transcript alterations were assessed by RNA-seq after ValRS-m knockdown in testes. (A) Volcan plot of differentially expressed genes (DEGs) from comparison of the control and ValRS-m knockdown groups. (B)Tissue expression profile of differentially expressed genes after ValRS-m knockdown in testis. (C) GO enrichment analysis of up-regulative DEGs in the testes of Bam Gal4>ValRS-m-hp relative to the control. (D) GO enrichment analysis of down-regulative DEGs in the testes of Bam Gal4>ValRS-m-hp relative to the control. (E) KEGG enrichment analysis of up-DEGs in the testes of ValRS-m RNAi relative to the control. (F) KEGG enrichment analysis of down-DEGs in the testes of ValRS-m RNAi relative to the control.
Figure 5.
Transcript alterations were assessed by RNA-seq after ValRS-m knockdown in testes. (A) Volcan plot of differentially expressed genes (DEGs) from comparison of the control and ValRS-m knockdown groups. (B)Tissue expression profile of differentially expressed genes after ValRS-m knockdown in testis. (C) GO enrichment analysis of up-regulative DEGs in the testes of Bam Gal4>ValRS-m-hp relative to the control. (D) GO enrichment analysis of down-regulative DEGs in the testes of Bam Gal4>ValRS-m-hp relative to the control. (E) KEGG enrichment analysis of up-DEGs in the testes of ValRS-m RNAi relative to the control. (F) KEGG enrichment analysis of down-DEGs in the testes of ValRS-m RNAi relative to the control.
Figure 6.
RT-qPCR validation of DEGs after ValRS-m knockdown in testes. (A) DEGs associated with mitochondria. (B) DEGs related to the meiosis. *P < 0.05; **P < 0.01; ***P < 0.001.
Figure 6.
RT-qPCR validation of DEGs after ValRS-m knockdown in testes. (A) DEGs associated with mitochondria. (B) DEGs related to the meiosis. *P < 0.05; **P < 0.01; ***P < 0.001.
Table 1.
Classification of DEGs that were downregulated (fold change ≥ 2, P value < 5%) in ValRS-m -knockdown males compared with controls by qRT-PCR.
Table 1.
Classification of DEGs that were downregulated (fold change ≥ 2, P value < 5%) in ValRS-m -knockdown males compared with controls by qRT-PCR.
Annotation ID |
Gene Name |
Description |
Male meiosis-related |
CG2075 |
Aly (always early) |
the onset of spermatid differentiation |
CG13493 |
Comr (cookie monster) |
involved in spermatogenesis and transcription regulation |
CG8484 |
Topi (matotopetli) |
male meiotic division and spermatid differentiation |
CG12442 |
Wuc (Wake-up-call) |
male meiotic nuclear division |
CG5970 |
Cbc (crowded by cid) |
required at the transition to meiosis in spermatogenesis |
CG5626 |
Tsen54 (tRNA splicing endonuclease subunit 54) |
tRNA processing |
Mitochondria-related
|
CG14128 |
Sprn (Spermitin) |
testis-specific mitochondrial lumen protein |
CG1287 |
mics1 (Mitochondrial morphology and cristae structure 1) |
maintenance of mitochondrial morphology |
CG11196 |
Dic3 (Dicarboxylate carrier 3) |
mitochondrial dicarboxylate carrier |
CG8330 |
tomboy40 |
import of protein precursors into mitochondria |
CG6691 |
ttm3 (tiny tim 3) |
mitochondrial protein translocation |
CG12201 |
GC2 (Glutamate Carrier 2) |
catalyzes the transport of L-glutamate across the inner mitochondrial membrane |
CG4706 |
mAcon2 (Mitochondrial aconitase 2) |
enable 4 irons, 4 sulfur cluster binding activity |
CG30354 |
UQCR-11L (Ubiquinol-cytochrome c reductase 11 kDa subunit-like) |
mitochondrial electron transport |
CG13263 |
cyt-c-d (Cytochrome c distal) |
Electron carrier protein, sperm individualization |
CG12736 |
mEFTu2 (mitochondrial translation elongation factor Tu 2) |
bring aminoacyl-tRNA to the ribosome during the elongation phase of mRNA translation |
CG14508 |
cyt-c1L (Cytochrome c1-like) |
enable ubiquinol-cytochrome-c reductase activity |
CG8257 |
CysRS-m (Cysteinyl-tRNA synthetase, mitochondrial) |
enable ATP binding activity and cysteine-tRNA ligase activity |
CG6404 |
OXA1L (OXA1L mitochondrial inner membrane protein) |
mitochondrial proton-transporting ATP synthase complex assembly |
CG7311 |
Gpo3 (Glycerophosphate oxidase 3) |
enable glycerol-3-phosphate dehydrogenase (quinone) activity |
CG6412 |
mEFTs (mitochondrial translation elongation factor Ts) |
recharge the products of mEFTu1 and mEFTu2 with GTP |
CG14290 |
Mpc1 (Mitochondrial pyruvate carrier) |
transports pyruvate across the mitochondrial inner membrane |
CG1134 |
Mul1 (Mitochondrial E3 ubiquitin protein ligase 1) |
in the control of mitochondrial morphology by promoting mitochondrial fission |
CG6914 |
ND-B14.5AL (NADH dehydrogenase (ubiquinone) B14.5 A subunit-like) |
mitochondrial electron transport, NADH to ubiquinone |
CG5718 |
SdhAL (Succinate dehydrogenase, subunit A (flavoprotein)-like) |
mitochondrial respiratory chain complex II |
CG7813 |
Knon (knotted onions) |
Nebenkern assembly |
CG14909 |
VhaM9.7-d (Vacuolar H+ ATPase M9.7 subunit d) |
proton-transporting ATPase activity |
CG4683 |
Tengl4 (Testis EndoG-Like 4) |
active in mitochondrial inner membrane and nucleus |
CG30329 |
Vha100-3 (Vacuolar H+ ATPase 100kD subunit 3) |
ATPase binding activity |
CG12162 |
Poldip2 (Polymerase (DNA-directed), delta interacting protein 2) |
active in mitochondrial nucleoid and nucleus |
CG7615 |
Fig (fos intronic gene) |
active in mitochondrion |
CG8937 |
Hsc70-1 (Heat shock protein 70 cognate 1) |
ATP hydrolysis activity |
CG3092 |
CG3092 |
protein insertion into mitochondrial inner membrane from matrix |
CG31913 |
CG31913 |
mitochondrial cytochrome c oxidase assembly |