Genome-scale analysis of in vivo spatiotemporal promoter activity in Caenorhabditis elegans

Proton Microbeam Targeted Irradiation of the Gonad Primordium Region Induces Developmental Alterations Associated with Heat Shock Responses and Cuticle Defense in Caenorhabditis elegans

We describe a methodology to manipulate Caenorhabditis elegans (C. elegans) and irradiate the stem progenitor gonad region using three MeV protons at a specific developmental stage (L1). The consequences of the targeted irradiation were first investigated by considering the organogenesis of the vulva and gonad, two well-defined and characterized developmental systems in C. elegans. In addition, we adapted high-throughput analysis protocols, using cell-sorting assays (COPAS) and whole transcriptome analysis, to the limited number of worms (>300) imposed by the selective irradiation approach. Here, the presented status report validated protocols to (i) deliver a controlled dose in specific regions of the worms; (ii) immobilize synchronized worm populations (>300); (iii) specifically target dedicated cells; (iv) study the radiation-induced developmental alterations and gene induction involved in cellular stress (heat shock protein) and cuticle injury responses that were found.

Lipid droplets modulate proteostasis, SQST-1/SQSTM1 dynamics, and lifespan in C. elegans

In several long-lived Caenorhabditis elegans strains, such as insulin/IGF-1 receptor daf-2 mutants, enhanced proteostatic mechanisms are accompanied by elevated intestinal lipid stores, but their role in longevity is unclear. Here, while determining the regulatory network of the selective autophagy receptor SQST-1/SQSTM1, we uncovered an important role for lipid droplets in proteostasis and longevity. Using genome-wide RNAi screening, we identified several SQST-1 modulators, including lipid droplets-associated and aggregation-prone proteins. Expansion of intestinal lipid droplets by silencing the conserved cytosolic triacylglycerol lipase gene atgl-1/ATGL enhanced autophagy, and extended lifespan. Notably, a substantial amount of ubiquitinated proteins were found on lipid droplets. Reducing lipid droplet levels exacerbated the proteostatic collapse when autophagy or proteasome function was compromised, and significantly reduced the lifespan of long-lived daf-2 animals. Altogether, our study uncovered a key role for lipid droplets in C. elegans as a proteostatic mediator that modulates ubiquitinated protein accumulation, facilitates autophagy, and promotes longevity.

VISTA: visualizing the spatial transcriptome of the C.elegans nervous system

Summary

Profiling the transcriptomes of single cells without sacrificing spatial information is a major goal of the field of spatial transcriptomics, but current technologies require tradeoffs between single-cell resolution and whole-transcriptome coverage. In one animal species, the nematode worm Caenorhabditis elegans, a comprehensive spatial transcriptome with single-cell resolution is attainable using existing datasets, thanks to the worm's invariant cell lineage and a series of recently generated single cell transcriptomes. Here we present VISTA, which leverages these datasets to provide a visualization of the worm spatial transcriptome, focusing specifically on the nervous system. VISTA allows users to input a query gene and visualize its expression across all neurons in the form of a "spatial heatmap" in which the color of a cell reports the expression level. Underlying gene expression values (in Transcripts Per Million) are displayed when an individual cell is selected. We provide examples of the utility of VISTA for identifying striking new gene expression patterns in specific neurons, and for resolving cellular identities of ambiguous expression patterns generated from in vivo reporter genes. The ability to easily obtain gene-level snapshots of the neuronal spatial transcriptome should facilitate studies on neuron-specific gene expression and regulation and provide a template for the high-resolution spatial transcriptomes the field hopes to obtain for various animal species in the future.

Availability and implementation

VISTA is freely available at the following URL: https://public.tableau.com/app/profile/smu.oit.data.insights/viz/VISTA_16814210566130/VISTA.

VISTA: Visualizing the Spatial Transcriptome of the C. elegans Nervous System

Profiling the transcriptomes of single cells without sacrificing spatial information is a major goal of the field of spatial transcriptomics, but current technologies require tradeoffs between single-cell resolution and whole-transcriptome coverage. In one animal species, the nematode worm C. elegans, a comprehensive spatial transcriptome with single-cell resolution is attainable using existing datasets, thanks to the worm's invariant cell lineage and a series of recently-generated single cell transcriptomes. Here we present VISTA, which leverages these datasets to provide a visualization of the worm spatial transcriptome, focusing specifically on the nervous system. VISTA allows users to input a query gene and visualize its expression across all neurons in the form of a "spatial heatmap" in which the color of a cell reports the expression level. Underlying gene expression values (in Transcripts Per Million) are displayed when an individual cell is selected. We provide examples of the utility of VISTA for identifying striking new gene expression patterns in specific neurons, and for resolving cellular identities of ambiguous expression patterns generated from in vivo reporter genes. The ability to easily obtain gene-level snapshots of the neuronal spatial transcriptome should facilitate studies on neuron-specific gene expression and regulation, and provide a template for the high-resolution spatial transcriptomes the field hopes to obtain for various animal species in the future.

An expandable FLP-ON::TIR1 system for precise spatiotemporal protein degradation in Caenorhabditis elegans

The auxin-inducible degradation system has been widely adopted in the Caenorhabditis elegans research community for its ability to empirically control the spatiotemporal expression of target proteins. This system can efficiently degrade auxin-inducible degron (AID)-tagged proteins via the expression of a ligand-activatable AtTIR1 protein derived from A. thaliana that adapts target proteins to the endogenous C. elegans proteasome. While broad expression of AtTIR1 using strong, ubiquitous promoters can lead to rapid degradation of AID-tagged proteins, cell type-specific expression of AtTIR1 using spatially restricted promoters often results in less efficient target protein degradation. To circumvent this limitation, we have developed an FLP/FRT3-based system that functions to reanimate a dormant, high-powered promoter that can drive sufficient AtTIR1 expression in a cell type-specific manner. We benchmark the utility of this system by generating a number of tissue-specific FLP-ON::TIR1 drivers to reveal genetically separable cell type-specific phenotypes for several target proteins. We also demonstrate that the FLP-ON::TIR1 system is compatible with enhanced degron epitopes. Finally, we provide an expandable toolkit utilizing the basic FLP-ON::TIR1 system that can be adapted to drive optimized AtTIR1 expression in any tissue or cell type of interest.

A preliminary analysis of prognostic genes in advanced laryngeal squamous cell carcinoma patients with postoperative radiotherapy

Advanced laryngeal squamous cell carcinoma (LSCC) has a high mortality rate, and the prognosis is poor. However, the underlying molecular biological mechanisms bringing about the development and progression of advanced LSCC are not entirely clarified. This study aimed to find out the potential biomarkers to predict the prognosis in advanced LSCC patients who had undergone postoperative radiotherapy alone. The next-generation sequencing of RNA was performed to detect the mRNAs expression profiling in 10 advanced LSCC samples, comprised of 5 samples from LSCC patients with favorable outcome and 5 samples from paired patients with poor outcome. Then bioinformatics analysis including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were used to find out functional core genes that were significantly different between the two groups. 1630 differentially expressed genes (DEGs) were confirmed to have significant differences between the two groups. 53 GO terms and 19 pathways which were closely related to the DEGs were identified. Finally, 52 intersection DEGs which were both related to the top three GO terms and pathways were identified. The expression of several core genes was confirmed with RT-qPCR in tissues from another 75 patients. RT-qPCR confirmed that the genes of c-JUN, LYN, PIK3R2, and TNFAIP3 were significantly differentially expressed between the two groups, which was in accordance with the RNA sequencing data. The DEGs identified above may be potential prognostic markers for advanced LSCC patients with postoperative radiotherapy, and may provide essential guidance for following-up.

Whole transcriptome expression profiles in kidney samples from rats with hyperuricaemic nephropathy

Hyperuricaemic nephropathy (HN) is a common clinical complication of hyperuricaemia (HUA) and poses a huge threat to human health. Hence, we aimed to prospectively investigate the dysregulated genes, pathways and networks involved in HN by performing whole transcriptome sequencing using RNA sequencing. Six kidney samples from HN group (n = 3) and a control group (n = 3) were obtained to conduct RNA sequencing. To disclose the relevant signalling pathways, we conducted the analysis of differentially expressed genes (DEGs), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. A competitive endogenous RNA (ceRNA) network was established to reveal the interactions between lncRNAs, circRNAs, mRNAs and miRNAs and investigate the potential mechanisms of HN. Ultimately, 2250 mRNAs, 306 lncRNAs, 5 circRNAs, and 70 miRNAs were determined to be significantly differentially expressed in the HN group relative to the control group. We further authenticated 8 differentially expressed (DE)-ncRNAs by quantitative real-time polymerase chain reaction, and these findings were in accordance with the sequencing results. The analysis results evidently showed that these DE-ncRNAs were significantly enriched in pathways related to inflammatory reaction. In conclusion, HUA may generate abnormal gene expression changes and regulate signalling pathways in kidney samples. Potentially related genes and pathways involved in HN were identified.

Unconventional Myosins from Caenorhabditis elegans as a Probe to Study Human Orthologues

Unconventional myosins are a superfamily of actin-based motor proteins that perform a number of roles in fundamental cellular processes, including (but not limited to) intracellular trafficking, cell motility, endocytosis, exocytosis and cytokinesis. 40 myosins genes have been identified in humans, which belong to different 12 classes based on their domain structure and organisation. These genes are widely expressed in different tissues, and mutations leading to loss of function are associated with a wide variety of pathologies while over-expression often results in cancer. Caenorhabditis elegans (C. elegans) is a small, free-living, non-parasitic nematode. ~38% of the genome of C. elegans has predicted orthologues in the human genome, making it a valuable tool to study the function of human counterparts and human diseases. To date, 8 unconventional myosin genes have been identified in the nematode, from 6 different classes with high homology to human paralogues. The hum-1 and hum-5 (heavy chain of an unconventional myosin) genes encode myosin of class I, hum-2 of class V, hum-3 and hum-8 of class VI, hum-6 of class VII and hum-7 of class IX. The hum-4 gene encodes a high molecular mass myosin (307 kDa) that is one of the most highly divergent myosins and is a member of class XII. Mutations in many of the human orthologues are lethal, indicating their essential properties. However, a functional characterisation for many of these genes in C. elegans has not yet been performed. This article reviews the current knowledge of unconventional myosin genes in C. elegans and explores the potential use of the nematode to study the function and regulation of myosin motors to provide valuable insights into their role in diseases.

Identification of Competing Endogenous RNAs (ceRNAs) Network Associated with Drought Tolerance in Medicago truncatula with Rhizobium Symbiosis

Drought, bringing the risks of agricultural production losses, is becoming a globally environmental stress. Previous results suggested that legumes with nodules exhibited superior drought tolerance compared with the non-nodule group. To investigate the molecular mechanism of rhizobium symbiosis impacting drought tolerance, transcriptome and sRNAome sequencing were performed to identify the potential mRNA-miRNA-ncRNA dynamic network. Our results revealed that seedlings with active nodules exhibited enhanced drought tolerance by reserving energy, synthesizing N-glycans, and medicating systemic acquired resistance due to the early effects of symbiotic nitrogen fixation (SNF) triggered in contrast to the drought susceptible with inactive nodules. The improved drought tolerance might be involved in the decreased expression levels of miRNA such as mtr_miR169l-5p, mtr_miR398b, and mtr_miR398c and its target genes in seedlings with active nodules. Based on the negative expression pattern between miRNA and its target genes, we constructed an mRNA-miR169l-ncRNA ceRNA network. During severe drought stress, the lncRNA alternative splicings TCONS_00049507 and TCONS_00049510 competitively interacted with mtr_miR169l-5p, which upregulated the expression of NUCLEAR FACTOR-Y (NF-Y) transcription factor subfamily NF-YA genes MtNF-YA2 and MtNF-YA3 to regulate their downstream drought-response genes. Our results emphasized the importance of SNF plants affecting drought tolerance. In conclusion, our work provides insight into ceRNA involvement in rhizobium symbiosis contributing to drought tolerance and provides molecular evidence for future study.

The Smart Soil Organism Detector: An instrument and machine learning pipeline for soil species identification

Understanding the diversity of soil organisms is complicated by both scale and substrate. Every footprint we leave in the soil covers hundreds to millions of organisms yet we cannot see them without extremely laborious extraction and microsopy endeavors. Studying them is also challenging. Keeping them alive so that we can understand their lifecycles and ecological roles ranges from difficult to impossible. Functional and taxonomic identification of soil organisms, while possible, is also challenging. Here we present the Smart Soil Organism Detector, an instrument and machine learning pipeline that combines high-resolution imaging, multi-spectral sensing, large-bore flow cytometry, and machine learning to extract, isolate, count, identify, and separate soil organisms in a high-throughput, high-resolution, non-destructive, and reproducible manner. This system is not only capable of separating alive nematodes, dead nematodes, and nematode cuticles from soil with 100% out-of-sample accuracy, but also capable of identifying nematode strains (sub-species) with 95.5% out-of-sample accuracy and 99.4% specificity. Soil micro-arthropods were identified to class with 96.1% out-of-sample accuracy. Broadly applicable across soil taxa, the Smart SOD system is a tool for understanding global soil biodiversity.

ACE2 pathway regulates thermogenesis and energy metabolism

Identification of key regulators of energy homeostasis holds important therapeutic promise for metabolic disorders, such as obesity and diabetes. ACE2 cleaves angiotensin II (Ang II) to generate Ang-(1-7) which acts mainly through the Mas1 receptor. Here, we identify ACE2 pathway as a critical regulator in the maintenance of thermogenesis and energy expenditure. We found that ACE2 is highly expressed in brown adipose tissue (BAT) and that cold stimulation increases ACE2 and Ang-(1-7) levels in BAT and serum. Ace2 knockout mice (Ace2^-/y) and Mas1 knockout mice (Mas1^-/-) displayed impaired thermogenesis. Mice transplanted with brown adipose tissue from Mas1^-/- display metabolic abnormalities consistent with those seen in the Ace2 and Mas1 knockout mice. In contrast, impaired thermogenesis of Lepr^db/db obese diabetic mice and high-fat diet-induced obese mice were ameliorated by overexpression of Ace2 or continuous infusion of Ang-(1-7). Activation of ACE2 pathway was associated with improvement of metabolic parameters, including blood glucose, lipids, and energy expenditure in multiple animal models. Consistently, ACE2 pathway remarkably enhanced the browning of white adipose tissue. Mechanistically, we showed that ACE2 pathway activated Akt/FoxO1 and PKA pathway, leading to induction of UCP1 and activation of mitochondrial function. Our data propose that adaptive thermogenesis requires regulation of ACE2 pathway and highlight novel potential therapeutic targets for the treatment of metabolic disorders.

Hippocampus chronic deep brain stimulation induces reversible transcript changes in a macaque model of mesial temporal lobe epilepsy

BACKGROUND

Deep brain stimulation (DBS) has seizure-suppressing effects but the molecular mechanisms underlying its therapeutic action remain unclear. This study aimed to systematically elucidate the mechanisms underlying DBS-induced seizure suppression at a molecular level.

METHODS

We established a macaque model of mesial temporal lobe epilepsy (mTLE), and continuous high-frequency hippocampus DBS (hip-DBS) was applied for 3 months. The effects of hip-DBS on hippocampus gene expression were examined using high-throughput microarray analysis followed by bioinformatics analysis. Moreover, the microarray results were validated using quantitative real-time polymerase chain reaction (qRT-PCR) and Western blot analyses.

RESULTS

The results showed that chronic hip-DBS modulated the hippocampal gene expression. We identified 4119 differentially expressed genes and assigned these genes to 16 model profiles. Series test of cluster analysis showed that profiles 5, 3, and 2 were the predominant expression profiles. Moreover, profile 5 was mainly involved in focal adhesion and extracellular matrix-receptor interaction pathway. Nine dysregulated genes (Arhgap5, Col1a2, Itgb1, Pik3r1, Lama4, Fn1, Col3a1, Itga9, and Shc4) and three genes (Col1a2, Itgb1, and Flna) in these two pathways were further validated by qRT-PCR and Western blot analyses, respectively, which showed a concordance.

CONCLUSION

Our findings suggest that hip-DBS could markedly reverse mTLE-induced abnormal gene expression. Findings from this study establish the basis for further investigation of the underlying regulatory mechanisms of DBS for mTLE.

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.

Aberrant lncRNA Profiles Are Associated With Chronic Benzene Poisoning and Acute Myelocytic Leukemia

OBJECTIVE

This study investigates the mechanisms of benzene hematotoxicity.

METHODS

We used microarray to detect expression profiles of long non-coding RNAs (lncRNAs) and mRNAs in peripheral lymphocytes from chronic benzene poisoning, acute myelocytic leukemia, and healthy controls. The lncRNAs and mRNAs were validated using real-time quantitative PCR (RT-qPCR). Cytokinesis-block micronucleus assay was used to analyze chromosomal aberration.

RESULTS

We found 173 upregulated and 258 downregulated lncRNAs, and 695 upregulated and 804 downregulated mRNAs. The lncRNA CUST_40243 and mRNA PDGFC and CDKN1A associated with chronic benzene poisoning. Relevant inflammatory response, hematopoietic cell lineage, and cell cycle may be important pathways for the sifted lncRNAs and mRNAs. Furthermore, micronuclei frequency was significantly higher in off-post chronic benzene poisoning patients.

CONCLUSIONS

Chromosomal aberration induced by benzene exposure is irreversible. The lncRNA CUST_40243 and mRNA PDGFC and CDKN1A are related to chronic benzene poisoning.

Quercetin regulates ERα mediated differentiation of BMSCs through circular RNA

Circular RNA (circRNA) participates in regulation of gene transcription, while estrogen receptor alpha (ERα) and quercetin (QUE) positively regulate bone formation, but little is known about the correlation among circRNA, ERα and QUE. In this experiment, we created an ERα-deficient rBMSC model treated with QUE and evaluated the effects of ERα or QUE on rBMSCs, then analyzed differentially-expressed circRNAs by RNA-Seq and bioinformatics. The results showed that ERα deficiency constrained osteogenic differentiation and stimulated adipocytic differentiation of rBMSCs, while QUE abrogated those effects. We identified 136 differentially-expressed circRNAs in the Lv-shERα group and 120 differentially-expressed circRNAs in the Lv-shERα + QUE group. Thirty-two circRNAs retroregulated by ERα and QUE were involved in Rap1 and Wnt signaling, and four of them together sponged miR-326-5p, the target genes of which are osteogenic and adipogenic differentiation factors. Further study showed that over-expressed miR-326-5p could stimulate osteogenic differentiation, while attenuating adipogenic differentiation of rBMSCs. Therefore, we concluded that ERα and QUE might regulate the differentiation of rBMSCs through the circRNA-miR-326-5p-mRNA axis.

Transcriptomic study of lipopolysaccharide-induced sepsis damage in a mouse heart model

Sepsis is an emergency systemic illness caused by pathogen infection and the combined result of the underactivity and overactivity of a patient's own immune system. However, the molecular mechanism of this illness remains largely unknown. Lipopolysaccharide (LPS) was injected to establish a sepsis model, and heart tissue was used to analyze transcriptome changes in mice. LPS injection was used to develop a sepsis model, which resulted in cardiac tissue rearrangement and inflammatory response activation. An RNA-sequencing-based transcriptome assay using mouse heart tissue with or without LPS injection showed that 3,326 and 1,769 genes were upregulated and downregulated, respectively (>2-fold changes; P<0.05). Furthermore, these differentially expressed genes were classified into 20 pathways, including ‘Wnt signaling pathway’, ‘VEGF signaling pathway’ and ‘TGF-β signaling pathway’, and these altered genes were enriched in 41 Gene Ontology terms. The application of Wnt3a inhibited the activation of the LPS-induced inflammatory response and activated Wnt signaling, as well as protecting against LPS-mediated cardiac tissue damage in mice. In contrast, inhibition of Wnt signaling by injection of its inhibitor IWR induced plasminogen activator inhibitor-1 expression and resulted in cardiac structure derangement, which was similar to the symptoms caused by injection of LPS, suggesting that LPS-induced damage to heart tissue may be via inhibition of Wnt signaling. The present analyses showed that Wnt signaling serves a pivotal role in sepsis development and may improve our understanding of the molecular basis underlying sepsis.

Regulatory systems that mediate the effects of temperature on the lifespan of Caenorhabditis elegans

Temperature affects animal physiology, including aging and lifespan. How temperature and biological systems interact to influence aging and lifespan has been investigated using model organisms, including the nematode Caenorhabditis elegans. In this review, we discuss mechanisms by which diverse cellular factors modulate the effects of ambient temperatures on aging and lifespan in C. elegans. C. elegans thermosensory neurons alleviate lifespan-shortening effects of high temperatures via sterol endocrine signaling and probably through systemic regulation of cytosolic proteostasis. At low temperatures, C. elegans displays a long lifespan by upregulating the cold-sensing TRPA channel, lipid homeostasis, germline-mediated prostaglandin signaling, and autophagy. In addition, co-chaperone p23 amplifies lifespan changes affected by high and low temperatures. Our review summarizes how external temperatures modulate C. elegans lifespan and provides information regarding responses of biological processes to temperature changes, which may affect health and aging at an organism level.

Fluorescence-based sorting of Caenorhabditis elegans via acoustofluidics

Effectively isolating and categorizing large quantities of Caenorhabditis elegans (C. elegans) based on different phenotypes is important for most worm research, especially genetics. Here we present an integrated acoustofluidic chip capable of identifying worms of interest based on expression of a fluorescent protein in a continuous flow and then separate them accordingly in a high-throughput manner. Utilizing planar fiber optics as the detection unit, our acoustofluidic device requires no temporary immobilization of worms for interrogation/detection, thereby improving the throughput. Implementing surface acoustic waves (SAW) as the sorting unit, our device provides a contact-free method to move worms of interest to the desired outlet, thus ensuring the biocompatibility for our chip. Our device can sort worms of different developmental stages (L3 and L4 stage worms) at high throughput and accuracy. For example, L3 worms can be processed at a throughput of around 70 worms per min with a sample purity over 99%, which remains over 90% when the throughput is increased to around 115 worms per min. In our acoustofluidic chip, the time period to complete the detection and sorting of one worm is only 50 ms, which outperforms nearly all existing microfluidics-based worm sorting devices and may be further reduced to achieve higher throughput.

A Persistence Detector for Metabolic Network Rewiring in an Animal

Biological systems must possess mechanisms that prevent inappropriate responses to spurious environmental inputs. Caenorhabditis elegans has two breakdown pathways for the short-chain fatty acid propionate: a canonical, vitamin B12-dependent pathway and a propionate shunt that is used when vitamin B12 levels are low. The shunt pathway is kept off when there is sufficient flux through the canonical pathway, likely to avoid generating shunt-specific toxic intermediates. Here, we discovered a transcriptional regulatory circuit that activates shunt gene expression upon propionate buildup. Nuclear hormone receptor 10 (NHR-10) and NHR-68 function together as a "persistence detector" in a type 1, coherent feed-forward loop with an AND-logic gate to delay shunt activation upon propionate accumulation and to avoid spurious shunt activation in response to a non-sustained pulse of propionate. Together, our findings identify a persistence detector in an animal, which transcriptionally rewires propionate metabolism to maintain homeostasis.

Transcriptome Analysis of Rice Roots in Response to Root-Knot Nematode Infection

Meloidogyne incognita and Meloidogyne graminicola are root-knot nematodes (RKNs) infecting rice (Oryza sativa L.) roots and severely decreasing yield, whose mechanisms of action remain unclear. We investigated RKN invasion and development in rice roots through RNA-seq transcriptome analysis. The results showed that 952 and 647 genes were differently expressed after 6 (invasion stage) and 18 (development stage) days post inoculation, respectively. Gene annotation showed that the differentially expressed genes were classified into diverse metabolic and stress response categories. Furthermore, phytohormone, transcription factor, redox signaling, and defense response pathways were enriched upon RKN infection. RNA-seq validation using qRT-PCR confirmed that CBL-interacting protein kinase (CIPK) genes (CIPK5, 8, 9, 11, 14, 23, 24, and 31) as well as brassinosteroid (BR)-related genes (OsBAK1, OsBRI1, D2, and D11) were altered by RKN infection. Analysis of the CIPK9 mutant and overexpressor indicated that the RKN populations were smaller in cipk9 and larger in CIPK9 OX, while more galls were produced in CIPK9 OX plant roots than the in wild-type roots. Significantly fewer numbers of second-stage infective juveniles (J2s) were observed in the plants expressing the BR biosynthesis gene D2 mutant and the BR receptor BRI1 activation-tagged mutant (bri1-D), and fewer galls were observed in bri1-D roots than in wild-type roots. The roots of plants expressing the regulator of ethylene signaling ERS1 (ethylene response sensor 1) mutant contained higher numbers of J2s and developed more galls compared with wild-type roots, suggesting that these signals function in RKN invasion or development. Our findings broaden our understanding of rice responses to RKN invasion and provide useful information for further research on RKN defense mechanisms.

Expression analysis of microRNAs and mRNAs in myofibroblast differentiation of lung resident mesenchymal stem cells

Idiopathic pulmonary fibrosis (IPF) is a serious lung disease that involved the myofibroblast differentiation of lung resident mesenchymal stem cells (LR-MSCs). However, the specific molecular mechanisms of myofibroblast differentiation of LR-MSCs still remain a mystery. In this study, a comprehensive analysis of miRNAs and mRNAs changes in LR-MSCs treated with TGF-β1 was performed. Through computational approaches, the pivotal roles of differentially expressed miRNAs that were associated with tight junction, pathways in cancer, focal adhesion, and cytokine-cytokine receptor interaction were shown. Kruppel-like factor 4 (Klf4) and inhibitor of growth family, member 5 (Ing5) may be the targets for the therapy of pulmonary fibrosis by inhibiting myofibroblast differentiation of LR-MSCs and EMT. Collectively, a molecular paradigm for understanding myofibroblast differentiation of LR-MSCs in IPF was provided by the integrated miRNA/mRNA analyses.

Large-scale analyses identify a cluster of novel long noncoding RNAs as potential competitive endogenous RNAs in progression of hepatocellular carcinoma

The abnormal expression of noncoding RNAs has attracted increasing interest in the field of hepatocellular carcinoma progression. However, the underlying molecular mechanisms mediated by noncoding RNAs in these processes are unclear. Here, we obtained the expression profiles of long noncoding RNAs, microRNAs, and mRNAs from the Gene Expression Omnibus database and identified hepatocarcinogenesis-specific differentially expressed transcripts. Next, we identified significant Gene Ontology and pathway terms that the differentially expressed transcripts involved in. Using functional analysis and target prediction, we constructed a hepatocellular carcinoma-associated deregulated competitive endogenous RNA network to reveal the potential mechanisms underlying tumor progression. By analyzing The Cancer Genome Atlas dataset, six key long noncoding RNAs showed significant association with overall survival as well as strong correlation with some microRNAs and mRNAs in the competitive endogenous RNA network. We further validated the above results and determined their diagnostic and prognostic value in clinical samples. Importantly, by large-scale analyses, we identified a cluster of long noncoding RNAs, GBAP1, MCM3AP-AS1, SLC16A1-AS1, C3P1, DIO3OS, and HNF4A-AS1 as candidate biomarkers for the diagnosis and prognosis of hepatocellular carcinoma, which will improve our understanding of competitive endogenous RNA-mediated regulatory mechanisms underlying hepatocellular carcinoma development and will provide novel therapeutic targets in the future.

Molecular genetic and biochemical characterization of a putative family of zinc metalloproteins in Caenorhabditis elegans

The first characterization of W08E12.2, W08E12.3, W08E12.4 and W08E12.5, four putative metalloproteins in C. elegans. (A) phase contrast microscopy, (B) fluorescence microscopy of PW08E12.3;W08E12.4::GFP.

Four highly similar genes (W08E12.2, W08E12.3, W08E12.4 and W08E12.5) which are consecutively aligned on chromosome IV of the C. elegans genome are predicted to code for small (120–141aa) yet cysteine rich (18–19 cysteines) proteins. Cloning and sequencing of the genomic regions of the isoforms confirmed the presence and order of all genes. The generation of transgenic worms strains with an integrated single copy or extrachromosomal multi-copy PW08E12.3;W08E12.4::GFP uncovered that W08E12.3 and W08E12.4 are constitutively expressed in the pharynx and significantly induced in worms exposed to 100 μM Zn. Knockdown by RNAi did not have a marked consequence on reproductive performance nor was a Zn-dependent effect on nematode growth observed. However, RNAi of these genes led to an accumulation of Zn in the intestinal cells. W08E12.3 was recombinantly expressed in E. coli and the purified protein was shown to be able to bind up to 6.5 Zn molecules at neutral pH. Zn-binding was acid-labile and the apo protein was observed at pH < 4.3. This characterization suggests W08E12.2, W08E12.3, W08E12.4 and W08E12.5 belong to a family of putative Metalloproteins which, akin to metallothioneins, may play an important role in Zn-sensing, homeostasis and/or detoxification.

Duplication of hsp-110 Is Implicated in Differential Success of Globodera Species under Climate Change

Managing the emergence and spread of crop pests and pathogens is essential for global food security. Understanding how organisms have adapted to their native climate is key to predicting the impact of climate change. The potato cyst nematodes Globodera pallida and G. rostochiensis are economically important plant pathogens that cause yield losses of up to 50% in potato. The two species have different thermal optima that may relate to differences in the altitude of their regions of origin in the Andes. Here, we demonstrate that juveniles of G. pallida are less able to recover from heat stress than those of G. rostochiensis. Genome-wide analysis revealed that while both Globodera species respond to heat stress by induction of various protective heat-inducible genes, G. pallida experiences heat stress at lower temperatures. We use C. elegans as a model to demonstrate the dependence of the heat stress response on expression of Heat Shock Factor-1 (HSF-1). Moreover, we show that hsp-110 is induced by heat stress in G. rostochiensis, but not in the less thermotolerant G. pallida. Sequence analysis revealed that this gene and its promoter was duplicated in G. rostochiensis and acquired thermoregulatory properties. We show that hsp-110 is required for recovery from acute thermal stress in both C. elegans and in G. rostochiensis. Our findings point towards an underlying molecular mechanism that allows the differential expansion of one species relative to another closely related species under current climate change scenarios. Similar mechanisms may be true of other invertebrate species with pest status.

lncRNA Gm10451 regulates PTIP to facilitate iPSCs-derived β-like cell differentiation by targeting miR-338-3p as a ceRNA

iPSCs-derived insulin-producing cell transplantation is a promising strategy for diabetes therapy. Although there have been many protocols of mature, glucose-responsive β cells induced in vitro over the past few years, many underlying problems remain to be resolved. As a crucial regulator, long noncoding RNAs (lncRNAs) participate in numerous biological processes, including the maintenance of pluripotency, and stem cell differentiation. In this study, we identified a novel lncRNA Gm10451 as a functional regulator for β-like cell differentiation. Localized to the cytoplasm, Gm10451 regulates histone H3K4 methyltransferase complex PTIP to facilitate Insulin⁺/Nkx6.1⁺ β-like cell differentiation by targeting miR-338-3p as a competing endogenous RNA (ceRNA). miR-338-3p has also been shown to suppress Nkx6.1⁺ early-stage β-like cell differentiation by targeting PTIP. Following transplantation into streptozotocin (STZ)-mice, Gm10451 loss in β-like cells prevented the expression of mature β-cell makers, such as Insulin, Nkx6.1, and Mafa. Accordingly, hyperglycemia in the mice was not resolved. Taken together, this study provides an efficient epigenetic target for generating more mature and functional iPSCs-derived β-like cells. We anticipate that pancreatic organoids, which are generated from human stem cells, biological materials, and epigenetic modifications, can be used in the future as a novel diabetes treatment option.

Long noncoding RNA uc.4 inhibits cell differentiation in heart development by altering DNA methylation

In previous studies, we have demonstrated that long noncoding RNA uc.4 may influence the cell differentiation through the TGF-β signaling pathway, suppressed the heart development of zebrafish and resulting cardiac malformation. DNA methylation plays a significant role in the heart development and disordered of DNA methylation may cause disruption of control of gene promoter. In this study, methylated DNA immunoprecipitation was performed to identify the different expression levels of methylation regions. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were also performed to identify the possible biological process and pathway that uc.4 may join, associated with Rap1 signaling pathway, gonadotropin-releasing hormone signaling pathway, and Calcium signaling pathway. We found that the distribution of differentially methylated regions peaks was mainly located in intergenic and intron regions. Altogether, our result showed that differentially methylated genes are significantly expressed in uc.4-overexpression cells, providing valuable data for further exploration of the role of uc.4 in heart development.

Screening common signaling pathways associated with drug resistance in non-small cell lung cancer via gene expression profile analysis

Lung cancer is the leading cause of cancer‐related deaths worldwide. Although several therapeutic strategies have been employed to curb lung cancer, the survival rate is still poor owing to the development of drug resistance. The mechanisms underlying drug resistance development are incompletely understood. Here, we aimed to identify the common signaling pathways involved in drug resistance in non‐small cell lung cancer (NSCLC). Three published transcriptome microarray data were downloaded from the Gene Expression Omnibus (GEO) database comprising different drug‐resistant cell lines and their parental cell lines. Differentially expressed genes (DEGs) were identified and used to perform Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. An overlapping analysis was performed for KEGG pathways enriched from all the three datasets to identify the common signaling pathways. As a result, we found that metabolic pathways, ubiquitin‐mediated proteolysis, and mitogen‐activated protein kinase (MAPK) signaling were the most aberrantly expressed signaling pathways. The knockdown of nicotinamide phosphoribosyltransferase (NAMPT), the gene involved in metabolic pathways and known to be upregulated in drug‐resistant tumor cells, was shown to increase the apoptosis of cisplatin‐resistant A549 cells following cisplatin treatment. Thus, our results provide an in‐depth analysis of the signaling pathways that are commonly altered in drug‐resistant NSCLC cell lines and highlight the potential strategy that facilitates the development of interventions to interfere with upregulated signaling pathways as well as to boost downregulated signaling pathways in drug‐resistant tumors for the elimination of multiple resistance of NSCLC.

Transcriptomic changes associated with PCK1 overexpression in hepatocellular carcinoma cells detected by RNA-seq

Phosphoenolpyruvate carboxykinase 1 (PCK1), a step limiting enzyme of gluconeogenesis, is downregulated in hepatocellular carcinoma (HCC). Overexpression of PCK1 has been shown to suppress hepatoma cell growth, but the underlying mechanism remains unclear. We used recombinant adenovirus overexpressing PCK1 or GFP in Huh7 cells, and the differentially expressed genes (DEGs) were identified by RNA-Seq. 180 were upregulated by PCK1 overexpression, whereas 316 were downregulated. Pathway analysis illustrated that PCK1 was closely correlated with Wnt signaling pathway and TGF-beta signaling pathway. Hence, Wnt signaling pathway and its downstream component, FZD2, FZD6, FZD7 and β-catenin were confirmed by qRT-PCR and Western blot. In vivo we also observed that PCK1 had restrained tumor growth as a result of decreasing expression of β-catenin. Whole-transcriptomic profile analysis discovered that overexpression of PCK1 downregulates several oncogenic signaling pathways in HCC, providing potential therapeutic targets for improving HCC therapy.

Integrated analysis of dysregulated long non-coding RNAs/microRNAs/mRNAs in metastasis of lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD), largely remains a primary cause of cancer-related death worldwide. The molecular mechanisms in LUAD metastasis have not been completely uncovered.

METHODS

In this study, we identified differentially expressed genes (DEGs), miRNAs (DEMs) and lncRNAs (DELs) underlying metastasis of LUAD from The Cancer Genome Atlas database. Intersection mRNAs were used to perform gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway and co-expression network analysis. In addition, survival analyses of intersection mRNAs were conducted. Finally, intersection mRNAs, miRNAs and lncRNAs were subjected to construct miRNA-mRNA-lncRNA network.

RESULTS

A total of 1015 DEGs, 54 DEMs and 22 DELs were identified in LUAD metastasis and non-metastasis samples. GO and KEGG pathway analysis had proven that the functions of intersection mRNAs were closely related with many important processes in cancer pathogenesis. Among the co-expression interactions network, 22 genes in the co-expression network were over the degree 20. These genes imply that they have connections with many other gene nodes. In addition, 14 target genes (ARHGAP11A, ASPM, HELLS, PRC1, TMPO, ARHGAP30, CD52, IL16, IRF8, P2RY13, PRKCB, PTPRC, SASH3 and TRAF3IP3) were found to be associated with survival in patients with LUAD significantly (log-rank P < 0.05). Two lncRNAs (LOC96610 and ADAM6) acting as ceRNAs were identified based on the miRNA-mRNA-lncRNA network.

CONCLUSIONS

Taken together, the results may provide a novel perspective to develop a multiple gene diagnostic tool for LUAD prognosis, which might also provide potential biomarkers or therapeutic targets for LUAD.

CRISPR-Cas9 human gene replacement and phenomic characterization in Caenorhabditis elegans to understand the functional conservation of human genes and decipher variants of uncertain significance

Our ability to sequence genomes has vastly surpassed our ability to interpret the genetic variation we discover. This presents a major challenge in the clinical setting, where the recent application of whole-exome and whole-genome sequencing has uncovered thousands of genetic variants of uncertain significance. Here, we present a strategy for targeted human gene replacement and phenomic characterization, based on CRISPR-Cas9 genome engineering in the genetic model organism Caenorhabditis elegans, that will facilitate assessment of the functional conservation of human genes and structure-function analysis of disease-associated variants with unprecedented precision. We validate our strategy by demonstrating that direct single-copy replacement of the C. elegans ortholog (daf-18) with the critical human disease-associated gene phosphatase and tensin homolog (PTEN) is sufficient to rescue multiple phenotypic abnormalities caused by complete deletion of daf-18, including complex chemosensory and mechanosensory impairments. In addition, we used our strategy to generate animals harboring a single copy of the known pathogenic lipid phosphatase inactive PTEN variant (PTEN-G129E), and showed that our automated in vivo phenotypic assays could accurately and efficiently classify this missense variant as loss of function. The integrated nature of the human transgenes allows for analysis of both homozygous and heterozygous variants and greatly facilitates high-throughput precision medicine drug screens. By combining genome engineering with rapid and automated phenotypic characterization, our strategy streamlines the identification of novel conserved gene functions in complex sensory and learning phenotypes that can be used as in vivo functional assays to decipher variants of uncertain significance.

Gene profiles to characterize the combined toxicity induced by low level co-exposure of silica nanoparticles and benzo[a]pyrene using whole genome microarrays in zebrafish embryos

Several studies have suggested that air pollutants combine exposure have greater adverse effects. However, limited studies were available on the combined toxicity of silica nanoparticles (SiNPs) and benzo[a]pyrene (B[a]P). The study was to evaluate the toxic effect and mechanisms of low-dose exposure of SiNPs, B[a]P and co-exposure in zebrafish embryos. In this study, zebrafish embryos received intravenous microinjection of SiNPs and B[a]P, and then was used to select differentially expressed genes by microarray analysis. Multiple bioinformatics analyses and STC analysis were done to identify key genes, pathways and biological processes and the expression trend of genes in each group. 1) 3065 differentially expressed genes were identified in zebrafish embryos. 2) These differentially expressed genes were involved in multiple biological processes and cellular processes such as immunity, response to stimuli, cell proliferation, adhesion, signaling transduction, and embryonic development. 3) Dynamic Gene Network analysis was used to identify a subgroup of 26 core genes that involved in multiple biological processes and cellular processes. 4) Pathway analysis and Signal-net analysis indicated that the MAPK signaling pathway, calcium signaling pathway, p53 signaling pathway, PI3k/Akt signaling pathway, and several pathways associated with immune response were the most prominent significant pathways induced by co-exposure of SiNPs and B[a]P in zebrafish embryos. Our study demonstrated that the molecular actions of co-treated with SiNPs and B[a]P on the immune system, inflammatory process and cardiovascular development had more severe toxicity than single exposure.

Gene expression profiling of human bone marrow mesenchymal stem cells during osteogenic differentiation

OBJECTIVE

Osteogenesis is a multiple-step process through which osteoblasts are derived from bone marrow mesenchymal stem cells (MSCs) with multilineage differentiation potential. This study aimed to analyze gene expression profiling during osteogenic differentiation of MSCs.

MATERIALS AND METHODS

Human MSCs were isolated and induced for differentiation in osteogenic medium. Full-genome gene expression microarrays and gene ontology analysis were performed.

RESULTS

A total of 1,680 differentially expressed genes in differentiated MSCs were identified including 430 upregulated and 1,250 downregulated. Moreover, pathway-act-network analysis showed that cell cycle, p53 signaling pathway and focal adhesion, had high degree (>5). The ribonucleotide reductase M1, thymidine kinase 1 and histone cluster 1 H3e also showed high degree (>10). Polymerase chain reaction analysis confirmed the differential expression of insulin-like growth factor binding protein 3, SMAD family member 3, transforming growth factor beta 2, and fibroblast growth factor 14 in differentiated MSCs.

CONCLUSIONS

Gene expression profiling provides a foundation to reveal the mechanisms that regulate osteogenic differentiation of MSCs.

miRNA array analysis of plasma miRNA alterations in rats exposed to a high altitude hypoxic environment

In the present study, the microRNA (miRNA) expression profiles of rats exposed to high altitude hypoxia and normal conditions were obtained from miRNA array analysis. Bioinformatics analyses, including the use of the Gene Oncology and Kyoto Encyclopedia of Genes and Genomes databases, were used to identify the genes and pathways, which were specifically associated with high altitude hypoxic environment‑associated miRNAs. A total of 26 miRNAs were differentially expressed in the two groups, comprising six upregulated and 20 downregulated miRNAs. In the present study, a novel pattern of upregulated miRNAs and their associated pathways were constructed, including proteoglycans in cancer, spliceosome, gluamatergic synapse, glycolysis/gluconeogenesis, Foxo, cGMP‑PKG and p53 signaling pathways, which may provide novel targets for diagnosing and understanding the mechanism of high altitude hypoxia‑induced disease.

Altered DNA Methylation Sites in Peripheral Blood Leukocytes from Patients with Simple Steatosis and Nonalcoholic Steatohepatitis (NASH)

Background

The aim of this study was to identify DNA methylation sites in peripheral blood leukocytes from patients with histologically confirmed nonalcoholic fatty liver disease (NAFLD) that included simple hepatic steatosis and nonalcoholic steatohepatitis (NASH).

Material/Methods

DNA was isolated from peripheral blood leukocytes from patients with histologically diagnosed NAFLD (n=35), including simple hepatic steatosis (n=18) and NASH (n=17). Healthy controls included individuals without liver disease (n=30). DNA was hybridized, and DNA methylation was interrogated in an epigenome-wide association study (EWAS). DNA methylation levels (β-values) were correlated with serum lipid profiles, liver enzymes, and liver histology.

Results

Circulating blood leukocytes from 35 patients with NAFLD (simple steatosis and NASH) contained 65 CpG sites, which represented 60 genes that were differentially methylated when compared with healthy controls. In the simple hepatic steatosis group (n=18), 42 methylated CpG sites were found to be associated with increased levels of serum alanine aminotransferase (ALT), and 32 methylated CpG sites were associated with increased serum lipid profiles. In the NASH group (n=17), when compared with the simple hepatic steatosis group, methylated CpG sites showed significant correlations with the presence of lobular inflammation compared with hepatic steatosis and fibrosis. Six differentially methylated CpG sites were identified in the ACSL4, CRLS1, CTP1A, SIGIRR, SSBP1 and ZNF622 genes, which were associated with histologically confirmed simple hepatic steatosis and NASH.

Conclusions

The study identified some key methylated CpG sites from peripheral blood leukocytes, which might be used as serum biomarkers to stratify NAFLD patients into simple hepatic steatosis and NASH.

Ovarian transcriptomic analysis reveals the alternative splicing events associated with fecundity in different sheep breeds

Alternative splicing (AS) is one of the most common mechanisms that accounts for the greater macromolecular and cellular complexity of higher eukaryotic organisms. This study focused on the splicing events in the ovaries of different sheep breeds, namely the Han and Dorset breeds. Of the groups studied, cassette splicing events accounted for the maximum number of the AS events with significant differences, whereas the splicing events that were mutually exclusive with introns accounted for the smallest proportion of splicing events. Greater than 1000 AS events with significant differences were identified between the Han BB and Dorset sheep. The number of AS events with significant differences between Han ++ and Dorset sheep, however, was fewer than that of the comparison of Han BB and Dorset sheep. Seven randomly selected genes with AS events were detected in this study and were validated by reverse transcription-polymerase chain reaction (RT-PCR). Furthermore, there are many genes which were common to the two genotype groups (Han BB and Dorset sheep, as well as the Han ++ and Dorset sheep). In addition, genes detected in the present study were involved in different pathways, including the pathways related with fertility or fecundity. The present study could provide the detailed understanding on the mechanisms of alternative splicing events associated with fecundity in different sheep breeds.

RAVL1 Activates Brassinosteroids and Ethylene Signaling to Modulate Response to Sheath Blight Disease in Rice

Rhizoctonia solani causes sheath blight disease in rice; however, the defense mechanism of rice plants against R. solani remains elusive. To analyze the roles of brassinosteroid (BR) and ethylene signaling on rice defense to R. solani, wild-type (WT) rice and several mutants and overexpressing (OX) lines were inoculated with R. solani. Mutants d61-1 and d2 were less susceptible to sheath blight disease, bri1-D was more susceptible, and ravl1 and d61-1/EIL1 Ri5 were similarly susceptible compared with WT. The double mutant ravl1/d61-1 was phenotypically similar to the ravl1 mutant. Transcriptome analysis, chromatin immunoprecipitation assay, electrophoretic mobility shift assay, and transient assays indicted that RAVL1 might directly activate Ethylene insensitive 3-like 1 (EIL1), a master regulator of ethylene signaling. Mutants ers1 and d61-1/RAVL1 OX were resistant to sheath blight disease, whereas EIL1 RNAi mutants and RAVL1 OX were more susceptible than WT. BRI1 and D2 expression in EIL1 Ri5/RAVL1 OX and EIL1 expression in d61-1/RAVL1 OX indicated that RAVL1 activates BRI1/D2 and EIL1, respectively, independent of BR and ethylene signaling. Our analyses provide information on how BR and ethylene signaling regulate sheath blight disease and on the regulatory function of RAVL1 in rice sheath blight disease.

System Pharmacology-Based Strategy to Decode the Synergistic Mechanism of Zhi-zhu Wan for Functional Dyspepsia

Functional dyspepsia (FD) is a widely prevalent gastrointestinal disorder throughout the world, whereas the efficacy of current treatment in the Western countries is limited. As the symptom is equivalent to the traditional Chinese medicine (TCM) term "stuffiness and fullness," FD can be treated with Zhi-zhu Wan (ZZW) which is a kind of Chinese patent medicine. However, the "multi-component" and "multi-target" feature of Chinese patent medicine makes it challenge to elucidate the potential therapeutic mechanisms of ZZW on FD. Presently, a novel system pharmacology model including pharmacokinetic parameters, pharmacological data, and component contribution score (CS) is constructed to decipher the potential therapeutic mechanism of ZZW on FD. Finally, 61 components with favorable pharmacokinetic profiles and biological activities were obtained through ADME (absorption, distribution, metabolism, and excretion) screening in silico. The related targets of these components are identified by component targeting process followed by GO analysis and pathway enrichment analysis. And systematic analysis found that through acting on the target related to inflammation, gastrointestinal peristalsis, and mental disorder, ZZW plays a synergistic and complementary effect on FD at the pathway level. Furthermore, the component CS showed that 29 components contributed 90.18% of the total CS values of ZZW for the FD treatment, which suggested that the effective therapeutic effects of ZZW for FD are derived from all active components, not a few components. This study proposes the system pharmacology method and discovers the potent combination therapeutic mechanisms of ZZW for FD. This strategy will provide a reference method for other TCM mechanism research.

Altered DNA Methylation of Long Noncoding RNA uc.167 Inhibits Cell Differentiation in Heart Development

In previous studies, we have demonstrated the function of uc.167 in the heart development. DNA methylation plays a crucial role in regulating the expression of developmental genes during embryonic development. In this study, the methylomic landscape was investigated in order to identify the DNA methylation alterations. Methylated DNA immunoprecipitation (MeDIP) was performed to examine the differences in methylation status of overexpressed uc.167 in P19 cells. GO and KEGG pathway analyses of differentially methylated genes were also conducted. We found that the distribution of differentially methylated regions (DMRs) peaks in different components of genome was mainly located in intergenic regions and intron. The biological process associated with uc.167 was focal adhesion and Rap1 signaling pathway. MEF2C was significantly decreased in uc.167 overexpressed group, suggesting that uc.167 may influence the P19 differentiation through MEF2C reduction. Taken together, our findings revealed that the effect of uc.167 on P19 differentiation may be attributed to the altered methylation of specific genes.

Dynamic transcriptome profiling towards understanding the morphogenesis and development of diverse feather in domestic duck

Background

Feathers with complex and fine structure are hallmark avian integument appendages, which have contributed significantly to the survival and breeding for birds. Here, we aimed to explore the differentiation, morphogenesis and development of diverse feathers in the domestic duck.

Results

Transcriptome profiles of skin owing feather follicle from two body parts at three physiological stages were constructed to understand the molecular network and excavate the candidate genes associated with the development of plumulaceous and flight feather structures. The venn analysis of differentially expressed genes (DEGs) between abdomen and wing skin tissues at three developmental stages showed that 38 genes owing identical differentially expression pattern. Together, our data suggest that feather morphological and structural diversity can be possibly related to the homeobox proteins. The key series-clusters, many candidate biological processes and genes were identified for the morphogenesis, growth and development of two feather types. Through comparing the results of developmental transcriptomes from plumulaceous and flight feather, we found that DEGs belonging to the family of WNT, FGF and BMP have certain differences; even the consistent DEGs of skin and feather follicle transcriptomes from abdomen and wing have the different expression patterns.

Conclusions

Overall, this study detected many functional genes and showed differences in the molecular mechanisms of diverse feather developments. The findings in WNT, FGF and BMP, which were consistent with biological experiments, showed more possible complex modulations. A correlative role of HOX genes was also suggested but future biological verification experiments are required. This work provided valuable information for subsequent research on the morphogenesis of feathers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4778-7) contains supplementary material, which is available to authorized users.

Combined flow cytometry and high-throughput image analysis for the study of essential genes in Caenorhabditis elegans

Background

Advances in automated image-based microscopy platforms coupled with high-throughput liquid workflows have facilitated the design of large-scale screens utilising multicellular model organisms such as Caenorhabditis elegans to identify genetic interactions, therapeutic drugs or disease modifiers. However, the analysis of essential genes has lagged behind because lethal or sterile mutations pose a bottleneck for high-throughput approaches, and a systematic way to analyse genetic interactions of essential genes in multicellular organisms has been lacking.

Results

In C. elegans, non-conditional lethal mutations can be maintained in heterozygosity using chromosome balancers, commonly expressing green fluorescent protein (GFP) in the pharynx. However, gene expression or function is typically monitored by the use of fluorescent reporters marked with the same fluorophore, presenting a challenge to sort worm populations of interest, particularly at early larval stages. Here, we develop a sorting strategy capable of selecting homozygous mutants carrying a GFP stress reporter from GFP-balanced animals at the second larval stage. Because sorting is not completely error-free, we develop an automated high-throughput image analysis protocol that identifies and discards animals carrying the chromosome balancer. We demonstrate the experimental usefulness of combining sorting of homozygous lethal mutants and automated image analysis in a functional genomic RNA interference (RNAi) screen for genes that genetically interact with mitochondrial prohibitin (PHB). Lack of PHB results in embryonic lethality, while homozygous PHB deletion mutants develop into sterile adults due to maternal contribution and strongly induce the mitochondrial unfolded protein response (UPR^mt). In a chromosome-wide RNAi screen for C. elegans genes having human orthologues, we uncover both known and new PHB genetic interactors affecting the UPR^mt and growth.

Conclusions

The method presented here allows the study of balanced lethal mutations in a high-throughput manner. It can be easily adapted depending on the user’s requirements and should serve as a useful resource for the C. elegans community for probing new biological aspects of essential nematode genes as well as the generation of more comprehensive genetic networks.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0496-5) contains supplementary material, which is available to authorized users.

A study of the key genes and inflammatory signaling pathways involved in HLA-B27-associated acute anterior uveitis families

The present study was conducted to investigate the key genes and the inflammatory signaling pathways involved in HLA-B27-associated acute anterior uveitis (AAU) families. Four families with HLA-B27-positive AAU patients and their HLA-B27-positive blood relatives were included in the study. Peripheral blood monocytes were isolated from the subjects and stimulated by lipopolysaccharides (LPS). Gene expression microarrays were used to identify the differentially expressed genes (DEGs), and the DEGs were analyzed by a range of bioinformatics-based techniques, including Gene Ontology (GO), Pathway analysis, Signal-Net analysis and Gene Relation Network (Gene-Rel-Net). Finally, ELISA was used to quantify cytokines in the supernatant. The gene expression microarrays identified 801 DEGs, including 349 upregulated and 452 downregulated genes. The GO analysis revealed several important functions, including metabolic, immune and inflammatory responses. The pathway analysis highlighted the enhanced activity of Staphylococcus aureus infection, chemokine and metabolic signaling pathways, as well as cytokine-to-cytokine receptor interactions. A total of 18 DEGs that were found to play critical roles by Signal-Net and Gene-Rel-Net and verified by quantitative polymerase chain reaction analysis were identified as key genes. In conclusion, monocytes from the AUU patients were more sensitive and exhibited a more prominent inflammatory response to stimulation by LPS compared with monocytes from healthy HLA-B27-positive blood relatives. These characterized DEGs may provide new evidence for the pathogenesis of AAU and help identify new therapeutic targets.

Genome-wide analysis of DNA methylation in human peripheral leukocytes identifies potential biomarkers of nonalcoholic fatty liver disease

The aim of the present study was to uncover the role of leukocytic DNA methylation in the evaluation of nonalcoholic fatty liver disease (NAFLD). Patients with biopsy-proven NAFLD (n=35) and normal controls (n=30) were recruited from Chinese Han population. Their DNA methylation in peripheral leukocytes was subjected to genome-wide profiling. The association between differential methylation of CpG sites and NAFLD was further investigated on the basis of histopathological classification, bioinformatics, and pyrosequencing. A panel of 863 differentially methylated CpG sites dominated by global hypomethylation, characterized the NAFLD patients. Hypomethylated CpG sites of Acyl-CoA synthetase long-chain family member 4 (ACSL4) (cg15536552) and carnitine palmitoyltransferase 1C (CPT1C) (cg21604803) associated with the increased risk of NAFLD [cg15536552, odds ratio (OR): 11.44, 95% confidence interval (CI): 1.04‑125.37, P=0.046; cg21604803, OR: 6.57, 95% CI: 1.02-42.15, P=0.047] at cut-off β-values of <3.36 (ACSL4 cg15536552) and <3.54 (CPT1C cg21604803), respectively, after the adjustment of age, sex, body mass index (BMI) and homeostasis model assessment of insulin resistant (HOMA-IR). Their methylation levels also served as biomarkers of NAFLD (ACSL4 cg15536552, AUC: 0.80, 95% CI: 0.62-0.98, P=0.009; CPT1C cg21604803, AUC: 0.78, 95% CI: 0.65-0.91, P=0.001). Pathologically, lowered methylation level (β-values <3.26) of ACSL4 (cg15536552) conferred susceptibility to nonalcoholic steatohepatitis (NASH). Taken together, genome-wide hypomethylation of peripheral leukocytes may differentiate NAFLD patients from normal controls. The leukocytic hypomethylated ACSL4 (cg15536552) was suggested to be a biomarker for the pathological characteristics of NAFLD.