Application of Transcriptional Gene Modules to Analysis of Caenorhabditis elegans' Gene Expression Data

The extracellular matrix integrates mitochondrial homeostasis

Cellular homeostasis is intricately influenced by stimuli from the microenvironment, including signaling molecules, metabolites, and pathogens. Functioning as a signaling hub within the cell, mitochondria integrate information from various intracellular compartments to regulate cellular signaling and metabolism. Multiple studies have shown that mitochondria may respond to various extracellular signaling events. However, it is less clear how changes in the extracellular matrix (ECM) can impact mitochondrial homeostasis to regulate animal physiology. We find that ECM remodeling alters mitochondrial homeostasis in an evolutionarily conserved manner. Mechanistically, ECM remodeling triggers a TGF-β response to induce mitochondrial fission and the unfolded protein response of the mitochondria (UPRMT). At the organismal level, ECM remodeling promotes defense of animals against pathogens through enhanced mitochondrial stress responses. We postulate that this ECM-mitochondria crosstalk represents an ancient immune pathway, which detects infection- or mechanical-stress-induced ECM damage, thereby initiating adaptive mitochondria-based immune and metabolic responses.

Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Spaceflight Induces Strength Decline in Caenorhabditis elegans

Background: Understanding and countering the well-established negative health consequences of spaceflight remains a primary challenge preventing safe deep space exploration. Targeted/personalized therapeutics are at the forefront of space medicine strategies, and cross-species molecular signatures now define the 'typical' spaceflight response. However, a lack of direct genotype-phenotype associations currently limits the robustness and, therefore, the therapeutic utility of putative mechanisms underpinning pathological changes in flight. Methods: We employed the worm Caenorhabditis elegans as a validated model of space biology, combined with 'NemaFlex-S' microfluidic devices for assessing animal strength production as one of the most reproducible physiological responses to spaceflight. Wild-type and dys-1 (BZ33) strains (a Duchenne muscular dystrophy (DMD) model for comparing predisposed muscle weak animals) were cultured on the International Space Station in chemically defined media before loading second-generation gravid adults into NemaFlex-S devices to assess individual animal strength. These same cultures were then frozen on orbit before returning to Earth for next-generation sequencing transcriptomic analysis. Results: Neuromuscular strength was lower in flight versus ground controls (16.6% decline, p < 0.05), with dys-1 significantly more (23% less strength, p < 0.01) affected than wild types. The transcriptional gene ontology signatures characterizing both strains of weaker animals in flight strongly corroborate previous results across species, enriched for upregulated stress response pathways and downregulated mitochondrial and cytoskeletal processes. Functional gene cluster analysis extended this to implicate decreased neuronal function, including abnormal calcium handling and acetylcholine signaling, in space-induced strength declines under the predicted control of UNC-89 and DAF-19 transcription factors. Finally, gene modules specifically altered in dys-1 animals in flight again cluster to neuronal/neuromuscular pathways, suggesting strength loss in DMD comprises a strong neuronal component that predisposes these animals to exacerbated strength loss in space. Conclusions: Highly reproducible gene signatures are strongly associated with space-induced neuromuscular strength loss across species and neuronal changes in calcium/acetylcholine signaling require further study. These results promote targeted medical efforts towards and provide an in vivo model for safely sending animals and people into deep space in the near future.

Individual cell types in C. elegans age differently and activate distinct cell-protective responses

Aging is characterized by a global decline in physiological function. However, by constructing a complete single-cell gene expression atlas, we find that Caenorhabditis elegans aging is not random in nature but instead is characterized by coordinated changes in functionally related metabolic, proteostasis, and stress-response genes in a cell-type-specific fashion, with downregulation of energy metabolism being the only nearly universal change. Similarly, the rates at which cells age differ significantly between cell types. In some cell types, aging is characterized by an increase in cell-to-cell variance, whereas in others, variance actually decreases. Remarkably, multiple resilience-enhancing transcription factors known to extend lifespan are activated across many cell types with age; we discovered new longevity candidates, such as GEI-3, among these. Together, our findings suggest that cells do not age passively but instead react strongly, and individualistically, to events that occur during aging. This atlas can be queried through a public interface.

Defining characteristics and conservation of poorly annotated genes in Caenorhabditis elegans using WormCat 2.0

Omics tools provide broad datasets for biological discovery. However, the computational tools for identifying important genes or pathways in RNA-seq, proteomics, or GWAS (Genome-Wide Association Study) data depend on Gene Ontogeny annotations and are biased toward well-described pathways. This limits their utility as poorly annotated genes, which could have novel functions, are often passed over. Recently, we developed an annotation and category enrichment tool for Caenorhabditis elegans genomic data, WormCat, which provides an intuitive visualization output. Unlike Gene Ontogeny-based enrichment tools, which exclude genes with no annotation information, WormCat 2.0 retains these genes as a special UNASSIGNED category. Here, we show that the UNASSIGNED gene category enrichment exhibits tissue-specific expression patterns and can include genes with biological functions identified in published datasets. Poorly annotated genes are often considered to be potentially species-specific and thus, of reduced interest to the biomedical community. Instead, we find that around 3% of the UNASSIGNED genes have human orthologs, including some linked to human diseases. These human orthologs themselves have little annotation information. A recently developed method that incorporates lineage relationships (abSENSE) indicates that the failure of BLAST to detect homology explains the apparent lineage specificity for many UNASSIGNED genes. This suggests that a larger subset could be related to human genes. WormCat provides an annotation strategy that allows the association of UNASSIGNED genes with specific phenotypes and known pathways. Building these associations in C. elegans, with its robust genetic tools, provides a path to further functional study and insight into these understudied genes.

Rotenone Modulates Caenorhabditis elegans Immunometabolism and Pathogen Susceptibility

Mitochondria are central players in host immunometabolism as they function not only as metabolic hubs but also as signaling platforms regulating innate immunity. Environmental exposures to mitochondrial toxicants occur widely and are increasingly frequent. Exposures to these mitotoxicants may pose a serious threat to organismal health and the onset of diseases by disrupting immunometabolic pathways. In this study, we investigated whether the Complex I inhibitor rotenone could alter C. elegans immunometabolism and disease susceptibility. C. elegans embryos were exposed to rotenone (0.5 µM) or DMSO (0.125%) until they reached the L4 larval stage. Inhibition of mitochondrial respiration by rotenone and disruption of mitochondrial metabolism were evidenced by rotenone-induced detrimental effects on mitochondrial efficiency and nematode growth and development. Next, through transcriptomic analysis, we investigated if this specific but mild mitochondrial stress that we detected would lead to the modulation of immunometabolic pathways. We found 179 differentially expressed genes (DEG), which were mostly involved in detoxification, energy metabolism, and pathogen defense. Interestingly, among the down-regulated DEG, most of the known genes were involved in immune defense, and most of these were identified as commonly upregulated during P. aeruginosa infection. Furthermore, rotenone increased susceptibility to the pathogen Pseudomonas aeruginosa (PA14). However, it increased resistance to Salmonella enterica (SL1344). To shed light on potential mechanisms related to these divergent effects on pathogen resistance, we assessed the activation of the mitochondrial unfolded protein response (UPRmt), a well-known immunometabolic pathway in C. elegans which links mitochondria and immunity and provides resistance to pathogen infection. The UPRmt pathway was activated in rotenone-treated nematodes further exposed for 24 h to the pathogenic bacteria P. aeruginosa and S. enterica or the common bacterial food source Escherichia coli (OP50). However, P. aeruginosa alone suppressed UPRmt activation and rotenone treatment rescued its activation only to the level of DMSO-exposed nematodes fed with E. coli. Module-weighted annotation bioinformatics analysis was also consistent with UPRmt activation in rotenone-exposed nematodes consistent with the UPR being involved in the increased resistance to S. enterica. Together, our results demonstrate that the mitotoxicant rotenone can disrupt C. elegans immunometabolism in ways likely protective against some pathogen species but sensitizing against others.

Transcriptomic profiling in canines and humans reveals cancer specific gene modules and biological mechanisms common to both species

Understanding relationships between spontaneous cancer in companion (pet) canines and humans can facilitate biomarker and drug development in both species. Towards this end we developed an experimental-bioinformatic protocol that analyzes canine transcriptomics data in the context of existing human data to evaluate comparative relevance of canine to human cancer. We used this protocol to characterize five canine cancers: melanoma, osteosarcoma, pulmonary carcinoma, B- and T-cell lymphoma, in 60 dogs. We applied an unsupervised, iterative clustering method that yielded five co-expression modules and found that each cancer exhibited a unique module expression profile. We constructed cancer models based on the co-expression modules and used the models to successfully classify the canine data. These canine-derived models also successfully classified human tumors representing the same cancers, indicating shared cancer biology between canines and humans. Annotation of the module genes identified cancer specific pathways relevant to cells-of-origin and tumor biology. For example, annotations associated with melanin production (PMEL, GPNMB, and BACE2), synthesis of bone material (COL5A2, COL6A3, and COL12A1), synthesis of pulmonary surfactant (CTSH, LPCAT1, and NAPSA), ribosomal proteins (RPL8, RPS7, and RPLP0), and epigenetic regulation (EDEM1, PTK2B, and JAK1) were unique to melanoma, osteosarcoma, pulmonary carcinoma, B- and T-cell lymphoma, respectively. In total, 152 biomarker candidates were selected from highly expressing modules for each cancer type. Many of these biomarker candidates are under-explored as drug discovery targets and warrant further study. The demonstrated transferability of classification models from canines to humans enforces the idea that tumor biology, biomarker targets, and associated therapeutics, discovered in canines, may translate to human medicine.

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.