Marília Körbes Rockenbach1
Lucas Rosa Fraga2,3,4
Thayne Woycinck Kowalski1,4,5,6,7*
Maria Teresa Vieira Sanseverino1,6,8- 1Postgraduate Program in Genetics and Molecular Biology, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil
- 2Department of Morphological Sciences, Institute of Health Sciences, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil
- 3Postgraduate Program in Medicine: Medical Sciences, Universidade Federal do Rio Grande do Sul (UFRGS), Porto Alegre, Brazil
- 4Laboratory of Genomic Medicine, Center of Experimental Research, Hospital de Clínicas de Porto Alegre, Porto Alegre, Brazil
- 5Bioinformatics Core, Hospital de Clínicas de Porto Alegre, Porto Alegre, Brazil
- 6Medical Genetics Service, Hospital de Clínicas de Porto Alegre, Porto Alegre, Brazil
- 7Centro Universitário CESUCA, Cachoeirinha, Brazil
- 8School of Medicine, Pontifícia Universidade Catolica do Rio Grande do Sul, Porto Alegre, Brazil
Gametes are specialized cells that, at fertilization, give rise to a totipotent zygote capable of generating an entire organism. Female and male germ cells undergo meiosis to produce mature gametes; however, sex-specific events of oogenesis and spermatogenesis contribute to specific roles of gametes in reproductive issues. We investigate the differential gene expression (DGE) of meiosis-related genes in human female and male gonads and gametes in normal and pathological conditions. The transcriptome data for the DGE analysis was obtained through the Gene Expression Omnibus repository, comprising human ovary and testicle samples of the prenatal period and adulthood, additionally to male (non-obstructive azoospermia (NOA) and teratozoospermia), and female (polycystic ovary syndrome (PCOS) and advanced maternal age) reproductive conditions. Gene ontology terms related to meiosis were associated with 678 genes, of which 17 genes in common were differentially expressed between the testicle and ovary during the prenatal period and adulthood. Except for SERPINA5 and SOX9, the 17 meiosis-related genes were downregulated in the testicle during the prenatal period and upregulated in adulthood compared to the ovary. No differences were observed in the oocytes of PCOS patients; however, meiosis-related genes were differentially expressed according to the patient’s age and maturity of the oocyte. In NOA and teratozoospermia, 145 meiosis-related genes were differentially expressed in comparison to the control, including OOEP; despite no recognized role in male reproduction, OOEP was co-expressed with genes related to male fertility. Taking together, these results shed light on potential genes that might be relevant to comprehend human fertility disorders.
1 Introduction
Germ cells undergo a series of mitotic and meiotic divisions, followed by a differentiation program to produce highly specialized haploid gametes that, at fertilization, fuse and give rise to a totipotent zygote, capable of generating all somatic lineages/tissues and the next-generation of gametes (Larose et al., 2019; Farini and De Felici, 2022). The proper completion of meiosis is essential for fertility and for ensuring the normal development of the offspring (Farini and De Felici, 2022). Therefore, revealing the molecular mechanisms and pathways underlying this unique process is pivotal for understanding human reproduction and fertility, as well as helping to identify therapeutic and diagnostic targets (Larose et al., 2019; Farini and De Felici, 2022).
Infertility affects 8%–12% of reproductive-age couples worldwide (Ombelet et al., 2008), reaching 30% in some regions of the world (Mascarenhas et al., 2012). In 33% of infertility cases, the inability to have children is due to female factors alone, in 20% to male factors, in 39% to both male and female factors, and in 8% the causes remain unknown, termed as unexplained infertility (Thonneau et al., 1991). Understanding meiosis may reveal possible mechanisms underlying the cases of unexplained infertility, since proper gametogenesis is important to the formation of mature gametes and to guarantee the first events of embryo development (Larose et al., 2019).
Both female and male germ cells undergo meiosis; however, the timing of meiosis and the quantity of gametes that are produced differs in males and females (Farini and De Felici, 2022). In mammals, meiosis begins during fetal life in females and continues until menopause, whilst in males it begins at puberty and is maintained throughout life (Farini and De Felici, 2022). Understanding the molecular mechanisms behind the differences in meiosis onset and progression in males and females may contribute to the elucidation of genes involved in conditions affecting human reproduction.
The current study aimed to investigate if the gene expression profile of meiosis-related genes in human male and female gonads and gametes could be associated with human reproductive disorders. To achieve that, the expression of 678 meiosis-related genes was evaluated in normal and pathological conditions affecting human reproduction, through publicly available datasets. The interesting pattern of expression demonstrated in this study sheds light on the importance of gene expression regulation programs in human reproductive disorders and elucidates potential biomarkers in the reproductive medicine field.
2 Methods
2.1 Search for ontologies and genes related to meiosis
The entire study was performed through bioinformatic analysis of publicly available data, divided in two stages: (1) Search for genes related to meiosis, and (2) Gene expression analysis of meiosis-related genes in human gonads and gametes (Figure 1). To obtain the genes related to meiosis for secondary analysis of gene expression, an ontology-based search was conducted in the R software (v.3.6.3). Gene ontology (GO) terms related to meiosis were selected using the GOterm function from the AnnotationDBI R package (v.1.48.0) and the following keywords: “meiosis OR meiotic”, “oogenesis”, “spermatogenesis”, and “gametogenesis”. A list containing all the GO terms retrieved (Supplementary Table S1) was evaluated by two authors (MKR and TWK), independently. Only GO terms associated with human gametogenesis were included. GO terms selected by both authors were included; when one GO term was selected by only one author, the selection was discussed and a consensus was reached. The genes associated with the selected ontologies were obtained through the biomaRt package (v.2.42.0), considering the human reference genome (hg38).
FIGURE 1. Study design and main results found. A search for gene ontology (GO) terms related to the meiosis process was conducted in the R software, resulting in 281 GO terms. Of these, 208 were selected, being further associated with 678 meiosis-related genes. Six datasets were selected in the Gene Expression Omnibus (GEO) repository for differential gene expression analysis. Considering the 678 meiosis-related genes, 17 were differentially expressed in the testicle compared to the ovary in both the prenatal period and adulthood. No statistically significant differences were observed in MII oocytes of PCOS patients in comparison to the Control group, nor in GV oocytes of AMA in comparison to YMA patients. However, 48 meiosis-related genes were differentially expressed in MII oocytes of AMA in comparison to YMA patients. Besides, when comparing GV vs. MII oocytes, 258 and 209 meiosis-related genes were differentially expressed in AMA and YMA patients, respectively. In male reproductive disorders’ scenario, 238 and 290 meiosis-related genes were differentially expressed in teratozoospermic and NOA patients in comparison to the Control groups, respectively, of which 145 genes were differentially expressed in both conditions. Green up arrows mean upregulated genes, whilst red down arrows downregulated genes. wpc (weeks post-conception); PCOS (Polycystic Ovary Syndrome); AMA (advanced maternal age); YMA (young maternal age); GV (germinal vesicle); MII (mature oocytes at metaphase II stage); NOA (non-obstructive azoospermia); OA (obstructive azoospermia).
Further information regarding the selected genes functions and clinical associations were obtained through the GeneCards and Hugo Gene Nomenclature Committee (HGNC) databases (Stelzer et al., 2016; Tweedie et al., 2021), as well as through the biomaRt package (v.2.42.0), filtering by human phenotype descriptions available in the Ensembl database (Cunningham et al., 2022).
2.2 Obtainment and comparisons of transcriptome data
To analyze the expression of meiosis-related genes in human reproductive disorders, we performed differential gene expression (DGE) analysis in transcriptomic data available in the Gene Expression Omnibus (GEO) repository (NCBI, USA) (Edgar et al., 2002; Barrett et al., 2013). The keywords “oocyte”, “oogonia”, “mature oocyte”, “immature oocyte”, “germinal vesicle”, “ovary”, “oogenesis”, “gametogenesis”, “sperm” “spermatozoa”, “spermatid”, “spermatocyte”, “spermatogonia”, “spermatogenesis”, and “testes OR testicle” were used for the datasets search, filtering by Entry type (Series), organism (Homo sapiens), and study type (Expression Profile by Array or Expression Profile by Throughput Sequencing). Only studies performed in consolidated platforms, containing the raw data, experimental design and sample groups available were included. For microarrays, customized libraries without probe correspondence were also excluded. Following these criteria, six datasets were selected for the DGE analysis, comprising human samples of male and female gonads, mature and immature oocytes, and sperm, under normal and reproductive pathological conditions. The datasets are publicly available in the GEO repository, registered as GSE7307 (Hever et al., 2007), GSE116278 (Lecluze et al., 2020), GSE5850 (Wood et al., 2007), GSE164371 (Ntostis et al., 2021), GSE6872 (Platts et al., 2007), and GSE145467 (Hodžić et al., 2021). Additional information about the selected datasets for DGE analysis is available in (Supplementary Table S2).
The meiosis-related genes obtained from the gene ontology-based search (678 genes) had their expression compared between the testicle and ovary samples in five different life stages (6, 9, 11–12, 17 weeks post-conception (wpc), and adult), and in conditions affecting human reproduction. The comparison of the expression in the testicle vs. ovary was performed separately for each of the five stages. Regarding conditions affecting female fertility, the DGE of in vivo matured metaphase II (MII) oocytes in Polycystic Ovary Syndrome (PCOS) versus Control patients was analyzed, as well as in MII and germinal vesicle (GV) oocytes of advanced maternal age (41–44 years old) and young maternal age (21–26 years old) women. Concerning male reproductive disorders, the DGE was analyzed in sperm samples of Teratozoospermic (≤3% of spermatozoa displaying the ideal form) versus Normospermic patients, as well as in testicle samples of patients with non-obstructive azoospermia (NOA–impaired spermatogenesis) versus obstructive azoospermia (OA–normal spermatogenesis).
2.3 Differential gene expression analysis
DGE analysis was conducted in the R environment (v.3.6.3). Firstly, for RNA-Seq data, the alignment was performed through the Galaxy Europe server (Jalili et al., 2020), using the HISAT2 (Kim et al., 2019) alignment tool against human reference genome hg38; transcript count was performed through featureCounts tool (Liao et al., 2014). The mapping with HISAT2 and transcript count with featureCounts were based on a standardized workflow from our Bioinformatics Core, which one of the authors (TWK) is part of. The parameters for RNA-Seq data alignment and transcript count were the default ones, and the alignment rate was above 80% for all the samples analyzed. The DGE was accessed in the aligned transcriptomes using the edgeR (v.3.28.1) (Robinson et al., 2010) package. Considering microarray data, the packages affy (v.1.64.0) (Gautier et al., 2004) and limma (v.3.42.2) (Ritchie et al., 2015) were used to normalize and evaluate the DGE. RNA-Seq data was normalized through the trimmed mean of M values (TMM) and microarray data by robust multi-array average (RMA).
2.4 Systems biology analysis
In the light of the DGE results, a gene-gene co-expression network for OOEP, a member of the Subcortical Maternal Complex (SCMC), was performed in the Cytoscape software (v.3.8) using the GeneMania application (Montojo et al., 2010). Only co-expression data was considered.
2.5 Statistical analysis
Statistical analysis were performed in the R environment (v.3.6.3). The results are demonstrated as values of log2 fold-change (logFC), and adjusted p-value for false discovery rate (FDR). The logFC measures the ratio between two quantities in a base 2 logarithmic scale (for instance, a doubling is equivalent to a logFC of 1, a quadrupling to a logFC of 2, and so on). Similarly, negative logFC values mean decreased amounts. The DGE was considered significant when identifying a gene with both a logFC ≥ |1| and an adjusted p-value ≤0.05.
3 Results
3.1 Ontologies and genes related to meiosis
The search terms were associated with 281 ontologies, being selected 208 for further analysis (Supplementary Table S1); Although GO includes three independent categories: molecular function, biological process, and cellular component, our keywords only returned GO terms related to biological processes. 678 genes were related to these ontologies (Supplementary Table S3), mostly involved in mechanisms of DNA repair, chromosome segregation, and meiotic regulation. Of the 678 meiosis-related genes, 243 had clinical associations available in the Ensembl database (Supplementary Table S4), including associations to spermatogenic failure, sperm motility disorder, oocyte meiotic arrest, and oocyte maturation defect.
3.2 Differential expression of meiosis-related genes under normal and pathological human reproductive conditions
The DGE of the 678 genes related to meiosis was assessed in human male and female gonads according to the life period (in adulthood and at 6, 9, 11–12, and 17 wpc). Of the 678 genes, 17 (ADAD1, BRDT, DAZL, DDX4, GTSF1, HORMAD1, MAEL, MAJIN, OOEP, SERPINA5, SOX9, SPATA22, SYCE1, SYCP3, TDRD1, TDRD12, and TEX11) were differentially expressed in the testicle in comparison to the ovary in all the analyzed developmental periods. In the adult gonads, more meiosis-related genes were upregulated in the testicle in comparison to the ovary, whilst in all the prenatal periods the meiosis genes were mostly upregulated in the ovary. The same pattern was observed for the 17 differentially expressed genes in all developmental periods, except for SERPINA5 and SOX9, which maintain the upregulation in the testicle in both prenatal and adults periods, as expected considering their roles in male reproduction (Yang and Geiger, 2017; Selvam et al., 2019; Jiménez et al., 2021) (Figure 2). These 17 genes are mainly related to processes involved in homologous chromosome pairing, recombination, and segregation (Supplementary Table S5).
FIGURE 2. Heatmap of the differential gene expression of the 17 meiosis-related genes differentially expressed in all the developmental periods analyzed. The X-axis indicates the analyzed genes, whilst the Y-axis represents the comparisons performed. Except for SERPINA5 and SOX9, which maintain the upregulation in the testicle in both prenatal and adult periods, all the other genes were downregulated in the testicle in comparison to the ovary during the prenatal periods (6 wpc, 9 wpc, 11–12 wpc, and 17 wpc), and upregulated in the adulthood. No statistically significant differences were observed in MII oocytes of PCOS patients in comparison to the Control group, neither in MII and GV oocytes of AMA in comparison to YMA patients. However, when comparing GV vs. MII oocytes, a downregulation of OOEP and TDRD1, as well as an upregulation of GTSF1 and SYCE1, was observed in both AMA and YMA patients. A downregulation of MAJIN in GV in comparison to MII oocytes was observed only in YMA, and a downregulated of BRDT and SPATA22 in AMA. A downregulation of ADAD1, BRDT, GTSF1, HORMAD1, OOEP, SPATA22, SYCP3, and TEX11 in both teratozoospermic and NOA patients in comparison to the Control groups. A downregulation of DAZL, MAEL, and SERPINA5 was demonstrated only in teratozoospermic patients, as well as a downregulation of DDX4, MAJIN, TDRD1, and TDRD12 in NOA patients, in comparison to the Control groups. wpc (weeks post-conception); PCOS (Polycystic Ovary Syndrome); AMA (advanced maternal age); YMA (young maternal age); GV (germinal vesicle); MII (mature oocytes at metaphase II stage); NOA (non-obstructive azoospermia–impaired spermatogenesis); OA (obstructive azoospermia–normal spermatogenesis); logFC (Log2 Fold Change); * means statistically significant differences between groups (logFC ≥ |1| and adjusted p-value ≤0.05). Warm colors (red) represent upregulated genes, whilst cool colors (blue) downregulated genes.
No statistically significant differences were observed in MII oocytes of PCOS patients in comparison to the Control group, nor in GV oocytes of advanced maternal age in comparison to young maternal age patients. However, 48 meiosis-related genes were differentially expressed in MII oocytes of advanced maternal age in comparison to young maternal age patients. Besides, when comparing GV vs. MII oocytes, 258 and 209 meiosis-related genes were differentially expressed in advanced maternal age and young maternal age patients, respectively.
In view of male reproductive disorders, 238 and 290 meiosis-related genes were differentially expressed in teratozoospermic and NOA patients in comparison to the control groups, respectively, of which 145 genes were differentially expressed in both conditions. Interestingly, the majority of the meiosis-related genes were downregulated in both teratozoospermic (213 genes) and NOA (254 genes) patients compared to the control samples.
Considering only the 17 differentially expressed genes between male and female gonads in all the developmental periods, a downregulation of OOEP and TDRD1, as well as an upregulation of GTSF1 and SYCE1, was observed in MII oocytes of both advanced maternal age and young maternal age patients. A downregulation of MAJIN in GV was observed only in young maternal age, and a downregulated of BRDT and SPATA22 in advanced maternal age. In male fertility disorders, it was observed downregulation of ADAD1, BRDT, GTSF1, HORMAD1, OOEP, SPATA22, SYCP3, and TEX11 in both teratozoospermic and NOA patients in comparison to the Control groups. A downregulation of DAZL, MAEL, and SERPINA5 was demonstrated only in teratozoospermic patients, as well as a downregulation of DDX4, MAJIN, TDRD1, and TDRD12 in NOA patients, in comparison to the Control groups (Figure 2).
The logFC and the adjusted p-Values of all the statistical analysis are available in (Supplementary Table S6).
3.3 Co-expression network for OOEP gene
According to the literature review performed by the completion of this manuscript, OOEP has no recognized role in male reproduction. Therefore, considering the downregulation of OOEP observed in male reproductive disorders, a co-expression network was performed for this gene to better understand the involvement of OOEP in male reproduction. It demonstrated the co-expression of OOEP with genes related to human reproduction, spermatogenesis, male infertility, and located in the Y chromosome (Figure 3). The main functions and reproductive disorders associated with these genes are shown in (Supplementary Table S7).
FIGURE 3. Oocyte expressed protein (OOEP) gene-gene co-expression network. The edges connecting OOEP with other genes represent co-expression, and larger circles mean stronger evidence of gene co-expression. Orange means genes related with human reproduction; gray: no associated with human reproduction; green: related with spermatogenesis and/or male infertility; blue: located in the Y chromosome. DPPA5 (Developmental Pluripotency Associated 5); KHDC1 (KH Domain Containing 1); KHDC1L (KH Domain Containing 1 Like); KHDC3L (KH Domain Containing 3 Like); FAM170A (Family With Sequence Similarity 170 Member A); DDX43 (DEAD-Box Helicase 43); MAGEB6 (MAGE Family Member B6); PTPN20 (Protein Tyrosine Phosphatase Non-Receptor Type 20); ARMC4 (Armadillo Repeat-Containing Protein 4); ESX1 (ESX Homeobox 1); KDM5D (Lysine Demethylase 5D); MAEL (Maelstrom Spermatogenic Transposon Silencer); SPESP1 (Sperm Equatorial Segment Protein 1); EIF1AY (Eukaryotic Translation Initiation Factor 1A Y-Linked); DDX3Y (DEAD-Box Helicase 3 Y-Linked); RPS4Y1 (Ribosomal Protein S4 Y-Linked 1); RPS4Y2 (Ribosomal Protein S4 Y-Linked 2); TMSB4Y (Thymosin Beta 4 Y-Linked); USP9Y (Ubiquitin Specific Peptidase 9 Y-Linked); UTY (Ubiquitously Transcribed Tetratricopeptide Repeat Containing, Y-Linked).
4 Discussion
Both female and male gametogenesis begin with the mitotic proliferation of oogonia and spermatogonial stem cells, respectively, followed by meiotic divisions and modifications until the formation of mature germ cells (Larose et al., 2019). However, the timing and quantity of gametes developed are very different in the two gonads, and gene expression patterns may govern, at least in part, these differences (Handel and Schimenti, 2010). In males, the spermatogonial stem cells remain quiescent during the prenatal period, starting to proliferate and entering meiosis after puberty, resulting in continuous sperm production throughout the reproductive lifespan (Handel and Schimenti, 2010). On the other hand, female meiosis initiates in fetal ovaries and is arrested in the prophase I stage before birth, being reestablished only after puberty in periodical intervals in response to hormonal stimuli (Handel and Schimenti, 2010). To accomplish these processes, intracellular signals, such as the cAMP levels in the oocyte, are necessaire for the maintenance of meiotic arrest or for the meiotic resumption, which is characterized by the germinal vesicle (GV) breakdown, followed by completion of meiosis I, and arrest at the metaphase of meiosis II until fertilization (Pan and Li, 2019). Hence, differences in gene expression in the GV (prophase I-arrested oocytes) and MII (mature oocytes), as demonstrated in the present study, may be related to molecular mechanisms involved in the regulation of meiotic arrest and resumption; these molecular mechanisms might also be related to the acquisition of oocyte competence, pivotal to support the early embryo development (Conti and Franciosi, 2018).
Complex biological processes, such as those involved in human reproduction, require different genes working in an orchestrated way to ensure regulatory networks (de la Fuente, 2010). Recent studies of global gene expression in oocytes and sperm have demonstrated important target genes involved in conditions affecting human reproduction (Hermann et al., 2018; Li et al., 2021; Ntostis et al., 2021; Santiago et al., 2021; Sayutti et al., 2022). In the present study, prior filtering of meiosis-related genes was performed, and these genes had their differential gene expression evaluated in female and male gonads and gametes under physiological and pathological conditions.
It was demonstrated that during adulthood, 17 meiosis-related genes were upregulated in the testicle in comparison to the ovary (ADAD1, BRDT, DAZL, DDX4, GTSF1, HORMAD1, MAEL, MAJIN, OOEP, SERPINA5, SOX9, SPATA22, SYCE1, SYCP3, TDRD1, TDRD12, and TEX11), whilst in all the prenatal periods analyzed, downregulation of these genes was observed in the testicle in comparison to the ovary, except for SERPINA5 and SOX9, that maintain the upregulation in the testicle. The orchestrated expression of these genes during meiosis, according to time and place, may be crucial for meiosis entry and normal gamete formation (Figure 4), considering their involvement in important processes such as: meiotic regulation, homologous chromosomes pairing, recombination and segregation, DNA double-strand breaks (DSBs) repair, and formation of the synaptonemal complex (GeneCards, 2022). In addition, considering the 678 meiosis-related genes, it was observed that more genes were upregulated in the ovary during the prenatal period compared to the testicle, whilst in the adult gonads more genes related to meiosis were upregulated in the testicle.
FIGURE 4. Illustration of the 17 meiosis-related genes expression in human ovary and testicle according to the life period and gametogenesis timing under physiological and pathological scenarios. Female meiosis initiates in fetal ovaries and is arrested in the prophase I stage before birth, being reestablished after puberty in periodical intervals in response to hormonal stimuli. In males, the spermatogonial stem cells remain quiescent during the prenatal period, starting to proliferate and entering meiosis after puberty, resulting in continuous sperm production throughout the reproductive lifespan. It was demonstrated that during adulthood, 17 meiosis-related genes were upregulated in the testicle in comparison to the ovary (ADAD1, BRDT, DAZL, DDX4, GTSF1, HORMAD1, MAEL, MAJIN, OOEP, SERPINA5, SOX9, SPATA22, SYCE1, SYCP3, TDRD1, TDRD12, and TEX11), whilst in all the prenatal periods analyzed, downregulation of these genes was observed in the testicle in comparison to the ovary, except by SERPINA5 and SOX9, that maintain the upregulation in the testicle. It was observed a downregulation of ADAD1, BRDT, GTSF1, HORMAD1, OOEP, SPATA22, SYCP3, and TEX11 in both teratozoospermic and NOA patients in comparison to the Control groups. A downregulation of DAZL, MAEL, and SERPINA5 was demonstrated only in teratozoospermic patients, as well as a downregulation of DDX4, MAJIN, TDRD1, and TDRD12 in NOA patients, in comparison to the Control groups. Comparing GV vs. MII oocytes, a downregulation of OOEP and TDRD1, as well as an upregulation of GTSF1 and SYCE1, was observed in MII oocytes of both AMA and YMA patients. A downregulation of MAJIN in GV was observed only in YMA, and a downregulated of BRDT and SPATA22 in AMA. NOA (non-obstructive azoospermia); GV (germinal vesicle); MII (mature oocytes at metaphase II stage); YMA (young maternal age); AMA (advanced maternal age); blue up arrows mean upregulated genes, whilst red down arrows mean downregulated genes.
During the meiosis prophase I, the synaptonemal complex mediates the homologous chromosomes pairing and synapsis, a fundamental step for the recombination, that is initiated by DNA double strand breaks, and the segregation events (Handel and Schimenti, 2010). In human fetal ovaries, meiotic entry is initiated asynchronously around 8-9 wpc, confirmed by the increased expression of meiotic proteins by this time, such as DAZL, DDX4, and SYCP3 (Farini and De Felici, 2022); this increased expression was corroborated in the present study. The conditions for meiotic beginning depend on intrinsic and extrinsic factors, being proposed that germ cells of both sexes receive at the same time stimuli for meiotic beginning, however, somatic cells in the developing gonads induce or prevent meiosis entry and progression (Farini and De Felici, 2022). Therefore, the distinct pattern of meiosis-related gene expression observed in the testicle and ovary according to the life period may reflect the differences in female and male meiosis onset, through differential gene expression of germ and somatic cells in fetal and adult gonads. The upregulation of the 17 meiosis-related genes in the fetal ovary and adult testicle, coinciding with the meiosis onset in both sexes, is in agreement with the roles of these genes in processes mainly related to prophase I, the first stage of meiosis.
SERPINA5 and SOX9 were upregulated in the testicle in both prenatal and adult periods. SOX9 is a transcription factor expressed in the Sertoli cells of the male gonad throughout life (Jiménez et al., 2021). At the time of sex determination, the expression of the Y-linked, sex determining-region - SRY - in the Sertoli cells leads to the upregulation of SOX9, which activates male-promoting genes necessary for the differentiation of Sertoli cells into the testicle cords (Jiménez et al., 2021). SERPINA5 encodes a glycoprotein with an inhibitory effect in several serine proteases; in seminal plasma, inactivates several serine proteases implicated in the reproductive system, and has a control on the sperm motility and fertilization (GeneCards, 2022). Although the roles of SERPINA5 are poorly understood (Yang and Geiger, 2017), studies have pointed to a role in male reproduction (Yang and Geiger, 2017; Selvam et al., 2019). Therefore, the constant upregulation of SOX9 and SERPINA5 in the testicle samples in comparison to the ovary is in agreement with their roles in processes related to male reproduction.
Chromosome segregation errors during meiosis, leading to aneuploidies, are one of the most common causes of infertility, pregnancy failure, and genetic disorders in the offspring (Mikwar et al., 2020). Although meiotic errors may occur in both male and female germ cells, the majority are of maternal origin, being the advanced maternal age (≥35 years) a risk factor for chromosome segregation errors (Mikwar et al., 2020). The present study demonstrated 48 meiosis-related genes differentially expressed in mature oocytes of advanced maternal age in comparison to young maternal age patients. Among them are TUBB8, which encodes the primary beta-tubulin subunit of oocytes and early embryos, and BUB3, which encodes a protein involved in spindle checkpoint function. TUBB8 and BUB3 were already associated with infertility due to oocyte meiotic arrest and mosaic variegated aneuploidy syndrome, respectively (GeneCards, 2022). Genes related to meiotic recombination, chromosome cohesion, spindle assembly checkpoint, post-translational modification of histones and tubulin, and mitochondrial functions, are great targets for research, since impairment in these mechanisms has been proposed to explain the higher incidence of oocyte aneuploidies in older females (Mikwar et al., 2020).
Aside from advanced maternal age, it was expected that alterations in meiosis-related gene expression could be associated with human reproductive disorders, considering their pattern of expression in male and female gonads in critical periods of gametogenesis. Nevertheless, no differences were observed in the oocytes’ gene expression of PCOS patients in comparison to the control group. PCOS is a multifactorial disorder, influenced by intrinsic and extrinsic factors, such as insulin resistance, hyperandrogenism, environmental factors, genetics, and epigenetics (Sadeghi et al., 2022). Therefore, as many factors are involved in the etiology of PCOS, the meiosis-related genes considered in the present study probably do not play a crucial role in the mechanisms behind PCOS etiology. Furthermore, a recent study on global gene expression of oocytes and cumulus cells from PCOS patients have not demonstrated the genes evaluated here as possible candidates (Li et al., 2021).
Contrastingly, when evaluating male reproductive disorders, 238 and 290 meiosis-related genes were differentially expressed in teratozoospermic and NOA patients in comparison to the control group, respectively. Most of these genes were downregulated in the conditions affecting human male reproduction, demonstrating that expression of meiosis-related genes may be necessaire in processes related to male fertility. Additionally, 145 genes were differentially expressed in both teratozoospermic and NOA patients compared to the control group, many of them already associated with spermatogenic failure and sperm maturation defect, according to our search for phenotype descriptions in the Ensembl database. Previous studies have demonstrated a downregulation profile in meiosis-related genes in both NOA and teratozoospermic patients (Fu et al., 2016; He et al., 2021). He et al. demonstrated 430 downregulated genes in NOA in comparison to the control group, many of them associated with spermatogenesis, multicellular organism development, cell differentiation, spermatid development, and sperm motility (He et al., 2021); Fu et al. demonstrated 2158 downregulated genes in teratozoospermic in comparison to normospermic patients, mostly related to reproduction, mitotic cell cycle, and spermatogenesis (Fu et al., 2016).
Considering the 17 meiosis-related genes differentially expressed in all the developmental periods analyzed, it was observed downregulation of ADAD1, BRDT, GTSF1, HORMAD1, OOEP, SPATA22, SYCP3, and TEX11 in both teratozoospermic sperm and NOA testicle samples, evidencing a role in male reproduction. Except for OOEP, all the other downregulated genes have been linked to male infertility disorders, including NOA and teratozoospermia (Aarabi et al., 2006; Miyamoto et al., 2012; Chen et al., 2020; Greither et al., 2020; Snyder et al., 2020; Her et al., 2021; Ji et al., 2021; Wu et al., 2021; Omolaoye et al., 2022). OOEP is a member of the SCMC, a maternal complex essential to the progression beyond the first embryonic cell divisions (Li et al., 2008; Zhu et al., 2014). Although mutations in the SCMC genes are associated with female reproductive disorders (Demond et al., 2019; Zhang et al., 2019), a relationship between SCMC components and male reproduction was only recently suggested by (Feng et al., 2020), evidencing Tle6 playing a role in the proliferation and cell cycle of the mouse spermatogonia (Feng et al., 2020).
Although no obvious association of OOEP with male reproduction the co-expression network performed for OOEP demonstrated the co-expression with genes located in the Y chromosome, related to human reproduction, spermatogenesis, and male infertility. Since co-expressed genes are simultaneously active (van Dam et al., 2018), downregulation of OOEP in patients with NOA and teratozoospermia could affect the entire gene expression regulation, even though the structure or even the protein expression might not be affected. Furthermore, Ooep has been related to processes involved in DNA repair (He et al., 2018; Zhao et al., 2018), and impairment in these mechanisms is associated with spermatogenesis arrest (Gunes et al., 2015), and teratozoospermia (Ammar et al., 2020).
This study has some limitations, such as not having the validation of gene expression results by qPCR and not demonstrating protein expression. However, all data used for analysis are from experimental studies of different datasets that pointed in the same direction. In this way, the results presented here contribute to the understanding of gene expression patterns involved in the differences between male and female meiosis, by (I) demonstrating that meiosis-related genes are differentially expressed in human male and female gonads throughout gametogenesis and (II) presenting altered gene expression in human male reproductive disorders and oocytes of advanced maternal age patients. Hence, concerning basic sciences and assisted reproduction application, future studies aiming to elucidate the molecular mechanisms of meiosis should consider the genes presented in this study. These results pinpoint to possible new candidate genes that could have clinical implications, therefore should be considered biomarkers in panels for patients with NOA or teratozoospermia.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
MR contributed to the study conception, acquisition of data, data analysis, and drafted the manuscript. LF participated in the interpretation of data and manuscript revisions. TK supervised the entire study, including conception, analysis, interpretation of data, and revisions. MS contributed to the study conception, interpretation of data, and revisions to the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the Fundo de Incentivo à Pesquisa e Evento (FIPE) do Hospital de Clínicas de Porto Alegre (grants numbers 2019-0173 and 2019-0479), as well as Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES–grant number 88887.503016/2020-00).
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/fgene.2023.1125097/full#supplementary-material
References
Aarabi, M., Modarressi, M. H., Soltanghoraee, H., Behjati, R., Amirjannati, N., and Akhondi, M. M. (2006). Testicular expression of synaptonemal complex protein 3 (SYCP3) messenger ribonucleic acid in 110 patients with nonobstructive azoospermia. Fertil. Steril. 86, 325–331. doi:10.1016/J.FERTNSTERT.2005.12.070
Ammar, O., Mehdi, M., and Muratori, M. (2020). Teratozoospermia: Its association with sperm DNA defects, apoptotic alterations, and oxidative stress. Andrology 8, 1095–1106. doi:10.1111/ANDR.12778
Barrett, T., Wilhite, S., Ledoux, P., Evangelista, C., Kim, I., Tomashevsky, M., et al. (2013). NCBI GEO: Archive for functional genomics data sets--update. Nucleic Acids Res. 41, 991–995. doi:10.1093/NAR/GKS1193
Chen, S., Wang, G., Zheng, X., Ge, S., Dai, Y., Ping, P., et al. (2020). Whole-exome sequencing of a large Chinese azoospermia and severe oligospermia cohort identifies novel infertility causative variants and genes. Hum. Mol. Genet. 29, 2451–2459. doi:10.1093/HMG/DDAA101
Conti, M., and Franciosi, F. (2018). Acquisition of oocyte competence to develop as an embryo: Integrated nuclear and cytoplasmic events. Hum. Reprod. Update 24, 245–266. doi:10.1093/humupd/dmx040
Cunningham, F., Allen, J. E., Allen, J., Alvarez-Jarreta, J., Amode, M. R., Armean, I. M., et al. (2022). Ensembl 2022. Nucleic Acids Res. 50, D988–D995. doi:10.1093/NAR/GKAB1049
de la Fuente, A. (2010). From “differential expression” to “differential networking” - identification of dysfunctional regulatory networks in diseases. Trends Genet. 26, 326–333. doi:10.1016/J.TIG.2010.05.001
Demond, H., Anvar, Z., Jahromi, B. N., Sparago, A., Verma, A., Davari, M., et al. (2019). A KHDC3L mutation resulting in recurrent hydatidiform mole causes genome-wide DNA methylation loss in oocytes and persistent imprinting defects post-fertilisation. Genome Med. 11, 84–97. doi:10.1186/s13073-019-0694-y
Edgar, R., Domrachev, M., and Lash, A. (2002). Gene expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 30, 207–210. doi:10.1093/NAR/30.1.207
Farini, D., and De Felici, M. (2022). The beginning of meiosis in mammalian female germ cells: A never-ending story of intrinsic and extrinsic factors. Int. J. Mol. Sci. 23, 12571. doi:10.3390/IJMS232012571
Feng, M., Bai, Y., Chen, Y., and Wang, K. (2020). Knockout of the transducin-like enhancer of split 6 gene affects the proliferation and cell cycle process of mouse spermatogonia. Int. J. Mol. Sci. 21, 5827–5839. doi:10.3390/IJMS21165827
Fu, G., Wei, Y., Wang, X., and Yu, L. (2016). Identification of candidate causal genes and their associated pathogenic mechanisms underlying teratozoospermia based on the spermatozoa transcript profiles. Andrologia 48, 576–583. doi:10.1111/AND.12484
Gautier, L., Cope, L., Bolstad, B. M., and Irizarry, R. A. (2004). Affy - analysis of affymetrix GeneChip data at the probe level. Bioinformatics 20, 307–315. doi:10.1093/bioinformatics/btg405
GeneCards (2022). Human genes. Gene Database. Available at: https://www.genecards.org/(Accessed November 29, 2022).
Greither, T., Schumacher, J., Dejung, M., Behre, H. M., Zischler, H., Butter, F., et al. (2020). Fertility relevance probability analysis shortlists genetic markers for male fertility impairment. Cytogenet. Genome Res. 160, 506–522. doi:10.1159/000511117
Gunes, S., Al-Sadaan, M., and Agarwal, A. (2015). Spermatogenesis, DNA damage and DNA repair mechanisms in male infertility. Reprod. Biomed. Online 31, 309–319. doi:10.1016/J.RBMO.2015.06.010
Handel, M. A., and Schimenti, J. C. (2010). Genetics of mammalian meiosis: Regulation, dynamics and impact on fertility. Nat. Rev. Genet. 11, 124–136. doi:10.1038/NRG2723
He, D., Wang, L., Zhang, Z., Guo, K., Li, J., He, X., et al. (2018). Maternal gene Ooep may participate in homologous recombination-mediated DNA double-strand break repair in mouse oocytes. Zool. Res. 39, 387–395. doi:10.24272/J.ISSN.2095-8137.2018.067
He, H., Yu, F., Shen, W., Chen, K., Zhang, L., Lou, S., et al. (2021). The novel key genes of non-obstructive azoospermia affect spermatogenesis: Transcriptomic analysis based on RNA-seq and scRNA-seq data. Front. Genet. 12, 608629. doi:10.3389/FGENE.2021.608629
Her, Y. R., Wang, L., Chepelev, I., Manterola, M., Berkovits, B., Cui, K., et al. (2021). Genome-wide chromatin occupancy of BRDT and gene expression analysis suggest transcriptional partners and specific epigenetic landscapes that regulate gene expression during spermatogenesis. Mol. Reprod. Dev. 88, 141–157. doi:10.1002/MRD.23449
Hermann, B. P., Cheng, K., Singh, A., Roa-De La Cruz, L., Mutoji, K. N., Chen, I. C., et al. (2018). The mammalian spermatogenesis single-cell transcriptome, from spermatogonial stem cells to spermatids. Cell Rep. 25, 1650–1667. doi:10.1016/J.CELREP.2018.10.026
Hever, A., Roth, R., Hevezi, P., Marin, M., Acosta, J., Acosta, H., et al. (2007). Human endometriosis is associated with plasma cells and overexpression of B lymphocyte stimulator. Proc. Natl. Acad. Sci. U. S. A. 104, 12451–12456. doi:10.1073/PNAS.0703451104
Hodžić, A., Maver, A., Plaseska-Karanfilska, D., Ristanović, M., Noveski, P., Zorn, B., et al. (2021). De novo mutations in idiopathic male infertility-A pilot study. Andrology 9, 212–220. doi:10.1111/ANDR.12897
Jalili, V., Afgan, E., Gu, Q., Clements, D., Blankenberg, D., Goecks, J., et al. (2020). The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2020 update. Nucleic Acids Res. 48, W395–W402. doi:10.1093/NAR/GKAA434
Ji, Z., Yao, C., Yang, C., Huang, C., Zhao, L., Han, X., et al. (2021). Novel hemizygous mutations of TEX11 cause meiotic arrest and non-obstructive azoospermia in Chinese han population. Front. Genet. 12, 741355. doi:10.3389/FGENE.2021.741355
Jiménez, R., Burgos, M., and Barrionuevo, F. J. (2021). Sex maintenance in mammals. Genes (Basel). 12, 999. doi:10.3390/GENES12070999
Kim, D., Paggi, J., Park, C., Bennett, C., and Salzberg, S. (2019). Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37, 907–915. doi:10.1038/S41587-019-0201-4
Larose, H., Shami, A. N., Abbott, H., Manske, G., Lei, L., and Hammoud, S. S. (2019). “Gametogenesis: A journey from inception to conception,” in Current topics in developmental biology (Academic Press), 257–310. doi:10.1016/bs.ctdb.2018.12.006
Lecluze, E., Rolland, A., Filis, P., Evrard, B., Leverrier-Penna, S., Maamar, M., et al. (2020). Dynamics of the transcriptional landscape during human fetal testis and ovary development. Hum. Reprod. 35, 1099–1119. doi:10.1093/HUMREP/DEAA041
Li, J., Chen, H., Gou, M., Tian, C., Wang, H., Song, X., et al. (2021). Molecular features of polycystic ovary syndrome revealed by transcriptome analysis of oocytes and cumulus cells. Front. Cell Dev. Biol. 9, 2378. doi:10.3389/fcell.2021.735684
Li, L., Baibakov, B., and Dean, J. (2008). A subcortical maternal complex essential for preimplantation mouse embryogenesis. Dev. Cell 15, 416–425. doi:10.1016/j.devcel.2008.07.010
Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi:10.1093/BIOINFORMATICS/BTT656
Mascarenhas, M. N., Flaxman, S. R., Boerma, T., Vanderpoel, S., and Stevens, G. A. (2012). National, regional, and global trends in infertility prevalence since 1990: A systematic analysis of 277 health surveys. PLoS Med. 9, e1001356. doi:10.1371/journal.pmed.1001356
Mikwar, M., MacFarlane, A. J., and Marchetti, F. (2020). Mechanisms of oocyte aneuploidy associated with advanced maternal age. Mutat. Res. Rev. Mutat. Res. 785, 108320. doi:10.1016/J.MRREV.2020.108320
Miyamoto, T., Tsujimura, A., Miyagawa, Y., Koh, E., Namiki, M., Horikawa, M., et al. (2012). Single-nucleotide polymorphisms in HORMAD1 may be a risk factor for azoospermia caused by meiotic arrest in Japanese patients. Asian J. Androl. 14, 580–583. doi:10.1038/AJA.2011.180
Montojo, J., Zuberi, K., Rodriguez, H., Kazi, F., Wright, G., Donaldson, S., et al. (2010). GeneMANIA Cytoscape plugin: Fast gene function predictions on the desktop. Bioinformatics 26, 2927–2928. doi:10.1093/BIOINFORMATICS/BTQ562
Ntostis, P., Iles, D., Kokkali, G., Vaxevanoglou, T., Kanavakis, E., Pantou, A., et al. (2021). The impact of maternal age on gene expression during the GV to MII transition in euploid human oocytes. Hum. Reprod. 37, 80–92. doi:10.1093/HUMREP/DEAB226
Ombelet, W., Cooke, I., Dyer, S., Serour, G., and Devroey, P. (2008). Infertility and the provision of infertility medical services in developing countries. Hum. Reprod. Update 14, 605–621. doi:10.1093/humupd/dmn042
Omolaoye, T. S., Hachim, M. Y., and du Plessis, S. S. (2022). Using publicly available transcriptomic data to identify mechanistic and diagnostic biomarkers in azoospermia and overall male infertility. Sci. Rep. 12, 2584. doi:10.1038/S41598-022-06476-1
Pan, B., and Li, J. (2019). The art of oocyte meiotic arrest regulation. Reprod. Biol. Endocrinol. 17, 8–12. doi:10.1186/S12958-018-0445-8
Platts, A., Dix, D., Chemes, H., Thompson, K., Goodrich, R., Rockett, J., et al. (2007). Success and failure in human spermatogenesis as revealed by teratozoospermic RNAs. Hum. Mol. Genet. 16, 763–773. doi:10.1093/HMG/DDM012
Ritchie, M., Phipson, B., Wu, D., Hu, Y., Law, C., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 43, 47–59. doi:10.1093/NAR/GKV007
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, 139–140. doi:10.1093/BIOINFORMATICS/BTP616
Sadeghi, H. M., Adeli, I., Calina, D., Docea, A. O., Mousavi, T., Daniali, M., et al. (2022). Polycystic ovary syndrome: A comprehensive review of pathogenesis, management, and drug repurposing. Int. J. Mol. Sci. 23, 583. doi:10.3390/IJMS23020583
Santiago, J., Silva, J. V., Howl, J., Santos, M. A. S., and Fardilha, M. (2021). All you need to know about sperm RNAs. Hum. Reprod. Update 28, 67–91. doi:10.1093/HUMUPD/DMAB034
Sayutti, N., Abu, M. A., and Ahmad, M. F. (2022). PCOS and role of cumulus gene expression in assessing oocytes quality. Front. Endocrinol. (Lausanne). 13, 768. doi:10.3389/fendo.2022.843867
Selvam, M. K. P., Agarwal, A., Pushparaj, P. N., Baskaran, S., and Bendou, H. (2019). Sperm proteome analysis and identification of fertility-associated biomarkers in unexplained male infertility. Genes (Basel). 10, 522. doi:10.3390/GENES10070522
Snyder, E., Chukrallah, L., Seltzer, K., Goodwin, L., and Braun, R. E. (2020). ADAD1 and ADAD2, testis-specific adenosine deaminase domain-containing proteins, are required for male fertility. Sci. Rep. 10, 11536. doi:10.1038/S41598-020-67834-5
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinforma. 54, 1.30.1–1.30.33. doi:10.1002/CPBI.5
Thonneau, P., Marchand, S., Tallec, A., Ferial, M. L., Ducot, B., Lansac, J., et al. (1991). Incidence and main causes of infertility in a resident population (1,850,000) of three French regions (1988-1989). Hum. Reprod. 6, 811–816. doi:10.1093/OXFORDJOURNALS.HUMREP.A137433
Tweedie, S., Braschi, B., Gray, K., Jones, T. E. M., Seal, R. L., Yates, B., et al. (2021). Genenames.org: The HGNC and VGNC resources in 2021. Nucleic Acids Res. 49, D939–D946. doi:10.1093/NAR/GKAA980
van Dam, S., Võsa, U., van der Graaf, A., Franke, L., and de Magalhães, J. (2018). Gene co-expression analysis for functional classification and gene-disease predictions. Brief. Bioinform. 19, 575–592. doi:10.1093/BIB/BBW139
Wood, J. R., Dumesic, D. A., Abbott, D. H., and Strauss, J. F. (2007). Molecular abnormalities in oocytes from women with polycystic ovary syndrome revealed by microarray analysis. J. Clin. Endocrinol. Metab. 92, 705–713. doi:10.1210/JC.2006-2123
Wu, Y., Li, Y., Murtaza, G., Zhou, J., Jiao, Y., Gong, C., et al. (2021). Whole-exome sequencing of consanguineous families with infertile men and women identifies homologous mutations in SPATA22 and MEIOB. Hum. Reprod. 36, 2793–2804. doi:10.1093/HUMREP/DEAB185
Yang, H., and Geiger, M. (2017). Cell penetrating SERPINA5 (ProteinC inhibitor, PCI): More questions than answers. Semin. Cell Dev. Biol. 62, 187–193. doi:10.1016/J.SEMCDB.2016.10.007
Zhang, W., Chen, Z., Zhang, D., Zhao, B., Liu, L., Xie, Z., et al. (2019). KHDC3L mutation causes recurrent pregnancy loss by inducing genomic instability of human early embryonic cells. PLoS Biol. 17, 3000468–3000498. doi:10.1371/JOURNAL.PBIO.3000468
Zhao, B., Zhang, W., Cun, Y., Li, J., Liu, Y., Gao, J., et al. (2018). Mouse embryonic stem cells have increased capacity for replication fork restart driven by the specific Filia-Floped protein complex. Cell Res. 28, 69–89. doi:10.1038/CR.2017.139
Keywords: oogenesis, spermatogenesis, germinal vesicle, mature oocyte, infertility
Citation: Rockenbach MK, Fraga LR, Kowalski TW and Sanseverino MTV (2023) Expression profiles of meiotic genes in male vs. female gonads and gametes: Insights into fertility issues. Front. Genet. 14:1125097. doi: 10.3389/fgene.2023.1125097
Received: 15 December 2022; Accepted: 02 March 2023;
Published: 14 March 2023.
Edited by:
Kenta Nakai, The University of Tokyo, JapanReviewed by:
Nélida Rodríguez-Osorio, Universidad de la República, UruguayHong Zhang, University of Basel, Switzerland
Copyright © 2023 Rockenbach, Fraga, Kowalski and Sanseverino. 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: Thayne Woycinck Kowalski, dGhheW5ld2tAZ21haWwuY29t, dGtvd2Fsc2tpQGhjcGEuZWR1LmJy