Search

Abstract While prostaglandin E2 (PGE2) is produced in human tumor microenvironment (TME), its role therein remains poorly understood. Here, we examine this issue by comparative single-cell RNA sequencing of immune cells infiltrating human cancers and syngeneic tumors in female mice. PGE receptors EP4 and EP2 are expressed in lymphocytes and myeloid cells, and their expression is associated with the downregulation of oxidative phosphorylation (OXPHOS) and MYC targets, glycolysis and ribosomal proteins (RPs). Mechanistically, CD8+ T cells express EP4 and EP2 upon TCR activation, and PGE2 blocks IL-2-STAT5 signaling by downregulating Il2ra, which downregulates c-Myc and PGC-1 to decrease OXPHOS, glycolysis, and RPs, impairing migration, expansion, survival, and antitumor activity. Similarly, EP4 and EP2 are induced upon macrophage activation, and PGE2 downregulates c-Myc and OXPHOS in M1-like macrophages. These results suggest that PGE2-EP4/EP2 signaling impairs both adaptive and innate immunity in TME by hampering bioenergetics and ribosome biogenesis of tumor-infiltrating immune cells.

The roles of antibodies secreted by subsets of B cells in dengue virus (DENV) infection have been extensively studied, yet, the contribution of tissue-homing B cells to antiviral immunity remains unclear. In this study, we performed a comprehensive analysis of B cell subpopulations in peripheral blood samples from DENV-infected patients using single-cell RNA-sequencing (scRNA-seq) datasets and flow cytometry. We showed that plasma cells (PCs) and plasmablasts (PBs) were the predominant B cell populations during the acute phase of secondary natural DENV infection, but not in convalescent phase nor in healthy controls. Interestingly, these cells expressed proliferation, adhesion, and tissue-homing genes, including SELPLG, a homing marker of the skin, the initial infected site of DENV. Flow cytometry analysis confirmed a significant upregulation of cell surface expression of a cutaneous lymphocyte-associated antigen (CLA) encoded by SELPLG in PCs and PBs, compared to naive and memory B cells from the same patients. The analysis of an independent single-cell B-cell receptor sequencing (scBCR-seq) dataset of DENV-infected patients revealed that the peripheral blood PCs and PBs exhibited the highest clonal expansion in secondary DENV infection compared to other B cell subsets. These clonally expanded cells also expressed the highest levels of tissue-homing genes, including SELPLG. In addition, by utilizing a public scRNA-seq dataset of SARS-CoV2 infection, we demonstrated the upregulation of several tissue-homing genes in PCs and PBs. Our study provides evidence for the potential roles of tissue-homing B cell subsets in the context of immune responses against viral infections in humans.

Abstract The binding of transcription factors at proximal promoters and distal enhancers is central to gene regulation. Identifying regulatory motifs and quantifying their impact on expression remains challenging. Using a convolutional neural network trained on single-cell data, we infer putative regulatory motifs and cell type-specific importance. Our model, scover, explains 29% of the variance in gene expression in multiple mouse tissues. Applying scover to distal enhancers identified using scATAC-seq from the developing human brain, we identify cell type-specific motif activities in distal enhancers. Scover can identify regulatory motifs and their importance from single-cell data where all parameters and outputs are easily interpretable.

Soil salinity is a complex abiotic stress that involves several biological pathways. Hence, focusing on a specific or a few salt-tolerant phenotypes is unlikely to provide comprehensive insights into the intricate and interwinding mechanisms that regulate salt responsiveness. In this study, we develop a heuristic framework for systematically integrating and comprehensively evaluating quantitative trait loci (QTL) analyses from multiple stress-related traits obtained by different studies. Making use of a combined set of 46 salinity-related traits from three independent studies that were based on the same chromosome segment substitution line (CSSL) population of rice (Oryza sativa), we demonstrate how our approach can address technical biases and limitations from different QTL studies and calling methods. This allows us to compile a comprehensive list of trait-specific and multi-trait QTLs, as well as salinity-related candidate genes. In doing so, we discover several novel relationships between traits that demonstrate similar trends of phenotype scores across the CSSLs, as well as the similarities between genomic locations that the traits were mapped to. Finally, we experimentally validate our findings by expression analyses and functional validations of several selected candidate genes from multiple pathways in rice and Arabidopsis orthologous genes, including OsKS7 (ENT-KAURENE SYNTHASE 7), OsNUC1 (NUCLEOLIN 1) and OsFRO1 (FERRIC REDUCTASE OXIDASE 1) to name a few. This work not only introduces a novel approach for conducting comparative analyses of multiple QTLs, but also provides a list of candidate genes and testable hypotheses for salinity-related mechanisms across several biological pathways.

Rice is a staple food for more than half of the world’s population. However, the pervasive problem of salinity is severely undermining rice production, especially in coastal and low-lying areas where soil salinization is widespread. This stress, exacerbated by climate change, necessitates the development of salt-tolerant rice varieties to ensure food security. In this study, an F2:3 population (n = 454) from a cross of KDML105 and its chromosome segment substitution line (CSSL) was used to identify genomic regions associated with salt tolerance at the seedling stage. Using the QTL-seq approach, a QTL significantly associated with salt tolerance was identified on chromosome 1. Annotation of candidate genes in this region revealed the potential regulators of salt tolerance, including MIKC-type MADS domain proteins, calmodulin-binding transcription factors, and NB-ARC domain-containing proteins. These and other identified genes provide insights into the genetic basis of salt tolerance. This study underscores the importance of using advanced genomics tools and CSSL populations in the study of complex traits such as salt tolerance in rice. Several candidate genes identified in this study could be used in further studies on molecular or physiological mechanisms related to the salt response and tolerance mechanism in rice. Additionally, these genes could also be utilized in plant breeding programs for salt tolerance.

Summary: Here, we present a computational approach for investigating highly variable genes (HVGs) associated with biological pathways of interest, across multiple time points and cell types in single-cell RNA-sequencing (scRNA-seq) data. Using public dengue virus and COVID-19 datasets, we describe steps for using the framework to characterize the dynamic expression levels of HVGs related to common and cell-type-specific biological pathways over multiple immune cell types.For complete details on the use and execution of this protocol, please refer to Arora et al.1 : Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Summary: Effective clinical management of acute dengue virus (DENV) infection relies on the timing of suitable treatments during the disease progression. We analyzed single-cell transcriptomic profiles of the peripheral blood mononuclear cell samples from two DENV patients, collected daily during acute phase and also at convalescence. Key immune cell types demonstrated different dynamic responses over the course of the infection. On the day before defervescence (Day −1), we observed the peak expression of several prominent genes in the adaptive immunological pathways. We also characterized unique effector T cell clusters that expressed skin-homing signature genes at Day −1, whereas upregulation of skin and gut homing genes was also observed in plasma cells and plasmablasts during the febrile period. This work provides an overview of unique molecular dynamics that signify the entry of the critical phase, and the findings could improve the patient management of DENV infection.

A comprehensive understanding of dental pulp cellular compositions and their molecular responses to infection are crucial for the advancement of regenerative dentistry. Here, we presented a pilot study of single-cell transcriptomic profiles of 6,810 pulpal cells isolated from a sound human maxillary third molar and three carious teeth with enamel and deep dental caries. We observed altered immune cell compositions of the dental pulp in deep, but not enamel ones. Differential expression analysis revealed up-regulation of several pro-inflammatory, anti-inflammatory, and mineralization-related genes in the immune and stromal cells of the deep dental caries. Making use of an algorithm for predicting cell-to-cell interactions from single-cell transcriptomic profiles, we showed an increase in cell-cell interactions between B cells, plasma cells and macrophages, and other cell types in deep dental caries, including those between TIMP1 (odontoblasts)—CD63 (myeloid cells), and CCL2 (macrophages)—ACKR1 (endothelial cells). Collectively, our work highlighted the single-cell level gene regulations and intercellular interactions in the dental pulps in health and disease.

Abstract Molecular mechanisms of how constant temperatures affect flowering time have been largely characterized in the model plant Arabidopsis thaliana; however, the effect of natural daily variable temperature outside laboratories is only partly explored. Several flowering genes have been shown to play important roles in temperature responses, including PHYTOCHROME‐INTERACTING FACTOR 4 (PIF4) and FLOWERING LOCUS C (FLC), the two genes encoding for the transcription factors (TFs) that act antagonistically to regulate flowering time by activating and repressing floral integrator FLOWERING LOCUS T (FT), respectively. In this study, we have taken a multidisciplinary approach to explore the contribution of PIF4 to the early flowering observed in the daily variable temperature (VAR) and to broaden its transcriptional network using publicly available transcriptomic data. We observed early flowering in the natural accessions Col‐0, C24 and their late flowering hybrid C24xCol grown under VAR, as compared with a constant temperature (CON). The loss‐of‐function mutation of PIF4 exhibits later flowering in VAR in both the Col‐0 parent and the C24xCol hybrid, suggesting that PIF4, at least in part, contributes to acceleration of flowering in the VAR condition. To investigate the interplay between PIF4 and its flowering regulator counterparts, FLC and FT, we performed transcriptional analyses and found that VAR increased PIF4 transcription at the end of the day when temperature peaked at 32°C, when FT transcription was also elevated. On the other hand, we observed a decrease in FLC transcription in the 4‐week‐old plants grown in VAR, as well as in the plants with PIF4 overexpression grown in CON. These results raise a possibility that PIF4 might also regulate FT indirectly through the repression of FLC, in addition to the well‐characterized direct control of PIF4 over FT. To further expand our view on the PIF4‐orientated flowering gene network in response to temperature changes, we have constructed a coexpression–transcriptional regulatory network by combining publicly available transcriptomic data and gene regulatory interactions of PIF4 and its closely related flowering genes, PIF5, FLC, and ELF3. The network model reveals conserved and tissue‐specific regulatory functions, which are useful for confirming as well as predicting the functions and regulatory interactions between these key flowering genes.

BackgroundDengue virus (DENV) infection has a global impact on public health. The clinical outcomes (of DENV) can vary from a flu-like illness called dengue fever (DF), to a more severe form, known as dengue hemorrhagic fever (DHF). The underlying innate immune mechanisms leading to protective or detrimental outcomes have not been fully elucidated. Helper innate lymphoid cells (hILCs), an innate lymphocyte recently discovered, functionally resemble T-helper cells and are important in inflammation and homeostasis. However, the role of hILCs in DENV infection had been unexplored.MethodsWe performed flow cytometry to investigate the frequency and phenotype of hILCs in peripheral blood mononuclear cells from DENV-infected patients of different disease severities (DF and DHF), and at different phases (febrile and convalescence) of infection. Intracellular cytokine staining of hILCs from DF and DHF were also evaluated by flow cytometry after ex vivo stimulation. Further, the hILCs were sorted and subjected to transcriptome analysis using RNA sequencing. Differential gene expression analysis was performed to compare the febrile and convalescent phase samples in DF and DHF. Selected differentially expressed genes were then validated by quantitative PCR.ResultsPhenotypic analysis showed marked activation of all three hILC subsets during the febrile phase as shown by higher CD69 expression when compared to paired convalescent samples, although the frequency of hILCs remained unchanged. Upon ex vivo stimulation, hILCs from febrile phase DHF produced significantly higher IFN-γ and IL-4 when compared to those of DF. Transcriptomic analysis showed unique hILCs gene expression in DF and DHF, suggesting that divergent functions of hILCs may be associated with different disease severities. Differential gene expression analysis indicated that hILCs function both in cytokine secretion and cytotoxicity during the febrile phase of DENV infection.ConclusionsHelper ILCs are activated in the febrile phase of DENV infection and display unique transcriptomic changes as well as cytokine production that correlate with severity. Targeting hILCs during early innate response to DENV might help shape subsequent immune responses and potentially lessen the disease severity in the future.

Spinal muscular atrophy (SMA) is a severe, debilitating neuromuscular condition characterised by loss of motor neurons and progressive muscle wasting. SMA is caused by a loss of expression of SMN1 that encodes the survival motor neuron (SMN) protein necessary for the survival of motor neurons. Restoration of SMN expression through increased inclusion of SMN2 exon 7 is known to ameliorate symptoms in SMA patients. As a consequence, regulation of pre-mRNA splicing of SMN2 could provide a potential molecular therapy for SMA. In this study, we explored if splice switching antisense oligonucleotides could redirect the splicing repressor hnRNPA1 to the hnRNPA1b isoform and restore SMN expression in fibroblasts from a type I SMA patient. Antisense oligonucleotides (AOs) were designed to promote exon 7b retention in the mature mRNA and induce the hnRNPA1b isoform. RT-PCR and western blot analysis were used to assess and monitor the efficiency of different AO combinations. A combination of AOs targeting multiple silencing motifs in hnRNPA1 pre-mRNA led to robust hnRNPA1b induction, which, in turn, significantly increased expression of full-length SMN (FL-SMN) protein. A combination of PMOs targeting the same motifs also strongly induced hnRNPA1b isoform, but surprisingly SMN2 exon 5 skipping was detected, and the PMO cocktail did not lead to a significant increase in expression of FL-SMN protein. We further performed RNA sequencing to assess the genome-wide effects of hnRNPA1b induction. Some 3244 genes were differentially expressed between the hnRNPA1b-induced and untreated SMA fibroblasts, which are functionally enriched in cell cycle and chromosome segregation processes. RT-PCR analysis demonstrated that expression of the master regulator of these enrichment pathways, MYBL2 and FOXM1B, were reduced in response to PMO treatment. These findings suggested that induction of hnRNPA1b can promote SMN protein expression, but not at sufficient levels to be clinically relevant.

Summary: Temperature is a key environmental variable influencing plant growth and survival. Protection against high temperature stress in eukaryotes is coordinated by heat shock factors (HSFs), transcription factors that activate the expression of protective chaperones such as HEAT SHOCK PROTEIN 70 (HSP70); however, the pathway by which temperature is sensed and integrated with other environmental signals into adaptive responses is not well understood. Plants are exposed to considerable diurnal variation in temperature, and we have found that there is diurnal variation in thermotolerance in Arabidopsis thaliana, with maximal thermotolerance coinciding with higher HSP70 expression during the day. In a forward genetic screen, we identified a key role for the chloroplast in controlling this response, suggesting that light-induced chloroplast signaling plays a key role. Consistent with this, we are able to globally activate binding of HSFA1a to its targets by altering redox status in planta independently of a heat shock.

Abstract With the wealth of data accumulated from completely sequenced genomes and other high-throughput experiments, global studies of biological systems, by simultaneously investigating multiple biological entities (e.g. genes, transcripts, proteins), has become a routine. Network representation is frequently used to capture the presence of these molecules as well as their relationship. Network biology has been widely used in molecular biology and genetics, where several network properties have been shown to be functionally important. Here, we discuss how such methodology can be useful to translational biomedical research, where scientists traditionally focus on one or a small set of genes, diseases, and drug candidates at any one time. We first give an overview of network representation frequently used in biology: what nodes and edges represent, and review its application in preclinical research to date. Using cancer as an example, we review how network biology can facilitate system-wide approaches to identify targeted small molecule inhibitors. These types of inhibitors have the potential to be more specific, resulting in high efficacy treatments with less side effects, compared to the conventional treatments such as chemotherapy. Global analysis may provide better insight into the overall picture of human diseases, as well as identify previously overlooked problems, leading to rapid advances in medicine. From the clinicians’ point of view, it is necessary to bridge the gap between theoretical network biology and practical biomedical research, in order to improve the diagnosis, prevention, and treatment of the world’s major diseases.

Altered cellular metabolism is a fundamental adaptation of cancer during rapid proliferation as a result of growth factor overstimulation. We review different pathways involving metabolic alterations in cancers including aerobic glycolysis, pentose phosphate pathway, de novo fatty acid synthesis, and serine and glycine metabolism. Although oncoproteins, c-MYC, HIF1α and p53 are the major drivers of this metabolic reprogramming, post-transcriptional regulation by microRNAs (miR) also plays an important role in finely adjusting the requirement of the key metabolic enzymes underlying this metabolic reprogramming. We also combine the literature data on the miRNAs that potentially regulate 40 metabolic enzymes responsible for metabolic reprogramming in cancers, with additional miRs from computational prediction. Our analyses show that: (1) a metabolic enzyme is frequently regulated by multiple miRs, (2) confidence scores from prediction algorithms might be useful to help narrow down functional miR-mRNA interaction, which might be worth further experimental validation. By combining known and predicted interactions of oncogenic transcription factors (TFs) (c-MYC, HIF1α and p53), sterol regulatory element binding protein 1 (SREBP1), 40 metabolic enzymes, and regulatory miRs we have established one of the first reference maps for miRs and oncogenic TFs that regulate metabolic reprogramming in cancers. The combined network shows that glycolytic enzymes are linked to miRs via p53, c-MYC, HIF1α, whereas the genes in serine, glycine and one carbon metabolism are regulated via the c-MYC, as well as other regulatory organization that cannot be observed by investigating individual miRs, TFs, and target genes.

Understanding specific biological processes involving in salt tolerance mechanisms is important for improving traits conferring tolerance to salinity, one of the most important abiotic stresses in plants. Under drought and salinity stresses, plants share overlapping responsive mechanisms such as physiological changes and activation of signaling molecules, which induce and transmit signals through regulator genes in a regulatory network. In this study, two near isogenic lines of rice carrying chromosome segments of drought tolerance QTL on chromosome 8 from IR68586-F2-CA-31 (DH103) in the genetic background of sensitive cultivar “Khao Dawk Mali 105; KDML105” (designated as CSSL8-94 and CSSL8-95) were used to investigate physiological responses to salt stress [namely growth, Na+/K+ ratio, water status, osmotic adjustment, photosynthetic parameters, electrolyte leakage (EL), malondialdehyde (MDA), proline and sugar accumulations], compared with the standard salt tolerant (Pokkali; PK) and their recurrent parent (KDML105) rice cultivars. Physiological examination indicated that both CSSLs showed superior salt-tolerant level to KDML105. Our results suggested that salt tolerance ability of these CSSL lines may be resulted from high performance photosynthesis, better osmotic adjustment, and less oxidative stress damage under salt conditions. Moreover, to explore new candidate genes that might take part in salt tolerance mechanisms, we performed co-expression network analysis for genes identified in the CSSL rice, and found that Os08g419090, the gene involved with tetrapyrrole and porphyrin biosynthetic process (chlorophyll biosynthetic process), Os08g43230 and Os08g43440 (encoded TraB family protein and cytochrome P450, respectively) might have unprecedented roles in salt stress tolerance.

Intrahepatic cholangiocarcinoma (ICC) is the cancer of the intrahepatic bile ducts, and together with hepatocellular carcinoma (HCC), constitute the majority of primary liver cancers. ICC is a rare disorder as its overall incidence is < 1/100,000 in the United States and Europe. However, it shows much higher incidence in particular geographical regions, such as northeastern Thailand, where liver fluke infection is the most common risk factor of ICC. Since the early stages of ICC are often asymptomatic, the patients are usually diagnosed at advanced stages with no effective treatments available, leading to the high mortality rate. In addition, unclear genetic mechanisms, heterogeneous nature, and various etiologies complicate the development of new efficient treatments. Recently, a number of studies have employed high-throughput approaches, including next-generation sequencing and mass spectrometry, in order to understand ICC in different biological aspects. In general, the majority of recurrent genetic alterations identified in ICC are enriched in known tumor suppressor genes and oncogenes, such as mutations in TP53, KRAS, BAP1, ARID1A, IDH1, IDH2, and novel FGFR2 fusion genes. Yet, there are no major driver genes with immediate clinical solutions characterized. Interestingly, recent studies utilized multi-omics data to classify ICC into two main subgroups, one with immune response genes as the main driving factor, while another is enriched with driver mutations in the genes associated with epigenetic regulations, such as IDH1 and IDH2. The two subgroups also show different hypermethylation patterns in the promoter regions. Additionally, the immune response induced by host-pathogen interactions, i.e., liver fluke infection, may further stimulate tumor growth through alterations of the tumor microenvironment. For in-depth functional studies, although many ICC cell lines have been globally established, these homogeneous cell lines may not fully explain the highly heterogeneous genetic contents of this disorder. Therefore, the advent of patient-derived xenograft and 3D patient-derived organoids as new disease models together with the understanding of evolution and genetic alterations of tumor cells at the single-cell resolution will likely become the main focus to fill the current translational research gaps of ICC in the future.

Abstract The expression level of a gene is often used as a proxy for determining whether the protein or RNA product is functional in a cell or tissue. Therefore, it is of fundamental importance to understand the global distribution of gene expression levels, and to be able to interpret it mechanistically and functionally. Here we use RNA sequencing (RNA‐seq) of mouse Th2 cells, coupled with a range of other techniques, to show that all genes can be separated, based on their expression abundance, into two distinct groups: one group comprised of lowly expressed and putatively non‐functional mRNAs, and the other of highly expressed mRNAs with active chromatin marks at their promoters. These observations are confirmed in many other microarray and RNA‐seq data sets of metazoan cell types.

Copy number variation (CNV) is a major genetic polymorphism contributing to genetic diversity and human evolution. Clinical application of CNVs for diagnostic purposes largely depends on sufficient population CNV data for accurate interpretation. CNVs from general population in currently available databases help classify CNVs of uncertain clinical significance, and benign CNVs. Earlier studies of CNV distribution in several populations worldwide showed that a significant fraction of CNVs are population specific. In this study, we characterized and analyzed CNVs in 3,017 unrelated Thai individuals genotyped with the Illumina Human610, Illumina HumanOmniexpress, or Illumina HapMap550v3 platform. We employed hidden Markov model and circular binary segmentation methods to identify CNVs, extracted 23,458 CNVs consistently identified by both algorithms, and cataloged these high confident CNVs into our publicly available Thai CNV database. Analysis of CNVs in the Thai population identified a median of eight autosomal CNVs per individual. Most CNVs (96.73%) did not overlap with any known chromosomal imbalance syndromes documented in the DECIPHER database. When compared with CNVs in the 11 HapMap3 populations, CNVs found in the Thai population shared several characteristics with CNVs characterized in HapMap3. Common CNVs in Thais had similar frequencies to those in the HapMap3 populations, and all high frequency CNVs (>20%) found in Thai individuals could also be identified in HapMap3. The majorities of CNVs discovered in the Thai population, however, were of low frequency, or uniquely identified in Thais. When performing hierarchical clustering using CNV frequencies, the CNV data were clustered into Africans, Europeans, and Asians, in line with the clustering performed with single nucleotide polymorphism (SNP) data. As CNV data are specific to origin of population, our population-specific reference database will serve as a valuable addition to the existing resources for the investigation of clinical significance of CNVs in Thais and related ethnicities.