Reut Shabtai
Yonatan B. Tzur- Department of Genetics, Institute of Life Sciences, The Hebrew University of Jerusalem, Jerusalem, Israel
The testis is the mammalian tissue with the highest expression levels of long intergenic non-coding RNAs (lincRNAs). However, most in vivo models have not found significant reductions in male fertility when highly expressed lincRNA genes were removed. This suggests that certain lincRNAs may act redundantly or lack functional roles. In the genome of the nematode Caenorhabditis elegans, there is an order of magnitude fewer lincRNA genes than in mammals. This characteristic lowers the potential for redundancy, making it an ideal model to test these possibilities. We identified five highly and dynamically expressed lincRNAs in male C. elegans gonads and quantified the fertility of worm strains in which these genes were removed. In contrast to the hermaphrodites of deletion strains, which exhibited no significant reductions in broods, smaller brood sizes were observed in the progeny of males of three of the lincRNA deleted strains. This demonstrates reduced male fertility in worms with those genes removed. Interestingly, reduced brood size was statistically significant only in the last days of egg laying in two of these strains. This suggests the effect is due to early deterioration and aging of the transferred sperm. We detected a mild increase in embryonic lethality in only one of the strains, supporting the possibility that these lincRNAs do not affect fertility through critical roles in essential meiotic processes. Together our results indicate a sexually dimorphic outcome on fertility when lincRNA are removed and show that, unlike mammals, individual lincRNAs in C. elegans do play significant roles in male fertility.
Introduction
Spermatogenesis is a complex developmental plan in which germ stem cells differentiate into mature spermatozoa. Spermatogenesis consists of several stages including primordial germ cell expansion, two meiotic divisions, and differentiation into spermatozoa (L'Hernault, 1997; Yoshida, 2010; Chu and Shakes, 2013; Ellis and Stanfield, 2014; Griswold, 2016). During the meiotic stage, the number of chromosomes is reduced by half to create haploid cells. This is achieved by unique chromosome interactions, including homologous chromosome pairing, synapsis, and recombination (Colaiacovo et al., 2003; Couteau and Zetka, 2005; Hayashi et al., 2010; Schild-Prufert et al., 2011; Reichman et al., 2017). Unlike many other stages in spermatogenesis, meiosis is also executed in a similar fashion in oogenesis. Failure to complete meiotic specific processes, in many cases, leads to apoptotic programmed cell death (Ye et al., 2014).
Long non-coding RNAs (lncRNAs) are transcribed in a similar process as mRNAs, and are often capped, spliced, and poly-adenylated, but not translated. LncRNAs are transcribed from tens of thousands of loci in the human genome (reviewed in (Ulitsky and Bartel, 2013; Fatica and Bozzoni, 2014; Marques and Ponting, 2014; Blythe et al., 2016; Melissari and Grote, 2016; Deniz and Erman, 2017; Kopp and Mendell, 2018; Fico et al., 2019; Shields et al., 2019). Approximately half of the human lncRNAs are transcribed from genomic loci that don’t overlap with coding genes and are denoted as long intergenic non-coding RNAs (lincRNAs) (Uszczynska-Ratajczak et al., 2018).
Testis has the most complex transcriptome and expresses the highest levels and the largest repertoire of lncRNAs compared with all other mammalian tissues (Soumillon et al., 2013; Necsulea et al., 2014; Washietl et al., 2014; Hezroni et al., 2015; Hong et al., 2018). Surprisingly, in almost all cases, lncRNA in vivo knockout or knockdown mammalian models failed to exhibit significant reductions in male fertility (e.g., (Mehta et al., 2021; Chadourne et al., 2021; Li et al., 2021; Zhang et al., 2021). Several hypotheses were suggested to explain this conundrum, including promiscuous transcription due to rapid changes in chromatin structure and functional redundancy of several lncRNAs (Reviewed in (Tzur, 2022).
Wild-type C. elegans nematodes exist as hermaphrodites and males. Hermaphrodites, which have two X chromosomes, produce sperm during larval stages, but switch to oogenesis in adulthood (L'Hernault, 1997; Chu and Shakes, 2013; Ellis and Stanfield, 2014; Griswold, 2016). Oocytes can be fertilized by self-sperm or by sperm transferred by males. Several genetic mutations lead to hermaphrodites that cannot self-fertilize, making them functionally females. Male C. elegans worms have only one X chromosome, produce only sperm, and are present at ∼0.1% of the population under normal laboratory conditions. In both male and hermaphrodite gonads, the germ cells are arranged in a spatio-temporal manner from proliferative stem cells to mature gametes (L'Hernault, 1997; Chu and Shakes, 2013; Ellis and Stanfield, 2014; Griswold, 2016). Humans and C. elegans have a similar number of coding genes, yet only a few hundred lincRNA genes were found in the latter’s genome (Nam and Bartel, 2012; Akay et al., 2019). This order of magnitude fewer lincRNA genes reduces the likelihood of redundancy in the lincRNA roles present during the worm’s spermatogenesis process.
The study of fertility and meiosis in C. elegans has uncovered many evolutionary conserved processes (reviewed in (Hubbard and Greenstein, 2005; Lui and Colaiacovo, 2013; Hillers et al., 2017; Hubbard and Schedl, 2019). However, most of this work has examined hermaphrodites. Far less is known about fertility in males [reviewed in (LHernault, 2009; Ellis and Stanfield, 2014)]. Similar to hermaphrodites, the distal side of the male gonad cells undergo proliferation and complete meiotic reductional division as they move proximally. This step includes pairing, synapsis and crossovers of homologous chromosomes (Chu and Shakes, 2013). Similar to other metazoans, in C. elegans males the meiotic divisions create four gametes (in contrast to oogenesis in which only one oocyte is formed from each progenitor cell), as well as shedding most of the cytoplasmic components and formation of a residual body (Shakes et al., 2009; Chu and Shakes, 2013) and chromatin compaction assisted by protamine proteins as spermatids develop into mature spermatozoa (Nishimura and LHernault, 2017).
We have previously analyzed the transcriptomic changes along the stages of oogenesis and spermatogenesis in C. elegans. We used laser capture microdissection to cut both hermaphrodite and male gonads into 10 sequential segments (Tzur et al., 2018). RNASeq analysis allowed us to quantitatively compare gene expression between the two gametogenesis processes and between different stages within oogenesis or spermatogenesis. To assess the roles of lincRNAs in oogenesis, we previously used these databases and found lincRNAs that are highly and dynamically expressed in the hermaphrodite gonad (Ishtayeh et al., 2021). We engineered full genomic homozygous deletion strains for these lincRNA genes, thus preventing expression of any part of the gene ((Ishtayeh et al., 2021) and Methods). Surprisingly, we found no change in hermaphrodite fertility without those lincRNAs (Ishtayeh et al., 2021). These worms also did not have higher than wild-type rates of embryonic lethality, germline apoptosis, defects in synapsis of homologous chromosomes, or bivalent structure of mature oocytes. Therefore, these lincRNAs are not required for normal spermatogenesis and oogenesis in hermaphrodites (Ishtayeh et al., 2021). However, the effect of these lincRNA gene deletions on male fertility has not been explored.
In this article we report our analysis of male fertility in worm strains in which we deleted the five lincRNA genes which we previously determined are redundant for hermaphrodite fertility (Ishtayeh et al., 2021). In three of those strains, we found a significant reduction in the brood size of mutant males, and in one we found a mild increase in embryonic lethality. This stands in contrast with hermaphrodites, despite the fact that they also produce sperm. Thus, our work indicates that some of the lincRNAs in worms are required specifically for male fertility.
Materials and methods
Strains and alleles
The fem-2 worms were cultured at 15°C and transferred to 25°C prior to progeny quantification experiments (see below). All other strains were cultured under standard conditions at 20°C (Brenner, 1974). The N2 Bristol strain was utilized as the wild-type background. Worms were grown on NGM plates with Escherichia coli OP50 (Brenner, 1974).
The following mutations were used in this study: LGI: linc-9(huj24), LGII: linc-4(huj25), LGIII: fem-2(b245), LGV: linc-168(huj10), LGX: linc-7(huj9), linc-20(huj21).
Engineering the lincRNA deletion strains was reported in (Ishtayeh et al., 2021). In short, we used the CRISPR genome editing method described in (Achache et al., 2019) with gRNAs directed to regions upstream and downstream of the gene. This strategy ensures that every single base of the gene is removed from the genome. Deletions were identified by PCR and complete deletion of the genes was verified by Sanger sequencing. Sequences and molecular data are detailed in Supplementary Table S1 in (Ishtayeh et al., 2021). After five times outcrossing, to minimize possible off-target mutations, we established homozygous strains which were used in this study.
LincRNA selection
To find lincRNAs which affect male but not hermaphrodite fertility we used lincRNAs reported in (Ishtayeh et al., 2021), in which full lincRNA deletion alleles were engineered and five times outcrossed strains were carefully tested for hermaphrodite’s fertility related phenotypes. The strains that showed no effect were selected and the level of the transcripts along the male gonad were evaluated from the datasets published in (Tzur et al., 2018). LincRNAs with dynamic expression and at least an average of six counts along the ten segments were selected for analysis in males.
RNA level analysis
RNASeq data of the lincRNAs and meiotic genes presented in Supplementary Figure S1 were extracted from the male analyses reported in (Tzur et al., 2018). Normalized values of the ten sequential segments from two male gonads were averaged, increased by one and log2 transformed. The resulting values were plotted as data points of the values vs. the 10 segments from (1) proliferating to the (10) spermatozoa stage.
Progeny quantification
To quantify the progeny, we used the well based method of Plate Phenotype (Fay, 2013). Worms were feminized by transferring fem-2(b245) L4 larvae from 15°C to 25°C. 48–72 h after the temperature shift, L4 fem-2(b245) F1 progeny were individually isolated to plates with five young adult males of the appropriate genotype (wild type or with specific lincRNA deletion) and maintained at 20°C until the end of the experiment. After 24 h, the female and male worms were transferred to a new plate. The female was transferred four additional times every 24 h.
The embryos and hatched progeny in each plate were counted at the time of the adult transfer and 24 h later. Only data from worms that survived until completion of the experiment were used. The n value was: linc-4 = 19, linc-7 = 17, linc-9 = 16, linc-20 = 11, linc-168 = 16, and wild type = 64 individual mated females. To verify full feminization, at least three plates with individual fem-2(b245) worms were isolated without males and monitored for zero F1 progeny. Mating feminized worms with wild-type males served as control for the effects of males with lincRNA deletion.
Embryonic lethality quantification
Embryonic lethality was quantified as in (Achache et al., 2019). In short, the number of larvae and embryos were counted in each plate and used for progeny quantification (see above) at two time points: After the transfer of the P0 individuals, and 24 h later. C. elegans L1 larvae hatch ∼800 min after fertilization (Sulston et al., 1983), therefore any embryo found on the plate during the second count (24 h after removal of the P0 worms) was termed “dead”. For each mated feminized worm, the total number of dead embryos was divided by the total progeny (larvae and dead embryos) to determine this worm’s progeny embryo lethality.
Statistics
For the data presented in Figure 4B we used the Fisher’s exact test. For all other analyses we used the two tailed Mann-Whitney U test.
Results
Analysis of lincRNA gene expression
Germ cells in the gonad of both male and hermaphrodite worms are arranged in spatiotemporal pattern from proliferative cells to mature gametes (L'Hernault, 1997; Chu and Shakes, 2013; Ellis and Stanfield, 2014; Griswold, 2016). We have utilized this setup in the past to determine the transcriptome of 10 sequential segments in along gonads of hermaphrodites and males to analyze gene expression during oogenesis and spermatogenesis respectively (Tzur et al., 2018). In that work, RNASeq analysis was used to quantify the transcriptome in 10 sequential segments of gonads from both sexes (Tzur et al., 2018). In a follow-up work we used the datasets to identify lincRNAs genes that are highly and dynamically expressed in the hermaphrodite gonads (Ishtayeh et al., 2021). Using CRISPR genome engineering we created homozygote worm strains with a full deletion of one of these lincRNAs. Hermaphrodite worms of these strains showed no change in fertility under normal laboratory conditions (Ishtayeh et al., 2021). The lack of fertility reduction in these strains indicates the lincRNAs are redundant for both oogenesis and spermatogenesis in hermaphrodites. Previous reports identified genes that play a role in male, but not hermaphrodite, spermatogenesis [e.g., (Stanfield and Villeneuve, 2006)]. To find lincRNA genes with fertility roles only in males we focused on lincRNAs found to be redundant for hermaphrodite fertility (Ishtayeh et al., 2021) and have significant and dynamic expression in the male gonad (see Methods).
We used the RNAseq datasets of male gonad sections published in (Tzur et al., 2018) to find the expression patterns of these lincRNAs in males. We found that five lincRNAs (linc-4, linc-7, linc-9, linc-20, and linc168) which were found to be redundant for hermaphrodite fertility (Ishtayeh et al., 2021) are highly and dynamically expressed in males (Figure 1). Two of these lincRNAs, linc-7 and linc-168, are among the highest expressed transcripts in the male gonad. The peak expression of linc-7 and linc-168 is higher than many genes known to be essential for successful spermatogenesis (e.g., mpk-1, syp-1, and htp-3; Supplementary Figure S1). Only the male-specific protein gene msp-51 has a peak expression higher than linc-168. In contrast, most lincRNAs are generally expressed at lower levels than coding genes in both mammals and C. elegans (Nam and Bartel, 2012; Soumillon et al., 2013; Akay et al., 2019).
FIGURE 1. Expression patterns of lincRNAs in the male gonad. Log2 of normalized expression values of five lincRNAs with high levels of expression along the male gonad from proliferative to mature sperm stages. X-axis numbers correspond to the segments used for the analysis and refer to the following stages: 1-2 proliferative, 2–4 leptotene/zygotene, 5-6 pachytene, 7-8 condensation, and division, 9-10 spermiogenesis. Adapted from (Tzur et al., 2018).
Our analysis indicates several similarities between oogenesis and spermatogenesis lincRNA expression. Ilinc-7 and linc-168 are the lincRNAs with the highest expression in both developmental processes [Figure 1 and (Ishtayeh et al., 2021)]. Additionally, linc-7 expression is higher at the first half of the gonad in both sexes, while linc-168 expression is mostly stable [Figure 1 and (Ishtayeh et al., 2021)]. linc-7 and linc-20, which are both transcribed from the X chromosome, are expressed during the early parts of the gonads [Figure 1 and (Ishtayeh et al., 2021)]. This stands in contrast to most X-linked genes which are silenced during this stage in both males and hermaphrodites (Kelly et al., 2002; Bean et al., 2004; Tzur et al., 2018). Some differences do exist in the expression pattern of these lincRNAs between spermatogenesis and oogenesis. For example, linc-9 is expressed at the highest levels during late pachytene and diplotene during oogenesis (Ishtayeh et al., 2021), but is primarily, present in premeiotic and leptotene/zygotene stages in spermatogenesis (Figure 1). We conclude that in male gonads, expression patterns of the lincRNAs discussed here are mostly high and share similar dynamics to their expression in hermaphrodite gonads.
Deletion of linc-7, linc-9, and linc-168 leads to reduced male fertility
To test male fertility without the interference of hermaphrodite self-fertilization we used fem-2(b245) worms which contain a temperature-sensitive mutation that causes XX worms to be functionally females at 25°C [see Methods (Hodgkin, 1986)]. We isolated L4 feminized larvae together with males for 2 days, removed the males, and continued to score the progeny for four more days (see Methods). Fertilized females mated with wild-type males laid an average of 543 ± 17 embryos during this 6 days period, in line with previous reports of mated wild-type females (Pickett et al., 2013). In contrast, females mated with male worms with lincRNA deletions showed a significant reduction in the number of laid embryos in three of the five strains tested (Figure 2; linc-7, linc-9, and linc-168, p-value by the Mann-Whitney test = 0.04, 0.02, and 0.008 respectively). We observed a similar, but not statistically significant, reduction for linc-20 males (Figure 2). Interestingly, linc-9, and linc-20 are highly paralogous (Ishtayeh et al., 2021), suggesting they may be partially redundant [albeit see (Rappaport et al., 2022)]. We note that linc-4 males had brood size larger than those of the wild-type males, however, this was not statistically significant (p-value >0.1 by the Mann-Whitney test). In conclusion three out of the five lincRNA deletion strains have a male-specific reduction in fertility.
FIGURE 2. Deletion of three lincRNA genes leads to reduced male fertility. Average progeny brood sizes for females mated with males of the indicated lines. Mann-Whitney p-value: * <0.05, ** <0.01.
No substantial increase in embryonic lethality of the progeny of lincRNAs deletion mutants
Several mutations that create aberrations in gametogenesis, especially those connected with meiotic processes, were shown to lead to embryonic lethality of the progeny. For example, mutations that hamper homologous chromosome pairing (Phillips and Dernburg, 2006), synapsis (MacQueen et al., 2002; Smolikov et al., 2007a; Smolikov et al., 2007b) and recombination, prevent correct transfer of chromosomes to the egg and sperm. This prevents successful completion of embryogenesis of the progeny. During oogenesis, many meiocytes that fail to correctly undergo synapsis and recombination are removed by apoptosis and do not end up as oocytes. In the worm’s male germline, these failures don’t lead to apoptosis, meaning even mild defects can increase embryonic lethality (Bhalla and Dernburg, 2005; Jaramillo-Lambert et al., 2010; Ye et al., 2014; Bohr et al., 2016). To find if the reduced male fertility we found in three of the lincRNA deletions (Figure 2) also lead to increased embryonic lethality we quantified the number of hatched and unhatched embryos. In the progeny of the lincRNA mutant males, we did not observe any substantial increase in embryonic lethality (Figure 3). Notably linc-9 and linc-20 did show a slight increase in embryonic lethality (0.4% vs. 0.1% compared to wild-type males), but only linc-9 showed a statistically significant increase. These results raise the possibility that these two extremely paralogous lincRNAs may play redundant roles. Taken together these mutations do not lead to any relevant increase in embryonic lethality, suggesting the reduction in fertility isn’t caused by meiotic failures.
FIGURE 3. Progeny of males with deletions in lincRNA genes do not undergo substantial embryonic lethality. Average embryonic lethality of the progeny of males of the indicated genotypes. Mann-Whitney p-value: * <0.05.
Most of the reduction in fertility in females mated with linc-7 and linc-168 males occurs during late stages of egg laying.
The egg-laying dynamics in C. elegans change under different conditions. For example, hermaphrodite worms reach their egg-laying peak during their second day of adulthood while fertilized females reach it on the third day [Figure 4 and (Pickett et al., 2013)]. Moreover, several mutations were shown to change the dynamics of egg-laying [e.g., (Kadandale and Singson, 2004)].
FIGURE 4. Dynamics of brood size along the reproductive term of lincRNA deleted males. (A) Average progeny brood size for females mated with males of the indicated genotypes along six 24 h intervals. X-axis timepoints refer to days post fertilization. (B) Percentage of embryos laid on day 2 and day 3 by progeny of WT and linc-168 males.
To find whether the deletion of the lincRNAs resulted in a change in the egg-laying dynamics we compared the brood size in each of the 24 h intervals over the 6 days period (Figure 4A). We found that the general dynamics of egg-laying of females mated with the lincRNA mutant males were similar to wild-type males, and most showed a peak of brood production in the third day of adulthood (Figure 4A). The only exception was linc-168 in which 56% of the mated worms reached a peak on the second day of adulthood compared to 24% of the worms mated with wild-type males (Figure 4B; p-value <0.05 by the Fisher exact test). In the three strains that demonstrated reduced fertility we found variation in the time points in which most of the reduction compared to wild type was measured. linc-9 displayed a non-significant increase in brood size the first day after mating but had significant reduction in the number of progeny compared to wild-type in days that followed (Figure 4A; p-value <0.05 by the Mann-Whitney test). Conversely, in linc-7 and linc-168, we found reduced brood size across the entire tested period, but we only found a statistically significant reduction in broods on days five and six after mating. Given that the total brood size of linc-9 is larger than linc-7 and linc-168 males, these results raise the possibility for a more substantial effect in the last days of egg laying after fertilization by linc-7 and linc-168. We conclude that the reduction in brood size of females mated by lincRNA deletion strain males occurs at different times along the fertility term depending on the removed gene.
Discussion
How can the same mutation lead to male-specific fertility phenotype in a process that occurs in both sexes?
Sexual dimorphism manifests as different phenotypes between males and females with a similar genetic background. In C. elegans it is possible to test this phenomenon for spermatogenesis since both sexes produce sperm. Stanfield and Villeneuve reported in the past that a mutation in swm-1 reduces sperm activation, and therefore fertility in male worms but not in hermaphrodites (Stanfield and Villeneuve, 2006). Conversely, mutations in genes from the spe-8 group are required for hermaphrodite self-fertility, but males are fully fertile (LHernault et al., 1988; Shakes and Ward, 1989). Mutant lincRNA genes tested in this work displayed effects on male fertility but did not change reproduction of hermaphrodites. Despite the fact that fertility in both sexes depends on sperm, surprisingly, we found a significant reduction in male fertility in three out of the five mutants tested.
When comparing the effects of these lincRNAs in hermaphrodites [as reported in (Ishtayeh et al., 2021)] to their effects in males we report here, several similarities exist. In both sexes two of these lincRNAs (linc-4 and linc-20) are redundant for fertility. Also, deletion of none of these lincRNA led to high embryonic lethality. On the other hand, here we show that three lincRNAs are required for normal fertility in males but not in hermaphrodites, although both sexes produce sperm.
What could lead to the difference between male and hermaphrodite when it comes to lack of lincRNA genes? Several hypotheses can be envisioned. Spermatozoa differs in males and hermaphrodites. For example, mature sperm in males are larger and faster than those produced by hermaphrodites (LaMunyon and Ward, 1998). Moreover, unlike hermaphrodites, male sperm must crawl from the vulva to the spermatheca. Therefore, subtle changes in crawling efficiency or orientation will be more critical for successful fertilization in male vs. hermaphrodite sperm. It is thus possible that the lincRNA mutations affect specific features in male sperm such as size or activation. Second, although the basic genetic program of spermatogenesis is similar between males and hermaphrodites, some variations in gene expression are present (L'Hernault, 1997; Ellis and Stanfield, 2014). Moreover, hermaphrodites have two X chromosomes whereas males have only one. Of note, two of the lincRNAs, linc-7, and linc-20, are coded on the X chromosome, and contrary to most X-linked genes they are expressed in male gonads and early stages of hermaphrodite oogenesis [Figure 1 and (Tzur et al., 2018; Ishtayeh et al., 2021)]. lincRNA deletion can affect gene expression within the worm’s germline [e.g., (Ishtayeh et al., 2021)], potentially leading to effects on male-specific pathways. Third, fertility could be reduced due to effects on the soma and not the germline. This could be manifested by specific sex-related roles of the soma, such as the requirement for males to interreact and transfer the sperm to the hermaphrodite, or specific interactions between somatic and germ cells within the same organism. Deletion of the lincRNA could result in a male-specific somatic change that will affect sperm quality and or transfer efficiency.
The reduced fertility of linc-7 and linc-168 is probably due to reduced quality or early aging of the transferred sperm.
Given these options, what could be the cause of the reduced fertility in the specific strains described here? Failure in chromosome segregation, whether it be meiotic or mitotic, leads to embryonic lethality. In hermaphrodites, but not in males, these are also accompanied by apoptosis. Therefore, meiotic failures are expected to lead to higher rates of embryonic lethality of the progeny if they occur in males. The only strain in which statistically significant embryonic lethality of the progeny is observed is linc-9. Yet, even in that case, the increase was extremely mild, below 1%, far less than the reduction in fertility. This makes the possibility of meiotic critical errors less likely to be the cause of the reduced fertility in this strain.
The analysis of brood size dynamics in the days following the mating is more informative. In our experiments the feminized worms were mated as L4 larvae for 48 h. The sperm which was transferred to them by the males had to stay functional for the entire period of the experiment (∼6 days). Consistent lower progeny (as compared to mating with wild type males) throughout the experiment period would suggest an overall lower sperm quality, whereas reduction only in the last days of the experiment would suggest that the sperm that was transferred deteriorated within the mated females. The significant reduction in brood from the second to the sixth day after mating in linc-9, correspond to the first scenario and suggests that it is lower quality sperm that leads to the reduced brood sizes. This lower quality could result in a reduced number of successful fertilization events already on the second day. On the other hand, the specific reduction in brood size of females mated with linc-7 and linc-168 during days five and six can correspond to the second option, and therefore the sperm that was transferred “aged” and deteriorated after a few days. It is also possible that the number of sperm transferred was lower. However, the fact that we introduced an excess number of males per female (see methods) suggests the former option is more likely. These results raise the question regarding the molecular mechanism by which linc-7, linc-9, and linc-168 promote normal male fertility. Very little is known about the molecular roles of lincRNAs in C. elegans. We showed that deletion of linc-4 leads to transcriptomic changes, with a significant enrichment of germline genes and genes involved in cuticle formation (Ishtayeh et al., 2021). It is therefore possible that the three lincRNAs that have roles in male fertility also work via transcription modulation. Additionally, the very high levels of linc-168 and linc-7 raise the possibility they act as competitive inhibitors for RNA, proteins and specific genomic sites, as was previously demonstrated for lincRNAs in mammals (Ulitsky and Bartel, 2013; Fatica and Bozzoni, 2014; Marques and Ponting, 2014; Blythe et al., 2016; Melissari and Grote, 2016; Deniz and Erman, 2017; Kopp and Mendell, 2018; Fico et al., 2019; Shields et al., 2019).
Is it possible that in invertebrates there is less lncRNA redundancy than in mammals?
LincRNAs have generally low sequence evolutionary conservation in both mammals and nematodes (Ulitsky and Bartel, 2013; Fatica and Bozzoni, 2014; Marques and Ponting, 2014; Blythe et al., 2016; Melissari and Grote, 2016; Deniz and Erman, 2017; Kopp and Mendell, 2018; Fico et al., 2019; Shields et al., 2019). Previous works suggested that lincRNAs with more prominent physiological roles will also show higher sequence conservation. However, linc-7, linc-9 and linc-168 have very low conservation even when compared to C. brenneri, C. briggsae, and C. remanei (Blast score lower than 45).
Our results stand in contrast to mice models in which many knock-outs of a single lincRNAs didn’t significantly reduce male fertility. We propose that similar to Drosophila, in which single mutations in many lncRNAs genes reduced male fertility (Wen et al., 2016), C. elegans highly expressed lincRNAs aren’t redundant. However, our results did raise the possibility that some level of redundancy may be present between linc-9 and linc-20. It will be interesting to revisit the roles of lncRNAs in mammalian spermatogenesis and find groups of lncRNAs that perform the same molecular role. Taken together we conclude that the reduced fertility in males could be the result of either the number, quality, or longevity of the sperm. Future studies will probably test these models and find how the lincRNAs affect mating, the number, and the quality of the transferred sperm.
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
RS performed experiments and analysed the results. YBT designed experiments, analysed the results, and wrote the manuscript.
Funding
This work was supported by Israel Science Foundation (Grant Number 979/21), Ministry of Science and Technology, Israel (Grant Number 100594), and the Shemesh fund to YBT.
Acknowledgments
We thank the Caenorhabditis Genetics Center for kindly providing strains. We thank the Tzur lab and especially Kayla Nennig-Kniaz for helpful discussions, proofreading and valuable comments.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2023.1115605/full#supplementary-material
Supplementary Figure S1 | Comparison between the expression patterns of linc-7, linc-168, and prominent meiotic genes. Log2 of normalized expression values of linc-7 and linc-168 and the indicated meiotic genes along the male gonad from proliferative (1) to mature sperm (10) stages. X-axis numbers correspond to the segments used for the analysis and refer to the following stages: 1-2 proliferative, 2-4 leptotene/zygotene, 5-6 pachytene, 7-8 condensation and division, 9-10 spermiogenesis. Note the maximum levels of the lincRNAs are higher than those of the meiotic genes. Adapted from (Tzur et al., 2018).
References
Achache, H., Laurent, L., Hecker-Mimoun, Y., Ishtayeh, H., Rappaport, Y., Kroizer, E., et al. (2019). Progression of meiosis is coordinated by the level and location of MAPK activation via OGR-2 in Caenorhabditis elegans. Genetics 212, 213–229. doi:10.1534/genetics.119.302080
Akay, A., Jordan, D., Navarro, I. C., Wrzesinski, T., Ponting, C. P., Miska, E. A., et al. (2019). Identification of functional long non-coding RNAs in C. elegans. BMC Biol. 17, 14. doi:10.1186/s12915-019-0635-7
Bean, C. J., Schaner, C. E., and Kelly, W. G. (2004). Meiotic pairing and imprinted X chromatin assembly in Caenorhabditis elegans. Nat. Genet. 36, 100–105. doi:10.1038/ng1283
Bhalla, N., and Dernburg, A. F. (2005). A conserved checkpoint monitors meiotic chromosome synapsis in Caenorhabditis elegans. Science 310, 1683–1686. doi:10.1126/science.1117468
Blythe, A. J., Fox, A. H., and Bond, C. S. (2016). The ins and outs of lncRNA structure: How, why and what comes next? Biochim. Biophys. Acta 1859, 46–58.
Bohr, T., Ashley, G., Eggleston, E., Firestone, K., and Bhalla, N. (2016). Synaptonemal complex components are required for meiotic checkpoint function in Caenorhabditis elegans. Genetics 204, 987–997. doi:10.1534/genetics.116.191494
Brenner, S. (1974). The genetics of Caenorhabditis elegans. Genetics 77, 71–94. doi:10.1093/genetics/77.1.71
Chadourne, M., Poumerol, E., Jouneau, L., Passet, B., Castille, J., Sellem, E., et al. (2021). Structural and functional characterization of a testicular long non-coding RNA (4930463O16Rik) identified in the meiotic arrest of the mouse topaz1 (-/-) testes. Front. Cell Dev. Biol. 9, 700290.
Chu, D. S., and Shakes, D. C. (2013). Spermatogenesis. Adv. Exp. Med. Biol. 757, 171–203. doi:10.1007/978-1-4614-4015-4_7
Colaiacovo, M. P., Macqueen, A. J., Martinez-Perez, E., Mcdonald, K., Adamo, A., La Volpe, A., et al. (2003). Synaptonemal complex assembly in C. elegans is dispensable for loading strand-exchange proteins but critical for proper completion of recombination. Dev. Cell 5, 463–474. doi:10.1016/s1534-5807(03)00232-6
Couteau, F., and Zetka, M. (2005). HTP-1 coordinates synaptonemal complex assembly with homolog alignment during meiosis in C. elegans. Genes Dev. 19, 2744–2756. doi:10.1101/gad.1348205
Deniz, E., and Erman, B. (2017). Long noncoding RNA (lincRNA), a new paradigm in gene expression control. Funct. Integr. Genomics 17, 135–143.
Ellis, R. E., and Stanfield, G. M. (2014). The regulation of spermatogenesis and sperm function in nematodes. Semin. Cell Dev. Biol. 29, 17–30. doi:10.1016/j.semcdb.2014.04.005
Fatica, A., and Bozzoni, I. (2014). Long non-coding RNAs: New players in cell differentiation and development. Nat. Rev. Genet. 15, 7–21.
Fico, A., Fiorenzano, A., Pascale, E., Patriarca, E. J., and Minchiotti, G. (2019). Long non-coding RNA in stem cell pluripotency and lineage commitment: Functions and evolutionary conservation. Cell Mol. Life Sci. 76, 1459–1471.
Griswold, M. D. (2016). Spermatogenesis: The commitment to meiosis. Physiol. Rev. 96, 1–17. doi:10.1152/physrev.00013.2015
Hayashi, M., Mlynarczyk-Evans, S., and Villeneuve, A. M. (2010). The synaptonemal complex shapes the crossover landscape through cooperative assembly, crossover promotion and crossover inhibition during Caenorhabditis elegans meiosis. Genetics 186, 45–58. doi:10.1534/genetics.110.115501
Hezroni, H., Koppstein, D., Schwartz, M. G., Avrutin, A., Bartel, D. P., and Ulitsky, I. (2015). Principles of long noncoding RNA evolution derived from direct comparison of transcriptomes in 17 species. Cell Rep. 11, 1110–1122. doi:10.1016/j.celrep.2015.04.023
Hillers, K. J., Jantsch, V., Martinez-Perez, E., and Yanowitz, J. L. (2017). Meiosis. WormBook, 1–43.
Hodgkin, J. (1986). Sex determination in the nematode C. elegans: Analysis of tra-3 suppressors and characterization of fem genes. Genetics 114, 15–52. doi:10.1093/genetics/114.1.15
Hong, S. H., Kwon, J. T., Kim, J., Jeong, J., Kim, J., Lee, S., et al. (2018). Profiling of testis-specific long noncoding RNAs in mice. BMC Genomics 19, 539. doi:10.1186/s12864-018-4931-3
Hubbard, E. J. A., and Schedl, T. (2019). Biology of the Caenorhabditis elegans germline stem cell system. Genetics 213, 1145–1188. doi:10.1534/genetics.119.300238
Ishtayeh, H., Achache, H., Kroizer, E., Rappaport, Y., Itskovits, E., Gingold, H., et al. (2021). Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth. RNA Biol. 18, 435–445. doi:10.1080/15476286.2020.1814549
Jaramillo-Lambert, A., Harigaya, Y., Vitt, J., Villeneuve, A., and Engebrecht, J. (2010). Meiotic errors activate checkpoints that improve gamete quality without triggering apoptosis in male germ cells. Curr. Biol. 20, 2078–2089. doi:10.1016/j.cub.2010.10.008
Kadandale, P., and Singson, A. (2004). Oocyte production and sperm utilization patterns in semi-fertile strains of Caenorhabditis elegans. BMC Dev. Biol. 4, 3. doi:10.1186/1471-213X-4-3
Kelly, W. G., Schaner, C. E., Dernburg, A. F., Lee, M. H., Kim, S. K., Villeneuve, A. M., et al. (2002). X-chromosome silencing in the germline of C. elegans. Development 129, 479–492. doi:10.1242/dev.129.2.479
Kopp, F., and Mendell, J. T. (2018). Functional classification and experimental dissection of long noncoding RNAs. Cell 172, 393–407.
L'hernault, S. W. (1997). “Spermatogenesis,” in C. elegans II. Editors D. L. Riddle, T. Blumenthal, B. J. Meyer, and J. R. Priess 2nd ed (Cold Spring Harbor). (NY)).
Lamunyon, C. W., and Ward, S. (1998). Larger sperm outcompete smaller sperm in the nematode Caenorhabditis elegans. Proc. Biol. Sci. 265, 1997–2002. doi:10.1098/rspb.1998.0531
Lhernault, S. W., Shakes, D. C., and Ward, S. (1988). Developmental genetics of chromosome I spermatogenesis-defective mutants in the nematode Caenorhabditis elegans. Genetics 120, 435–452. doi:10.1093/genetics/120.2.435
Lhernault, S. W. (2009). The genetics and cell biology of spermatogenesis in the nematode C. elegans. Mol. Cell Endocrinol. 306, 59–65. doi:10.1016/j.mce.2009.01.008
Li, P. Y., Li, S. Q., Gao, S. G., and Dong, D. Y. (2021). CRISPR/Cas9-mediated gene editing on Sox2ot promoter leads to its truncated expression and does not influence neural tube closure and embryonic development in mice. Biochem. Biophys. Res. Commun. 573, 107–111.
Lui, D. Y., and Colaiacovo, M. P. (2013). Meiotic development in Caenorhabditis elegans. Adv. Exp. Med. Biol. 757, 133–170. doi:10.1007/978-1-4614-4015-4_6
Macqueen, A. J., Colaiacovo, M. P., Mcdonald, K., and Villeneuve, A. M. (2002). Synapsis-dependent and -independent mechanisms stabilize homolog pairing during meiotic prophase in C. elegans. Genes Dev. 16, 2428–2442. doi:10.1101/gad.1011602
Marques, A. C., and Ponting, C. P. (2014). Intergenic lncRNAs and the evolution of gene expression. Curr. Opin. Genet. Dev. 27, 48–53.
Mehta, S. L., Chokkalla, A. K., Kim, T., Bathula, S., Chelluboina, B., Morris-Blanco, K. C., et al. (2021). Long noncoding RNA fos downstream transcript is developmentally dispensable but vital for shaping the poststroke functional outcome. Stroke 52, 2381–2392.
Melissari, M. T., and Grote, P. (2016). Roles for long non-coding RNAs in physiology and disease. Pflugers Arch. 468, 945–958.
Nam, J. W., and Bartel, D. P. (2012). Long noncoding RNAs in C. elegans. Genome Res. 22, 2529–2540. doi:10.1101/gr.140475.112
Necsulea, A., Soumillon, M., Warnefors, M., Liechti, A., Daish, T., Zeller, U., et al. (2014). The evolution of lncRNA repertoires and expression patterns in tetrapods. Nature 505, 635–640. doi:10.1038/nature12943
Nishimura, H., and Lhernault, S. W. (2017). Spermatogenesis. Curr. Biol. 27, R988–R994. doi:10.1016/j.cub.2017.07.067
Phillips, C. M., and Dernburg, A. F. (2006). A family of zinc-finger proteins is required for chromosome-specific pairing and synapsis during meiosis in C. elegans. Dev. Cell 11, 817–829. doi:10.1016/j.devcel.2006.09.020
Pickett, C. L., Dietrich, N., Chen, J., Xiong, C., and Kornfeld, K. (2013). Mated progeny production is a biomarker of aging in Caenorhabditis elegans. G3 (Bethesda) 3, 2219–2232. doi:10.1534/g3.113.008664
Rappaport, Y., Falk, R., Achache, H., and Tzur, Y. B. (2022). linc-20 and linc-9 do not have compensatory fertility roles in C. elegans. Micropubl. Biol. 2022. doi:10.17912/micropub.biology.000524
Reichman, R., Alleva, B., and Smolikove, S. (2017). Prophase I: Preparing chromosomes for segregation in the developing oocyte. Results Probl. Cell Differ. 59, 125–173. doi:10.1007/978-3-319-44820-6_5
Schild-Prufert, K., Saito, T. T., Smolikov, S., Gu, Y., Hincapie, M., Hill, D. E., et al. (2011). Organization of the synaptonemal complex during meiosis in Caenorhabditis elegans. Genetics 189, 411–421. doi:10.1534/genetics.111.132431
Shakes, D. C., and Ward, S. (1989). Initiation of spermiogenesis in C. elegans: A pharmacological and genetic analysis. Dev. Biol. 134, 189–200. doi:10.1016/0012-1606(89)90088-2
Shakes, D. C., Wu, J. C., Sadler, P. L., Laprade, K., Moore, L. L., Noritake, A., et al. (2009). Spermatogenesis-specific features of the meiotic program in Caenorhabditis elegans. PLoS Genet. 5, e1000611. doi:10.1371/journal.pgen.1000611
Shields, E. J., Petracovici, A. F., and Bonasio, R. (2019). lncRedibly versatile: biochemical and biological functions of long noncoding RNAs. Biochem. J. 476, 1083–1104.
Smolikov, S., Eizinger, A., Hurlburt, A., Rogers, E., Villeneuve, A. M., and Colaiacovo, M. P. (2007a). Synapsis-defective mutants reveal a correlation between chromosome conformation and the mode of double-strand break repair during Caenorhabditis elegans meiosis. Genetics 176, 2027–2033. doi:10.1534/genetics.107.076968
Smolikov, S., Eizinger, A., Schild-Prufert, K., Hurlburt, A., Mcdonald, K., Engebrecht, J., et al. (2007b). SYP-3 restricts synaptonemal complex assembly to bridge paired chromosome axes during meiosis in Caenorhabditis elegans. Genetics 176, 2015–2025. doi:10.1534/genetics.107.072413
Soumillon, M., Necsulea, A., Weier, M., Brawand, D., Zhang, X., Gu, H., et al. (2013). Cellular source and mechanisms of high transcriptome complexity in the mammalian testis. Cell Rep. 3, 2179–2190. doi:10.1016/j.celrep.2013.05.031
Stanfield, G. M., and Villeneuve, A. M. (2006). Regulation of sperm activation by SWM-1 is required for reproductive success of C. elegans males. Curr. Biol. 16, 252–263. doi:10.1016/j.cub.2005.12.041
Sulston, J. E., Schierenberg, E., White, J. G., and Thomson, J. N. (1983). The embryonic cell lineage of the nematode Caenorhabditis elegans. Dev. Biol. 100, 64–119. doi:10.1016/0012-1606(83)90201-4
Tzur, Y. B. (2022). lncRNAs in fertility: redefining the gene expression paradigm? Trends Genet. 38, 1170–1179. doi:10.1016/j.tig.2022.05.013
Tzur, Y. B., Winter, E., Gao, J., Hashimshony, T., Yanai, I., and Colaiacovo, M. P. (2018). Spatiotemporal gene expression analysis of the Caenorhabditis elegans germline uncovers a syncytial expression switch. Genetics 210, 587–605. doi:10.1534/genetics.118.301315
Ulitsky, I., and Bartel, D. P. (2013). lincRNAs: genomics, evolution, and mechanisms. Cell 154, 26–46.
Uszczynska-Ratajczak, B., Lagarde, J., Frankish, A., Guigo, R., and Johnson, R. (2018). Towards a complete map of the human long non-coding RNA transcriptome. Nat. Rev. Genet. 19, 535–548. doi:10.1038/s41576-018-0017-y
Washietl, S., Kellis, M., and Garber, M. (2014). Evolutionary dynamics and tissue specificity of human long noncoding RNAs in six mammals. Genome Res. 24, 616–628. doi:10.1101/gr.165035.113
Wen, K., Yang, L., Xiong, T., Di, C., Ma, D., Wu, M., et al. (2016). Critical roles of long noncoding RNAs in Drosophila spermatogenesis. Genome Res. 26, 1233–1244. doi:10.1101/gr.199547.115
Ye, A. L., Ragle, J. M., Conradt, B., and Bhalla, N. (2014). Differential regulation of germline apoptosis in response to meiotic checkpoint activation. Genetics 198, 995–1000. doi:10.1534/genetics.114.170241
Yoshida, S. (2010). Stem cells in mammalian spermatogenesis. Dev. Growth Differ. 52, 311–317. doi:10.1111/j.1440-169X.2010.01174.x
Keywords: lncRNA, lincRNA, long intergenic non-coding RNA, spermatogenesis, fertility, C. elegans
Citation: Shabtai R and Tzur YB (2023) Male-specific roles of lincRNA in C. elegans fertility. Front. Cell Dev. Biol. 11:1115605. doi: 10.3389/fcell.2023.1115605
Received: 04 December 2022; Accepted: 08 March 2023;
Published: 23 March 2023.
Edited by:
Barak Rotblat, Ben-Gurion University of the Negev, IsraelReviewed by:
Gabriel Leprivier, Heinrich Heine University of Düsseldorf, GermanyMarcia Manterola, University of Chile, Chile
Copyright © 2023 Shabtai and Tzur. 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: Yonatan B. Tzur, tzur@mail.huji.ac.il