Spatiotemporal Gene Expression Analysis of the Caenorhabditis elegans Germline Uncovers a Syncytial Expression Switch

Spatiotemporal single-cell architecture of gene expression in the Caenorhabditis elegans germ cells

Spermatogenesis is an intricate and tightly controlled process encompassing various layers of gene expression regulation. Despite the advance of our current understanding, the developmental trajectory and regulatory mechanisms dictating spermatogenesis remain elusive. In this study, we have generated single-cell gene expression profiles for Caenorhabditis elegans sperm cells and constructed gene regulatory networks alongside the developmental trajectories of these cells. Our findings indicate that each pre- and post-developmental stage is closely linked by co-expressed genes, while simultaneously being uniquely identified by the combined expression of specific gene families. To illustrate the applicability of this exhaustive gene expression catalog, we used gene regulatory networks to uncover potential transcription factors for (1) the expression of genes in the phosphorylation pathway, identifying NHR-23-to-phosphatase regulation for the meiotic cell division process; and (2) the expression of constituent components of small RNA pathways, identifying ELT-1-to-Argonaute protein regulation for siRNA maintenance and sperm activation. We expect that this sperm cell-specific gene expression directory will prompt investigations into the underlying mechanisms determining anatomy, differentiation, and function across the reproductive system. Finally, our expression data can be explored using the web application CelegansGermAtlas ( https://scgerm-atlas.sjtu.edu.cn/website/#/home ).

Parental-effect gene-drive elements under partial selfing, or why do Caenorhabditis genomes have hyperdivergent regions?

Self-fertile Caenorhabditis nematodes carry a surprising number of Medea elements, alleles that act in heterozygous mothers and cause death or developmental delay in offspring that do not inherit them. At some loci, both alleles in a cross operate as independent Medeas, affecting all the homozygous progeny of a selfing heterozygote. The genomic coincidence of Medea elements and ancient, deeply coalescing haplotypes, which pepper the otherwise homogeneous genomes of these animals, raises questions about how these apparent gene-drive elements persist for long periods of time. Here, I investigate how mating system affects the evolution of Medeas, and their paternal-effect counterparts, peels. Despite an intuition that antagonistic alleles should induce balancing selection by killing homozygotes, models show that, under partial selfing, antagonistic elements experience positive frequency dependence: the common allele drives the rare one extinct, even if the rare one is more penetrant. Analytical results for the threshold frequency required for one allele to invade a population show that a very weakly penetrant allele, one whose effects would escape laboratory detection, could nevertheless prevent a much more penetrant allele from invading under high rates of selfing. Ubiquitous weak antagonistic Medeas and peels could then act as localized barriers to gene flow between populations, generating genomic islands of deep coalescence. Analysis of gene expression data, however, suggests that this cannot be the whole story. A complementary explanation is that ordinary ecological balancing selection generates ancient haplotypes on which Medeas can evolve, while high homozygosity in these selfers minimizes the role of gene drive in their evolution.

Role of male gonad-enriched microRNAs in sperm production in Caenorhabditis elegans

Germ cell development and gamete production in animals require small RNA pathways. While studies indicate that microRNAs (miRNAs) are necessary for normal sperm production and function, the specific roles for individual miRNAs are largely unknown. Here, we use small RNA sequencing (RNA-seq) of dissected gonads and functional analysis of new loss-of-function alleles to identify functions for miRNAs in the control of fecundity and sperm production in Caenorhabditis elegans males and hermaphrodites. We describe a set of 29 male gonad-enriched miRNAs and identify a set of individual miRNAs (mir-58.1 and mir-235) and a miRNA cluster (mir-4807-4810.1) that are required for optimal sperm production at 20°C and a set of miRNAs (mir-49, mir-57, mir-83, mir-261, and mir-357/358) that are required for sperm production at 25°C. We observed defects in meiotic progression in mutants missing mir-58.1, mir-83, mir-235, and mir-4807-4810.1, which may contribute to the observed defects in sperm production. Further, analysis of multiple mutants of these miRNAs suggested genetic interactions between these miRNAs. This study provides insights on the regulatory roles of miRNAs that promote optimal sperm production and fecundity in males and hermaphrodites.

ADBP-1 regulates ADR-2 nuclear localization to control editing substrate selection

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is a prevalent and conserved RNA modification. While A-to-I RNA editing is essential in mammals, in Caenorhabditis elegans, it is not, making them invaluable for RNA editing research. In C. elegans, ADR-2 is the sole catalytic A-to-I editing enzyme, and ADR-1 is an RNA editing regulator. ADAR localization is well-studied in humans but not well-established in C. elegans. In this study, we examine the cellular and tissue-specific localization of ADR-2. We show that while ADR-2 is present in most cells in the embryo, at later developmental stages, its expression is both tissue- and cell-type-specific. Additionally, both ADARs are mainly in the nucleus. ADR-2 is adjacent to the chromosomes during the cell cycle. We show that the nuclear localization of endogenous ADR-2 depends on ADBP-1, not ADR-1. In adbp-1 mutant worms, ADR-2 is mislocalized, while ADR-1 is not, leading to decreased editing levels and de-novo editing, mostly in exons, suggesting that ADR-2 is also functional in the cytoplasm. Besides, mutated ADBP-1 affects gene expression. Furthermore, we show that ADR-2 targets adenosines with different surrounding nucleotides in exons and introns. Our findings indicate that ADR-2 cellular localization is highly regulated and affects its function.

Transcriptomic analysis of the spatiotemporal axis of oogenesis and fertilization in C. elegans

Caenorhabditis elegans hermaphrodite presents a unique model to study the formation of oocytes. However, the size of the model animal and difficulties in retrieval of specific stages of the germline have obviated closer systematic studies of this process throughout the years. Here, we present a transcriptomic level analysis into the oogenesis of C. elegans hermaphrodites. We dissected a hermaphrodite gonad into seven sections corresponding to the mitotic distal region, the pachytene region, the diplotene region, the early diakinesis region and the 3 most proximal oocytes, and deeply sequenced the transcriptome of each of them along with that of the fertilized egg using a single-cell RNA-seq (scRNA-seq) protocol. We identified specific gene expression events as well as gene splicing events in finer detail along the gonad and provided novel insights into underlying mechanisms of the oogenesis process. Furthermore, through careful review of relevant research literature coupled with patterns observed in our analysis, we delineate transcripts that may serve functions in the interactions between the germline and cells of the somatic gonad. These results expand our knowledge of the transcriptomic space of the C. elegans germline and lay a foundation on which future studies of the germline can be based upon.

A nucleic acid binding protein map of germline regulation in Caenorhabditis elegans

Fertility requires the faithful proliferation of germ cells and their differentiation into gametes. Controlling these cellular states demands precise timing and expression of gene networks. Nucleic acid binding proteins (NBPs) play critical roles in gene expression networks that influence germ cell development. There has, however, been no functional analysis of the entire NBP repertoire in controlling in vivo germ cell development. Here, we analyzed germ cell states and germline architecture to systematically investigate the function of 364 germline-expressed NBPs in the Caenorhabditis elegans germ line. Using germline-specific knockdown, automated germ cell counting, and high-content analysis of germ cell nuclei and plasma membrane organization, we identify 156 NBPs with discrete autonomous germline functions. By identifying NBPs that control the germ cell cycle, proliferation, differentiation, germline structure and fertility, we have created an atlas for mechanistic dissection of germ cell behavior and gamete production.

BRCA1/BRC-1 and SMC-5/6 regulate DNA repair pathway engagement during Caenorhabditis elegans meiosis

The preservation of genome integrity during sperm and egg development is vital for reproductive success. During meiosis, the tumor suppressor BRCA1/BRC-1 and structural maintenance of chromosomes 5/6 (SMC-5/6) complex genetically interact to promote high fidelity DNA double strand break (DSB) repair, but the specific DSB repair outcomes these proteins regulate remain unknown. Using genetic and cytological methods to monitor resolution of DSBs with different repair partners in Caenorhabditis elegans, we demonstrate that both BRC-1 and SMC-5 repress intersister crossover recombination events. Sequencing analysis of conversion tracts from homolog-independent DSB repair events further indicates that BRC-1 regulates intersister/intrachromatid noncrossover conversion tract length. Moreover, we find that BRC-1 specifically inhibits error prone repair of DSBs induced at mid-pachytene. Finally, we reveal functional interactions of BRC-1 and SMC-5/6 in regulating repair pathway engagement: BRC-1 is required for localization of recombinase proteins to DSBs in smc-5 mutants and enhances DSB repair defects in smc-5 mutants by repressing theta-mediated end joining (TMEJ). These results are consistent with a model in which some functions of BRC-1 act upstream of SMC-5/6 to promote recombination and inhibit error-prone DSB repair, while SMC-5/6 acts downstream of BRC-1 to regulate the formation or resolution of recombination intermediates. Taken together, our study illuminates the coordinated interplay of BRC-1 and SMC-5/6 to regulate DSB repair outcomes in the germline.

Parental-effect gene-drive elements under partial selfing, or why do Caenorhabditis genomes have hyperdivergent regions?

Self-fertile Caenorhabditis nematodes carry a surprising number of Medea elements, alleles that act in heterozygous mothers and cause death or developmental delay in offspring that don't inherit them. At some loci, both alleles in a cross operate as independent Medeas, affecting all the homozygous progeny of a selfing heterozygote. The genomic coincidence of Medea elements and ancient, deeply coalescing haplotypes, which pepper the otherwise homogeneous genomes of these animals, raises questions about how these apparent gene-drive elements persist for long periods of time. Here I investigate how mating system affects the evolution of Medeas, and their paternal-effect counterparts, peels. Despite an intuition that antagonistic alleles should induce balancing selection by killing homozygotes, models show that, under partial selfing, antagonistic elements experience positive frequency dependence: the common allele drives the rare one extinct, even if the rare one is more penetrant. Analytical results for the threshold frequency required for one allele to invade a population show that a very weakly penetrant allele, one whose effects would escape laboratory detection, could nevertheless prevent a much more penetrant allele from invading under high rates of selfing. Ubiquitous weak antagonistic Medeas and peels could then act as localized barriers to gene flow between populations, generating genomic islands of deep coalescence. Analysis of gene expression data, however, suggest that this cannot be the whole story. A complementary explanation is that ordinary ecological balancing selection generates ancient haplotypes on which Medeas can evolve, while high homozygosity in these selfers minimizes the role of gene drive in their evolution.

ADBP-1 regulates ADR-2 nuclear localization to control editing substrate selection

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by ADAR enzymes, is a prevalent and conserved RNA modification. While A-to-I RNA editing is essential in mammals, in Caenorhabditis elegans , it is not, making them invaluable for RNA editing research. In C. elegans , ADR-2 is the sole catalytic A-to-I editing enzyme, and ADR-1 is an RNA editing regulator. ADAR localization is well-studied in humans but not well-established in C. elegans . In this study, we examine the cellular and tissue-specific localization of ADR-2. We show that while ADR-2 is present in most cells in the embryo, at later developmental stages, its expression is both tissue- and cell-type-specific. Additionally, both ADARs are mainly in the nucleus. ADR-2 is adjacent to the chromosomes during the cell cycle. We show that the nuclear localization of endogenous ADR-2 depends on ADBP-1, not ADR-1. In adbp-1 mutant worms, ADR-2 is mislocalized, while ADR-1 is not, leading to decreased editing levels and de-novo editing, mostly in exons, suggesting that ADR-2 is also functional in the cytoplasm. Besides, mutated ADBP-1 affects gene expression. Furthermore, we show that ADR-2 targets adenosines with different surrounding nucleotides in exons and introns. Our findings indicate that ADR-2 cellular localization is highly regulated and affects its function.

Transcriptomic Analysis of the Spatiotemporal Axis of Oogenesis and Fertilization in C. elegans

The oocyte germline of the C. elegans hermaphrodite presents a unique model to study the formation of oocytes. However, the size of the model animal and difficulties in retrieval of specific stages of the germline have obviated closer systematic studies of this process throughout the years. Here, we present a transcriptomic level analysis into the oogenesis of C. elegans hermaphrodites. We dissected a hermaphrodite gonad into seven sections corresponding to the mitotic distal region, the pachytene, the diplotene, the early diakinesis region and the 3 most proximal oocytes, and deeply sequenced the transcriptome of each of them along with that of the fertilized egg using a single-cell RNA-seq protocol. We identified specific gene expression events as well as gene splicing events in finer detail along the oocyte germline and provided novel insights into underlying mechanisms of the oogenesis process. Furthermore, through careful review of relevant research literature coupled with patterns observed in our analysis, we attempt to delineate transcripts that may serve functions in the interaction between the germline and cells of the somatic gonad. These results expand our knowledge of the transcriptomic space of the C. elegans germline and lay a foundation on which future studies of the germline can be based upon.

The matrisome landscape controlling in vivo germ cell fates

The developmental fate of cells is regulated by intrinsic factors and the extracellular environment. The extracellular matrix (matrisome) delivers chemical and mechanical cues that can modify cellular development. However, comprehensive understanding of how matrisome factors control cells in vivo is lacking. Here we show that specific matrisome factors act individually and collectively to control germ cell development. Surveying development of undifferentiated germline stem cells through to mature oocytes in the Caenorhabditis elegans germ line enabled holistic functional analysis of 443 conserved matrisome-coding genes. Using high-content imaging, 3D reconstruction, and cell behavior analysis, we identify 321 matrisome genes that impact germ cell development, the majority of which (>80%) are undescribed. Our analysis identifies key matrisome networks acting autonomously and non-autonomously to coordinate germ cell behavior. Further, our results demonstrate that germ cell development requires continual remodeling of the matrisome landscape. Together, this study provides a comprehensive platform for deciphering how extracellular signaling controls cellular development and anticipate this will establish new opportunities for manipulating cell fates.

SIN-3 acts in distinct complexes to regulate the germline transcriptional program in Caenorhabditis elegans

The transcriptional co-regulator SIN3 influences gene expression through multiple interactions that include histone deacetylases. Haploinsufficiency and mutations in SIN3 are the underlying cause of Witteveen-Kolk syndrome and related intellectual disability and autism syndromes, emphasizing its key role in development. However, little is known about the diversity of its interactions and functions in developmental processes. Here, we show that loss of SIN-3, the single SIN3 homolog in Caenorhabditis elegans, results in maternal-effect sterility associated with de-regulation of the germline transcriptome, including de-silencing of X-linked genes. We identify at least two distinct SIN3 complexes containing specific histone deacetylases and show that they differentially contribute to fertility. Single-cell, single-molecule fluorescence in situ hybridization reveals that in sin-3 mutants the X chromosome becomes re-expressed prematurely and in a stochastic manner in individual germ cells, suggesting a role for SIN-3 in its silencing. Furthermore, we identify histone residues whose acetylation increases in the absence of SIN-3. Together, this work provides a powerful framework for the in vivo study of SIN3 and associated proteins.

Male gonad-enriched microRNAs function to control sperm production in C. elegans

Germ cell development and gamete production in animals require small RNA pathways. While studies indicate that microRNAs (miRNAs) are necessary for normal sperm production and function, the specific roles for individual miRNAs are largely unknown. Here, we use small RNA sequencing of dissected gonads and functional analysis of new loss of function alleles to identify functions for miRNAs in the control of fecundity and sperm production in Caenorhabditis elegans males and hermaphrodites. We describe a set of 29 male gonad-enriched miRNAs and identify a set of 3 individual miRNAs (mir-58.1, mir-83, and mir-235) and a miRNA cluster (mir-4807-4810.1) that are required for optimal sperm production at 20°C and 5 additional miRNAs (mir-49, mir-57, mir-261, and mir-357/358) that are required for sperm production at 25°C. We observed defects in meiotic progression in mir-58.1, mir-83, mir-235, and mir-4807-4810.1 mutants that may contribute to the reduced number of sperm. Further, analysis of multiple mutants of these miRNAs suggested complex genetic interactions between these miRNAs for sperm production. This study provides insights on the regulatory roles of miRNAs that promote optimal sperm production and fecundity in males and hermaphrodites.

Toward Gene-Correlated Spatially Resolved Metabolomics with Fingerprint Coherent Raman Imaging

Raman spectroscopy has long been known to provide sufficient information to discriminate distinct cell phenotypes. Underlying this discriminating capability is that Raman spectra provide an overall readout of the metabolic profiles that change with transcriptomic activity. Robustly associating Raman spectral changes with the regulation of specific signaling pathways may be possible, but the spectral signals of interest may be weak and vary somewhat among individuals. Establishing a Raman-to-transcriptome mapping will thus require tightly controlled and easily manipulated biological systems and high-throughput spectral acquisition. We attempt to meet these requirements using broadband coherent anti-Stokes Raman scattering (BCARS) microscopy to spatio-spectrally map the C. elegans hermaphrodite gonad in vivo at subcellular resolution. The C. elegans hermaphrodite gonad is an ideal model system with a sequential, continuous process of highly regulated spatiotemporal cellular events. We demonstrate that the BCARS spatio-spectral signatures correlate with gene expression profiles in the gonad, evincing that BCARS has potential as a spatially resolved omics surrogate.

Spatial single-cell sequencing of meiosis I arrested oocytes indicates acquisition of maternal transcripts from the soma

Maternal RNAs are stored from minutes to decades in oocytes throughout meiosis I arrest in a transcriptionally quiescent state. Recent reports, however, propose a role for nascent transcription in arrested oocytes. Whether arrested oocytes launch nascent transcription in response to environmental or hormonal signals while maintaining the meiosis I arrest remains undetermined. We test this by integrating single-cell RNA sequencing, RNA velocity, and RNA fluorescence in situ hybridization on C. elegans meiosis I arrested oocytes. We identify transcripts that increase as the arrested meiosis I oocyte ages, but rule out extracellular signaling through ERK MAPK and nascent transcription as a mechanism for this increase. We report transcript acquisition from neighboring somatic cells as a mechanism of transcript increase during meiosis I arrest. These analyses provide a deeper view at single-cell resolution of the RNA landscape of a meiosis I arrested oocyte and as it prepares for oocyte maturation and fertilization.

A comprehensive survey of C. elegans argonaute proteins reveals organism-wide gene regulatory networks and functions

Argonaute (AGO) proteins associate with small RNAs to direct their effector function on complementary transcripts. The nematode Caenorhabditis elegans contains an expanded family of 19 functional AGO proteins, many of which have not been fully characterized. In this work, we systematically analyzed every C. elegans AGO using CRISPR-Cas9 genome editing to introduce GFP::3xFLAG tags. We have characterized the expression patterns of each AGO throughout development, identified small RNA binding complements, and determined the effects of ago loss on small RNA populations and developmental phenotypes. Our analysis indicates stratification of subsets of AGOs into distinct regulatory modules, and integration of our data led us to uncover novel stress-induced fertility and pathogen response phenotypes due to ago loss.

GRAS-1 is a novel regulator of early meiotic chromosome dynamics in C. elegans

Chromosome movements and licensing of synapsis must be tightly regulated during early meiosis to ensure accurate chromosome segregation and avoid aneuploidy, although how these steps are coordinated is not fully understood. Here we show that GRAS-1, the worm homolog of mammalian GRASP/Tamalin and CYTIP, coordinates early meiotic events with cytoskeletal forces outside the nucleus. GRAS-1 localizes close to the nuclear envelope (NE) in early prophase I and interacts with NE and cytoskeleton proteins. Delayed homologous chromosome pairing, synaptonemal complex (SC) assembly, and DNA double-strand break repair progression are partially rescued by the expression of human CYTIP in gras-1 mutants, supporting functional conservation. However, Tamalin, Cytip double knockout mice do not exhibit obvious fertility or meiotic defects, suggesting evolutionary differences between mammals. gras-1 mutants show accelerated chromosome movement during early prophase I, implicating GRAS-1 in regulating chromosome dynamics. GRAS-1-mediated regulation of chromosome movement is DHC-1-dependent, placing it acting within the LINC-controlled pathway, and depends on GRAS-1 phosphorylation at a C-terminal S/T cluster. We propose that GRAS-1 coordinates the early steps of homology search and licensing of SC assembly by regulating the pace of chromosome movement in early prophase I.

Male-specific roles of lincRNA in C. elegans fertility

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.

Chromatin landscape, DSB levels, and cKU-70/80 contribute to patterning of meiotic DSB processing along chromosomes in C. elegans

Programmed DNA double-strand break (DSB) formation is essential for achieving accurate chromosome segregation during meiosis. DSB repair timing and template choice are tightly regulated. However, little is known about how DSB distribution and the choice of repair pathway are regulated along the length of chromosomes, which has direct effects on the recombination landscape and chromosome remodeling at late prophase I. Here, we use the spatiotemporal resolution of meiosis in the Caenorhabditis elegans germline along with genetic approaches to study distribution of DSB processing and its regulation. High-resolution imaging of computationally straightened chromosomes immunostained for the RAD-51 recombinase marking DSB repair sites reveals that the pattern of RAD-51 foci throughout pachytene resembles crossover distribution in wild type. Specifically, RAD-51 foci occur primarily along the gene-poor distal thirds of the chromosomes in both early and late pachytene, and on both the X and the autosomes. However, this biased off-center distribution can be abrogated by the formation of excess DSBs. Reduced condensin function, but not an increase in total physical axial length, results in a homogeneous distribution of RAD-51 foci, whereas regulation of H3K9 methylation is required for the enrichment of RAD-51 at off-center positions. Finally, the DSB recognition heterodimer cKU-70/80, but not the non-homologous end-joining canonical ligase LIG-4, contributes to the enriched off-center distribution of RAD-51 foci. Taken together, our data supports a model by which regulation of the chromatin landscape, DSB levels, and DSB detection by cKU-70/80 collaborate to promote DSB processing by homologous recombination at off-center regions of the chromosomes in C. elegans.

lncRNAs in fertility: redefining the gene expression paradigm?

Comparative transcriptome approaches assume that highly or dynamically expressed genes are important. This has led to the identification of many genes critical for cellular activity and organism development. However, while testes express the highest levels of long noncoding RNAs (lncRNAs), there is scarcely any evidence for lncRNAs with significant roles in fertility. This was explained by changes in chromatin structure during spermatogenesis that lead to 'promiscuous transcription' with no functional roles for the transcripts. Recent discoveries offer novel and surprising alternatives. Here, I review the current knowledge regarding the involvement of lncRNAs in fertility, why I find gametogenesis different from other developmental processes, offer models to explain why the experimental evidence did not meet theoretical predictions, and suggest possible approaches to test the models.

Sperm-inherited H3K27me3 epialleles are transmitted transgenerationally in cis

The transmission of chromatin states from parent cells to daughter cells preserves cell-specific transcriptional states and thus cell identity through cell division. The mechanism that underpins this process is not fully understood. The role that chromatin states serve in transmitting gene expression information across generations via sperm and oocytes is even less understood. Here, we utilized a model in which Caenorhabditis elegans sperm and oocyte alleles were inherited in different states of the repressive mark H3K27me3. This resulted in the alleles achieving different transcriptional states within the nuclei of offspring. Using this model, we showed that sperm alleles inherited without H3K27me3 were sensitive to up-regulation in offspring somatic and germline tissues, and tissue context determined which genes were up-regulated. We found that the subset of sperm alleles that were up-regulated in offspring germlines retained the H3K27me3(-) state and were transmitted to grandoffspring as H3K27me3(-) and up-regulated epialleles, demonstrating that H3K27me3 can serve as a transgenerational epigenetic carrier in C. elegans.

Maternal H3K36 and H3K27 HMTs protect germline development via regulation of the transcription factor LIN-15B

Maternally synthesized products play critical roles in the development of offspring. A premier example is the Caenorhabditis elegans H3K36 methyltransferase MES-4, which is essential for germline survival and development in offspring. How maternal MES-4 protects the germline is not well understood, but its role in H3K36 methylation hinted that it may regulate gene expression in primordial germ cells (PGCs). We tested this hypothesis by profiling transcripts from nascent germlines (PGCs and their descendants) dissected from wild-type and mes-4 mutant (lacking maternal and zygotic MES-4) larvae. mes-4 nascent germlines displayed downregulation of some germline genes, upregulation of some somatic genes, and dramatic upregulation of hundreds of genes on the X chromosome. We demonstrated that upregulation of one or more genes on the X is the cause of germline death by generating and analyzing mes-4 mutants that inherited different endowments of X chromosome(s). Intriguingly, removal of the THAP transcription factor LIN-15B from mes-4 mutants reduced X misexpression and prevented germline death. lin-15B is X-linked and misexpressed in mes-4 PGCs, identifying it as a critical target for MES-4 repression. The above findings extend to the H3K27 methyltransferase MES-2/3/6, the C. elegans version of polycomb repressive complex 2. We propose that maternal MES-4 and PRC2 cooperate to protect germline survival by preventing synthesis of germline-toxic products encoded by genes on the X chromosome, including the key transcription factor LIN-15B.

Small RNAs and their protein partners in animal meiosis

Meiosis is characterized by highly regulated transitions in gene expression that require diverse mechanisms of gene regulation. For example, in male mammals, transcription undergoes a global shut-down in early prophase I of meiosis, followed by increasing transcriptional activity into pachynema. Later, as spermiogenesis proceeds, the histones bound to DNA are replaced with transition proteins, which are themselves replaced with protamines, resulting in a highly condensed nucleus with repressed transcriptional activity. In addition, two specialized gene silencing events take place during prophase I: meiotic silencing of unsynapsed chromatin (MSUC), and the sex chromatin specific mechanism, meiotic sex chromosome inactivation (MSCI). Notably, conserved roles for the RNA binding protein (RBP) machinery that functions with small non-coding RNAs have been described as participating in these meiosis-specific mechanisms, suggesting that RNA-mediated gene regulation is critical for fertility in many species. Here, we review roles of small RNAs and their associated RBPs in meiosis-related processes such as centromere function, silencing of unpaired chromatin and meiotic recombination. We will discuss the emerging evidence of non-canonical functions of these components in meiosis.

R-loop-induced irreparable DNA damage evades checkpoint detection in the C. elegans germline

Accumulation of DNA–RNA hybrids in the form of R-loops can result in replication–transcription conflict that leads to the formation of DNA double strand breaks (DSBs). Using null mutants for the two Caenorhabditis elegans genes encoding for RNaseH1 and RNaseH2, we identify novel effects of R-loop accumulation in the germline. R-loop accumulation leads, as expected, to replication stress, followed by the formation of DSBs. A subset of these DSBs are irreparable. However, unlike irreparable DSBs generated in other systems, which trigger permanent cell cycle arrest, germline irreparable DSBs are propagated to oocytes. Despite DNA damage checkpoint activation in the stem cell niche, the signaling cannot be sustained and nuclei with irreparable DNA damage progress into meiosis. Moreover, unlike other forms of DNA damage that increase germline apoptosis, R-loop-generated DSBs remain undetected by the apoptotic checkpoint. This coincides with attenuation of ATM/ATR signaling in mid-to-late meiotic prophase I. These data altogether indicate that in the germline, DSBs that are generated by R-loops can lead to irreparable DSBs that evade cellular machineries designed for damage recognition. These studies implicate germline R-loops as an especially dangerous driver of germline mutagenesis.

Sex Determination in Nematode Germ Cells

BACKGROUND

Animal germ cells differentiate as sperm or as oocytes. These sexual fates are controlled by complex regulatory pathways to ensure that the proper gametes are made at the appropriate times.

SUMMARY

Nematodes like Caenorhabditis elegans and its close relatives are ideal models for studying how this regulation works, because the XX animals are self-fertile hermaphrodites that produce both sperm and oocytes. In these worms, germ cells use the same signal transduction pathway that functions in somatic cells. This pathway determines the activity of the transcription factor TRA-1, a Gli protein that can repress male genes. However, the pathway is extensively modified in germ cells, largely by the action of translational regulators like the PUF proteins. Many of these modifications play critical roles in allowing the XX hermaphrodites to make sperm in an otherwise female body. Finally, TRA-1 cooperates with chromatin regulators in the germ line to control the activity of fog-1 and fog-3, which are essential for spermatogenesis. FOG-1 and FOG-3 work together to determine germ cell fates by blocking the translation of oogenic transcripts. Key Messages: Although there is great diversity in how germ cell fates are controlled in other animals, many of the key nematode genes are conserved, and the critical role of translational regulators may be universal.

linc-20 and linc-9 do not have compensatory fertility roles in C. elegans

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides which are transcribed from regions that do not overlap with protein coding sequences. Reproductive organs express high levels of lincRNAs, yet removal of many lincRNA genes with high and dynamic germline expression did not lead to fertility defects. It was previously suggested this stems from redundant roles of different lincRNA genes. We previously reported engineering C. elegans strains in which we deleted lincRNA genes with high and dynamic expression in the gonad. The individual mutations did not lead to major effects on fertility. Two of those lincRNA genes, linc-9 and linc-20, are highly homologous, suggesting they could perform redundant roles. Here we report that in the double mutant linc-9; linc-20 the brood size and embryonic lethality do not significantly differ from wild-type worms. This could be explained by either lack of fertility roles, or redundancy with other lincRNA genes.

Small RNA pathways in the nematode Ascaris in the absence of piRNAs

Small RNA pathways play key and diverse regulatory roles in C. elegans, but our understanding of their conservation and contributions in other nematodes is limited. We analyzed small RNA pathways in the divergent parasitic nematode Ascaris. Ascaris has ten Argonautes with five worm-specific Argonautes (WAGOs) that associate with secondary 5’-triphosphate 22-24G-RNAs. These small RNAs target repetitive sequences or mature mRNAs and are similar to the C. elegans mutator, nuclear, and CSR-1 small RNA pathways. Even in the absence of a piRNA pathway, Ascaris CSR-1 may still function to “license” as well as fine-tune or repress gene expression. Ascaris ALG-4 and its associated 26G-RNAs target and likely repress specific mRNAs during testis meiosis. Ascaris WAGO small RNAs demonstrate target plasticity changing their targets between repeats and mRNAs during development. We provide a unique and comprehensive view of mRNA and small RNA expression throughout spermatogenesis. Overall, our study illustrates the conservation, divergence, dynamics, and flexibility of small RNA pathways in nematodes.

The parasitic nematode Ascaris lacks piRNAs. Here the authors compare Argonaute proteins and small RNAs from C. elegans and Ascaris, expanding our understanding of the conservation, divergence, and flexibility of Argonautes and small RNA pathways in nematodes.

mDrop-Seq: Massively Parallel Single-Cell RNA-Seq of Saccharomyces cerevisiae and Candida albicans

Advances in high-throughput single-cell RNA sequencing (scRNA-seq) have been limited by technical challenges such as tough cell walls and low RNA quantity that prevent transcriptomic profiling of microbial species at throughput. We present microbial Drop-seq or mDrop-seq, a high-throughput scRNA-seq technique that is demonstrated on two yeast species, Saccharomyces cerevisiae, a popular model organism, and Candida albicans, a common opportunistic pathogen. We benchmarked mDrop-seq for sensitivity and specificity and used it to profile 35,109 S. cerevisiae cells to detect variation in mRNA levels between them. As a proof of concept, we quantified expression differences in heat shock S. cerevisiae using mDrop-seq. We detected differential activation of stress response genes within a seemingly homogenous population of S. cerevisiae under heat shock. We also applied mDrop-seq to C. albicans cells, a polymorphic and clinically relevant species of yeast with a thicker cell wall compared to S. cerevisiae. Single-cell transcriptomes in 39,705 C. albicans cells were characterized using mDrop-seq under different conditions, including exposure to fluconazole, a common anti-fungal drug. We noted differential regulation in stress response and drug target pathways between C. albicans cells, changes in cell cycle patterns and marked increases in histone activity when treated with fluconazole. We demonstrate mDrop-seq to be an affordable and scalable technique that can quantify the variability in gene expression in different yeast species. We hope that mDrop-seq will lead to a better understanding of genetic variation in pathogens in response to stimuli and find immediate applications in investigating drug resistance, infection outcome and developing new drugs and treatment strategies.

A somatic proteoglycan controls Notch-directed germ cell fate

Communication between the soma and germline optimizes germ cell fate programs. Notch receptors are key determinants of germ cell fate but how somatic signals direct Notch-dependent germ cell behavior is undefined. Here we demonstrate that SDN-1 (syndecan-1), a somatic transmembrane proteoglycan, controls expression of the GLP-1 (germline proliferation-1) Notch receptor in the Caenorhabditis elegans germline. We find that SDN-1 control of a somatic TRP calcium channel governs calcium-dependent binding of an AP-2 transcription factor (APTF-2) to the glp-1 promoter. Hence, SDN-1 signaling promotes GLP-1 expression and mitotic germ cell fate. Together, these data reveal SDN-1 as a putative communication nexus between the germline and its somatic environment to control germ cell fate decisions.

Bisection of the X chromosome disrupts the initiation of chromosome silencing during meiosis in Caenorhabditis elegans

During meiosis, gene expression is silenced in aberrantly unsynapsed chromatin and in heterogametic sex chromosomes. Initiation of sex chromosome silencing is disrupted in meiocytes with sex chromosome-autosome translocations. To determine whether this is due to aberrant synapsis or loss of continuity of sex chromosomes, we engineered Caenorhabditis elegans nematodes with non-translocated, bisected X chromosomes. In early meiocytes of mutant males and hermaphrodites, X segments are enriched with euchromatin assembly markers and active RNA polymerase II staining, indicating active transcription. Analysis of RNA-seq data showed that genes from the X chromosome are upregulated in gonads of mutant worms. Contrary to previous models, which predicted that any unsynapsed chromatin is silenced during meiosis, our data indicate that unsynapsed X segments are transcribed. Therefore, our results suggest that sex chromosome chromatin has a unique character that facilitates its meiotic expression when its continuity is lost, regardless of whether or not it is synapsed.

HIM-17 regulates the position of recombination events and GSP-1/2 localization to establish short arm identity on bivalents in meiosis

The position of recombination events established along chromosomes in early prophase I and the chromosome remodeling that takes place in late prophase I are intrinsically linked steps of meiosis that need to be tightly regulated to ensure accurate chromosome segregation and haploid gamete formation. Here, we show that RAD-51 foci, which form at the sites of programmed meiotic DNA double-strand breaks (DSBs), exhibit a biased distribution toward off-centered positions along the chromosomes in wild-type Caenorhabditis elegans, and we identify two meiotic roles for chromatin-associated protein HIM-17 that ensure normal chromosome remodeling in late prophase I. During early prophase I, HIM-17 regulates the distribution of DSB-dependent RAD-51 foci and crossovers on chromosomes, which is critical for the formation of distinct chromosome subdomains (short and long arms of the bivalents) later during chromosome remodeling. During late prophase I, HIM-17 promotes the normal expression and localization of protein phosphatases GSP-1/2 to the surface of the bivalent chromosomes and may promote GSP-1 phosphorylation, thereby antagonizing Aurora B kinase AIR-2 loading on the long arms and preventing premature loss of sister chromatid cohesion. We propose that HIM-17 plays distinct roles at different stages during meiotic progression that converge to promote normal chromosome remodeling and accurate chromosome segregation.

Sperm fate is promoted by the mir-44 microRNA family in the Caenorhabditis elegans hermaphrodite germline

Posttranscriptional regulation of gene expression, typically effected by RNA-binding proteins, microRNAs (miRNAs), and translation initiation factors, is essential for normal germ cell function. Numerous miRNAs have been detected in the germline; however, the functions of specific miRNAs remain largely unknown. Functions of miRNAs have been difficult to determine as miRNAs often modestly repress target mRNAs and are suggested to sculpt or fine tune gene expression to allow for the robust expression of cell fates. In Caenorhabditis elegans hermaphrodites, cell fate decisions are made for germline sex determination during larval development when sperm are generated in a short window before the switch to oocyte production. Here, analysis of newly generated mir-44 family mutants has identified a family of miRNAs that modulate the germline sex determination pathway in C. elegans. Mutants with the loss of mir-44 and mir-45 produce fewer sperm, showing both a delay in the specification and formation of sperm as well as an early termination of sperm specification accompanied by a premature switch to oocyte production. mir-44 and mir-45 are necessary for the normal period of fog-1 expression in larval development. Through genetic analysis, we find that mir-44 and mir-45 may act upstream of fbf-1 and fem-3 to promote sperm specification. Our research indicates that the mir-44 family promotes sperm cell fate specification during larval development and identifies an additional posttranscriptional regulator of the germline sex determination pathway.

Multivalent weak interactions between assembly units drive synaptonemal complex formation

The synaptonemal complex (SC) is an ordered but highly dynamic structure assembled between homologous chromosomes to control interhomologous crossover formation, ensuring accurate meiotic chromosome segregation. However, the mechanisms regulating SC assembly and dynamics remain unclear. Here, we identified two new SC components, SYP-5 and SYP-6, in Caenorhabditis elegans that have distinct expression patterns and form distinct SC assembly units with other SYPs through stable interactions. SYP-5 and SYP-6 exhibit diverse in vivo SC regulatory functions and distinct phase separation properties in cells. Charge-interacting elements (CIEs) are enriched in SC intrinsically disordered regions (IDRs), and IDR deletion or CIE removal confirmed a requirement for these elements in SC regulation. Our data support the theory that multivalent weak interactions between the SC units drive SC formation and that CIEs confer multivalency to the assembly units.

A Living Organism in your CRISPR Toolbox: Caenorhabditis elegans Is a Rapid and Efficient Model for Developing CRISPR-Cas Technologies

The Cas9 nuclease from Streptococcus pyogenes (SpCas9) is the most popular enzyme for CRISPR technologies. However, considering the wide diversity of microorganisms (discovered and still unknown), a massive number of CRISPR effectors are being and will be identified and characterized in the search of optimal Cas variants for each of the many applications of CRISPR. In this context, a versatile and efficient multicellular system for CRISPR editing such as Caenorhabditis elegans would be of great help in the development of these effectors. Here, we highlight the benefits of using C. elegans for the rapid evaluation of new CRISPR effectors, and for optimizing CRISPR efficiency in animals in several ways such as by modulating the balance between repair pathways, modifying chromatin accessibility, or controlling the expression and activity of nucleases and guide RNAs.

Spatial Transcriptomics of Nematodes Identifies Sperm Cells as a Source of Genomic Novelty and Rapid Evolution

Divergence of gene function and expression during development can give rise to phenotypic differences at the level of cells, tissues, organs, and ultimately whole organisms. To gain insights into the evolution of gene expression and novel genes at spatial resolution, we compared the spatially resolved transcriptomes of two distantly related nematodes, Caenorhabditis elegans and Pristionchus pacificus, that diverged 60-90 Ma. The spatial transcriptomes of adult worms show little evidence for strong conservation at the level of single genes. Instead, regional expression is largely driven by recent duplication and emergence of novel genes. Estimation of gene ages across anatomical structures revealed an enrichment of novel genes in sperm-related regions. This provides first evidence in nematodes for the "out of testis" hypothesis that has been previously postulated based on studies in Drosophila and mammals. "Out of testis" genes represent a mix of products of pervasive transcription as well as fast evolving members of ancient gene families. Strikingly, numerous novel genes have known functions during meiosis in Caenorhabditis elegans indicating that even universal processes such as meiosis may be targets of rapid evolution. Our study highlights the importance of novel genes in generating phenotypic diversity and explicitly characterizes gene origination in sperm-related regions. Furthermore, it proposes new functions for previously uncharacterized genes and establishes the spatial transcriptome of Pristionchus pacificus as a catalog for future studies on the evolution of gene expression and function.

The conserved molting/circadian rhythm regulator NHR-23/NR1F1 serves as an essential co-regulator of C. elegans spermatogenesis

In sexually reproducing metazoans, spermatogenesis is the process by which uncommitted germ cells give rise to haploid sperm. Work in model systems has revealed mechanisms controlling commitment to the sperm fate, but how this fate is subsequently executed remains less clear. While studying the well-established role of the conserved nuclear hormone receptor transcription factor, NHR-23/NR1F1, in regulating C. elegans molting, we discovered that NHR-23/NR1F1 is also constitutively expressed in developing primary spermatocytes and is a critical regulator of spermatogenesis. In this novel role, NHR-23/NR1F1 functions downstream of the canonical sex-determination pathway. Degron-mediated depletion of NHR-23/NR1F1 within hermaphrodite or male germlines causes sterility due to an absence of functional sperm, as depleted animals produce arrested primary spermatocytes rather than haploid sperm. These spermatocytes arrest in prometaphase I and fail to either progress to anaphase or attempt spermatid-residual body partitioning. They make sperm-specific membranous organelles but fail to assemble their major sperm protein into fibrous bodies. NHR-23/NR1F1 appears to function independently of the known SPE-44 gene regulatory network, revealing the existence of an NHR-23/NR1F1-mediated module that regulates the spermatogenesis program.

The molecular underpinnings of fertility: Genetic approaches in Caenorhabditis elegans

The study of mutations that impact fertility has a catch-22. Fertility mutants are often lost since they cannot simply be propagated and maintained. This has hindered progress in understanding the genetics of fertility. In mice, several molecules are found to be required for the interactions between the sperm and egg, with JUNO and IZUMO1 being the only known receptor pair on the egg and sperm surface, respectively. In Caenorhabditis elegans, a total of 12 proteins on the sperm or oocyte have been identified to mediate gamete interactions. Majority of these genes were identified through mutants isolated from genetic screens. In this review, we summarize the several key screening strategies that led to the identification of fertility mutants in C. elegans and provide a perspective about future research using genetic approaches. Recently, advancements in new technologies such as high-throughput sequencing and Crispr-based genome editing tools have accelerated the molecular, cell biological, and mechanistic analysis of fertility genes. We review how these valuable tools advance our understanding of the molecular underpinnings of fertilization. We draw parallels of the molecular mechanisms of fertilization between worms and mammals and argue that our work in C. elegans complements fertility research in humans and other species.

Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides that are transcribed from non-coding loci yet undergo biosynthesis similar to coding mRNAs. The disproportional number of lincRNAs expressed in testes suggests that lincRNAs are important during gametogenesis, but experimental evidence has implicated very few lincRNAs in this process. We took advantage of the relatively limited number of lincRNAs in the genome of the nematode Caenorhabditis elegans to systematically analyse the functions of lincRNAs during meiosis. We deleted six lincRNA genes that are highly and dynamically expressed in the C. elegans gonad and tested the effects on central meiotic processes. Surprisingly, whereas the lincRNA deletions did not strongly impact fertility, germline apoptosis, crossovers, or synapsis, linc-4 was required for somatic growth. Slower growth was observed in linc-4-deletion mutants and in worms depleted of linc-4 using RNAi, indicating that linc-4 transcripts are required for this post-embryonic process. Unexpectedly, analysis of worms depleted of linc-4 in soma versus germline showed that the somatic role stems from linc-4 expression in germline cells. This unique feature suggests that some lincRNAs, like some small non-coding RNAs, are required for germ-soma interactions.

Identification of transcriptome differences in goat ovaries at the follicular phase and the luteal phase using an RNA-Seq method

The ovaries, the main female reproductive organs, directly mediate ovulation and reproductive hormone secretion. These complex physiological processes are regulated by multiple genes and pathways. However, there is a lack of research on goat ovaries, and the molecular mechanisms underlying the signaling pathways remain unclear. In this study, Illumina HiSeq 4000 sequencing was used to sequence the transcriptomes of goat ovaries. The expression patterns of differentially expressed mRNAs in goat ovaries at both the follicular and luteal phases were determined by bioinformatics analysis. A total of 1,122, 014, 112 clean reads were obtained, and 3770 differentially expressed mRNAs were identified for further analysis. There were 1727 and 2043 upregulated mRNAs in the luteal phase and follicular phase, respectively. According to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, some mRNAs that were highly expressed in ovaries during the luteal phase, such as HSD17B7, 3BHSD, and SRD5A2, may be related to the synthesis of progesterone. In addition, some mRNAs that were highly expressed in ovaries during the follicular phase, such as RPL12, RPS13 and RPL10, are related to the growth and maturation of oocytes. Taken together, the findings of this study provide genome-wide mRNA expression profiles for goat ovaries at the follicular and luteal phases and identify mRNAs associated with goat hormone secretion and follicular development. In addition, this study provides a theoretical basis for further investigation of goat reproductive regulation.

X Chromosome Domain Architecture Regulates Caenorhabditis elegans Lifespan but Not Dosage Compensation

Mechanisms establishing higher-order chromosome structures and their roles in gene regulation are elusive. We analyzed chromosome architecture during nematode X chromosome dosage compensation, which represses transcription via a dosage-compensation condensin complex (DCC) that binds hermaphrodite Xs and establishes megabase-sized topologically associating domains (TADs). We show that DCC binding at high-occupancy sites (rex sites) defines eight TAD boundaries. Single rex deletions disrupted boundaries, and single insertions created new boundaries, demonstrating that a rex site is necessary and sufficient to define DCC-dependent boundary locations. Deleting eight rex sites (8rexΔ) recapitulated TAD structure of DCC mutants, permitting analysis when chromosome-wide domain architecture was disrupted but most DCC binding remained. 8rexΔ animals exhibited no changes in X expression and lacked dosage-compensation mutant phenotypes. Hence, TAD boundaries are neither the cause nor the consequence of DCC-mediated gene repression. Abrogating TAD structure did, however, reduce thermotolerance, accelerate aging, and shorten lifespan, implicating chromosome architecture in stress responses and aging.

Progression of Meiosis Is Coordinated by the Level and Location of MAPK Activation Via OGR-2 in Caenorhabditis elegans

During meiosis, a series of evolutionarily conserved events allow for reductional chromosome division, which is required for sexual reproduction. Although individual meiotic processes have been extensively studied, we currently know far less about how meiosis is regulated and coordinated. In the Caenorhabditis elegans gonad, mitogen-activated protein kinase (MAPK) signaling drives oogenesis while undergoing spatial activation and deactivation waves. However, it is currently unclear how MAPK activation is governed and how it facilitates the progression of oogenesis. Here, we show that the oocyte and germline-related 2 (ogr-2) gene affects proper progression of oogenesis. Complete deletion of ogr-2 results in delayed meiotic entry and late spatial onset of double-strand break repair. Elevated levels of apoptosis are observed in this mutant, independent of the meiotic canonical checkpoints; however, they are dependent on the MAPK terminal member MPK-1/ERK. MPK-1 activation is elevated in diplotene in ogr-2 mutants and its aberrant spatial activation correlates with stages where meiotic progression defects are evident. Deletion of ogr-2 significantly reduces the expression of lip-1, a phosphatase reported to repress MPK-1, which is consistent with OGR-2 localization at chromatin in germ cells. We suggest that OGR-2 modulates the expression of lip-1 to promote the timely progression of meiosis through MPK-1 spatial deactivation.