A Caenorhabditis elegans protein with a PRDM9-like SET domain localizes to chromatin-associated foci and promotes spermatocyte gene expression, sperm production and fertility

Two H3K23 histone methyltransferases, SET-32 and SET-21, function synergistically to promote nuclear RNAi-mediated transgenerational epigenetic inheritance in Caenorhabditis elegans

Nuclear RNAi in Caenorhabditis elegans induces a set of transgenerationally heritable marks of H3K9me3, H3K23me3, and H3K27me3 at the target genes. The function of H3K23me3 in the nuclear RNAi pathway is largely unknown due to the limited knowledge of H3K23 histone methyltransferase (HMT). In this study we identified SET-21 as a novel H3K23 HMT. By taking combined genetic, biochemical, imaging, and genomic approaches, we found that SET-21 functions synergistically with a previously reported H3K23 HMT SET-32 to deposit H3K23me3 at the native targets of germline nuclear RNAi. We identified a subset of native nuclear RNAi targets that are transcriptionally activated in the set-21;set-32 double mutant. SET-21 and SET-32 are also required for robust transgenerational gene silencing induced by exogenous dsRNA. The set-21;set-32 double mutant strain exhibits an enhanced temperature-sensitive mortal germline phenotype compared to the set-32 single mutant, while the set-21 single mutant animals are fertile. We also found that HRDE-1 and SET-32 are required for cosuppression, a transgene-induced gene silencing phenomenon, in C. elegans germline. Together, these results support a model in which H3K23 HMTs SET-21 and SET-32 function cooperatively as germline nuclear RNAi factors and promote the germline immortality under the heat stress.

Transcriptional responses to in vitro macrocyclic lactone exposure in Toxocara canis larvae using RNA-seq

Toxocara canis, the causative agent of zoonotic toxocariasis in humans, is a parasitic roundworm of canids with a complex lifecycle. While macrocyclic lactones (MLs) are successful at treating adult T. canis infections when used at FDA-approved doses in dogs, they fail to kill somatic third-stage larvae. In this study, we profiled the transcriptome of third-stage larvae derived from larvated eggs and treated in vitro with 10 μM of the MLs - ivermectin and moxidectin with Illumina sequencing. We analyzed transcriptional changes in comparison with untreated control larvae. In ivermectin-treated larvae, we identified 608 differentially expressed genes (DEGs), of which 453 were upregulated and 155 were downregulated. In moxidectin-treated larvae, we identified 1,413 DEGs, of which 902 were upregulated and 511 were downregulated. Notably, many DEGs were involved in critical biological processes and pathways including transcriptional regulation, energy metabolism, neuronal structure and function, physiological processes such as reproduction, excretory/secretory molecule production, host-parasite response mechanisms, and parasite elimination. We also assessed the expression of known ML targets and transporters, including glutamate-gated chloride channels (GluCls), and ATP-binding cassette (ABC) transporters, subfamily B, with a particular focus on P-glycoproteins (P-gps). We present gene names for previously uncharacterized T. canis GluCl genes using phylogenetic analysis of nematode orthologs to provide uniform gene nomenclature. Our study revealed that the expression of Tca-glc-3 and six ABCB genes, particularly four P-gps, were significantly altered in response to ML treatment. Compared to controls, Tca-glc-3, Tca-Pgp-11.2, and Tca-Pgp-13.2 were downregulated in ivermectin-treated larvae, while Tca-abcb1, Tca-abcb7, Tca-Pgp-11.2, and Tca-Pgp-13.2 were downregulated in moxidectin-treated larvae. Conversely, Tca-abcb9.1 and Tca-Pgp-11.3 were upregulated in moxidectin-treated larvae. These findings suggest that MLs broadly impact transcriptional regulation in T. canis larvae.

Two H3K23 histone methyltransferases, SET-32 and SET-21, function synergistically to promote nuclear RNAi-mediated transgenerational epigenetic inheritance in Caenorhabditis elegans

Nuclear RNAi in C. elegans induces a set of transgenerationally heritable marks of H3K9me3, H3K23me3, and H3K27me3 at the target genes. The function of H3K23me3 in the nuclear RNAi pathway is largely unknown due to the limited knowledge of H3K23 histone methyltransferase (HMT). In this study we identified SET-21 as a novel H3K23 HMT. By taking combined genetic, biochemical, imaging, and genomic approaches, we found that SET-21 functions synergistically with a previously reported H3K23 HMT SET-32 to deposit H3K23me3 at the native targets of germline nuclear RNAi. We identified a subset of native nuclear RNAi targets that are transcriptionally activated in the set-21;set-32 double mutant. SET-21 and SET-32 are also required for robust transgenerational gene silencing induced by exogenous dsRNA. The set-21;set-32 double mutant strain exhibits an enhanced temperature-sensitive mortal germline phenotype compared to the set-32 single mutant, while the set-21 single mutant animals are fertile. We also found that HRDE-1 and SET-32 are required for cosuppression, a transgene-induced gene silencing phenomenon, in C. elegans germline. Together, these results support a model in which H3K23 HMTs SET-21 and SET-32 function cooperatively to ensure the robustness of germline nuclear RNAi and promotes the germline immortality under the heat stress.

RNAi-dependent expression of sperm genes in ADL chemosensory neurons is required for olfactory responses in Caenorhabditis elegans

Environmental conditions experienced early in the life of an animal can result in gene expression changes later in its life history. We have previously shown that C. elegans animals that experienced the developmentally arrested and stress resistant dauer stage (postdauers) retain a cellular memory of early-life stress that manifests during adulthood as genome-wide changes in gene expression, chromatin states, and altered life history traits. One consequence of developmental reprogramming in C. elegans postdauer adults is the downregulation of osm-9 TRPV channel gene expression in the ADL chemosensory neurons resulting in reduced avoidance to a pheromone component, ascr#3. This altered response to ascr#3 requires the principal effector of the somatic nuclear RNAi pathway, the Argonaute (AGO) NRDE-3. To investigate the role of the somatic nuclear RNAi pathway in regulating the developmental reprogramming of ADL due to early-life stress, we profiled the mRNA transcriptome of control and postdauer ADL in wild-type and nrde-3 mutant adults. We found 711 differentially expressed (DE) genes between control and postdauer ADL neurons, 90% of which are dependent upon NRDE-3. Additionally, we identified a conserved sequence that is enriched in the upstream regulatory sequences of the NRDE-3-dependent differentially expressed genes. Surprisingly, 214 of the ADL DE genes are considered "germline-expressed", including 21 genes encoding the Major Sperm Proteins and two genes encoding the sperm-specific PP1 phosphatases, GSP-3 and GSP-4. Loss of function mutations in gsp-3 resulted in both aberrant avoidance and attraction behaviors. We also show that an AGO pseudogene, Y49F6A.1 (wago-11), is expressed in ADL and is required for ascr#3 avoidance. Overall, our results suggest that small RNAs and reproductive genes program the ADL mRNA transcriptome during their developmental history and highlight a nexus between neuronal and reproductive networks in calibrating animal neuroplasticity.

Characterization of the Pristionchus pacificus "epigenetic toolkit" reveals the evolutionary loss of the histone methyltransferase complex PRC2

Comparative approaches have revealed both divergent and convergent paths to achieving shared developmental outcomes. Thus, only through assembling multiple case studies can we understand biological principles. Yet, despite appreciating the conservation-or lack thereof-of developmental networks, the conservation of epigenetic mechanisms regulating these networks is poorly understood. The nematode Pristionchus pacificus has emerged as a model system of plasticity and epigenetic regulation as it exhibits a bacterivorous or omnivorous morph depending on its environment. Here, we determined the "epigenetic toolkit" available to P. pacificus as a resource for future functional work on plasticity, and as a comparison with Caenorhabditis elegans to investigate the conservation of epigenetic mechanisms. Broadly, we observed a similar cast of genes with putative epigenetic function between C. elegans and P. pacificus. However, we also found striking differences. Most notably, the histone methyltransferase complex PRC2 appears to be missing in P. pacificus. We described the deletion/pseudogenization of the PRC2 genes mes-2 and mes-6 and concluded that both were lost in the last common ancestor of P. pacificus and a related species P. arcanus. Interestingly, we observed the enzymatic product of PRC2 (H3K27me3) by mass spectrometry and immunofluorescence, suggesting that a currently unknown methyltransferase has been co-opted for heterochromatin silencing. Altogether, we have provided an inventory of epigenetic genes in P. pacificus to compare with C. elegans. This inventory will enable reverse-genetic experiments related to plasticity and has revealed the first loss of PRC2 in a multicellular organism.

Bacterial N4-methylcytosine as an epigenetic mark in eukaryotic DNA

DNA modifications are used to regulate gene expression and defend against invading genetic elements. In eukaryotes, modifications predominantly involve C5-methylcytosine (5mC) and occasionally N6-methyladenine (6mA), while bacteria frequently use N4-methylcytosine (4mC) in addition to 5mC and 6mA. Here we report that 4mC can serve as an epigenetic mark in eukaryotes. Bdelloid rotifers, tiny freshwater invertebrates with transposon-poor genomes rich in foreign genes, lack canonical eukaryotic C5-methyltransferases for 5mC addition, but encode an amino-methyltransferase, N4CMT, captured from bacteria >60 Mya. N4CMT deposits 4mC at active transposons and certain tandem repeats, and fusion to a chromodomain shapes its “histone-read-DNA-write” architecture recognizing silent chromatin marks. Furthermore, amplification of SETDB1 H3K9me3 histone methyltransferases yields variants preferentially binding 4mC-DNA, suggesting “DNA-read-histone-write” partnership to maintain chromatin-based silencing. Our results show how non-native DNA methyl groups can reshape epigenetic systems to silence transposons and demonstrate the potential of horizontal gene transfer to drive regulatory innovation in eukaryotes.

Eukaryotic DNA can be methylated as 5-methylcytosine and N6-methyladenine, but whether other forms of DNA methylation occur has been controversial. Here the authors show that a bacterial DNA methyltransferase was acquired >60 Mya in bdelloid rotifers that catalyzes N4-methylcytosine addition and is involved in suppression of transposon proliferation.

A systems biology analysis of reproductive toxicity effects induced by multigenerational exposure to ionizing radiation in C. elegans

Understanding the effects of chronic exposure to pollutants over generations is of primary importance for the protection of humans and the environment; however, to date, knowledge on the molecular mechanisms underlying multigenerational adverse effects is scarce. We employed a systems biology approach to analyze effects of chronic exposure to gamma radiation at molecular, tissue and individual levels in the nematode Caenorhabditis elegans. Our data show a decrease of 23% in the number of offspring on the first generation F0 and more than 40% in subsequent generations F1, F2 and F3. To unveil the impact on the germline, an in-depth analysis of reproductive processes involved in gametes formation was performed for all four generations. We measured a decrease in the number of mitotic germ cells accompanied by increased cell-cycle arrest in the distal part of the gonad. Further impact on the germline was manifested by decreased sperm quantity and quality. In order to obtain insight in the molecular mechanisms leading to decreased fecundity, gene expression was investigated via whole genome RNA sequencing. The transcriptomic analysis revealed modulation of transcription factors, as well as genes involved in stress response, unfolded protein response, lipid metabolism and reproduction. Furthermore, a drastic increase in the number of differentially expressed genes involved in defense response was measured in the last two generations, suggesting a cumulative stress effect of ionizing radiation exposure. Transcription factor binding site enrichment analysis and the use of transgenic strain identified daf-16/FOXO as a master regulator of genes differentially expressed in response to radiation. The presented data provide new knowledge with respect to the molecular mechanisms involved in reproductive toxic effects and accumulated stress resulting from multigenerational exposure to ionizing radiation.

hsp-90 and unc-45 depletion induce characteristic transcriptional signatures in coexpression cliques of C. elegans

Nematode development is characterized by progression through several larval stages. Thousands of genes were found in large scale RNAi-experiments to block this development at certain steps, two of which target the molecular chaperone HSP-90 and its cofactor UNC-45. Aiming to define the cause of arrest, we here investigate the status of nematodes after treatment with RNAi against hsp-90 and unc-45 by employing an in-depth transcriptional analysis of the arrested larvae. To identify misregulated transcriptional units, we calculate and validate genome-wide coexpression cliques covering the entire nematode genome. We define 307 coexpression cliques and more than half of these can be related to organismal functions by GO-term enrichment, phenotype enrichment or tissue enrichment analysis. Importantly, hsp-90 and unc-45 RNAi induce or repress many of these cliques in a coordinated manner, and then several specifically regulated cliques are observed. To map the developmental state of the arrested nematodes we define the expression behaviour of each of the cliques during development from embryo to adult nematode. hsp-90 RNAi can be seen to arrest development close to the L4 larval stage with further deviations in daf-16 regulated genes. unc-45 RNAi instead leads to arrested development at young adult stage prior to the programmatic downregulation of sperm-cell specific genes. In both cases processes can be defined to be misregulated upon depletion of the respective chaperone. With most of the defined gene cliques showing concerted behaviour at some stage of development from embryo to late adult, the “clique map” together with the clique-specific GO-terms, tissue and phenotype assignments will be a valuable tool in understanding concerted responses on the genome-wide level in Caenorhabditis elegans.

Investigating autism associated genes in C. elegans reveals candidates with a role in social behaviour

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterised by a triad of behavioural impairments and includes disruption in social behaviour. ASD has a clear genetic underpinning and hundreds of genes are implicated in its aetiology. However, how single penetrant genes disrupt activity of neural circuits which lead to affected behaviours is only beginning to be understood and less is known about how low penetrant genes interact to disrupt emergent behaviours. Investigations are well served by experimental approaches that allow tractable investigation of the underpinning genetic basis of circuits that control behaviours that operate in the biological domains that are neuro-atypical in autism. The model organism C. elegans provides an experimental platform to investigate the effect of genetic mutations on behavioural outputs including those that impact social biology. Here we use progeny-derived social cues that modulate C. elegans food leaving to assay genetic determinants of social behaviour. We used the SAFRI Gene database to identify C. elegans orthologues of human ASD associated genes. We identified a number of mutants that displayed selective deficits in response to progeny. The genetic determinants of this complex social behaviour highlight the important contribution of synaptopathy and implicates genes within cell signalling, epigenetics and phospholipid metabolism functional domains. The approach overlaps with a growing number of studies that investigate potential molecular determinants of autism in C. elegans. However, our use of a complex, sensory integrative, emergent behaviour provides routes to enrich new or underexplored biology with the identification of novel candidate genes with a definable role in social behaviour.

A nervous system-specific subnuclear organelle in Caenorhabditis elegans

We describe here phase-separated subnuclear organelles in the nematode Caenorhabditis elegans, which we term NUN (NUclear Nervous system-specific) bodies. Unlike other previously described subnuclear organelles, NUN bodies are highly cell type specific. In fully mature animals, 4-10 NUN bodies are observed exclusively in the nucleus of neuronal, glial and neuron-like cells, but not in other somatic cell types. Based on co-localization and genetic loss of function studies, NUN bodies are not related to other previously described subnuclear organelles, such as nucleoli, splicing speckles, paraspeckles, Polycomb bodies, promyelocytic leukemia bodies, gems, stress-induced nuclear bodies, or clastosomes. NUN bodies form immediately after cell cycle exit, before other signs of overt neuronal differentiation and are unaffected by the genetic elimination of transcription factors that control many other aspects of neuronal identity. In one unusual neuron class, the canal-associated neurons, NUN bodies remodel during larval development, and this remodeling depends on the Prd-type homeobox gene ceh-10. In conclusion, we have characterized here a novel subnuclear organelle whose cell type specificity poses the intriguing question of what biochemical process in the nucleus makes all nervous system-associated cells different from cells outside the nervous system.

The KMT2A gene: mRNA differential expression in the ovary and a novel 13-nt nucleotide sequence variant associated with litter size in cashmere goats

A genome-wide association study had shown that lysine methyltransferase 2A (KMT2A), which encodes the histone 3 lysine 4 methyltransferase and reportedly can regulate gametogenesis, steroidogenesis, and development as well as other biological processes, is a potential candidate gene influencing litter size in the dairy goat, suggesting its key function in animal reproduction. Here, we aimed to explore the genetic effects of the KMT2A gene on litter size in females of the Chinese indigenous cashmere goat, using a large sample size (n > 1,000), based on their levels of RNA transcription and DNA variation. First, mRNA expression levels of this gene in ovarian tissues between the low-prolific group (first-born litter size = 1) and high-prolific group (first-born litter size ≥2) were significantly different, revealing the potential functioning of KMT2A in goat prolific. Moreover, a novel 13-nt nucleotide sequence variant was identified in Shaanbei white cashmere goats (n = 1,616). In accordance with the independent chi-square (χ²) analysis, the distribution of genotypes (P = 2.57 × 10^-9) and allelotypes (P = 3.00 × 10^-7) between the low- and high-prolific groups differed significantly, indicating the 13-nt mutation was associated with litter size. Further analysis showed that the insertion/insertion (II) genotype was significantly different with insertion/deletion (ID) (P = 1.76 × 10^-9) and deletion/deletion (DD) (P = 7.00 × 10^-6), with goats having the DD genotype producing an average litter size larger than the other genotypes. Taken together, these findings suggest KMT2A can serve as a candidate gene for breeding goats, which may have implications for improving the future development of the goat industry.

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 Effects of Age, Cigarette Smoking, Sex, and Race on the Qualitative Characteristics of Lung Transcriptome

The within-sample relative expression orderings (REOs) of genes, which are stable qualitative transcriptional characteristics, can provide abundant information for a disease. Methods based on REO comparisons have been proposed for identifying differentially expressed genes (DEGs) at the individual level and for detecting disease-associated genes based on one-phenotype disease data by reusing data of normal samples from other sources. Here, we evaluated the effects of common potential confounding factors, including age, cigarette smoking, sex, and race, on the REOs of gene pairs within normal lung tissues transcriptome. Our results showed that age has little effect on REOs within lung tissues. We found that about 0.23% of the significantly stable REOs of gene pairs in nonsmokers' lung tissues are reversed in smokers' lung tissues, introduced by 344 DEGs between the two groups of samples (RankCompV2, FDR <0.05), which are enriched in metabolism of xenobiotics by cytochrome P450, glutathione metabolism, and other pathways (hypergeometric test, FDR <0.05). Comparison between the normal lung tissue samples of males and females revealed fewer reversal REOs introduced by 24 DEGs between the sex groups, among which 19 DEGs are located on sex chromosomes and 5 DEGs involving in spermatogenesis and regulation of oocyte are located on autosomes. Between the normal lung tissue samples of white and black people, we identified 22 DEGs (RankCompV2, FDR <0.05) which introduced a few reversal REOs between the two races. In summary, the REO-based study should take into account the confounding factors of cigarette smoking, sex, and race.

Chromatin Modifiers SET-25 and SET-32 Are Required for Establishment but Not Long-Term Maintenance of Transgenerational Epigenetic Inheritance

Some epigenetic modifications are inherited from one generation to the next, providing a potential mechanism for the inheritance of environmentally acquired traits. Transgenerational inheritance of RNAi phenotypes in Caenorhabditis elegans provides an excellent model to study this phenomenon, and although studies have implicated both chromatin modifications and small RNA pathways in heritable silencing, their relative contributions remain unclear. Here, we demonstrate that the putative histone methyltransferases SET-25 and SET-32 are required for establishment of a transgenerational silencing signal but not for long-term maintenance of this signal between subsequent generations, suggesting that transgenerational epigenetic inheritance is a multi-step process with distinct genetic requirements for establishment and maintenance of heritable silencing. Furthermore, small RNA sequencing reveals that the abundance of secondary siRNAs (thought to be the effector molecules of heritable silencing) does not correlate with silencing phenotypes. Together, our results suggest that the current mechanistic models of epigenetic inheritance are incomplete.

Gamma radiation induces life stage-dependent reprotoxicity in Caenorhabditis elegans via impairment of spermatogenesis

The current study investigated life stage, tissue and cell dependent sensitivity to ionizing radiation of the nematode Caenorhabditis elegans. Results showed that irradiation of post mitotic L4 stage larvae induced no significant effects with respect to mortality, morbidity or reproduction at either acute dose ≤6 Gy (1500 mGy·h^-1) or chronic exposure ≤15 Gy (≤100 mGy·h^-1). In contrast, chronic exposure from the embryo to the L4-young adult stage caused a dose and dose-rate dependent reprotoxicity with 43% reduction in total brood size at 6.7 Gy (108 mGy·h^-1). Systematic irradiation of the different developmental stages showed that the most sensitive life stage was L1 to young L4. Exposure during these stages was associated with dose-rate dependent genotoxic effects, resulting in a 1.8 to 2 fold increase in germ cell apoptosis in larvae subjected to 40 or 100 mGy·h^-1, respectively. This was accompanied by a dose-rate dependent reduction in the number of spermatids, which was positively correlated to the reprotoxic effect (0.99, PCC). RNAseq analysis of nematodes irradiated from L1 to L4 stage revealed a significant enrichment of differentially expressed genes related to both male and hermaphrodite reproductive processes. Gene network analysis revealed effects related to down-regulation of genes required for spindle formation and sperm meiosis/maturation, including smz-1, smz-2 and htas-1. Furthermore, the expression of a subset of 28 set-17 regulated Major Sperm Proteins (MSP) required for spermatid production was correlated (R² 0.80) to the reduction in reproduction and the number of spermatids. Collectively these observations corroborate the impairment of spermatogenesis as the major cause of gamma radiation induced life-stage dependent reprotoxic effect. Furthermore, the progeny of irradiated nematodes showed significant embryonal DNA damage that was associated with persistent effect on somatic growth. Unexpectedly, these nematodes maintained much of their reproductive capacity in spite of the reduced growth.