Temporal Changes In Gene Expression Research Articles

Story of an infection: Viral dynamics and host responses in the Caenorhabditis elegans-Orsay virus pathosystem.

Orsay virus (OrV) is the only known natural virus affecting Caenorhabditis elegans, with minimal impact on the animal's fitness due to its robust innate immune response. This study aimed to understand the interactions between C. elegans and OrV by tracking the infection's progression during larval development. Four distinct stages of infection were identified on the basis of viral load, with a peak in capsid-encoding RNA2 coinciding with the first signs of viral egression. Transcriptomic analysis revealed temporal changes in gene expression and functions induced by the infection. A specific set of up-regulated genes remained active throughout the infection, and genes correlated and anticorrelated with virus accumulation were identified. Responses to OrV mirrored reactions to other biotic stressors, distinguishing between virus-specific responses and broader immune responses. Moreover, mutants of early response genes and defense-related processes showed altered viral load progression, uncovering additional players in the antiviral defense response.

Chemically defined production of engineered cardiac tissue microspheres from hydrogel-encapsulated pluripotent stem cells.

Chemically defined, suspension culture conditions are a key requirement in realizing clinical translation of engineered cardiac tissues (ECTs). Building on our previous work producing functional ECT microspheres through differentiation of biomaterial encapsulated human induced pluripotent stem cells (hiPSCs), here we establish the ability to use chemically defined culture conditions, including stem cell media (E8) and cardiac differentiation media (chemically defined differentiation media with three components, CDM3). A custom microfluidic cell encapsulation system was used to encapsulate hiPSCs at a range of initial cell concentrations and diameters in the hybrid biomaterial, poly(ethylene glycol)-fibrinogen (PF), for the formation of highly spherical and uniform ECT microspheres for subsequent cardiac differentiation. Initial microsphere diameter could be tightly controlled, and microspheres could be produced with an initial diameter between 400 and 800 µm. Three days after encapsulation, cardiac differentiation was initiated through small molecule modulation of Wnt signaling in CDM3. Cardiac differentiation occurred resulting in in situ ECTformation; results showed that this differentiation protocol could be used to achieve cardiomyocyte (CM) contents greater than 90%, although there was relatively high variability in CM content and yield between differentiation batches. Spontaneous contraction of ECT microspheres initiated between Days 7 and 10 of differentiation and ECT microspheres responded to electrical pacing up to 1.5 Hz. Resulting CMs had well-defined sarcomeres and the gap junction protein, connexin 43, and had appropriate temporal changes in gene expression. In summary, this study demonstrated the proof-of-concept to produce functional ECT microspheres with chemically defined media in suspension culture in combination with biomaterial support of microsphere encapsulated hiPSCs.

Breast Cancer Cells Exhibit Mesenchymal-Epithelial Plasticity Following Dynamic Modulation of Matrix Stiffness.

Mesenchymal-epithelial transition (MET) is essential for tissue and organ development and is thought to contribute to cancer by enabling the establishment of metastatic lesions. Despite its importance in both health and disease, there is a lack of in vitro platforms to study MET and little is known about the regulation of MET by mechanical cues. Here, hyaluronic acid-based hydrogels with dynamic and tunable stiffnesses mimicking that of normal and tumorigenic mammary tissue are synthesized. The platform is then utilized to examine the response of mammary epithelial cells and breast cancer cells to dynamic modulation of matrix stiffness. Gradual softening of the hydrogels reduces proliferation and increases apoptosis of breast cancer cells. Moreover, breast cancer cells exhibit temporal changes in cell morphology, cytoskeletal organization, and gene expression that are consistent with mesenchymal-epithelial plasticity as the stiffness of the matrix is reduced. A reduction in matrix stiffness attenuates the expression of integrin-linked kinase, and inhibition of integrin-linked kinase impacts proliferation, apoptosis, and gene expression in cells cultured on stiff and dynamic hydrogels. Overall, these findings reveal intermediate epithelial/mesenchymal states as cells move along a matrix stiffness-mediated MET trajectory and suggest an important role for matrix mechanics in regulating mesenchymal-epithelial plasticity.

The link between gene duplication and divergent patterns of gene expression across a complex life cycle.

The diversification of many lineages throughout natural history has frequently been associated with evolutionary changes in life cycle complexity. However, our understanding of the processes that facilitate differentiation in the morphologies and functions expressed by organisms throughout their life cycles is limited. Theory suggests that the expression of traits is decoupled across life stages, thus allowing for their evolutionary independence. Although trait decoupling between stages is well established, explanations of how said decoupling evolves have seldom been considered. Because the different phenotypes expressed by organisms throughout their life cycles are coded for by the same genome, trait decoupling must be mediated through divergence in gene expression between stages. Gene duplication has been identified as an important mechanism that enables divergence in gene function and expression between cells and tissues. Because stage transitions across life cycles require changes in tissue types and functions, we investigated the potential link between gene duplication and expression divergence between life stages. To explore this idea, we examined the temporal changes in gene expression across the monarch butterfly (Danaus plexippus) metamorphosis. We found that within homologous groups, more phylogenetically diverged genes exhibited more distinct temporal expression patterns. This relationship scaled such that more phylogenetically diverse homologous groups showed more diverse patterns of gene expression. Furthermore, we found that duplicate genes showed increased stage-specificity relative to singleton genes. Overall, our findings suggest an important link between gene duplication and the evolution of complex life cycles.

Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus.

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

Evaluating tumor inflammation and PD-L1 expression in sequential biopsies of real-world primary and metastatic breast cancer.

e13035 Background: Biomarkers for measuring immunogenicity of tumors such as PD-L1 immunohistochemistry (IHC), tumor mutational burden (TMB), and immune gene expression have been studied in many cancers including breast cancer. How these biomarkers change over time and how their detection differs with sequential testing in the real-world clinical setting has been studied to a lesser extent but is an important consideration for optimal patient care. Methods: We performed an analysis of 30 patients with primary and metastatic breast cancer sequentially tested by comprehensive genomic and immune profiling (CGIP) at two time points during routine clinical care to observe temporal changes in immunotherapy markers and immune gene expression. Genomic variants and TMB were assessed with DNA-seq, while immune gene expression was assessed with RNA-seq. PD-L1 expression was determined using IHC. Immune gene expression signatures characterizing the tumor microenvironment (TME) were calculated, including tumor inflammation (TIGS), cell proliferation (CP), and cancer testis antigen burden (CTAB). Differences between sequential tests were tested for using linear mixed effect models. Results: Time between biopsies for sequential tests ranged from the same day to 4.5 years, with 40% of specimens being collected within 30 days of the first. Sequential specimens were collected from breast (27%), non-breast (40%), or a mix of breast and non-breast (33%) tissue. Non-breast tissue types included sites from lymph node (18%), GI tract (15%), skin (10%), thorax (7%), and musculoskeletal sites (7%). Overall, TMB, PD-L1 IHC, and gene expression signatures were relatively stable between sequential tests (P > 0.05); however, trends were noted for second tests detecting lower tumor inflammation (TIGS score) and PD-L1 IHC (P = 0.1). Further, second tests detected significantly lower levels of tumor inflammation when both were performed on non-breast tissue (P = 0.04), but not when performed on all breast (P = 0.7) or a mix of breast and non-breast (P = 0.6) tissues. An opposite trend was noted for PD-L1 where second tests detected lower levels of PD-L1 IHC when both tests were performed on breast tissue (P = 0.1), but not when performed on all non-breast (P = 0.4) or a mix of breast and non-breast (P = 0.5) tissue. Tissue-dependent trends were still present even when excluding lymph node specimens (P = 0.1). Conclusions: Changes in immunotherapy markers and immune gene expression signatures, on average, remained relatively stable over time, except for tumor inflammation and PD-L1 IHC when sequentially testing all non-breast or breast tissue, respectively. These observations highlight the utility of clinical CGIP testing in revealing temporal changes of the TME, which may help inform clinical trial or treatment strategies occurring over time for primary and metastatic breast cancer.

Elucidating the ecophysiology of soybean pod-sucking stinkbug Riptortus pedestris (Hemiptera: Alydidae) based on de novo genome assembly and transcriptome analysis

Food security is important for the ever-growing global population. Soybean, Glycine max (L.) Merr., is cultivated worldwide providing a key source of food, protein and oil. Hence, it is imperative to maintain or to increase its yield under different conditions including challenges caused by abiotic and biotic stresses. In recent years, the soybean pod-sucking stinkbug Riptortus pedestris has emerged as an important agricultural insect pest in East, South and Southeast Asia. Here, we present a genomics resource for R. pedestris including its genome assembly, messenger RNA (mRNA) and microRNA (miRNA) transcriptomes at different developmental stages and from different organs. As insect hormone biosynthesis genes (genes involved in metamorphosis) and their regulators such as miRNAs are potential targets for pest control, we analyzed the sesquiterpenoid (juvenile) and ecdysteroid (molting) hormone biosynthesis pathway genes including their miRNAs and relevant neuropeptides. Temporal gene expression changes of these insect hormone biosynthesis pathways were observed at different developmental stages. Similarly, a diet-specific response in gene expression was also observed in both head and salivary glands. Furthermore, we observed that microRNAs (bantam, miR-14, miR-316, and miR-263) of R. pedestris fed with different types of soybeans were differentially expressed in the salivary glands indicating a diet-specific response. Interestingly, the opposite arms of miR-281 (-5p and -3p), a miRNA involved in regulating development, were predicted to target Hmgs genes of R. pedestris and soybean, respectively. These observations among others highlight stinkbug’s responses as a function of its interaction with soybean. In brief, the results of this study not only present salient findings that could be of potential use in pest management and mitigation but also provide an invaluable resource for R. pedestris as an insect model to facilitate studies on plant-pest interactions.

Temporal changes of gene expression in health, schizophrenia, bipolar disorder, and major depressive disorder

The molecular events underlying the development, manifestation, and course of schizophrenia, bipolar disorder, and major depressive disorder span from embryonic life to advanced age. However, little is known about the early dynamics of gene expression in these disorders due to their relatively late manifestation. To address this, we conducted a secondary analysis of post-mortem prefrontal cortex datasets using bioinformatics and machine learning techniques to identify differentially expressed gene modules associated with aging and the diseases, determine their time-perturbation points, and assess enrichment with expression quantitative trait loci (eQTL) genes. Our findings revealed early, mid, and late deregulation of expression of functional gene modules involved in neurodevelopment, plasticity, homeostasis, and immune response. This supports the hypothesis that multiple hits throughout life contribute to disease manifestation rather than a single early-life event. Moreover, the time-perturbed functional gene modules were associated with genetic loci affecting gene expression, highlighting the role of genetic factors in gene expression dynamics and the development of disease phenotypes. Our findings emphasize the importance of investigating time-dependent perturbations in gene expression before the age of onset in elucidating the molecular mechanisms of psychiatric disorders.

Spatiotemporal transcriptomic plasticity in barley roots: unravelling water deficit responses in distinct root zones

BackgroundDrought poses a major threat to agricultural production and thus food security. Understanding the processes shaping plant responses to water deficit is essential for global food safety. Though many studies examined the effect of water deficit on the whole-root level, the distinct functions of each root zone and their specific stress responses remain masked by this approach.ResultsIn this study, we investigated the effect of water deficit on root development of the spring barley (Hordeum vulgare L.) cultivar Morex and examined transcriptomic responses at the level of longitudinal root zones. Water deficit significantly reduced root growth rates after two days of treatment. RNA-sequencing revealed root zone and temporal gene expression changes depending on the duration of water deficit treatment. The majority of water deficit-regulated genes were unique for their respective root zone-by-treatment combination, though they were associated with commonly enriched gene ontology terms. Among these, we found terms associated with transport, detoxification, or cell wall formation affected by water deficit. Integration of weighted gene co-expression analyses identified differential hub genes, that highlighted the importance of modulating energy and protein metabolism and stress response.ConclusionOur findings provide new insights into the highly dynamic and spatiotemporal response cascade triggered by water deficit and the underlying genetic regulations on the level of root zones in the barley cultivar Morex, providing potential targets to enhance plant resilience against environmental constraints. This study further emphasizes the importance of considering spatial and temporal resolution when examining stress responses.

AGE-ASSOCIATED DECLINE IN TISSUE CIRCADIAN CLOCK FUNCTION

Abstract Cellular circadian clocks direct a daily transcriptional program that supports homeostasis and resilience. Emerging evidence has demonstrated age-associated changes in circadian functions. To define age-dependent changes at the systems level, we profile the circadian transcriptome in the hypothalamus, lung, heart, kidney, skeletal muscle, and adrenal gland in three age groups. We find age-dependent and tissue-specific clock output changes but these changes exhibit tissue specific patterns. In general, aging reduces the number of rhythmically expressed genes (REGs), indicative of weakened circadian control. These REGs are enriched for the hallmarks of aging, adding another dimension to our understanding of aging. We determined there were tissue specific changes in the phase distribution of REGs with the hypothalamus and skeletal muscle tissues of the old mice showing a large reduction with REG peak expression at only one phase of the day. Analyzing all expressed genes within a tissue for differences at four distinct times of day identified unique clusters of differentially expressed genes (DEGs). These outcomes identify non-circadian but temporally distinct age associated changes in gene expression. This analysis extends the landscape for understanding aging and highlights the impact of aging on circadian clock function and temporal changes in gene expression.

Transcriptional changes during isoproterenol-induced cardiac fibrosis in mice.

Introduction: β-adrenergic stimulation using β-agonists such as isoproterenol has been routinely used to induce cardiac fibrosis in experimental animal models. Although transcriptome changes in surgical models of cardiac fibrosis such as transverse aortic constriction (TAC) and coronary artery ligation (CAL) are well-studied, transcriptional changes during isoproterenol-induced cardiac fibrosis are not well-explored. Methods: Cardiac fibrosis was induced in male C57BL6 mice by administration of isoproterenol for 4, 8, or 11days at 50mg/kg/day dose. Temporal changes in gene expression were studied by RNA sequencing. Results and discussion: We observed a significant alteration in the transcriptome profile across the different experimental groups compared to the saline group. Isoproterenol treatment caused upregulation of genes associated with ECM organization, cell-cell contact, three-dimensional structure, and cell growth, while genes associated with fatty acid oxidation, sarcoplasmic reticulum calcium ion transport, and cardiac muscle contraction are downregulated. A number of known long non-coding RNAs (lncRNAs) and putative novel lncRNAs exhibited differential regulation. In conclusion, our study shows that isoproterenol administration leads to the dysregulation of genes relevant to ECM deposition and cardiac contraction, and serves as an excellent alternate model to the surgical models of heart failure.

Akaluc/AkaLumine bioluminescence system enables highly sensitive, non-invasive and temporal monitoring of gene expression in Drosophila

Bioluminescence generated by luciferase and luciferin has been extensively used in biological research. However, detecting signals from deep tissues in vivo poses a challenge to traditional methods. To overcome this, the Akaluc and AkaLumine bioluminescent systems were developed, resulting in improved signal detection. We evaluate the potential of Akaluc/AkaLumine in Drosophila melanogaster to establish a highly sensitive, non-invasive, and temporal detection method for gene expression. Our results show that oral administration of AkaLumine to flies expressing Akaluc provided a higher luminescence signal than Luc/D-luciferin, with no observed harmful effects on flies. The Akaluc/AkaLumine system allows for monitoring of dynamic temporal changes in gene expression. Additionally, using the Akaluc fusion gene allows for mRNA splicing monitoring. Our findings indicate that the Akaluc/AkaLumine system is a powerful bioluminescence tool for analyzing gene expression in deep tissues and small numbers of cells in Drosophila.

LP-48: Temporal gene expression changes in xenobiotic-handling genes after constitutive androstane receptor activation in a model of 3D primary human hepatocytes

LP-48: Temporal gene expression changes in xenobiotic-handling genes after constitutive androstane receptor activation in a model of 3D primary human hepatocytes

Hippocampal Inflammation and Gene Expression Changes in Peripheral Lipopolysaccharide Challenged Mice Showing Sickness and Anxiety-Like Behaviors.

Neuroinflammation is often associated with the development of depressive and anxiety disorders. The hippocampus is one of the brain regions affected by inflammation that is associated with these symptoms. However, the mechanism of hippocampal inflammation-induced emotional behavior remains unknown. The aim of this study was to clarify temporal changes in the neuroinflammatory responses in the hippocampus and the response of dentate gyrus (DG) neurons using peripheral lipopolysaccharide (LPS)-challenged mice. LPS administration induced anxiety-like activity in the elevated plus maze test and social interaction test after 24 h, at which time the mice had recovered from sickness behavior. We examined the hippocampal inflammation-related gene expression changes over time. The expression of interleukin-1β (Il1b) and tumor necrosis factor α (Tnfa) was rapidly enhanced and sustained until 24 h after LPS administration, whereas the expression of Il6 was transiently induced at approx. 6 h. IL-6-dependent downstream signaling of transducer and activator of transcription 3 (STAT3) was also activated approx. 3-6 h after LPS treatment. The expression of innate immune genes including interferon-induced transmembrane proteins such as interferon-induced transmembrane protein 1 (Ifitm1) and Ifitm3 and complement factors such as C1qa and C1qb started to increase approx. 6 h and showed sustained or further increase at 24 h. We also examined changes in the expression of several maturation markers in the DG and found that LPS enhanced the expression of calbindin 1 (Calb1), tryptophan-2,3-dioxigenase 2 (Tdo2), Il1rl, and neurotrophin-3 (Ntf3) at 24 h after LPS treatment. Collectively, these results demonstrate temporal changes of inflammation and gene expression in the hippocampus in LPS-induced sickness and anxiety-like behaviors.

Dual RNA-seq identifies genes and pathways modulated during Clostridioides difficile colonization.

The gastrointestinal pathogen Clostridioides difficile is the most common cause of hospital-acquired diarrhea. Bacterial interactions with the gut mucosa are crucial for the establishment of C. difficile infection; however, key infection events like bacterial attachment and gut penetration are still poorly defined. To better understand the initial events that occur when this anaerobe interacts with the human gut epithelium, we employed a dual RNA-sequencing approach to study the bacterial and host transcriptomic profiles during C. difficile infection in a dual environment in vitro human gut model. Temporal changes in gene expression during infection were studied in bacterial and epithelial cells over 3-24 hours. While there were several common differentially expressed bacterial genes across different timepoints after infection, mammalian transcriptional profiles were quite distinct, with little overlap. Interestingly, an induction of colonic receptors for C. difficile toxins was observed, along with the downregulation of genes encoding immune response markers. Several cell wall-associated proteins were downregulated in C. difficile when associated with host cells, including slpA, which encodes the main S-layer protein. Gene function and pathway enrichment analyses revealed a potential modulation of the purine/pyrimidine synthesis pathways both in the mammalian and bacterial cells. We observed that proline-proline endopeptidase, a secreted metalloprotease responsible for cell surface protein cleavage, is downregulated during infection, and a mutant lacking this enzyme demonstrated enhanced adhesion to epithelial cells during infection. This study provides new insight into the host and bacterial pathways based on gene expression modulation during the initial contact of C. difficile with gut cells. IMPORTANCE The initial interactions between the colonic epithelium and the bacterium are likely critical in the establishment of Clostridioides difficile infection, one of the major causes of hospital-acquired diarrhea worldwide. Molecular interactions between C. difficile and human gut cells have not been well defined mainly due to the technical challenges of studying cellular host-pathogen interactions with this anaerobe. Here we have examined transcriptional changes occurring in the pathogen and host cells during the initial 24 hours of infection. Our data indicate several changes in metabolic pathways and virulence-associated factors during the initial bacterium-host cell contact and early stages of infection. We describe canonical pathways enriched based on the expression profiles of a dual RNA sequencing in the host and bacterium, and functions of bacterial factors that are modulated during infection. This study thus provides fresh insight into the early C. difficile infection process.

Post-COVID symptoms are associated with endotypes reflecting poor inflammatory and hemostatic modulation.

Persistent symptoms after COVID-19 infection ("long COVID") negatively affects almost half of COVID-19 survivors. Despite its prevalence, its pathophysiology is poorly understood, with multiple host systems likely affected. Here, we followed patients from hospital to discharge and used a systems-biology approach to identify mechanisms of long COVID. RNA-seq was performed on whole blood collected early in hospital and 4-12 weeks after discharge from 24 adult COVID-19 patients (10 reported post-COVID symptoms after discharge). Differential gene expression analysis, pathway enrichment, and machine learning methods were used to identify underlying mechanisms for post-COVID symptom development. Compared to patients with post-COVID symptoms, patients without post-COVID symptoms had larger temporal gene expression changes associated with downregulation of inflammatory and coagulation genes over time. Patients could also be separated into three patient endotypes with differing mechanistic trajectories, which was validated in another published patient cohort. The "Resolved" endotype (lowest rate of post-COVID symptoms) had robust inflammatory and hemostatic responses in hospital that resolved after discharge. Conversely, the inflammatory/hemostatic responses of "Suppressive" and "Unresolved" endotypes (higher rates of patients with post-COVID symptoms) were persistently dampened and activated, respectively. These endotypes were accurately defined by specific blood gene expression signatures (6-7 genes) for potential clinical stratification. This study allowed analysis of long COVID whole blood transcriptomics trajectories while accounting for the issue of patient heterogeneity. Two of the three identified and externally validated endotypes ("Unresolved" and "Suppressive") were associated with higher rates of post-COVID symptoms and either persistently activated or suppressed inflammation and coagulation processes. Gene biomarkers in blood could potentially be used clinically to stratify patients into different endotypes, paving the way for personalized long COVID treatment.

A TRANSCRIPTOMIC APPRECIATION OF CHILDHOOD MENINGOCOCCAL AND POLYMICROBIAL SEPSIS FROM A PRO-INFLAMMATORY AND TRAJECTORIAL PERSPECTIVE, A ROLE FOR VASCULAR ENDOTHELIAL GROWTH FACTOR A AND B MODULATION?

This study investigated the temporal dynamics of childhood sepsis by analyzing gene expression changes associated with proinflammatory processes. Five datasets, including four meningococcal sepsis shock (MSS) datasets (two temporal and two longitudinal) and one polymicrobial sepsis dataset, were selected to track temporal changes in gene expression. Hierarchical clustering revealed three temporal phases: early, intermediate, and late, providing a framework for understanding sepsis progression. Principal component analysis supported the identification of gene expression trajectories. Differential gene analysis highlighted consistent upregulation of vascular endothelial growth factor A (VEGF-A) and nuclear factor κB1 (NFKB1), genes involved in inflammation, across the sepsis datasets. NFKB1 gene expression also showed temporal changes in the MSS datasets. In the postmortem dataset comparing MSS cases to controls, VEGF-A was upregulated and VEGF-B downregulated. Renal tissue exhibited higher VEGF-A expression compared with other tissues. Similar VEGF-A upregulation and VEGF-B downregulation patterns were observed in the cross-sectional MSS datasets and the polymicrobial sepsis dataset. Hexagonal plots confirmed VEGF-R (VEGF receptor)-VEGF-R2 signaling pathway enrichment in the MSS cross-sectional studies. The polymicrobial sepsis dataset also showed enrichment of the VEGF pathway in septic shock day 3 and sepsis day 3 samples compared with controls. These findings provide unique insights into the dynamic nature of sepsis from a transcriptomic perspective and suggest potential implications for biomarker development. Future research should focus on larger-scale temporal transcriptomic studies with appropriate control groups and validate the identified gene combination as a potential biomarker panel for sepsis.

Dynamic nucleosome remodeling mediated by YY1 underlies early mouse development.

Nucleosome positioning can alter the accessibility of DNA-binding proteins to their cognate DNA elements, and thus its precise control is essential for cell identity and function. Mammalian preimplantation embryos undergo temporal changes in gene expression and cell potency, suggesting the involvement of dynamic epigenetic control during this developmental phase. However, the dynamics of nucleosome organization during early development are poorly understood. In this study, using a low-input MNase-seq method, we show that nucleosome positioning is globally obscure in zygotes but becomes well defined during subsequent development. Down-regulation of the chromatin assembly in embryonic stem cells can partially reverse nucleosome organization into a zygote-like pattern, suggesting a possible link between the chromatin assembly pathway and fuzzy nucleosomes in zygotes. We also reveal that YY1, a zinc finger-containing transcription factor expressed upon zygotic genome activation, regulates the de novo formation of well-positioned nucleosome arrays at the regulatory elements through identifying YY1-binding sites in eight-cell embryos. The YY1-binding regions acquire H3K27ac enrichment around the eight-cell and morula stages, and YY1 depletion impairs the morula-to-blastocyst transition. Thus, our study delineates the remodeling of nucleosome organization and its underlying mechanism during early mouse development.

Transcriptomic forecasting with neural ordinary differential equations

Single-cell transcriptomics technologies can uncover changes in the molecular states that underlie cellular phenotypes. However, understanding the dynamic cellular processes requires extending from inferring trajectories from snapshots of cellular states to estimating temporal changes in cellular gene expression. To address this challenge, we have developed a neural ordinary differential-equation-based method, RNAForecaster, for predicting gene expression states in single cells for multiple future time steps in an embedding-independent manner. We demonstrate that RNAForecaster can accurately predict future expression states in simulated single-cell transcriptomic data with cellular tracking over time. We then show that by using metabolic labeling single-cell RNA sequencing (scRNA-seq) data from constitutively dividing cells, RNAForecaster accurately recapitulates many of the expected changes in gene expression during progression through the cell cycle over a 3-day period. Thus, RNAForecaster enables short-term estimation of future expression states in biological systems from high-throughput datasets with temporal information.

Plasma corticosterone rhythmicity loss dysregulates the renal artery circadian transcriptome and abolishes daily vasodilation variation

Molecular clocks control a daily transcriptional schedule that maintains tissue homeostasis. Glucocorticoids synchronize peripheral clocks with the master clock in the suprachiasmatic nucleus in response to ambient light. In humans and mice, arrhythmic glucocorticoids induce non-dipping blood pressure and vascular dysfunction. The mechanisms of this are poorly understood. We hypothesize that arrhythmic activation of the glucocorticoid receptor desynchronizes the molecular clock mechanism in the renal artery, a major component of blood pressure control, dysregulating normal function. Our aim was to dissect renal artery circadian physiology transcriptionally and functionally and interrogate the effect of arrhythmic glucocorticoids. Male mice, kept on a 12:12 light:dark cycle, were anesthetised with 4% isolflurane and implanted with a hypodermal pellet containing either vehicle or corticosterone (≈3.7 mg/kg per day), to flatten the endogenous glucocorticoid rhythm. After 7 days, mice were culled by cervical dislocation at 2hr intervals for 48 hours and the renal arteries taken. The renal artery circadian transcriptome was profiled using Illumina Truseq Stranded RNA-sequencing. Cosinor regression (Limorhyde) was applied to the control group and identified 465 of 14425 (~3%) protein coding transcripts with rhythmic expression (p<0.005 for fit to a cosinor waveform). Of these, 265 (62%) lost their rhythm in the corticosterone-treated group. These arrhythmic transcripts were interrogated by gene set enrichment analysis (Biological Process, FDR<0.05) and enriched gene sets were related to “Circadian Rhythm,” “Transforming growth factor beta” and “Smooth muscle proliferation.” Paradoxically, 919 transcripts gained de novo rhythmicity in the corticosterone group and were related to mitochondrial respiration and ATP synthesis. In parallel, we used wire myography to determine the diurnal variation in endothelium-independent vasodilation via sodium nitroprusside dose response. Renal arteries were taken from corticosterone and vehicle treated mice as described above, either during the active phase (subjective night) or inactive phase (subjective day). Data are mean ± SD and analysed by 2-way ANOVA. Vasodilation was higher during the night compared to the day in the control group, (night logEC50 =8.5±0.4 versus day logEC50 =7.5±0.4 p=0.01, n=12). In the corticosterone treated group circadian variability in vasodilation was blunted (night logEC50 =7.7±0.9 versus day logEC50=7.7±0.4, p=0.9, n=12). We speculate that temporal misalignment of the vascular molecular circadian machinery and the subsequent temporal changes in gene expression may have a direct effect on dilatory capacity. This novel study extends our understanding of circadian vascular physiology, undelining the impact of glucocorticoids on temporal gene expression.This may be clinically relevant in the pathogenesis of vascular dysfunction linked to elevated glucocorticoids in metabolic syndrome. This study was funded by the British Heart Foundation (PhD studentship code: FS/18/57/34178) and the Kidney Research UK (Intermediate Grant Code: INT_001_20181115 and IN_002_2021 0729). Analyses were carried out by Edinburgh Genomics at the University of Edinburgh. Edinburgh Genomics is partly supported through core grants from NERC (R8/H10/56), MRC (MR/K001744/1) and BBSRC (BB/J004243/1). This is the full abstract presented at the American Physiology Summit 2023 meeting and is only available in HTML format. There are no additional versions or additional content available for this abstract. Physiology was not involved in the peer review process.