- 1Department of Genome Editing, Institute of Biomedical Science, Kansai Medical University, Osaka, Japan
- 2Research Institute for Microbial Diseases, Osaka University, Osaka, Japan
- 3Graduate School of Pharmaceutical Sciences, Osaka University, Osaka, Japan
- 4Precursory Research for Embryonic Science and Technology (PRESTO), Japan Science and Technology Agency, Tokyo, Japan
- 5National Institute for Basic Biology, National Institutes of Natural Sciences, Okazaki, Japan
- 6Center for Infectious Disease Education and Research, Osaka University, Osaka, Japan
- 7The Institute of Medical Science, The University of Tokyo, Tokyo, Japan
Since the advent of gene-targeting technology in embryonic stem cells, mice have become a primary model organism for investigating human gene function due to the striking genomic similarities between the two species. With the introduction of the CRISPR/Cas9 system for genome editing in mice, the pace of loss-of-function analysis has accelerated significantly. This has led to the identification of numerous genes that play crucial roles in male reproductive processes, including meiosis, chromatin condensation, flagellum formation in the testis, sperm maturation in the epididymis, and fertilization in the oviduct. Despite the advancements, the functions of many genes, particularly those enriched in male reproductive tissues, remain largely unknown. In our study, we focused on 15 genes and generated 13 gene-deficient mice [4933411K16Rik, Adam triple (Adam20, Adam25, and Adam39), BC048671, Cfap68, Gm4846, Gm4984, Gm13570, Nt5c1b, Ppp1r42, Saxo4, Sh3d21, Spz1, and Tektl1] to elucidate their roles in male fertility. Surprisingly, all 13 gene-deficient mice exhibited normal fertility in natural breeding experiments, indicating that these genes are not essential for male fertility. These findings have important implications as they may help prevent other research laboratories from duplicating efforts to generate knockout mice for genes that do not demonstrate an apparent phenotype related to male fertility. By shedding light on the dispensability of these genes, our study contributes to a more efficient allocation of research resources in the exploration of male reproductive biology.
1 Introduction
Approximately half of all infertility cases stem from male factors, affecting one in six couples globally, as reported by the World Health Organization (Agarwal et al., 2021). Male infertility, recognized as a complex clinical disorder with a wide range of phenotypic manifestations, has traditionally been categorized into pre-testicular, testicular, and post-testicular factors. However, a newer clinical classification now divides male infertility into four main categories: reproductive tract obstruction or dysfunction, hypothalamic-pituitary axis dysfunction, semen abnormalities, and reduced sperm quality and/or quantity (Tournaye et al., 2017; Krausz and Riera-Escamilla, 2018). Due to this multifaceted nature, diagnosing male infertility proves challenging, with the underlying causes of reproductive abnormalities remaining unclear in nearly half of cases.
In mammals, male fertility hinges on the consistent production of motile sperm capable of fertilizing an egg. Within testes, spermatogonial stem cells differentiate into spermatocytes, which then undergo meiosis to become haploid spermatids. These spermatids undergo significant changes, including acrosome formation, chromatin condensation, mitochondrial rearrangement, and flagellum formation, to develop into testicular spermatozoa (Yogo, 2022). Following this, sperm undergo additional maturation while traveling through the epididymis, acquiring the capability to reach the oviduct where they encounter ovulated eggs. Sperm passes through cumulus cells surrounding zona pellucida; only acrosome-reacted sperm can penetrate this layer. Finally, sperm fuse with the egg’s plasma membrane for fertilization (Ikawa et al., 2010). In testis, previous reports suggest that 1000–2000 germ cell-enriched genes may have the potential to regulate and facilitate these processes. (Schultz et al., 2003). Publications have discussed the expression of these genes, yet their in vivo functions in both humans and mice remain largely unknown. This is primarily attributed to the cost associated with knock-out techniques. However, ongoing exploration of whole gene functions in mice has been improving this situation (Miyata et al., 2016; Fujihara et al., 2018; Okabe, 2018).
Our research focuses on uncovering functions of genes enriched in the male reproductive system, elucidating their roles through detailed studies in gene-deficient mice (Miyata et al., 2016; Lu et al., 2019; Noda et al., 2019; Park et al., 2020; Sun et al., 2020). Here, we present loss-of-function analyses of 15 such genes (4933411K16Rik, Adam20, Adam25, Adam39, BC048671, Cfap68, Gm4846, Gm4984, Gm13570, Nt5c1b, Ppp1r42, Saxo4, Sh3d21, Spz1, and Tektl1), revealing that each is dispensable for male fertility.
2 Materials and methods
2.1 Animals
The animal experiments conducted in this study received approval from the Institutional Animal Care and Use Committees of Osaka University (Osaka, Japan) and Kansai Medical University (Osaka, Japan), adhering strictly to guidelines and regulations for animal experiments. B6D2F1 and ICR mice utilized in this research were purchased from CLEA Japan, Inc. (Tokyo, Japan), Japan SLC, Inc. (Shizuoka, Japan), or Shimizu Laboratory Supplies (Kyoto, Japan). In this study, all gene-modified mice generated on a hybrid genetic background of B6D2F1 and these mutant mice will be accessible to other researchers via either the RIKEN BioResource Research Center, Ibaraki, Japan, or the Center for Animal Resources and Development (CARD), Kumamoto University, Kumamoto, Japan (Supplementary Table S1).
2.2 Digital PCR
We generated digital PCRs (heatmaps) illustrating the average transcript per million values across a range of mouse and human reproductive and non-reproductive tissues, employing methodologies consistent with our previously published work (Robertson et al., 2020). In summary, sequences corresponding to various tissues were retrieved from Sequence Read Archives, subjected to trimming via TrimGalore, and aligned against the human genome (GRCh38) or mouse genome (GRCm38) utilizing HISAT2 (Kim et al., 2015). Gene expression levels in each tissue were quantified using Feature Counts, while RUVr was utilized for batch correction, eliminating unwanted variation (Risso et al., 2014). Differential gene expression between each non-reproductive tissue and reproductive tissue was determined using EdgeR (Robinson et al., 2010). To visualize the expression heatmap, we employed software accessible at https://orit.research.bcm.edu/MRGDv2.
2.3 Reverse transcription PCR (RT-PCR)
Total RNA from various tissues of wild-type mice were extracted using the TRIzol reagent (Thermo Fisher Scientific, Waltham, MA, United States) according to the standard procedure. Total RNA was reverse transcribed into cDNA using the SuperScript IV First-Strand Synthesis System (Thermo Fisher Scientific) and deoxyribonuclease (RT Grade) (Nippon Gene). PCR primers for each gene are listed in Supplementary Table S2.
2.4 Generation of KO mice with the CRISPR/Cas9 system
All gene-deficient mouse lines in this study were produced using the CRISPR/Cas9 system. CRISPRdirect software was used to find off-target sequences when we selected guide RNAs (Supplementary Table S3) (Naito et al., 2015). To prevent truncated protein expression and disrupt the functional role of proteins derived from mutated genes, we deliberately targeted the nearly complete coding regions of the gene locus. To generate gene-deficient mice, we performed electroporation using zygotes as described previously (Abbasi et al., 2018; Noda et al., 2019). The crRNA and tracrRNA (Integrated DNA Technologies, Coralville, IA, United States or Merck, Darmstadt, Germany) were mixed and hybridized by heating to 95°C and cooling to room temperature in the duplex buffer (Integrated DNA Technologies). The Cas9 protein (Thermo Fisher Scientific or Integrated DNA Technologies) was mixed with hybridized RNAs in Opti-MEM (Thermo Fisher Scientific). This solution was incubated at 37°C to generate the gRNA/Cas9 ribonucleoprotein (RNP) complex, and the resulting complex was electroporated into fertilized eggs using NEPA21 Super Electroporator (Nepagene, Chiba, Japan). The treated embryos that developed to the 2-cell stage were transplanted into the oviducts of pseudopregnant ICR females at 0.5 days after mating with vasectomized males. Pups (Founder 0) were obtained by natural or cesarean section and F0 mice that harbor each mutation are mated to wild-type B6D2F1 mice from a commercial source to produce heterozygous F1 offspring. Subsequently, F1 heterozygotes are interbred to produce F2 homozygous mutated offsprings. All experiments were conducted using F2 or F3-generation mice for analysis. Genotyping was performed by PCR using the primers listed in Supplementary Table S4. To examine the deleted region in the mutated allele of each line, we purified the amplified PCR products using the respective primers that detected the mutated allele. Subsequently, we subjected them to Sanger sequencing.
2.5 Morphological and histological analysis of testes
After measuring the whole body and testicular weight, testes were fixed in Bouin’s solution (Polysciences, Warrington, PA, United States or FUJIFILM Wako, Osaka, Japan), embedded in paraffin, sectioned at a thickness of 5 μm on a Microm HM325 microtome (Microm, Walldorf, Germany) or RM2125RT microtome (Leica, Nussloch, Germany), rehydrated, and treated with 1% periodic acid for 10 min. Schiff’s reagent (FUJIFILM Wako) was then applied for 20 min. The Periodic acid-Shiff staining sections were then stained with Mayer’s hematoxylin solution (FUJIFILM Wako). Testicular sections of Adam triple- and Sh3d21-deficient mice were stained with Mayer’s hematoxylin and EosinY solution (FUJIFILM Wako) before imaging. These sections were observed using an Olympus BX53 microscope equipped with an Olympus DP74 color camera (Olympus, Tokyo, Japan) or BZ-9000 (Keyence, Osaka, Japan).
2.6 Analysis of morphology and motility of spermatozoa
Spermatozoa from cauda epididymis were suspended in TYH medium (Muro et al., 2016), incubated at 37°C for either 10–30 min or 120 min, diluted, and then placed on MAS-coated glass slides (Matsunami Glass, Osaka, Japan) for morphology assessment using an Olympus BX53 microscope or placed in Leja glass chambers (imv Technologies, Orne, France) for sperm motility analysis using the CEROS II sperm analysis system (software version 1.5; Hamilton Thorne Biosciences, Beverly, MA, United States) or SMAS (DITECT, Tokyo, Japan).
2.7 Sperm-egg fusion assay
The fusion assay was performed as previously described (Inoue et al., 2005). Specifically, cauda epididymal spermatozoa were incubated in TYH drops for 2 h before the fusion assay. Unfertilized eggs were collected from hormone-treated females 14 h after hCG injection. To remove the ZP, eggs were treated with 1 mg/mL collagenase (038–22361, FUJIFILM Wako) for 30 min. After washing in TYH, eggs were stained with 1 μg/mL Hoechst 33,342 (DOJINDO, Kumamoto, Japan) for 10 min and washed repeatedly in TYH drops without Hoechst dye. For insemination, 2 × 105 spermatozoa per mL were used. After incubation for 30 min, the eggs were fixed with 0.25% glutaraldehyde (FUJIFILM Wako), and we counted the percentage of fused eggs.
2.8 Fertility analysis of KO lines
Sexually mature control or null male mice were individually housed with two to three 8-week-old wild-type B6D2F1 female mice in the same cage for at least 8 weeks (except for Adam20, Adam25, Adam39-triple and Sh3d21-deficient mice). To ensure statistical validity, three or four males were tested for each knockout line. At the end of the mating period, male mice were removed from the cages, and the wild-type females were kept for an additional 3 weeks to allow for the delivery of their final litters. During the fertility test, the number of pups and copulation plugs were counted weekday mornings.
2.9 Statistical analyses
Statistical difference was determined using the Student’s t-test by Microsoft Office Excel (Microsoft Corporation, Redmond, WA, United States). Differences were considered statically significant if the p < 0.05. Data represent the mean ± standard deviation (s.d.).
3 Results
3.1 In silico expression analyses of 15 reproductive organ-enriched genes
We attempted to confirm the expression patterns of 13 genes of interest by digital PCR using Mammalian Reproductive Genetics Database V2 (Figures 1A,B). As shown in Figures 1A,B, all genes except 4933411K16Rik, Cfap68, Saxo4, and Sh3d21 demonstrated predominant expression in the testes or epididymides of mice and human. 4933411K16Rik, Cfap68, Saxo4, and Sh3d21 exhibited high expression levels in the testis and epididymis, with slightly weaker signals detected in somatic tissues. In a study by Benoît Vanderperre et al., additional members of the mitochondrial pyruvate carrier (MPC) family in placental mammals were identified using MPC1 and MPC2 protein sequences as input for analysis via the tblastn algorithm (Vanderperre et al., 2016). The additional MPC family genes, Gm4984 and Gm13570, were not identified in these databases. We conducted reverse transcription-polymerase chain reaction (RT-PCR) analysis using mouse tissues, revealing that somatic-type genes encoding a mitochondrial pyruvate carrier protein, Mpc1, and Mpc2, were ubiquitously expressed. Meanwhile, Gm4984 and Gm13570, predicted to be MPC family genes, exhibited specific expression in the testis (Figure 1C).
Figure 1. In silico analyses of the expression patterns of specific genes in multiple tissues and spermatogenic cells. (A) Digital PCR indicates the expression patterns of 13 genes of interest in mouse tissues. The average transcript per million (TPM) value per tissue per gene was generated from 77 published mouse RNA-seq data sets and 147 human RNA-seq data sets, respectively. Only Cfap68 is detected slightly weak signals in ubiquitous tissues and highly expressed in testis and epididymis. The other 12 genes showed predominant or restricted expression in mouse testes or epididymides. White = 0 TPM, Black ≥30 TPM. Eif3l is used as an internal control. (B) Digital PCR indicates the expression patterns of 10 genes of interest in human tissues. White = 0 TPM, Black ≥30 TPM. GAPDH is used as an internal control. ADAM25, ADAM39, and GM4846 are not found in human genomes. (C) RT-PCR indicates the expression patterns of Mpc1, Mpc2, Gm4984, and Gm13570 in mouse tissues. Actb is used as an internal control.
3.2 Production and phenotypic analyses of gene-deficient mice
Individual gene-deficient mouse lines were generated by CRISPR/Cas9 system, and phenotypic analyses were conducted to investigate the essentiality of these genes in male reproduction. The analyses of Adam20, Adam25, Adam39-triple, Sh3d21, Gm4984, and Gm13570 deficient males will be described herein, whereas those of 4933411K16Rik, BC048671, Cfap68, Gm4846, Nt5c1b, Ppp1r42, Saxo4, Spz1, and Tektl1 deficient males are presented in Table 1 and Supplementary Figures S1–S3.
Table 1. Results of the fertility tests for the 13 mutant mouse lines. Each wild-type (WT) or mutant males were individually caged with WT female mice. NE indicates not examined.
In 2017, Yan-Wei Sha et al., 2018 reported that the heterozygous variant in ADAM20 mutation (NM_003814:exon2:c.641A>C:p.D214A) was found in an infertile male patient in China, who was concluded to be associated with sperm-egg fusion disorder. In mice, Adam25 and Adam39, which are adjacent to the Adam20 locus, are considered homologs of human ADAM20, as neither gene is identified in the human genome. These proteins might have a redundant role in mice. To assess the function of Adam20 in mice, we established Adam20, Adam25, and Adam39 triple deficient mice. Adam20, Adam25, and Adam39 are positioned consecutively on the forward strand of mouse chromosome 8. These genes encode a single-pass type I membrane protein respectively and contain metalloproteinase, disintegrin, and EGF-like domains. Using two guide RNAs targeting the 5′ region of Adam25 and the 3′ region of Adam39, we obtained mutant mice with a 78,899 bp deletion containing the coding region of the three genes (Figures 2A,B). Control and Adam triple null males displayed normal testicular appearance and comparable testis weights relative to body weights (Figure 2C). Histological analyses revealed no abnormalities in the testes and epididymides of control and Adam triple null males (Figure 2D). Additionally, both groups exhibited normal sperm morphology (Figure 2E). Despite reports in humans suggesting the importance of ADAM20 in sperm-egg fusion (Sha et al., 2018), our analysis using Adam triple deficient sperm revealed normal fusion ability compared to normal sperm (Figure 2F). Moreover, Adam triple deficient mice showed normal litter size mating with wild-type female mice (Figure 2G).
Figure 2. Phenotypic analysis of Adam20, Adam25, and Adam39 triple (Adam-triple) -deficient male mice. (A) Genomic structure and knockout strategy of Adam20, Adam25, and Adam39. Two guide RNAs were designed to target the upstream of exon 1 of Adam25 and downstream of exon 2 of Adam39. Eight primers (Fw1, 2, 3, 4 and Rv1, 2, 3, 4) were designed for genotyping the mutant mice. (B) Mutant and wild-type alleles were detected by genomic PCR using primer sets Fw1-Rv1 (mutant allele), Fw2-Rv2 (Adam25 WT allele), Fw3-Rv3 (Adam20 WT allele), and Fw4-Rv4 (Adam39 WT allele). DNA sequence of the mutant allele by Sanger sequencing was shown in the lower panel. (C) Testis appearance and testis to body weight ratios of control (Wild type) and Adam-triple homozygous mutated mice. Scale bar = 3 mm. (D) Histological analysis of testes and epididymides in control and Adam-triple deficient mice. Scale bar = 50 μm. (E) Morphology of cauda epididymal spermatozoa in control and Adam-triple deficient mice. Scale bar = 20 μm. (F) Analysis of sperm and egg fusion ability. (G) Average litter size of control and Adam triple-deficient male mice. NS indicates not significant.
Protein phosphorylation plays a crucial role in germ cell development in mammals, spanning from spermatogonia to spermatids, and including key processes like sperm capacitation, motility, the acrosome reaction, and fertilization (Zhang et al., 2023). Several protein kinases, such as mitogen-activated protein kinases (MAPKs), aurora kinase C (AURKC), and Src family kinases (SFKs), are integral to various stages of spermatogenesis and are essential for mammalian fertility (Li et al., 2009; Coutton et al., 2015; Xiao et al., 2019). SH3D21 is identified as a candidate gene in male infertility, linked to a maturation arrest of spermatogenesis in infertile patient groups (Stouffs et al., 2012). Additionally, it has been characterized as a protein that might have effects on other metabolic disorders and a sensitizer for gemcitabine in the treatment of pancreatic cancer (Hellwege et al., 2017; Masoudi et al., 2019; Tan et al., 2022). SH3D21, bearing the SRC homology 3 (SH3) domain, is found in signaling proteins and serves as a classical protein-protein interaction module. SH3 domain-containing proteins regulate a multitude of processes including cell-environment and cell-cell communication, cytoskeletal rearrangement, cell migration and growth, protein trafficking and degradation, and immune responses (Kaneko et al., 2008; Tatarova et al., 2012; Mayer, 2015). However, the specific functions of SH3D21 in male reproductive systems remain poorly understood. This protein is encoded by a gene located on the reverse strand of mouse chromosome 4. Similarly, Sh3d21 deficient mice were produced using two guide RNAs targeting the region containing between exon 1 and exon 14, and we obtained mutant offspring bearing a 12,262 bp deletion were identified by genomic PCR and Sanger sequencing (Figures 3A,B). Phenotypic analyses demonstrated that Sh3d21 null males displayed normal testicular appearance and weights (Figure 3C), as well as normal testicular and epididymal histology (Figure 3D). Moreover, they exhibited normal sperm morphology and motility (Figures 3E,F) and maintained normal in vivo fertility (Figure 3G).
Figure 3. Phenotypic analysis of Sh3d21 deficient male mice. (A) Genomic structure and knockout strategy of Sh3d21. Two guide RNAs were designed to target the first coding exon (Exon 1) and Exon 14. Four primers (Fw1, 2 and Rv1, 2) were designed for genotyping the mutant mice. (B) Mutant and wild-type alleles were detected by genomic PCR using primer sets Fw1-Rv1 (mutant allele) and Fw2-Rv2 (WT allele). DNA sequence and chromatograph of the mutant allele by Sanger sequencing were shown in the lower panel. (C) Testis appearance and testis to body weight ratios of Sh3d21 heterozygous and homozygous mutated mice. Scale bar = 3 mm. (D) Histological analysis of testes and epididymides in control and Sh3d21 deficient mice. Scale bar = 50 μm. (E) Morphology of cauda epididymal spermatozoa in control and Sh3d21 deficient mice. Scale bar = 20 μm. (F) Proportion of motile sperm and progressively motile sperm after 120 min of incubation using computer-assisted sperm analysis (CASA). (G) Average litter size of control and Sh3d21 deficient male mice. NS indicates not significant.
Mitochondria play a vital role in cellular energy production, generating ATP and metabolites essential for cell survival and growth. Pyruvate, a key metabolite, serves as a major source of acetyl-CoA required for the tricarboxylic acid (TCA) cycle, a critical component of aerobic respiration (Martínez-Reyes and Chandel, 2020). Acetyl-CoA derived from pyruvate also contributes to synthesizing lipids and amino acids, supporting anabolic processes. Thus, understanding the molecular mechanisms governing pyruvate fate within cells is crucial for comprehending cellular metabolism regulation and overall cell function and fate. Central to this process is the transport of pyruvate from the cytosol to the mitochondrial matrix. It is believed that cytosolic pyruvate transport to the outer mitochondrial membrane via voltage-dependent anion channel/porin channels (Palmieri et al., 2011) and the inner mitochondrial membrane via a specific transporter called the mitochondrial pyruvate carrier (MPC) complexes (Bricker et al., 2012; Herzig et al., 2012). Gm4984 (MPC1L-mitochondrial pyruvate carrier 1 like) and Gm13570 (MPCLP- mitochondrial pyruvate carrier like protein) were identified as candidates of MPC1 family based on protein sequence similarity and are expressed at very high concentrations in mice, particularly in the mitochondria of adult reproductive tissues (Vanderperre et al., 2016). Gm4984, located as a single-exon gene on the forward strand of mouse chromosome X, and Gm13570, also a single-exon gene positioned on the reverse strand of mouse chromosome 2, are predicted to encode mitochondria-localized proteins with a role as pyruvate carriers. (Vanderperre et al., 2016) (Figure 4A). Using two guide RNAs in each gene that flanking the single exon, we introduced a 1,558 bp deletion (+2 bp insertion) into the locus of Gm4984 and a 1,193 bp deletion (+4 bp insertion) into the locus of Gm13570 (Figures 4A,B). Gm4984 null mice displayed slightly higher ratio compare to control and Gm13570 null mice displayed normal testis weights (+/+: 3.2 ± 0.3, Gm4984 null: 3.7 ± 0.4, Gm13570 null: 3.0 ± 0.3, Figure 4C). Both Gm4984 null and Gm13570 null mice displayed normal testicular histology (Figure 4D). Furthermore, both null mice exhibited normal sperm morphology (Figure 4E), with computer-assisted sperm analysis indicating comparable sperm motility and progressive motility (Figure 4F). In terms of sperm swimming velocity comparison between control and null males, there was a slight decrease in the VAP and VCL parameters in Gm4984 null at the 10 min incubation. However, after 120 min of incubation, which induced capacitation, there were no significant differences observed in any parameter (Figure 4G). Importantly, both null mice maintained normal in vivo fertility and slightly increased the litter size in Gm4984 null mice (+/+: 8.4 ± 2.0, Gm4984 null: 10.7 ± 1.9, Gm13570 null: 8.9 ± 2.6, Figure 4H).
Figure 4. Phenotypic analysis of Gm4984 and Gm13570 deficient male mice. (A) Genomic structure and knockout strategy of Gm4984 and Gm13570. Two guide RNAs were designed respectively to target the first coding exon (Exon 1) of both genes. Four primers (Fw1, 2 and Rv1, 2) were designed for genotyping the mutant mice. (B) Mutant and wild-type alleles were detected by genomic PCR using primer sets Fw1-Rv1 (wild-type and mutant allele of Gm4984) and Fw2-Rv2 (wild-type and mutant allele of Gm13570). DNA sequence and chromatograph of the mutant allele by Sanger sequencing were shown in the right panel. (C) Testis to body weight ratios of control, Gm4984, and Gm13570 deficient mice. (D) Histological analysis of testes in control, Gm4984, and Gm13570 deficient mice. Scale bar = 50 μm. (E) Morphology of cauda epididymal spermatozoa in control, Gm4984, and Gm13570 deficient mice. Scale bar = 20 μm. (F) Proportion of motile sperm and progressively motile sperm after 10 and 120 min of incubation using computer-assisted sperm analysis (CASA). (G) Motility of cauda epididymal spermatozoa in control, Gm4984, and Gm13570 deficient mice. Sperm motility and kinetic parameters were measured after incubation in TYH media for 10 and 120 min. VAP, average path velocity; VSL, straight-line velocity; VCL, curvilinear velocity. (H) Average litter size of control, Gm4984, and Gm13570 deficient male mice. The reduction in litter size within each deficient line isn't notably significant when compared to the control group. * indicates p < 0.05 and NS indicates not significant.
3.3 Fertility results
Through housing the gene-deficient males individually with wild-type females, we discovered that these males were able to sire offspring with average litter sizes comparable to those sired by wild-type males (Table 1). Hence, all mutant males examined in this study displayed normal fertility.
4 Discussion
About 15 out of 100 couples have trouble getting pregnant. Roughly half of infertility cases are attributed to males and not fully understood, and there may be thousands of genes that can be harbored for this issue (Esteves, 2013; Lotti and Maggi, 2018). The application of the CRISPR-Cas9 system has helped clarify certain genes related to male reproduction (Miyata et al., 2016; Lu et al., 2019; Park et al., 2020; Sun et al., 2020). In this study, we continue utilizing CRISPR/Cas9 to analyze the role of these 15 genes in male reproduction. All of the genes under analysis exhibited higher expression levels in the human and mouse testes or epididymis. Despite the limited expression of some of these genes in other tissues, we did not observe any aberrant behavior or appearance in the mutant mice. This suggests that these genes may not be important for human male reproduction either.
The A Disintegrin and Metalloproteinases (ADAMs) belong to the metalloproteinase superfamily, which comprises a diverse group of multi-domain transmembrane and secreted proteins with varied biological functions such as inflammation and cancers, and reproductive functions (Cho, 2012; Arai et al., 2023). In 2017, Yan-Wei Sha et al., 2018 identified a heterozygous variant (NM_003814:exon2:c.641A>C:p.D214A) in ADAM20 in an infertile male patient in China, associated with a sperm-egg fusion disorder. Adam25 and Adam39, neighboring Adam20 locus in mice, are considered human ADAM20 homologs. To confirm Adam20 in vivo function in mammals, we developed triple deficient mice lacking Adam20, Adam25, and Adam39. Surprisingly, the Adam triple deficient mice showed normal spermatogenesis, sperm morphology, fusion ability, and in vivo fertility compared to control. In mice, many ADAM proteins that are specifically expressed in reproductive tissues were identified and remained unknown function in mammalian fertility (Cho, 2012). These proteins might work as a key molecule in the fusion ability of sperm or compensate for the function of ADAM20 in mice. The advancements in next-generation sequencing have revolutionized genome characterization, facilitating precise identification of human genome mutations that may underlie disease. However, validating whether gene deficiency or abnormalities cause the disease in vivo is essential. Utilizing gene-deficient mice models is a potent strategy for pinpointing disease-causing genes. These models are pivotal in unraveling critical functions, not only in reproductive biology but also in other areas.
Protein phosphorylation is pivotal in mammalian germ cell development, encompassing processes from spermatogonia to spermatids, including sperm capacitation, motility, the acrosome reaction, and fertilization (Zhang et al., 2023). SH3D21 emerges as a candidate gene in male infertility, associated with spermatogenesis maturation arrest in infertile patients (Stouffs et al., 2012). However, its specific roles in male reproductive systems remain unclear. Digital PCR analyses indicate high expression of Sh3d21 in the testis, implying potential significant roles during spermatogenesis. Surprisingly, Sh3d21 deficient mice displayed normal spermatogenesis and in vivo fertility. The SRC homology 3 family encompasses approximately 300 annotated proteins containing SH3 domains in humans (Mayer, 2015), suggesting that other family genes enriched in the testis might compensate for the function of SH3D21.
Gm4984 (MPC1L-mitochondrial pyruvate carrier 1 like) and Gm13570 (MPCLP- mitochondrial pyruvate carrier like protein) were identified as MPC1 family candidates based on protein sequence similarity. They are highly expressed in mice, particularly in adult reproductive tissue mitochondria (Vanderperre et al., 2016). While GM4984 in humans and mice are highly conserved in the first 100 amino acids and are exclusively expressed in round and elongated haploid spermatids during their differentiation into mature sperm cells (Vanderperre et al., 2016), our observations did not reveal any abnormalities in the male reproductive system of both gene-deficient mice. Testicular germ cells and mature sperm possess MPC1 and MPC2, which are somatic types of mitochondrial pyruvate carrier complex subunits. It is plausible that the functions of these genes are sufficient to complete the fertilization process under natural conditions, despite the potential importance of GM4984 and GM13570 in the male reproductive system.
PPP1R42, localized at the basal body primary cilia in retinal epithelial cells, has been implicated in maintaining centrosome number. Inhibition of PPP1R42 expression results in increased centrosome number, leading to mis-segregation and cell death (DeVaul et al., 2013). Cfap68 gene encodes cilia and flagella-associated protein 68, which is involved in the structural specialization of the sperm tail (Leung et al., 2023). SAXO4 protein is a member of the stabilizer of axonemal microtubules (SAXO) group which decorates sperm doublet microtubules and likely stabilizes the tubulin lattice during sperm movement. TEKTL1 (CCDC105) which is a tektin-like sperm-microtubule inner protein, has been described as a conserved, sperm-specific, tektin-like molecule that forms a filament in the cleft between protofilaments and binds several core- and sperm-microtubule inner proteins. It may help stabilize the outer junction which suggests that deletion of Tektl1 may lead to fertility problems. These mice deficient in each protein showed normal spermatogenesis, sperm morphology, and in vivo fertility. Utilizing in situ cryoelectron tomography (Cryo-ET), three research groups have identified proteins associated with microtubule doublets (DMTs) and microtubule inner proteins (MIPs) within sperm flagella (Chen et al., 2023; Leung et al., 2023; Tai et al., 2023). Among these, Tekt5 is identified as a gene expressed specifically in the testis and is expected that play essential roles in DMTs structure. Surprisingly, Tekt5 null mice exhibit normal fertility despite an increased ratio of sperm with defective flagella (Chen et al., 2023). Cryo-ET analysis of Tekt5 null sperm reveals low occupancies in 3-helix bundles and this analysis may reveal abnormalities in the flagellum of Tektl1 null sperm. These findings suggest that while structural biology provides intricate insights into biomolecular architecture, functional studies are indispensable for unraveling their biological significance and contributions to physiological processes and diseases. This synergy between structural and functional approaches is vital for a comprehensive understanding of biology. Other similar proteins, CFAP proteins, SAXO proteins, and TEKTIN proteins were identified in sperm flagellar components in sea urchin and bovine and these proteins might have duplicate roles and compensate for the function of deficient proteins (Leung et al., 2023). While Cfap68 exhibits weak expression in various tissues such as cerebral cortex, kidney, lung, and small intestine of mice, its deficiency does not seem to lead to obvious abnormalities or lethality. Delving deeper into the analysis of these tissues could potentially unveil additional functions of Cfap68 beyond what is currently investigated in this study.
Similarly, we identified four other relatively unknown genes—4933411K16Rik, BC048671, N15c1b, and Spz1 that are highly expressed in the testis and Gm4846 specifically expressed in the epididymis (Hammarström et al., 1995; Hsu et al., 2001; Kiyozumi et al., 2020; Kiyozumi et al., 2023; Kiyozumi, 2024). Despite their significant expression patterns, our study revealed that these genes are not essential for male fertility, as no differences were observed between control and gene-deficient mice in terms of fertility.
In summary, our study utilized the CRISPR/Cas9 system to disrupt 15 genes enriched in reproductive tissues in mice. Through in vivo functional analysis of 13 established gene-deficient mouse lines, we demonstrated that none of these genes are essential for normal male fertility. These gene-deficient mouse lines will be deposited at public institutions in Japan, facilitating future research to explore overlapping gene functions in mammals and understand the mechanisms underlying male reproductive function. The evolution of the CRISPR/Cas9 system has enabled the creation of gene-deficient animals in various experimental species, offering valuable insights into reproductive biology (Lee et al., 2018). Furthermore, exploring proteins like acrosin may reveal significant implications in other animal species, driving forward our understanding of reproductive processes and unveiling novel perspectives (Hirose et al., 2020).
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The animal study was approved by the Institutional Animal Care and Use Committees of Osaka University The Institutional Animal Care and Use Committees of Kansai Medical University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
TN: Writing–original draft, Writing–review and editing. KT: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing. KS: Data curation, Funding acquisition, Investigation, Writing–review and editing. HW: Data curation, Investigation, Writing–review and editing. DM: Data curation, Investigation, Writing–review and editing. ST: Data curation, Investigation, Writing–review and editing. DK: Data curation, Funding acquisition, Investigation, Writing–review and editing. MI: Conceptualization, Funding acquisition, Project administration, Supervision, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by Ministry of Education, Culture, Sports, Science and Technology/Japan Society for the Promotion of Science KAKENHI (Grants-in-Aid for Scientific Research) Grants JP21H02487 and JP21K19263 to DK, JP23K05831 to KS, JP20K06476 to KT, and JP25112007, JP17H01394, and JP19H05750 to MI; Japan Science and Technology Agency Grant JPMJPR2143 to DK, Japan Agency for Medical Research and Development Grant JP19gm5010001 to MI, grants from Kansai Medical University Molecular Imaging Center of Diseases to KT and Takeda Science Foundation grants to KT and MI.
Acknowledgments
We thank Saki Nishioka, Kotone Kawamura, Kaito Yamamoto, Yonggang Lu, and Chen Pan (Research Institue for Microbial Diseases, Osaka University), and the Biotechnology Research and Development (nonprofit organization) for excellent technical assistance.
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.2024.1411162/full#supplementary-material
References
Abbasi, F., Miyata, H., Shimada, K., Morohoshi, A., Nozawa, K., Matsumura, T., et al. (2018). RSPH6A is required for sperm flagellum formation and male fertility in mice. J. Cell Sci. 131 (19), jcs221648. doi:10.1242/jcs.221648
Agarwal, A., Baskaran, S., Parekh, N., Cho, C. L., Henkel, R., Vij, S., et al. (2021). Male infertility. Lancet 397 (10271), 319–333. doi:10.1016/S0140-6736(20)32667-2
Arai, J., Otoyama, Y., Nozawa, H., Kato, N., and Yoshida, H. (2023). The immunological role of ADAMs in the field of gastroenterological chronic inflammatory diseases and cancers: a review. Oncogene 42 (8), 549–558. doi:10.1038/s41388-022-02583-5
Bricker, D. K., Taylor, E. B., Schell, J. C., Orsak, T., Boutron, A., Chen, Y. C., et al. (2012). A mitochondrial pyruvate carrier required for pyruvate uptake in yeast, Drosophila, and humans. Science 337 (6090), 96–100. doi:10.1126/science.1218099
Chen, Z., Shiozaki, M., Haas, K. M., Skinner, W. M., Zhao, S., Guo, C., et al. (2023). De novo protein identification in mammalian sperm using in situ cryoelectron tomography and AlphaFold2 docking. Cell 186 (23), 5041–5053.e19. doi:10.1016/j.cell.2023.09.017
Cho, C. (2012). Testicular and epididymal ADAMs: expression and function during fertilization. Nat. Rev. Urol. 9 (10), 550–560. doi:10.1038/nrurol.2012.167
Coutton, C., Escoffier, J., Martinez, G., Arnoult, C., and Ray, P. F. (2015). Teratozoospermia: spotlight on the main genetic actors in the human. Hum. Reprod. Update 21 (4), 455–485. doi:10.1093/humupd/dmv020
DeVaul, N., Wang, R., and Sperry, A. O. (2013). PPP1R42, a PP1 binding protein, regulates centrosome dynamics in ARPE-19 cells. Biol. Cell 105 (8), 359–371. doi:10.1111/boc.201300019
Esteves, S. C. (2013). A clinical appraisal of the genetic basis in unexplained male infertility. J. Hum. Reprod. Sci. 6 (3), 176–182. doi:10.4103/0974-1208.121419
Fujihara, Y., Miyata, H., and Ikawa, M. (2018). Factors controlling sperm migration through the oviduct revealed by gene-modified mouse models. Exp. Anim. 67 (2), 91–104. doi:10.1538/expanim.17-0153
Hammarström, L., Samuelsson, J., and Grimfors, G. (1995). Subcutaneous gammaglobulin for patients with secondary hypogammaglobulinaemia. Lancet. 345 (8946), 382–383. doi:10.1016/s0140-6736(95)90367-4
Hellwege, J. N., Velez Edwards, D. R., Acra, S., Chen, K., Buchowski, M. S., and Edwards, T. L. (2017). Association of gene coding variation and resting metabolic rate in a multi-ethnic sample of children and adults. BMC Obes. 4, 12. doi:10.1186/s40608-017-0145-5
Herzig, S., Raemy, E., Montessuit, S., Veuthey, J. L., Zamboni, N., Westermann, B., et al. (2012). Identification and functional expression of the mitochondrial pyruvate carrier. Science 337 (6090), 93–96. doi:10.1126/science.1218530
Hirose, M., Honda, A., Fulka, H., Tamura-Nakano, M., Matoba, S., Tomishima, T., et al. (2020). Acrosin is essential for sperm penetration through the zona pellucida in hamsters. Proc. Natl. Acad. Sci. U. S. A. 117 (5), 2513–2518. doi:10.1073/pnas.1917595117
Hsu, S. H., Shyu, H. W., Hsieh-Li, H. M., and Li, H. (2001). Spz1, a novel bHLH-Zip protein, is specifically expressed in testis. Mech. Dev. 100 (2), 177–187. doi:10.1016/s0925-4773(00)00513-x
Ikawa, M., Inoue, N., Benham, A. M., and Okabe, M. (2010). Fertilization: a sperm's journey to and interaction with the oocyte. J. Clin. Invest. 120 (4), 984–994. doi:10.1172/JCI41585
Inoue, N., Ikawa, M., Isotani, A., and Okabe, M. (2005). The immunoglobulin superfamily protein Izumo is required for sperm to fuse with eggs. Nature 434 (7030), 234–238. doi:10.1038/nature03362
Kaneko, T., Li, L., and Li, S. S. (2008). The SH3 domain--a family of versatile peptide- and protein-recognition module. Front. Biosci. 13, 4938–4952. doi:10.2741/3053
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12 (4), 357–360. doi:10.1038/nmeth.3317
Kiyozumi, D. (2024). Busulfan administration replicated the characteristics of the epididymal initial segment observed in mice lacking testis-epididymis lumicrine signaling. J. Reprod. Dev. 70, 104–114. doi:10.1262/jrd.2023-102
Kiyozumi, D., Noda, T., Yamaguchi, R., Tobita, T., Matsumura, T., Shimada, K., et al. (2020). NELL2-mediated lumicrine signaling through OVCH2 is required for male fertility. Science 368 (6495), 1132–1135. doi:10.1126/science.aay5134
Kiyozumi, D., Shimada, K., Chalick, M., Emori, C., Kodani, M., Oura, S., et al. (2023). A small secreted protein NICOL regulates lumicrine-mediated sperm maturation and male fertility. Nat. Commun. 14 (1), 2354. doi:10.1038/s41467-023-37984-x
Krausz, C., and Riera-Escamilla, A. (2018). Genetics of male infertility. Nat. Rev. Urol. 15 (6), 369–384. doi:10.1038/s41585-018-0003-3
Lee, J. G., Sung, Y. H., and Baek, I. J. (2018). Generation of genetically-engineered animals using engineered endonucleases. Arch. Pharm. Res. 41 (9), 885–897. doi:10.1007/s12272-018-1037-z
Leung, M. R., Zeng, J., Wang, X., Roelofs, M. C., Huang, W., Zenezini Chiozzi, R., et al. (2023). Structural specializations of the sperm tail. Cell 186 (13), 2880–2896.e17. doi:10.1016/j.cell.2023.05.026
Li, M. W., Mruk, D. D., and Cheng, C. Y. (2009). Mitogen-activated protein kinases in male reproductive function. Trends Mol. Med. 15 (4), 159–168. doi:10.1016/j.molmed.2009.02.002
Lotti, F., and Maggi, M. (2018). Sexual dysfunction and male infertility. Nat. Rev. Urol. 15 (5), 287–307. doi:10.1038/nrurol.2018.20
Lu, Y., Oura, S., Matsumura, T., Oji, A., Sakurai, N., Fujihara, Y., et al. (2019). CRISPR/Cas9-mediated genome editing reveals 30 testis-enriched genes dispensable for male fertility in mice†. Biol. Reprod. 101 (2), 501–511. doi:10.1093/biolre/ioz103
Martínez-Reyes, I., and Chandel, N. S. (2020). Mitochondrial TCA cycle metabolites control physiology and disease. Nat. Commun. 11 (1), 102. doi:10.1038/s41467-019-13668-3
Masoudi, M., Seki, M., Yazdanparast, R., Yachie, N., and Aburatani, H. (2019). A genome-scale CRISPR/Cas9 knockout screening reveals SH3D21 as a sensitizer for gemcitabine. Sci. Rep. 9 (1), 19188. doi:10.1038/s41598-019-55893-2
Mayer, B. J. (2015). The discovery of modular binding domains: building blocks of cell signalling. Nat. Rev. Mol. Cell Biol. 16 (11), 691–698. doi:10.1038/nrm4068
Miyata, H., Castaneda, J. M., Fujihara, Y., Yu, Z., Archambeault, D. R., Isotani, A., et al. (2016). Genome engineering uncovers 54 evolutionarily conserved and testis-enriched genes that are not required for male fertility in mice. Proc. Natl. Acad. Sci. U. S. A. 113 (28), 7704–7710. doi:10.1073/pnas.1608458113
Muro, Y., Hasuwa, H., Isotani, A., Miyata, H., Yamagata, K., Ikawa, M., et al. (2016). Behavior of mouse spermatozoa in the female reproductive tract from soon after mating to the beginning of fertilization. Biol. Reprod. 94 (4), 80. doi:10.1095/biolreprod.115.135368
Naito, Y., Hino, K., Bono, H., and Ui-Tei, K. (2015). CRISPRdirect: software for designing CRISPR/Cas guide RNA with reduced off-target sites. Bioinformatics 31 (7), 1120–1123. doi:10.1093/bioinformatics/btu743
Noda, T., Sakurai, N., Nozawa, K., Kobayashi, S., Devlin, D. J., Matzuk, M. M., et al. (2019). Nine genes abundantly expressed in the epididymis are not essential for male fecundity in mice. Andrology 7 (5), 644–653. doi:10.1111/andr.12621
Okabe, M. (2018). Sperm-egg interaction and fertilization: past, present, and future. Biol. Reprod. 99 (1), 134–146. doi:10.1093/biolre/ioy028
Palmieri, F., Pierri, C. L., De Grassi, A., Nunes-Nesi, A., and Fernie, A. R. (2011). Evolution, structure and function of mitochondrial carriers: a review with new insights. Plant J. 66 (1), 161–181. doi:10.1111/j.1365-313X.2011.04516.x
Park, S., Shimada, K., Fujihara, Y., Xu, Z., Larasati, T., Pratiwi, P., et al. (2020). CRISPR/Cas9-mediated genome-edited mice reveal 10 testis-enriched genes are dispensable for male fecundity. Biol. Reprod. 103 (2), 195–204. doi:10.1093/biolre/ioaa084
Risso, D., Ngai, J., Speed, T. P., and Dudoit, S. (2014). Normalization of RNA-seq data using factor analysis of control genes or samples. Nat. Biotechnol. 32 (9), 896–902. doi:10.1038/nbt.2931
Robertson, M. J., Kent, K., Tharp, N., Nozawa, K., Dean, L., Mathew, M., et al. (2020). Large-scale discovery of male reproductive tract-specific genes through analysis of RNA-seq datasets. BMC Biol. 18 (1), 103. doi:10.1186/s12915-020-00826-z
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26 (1), 139–140. doi:10.1093/bioinformatics/btp616
Schultz, N., Hamra, F. K., and Garbers, D. L. (2003). A multitude of genes expressed solely in meiotic or postmeiotic spermatogenic cells offers a myriad of contraceptive targets. Proc. Natl. Acad. Sci. U. S. A. 100 (21), 12201–12206. doi:10.1073/pnas.1635054100
Sha, Y. W., Xu, X., Ji, Z. Y., Mei, L. B., Qiu, P. P., Ji, H., et al. (2018). Sperm-egg fusion disorder in a Chinese male patient was associated with a rare ADAM20 variant. Oncotarget 9 (2), 2086–2091. doi:10.18632/oncotarget.23331
Stouffs, K., Vandermaelen, D., Massart, A., Menten, B., Vergult, S., Tournaye, H., et al. (2012). Array comparative genomic hybridization in male infertility. Hum. Reprod. 27 (3), 921–929. doi:10.1093/humrep/der440
Sun, J., Lu, Y., Nozawa, K., Xu, Z., Morohoshi, A., Castaneda, J. M., et al. (2020). CRISPR/Cas9-based genome editing in mice uncovers 13 testis- or epididymis-enriched genes individually dispensable for male reproduction†. Biol. Reprod. 103 (2), 183–194. doi:10.1093/biolre/ioaa083
Tai, L., Yin, G., Huang, X., Sun, F., and Zhu, Y. (2023). In-cell structural insight into the stability of sperm microtubule doublet. Cell Discov. 9 (1), 116. doi:10.1038/s41421-023-00606-3
Tan, X., Zhu, T., Zhang, L., Fu, L., Hu, Y., Li, H., et al. (2022). miR-669a-5p promotes adipogenic differentiation and induces browning in preadipocytes. Adipocyte 11 (1), 120–132. doi:10.1080/21623945.2022.2030570
Tatarova, Z., Brabek, J., Rosel, D., and Novotny, M. (2012). SH3 domain tyrosine phosphorylation--sites, role and evolution. PLoS One 7 (5), e36310. doi:10.1371/journal.pone.0036310
Tournaye, H., Krausz, C., and Oates, R. D. (2017). Novel concepts in the aetiology of male reproductive impairment. Lancet Diabetes Endocrinol. 5 (7), 544–553. doi:10.1016/S2213-8587(16)30040-7
Vanderperre, B., Cermakova, K., Escoffier, J., Kaba, M., Bender, T., Nef, S., et al. (2016). MPC1-like is a placental mammal-specific mitochondrial pyruvate carrier subunit expressed in postmeiotic male germ cells. J. Biol. Chem. 291 (32), 16448–16461. doi:10.1074/jbc.M116.733840
Xiao, X., Yang, Y., Mao, B., Cheng, C. Y., and Ni, Y. (2019). Emerging role for SRC family kinases in junction dynamics during spermatogenesis. Reproduction 157 (3), R85–R94. doi:10.1530/REP-18-0440
Yogo, K. (2022). Molecular basis of the morphogenesis of sperm head and tail in mice. Reprod. Med. Biol. 21 (1), e12466. doi:10.1002/rmb2.12466
Keywords: CRISPR/Cas9, knockout mice, male infertility, spermatozoa, testis
Citation: Nguyen TTT, Tokuhiro K, Shimada K, Wang H, Mashiko D, Tonai S, Kiyozumi D and Ikawa M (2024) Gene-deficient mouse model established by CRISPR/Cas9 system reveals 15 reproductive organ-enriched genes dispensable for male fertility. Front. Cell Dev. Biol. 12:1411162. doi: 10.3389/fcell.2024.1411162
Received: 03 April 2024; Accepted: 02 May 2024;
Published: 21 May 2024.
Edited by:
Vanina Gabriela Da Ros, CONICET Institute of Biology and Experimental Medicine (IBYME), ArgentinaReviewed by:
Hubert Schorle, University Hospital Bonn, GermanySu-Ren Chen, Beijing Normal University, China
Brendan J. Houston, The University of Melbourne, Australia
Copyright © 2024 Nguyen, Tokuhiro, Shimada, Wang, Mashiko, Tonai, Kiyozumi and Ikawa. 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: Keizo Tokuhiro, tokuhirk@hirakata.kmu.ac.jp; Masahito Ikawa, ikawa@biken.osaka-u.ac.jp
†These authors have contributed equally to this work