Systematic Identification of Housekeeping Genes Possibly Used as References in Caenorhabditis elegans by Large-Scale Data Integration

The Molecular Basis of Differentiation Wave Activity in Embryogenesis

As development varies greatly across the tree of life, it may seem difficult to suggest a model that proposes a single mechanism for understanding collective cell behaviors and the coordination of tissue formation. Here we propose a mechanism called differentiation waves, which unify many disparate results involving developmental systems from across the tree of life. We demonstrate how a relatively simple model of differentiation proceeds not from function-related molecular mechanisms, but from so-called differentiation waves. A phenotypic model of differentiation waves is introduced, and its relation to molecular mechanisms is proposed. These waves contribute to a differentiation tree, which is an alternate way of viewing cell lineage and local action of the molecular factors. We construct a model of differentiation wave-related molecular mechanisms (genome, epigenome, and proteome) based on bioinformatic data from the nematode Caenorhabditis elegans. To validate this approach across different modes of development, we evaluate protein expression across different types of development by comparing Caenorhabditis elegans with several model organisms: fruit flies (Drosophila melanogaster), yeast (Saccharomyces cerevisiae), and mouse (Mus musculus). Inspired by gene regulatory networks, two Models of Interactive Contributions (fully-connected MICs and ordered MICs) are used to suggest potential genomic contributions to differentiation wave-related proteins. This, in turn, provides a framework for understanding differentiation and development.

Selection and Validation of Reference Genes for Gene Expression in Bactericera gobica Loginova under Different Insecticide Stresses

Insecticide resistance has long been a problem in crop pest control. Bactericera gobica is a major pest on the well-known medicinal plants Lycium barbarum L. Investigating insecticide resistance mechanisms of B. gobica will help to identify pesticide reduction strategies to control the pest. Gene expression normalization by RT-qPCR requires the selection and validation of appropriate reference genes (RGs). Here, 15 candidate RGs were selected from transcriptome data of B. gobica. Their expression stability was evaluated with five algorithms (Delta Ct, GeNorm, Normfinder, BestKeeper and RefFinder) for sample types differing in response to five insecticide stresses and in four other experimental conditions. Our results indicated that the RGs RPL10 + RPS15 for Imidacloprid and Abamectin; RPL10 + AK for Thiamethoxam; RPL32 + RPL10 for λ-cyhalothrin; RPL10 + RPL8 for Matrine; and EF2 + RPL32 under different insecticide stresses were the most suitable RGs for RT-qPCR normalization. EF1α + RPL8, EF1α + β-actin, β-actin + EF2 and β-actin + RPS15 were the optimal combination of RGs under odor stimulation, temperature, developmental stages and both sexes, respectively. Overall, EF2 and RPL8 were the two most stable RGs in all conditions, while α-TUB and RPL32 were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target cytochrome P450 CYP6a1 gene between adult and nymph stages and under imidacloprid stress. The results of CYP6a1 expression were consistent with transcriptome data. This study is the first research on the most stable RG selection in B. gobica nymphs exposed to different insecticides, which will contribute to further research on insecticide resistance mechanisms in B. gobica.

Integration of RNA-seq and ATAC-seq analyzes the effect of low dose neutron-γ radiation on gene expression of lymphocytes from oilfield logging workers

Objective: Although radiation workers are exposed to much lower doses of neutron-γ rays than those suffered in nuclear explosions and accidents, it does not mean that their health is not affected by radiation. Lower doses of radiation do not always cause morphological aberrations in chromosomes, so more sophisticated tests must be sought to specific alterations in the exposed cells. Our goal was to characterize the specific gene expression in lymphocytes from logging workers who were continuously exposed to low doses of neutron-γ radiation. We hypothesized that the combination of cell type-specific transcriptomes and open chromatin profiles would identify lymphocyte-specific gene alterations induced by long-term radiation with low-dose neutron-γ-rays and discover new regulatory pathways and transcriptional regulatory elements. Methods: Lymphocytes were extracted from workers who have been occupationally exposed to neutron-γ and workers unexposed to radiation in the same company. mRNA-seq and ATAC-seq (Assay for Transposase-Accessible Chromatin with high-throughput sequencing) were performed, followed integrative analysis to identify specific gene regulatory regions induced by neutron-γ radiation. A qPCR assay was then performed to verify the downregulation of RNA coding for ribosomal proteins and flow cytometry was used to detect ribosomal protein expression and cell cycle alterations. Results: We identified transcripts that were specifically induced by neutron-γ radiation and discovered differential open chromatin regions that correlated with these gene activation patterns. Notably, we observed a downward trend in the expression of both differentially expressed genes and open chromatin peaks. Our most significant finding was that the differential peak upregulated in ATAC-seq, while the differential gene was downregulated in the ribosome pathway. We confirmed that neutron-γ radiation leads to transcriptional inhibition by analyzing the most enriched promoters, examining RPS18 and RPS27A expression by qPCR, and analyzing protein-protein interactions of the differential genes. Ribosomal protein expression and cell cycle were also affected by neutron-γ as detected by flow cytometry. Conclusion: We have comprehensively analyzed the genetic landscape of human lymphocytes based on chromatin accessibility and transcript levels, enabling the identification of novel neutron-γ induced signature genes not previously known. By comparing fine-mapping of open chromatin and RNA reads, we have determined that neutron-γ specifically leads to downregulation of genes in the ribosome pathway, with pseudogenes potentially playing a crucial role.

Eicosapentaenoic acid influences the pathogenesis of Candida albicans in Caenorhabditis elegans via inhibition of hyphal formation and stimulation of the host immune response

The intake of omega-3 polyunsaturated fatty acids, including eicosapentaenoic acid (EPA), is associated with health benefits due to its anti-inflammatory properties. This fatty acid also exhibits antifungal properties in vitro. In order to determine if this antifungal property is valid in vivo, we examined how EPA affects Candida albicans pathogenesis in the Caenorhabditis elegans infection model, an alternative to mammalian host models. The nematodes were supplemented with EPA prior to infection, and the influence of EPA on C. elegans lipid metabolism, survival and immune response was studied. In addition, the influence of EPA on hyphal formation in C. albicans was investigated. It was discovered that EPA supplementation changed the lipid composition, but not the unsaturation index of C. elegans by regulating genes involved in fatty acid and eicosanoid production. EPA supplementation also delayed killing of C. elegans by C. albicans due to the inhibition of hyphal formation in vivo, via the action of the eicosanoid metabolite of EPA, 17,18-epoxyeicosatetraenoic acid. Moreover, EPA supplementation also caused differential expression of biofilm-related gene expression in C. albicans and stimulated the immune response of C. elegans. This provides a link between EPA and host susceptibility to microbial infection in this model.

Whole-body gene expression atlas of an adult metazoan

Gene activity defines cell identity, drives intercellular communication, and underlies the functioning of multicellular organisms. We present the single-cell resolution atlas of gene activity of a fertile adult metazoan: Caenorhabditis elegans. This compendium comprises 180 distinct cell types and 19,657 expressed genes. We predict 7541 transcription factor expression profile associations likely responsible for defining cellular identity. We predict thousands of intercellular interactions across the C. elegans body and the ligand-receptor pairs that mediate them, some of which we experimentally validate. We identify 172 genes that show consistent expression across cell types, are involved in basic and essential functions, and are conserved across phyla; therefore, we present them as experimentally validated housekeeping genes. We developed the WormSeq application to explore these data. In addition to the integrated gene-to-systems biology, we present genome-scale single-cell resolution testable hypotheses that we anticipate will advance our understanding of the molecular mechanisms, underlying the functioning of a multicellular organism and the perturbations that lead to its malfunction.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

Transcriptome-Based Selection and Validation of Reference Genes for Gene Expression in Goji Fruit Fly (Neoceratitis asiatica Becker) under Developmental Stages and Five Abiotic Stresses

Goji fruit fly, Neoceratitis asiatica, is a major pest on the well-known medicinal plant Lycium barbarum. Dissecting molecular mechanisms of infestation and host selection of N. asiatica will contribute to the determination of best management practices for pest fly control. Gene expression normalization by Real-time quantitative PCR (qPCR) requires the selection and validation of appropriate reference genes (RGs). Hence, 15 candidate RGs were selected from transcriptome data of N. asiatica. Their expression stability was evaluated with five algorithms (∆Ct, Normfinder, GeNorm, BestKeeper, and RefFinder) for sample types differing in the developmental stage, sex, tissue type, and in response to five different abiotic stresses. Our results indicated that the RGs β-Actin + GST for sex, RPL32 + EF1α for tissue type, RPS13+ EF1α for developmental stages along with odor stimulation, color induction, and starvation-refeeding stresses, EF1α + GAPDH under insecticide stress, RPS13 + RPS18 under temperature stress, respectively, were selected as the most suitable RGs for qPCR normalization. Overall, RPS18 and EF1α were the two most stable RGs in all conditions, while RPS15 and EF1β were the least stable RGs. The corresponding suitable RGs and one unstable RG were used to normalize a target odorant-binding protein OBP56a gene in male and female antennae, different tissues, and under odor stimulation. The results of OBP56a expression were consistent with transcriptome data. Our study is the first research on the most stable RGs selection in N. asiatica, which will facilitate further studies on the mechanisms of host selection and insecticide resistance in N. asiatica.

Multiple reaction monitoring mass spectrometry for the discovery of environmentally modulated proteins in an aquatic invertebrate sentinel species, Gammarus fossarum

Multiple reaction monitoring (MRM) mass spectrometry is emerging as a relevant tool for measuring customized molecular markers in freshwater sentinel species. While this technique is typically used for the validation of protein molecular markers preselected from shotgun experiments, recent gains of MRM multiplexing capacity offer new possibilities to conduct large-scale screening of animal proteomes. By combining the strength of active biomonitoring strategies and MRM technologies, this study aims to propose a new strategy for the discovery of candidate proteins that respond to environmental variability. For this purpose, 249 peptides derived from 147 proteins were monitored by MRM in 273 male gammarids caged in 56 environmental sites, representative of the diversity of French water bodies. A methodology is here proposed to identify a set of customized housekeeping peptides (HKPs) used to correct analytical batch effects and allow proper comparison of peptide levels in gammarids. A comparative analysis performed on HKPs-normalized data resulted in the identification of peptides highly modulated in the environment and derived from proteins likely involved in the environmental stress response. Overall, this study proposes a breakthrough approach to screen and identify potential proteins responding to relevant environmental conditions in sentinel species.

TENT5 cytoplasmic noncanonical poly(A) polymerases regulate the innate immune response in animals

Innate immunity is the first line of host defense against pathogens. Here, through global transcriptome and proteome analyses, we uncover that newly described cytoplasmic poly(A) polymerase TENT-5 (terminal nucleotidyltransferase 5) enhances the expression of secreted innate immunity effector proteins in Caenorhabditis elegans. Direct RNA sequencing revealed that multiple mRNAs with signal peptide-encoding sequences have shorter poly(A) tails in tent-5-deficient worms. Those mRNAs are translated at the endoplasmic reticulum where a fraction of TENT-5 is present, implying that they represent its direct substrates. Loss of tent-5 makes worms more susceptible to bacterial infection. Notably, the role of TENT-5 in innate immunity is evolutionarily conserved. Its orthologs, TENT5A and TENT5C, are expressed in macrophages and induced during their activation. Analysis of macrophages devoid of TENT5A/C revealed their role in the regulation of secreted proteins involved in defense response. In summary, our study reveals cytoplasmic polyadenylation to be a previously unknown component of the posttranscriptional regulation of innate immunity in animals.

Cutting-Edge Methods for Better Understanding Cells

Cells are microscopic yet fundamental elements of life [...].

An Integrated Bioinformatics Approach to Identify Network-Derived Hub Genes in Starving Zebrafish

The present study was aimed at identifying causative hub genes within modules formed by co-expression and protein-protein interaction (PPI) networks, followed by Bayesian network (BN) construction in the liver transcriptome of starved zebrafish. To this end, the GSE11107 and GSE112272 datasets from the GEO databases were downloaded and meta-analyzed using the MetaDE package, an add-on R package. Differentially expressed genes (DEGs) were identified based upon expression intensity N(µ = 0.2, σ² = 0.4). Reconstruction of BNs was performed by the bnlearn R package on genes within modules using STRINGdb and CEMiTool. ndufs5 (shared among PPI, BN and COEX), rps26, rpl10, sdhc (shared between PPI and BN), ndufa6, ndufa10, ndufb8 (shared between PPI and COEX), skp1, atp5h, ndufb10, rpl5b, zgc:193613, zgc:123327, zgc:123178, wu:fc58f10, zgc:111986, wu:fc37b12, taldo1, wu:fb62f08, zgc:64133 and acp5a (shared between COEX and BN) were identified as causative hub genes affecting gene expression in the liver of starving zebrafish. Future work will shed light on using integrative analyses of miRNA and DNA microarrays simultaneously, and performing in silico and experimental validation of these hub-causative (CST) genes affecting starvation in zebrafish.

Identifying optimal reference genes for gene expression studies in Eurasian spruce bark beetle, Ips typographus (Coleoptera: Curculionidae: Scolytinae)

Eurasian spruce bark beetle (Ips typographus [L.]) causes substantial damage to spruce forests worldwide. Undoubtedly, more aggressive measures are necessary to restrict the enduring loss. Finishing genome sequencing is a landmark achievement for deploying molecular techniques (i.e., RNA interference) to manage this pest. Gene expression studies assist in understanding insect physiology and deployment of molecular approaches for pest management. RT-qPCR is a valuable technique for such studies. However, accuracy and reliability depend on suitable reference genes. With the genome sequence available and the growing requirement of molecular tools for aggressive forest pest management, it is crucial to find suitable reference genes in Ips typographus under different experimental conditions. Hence, we evaluated the stability of twelve candidate reference genes under diverse experimental conditions such as biotic (developmental, sex and tissues) and abiotic factors (i.e., temperature and juvenile hormone treatment) to identify the reference genes. Our results revealed that ribosomal protein 3a (RPS3-a) was the best reference gene across all the experimental conditions, with minor exceptions. However, the stability of the reference gene can differ based on experiments. Nevertheless, present study provides a comprehensive list of reference genes under different experimental conditions for Ips typographus and contributes to "future genomic and functional genomic research".

K. S, T. S, Sachdeva K, Bathula C, K. V, K. S. N, Damodar S, Dhar SK, Das M. Relative quantification of BCL2 mRNA for diagnostic usage needs stable uncontrolled genes as reference

Dysregulation of BCL2 is a pathophysiology observed in haematological malignancies. For implementation of available treatment-options it is preferred to know the relative quantification of BCL2 mRNA with appropriate reference genes. For the choice of reference genes-(i) Reference Genes were selected by assessing variation of >60,000 genes from 4 RNA-seq datasets of haematological malignancies followed by filtering based on their GO biological process annotations and proximity of their chromosomal locations to known disease translocations. Selected genes were experimentally validated across various haematological malignancy samples followed by stability comparison using geNorm, NormFinder, BestKeeper and RefFinder. (ii) 43 commonly used Reference Genes were obtained from literature through extensive systematic review. Levels of BCL2 mRNA was assessed by qPCR normalized either by novel reference genes from this study or GAPDH, the most cited reference gene in literature and compared. The analysis showed PTCD2, PPP1R3B and FBXW9 to be the most unregulated genes across lymph-nodes, bone marrow and PBMC samples unlike the Reference Genes used in literature. BCL2 mRNA level shows a consistent higher expression in haematological malignancy patients when normalized by these novel Reference Genes as opposed to GAPDH, the most cited Reference Gene. These reference genes should also be applicable in qPCR platforms using Taqman probes and other model systems including cell lines and rodent models. Absence of sample from healthy-normal individual in diagnostic cases call for careful selection of Reference Genes for relative quantification of a biomarker by qPCR.BCL2 can be used as molecular diagnostics only if normalized with a set of reference genes with stable yet low levels of expression across different types of haematological malignancies.