ViSEAGO: a Bioconductor package for clustering biological functions using Gene Ontology and semantic similarity

Implication of digestive functions and microbiota in the establishment of muscle glycogen differences between divergent lines for ultimate pH

Both the quality of chicken meat and the quality of chicks are influenced by the level of breast muscle glycogen reserves. In order to study the role of digestive metabolism in establishing this muscular phenotype, we compared two divergent chicken lines for the ultimate pH (pHu) of the breast meat, a proxy for glycogen reserves. Males aged 4 weeks had twice the breast muscle glycogen content in the pHu- line (low pHu) than in the pHu +  line (high pHu). The increase in glycogen reserves (pHu-) was associated with a higher relative weight of the proventriculus and gizzard, as well as better apparent ileal digestibility of nitrogen and calcium. The diversity of the cecal microbiota was comparable, but three bacterial genera (Lachnospira, Lachnospiraceae UCG-010, Caproiciproducens) varied between the lines. The differences observed could lead to down-regulation of carbon fixation in prokaryotes and of the citrate cycle in the pHu + line. RNA-seq analysis of the jejunum, the major site of nutrient absorption, revealed 149 genes differentially expressed (DE) between the lines, including several genes linked to immunity, hormonal response and circadian rhythms that are less expressed in pHu + animals. Others involved in cell migration and proliferation, and more generally tissue morphogenesis, also differed between the lines. Among the DE genes, several co-localized with Quantitative Trait Loci (QTL) controlling pHu and selection signatures identified in the divergent lines, such as the gene coding for ghrelin, a hormone regulating appetite.

RNA-seq dataset of the chorioallantoic membrane of male and female chicken embryos, after 11 and 15 days of incubation

The chicken chorioallantoic membrane (CAM) is an extraembryonic structure that exhibits many vital functions to support the development of the chicken embryo (gaseous exchange, innate defence, calcium transport from the eggshell to the embryo skeleton, homeostasis). Developing from day 6 of incubation, the CAM progressively differentiates into three functional layers (the chorionic epithelium in contact with the inner eggshell, the highly vascularized mesoderm, and the allantoic epithelium), between 11 and 15 days of incubation. This article describes the RNASeq dataset and the analyses performed on total CAMs collected from male and female embryos after 11 and 15 days of incubation. The datasets are available at the NCBI Gene Expression Omnibus (GEO) repository (http://www.ncbi.nlm.nih.gov/geo) using GSE199780 as the accession number. The statistical analysis of the data allowed identifying genes differentially expressed depending on the sex of the embryo at two time points of CAM differentiation. Knowing that the CAM is widely used as a model to study tumour growth, metastasis or wound healing, the resulting analysis highlights the necessity to include this sex variable in experimental assays to avoid any bias of interpretation. Indeed, the functional annotation of genes that are differentially expressed between male and female CAMs revealed an enrichment of activities and functions related to lipid metabolism, bone formation, and morphogenesis suggesting that the response of the CAM to external and experimental stimuli might be different depending on the sex of the embryo.

Decoded cardiopoietic cell secretome linkage to heart repair biosignature

Cardiopoiesis-primed human stem cells exert sustained benefit in treating heart failure despite limited retention following myocardial delivery. To assess potential paracrine contribution, the secretome of cardiopoiesis conditioned versus naïve human mesenchymal stromal cells was decoded by directed proteomics augmented with machine learning and systems interrogation. Cardiopoiesis doubled cellular protein output generating a distinct secretome that segregated the conditioned state. Altering the expression of 1035 secreted proteins, cardiopoiesis reshaped the secretome across functional classes. The resolved differential cardiopoietic secretome was enriched in mesoderm development and cardiac progenitor signaling processes, yielding a cardiovasculogenic profile bolstered by upregulated cardiogenic proteins. In tandem, cardiopoiesis enhanced the secretion of immunomodulatory proteins associated with cytokine signaling, leukocyte migration, and chemotaxis. Network analysis integrated the differential secretome within an interactome of 1745 molecules featuring prioritized regenerative processes. Secretome contribution to the repair signature of cardiopoietic cell-treated infarcted hearts was assessed in a murine coronary ligation model. Intramyocardial delivery of cardiopoietic cells improved the performance of failing hearts, with undirected proteomics revealing 50 myocardial proteins responsive to cell therapy. Pathway analysis linked the secretome to cardiac proteome remodeling, pinpointing 17 cardiopoiesis-upregulated secretome proteins directly upstream of 44% of the cell therapy-responsive cardiac proteome. Knockout, in silico, of this 22-protein secretome-dependent myocardial ensemble eliminated indices of the repair signature. Accordingly, in vivo, cell therapy rendered the secretome-dependent myocardial proteome of an infarcted heart indiscernible from healthy counterparts. Thus, the secretagogue effect of cardiopoiesis transforms the human stem cell secretome, endows regenerative competency, and upregulates candidate paracrine effectors of cell therapy-mediated molecular restitution.

Impact of tRNA-induced proline-to-serine mistranslation on the transcriptome of Drosophila melanogaster

Mistranslation is the misincorporation of an amino acid into a polypeptide. Mistranslation has diverse effects on multicellular eukaryotes and is implicated in several human diseases. In Drosophila melanogaster, a serine transfer RNA (tRNA) that misincorporates serine at proline codons (P→S) affects male and female flies differently. The mechanisms behind this discrepancy are currently unknown. Here, we compare the transcriptional response of male and female flies to P→S mistranslation to identify genes and cellular processes that underlie sex-specific differences. Both males and females downregulate genes associated with various metabolic processes in response to P→S mistranslation. Males downregulate genes associated with extracellular matrix organization and response to negative stimuli such as wounding, whereas females downregulate aerobic respiration and ATP synthesis genes. Both sexes upregulate genes associated with gametogenesis, but females also upregulate cell cycle and DNA repair genes. These observed differences in the transcriptional response of male and female flies to P→S mistranslation have important implications for the sex-specific impact of mistranslation on disease and tRNA therapeutics.

Jejunal transcriptomic profiling of carnosine synthesis precursor-related genes and pathways in slow-growing Korat chicken

Carnosine is a physiologically important molecule in normal human body functions. Chicken meat is an excellent source of carnosine; especially slow-growing Korat chicken (KR) females have a high carnosine content in their meat. The carnosine content of chicken meat can be increased by dietary supplementation of β-alanine (βA) and L-histidine (L-His). Our objective was to reveal the pathways and genes through jejunal transcriptomic profiling related to βA and L-His absorption and transportation. We collected whole jejunum samples from 5 control and 5 experimental KR chicken, fed with 1% βA and 0.5% L-His supplementation. A total of 407 differentially expressed genes (P < 0.05, log2 fold change ≥2) were identified, 272 of which were down-regulated and 135 up-regulated in the group with dietary supplementation compared to the control group. Based on the integrated analysis of the protein-protein interaction network and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway maps, 87 gene ontology terms were identified and 6 KEGG pathways were significantly (P < 0.05) enriched in the jejunum. The analyses revealed 6 key genes, KCND3, OPRM1, CCK, GCG, TRH, and GABBR2, that are related to neuroactive ligand-receptor interaction and the calcium signaling pathway. These findings give insight regarding the molecular mechanism related to carnosine precursor absorption and transportation in the jejunum and help to identify useful molecular markers for improving the carnosine content in slow-growing KR chicken meat.

The activity of early-life gene regulatory elements is hijacked in aging through pervasive AP-1-linked chromatin opening

A mechanistic connection between aging and development is largely unexplored. Through profiling age-related chromatin and transcriptional changes across 22 murine cell types, analyzed alongside previous mouse and human organismal maturation datasets, we uncovered a transcription factor binding site (TFBS) signature common to both processes. Early-life candidate cis-regulatory elements (cCREs), progressively losing accessibility during maturation and aging, are enriched for cell-type identity TFBSs. Conversely, cCREs gaining accessibility throughout life have a lower abundance of cell identity TFBSs but elevated activator protein 1 (AP-1) levels. We implicate TF redistribution toward these AP-1 TFBS-rich cCREs, in synergy with mild downregulation of cell identity TFs, as driving early-life cCRE accessibility loss and altering developmental and metabolic gene expression. Such remodeling can be triggered by elevating AP-1 or depleting repressive H3K27me3. We propose that AP-1-linked chromatin opening drives organismal maturation by disrupting cell identity TFBS-rich cCREs, thereby reprogramming transcriptome and cell function, a mechanism hijacked in aging through ongoing chromatin opening.

Heat tolerance, oxidative stress response tuning and robust gene activation in early-stage Drosophila melanogaster embryos

In organisms with complex life cycles, life stages that are most susceptible to environmental stress may determine species persistence in the face of climate change. Early embryos of Drosophila melanogaster are particularly sensitive to acute heat stress, yet tropical embryos have higher heat tolerance than temperate embryos, suggesting adaptive variation in embryonic heat tolerance. We compared transcriptomic responses to heat stress among tropical and temperate embryos to elucidate the gene regulatory basis of divergence in embryonic heat tolerance. The transcriptomes of tropical and temperate embryos differed in both constitutive and heat-stress-induced responses of the expression of relatively few genes, including genes involved in oxidative stress. Most of the transcriptomic response to heat stress was shared among all embryos. Embryos shifted the expression of thousands of genes, including increases in the expression of heat shock genes, suggesting robust zygotic gene activation and demonstrating that, contrary to previous reports, early embryos are not transcriptionally silent. The involvement of oxidative stress genes corroborates recent reports on the critical role of redox homeostasis in coordinating developmental transitions. By characterizing adaptive variation in the transcriptomic basis of embryonic heat tolerance, this study is a novel contribution to the literature on developmental physiology and developmental genetics.

Evolutionary and biomedical implications of sex differences in the primate brain transcriptome

Humans exhibit sex differences in the prevalence of many neurodevelopmental disorders and neurodegenerative diseases. Here, we generated one of the largest multi-brain-region bulk transcriptional datasets for the rhesus macaque and characterized sex-biased gene expression patterns to investigate the translatability of this species for sex-biased neurological conditions. We identify patterns similar to those in humans, which are associated with overlapping regulatory mechanisms, biological processes, and genes implicated in sex-biased human disorders, including autism. We also show that sex-biased genes exhibit greater genetic variance for expression and more tissue-specific expression patterns, which may facilitate rapid evolution of sex-biased genes. Our findings provide insights into the biological mechanisms underlying sex-biased disease and support the rhesus macaque model for the translational study of these conditions.

Sparse spectral graph analysis and its application to gastric cancer drug resistance-specific molecular interplays identification

Uncovering acquired drug resistance mechanisms has garnered considerable attention as drug resistance leads to treatment failure and death in patients with cancer. Although several bioinformatics studies developed various computational methodologies to uncover the drug resistance mechanisms in cancer chemotherapy, most studies were based on individual or differential gene expression analysis. However the single gene-based analysis is not enough, because perturbations in complex molecular networks are involved in anti-cancer drug resistance mechanisms. The main goal of this study is to reveal crucial molecular interplay that plays key roles in mechanism underlying acquired gastric cancer drug resistance. To uncover the mechanism and molecular characteristics of drug resistance, we propose a novel computational strategy that identified the differentially regulated gene networks. Our method measures dissimilarity of networks based on the eigenvalues of the Laplacian matrix. Especially, our strategy determined the networks' eigenstructure based on sparse eigen loadings, thus, the only crucial features to describe the graph structure are involved in the eigenanalysis without noise disturbance. We incorporated the network biology knowledge into eigenanalysis based on the network-constrained regularization. Therefore, we can achieve a biologically reliable interpretation of the differentially regulated gene network identification. Monte Carlo simulations show the outstanding performances of the proposed methodology for differentially regulated gene network identification. We applied our strategy to gastric cancer drug-resistant-specific molecular interplays and related markers. The identified drug resistance markers are verified through the literature. Our results suggest that the suppression and/or induction of COL4A1, PXDN and TGFBI and their molecular interplays enriched in the Extracellular-related pathways may provide crucial clues to enhance the chemosensitivity of gastric cancer. The developed strategy will be a useful tool to identify phenotype-specific molecular characteristics that can provide essential clues to uncover the complex cancer mechanism.

Integration mapping of cardiac fibroblast single-cell transcriptomes elucidates cellular principles of fibrosis in diverse pathologies

Single-cell technology has allowed researchers to probe tissue complexity and dynamics at unprecedented depth in health and disease. However, the generation of high-dimensionality single-cell atlases and virtual three-dimensional tissues requires integrated reference maps that harmonize disparate experimental designs, analytical pipelines, and taxonomies. Here, we present a comprehensive single-cell transcriptome integration map of cardiac fibrosis, which underpins pathophysiology in most cardiovascular diseases. Our findings reveal similarity between cardiac fibroblast (CF) identities and dynamics in ischemic versus pressure overload models of cardiomyopathy. We also describe timelines for commitment of activated CFs to proliferation and myofibrogenesis, profibrotic and antifibrotic polarization of myofibroblasts and matrifibrocytes, and CF conservation across mouse and human healthy and diseased hearts. These insights have the potential to inform knowledge-based therapies.

Orco-dependent survival of odorant receptor neurons in ants

Olfaction is essential for complex social behavior in insects. To discriminate complex social cues, ants evolved an expanded number of odorant receptor (Or) genes. Mutations in the obligate odorant co-receptor gene orco lead to the loss of ~80% of the antennal lobe glomeruli in the jumping ant Harpegnathos saltator. However, the cellular mechanism remains unclear. Here, we demonstrate massive apoptosis of odorant receptor neurons (ORNs) in the mid to late stages of pupal development, possibly due to ER stress in the absence of Orco. Further bulk and single-nucleus transcriptome analysis shows that, although most orco-expressing ORNs die in orco mutants, a small proportion of them survive: They express ionotropic receptor (Ir) genes that form IR complexes. In addition, we found that some Or genes are expressed in mechanosensory neurons and nonneuronal cells, possibly due to leaky regulation from nearby non-Or genes. Our findings provide a comprehensive overview of ORN development and Or expression in H. saltator.

Impact of tRNA-induced proline-to-serine mistranslation on the transcriptome of Drosophila melanogaster

Mistranslation is the misincorporation of an amino acid into a polypeptide. Mistranslation has diverse effects on multicellular eukaryotes and is implicated in several human diseases. In Drosophila melanogaster, a serine transfer RNA (tRNA) that misincorporates serine at proline codons (P→S) affects male and female flies differently. The mechanisms behind this discrepancy are currently unknown. Here, we compare the transcriptional response of male and female flies to P→S mistranslation to identify genes and cellular processes that underlie sex-specific differences. Both males and females downregulate genes associated with various metabolic processes in response to P→S mistranslation. Males downregulate genes associated with extracellular matrix organization and response to negative stimuli such as wounding, whereas females downregulate aerobic respiration and ATP synthesis genes. Both sexes upregulate genes associated with gametogenesis, but females also upregulate cell cycle and DNA repair genes. These observed differences in the transcriptional response of male and female flies to P→S mistranslation have important implications for the sex-specific impact of mistranslation on disease and tRNA therapeutics.

In-depth characterization of a selection of gut commensal bacteria reveals their functional capacities to metabolize dietary carbohydrates with prebiotic potential

The microbial utilization of dietary carbohydrates is closely linked to the pivotal role of the gut microbiome in human health. Inherent to the modulation of complex microbial communities, a prebiotic implies the selective utilization of a specific substrate, relying on the metabolic capacities of targeted microbes. In this study, we investigated the metabolic capacities of 17 commensal bacteria of the human gut microbiome toward dietary carbohydrates with prebiotic potential. First, in vitro experiments allowed the classification of bacterial growth and fermentation profiles in response to various carbon sources, including agave inulin, corn fiber, polydextrose, and citrus pectin. The influence of phylogenetic affiliation appeared to statistically outweigh carbon sources in determining the degree of carbohydrate utilization. Second, we narrowed our focus on six commensal bacteria representative of the Bacteroidetes and Firmicutes phyla to perform an untargeted high-resolution liquid chromatography-mass spectrometry metabolomic analysis: Bacteroides xylanisolvens, Bacteroides thetaiotaomicron, Bacteroides intestinalis, Subdoligranulum variabile, Roseburia intestinalis, and Eubacterium rectale exhibited distinct metabolomic profiles in response to different carbon sources. The relative abundance of bacterial metabolites was significantly influenced by dietary carbohydrates, with these effects being strain-specific and/or carbohydrate-specific. Particularly, the findings indicated an elevation in short-chain fatty acids and other metabolites, including succinate, gamma-aminobutyric acid, and nicotinic acid. These metabolites were associated with putative health benefits. Finally, an RNA-Seq transcriptomic approach provided deeper insights into the underlying mechanisms of carbohydrate metabolization. Restricting our focus on four commensal bacteria, including B. xylanisolvens, B. thetaiotaomicron, S. variabile, and R. intestinalis, carbon sources did significantly modulate the level of bacterial genes related to the enzymatic machinery involved in the metabolization of dietary carbohydrates. This study provides a holistic view of the molecular strategies induced during the dynamic interplay between dietary carbohydrates with prebiotic potential and gut commensal bacteria.

IMPORTANCE

This study explores at a molecular level the interactions between commensal health-relevant bacteria and dietary carbohydrates holding prebiotic potential. We showed that prebiotic breakdown involves the specific activation of gene expression related to carbohydrate metabolism. We also identified metabolites produced by each bacteria that are potentially related to our digestive health. The characterization of the functional activities of health-relevant bacteria toward prebiotic substances can yield a better application of prebiotics in clinical interventions and personalized nutrition. Overall, this study highlights the importance of identifying the impact of prebiotics at a low resolution of the gut microbiota to characterize the activities of targeted bacteria that can play a crucial role in our health.

Silencing Apoe with divalent-siRNAs improves amyloid burden and activates immune response pathways in Alzheimer's disease

INTRODUCTION

The most significant genetic risk factor for late-onset Alzheimer's disease (AD) is APOE4, with evidence for gain- and loss-of-function mechanisms. A clinical need remains for therapeutically relevant tools that potently modulate APOE expression.

METHODS

We optimized small interfering RNAs (di-siRNA, GalNAc) to potently silence brain or liver Apoe and evaluated the impact of each pool of Apoe on pathology.

RESULTS

In adult 5xFAD mice, siRNAs targeting CNS Apoe efficiently silenced Apoe expression and reduced amyloid burden without affecting systemic cholesterol, confirming that potent silencing of brain Apoe is sufficient to slow disease progression. Mechanistically, silencing Apoe reduced APOE-rich amyloid cores and activated immune system responses.

DISCUSSION

These results establish siRNA-based modulation of Apoe as a viable therapeutic approach, highlight immune activation as a key pathway affected by Apoe modulation, and provide the technology to further evaluate the impact of APOE silencing on neurodegeneration.

Licensing effects of inflammatory factors and TLR ligands on the regenerative capacity of adipose-derived mesenchymal stem cells

Introduction: Adipose tissue-derived mesenchymal stem cells are promising contributors to regenerative medicine, exhibiting the ability to regenerate tissues and modulate the immune system, which is particularly beneficial for addressing chronic inflammatory ulcers and wounds. Despite their inherent capabilities, research suggests that pretreatment amplifies therapeutic effectiveness. Methods: Our experimental design exposed adipose-derived mesenchymal stem cells to six inflammatory factors for 24 h. We subsequently evaluated gene expression and proteome profile alterations and observed the wound closure rate post-treatment. Results: Specific pretreatments, such as IL-1β, notably demonstrated an accelerated wound-healing process. Analysis of gene and protein expression profiles revealed alterations in pathways associated with tissue regeneration. Discussion: This suggests that licensed cells exhibit potentially higher therapeutic efficiency than untreated cells, shedding light on optimizing regenerative strategies using adipose tissue-derived stem cells.

Single-cell RNA sequencing unveils unique transcriptomic signatures of endothelial cells and role of ENO1 in response to disturbed flow

Flow patterns exert significant effects on vascular endothelial cells (ECs) to lead to the focal nature of atherosclerosis. Using a step flow chamber to investigate the effects of disturbed shear (DS) and pulsatile shear (PS) on ECs in the same flow channel, we conducted single-cell RNA sequencing analyses to explore the distinct transcriptomic profiles regulated by DS vs. PS. Integrated analysis identified eight cell clusters and demonstrated that DS induces EC transition from atheroprotective to proatherogenic phenotypes. Using an automated cell type annotation algorithm (SingleR), we showed that DS promoted endothelial-to-mesenchymal transition (EndMT) by inducing the transcriptional phenotypes for inflammation, hypoxia responses, transforming growth factor-beta (TGF-β) signaling, glycolysis, and fatty acid synthesis. Enolase 1 (ENO1), a key gene in glycolysis, was one of the top-ranked genes in the DS-induced EndMT cluster. Pseudotime trajectory analysis revealed that the kinetic expression of ENO1 was significantly associated with EndMT and that ENO1 silencing repressed the DS- and TGF-β-induced EC inflammation and EndMT. Consistent with these findings, ENO1 was highly expressed in ECs at the inner curvature of the mouse aortic arch (which is exposed to DS) and atherosclerotic lesions, suggesting its proatherogenic role in vivo. In summary, we present a comprehensive single-cell atlas of ECs in response to different flow patterns within the same flow channel. Among the DS-regulated genes, ENO1 plays an important role in DS-induced EC inflammation and EndMT. These results provide insights into how hemodynamic forces regulate vascular endothelium in health and disease.

EGFR activation differentially affects the inflammatory profiles of female human aortic and coronary artery endothelial cells

Endothelial cells (EC) are key players in vascular function, homeostasis and inflammation. EC show substantial heterogeneity due to inter-individual variability (e.g. sex-differences) and intra-individual differences as they originate from different organs or vessels. This variability may lead to different responsiveness to external stimuli. Here we compared the responsiveness of female human primary EC from the aorta (HAoEC) and coronary arteries (HCAEC) to Epidermal Growth Factor Receptor (EGFR) activation. EGFR is an important signal integration hub for vascular active substances with physiological and pathophysiological relevance. Our transcriptomic analysis suggested that EGFR activation differentially affects the inflammatory profiles of HAoEC and HCAEC, particularly by inducing a HCAEC-driven leukocyte attraction but a downregulation of adhesion molecule and chemoattractant expression in HAoEC. Experimental assessments of selected inflammation markers were performed to validate these predictions and the results confirmed a dual role of EGFR in these cells: its activation initiated an anti-inflammatory response in HAoEC but a pro-inflammatory one in HCAEC. Our study highlights that, although they are both arterial EC, female HAoEC and HCAEC are distinguishable with regard to the role of EGFR and its involvement in inflammation regulation, what may be relevant for vascular maintenance but also the pathogenesis of endothelial dysfunction.

K6-linked ubiquitylation marks formaldehyde-induced RNA-protein crosslinks for resolution

Reactive aldehydes are produced by normal cellular metabolism or after alcohol consumption, and they accumulate in human tissues if aldehyde clearance mechanisms are impaired. Their toxicity has been attributed to the damage they cause to genomic DNA and the subsequent inhibition of transcription and replication. However, whether interference with other cellular processes contributes to aldehyde toxicity has not been investigated. We demonstrate that formaldehyde induces RNA-protein crosslinks (RPCs) that stall the ribosome and inhibit translation in human cells. RPCs in the messenger RNA (mRNA) are recognized by the translating ribosomes, marked by atypical K6-linked ubiquitylation catalyzed by the RING-in-between-RING (RBR) E3 ligase RNF14, and subsequently resolved by the ubiquitin- and ATP-dependent unfoldase VCP. Our findings uncover an evolutionary conserved formaldehyde-induced stress response pathway that protects cells against RPC accumulation in the cytoplasm, and they suggest that RPCs contribute to the cellular and tissue toxicity of reactive aldehydes.

Differential Expression of RNAseq Imprinted Genes from Bovine Females Before and After Puberty

The productivity of beef cows depends on early reproduction traits such as puberty and has an economic impact on the efficiency of production system. Imprinted genes modulate many important endocrine processes such as growth, the onset of puberty and maternal reproductive and behavior. The role of imprinted genes in puberty is a challenging subject since they show the reciprocal role of maternal and paternal genomes in progeny. Although, there are evidences of the involvement of imprint genes in puberty in human, the role of this type of genes in the onset of puberty in cattle has not been studied yet. Here we examined the expression of 27 imprinted genes in pre and post puberty in a bovine model to find differentially expressed imprinted genes in maternal-paternal purebreds and reciprocal crosses across eight tissues and discussed the task of these genes in this crucial process of development and in onset of puberty. DLK1 and MKRN3 that previously described as cause of the central precocious puberty (CPP) in human were differentially expressed in this study. Functional annotation analysis of differentially imprinted genes in different tissues showed significant biological processes of cellular response to growth factor stimulus, response to growth factor, response to parathyroid hormone, developmental growth and the importance of alternative splicing. The results of this study have implications in understanding the role of imprinted genes in the onset of puberty in cattle.

Epithelial cell invasion by salmonella typhimurium induces modulation of genes controlled by aryl hydrocarbon receptor signaling and involved in extracellular matrix biogenesis

Salmonella is the only bacterium able to enter a host cell by the two known mechanisms: trigger and zipper. The trigger mechanism relies on the injection of bacterial effectors into the host cell through the Salmonella type III secretion system 1. In the zipper mechanism, mediated by the invasins Rck and PagN, the bacterium takes advantage of a cellular receptor for invasion. This study describes the transcriptomic reprogramming of the IEC-6 intestinal epithelial cell line to Salmonella Typhimurium strains that invaded cells by a trigger, a zipper, or both mechanisms. Using S. Typhimurium strains invalidated for one or other entry mechanism, we have shown that IEC-6 cells could support both entries. Comparison of the gene expression profiles of exposed cells showed that irrespective of the mechanism used for entry, the transcriptomic reprogramming of the cell was nearly the same. On the other hand, when gene expression was compared between cells unexposed or exposed to the bacterium, the transcriptomic reprogramming of exposed cells was significantly different. It is particularly interesting to note the modulation of expression of numerous target genes of the aryl hydrocarbon receptor showing that this transcription factor was activated by S. Typhimurium infection. Numerous genes associated with the extracellular matrix were also modified. This was confirmed at the protein level by western-blotting showing a dramatic modification in some extracellular matrix proteins. Analysis of a selected set of modulated genes showed that the expression of the majority of these genes was modulated during the intracellular life of S. Typhimurium.

Transcriptional impact of EGFR activation in human female vascular smooth muscle cells

Vascular smooth muscle cells (VSMC) are critical for the vascular tone, but they can also drive the development of vascular diseases when they lose their contractile phenotype and de-differentiate. Previous studies showed that the epidermal growth factor receptor (EGFR) of VSMC is critical for vascular health, but most of the underlying mechanisms by which VSMC-EGFR controls vascular fate have remained unknown. We combined RNA-sequencing and bioinformatics analysis to characterize the effect of EGFR-activation on the transcriptome of human primary VSMC (from different female donors) and to identify potentially affected cellular processes. Our results indicate that the activation of human VSMC-EGFR is sufficient to trigger a phenotypical switch toward a proliferative and inflammatory phenotype. The extent of this effect is nonetheless partly donor-dependent. Our hypothesis-generating study thus provides a first insight into mechanisms that could partly explain variable susceptibilities to vascular diseases in between individuals.

TLR7 promotes smoke-induced experimental lung damage through the activity of mast cell tryptase

Toll-like receptor 7 (TLR7) is known for eliciting immunity against single-stranded RNA viruses, and is increased in both human and cigarette smoke (CS)-induced, experimental chronic obstructive pulmonary disease (COPD). Here we show that the severity of CS-induced emphysema and COPD is reduced in TLR7-deficient mice, while inhalation of imiquimod, a TLR7-agonist, induces emphysema without CS exposure. This imiquimod-induced emphysema is reduced in mice deficient in mast cell protease-6, or when wild-type mice are treated with the mast cell stabilizer, cromolyn. Furthermore, therapeutic treatment with anti-TLR7 monoclonal antibody suppresses CS-induced emphysema, experimental COPD and accumulation of pulmonary mast cells in mice. Lastly, TLR7 mRNA is increased in pre-existing datasets from patients with COPD, while TLR7+ mast cells are increased in COPD lungs and associated with severity of COPD. Our results thus support roles for TLR7 in mediating emphysema and COPD through mast cell activity, and may implicate TLR7 as a potential therapeutic target.

YWHAG Deficiency Disrupts the EMT-Associated Network to Induce Oxidative Cell Death and Prevent Metastasis

Metastasis involves epithelial-to-mesenchymal transition (EMT), a process that is regulated by complex gene networks, where their deliberate disruption may yield a promising outcome. However, little is known about mechanisms that coordinate these metastasis-associated networks. To address this gap, hub genes with broad engagement across various human cancers by analyzing the transcriptomes of different cancer cell types undergoing EMT are identified. The oncogenic signaling adaptor protein tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein gamma (YWHAG) is ranked top for its clinical relevance and impact. The cellular kinome and transcriptome data are surveyed to construct the regulome of YWHAG, revealing stress responses and metabolic processes during cancer EMT. It is demonstrated that a YWHAG-dependent cytoprotective mechanism in the regulome is embedded in EMT-associated networks to protect cancer cells from oxidative catastrophe through enhanced autophagy during EMT. YWHAG deficiency results in a rapid accumulation of reactive oxygen species (ROS), delayed EMT, and cell death. Tumor allografts show that metastasis potential and overall survival time are correlated with the YWHAG expression level of cancer cell lines. Metastasized tumors have higher expression of YWHAG and autophagy-related genes than primary tumors. Silencing YWHAG diminishes primary tumor volumes, prevents metastasis, and prolongs the median survival period of the mice.

Candida albicans exhibits heterogeneous and adaptive cytoprotective responses to antifungal compounds

Candida albicans, an opportunistic human pathogen, poses a significant threat to human health and is associated with significant socio-economic burden. Current antifungal treatments fail, at least in part, because C. albicans can initiate a strong drug tolerance response that allows some cells to grow at drug concentrations above their minimal inhibitory concentration. To better characterize this cytoprotective tolerance program at the molecular single-cell level, we used a nanoliter droplet-based transcriptomics platform to profile thousands of individual fungal cells and establish their subpopulation characteristics in the absence and presence of antifungal drugs. Profiles of untreated cells exhibit heterogeneous expression that correlates with cell cycle stage with distinct metabolic and stress responses. At 2 days post-fluconazole exposure (a time when tolerance is measurable), surviving cells bifurcate into two major subpopulations: one characterized by the upregulation of genes encoding ribosomal proteins, rRNA processing machinery, and mitochondrial cellular respiration capacity, termed the Ribo-dominant (Rd) state; and the other enriched for genes encoding stress responses and related processes, termed the Stress-dominant (Sd) state. This bifurcation persists at 3 and 6 days post-treatment. We provide evidence that the ribosome assembly stress response (RASTR) is activated in these subpopulations and may facilitate cell survival.

Single-cell census of human tooth development enables generation of human enamel

Tooth enamel secreted by ameloblasts (AMs) is the hardest material in the human body, acting as a shield to protect the teeth. However, the enamel is gradually damaged or partially lost in over 90% of adults and cannot be regenerated due to a lack of ameloblasts in erupted teeth. Here, we use single-cell combinatorial indexing RNA sequencing (sci-RNA-seq) to establish a spatiotemporal single-cell census for the developing human tooth and identify regulatory mechanisms controlling the differentiation process of human ameloblasts. We identify key signaling pathways involved between the support cells and ameloblasts during fetal development and recapitulate those findings in human ameloblast in vitro differentiation from induced pluripotent stem cells (iPSCs). We furthermore develop a disease model of amelogenesis imperfecta in a three-dimensional (3D) organoid system and show AM maturation to mineralized structure in vivo. These studies pave the way for future regenerative dentistry.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Transcriptomic Analysis of Light-Induced Genes in Nasonia vitripennis: Possible Implications for Circadian Light Entrainment Pathways

Circadian entrainment to the environmental day-night cycle is essential for the optimal use of environmental resources. In insects, opsin-based photoreception in the compound eye and ocelli and CRYPTOCHROME1 (CRY1) in circadian clock neurons are thought to be involved in sensing photic information, but the genetic regulation of circadian light entrainment in species without light-sensitive CRY1 remains unclear. To elucidate a possible CRY1-independent light transduction cascade, we analyzed light-induced gene expression through RNA-sequencing in Nasonia vitripennis. Entrained wasps were subjected to a light pulse in the subjective night to reset the circadian clock, and light-induced changes in gene expression were characterized at four different time points in wasp heads. We used co-expression, functional annotation, and transcription factor binding motif analyses to gain insight into the molecular pathways in response to acute light stimulus and to form hypotheses about the circadian light-resetting pathway. Maximal gene induction was found after 2 h of light stimulation (1432 genes), and this included the opsin gene opblue and the core clock genes cry2 and npas2. Pathway and cluster analyses revealed light activation of glutamatergic and GABA-ergic neurotransmission, including CREB and AP-1 transcription pathway signaling. This suggests that circadian photic entrainment in Nasonia may require pathways that are similar to those in mammals. We propose a model for hymenopteran circadian light-resetting that involves opsin-based photoreception, glutamatergic neurotransmission, and gene induction of cry2 and npas2 to reset the circadian clock.

The CUL4B-based E3 ubiquitin ligase regulates mitosis and brain development by recruiting phospho-specific DCAFs

The paralogs CUL4A and CUL4B assemble cullin-RING E3 ubiquitin ligase (CRL) complexes regulating multiple chromatin-associated cellular functions. Although they are structurally similar, we found that the unique N-terminal extension of CUL4B is heavily phosphorylated during mitosis, and the phosphorylation pattern is perturbed in the CUL4B-P50L mutation causing X-linked intellectual disability (XLID). Phenotypic characterization and mutational analysis revealed that CUL4B phosphorylation is required for efficient progression through mitosis, controlling spindle positioning and cortical tension. While CUL4B phosphorylation triggers chromatin exclusion, it promotes binding to actin regulators and to two previously unrecognized CUL4B-specific substrate receptors (DCAFs), LIS1 and WDR1. Indeed, co-immunoprecipitation experiments and biochemical analysis revealed that LIS1 and WDR1 interact with DDB1, and their binding is enhanced by the phosphorylated N-terminal domain of CUL4B. Finally, a human forebrain organoid model demonstrated that CUL4B is required to develop stable ventricular structures that correlate with onset of forebrain differentiation. Together, our study uncovers previously unrecognized DCAFs relevant for mitosis and brain development that specifically bind CUL4B, but not the CUL4B-P50L patient mutant, by a phosphorylation-dependent mechanism.

Viruses in Laboratory Drosophila and Their Impact on Host Gene Expression

Drosophila melanogaster has one of the best characterized antiviral immune responses among invertebrates. However, relatively few easily transmitted natural virus isolates are available, and so many Drosophila experiments have been performed using artificial infection routes and artificial host-virus combinations. These may not reflect natural infections, especially for subtle phenotypes such as gene expression. Here, to explore the laboratory virus community and to better understand how natural virus infections induce changes in gene expression, we have analysed seven publicly available D. melanogaster transcriptomic sequencing datasets that were originally sequenced for projects unrelated to virus infection. We have found ten known viruses-including five that have not been experimentally isolated-but no previously unknown viruses. Our analysis of host gene expression revealed that numerous genes were differentially expressed in flies that were naturally infected with a virus. For example, flies infected with nora virus showed patterns of gene expression consistent with intestinal vacuolization and possible host repair via the upd3 JAK/STAT pathway. We also found marked sex differences in virus-induced differential gene expression. Our results show that natural virus infection in laboratory Drosophila does indeed induce detectable changes in gene expression, suggesting that this may form an important background condition for experimental studies in the laboratory.

Large scale microfluidic CRISPR screening for increased amylase secretion in yeast

Key to our ability to increase recombinant protein production through secretion is a better understanding of the pathways that interact to translate, process and export mature proteins to the surrounding environment, including the supporting cellular machinery that supplies necessary energy and building blocks. By combining droplet microfluidic screening with large-scale CRISPR libraries that perturb the expression of the majority of coding and non-coding genes in S. cerevisiae, we identified 345 genes for which an increase or decrease in gene expression resulted in increased secretion of α-amylase. Our results show that modulating the expression of genes involved in the trafficking of vesicles, endosome to Golgi transport, the phagophore assembly site, the cell cycle and energy supply improve α-amylase secretion. Besides protein-coding genes, we also find multiple long non-coding RNAs enriched in the vicinity of genes associated with endosomal, Golgi and vacuolar processes. We validated our results by overexpressing or deleting selected genes, which resulted in significant improvements in α-amylase secretion. The advantages, in terms of precision and speed, inherent to CRISPR based perturbations, enables iterative testing of new strains for increased protein secretion.

Y-Linked Copy Number Polymorphism of Target of Rapamycin Is Associated with Sexual Size Dimorphism in Seed Beetles

The Y chromosome is theorized to facilitate evolution of sexual dimorphism by accumulating sexually antagonistic loci, but empirical support is scarce. Due to the lack of recombination, Y chromosomes are prone to degenerative processes, which poses a constraint on their adaptive potential. Yet, in the seed beetle, Callosobruchus maculatus segregating Y linked variation affects male body size and thereby sexual size dimorphism (SSD). Here, we assemble C. maculatus sex chromosome sequences and identify molecular differences associated with Y-linked SSD variation. The assembled Y chromosome is largely euchromatic and contains over 400 genes, many of which are ampliconic with a mixed autosomal and X chromosome ancestry. Functional annotation suggests that the Y chromosome plays important roles in males beyond primary reproductive functions. Crucially, we find that, besides an autosomal copy of the gene target of rapamycin (TOR), males carry an additional TOR copy on the Y chromosome. TOR is a conserved regulator of growth across taxa, and our results suggest that a Y-linked TOR provides a male specific opportunity to alter body size. A comparison of Y haplotypes associated with male size difference uncovers a copy number variation for TOR, where the haplotype associated with decreased male size, and thereby increased sexual dimorphism, has two additional TOR copies. This suggests that sexual conflict over growth has been mitigated by autosome to Y translocation of TOR followed by gene duplications. Our results reveal that despite of suppressed recombination, the Y chromosome can harbor adaptive potential as a male-limited supergene.

In vitro analysis of PI3K pathway activation genes for exploring novel biomarkers and therapeutic targets in clear cell renal carcinoma

OBJECTIVES

The regulation of various cellular functions such as growth, proliferation, metabolism, and angiogenesis, is dependent on the PI3K pathway. Recent evidence has indicated that kidney renal clear cell carcinoma (KIRC) can be triggered by the deregulation of this pathway. The objective of this research was to investigate 25 genes associated with activation of the PI3K pathway in KIRC and control samples to identify four hub genes that might serve as novel molecular biomarkers and therapeutic targets for treating KIRC.

METHODS

Multi-omics in silico and in vitro analysis was employed to find hub genes related to the PI3K pathway that may be biomarkers and therapeutic targets for KIRC.

RESULTS

Using STRING software, a protein-protein interaction (PPI) network of 25 PI3K pathway-related genes was developed. Based on the degree scoring method, the top four hub genes were identified using Cytoscape's Cytohubba plug-in. TCGA datasets, KIRC (786-O and A-498), and normal (HK2) cells were used to validate the expression of hub genes. Additionally, further bioinformatic analyses were performed to investigate the mechanisms by which hub genes are involved in the development of KIRC. Out of a total of 25 PI3K pathway-related genes, we developed and validated a diagnostic and prognostic model based on the up-regulation of TP53 (tumor protein 53) and CCND1 (Cyclin D1) and the down-regulation of PTEN (Phosphatase and TENsin homolog deleted on chromosome 10), and GSK3B (Glycogen synthase kinase-3 beta) hub genes. The hub genes included in our model may be a novel therapeutic target for KIRC treatment. Additionally, associations between hub genes and infiltration of immune cells can enhance comprehension of immunotherapy for KIRC.

CONCLUSION

We have created a new diagnostic and prognostic model for KIRC patients that uses PI3K pathway-related hub genes (TP53, PTEN, CCND1, and GSK3B). Nevertheless, further experimental studies are required to ascertain the efficacy of our model.

The making of the oral microbiome in Agta hunter-gatherers

Ecological and genetic factors have influenced the composition of the human microbiome during our evolutionary history. We analysed the oral microbiota of the Agta, a hunter-gatherer population where some members have adopted an agricultural diet. We show that age is the strongest factor modulating the microbiome, probably through immunosenescence since we identified an increase in the number of species classified as pathogens with age. We also characterised biological and cultural processes generating sexual dimorphism in the oral microbiome. A small subset of oral bacteria is influenced by the host genome, linking host collagen genes to bacterial biofilm formation. Our data also suggest that shifting from a fish/meat diet to a rice-rich diet transforms their microbiome, mirroring the Neolithic transition. All of these factors have implications in the epidemiology of oral diseases. Thus, the human oral microbiome is multifactorial and shaped by various ecological and social factors that modify the oral environment.

Differential Methylation of Telomere-Related Genes Is Associated with Kidney Disease in Individuals with Type 1 Diabetes

Diabetic kidney disease (DKD) represents a major global health problem. Accelerated ageing is a key feature of DKD and, therefore, characteristics of accelerated ageing may provide useful biomarkers or therapeutic targets. Harnessing multi-omics, features affecting telomere biology and any associated methylome dysregulation in DKD were explored. Genotype data for nuclear genome polymorphisms in telomere-related genes were extracted from genome-wide case-control association data (n = 823 DKD/903 controls; n = 247 end-stage kidney disease (ESKD)/1479 controls). Telomere length was established using quantitative polymerase chain reaction. Quantitative methylation values for 1091 CpG sites in telomere-related genes were extracted from epigenome-wide case-control association data (n = 150 DKD/100 controls). Telomere length was significantly shorter in older age groups (p = 7.6 × 10^-6). Telomere length was also significantly reduced (p = 6.6 × 10^-5) in DKD versus control individuals, with significance remaining after covariate adjustment (p = 0.028). DKD and ESKD were nominally associated with telomere-related genetic variation, with Mendelian randomisation highlighting no significant association between genetically predicted telomere length and kidney disease. A total of 496 CpG sites in 212 genes reached epigenome-wide significance (p ≤ 10^-8) for DKD association, and 412 CpG sites in 193 genes for ESKD. Functional prediction revealed differentially methylated genes were enriched for Wnt signalling involvement. Harnessing previously published RNA-sequencing datasets, potential targets where epigenetic dysregulation may result in altered gene expression were revealed, useful as potential diagnostic and therapeutic targets for intervention.

Temporal tracking of microglial and monocyte single-cell transcriptomics in lethal flavivirus infection

As the resident parenchymal myeloid population in the central nervous system (CNS), microglia are strategically positioned to respond to neurotropic virus invasion and have been implicated in promoting both disease resolution and progression in the acute and post-infectious phase of virus encephalitis. In a mouse model of West Nile virus encephalitis (WNE), infection of the CNS results in recruitment of large numbers of peripheral immune cells into the brain, the majority being nitric oxide (NO)-producing Ly6C^hi inflammatory monocyte-derived cells (MCs). In this model, these cells enhance immunopathology and mortality. However, the contribution of microglia to this response is currently undefined. Here we used a combination of experimental tools, including single-cell RNA sequencing (scRNA-seq), microglia and MC depletion reagents, high-dimensional spectral cytometry and computational algorithms to dissect the differential contribution of microglia and MCs to the anti-viral immune response in severe neuroinflammation seen in WNE. Intriguingly, analysis of scRNA-seq data revealed 6 unique microglia and 3 unique MC clusters that were predominantly timepoint-specific, demonstrating substantial transcriptional adaptation with disease progression over the course of WNE. While microglia and MC adopted unique gene expression profiles, gene ontology enrichment analysis, coupled with microglia and MC depletion studies, demonstrated a role for both of these cells in the trafficking of peripheral immune cells into the CNS, T cell responses and viral clearance. Over the course of infection, microglia transitioned from a homeostatic to an anti-viral and then into an immune cell-recruiting phenotype. Conversely, MC adopted antigen-presenting, immune cell-recruiting and NO-producing phenotypes, which all had anti-viral function. Overall, this study defines for the first time the single-cell transcriptomic responses of microglia and MCs over the course of WNE, demonstrating both protective and pathological roles of these cells that could potentially be targeted for differential therapeutic intervention to dampen immune-mediated pathology, while maintaining viral clearance functions.

Integrating transcriptomic datasets across neurological disease identifies unique myeloid subpopulations driving disease-specific signatures

Microglia and bone marrow-derived monocytes are key elements of central nervous system (CNS) inflammation, both capable of enhancing and dampening immune-mediated pathology. However, the study-specific focus on individual cell types, disease models or experimental approaches has limited our ability to infer common and disease-specific responses. This meta-analysis integrates bulk and single-cell transcriptomic datasets of microglia and monocytes from disease models of autoimmunity, neurodegeneration, sterile injury, and infection to build a comprehensive resource connecting myeloid responses across CNS disease. We demonstrate that the bulk microglial and monocyte program is highly contingent on the disease environment, challenging the notion of a universal microglial disease signature. Integration of six single-cell RNA-sequencing datasets revealed that these disease-specific signatures are likely driven by differing proportions of unique myeloid subpopulations that were individually expanded in different disease settings. These subsets were functionally-defined as neurodegeneration-associated, inflammatory, interferon-responsive, phagocytic, antigen-presenting, and lipopolysaccharide-responsive cellular states, revealing a core set of myeloid responses at the single-cell level that are conserved across CNS pathology. Showcasing the predictive and practical value of this resource, we performed differential expression analysis on microglia and monocytes across disease and identified Cd81 as a new neuroinflammatory-stable gene that accurately identified microglia and distinguished them from monocyte-derived cells across all experimental models at both the bulk and single-cell level. Together, this resource dissects the influence of disease environment on shared immune response programmes to build a unified perspective of myeloid behavior across CNS pathology.

Transcriptome changes in chlorsulfuron-treated plants are caused by acetolactate synthase inhibition and not induction of a herbicide detoxification system in Marchantia polymorpha

A sensing mechanism in mammals perceives xenobiotics and induces the transcription of genes encoding proteins that detoxify these molecules. However, it is unclear if plants sense xenobiotics, and activate an analogous signalling system leading to their detoxification. Using the liverwort Marchantia polymorpha, we tested the hypothesis that there is a sensing system in plants that perceives herbicides resulting in the increased transcription of genes encoding proteins that detoxify these herbicides. Consistent with the hypothesis, we show that chlorsulfuron-treatment induces changes in the M. polymorpha transcriptome. However, these transcriptome changes do not occur in chlorsulfuron (CS)-treated target site resistant mutants, where the gene encoding the target carries a mutation that confers resistance to chlorsulfuron. Instead, we show that inactivation of the chlorsulfuron target, acetolactate synthase (ALS) (also known as acetohydroxyacid synthase (AHAS)), is required for the transcriptome response. These data demonstrate that the transcriptome changes in chlorsulfuron-treated plants are caused by disrupted amino acid synthesis and metabolism resulting from acetolactate synthase inhibition, and indicate that the transcriptome changes are not caused by a herbicide sensing mechanism.

Combination of bacterial N-acyl homoserine lactones primes Arabidopsis defenses via jasmonate metabolism

N-acyl homoserine lactones (AHLs) are important players in plant-bacteria interactions. Different AHL-producing bacteria can improve plant growth and resistance against plant pathogens. In nature, plants may host a variety of AHL-producing bacteria and frequently experience numerous AHLs at the same time. Therefore, a coordinated response to combined AHL molecules is necessary. The purpose of this study was to explore the mechanism of AHL-priming using combined AHL molecules including N-(3-oxo-hexanoyl)-L-homoserine lactone, N-3-oxo-octanoyl-L-homoserine lactone, N-3-oxo-dodecanoyl-L-homoserine lactone, and N-3-oxo-tetradecanoyl-L-homoserine lactone and AHL-producing bacteria including Serratia plymuthica HRO-C48, Rhizobium etli CFN42, Burkholderia graminis DSM17151, and Ensifer meliloti (Sinorhizobium meliloti) Rm2011. We used transcriptome analysis, phytohormone measurements, as well as genetic and microbiological approaches to assess how the combination of structurally diverse AHL molecules influence Arabidopsis (Arabidopsis thaliana). Our findings revealed a particular response to a mixture of AHL molecules (AHL mix). Different expression patterns indicated that the reaction of plants exposed to AHL mix differs from that of plants exposed to single AHL molecules. In addition, different content of jasmonic acid (JA) and derivatives revealed that jasmonates play an important role in AHL mix-induced priming. The fast and stable decreased concentration of COOH-JA-Ile after challenge with the flagellin-derived peptide flg22 indicated that AHL mix modifies the metabolism of jasmonates. Study of various JA- and salicylic acid-related Arabidopsis mutants strengthened the notion that JA homeostasis is involved in AHL-priming. Understanding how the combination of AHLs primes plants for enhanced resistance has the potential to broaden our approaches in sustainable agriculture and will help to effectively protect plants against pathogens.

TiSA: TimeSeriesAnalysis-a pipeline for the analysis of longitudinal transcriptomics data

Improved transcriptomic sequencing technologies now make it possible to perform longitudinal experiments, thus generating a large amount of data. Currently, there are no dedicated or comprehensive methods for the analysis of these experiments. In this article, we describe our TimeSeries Analysis pipeline (TiSA) which combines differential gene expression, clustering based on recursive thresholding, and a functional enrichment analysis. Differential gene expression is performed for both the temporal and conditional axes. Clustering is performed on the identified differentially expressed genes, with each cluster being evaluated using a functional enrichment analysis. We show that TiSA can be used to analyse longitudinal transcriptomic data from both microarrays and RNA-seq, as well as small, large, and/or datasets with missing data points. The tested datasets ranged in complexity, some originating from cell lines while another was from a longitudinal experiment of severity in COVID-19 patients. We have also included custom figures to aid with the biological interpretation of the data, these plots include Principal Component Analyses, Multi Dimensional Scaling plots, functional enrichment dotplots, trajectory plots, and complex heatmaps showing the broad overview of results. To date, TiSA is the first pipeline to provide an easy solution to the analysis of longitudinal transcriptomics experiments.

Bisphenol A and bisphenol S both disrupt ovine granulosa cell steroidogenesis but through different molecular pathways

BACKGROUND

Ovarian granulosa cells (GC) are essential for the development and maturation of a proper oocyte. GC are sensitive to endocrine disruptors, including bisphenol A (BPA) and its analogue bisphenol S (BPS), plasticisers present in everyday consumer products. BPA exhibits greater binding affinity for the membrane oestrogen receptor (GPER) than for the nuclear oestrogen receptors (ERα and ERβ). Here, we analysed the effects of BPA and BPS on the steroidogenesis of ovine GC in vitro, as well as their early mechanisms of action, the ovine being a relevant model to study human reproductive impairment. Disruption of GC steroidogenesis might alter oocyte quality and consequently fertility rate. In addition, we compared the effects of a specific GPER agonist (G-1) and antagonist (G-15) to those of BPA and BPS. Ewe GC were cultured with BPA or BPS (10 or 50 µM) or G-1 (1 µM) and/or G-15 (10 µM) for 48 h to study steroidogenesis.

RESULTS

Both BPA and BPS (10 µM) altered the secretion of progesterone, however, only BPS (10 µM) affected oestradiol secretion. RNA-seq was performed on GC after 1 h of culture with BPA or BPS (50 µM) or G-1 (10 µM), followed by real-time PCR analyses of differentially expressed genes after 12, 24 and 48 h of culture. The absence of induced GPER target genes showed that BPA and BPS did not activate GPER in GC after 1 h of treatment. These molecules exhibited mainly independent early mechanisms of action. Gene ontology analysis showed that after 1 h of treatment, BPA mainly disrupted the expression of the genes involved in metabolism and transcription, while BPS had a smaller effect and impaired cellular communications. BPA had a transient effect on the expression of CHAC1 (NOTCH signalling and oxidative balance), JUN (linked to MAPK pathway), NR4A1 (oestradiol secretion inhibition), ARRDC4 (endocytose of GPCR) and KLF10 (cell growth, differentiation and apoptosis), while expression changes were maintained over time for the genes LSMEM1 (linked to MAPK pathway), TXNIP (oxidative stress) and LIF (cell cycle regulation) after 12 and 48 h, respectively.

CONCLUSION

In conclusion, although they exhibited similar effects, BPA and BPS impaired different molecular pathways in GC in vitro. New investigations will be necessary to follow the temporal changes of these genes over time, as well as the biological processes involved.

Protein network analysis and functional enrichment via computational biotechnology unravel molecular and pathogenic mechanisms of kidney stone disease

Mass spectrometry-based proteomics has been extensively applied to current biomedical research. From such large-scale identification of proteins, several computational tools have been developed for determining protein-protein interactions (PPI) network and functional significance of the identified proteins and their complex. Analyses of PPI network and functional enrichment have been widely applied to various fields of biomedical research. Herein, we summarize commonly used tools for PPI network analysis and functional enrichment in kidney stone research and discuss their applications to kidney stone disease (KSD). Such computational approach has been used mainly to investigate PPI networks and functional significance of the proteins derived from urine of patients with kidney stone (stone formers), stone matrix, Randall's plaque, renal papilla, renal tubular cells, mitochondria and immune cells. The data obtained from computational biotechnology leads to experimental validation and investigations that offer new knowledge on kidney stone formation processes. Moreover, the computational approach may also lead to defining new therapeutic targets and preventive strategies for better outcome in KSD management.

rrvgo: a Bioconductor package for interpreting lists of Gene Ontology terms

Gene Ontology (GO) annotation is often used to guide the biological interpretation of high-throughput omics experiments, e.g. by analysing lists of differentially regulated genes for enriched GO terms. Due to the hierarchical nature of GOs, the resulting lists of enriched terms are usually redundant and difficult to summarise and interpret. To facilitate the interpretation of large lists of GO terms, I developed rrvgo, a Bioconductor package that aims at simplifying the redundancy of GO lists by grouping similar terms based on their semantic similarity. rrvgo also provides different visualization options to guide the interpretation of the summarized GO terms. Considering that several software tools have been developed for this purpose, rrvgo is unique at combining powerful visualizations in a programmatic interface coupled with up-to-date GO gene annotation provided by the Bioconductor project.

Asymmetric gene expression in grain development of reciprocal crosses between tetraploid and hexaploid wheats

Production of viable progeny from interploid crosses requires precise regulation of gene expression from maternal and paternal chromosomes, yet the transcripts contributed to hybrid seeds from polyploid parent species have rarely been explored. To investigate the genome-wide maternal and paternal contributions to polyploid grain development, we analyzed the transcriptomes of developing embryos, from zygote to maturity, alongside endosperm in two stages of development, using reciprocal crosses between tetraploid and hexaploid wheats. Reciprocal crosses between species with varied levels of ploidy displayed broad impacts on gene expression, including shifts in alternative splicing events in select crosses, as illustrated by active splicing events, enhanced protein synthesis and chromatin remodeling. Homoeologous gene expression was repressed on the univalent D genome in pentaploids, but this suppression was attenuated in crosses with a higher ploidy maternal parent. Imprinted genes were identified in endosperm and early embryo tissues, supporting predominant maternal effects on early embryogenesis. By systematically investigating the complex transcriptional networks in reciprocal-cross hybrids, this study presents a framework for understanding the genomic incompatibility and transcriptome shock that results from interspecific hybridization and uncovers the transcriptional impacts on hybrid seeds created from agriculturally-relevant polyploid species.

The establishment of COPD organoids to study host-pathogen interaction reveals enhanced viral fitness of SARS-CoV-2 in bronchi

Chronic obstructive pulmonary disease (COPD) is characterised by airflow limitation and infective exacerbations, however, in-vitro model systems for the study of host-pathogen interaction at the individual level are lacking. Here, we describe the establishment of nasopharyngeal and bronchial organoids from healthy individuals and COPD that recapitulate disease at the individual level. In contrast to healthy organoids, goblet cell hyperplasia and reduced ciliary beat frequency were observed in COPD organoids, hallmark features of the disease. Single-cell transcriptomics uncovered evidence for altered cellular differentiation trajectories in COPD organoids. SARS-CoV-2 infection of COPD organoids revealed more productive replication in bronchi, the key site of infection in severe COVID-19. Viral and bacterial exposure of organoids induced greater pro-inflammatory responses in COPD organoids. In summary, we present an organoid model that recapitulates the in vivo physiological lung microenvironment at the individual level and is amenable to the study of host-pathogen interaction and emerging infectious disease.

Transcriptomic FHITlow/pHER2high signature as a predictive factor of outcome and immunotherapy response in non-small cell lung cancer

Introduction

In recent decades, the development of immunotherapy and targeted therapies has considerably improved the outcome of non-small cell lung cancer (NSCLC) patients. Despite these impressive clinical benefits, new biomarkers are needed for an accurate stratification of NSCLC patients and a more personalized management. We recently showed that the tumor suppressor fragile histidine triad (FHIT), frequently lost in NSCLC, controls HER2 receptor activity in lung tumor cells and that tumor cells from NSCLC patients harboring a FHIT^low/pHER2^high phenotype are sensitive to anti-HER2 drugs. Here, we sought to identify the transcriptomic signature of this phenotype and evaluate its clinical significance.

Materials and methods

We performed RNA sequencing analysis on tumor cells isolated from NSCLC (n=12) according to FHIT/pHER2 status and a functional analysis of differentially regulated genes. We also investigated the FHIT^low/pHER2^high signature in The Cancer Genome Atlas (TCGA) lung adenocarcinoma (LUAD) (n=489) and lung squamous cell carcinoma (LUSC) (n=493) cohorts and used the tumor immune dysfunction and exclusion (TIDE) model to test the ability of this signature to predict response to immune checkpoint inhibitors (ICI).

Results

We showed that up-regulated genes in FHIT^low/pHER2^high tumors were associated with cell proliferation, metabolism and metastasis, whereas down-regulated genes were related to immune response. The FHIT^low/pHER2^high signature was associated with the higher size of tumors, lymph node involvement, and late TNM stages in LUAD and LUSC cohorts. It was identified as an independent predictor of overall survival (OS) in LUAD cohort. FHIT^low/pHER2^high tumors were also predictive of poor response to ICI in both LUAD and LUSC cohorts.

Conclusion

These data suggest that ICI might not be a relevant option for NSCLC patients with FHIT^low/pHER2^high tumors and that anti-HER2 targeted therapy could be a good therapeutic alternative for this molecular subclass with poorer prognosis.

orsum: a Python package for filtering and comparing enrichment analyses using a simple principle

Background

Enrichment analyses are widely applied to investigate lists of genes of interest. However, such analyses often result in long lists of annotation terms with high redundancy, making the interpretation and reporting difficult. Long annotation lists and redundancy also complicate the comparison of results obtained from different enrichment analyses. An approach to overcome these issues is using down-sized annotation collections composed of non-redundant terms. However, down-sized collections are generic and the level of detail may not fit the user’s study. Other available approaches include clustering and filtering tools, which are based on similarity measures and thresholds that can be complicated to comprehend and set.

Result

We propose orsum, a Python package to filter enrichment results. orsum can filter multiple enrichment results collectively and highlight common and specific annotation terms. Filtering in orsum is based on a simple principle: a term is discarded if there is a more significant term that annotates at least the same genes; the remaining more significant term becomes the representative term for the discarded term. This principle ensures that the main biological information is preserved in the filtered results while reducing redundancy. In addition, as the representative terms are selected from the original enrichment results, orsum outputs filtered terms tailored to the study. As a use case, we applied orsum to the enrichment analyses of four lists of genes, each associated with a neurodegenerative disease.

Conclusion

orsum provides a comprehensible and effective way of filtering and comparing enrichment results. It is available at https://anaconda.org/bioconda/orsum.

PARP1 proximity proteomics reveals interaction partners at stressed replication forks

PARP1 mediates poly-ADP-ribosylation of proteins on chromatin in response to different types of DNA lesions. PARP inhibitors are used for the treatment of BRCA1/2-deficient breast, ovarian, and prostate cancer. Loss of DNA replication fork protection is proposed as one mechanism that contributes to the vulnerability of BRCA1/2-deficient cells to PARP inhibitors. However, the mechanisms that regulate PARP1 activity at stressed replication forks remain poorly understood. Here, we performed proximity proteomics of PARP1 and isolation of proteins on stressed replication forks to map putative PARP1 regulators. We identified TPX2 as a direct PARP1-binding protein that regulates the auto-ADP-ribosylation activity of PARP1. TPX2 interacts with DNA damage response proteins and promotes homology-directed repair of DNA double-strand breaks. Moreover, TPX2 mRNA levels are increased in BRCA1/2-mutated breast and prostate cancers, and high TPX2 expression levels correlate with the sensitivity of cancer cells to PARP-trapping inhibitors. We propose that TPX2 confers a mitosis-independent function in the cellular response to replication stress by interacting with PARP1.

Ecological Speciation Promoted by Divergent Regulation of Functional Genes Within African Cichlid Fishes

Rapid ecological speciation along depth gradients has taken place repeatedly in freshwater fishes, yet molecular mechanisms facilitating such diversification are typically unclear. In Lake Masoko, an African crater lake, the cichlid Astatotilapia calliptera has diverged into shallow-littoral and deep-benthic ecomorphs with strikingly different jaw structures within the last 1,000 years. Using genome-wide transcriptome data, we explore two major regulatory transcriptional mechanisms, expression and splicing-QTL variants, and examine their contributions to differential gene expression underpinning functional phenotypes. We identified 7,550 genes with significant differential expression between ecomorphs, of which 5.4% were regulated by cis-regulatory expression QTLs, and 9.2% were regulated by cis-regulatory splicing QTLs. We also found strong signals of divergent selection on differentially expressed genes associated with craniofacial development. These results suggest that large-scale transcriptome modification plays an important role during early-stage speciation. We conclude that regulatory variants are important targets of selection driving ecologically relevant divergence in gene expression during adaptive diversification.

Transcriptional regulation underlying the temperature response of embryonic development rate in the winter moth

Climate change will strongly affect the developmental timing of insects, as their development rate depends largely on ambient temperature. However, we know little about the genetic mechanisms underlying the temperature sensitivity of embryonic development in insects. We investigated embryonic development rate in the winter moth (Operophtera brumata), a species with egg dormancy which has been under selection due to climate change. We used RNA sequencing to investigate which genes are involved in the regulation of winter moth embryonic development rate in response to temperature. Over the course of development, we sampled eggs before and after an experimental change in ambient temperature, including two early development weeks when the temperature sensitivity of eggs is low and two late development weeks when temperature sensitivity is high. We found temperature-responsive genes that responded in a similar way across development, as well as genes with a temperature response specific to a particular development week. Moreover, we identified genes whose temperature effect size changed around the switch in temperature sensitivity of development rate. Interesting candidate genes for regulating the temperature sensitivity of egg development rate included genes involved in histone modification, hormonal signalling, nervous system development and circadian clock genes. The diverse sets of temperature-responsive genes we found here indicate that there are many potential targets of selection to change the temperature sensitivity of embryonic development rate. Identifying for which of these genes there is genetic variation in wild insect populations will give insight into their adaptive potential in the face of climate change.

Alternative developmental and transcriptomic responses to host plant water limitation in a butterfly metapopulation

Predicting how climate change affects biotic interactions poses a challenge. Plant-insect herbivore interactions are particularly sensitive to climate change, as climate-induced changes in plant quality cascade into the performance of insect herbivores. Whereas the immediate survival of herbivore individuals depends on plastic responses to climate change-induced nutritional stress, long-term population persistence via evolutionary adaptation requires genetic variation for these responses. To assess the prospects for population persistence under climate change, it is therefore crucial to characterize response mechanisms to climate change-induced stressors, and quantify their variability in natural populations. Here, we test developmental and transcriptomic responses to water limitation-induced host plant quality change in a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We combine nuclear magnetic resonance spectroscopy on the plant metabolome, larval developmental assays and an RNA sequencing analysis of the larval transcriptome. We observed that responses to feeding on water-limited plants, in which amino acids and aromatic compounds are enriched, showed marked variation within the metapopulation, with individuals of some families performing better on control and others on water-limited plants. The transcriptomic responses were concordant with the developmental responses: families exhibiting opposite developmental responses also produced opposite transcriptomic responses (e.g. in growth-associated transcripts). The divergent responses in both larval development and transcriptome are associated with differences between families in amino acid catabolism and storage protein production. The results reveal intrapopulation variability in plasticity, suggesting that the Finnish M. cinxia metapopulation harbours potential for buffering against drought-induced changes in host plant quality.