Catalytically distinct IDH1 mutants tune phenotype severity in tumor models.

Mutations in isocitrate dehydrogenase 1 (IDH1) impart a neomorphic reaction that produces the oncometabolite D-2-hydroxyglutarate (D2HG), which can inhibit DNA and histone demethylases to drive tumorigenesis via epigenetic changes. Though heterozygous point mutations in patients primarily affect residue R132, there are myriad D2HG-producing mutants that display unique catalytic efficiency of D2HG production. Here, we show that catalytic efficiency of D2HG production is greater in IDH1 R132Q than R132H mutants, and expression of IDH1 R132Q in cellular and mouse xenograft models leads to higher D2HG concentrations in cells, tumors, and sera compared to R132H-expressing models. Reduced representation bisulfite sequencing (RRBS) analysis of xenograft tumors shows expression of IDH1 R132Q relative to R132H leads to hypermethylation patterns in pathways associated with DNA damage. Transcriptome analysis indicates that the IDH1 R132Q mutation has a more aggressive pro-tumor phenotype, with members of EGFR, Wnt, and PI3K signaling pathways differentially expressed, perhaps through non-epigenetic routes. Together, these data suggest that the catalytic efficiency of IDH1 mutants modulate D2HG levels in cellular and in vivo models, resulting in unique epigenetic and transcriptomic consequences where higher D2HG levels appear to be associated with more aggressive tumors.

danRerLib: a Python package for zebrafish transcriptomics.

Summary: Understanding the pathways and biological processes underlying differential gene expression is fundamental for characterizing gene expression changes in response to an experimental condition. Zebrafish, with a transcriptome closely mirroring that of humans, are frequently utilized as a model for human development and disease. However, a challenge arises due to the incomplete annotations of zebrafish pathways and biological processes, with more comprehensive annotations existing in humans. This incompleteness may result in biased functional enrichment findings and loss of knowledge. danRerLib, a versatile Python package for zebrafish transcriptomics researchers, overcomes this challenge and provides a suite of tools to be executed in Python including gene ID mapping, orthology mapping for the zebrafish and human taxonomy, and functional enrichment analysis utilizing the latest updated Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. danRerLib enables functional enrichment analysis for GO and KEGG pathways, even when they lack direct zebrafish annotations through the orthology of human-annotated functional annotations. This approach enables researchers to extend their analysis to a wider range of pathways, elucidating additional mechanisms of interest and greater insight into experimental results. Availability and implementation: danRerLib, along with comprehensive documentation and tutorials, is freely available. The source code is available at https://github.com/sdsucomptox/danrerlib/ with associated documentation and tutorials at https://sdsucomptox.github.io/danrerlib/. The package has been developed with Python 3.9 and is available for installation on the package management systems PIP (https://pypi.org/project/danrerlib/) and Conda (https://anaconda.org/sdsu_comptox/danrerlib) with additional installation instructions on the documentation website.

Development of a Dynamic Network Model to Identify Temporal Patterns of Structural Malformations in Zebrafish Embryos Exposed to a Model Toxicant, Tris(4-chlorophenyl)methanol.

Embryogenesis is a well-coordinated process relying on precise cues and environmental signals that direct spatiotemporal embryonic patterning. Quite often, when one error in this process occurs, others tend to co-occur. We posit that investigating the co-occurrence of these abnormalities over time would yield additional information about the mode of toxicity for chemicals. Here, we use the environmental contaminant tris(4-chlorophenyl)methanol (TCPMOH) as a model toxicant to assess the relationship between exposures and co-occurrence of developmental abnormalities in zebrafish embryos. We propose a dynamic network modeling approach to study the co-occurrence of abnormalities, including pericardial edema, yolk sac edema, cranial malformation, spinal deformity, delayed/failed swim bladder inflation, and mortality induced by TCPMOH exposure. TCPMOH-exposed samples revealed increased abnormality co-occurrence when compared to controls. The abnormalities were represented as nodes in the dynamic network model. Abnormalities with high co-occurrence over time were identified using network centrality scores. We found that the temporal patterns of abnormality co-occurrence varied between exposure groups. In particular, the high TCPMOH exposure group experienced abnormality co-occurrence earlier than the low exposure group. The network model also revealed that pericardial and yolk sac edema are the most common critical nodes among all TCPMOH exposure levels, preceding further abnormalities. Overall, this study introduces a dynamic network model as a tool for assessing developmental toxicology, integrating structural and temporal features with a concentration response.

Effect of circadian system disruption on the concentration and daily oscillations of cortisol, progesterone, melatonin, serotonin, growth hormone, and core body temperature in periparturient dairy cattle.

Metabolic, circadian, sleep, and reproductive systems are integrated and reciprocally regulated, but the understanding of the mechanism is limited. To study this integrated regulation, the circadian timing system was disrupted by exposing late pregnant nonlactating (dry) cows to chronic shifts in the light-dark phase, and rhythms of body temperature and circulating cortisol (CORT), progesterone (P4), serotonin (5HT), melatonin (MEL), and growth hormone (GH) concentrations were measured. Specifically, across 2 identical studies (1 and 2), at 35 d before expected calving (BEC) multiparous cows were assigned to control (CON; n = 24) and exposed to 16 h light and 8 h dark or phase shift (PS; n = 24) treatments and exposed to 6-h light-dark phase shifts every 3 d until parturition. All cows were exposed to control lighting after calving. Blood samples were collected in the first study at 0600 h on d 35 BEC, d 21 BEC, and 2 d before calving, and d 0, 2, 9, 15, and 22 postpartum (PP). A subset of cows (n = 6/group) in study 1 was blood sampled every 4 h over 48 h beginning on d 23 BEC, 9 BEC, and 5 PP. Body temperature was measured every 30 min (n = 8-16/treatment) for 48 h at 23 BEC and 9 BEC in both studies; and at 14 PP and 60 PP only in study 2. Treatment did not affect levels of CORT, GH, or P4 at 0600 h, but overall level of 5HT was lower and MEL higher in PS cows across days sampled. A 2-component versus single-component cosinor model better described [>coefficient of determination (R2);

The Axes of Life: A Roadmap for Understanding Dynamic Multiscale Systems.

The biological challenges facing humanity are complex, multi-factorial, and are intimately tied to the future of our health, welfare, and stewardship of the Earth. Tackling problems in diverse areas, such as agriculture, ecology, and health care require linking vast datasets that encompass numerous components and spatio-temporal scales. Here, we provide a new framework and a road map for using experiments and computation to understand dynamic biological systems that span multiple scales. We discuss theories that can help understand complex biological systems and highlight the limitations of existing methodologies and recommend data generation practices. The advent of new technologies such as big data analytics and artificial intelligence can help bridge different scales and data types. We recommend ways to make such models transparent, compatible with existing theories of biological function, and to make biological data sets readable by advanced machine learning algorithms. Overall, the barriers for tackling pressing biological challenges are not only technological, but also sociological. Hence, we also provide recommendations for promoting interdisciplinary interactions between scientists.

Mathematical modeling of the interaction between yolk utilization and fish growth in zebrafish, Danio rerio.

Optimal embryonic development plays a major role in the health of an individual beyond the developmental stage. Nutritional perturbation during development is associated with cardiovascular and metabolic disease later in life. With both nutritional uptake and overall growth being risk factors for eventual health, it is necessary to understand not only the behavior of the processes during development but also their interactions. In this study, we used differential equations, image analyses, curve fittings, parameter estimation and laboratory experiments to quantify the rate of yolk absorption and its effect on early development of a vertebrate model (Danio rerio). Findings from this study establish a nonlinear functional relationship between nutrient absorption and early fish growth. We found that the rate of change in fish length and yolk utilization is logistic, that is the yolk decays rapidly for a period of time before leveling out. An interesting finding from this study is that yolk utilization reaches its maximum at 84 h post-fertilization. We validated our mathematical models against experimental observations, making them powerful tools for replication and future simulations.

The ecotoxicological contaminant tris(4-chlorophenyl)methanol (TCPMOH) impacts embryonic development in zebrafish (Danio rerio).

Tris(4-chlorophenyl)methanol (TCPMOH) is a water contaminant with unknown etiology, but is believed to be a byproduct of DDT manufacturing. It is highly persistent in the environment, and bioaccumulates in marine species. TCPMOH has also been measured in human breast milk, which poses a risk for developing infants. However, almost no toxicity data is currently available. In this study, we investigate the hazard posed by developmental TCPMOH exposures using the zebrafish model (Danio rerio). Zebrafish (Danio rerio) embryos were exposed to 0, 0.1, 0.5, 1, or 5 µM TCPMOH beginning at 24 h post fertilization (hpf). Embryonic mortality and incidence of morphological deformities increased in a concentration-dependent manner with TCPMOH exposure. RNA sequencing assessed changes in gene expression associated with acute (4 hour) exposures to 50 nM TCPMOH. Developmental exposure to TCPMOH decreased expression of ahr2, as well as metabolic enzymes cyp1a1, cyp1b1, cyp1c1, cyp1c2, and cyp2y3 (p<0.05). These findings were concordant with decreased Cyp1a1 induction measured by the ethoxyresorufin-O-deethylase (EROD) assay (p<0.05). Pathways associated with xenobiotic metabolism, lipid metabolism, and transcriptional and translational regulation were decreased. Pathways involved in DNA replication and repair, carbohydrate metabolism, and endocrine function were upregulated. Overall, this study demonstrates that TCPMOH is acutely toxic to zebrafish embryos at elevated concentrations.

Extraction and Analysis of Respiratory Motion Using a Comprehensive Wearable Health Monitoring System.

Respiratory activity is an important vital sign of life that can indicate health status. Diseases such as bronchitis, emphysema, pneumonia and coronavirus cause respiratory disorders that affect the respiratory systems. Typically, the diagnosis of these diseases is facilitated by pulmonary auscultation using a stethoscope. We present a new attempt to develop a lightweight, comprehensive wearable sensor system to monitor respiration using a multi-sensor approach. We employed new wearable sensor technology using a novel integration of acoustics and biopotentials to monitor various vital signs on two volunteers. In this study, a new method to monitor lung function, such as respiration rate and tidal volume, is presented using the multi-sensor approach. Using the new sensor, we obtained lung sound, electrocardiogram (ECG), and electromyogram (EMG) measurements at the external intercostal muscles (EIM) and at the diaphragm during breathing cycles with 500 mL, 625 mL, 750 mL, 875 mL, and 1000 mL tidal volume. The tidal volumes were controlled with a spirometer. The duration of each breathing cycle was 8 s and was timed using a metronome. For each of the different tidal volumes, the EMG data was plotted against time and the area under the curve (AUC) was calculated. The AUC calculated from EMG data obtained at the diaphragm and EIM represent the expansion of the diaphragm and EIM respectively. AUC obtained from EMG data collected at the diaphragm had a lower variance between samples per tidal volume compared to those monitored at the EIM. Using cubic spline interpolation, we built a model for computing tidal volume from EMG data at the diaphragm. Our findings show that the new sensor can be used to measure respiration rate and variations thereof and holds potential to estimate tidal lung volume from EMG measurements obtained from the diaphragm.

Acetylcholine Regulates Pulmonary Pathology During Viral Infection and Recovery.

INTRODUCTION: This study was designed to explore the role of acetylcholine (ACh) in pulmonary viral infection and recovery. Inflammatory control is critical to recovery from respiratory viral infection. ACh secreted from non-neuronal sources, including lymphocytes, plays an important, albeit underappreciated, role in regulating immune-mediated inflammation. METHODS: ACh and lymphocyte cholinergic status in the lungs were measured over the course of influenza infection and recovery. The role of ACh was examined by inhibiting ACh synthesis in vivo. Pulmonary inflammation was monitored by Iba1 immunofluorescence, using a novel automated algorithm. Tissue repair was monitored histologically. RESULTS: Pulmonary ACh remained constant through the early stage of infection and increased during the peak of the acquired immune response. As the concentration of ACh increased, cholinergic lymphocytes appeared in the BAL and lungs. Cholinergic capacity was found primarily in CD4 T cells, but also in B cells and CD8 T cells. The cholinergic CD4+ T cells bound to influenza-specific tetramers and were retained in the resident memory regions of the lung up to 2 months after infection. Histologically, cholinergic lymphocytes were found in direct physical contact with activated macrophages throughout the lung. Inflammation was monitored by ionized calcium-binding adapter molecule 1 (Iba1) immunofluorescence, using a novel automated algorithm. When ACh production was inhibited, mice exhibited increased tissue inflammation and delayed recovery. Histologic examination revealed abnormal tissue repair when ACh was limited. CONCLUSION: These findings point to a previously unrecognized role for ACh in the transition from active immunity to recovery and pulmonary repair following respiratory viral infection.

Tissue geometry may govern lung branching mode selection.

Lung branching morphogenesis proceeds in three stereotyped modes (domain, planar, and orthogonal branching). Much is known about the molecular players, including growth factors such as fibroblast growth factor 10 but it is unknown how these signals could actuate the different branching patterns. With the aim of identifying mechanisms that may determine the different branching modes, we developed a computational model of the epithelial lung bud and its surrounding mesenchyme. We studied transport of morphogens and localization of morphogen flux at lobe surfaces and lobe edges. We find that a single simple mechanism is theoretically capable of directing an epithelial tubule to elongate, bend, flatten, or bifurcate, depending solely on geometric ratios of the tissues in the vicinity of a growing tubule tip. Furthermore, the same simple mechanism is capable of generating orthogonal or planar branching, depending only on the same geometric ratios.

Quantifying stretch and secretion in the embryonic lung: Implications for morphogenesis.

Branching in the embryonic lung is controlled by a variety of morphogens. Mechanics is also believed to play a significant role in lung branching. The relative roles and interactions of these two broad factors are challenging to determine. We considered three hypotheses for explaining why tracheal occlusion triples branching with no overall increase in size. Both hypotheses are based on tracheal occlusion blocking the exit of secretions. (H1) Increased lumen pressure stretches tissues; stretch receptors at shoulders of growing tips increase local rate of branching. (H2) Blocking exit of secretions blocks advective transport of morphogens, leading to (H2a) increased overall concentration of morphogens or (H2b) increased flux of morphogens at specific locations. We constructed and analyzed computational models of tissue stretch and solute transport in a 3D lung geometry. Observed tissue stresses and stretches were predominantly in locations unrelated to subsequent branch locations, suggesting that tissue stretch (H1) is not the mechanism of enhancement of branching. Morphogen concentration in the mesenchyme (H2a) increased with tracheal occlusion, consistent with previously reported results. Morphogen flux at the epithelial surface (H2b) completely changed its distribution pattern when the trachea was occluded, tripling the number of locations at which it was elevated. Our results are consistent with the hypothesis that tracheal occlusion blocks outflow of secretions, leading to a higher number of high-flux locations at branching tips, in turn leading to a large increase in number of branching locations.

Numerical analysis of the effect of T-tubule location on calcium transient in ventricular myocytes.

Intracellular calcium (Ca2+) signaling in cardiac myocytes is vital for proper functioning of the heart. Understanding the intracellular Ca2+ dynamics would give an insight into the functions of normal and diseased hearts. In the current study, spatiotemporal Ca2+ dynamics is investigated in ventricular myocytes by considering Ca2+ release and re-uptake via sarcolemma and transverse tubules (T-tubules), Ca2+ diffusion and buffering in the cytosol, and the blockade of Ca2+ activities associated with the sarcoplasmic reticulum. This study is carried out using a three dimensional (3D) geometric model of a branch of T-tubule extracted from the electron microscopy (EM) images of a partial ventricular myocyte. Mathematical modeling is done by using a system of partial differential equations involving Ca2+, buffers, and membrane channels. Numerical simulation results suggest that a lack of T-tubule structure at the vicinity of the cell surface could increase the peak time of Ca2+ concentration in myocytes. The results also show that T-tubules and mobile buffers play an important role in the regulation of Ca2+ transient in ventricular myocytes.

Mathematical modelling and numerical simulations of actin dynamics in the eukaryotic cell.

The aim of this article is to study cell deformation and cell movement by considering both the mechanical and biochemical properties of the cortical network of actin filaments and its concentration. Actin is a polymer that can exist either in filamentous form (F-actin) or in monometric form (G-actin) (Chen et al. in Trends Biochem Sci 25:19-23, 2000) and the filamentous form is arranged in a paired helix of two protofilaments (Ananthakrishnan et al. in Recent Res Devel Biophys 5:39-69, 2006). By assuming that cell deformations are a result of the cortical actin dynamics in the cell cytoskeleton, we consider a continuum mathematical model that couples the mechanics of the network of actin filaments with its bio-chemical dynamics. Numerical treatment of the model is carried out using the moving grid finite element method (Madzvamuse et al. in J Comput Phys 190:478-500, 2003). Furthermore, by assuming slow deformations of the cell, we use linear stability theory to validate the numerical simulation results close to bifurcation points. Far from bifurcation points, we show that the mathematical model is able to describe the complex cell deformations typically observed in experimental results. Our numerical results illustrate cell expansion, cell contraction, cell translation and cell relocation as well as cell protrusions. In all these results, the contractile tonicity formed by the association of actin filaments to the myosin II motor proteins is identified as a key bifurcation parameter.