With the development of transcriptomic technologies, we are able to quantify precise changes in gene expression profiles from astronauts and other organisms exposed to spaceflight. Members of NASA GeneLab and GeneLab-associated analysis working groups (AWGs) have developed a consensus pipeline for analyzing short-read RNA-sequencing data from spaceflight-associated experiments. The pipeline includes quality control, read trimming, mapping, and gene quantification steps, culminating in the detection of differentially expressed genes. This data analysis pipeline and the results of its execution using data submitted to GeneLab are now all publicly available through the GeneLab database. We present here the full details and rationale for the construction of this pipeline in order to promote transparency, reproducibility, and reusability of pipeline data; to provide a template for data processing of future spaceflight-relevant datasets; and to encourage cross-analysis of data from other databases with the data available in GeneLab.

To understand the physiological changes that occur in response to spaceflight, mice are transported to the International Space Station (ISS) and housed for variable periods of time before euthanasia on-orbit or return to Earth. Sample collection under such difficult conditions introduces confounding factors that need to be identified and addressed. We found large changes in the transcriptome of mouse tissues dissected and preserved on-orbit compared with tissues from mice euthanized on-orbit, preserved, and dissected after return to Earth. Changes due to preservation method eclipsed those between flight and ground samples, making it difficult to identify spaceflight-specific changes. Follow-on experiments to interrogate the roles of euthanasia methods, tissue and carcass preservation protocols, and library preparation methods suggested that differences due to preservation protocols are exacerbated when coupled with polyA selection. This has important implications for the interpretation of existing datasets and the design of future experiments.

Although antisense oligonucleotides (ASOs) are well tolerated preclinically and in the clinic, some sequences of ASOs can trigger an inflammatory response leading to B cell and macrophage activation in rodents. This prompted our investigation into the contribution of genetic architecture to the ASO-mediated inflammatory response. Genome-wide association (GWA) and transcriptomic analysis in a hybrid mouse diversity panel (HMDP) were used to identify and validate novel genes involved in the acute and delayed inflammatory response to a single 75 mg/kg dose of an inflammatory 2'-O-methoxyethyl (2'MOE) modified ASO. The acute response was measured 6 h after ASO administration, via evaluation for increased plasma production of interleukin 6 (IL6), IL10, monocyte chemoattractant protein 1 (MCP-1) and macrophage inflammatory protein-1β (MIP-1β). Delayed inflammation was evaluated by spleen weight increases after 96 h. We identified single nucleotide polymorphisms (SNPs) on chromosomes 16 and 17 associated with plasma MIP-1β, IL6, and MCP-1 levels, and one on chromosome 8 associated with increases in spleen weight. Systems genetic analysis utilizing transcriptomic data from HMDP strain macrophages determined that the acute inflammatory SNPs were expression quantitative trait locis (eQTLs) for CCAAT/enhancer-binding protein beta (Cebpb) and salt inducible kinase 1 (Sik1). The delayed inflammatory SNP was an eQTL for Rho guanine nucleotide exchange factor 10 (Arhgef10). In vitro assays in mouse primary cells and human cell lines have confirmed the HMDP finding that lower Sik1 expression increases the acute inflammatory response. Our results demonstrate the utility of using mouse GWA study (GWAS) and the HMDP for detecting genes modulating the inflammatory response to pro-inflammatory ASOs in a pharmacological setting.

Abstract Tumor irradiation during cancer treatment induces cell death by apoptosis, mitotic death, and clonogenic death. With expanded global access to cancer radiotherapy, there is a need to discover agents with anti-cancer and radiation sensitizing properties. in vitro high-throughput screening (HTS) systems routinely used to screen agents for cytotoxicity assays unfortunately do not directly measure clonogenic potential and fail to accurately predict the efficacy of an agent either in subsequent preclinical animal model testing or in clinical trials. Clonogenic assays unlike colorimetric and metabolic assays for determining cell viability and apoptosis, are less susceptible to artifacts. We are addressing these requirements through an integrated I-PARTS platform development and optimization. Here, the design and the prototype of a fully enclosed imaging/X-ray irradiation system will be presented. The proof of principle experimentation to obtain high content information on the cellular response to specific primary human cells (normal or tumors) to a variety of radiation doses and fraction modality will be described. The initial approaches taken towards the establishment of an integrated platform to screen anti-cancer drugs will be discussed. Citation Format: Rao V. Papineni, Sylvain Costes, Shayoni Ray, William McLaughlin, Brendan Lukaniec, Kurt Kruger, Keith Kanugi. Design and modular development of I-PARTS: An integrated platform for anti-cancer radiation therapeutics screening [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr LB-070.

MOTIVATION To curate and organize expensive spaceflight experiments conducted aboard space stations and maximize the scientific return of investment, while democratizing access to vast amounts of spaceflight related omics data generated from several model organisms. RESULTS The GeneLab Data System (GLDS) is an open access database containing fully coordinated and curated 'omics' (genomics, transcriptomics, proteomics, metabolomics) data, detailed metadata and radiation dosimetry for a variety of model organisms. GLDS is supported by an integrated data system allowing federated search across several public bioinformatics repositories. Archived datasets can be queried using full-text search (e.g. keywords, Boolean and wildcards) and results can be sorted in multifactorial manner using assistive filters. GLDS also provides a collaborative platform built on GenomeSpace for sharing files and analyses with collaborators. It currently houses 172 datasets and supports standard guidelines for submission of datasets, MIAME (for microarray), ENCODE Consortium Guidelines (for RNA-seq) and MIAPE Guidelines (for proteomics). AVAILABILITY AND IMPLEMENTATION https://genelab.nasa.gov/.

In this manuscript we address the challenges associated with the ability to predict radiation sensitivity associated with exposure to either cosmic radiation or X-rays in a population study, by monitoring DNA damage sensing protein 53BP1 forming small nuclear radiation-induced foci (RIF) as a surrogate biomarker of DNA double strand breaks (DSB). 76 primary skin fibroblasts were isolated from 10 collaborative cross strains and five reference inbred mice (C57Bl/6, BALB/CByJ, B6C3, C3H and CBA/CaJ) and exposed to three different charged nuclei of increasing LET (350 MeV/n Si, 350 MeV/n Ar and 600 MeV/n Fe) and X-ray. Our data brings strong evidence against the classic contact-first model where DSBs are assumed to be immobile and repaired at the lesion site. In contrast, our model suggests nearby DSBs move into single repair unit characterized by large RIF before the repair machinery kicks in. Such model has the advantage of being much more efficient molecularly but is poorly suited to deal with cosmic radiation, where energy is concentrated along the particle trajectory, inducing a large density of DSBs along each particle track. In accordance with this model, RIF quantification after X-ray exposition showed a saturated dose response for early time points post-irradiation for all strains. Similarly, the high-LET response showed that RIF number matched the number of track per cell, not the number of expected DSB per cell (1). At the temporal level, we noted that the percentage of unrepaired high-LET tracks over a 48 hour time-course increased with LET, confirming that the DNA repair process becomes more difficult as more DSB coalesce into single RIF. There was also good agreement between persistent RIF levels measured in-vitro in the primary skin cultures and survival levels of T-cells and B-cells collected in blood samples from 10 CC strains 24 hours after 0.1 Gy whole-body dose of X-ray. This suggests that persistent RIF 24 hour post-IR is a good surrogate in-vitro biomarker for in-vivo radiation toxicity. Finally, at the genomic level, large differences in repair rates between strains for high-LET allowed us to identify suggestive genetic loci associated with radiation sensitivity. Interestingly, the two highest LETs provided the most strain variation with a common locus on Chromosome 10 highly enriched for DNA repair associated genes we discussed in detail.

Translating fundamental biological discoveries from NASA Space Biology program into health risk from space flights has been an ongoing challenge. We propose to use NASA GeneLab database to gain new knowledge on potential systemic responses to space. Unbiased systems biology analysis of transcriptomic data from seven different rodent datasets reveals for the first time the existence of potential "master regulators" coordinating a systemic response to microgravity and/or space radiation with TGF-β1 being the most common regulator. We hypothesized the space environment leads to the release of biomolecules circulating inside the blood stream. Through datamining we identified 13 candidate microRNAs (miRNA) which are common in all studies and directly interact with TGF-β1 that can be potential circulating factors impacting space biology. This study exemplifies the utility of the GeneLab data repository to aid in the process of performing novel hypothesis-based research.

Antisense oligonucleotides (ASOs) have demonstrated variation of efficacy in patient populations. This has prompted our investigation into the contribution of genetic architecture to ASO pharmacokinetics (PK) and pharmacodynamics (PD). Genome wide association (GWA) and transcriptomic analysis in a hybrid mouse diversity panel (HMDP) were used to identify and validate novel genes involved in the uptake and efficacy of a single dose of a Malat1 constrained ethyl (cEt) modified ASO. The GWA of the HMDP identified two significant associations on chromosomes 4 and 10 with hepatic Malat1 ASO concentrations. Stabilin 2 (Stab2) and vesicle associated membrane protein 3 (Vamp3) were identified by cis-eQTL analysis. HMDP strains with lower Stab2 expression and Stab2 KO mice displayed significantly lower PK than strains with higher Stab2 expression and the wild type (WT) animals respectively, confirming the role of Stab2 in regulating hepatic Malat1 ASO uptake. GWA examining ASO efficacy uncovered three loci associated with Malat1 potency: Small Subunit Processome Component (Utp11l) on chromosome 4, Rho associated coiled-coil containing protein kinase 2 (Rock2) and Aci-reductone dioxygenase (Adi1) on chromosome 12. Our results demonstrate the utility of mouse GWAS using the HMDP in detecting genes capable of impacting the uptake of ASOs, and identifies genes critical for the activity of ASOs in vivo.

Treatment modalities for cancer radiation therapy have become increasingly diversified given the growing number of facilities providing proton and carbon-ion therapy in addition to the more historically accepted photon therapy.An understanding of high-LET radiobiology is critical for optimization of charged particle radiation therapy and potential DNA damage response.In this review, we present a comprehensive summary and comparison of these types of therapy monitored primarily by using DNA damage biomarkers.We focus on their relative profiles of dose distribution and mechanisms of action from the level of nucleic acid to tumor cell death.

We hypothesize that DNA damage induced by high local energy deposition, occurring when cells are traversed by high-LET (Linear Energy Transfer) particles, can be experimentally modeled by exposing cells to high doses of low-LET. In this work, we validate such hypothesis by characterizing and correlating the time dependence of 53BP1 radiation-induced foci (RIF) for various doses and LET across 72 primary skin fibroblast from mice. This genetically diverse population allows us to understand how genetic may modulate the dose and LET relationship. The cohort was made on average from 3 males and 3 females belonging to 15 different strains of mice with various genetic backgrounds, including the collaborative cross (CC) genetic model (10 strains) and 5 reference mice strains. Cells were exposed to two fluences of three HZE (High Atomic Energy) particles (Si 350 megaelectronvolts per nucleon, Ar 350 megaelectronvolts per nucleon and Fe 600 megaelectronvolts per nucleon) and to 0.1, 1 and 4 grays from a 160 kilovolt X-ray. Individual radiation sensitivity was investigated by high throughput measurements of DNA repair kinetics for different doses of each radiation type. The 53BP1 RIF dose response to high-LET particles showed a linear dependency that matched the expected number of tracks per cell, clearly illustrating the fact that close-by DNA double strand breaks along tracks cluster within one single RIF. By comparing the slope of the high-LET dose curve to the expected number of tracks per cell we computed the number of remaining unrepaired tracks as a function of time post-irradiation. Results show that the percentage of unrepaired track over a 48 hours follow-up is higher as the LET increases across all strains. We also observe a strong correlation between the high dose repair kinetics following exposure to 160 kilovolts X-ray and the repair kinetics of high-LET tracks, with higher correlation with higher LET. At the in-vivo level for the 10-CC strains, we observe that drops in the number of T-cells and B-cells found in the blood of mice 24 hours after exposure to 0.1 gray of 320 kilovolts X-ray correlate well with slower DNA repair kinetics in skin cells exposed to X-ray. Overall, our results suggest that repair kinetics found in skin is a surrogate marker for in-vivo radiation sensitivity in other tissue, such as blood cells, and that such response is modulated by genetic variability.

This study considers the problem of describing and predicting cleft formation during the early stages of branching morphogenesis in mouse submandibular salivary glands (SMG) under the influence of varied concentrations of epidermal growth factors (EGF). Given a time-lapse video of a growing SMG, first we build a descriptive model that captures the underlying biological process and quantifies the ground truth. Tissue-scale (global) and morphological features related to regions of interest (local features) are used to characterize the biological ground truth. Second, we devise a predictive growth model that simulates EGF-modulated branching morphogenesis using a dynamic graph algorithm, which is driven by biological parameters such as EGF concentration, mitosis rate, and cleft progression rate. Given the initial configuration of the SMG, the evolution of the dynamic graph predicts the cleft formation, while maintaining the local structural characteristics of the SMG. We determined that higher EGF concentrations cause the formation of higher number of buds and comparatively shallow cleft depths. Third, we compared the prediction accuracy of our model to the Glazier-Graner-Hogeweg (GGH) model, an on-lattice Monte-Carlo simulation model, under a specific energy function parameter set that allows new rounds of de novo cleft formation. The results demonstrate that the dynamic graph model yields comparable simulations of gland growth to that of the GGH model with a significantly lower computational complexity. Fourth, we enhanced this model to predict the SMG morphology for an EGF concentration without the assistance of a ground truth time-lapse biological video data; this is a substantial benefit of our model over other similar models that are guided and terminated by information regarding the final SMG morphology. Hence, our model is suitable for testing the impact of different biological parameters involved with the process of branching morphogenesis in silico, while reducing the requirement of in vivo experiments.

Coordinated actin microfilament and microtubule dynamics is required for salivary gland development, although the mechanisms by which they contribute to branching morphogenesis are not defined. Because LIM kinase (LIMK) regulates both actin and microtubule organization, we investigated the role of LIMK signaling in mouse embryonic submandibular salivary glands using ex vivo organ cultures. Both LIMK 1 and 2 were necessary for branching morphogenesis and functioned to promote epithelial early- and late-stage cleft progression through regulation of both microfilaments and microtubules. LIMK-dependent regulation of these cytoskeletal systems was required to control focal adhesion protein–dependent fibronectin assembly and integrin β1 activation, involving the LIMK effectors cofilin and TPPP/p25, for assembly of the actin- and tubulin-based cytoskeletal systems, respectively. We demonstrate that LIMK regulates the early stages of cleft formation—cleft initiation, stabilization, and progression—via establishment of actin stability. Further, we reveal a novel role for the microtubule assembly factor p25 in regulating stabilization and elongation of late-stage progressing clefts. This study demonstrates the existence of multiple actin- and microtubule-dependent stabilization steps that are controlled by LIMK and are required in cleft progression during branching morphogenesis.

Normal organ development, function, and repair are coordinated by interactions between the epithelium and the surrounding stromal cell populations. Cellular function and homeostasis are controlled by an array of chemical and physical cues originating from the cells themselves and from the surrounding extracellular matrix (ECM). Both the endogenous cell population and ECM modulate and rely on the maintenance of basal level of tension within the tissue as a cue for growth and differentiation [1]. Furthermore, the loss of this tensional homeostasis is synonymous with many pathological conditions including; cancer, wound healing, and degenerative diseases [2].

Branching morphogenesis is a complex developmental process requiring coordination of cell proliferation, cell migration, actomyosin‐mediated cell contractility, and adhesion formation, which is not understood. We utilized computer modeling to simulate cleft progression in the primary bud during salivary gland development. Using computer simulations where individual or multiple parameters were varied in specific subsets of cells and each repeated 100 times, we examined the function of multiple cellular processes in cleft formation. We developed quantitative cleft detection algorithms to report on cleft quality and used confocal imaging of salivary gland organ explant cultures to confirm predictions of the simulations, when possible. Manipulations of the percentage of cells undergoing mitosis and the location of those cells are not currently possible to perform experimentally. Of multiple conditions, computer simulations predicted that a low proliferation rate and equal division of proliferative events between the outer and inner cells in the epithelium produces ideal clefts more than other conditions. This work demonstrates that computer modeling can be used to make predictions regarding cellular mechanisms of branching morphogenesis when experimental manipulation is not technically feasible. Supported by NIH R01DE019244 (M.L. and B.Y.)

Branching morphogenesis occurs during the development of many organs, and the embryonic mouse submandibular gland (SMG) is a classical model for the study of branching morphogenesis. In the developing SMG, this process involves iterative steps of epithelial bud and duct formation, to ultimately give rise to a complex branched network of acini and ducts, which serve to produce and modify/transport the saliva, respectively, into the oral cavity(1-3). The epithelial-associated basement membrane and aspects of the mesenchymal compartment, including the mesenchyme cells, growth factors and the extracellular matrix, produced by these cells, are critical to the branching mechanism, although how the cellular and molecular events are coordinated remains poorly understood (4). The study of the molecular mechanisms driving epithelial morphogenesis advances our understanding of developmental mechanisms and provides insight into possible regenerative medicine approaches. Such studies have been hampered due to the lack of effective methods for genetic manipulation of the salivary epithelium. Currently, adenoviral transduction represents the most effective method for targeting epithelial cells in adult glands in vivo(5). However, in embryonic explants, dense mesenchyme and the basement membrane surrounding the epithelial cells impedes viral access to the epithelial cells. If the mesenchyme is removed, the epithelium can be transfected using adenoviruses, and epithelial rudiments can resume branching morphogenesis in the presence of Matrigel or laminin-111(6,7). Mesenchyme-free epithelial rudiment growth also requires additional supplementation with soluble growth factors and does not fully recapitulate branching morphogenesis as it occurs in intact glands(8). Here we describe a technique which facilitates adenoviral transduction of epithelial cells and culture of the transfected epithelium with associated mesenchyme. Following microdissection of the embryonic SMGs, removal of the mesenchyme, and viral infection of the epithelium with a GFP-containing adenovirus, we show that the epithelium spontaneously recombines with uninfected mesenchyme, recapitulating intact SMG glandular structure and branching morphogenesis. The genetically modified epithelial cell population can be easily monitored using standard fluorescence microscopy methods, if fluorescently-tagged adenoviral constructs are used. The tissue recombination method described here is currently the most effective and accessible method for transfection of epithelial cells with a wild-type or mutant vector within a complex 3D tissue construct that does not require generation of transgenic animals.

Cleft formation during submandibular salivary gland branching morphogenesis is the critical step initiating the growth and development of the complex adult organ. Previous experimental studies indicated requirements for several epithelial cellular processes, such as proliferation, migration, cell-cell adhesion, cell-extracellular matrix (matrix) adhesion, and cellular contraction in cleft formation; however, the relative contribution of each of these processes is not fully understood since it is not possible to experimentally manipulate each factor independently. We present here a comprehensive analysis of several cellular parameters regulating cleft progression during branching morphogenesis in the epithelial tissue of an early embryonic salivary gland at a local scale using an on lattice Monte-Carlo simulation model, the Glazier-Graner-Hogeweg model. We utilized measurements from time-lapse images of mouse submandibular gland organ explants to construct a temporally and spatially relevant cell-based 2D model. Our model simulates the effect of cellular proliferation, actomyosin contractility, cell-cell and cell-matrix adhesions on cleft progression, and it was used to test specific hypotheses regarding the function of these parameters in branching morphogenesis. We use innovative features capturing several aspects of cleft morphology and quantitatively analyze clefts formed during functional modification of the cellular parameters. Our simulations predict that a low epithelial mitosis rate and moderate level of actomyosin contractility in the cleft cells promote cleft progression. Raising or lowering levels of contractility and mitosis rate resulted in non-progressive clefts. We also show that lowered cell-cell adhesion in the cleft region and increased cleft cell-matrix adhesions are required for cleft progression. Using a classifier-based analysis, the relative importance of these four contributing cellular factors for effective cleft progression was determined as follows: cleft cell contractility, cleft region cell-cell adhesion strength, epithelial cell mitosis rate, and cell-matrix adhesion strength.

In this paper, we introduce a biologically motivated dynamic graph-based growth model to describe and predict the stages of cleft formation during the process of branching morphogenesis in the submandibular mouse gland (SMG) from 3 hrs after embryonic day E12 to 8 hrs after embryonic day E12, which can be considered as E12.5. Branching morphogenesis is the process by which many mammalian exocrine and endocrine glands undergo significant morphological transformations, from a primary bud to an adult organ. Although many studies have investigated the cellular and molecular mechanisms driving branching morphogenesis, it is not clear how the shape changes that are inherent to establishing organ structure are produced. Using morphological features extracted from sequential images of SMG organ cultures we were able to develop a dynamic graph-based predictive model that is able to mimic the process of cleft formation and predict the final state. In addition, we compare our model to a state-of-the-art Glazier-Graner-Hogeweg (GGH) simulative tool, and demonstrate that the dynamic graph-based predictive model has comparable accuracy in modeling growth of clefts across SMG developmental stages, as well as faster convergence to the target SMG morphology.

Pattern formation in developing tissues involves dynamic spatio-temporal changes in cellular organization and subsequent evolution of functional adult structures. Branching morphogenesis is a developmental mechanism by which patterns are generated in many developing organs, which is controlled by underlying molecular pathways. Understanding the relationship between molecular signaling, cellular behavior and resulting morphological change requires quantification and categorization of the cellular behavior. In this study, tissue-level and cellular changes in developing salivary gland in response to disruption of ROCK-mediated signaling by are modeled by building cell-graphs to compute mathematical features capturing structural properties at multiple scales. These features were used to generate multiscale cell-graph signatures of untreated and ROCK signaling disrupted salivary gland organ explants. From confocal images of mouse submandibular salivary gland organ explants in which epithelial and mesenchymal nuclei were marked, a multiscale feature set capturing global structural properties, local structural properties, spectral, and morphological properties of the tissues was derived. Six feature selection algorithms and multiway modeling of the data was performed to identify distinct subsets of cell graph features that can uniquely classify and differentiate between different cell populations. Multiscale cell-graph analysis was most effective in classification of the tissue state. Cellular and tissue organization, as defined by a multiscale subset of cell-graph features, are both quantitatively distinct in epithelial and mesenchymal cell types both in the presence and absence of ROCK inhibitors. Whereas tensor analysis demonstrate that epithelial tissue was affected the most by inhibition of ROCK signaling, significant multiscale changes in mesenchymal tissue organization were identified with this analysis that were not identified in previous biological studies. We here show how to define and calculate a multiscale feature set as an effective computational approach to identify and quantify changes at multiple biological scales and to distinguish between different states in developing tissues.

Branching morphogenesis is a developmental process shared by many organs, including the submandibular salivary gland. During morphogenesis, cells within the gland undergo rearrangements to cause changes in the overall tissue morphology. This work presents a methodology based on cell-graphs to quantify these changes in cellular arrangements. Multiple confocal images of developing salivary gland organ cultures are captured. These cultures are immunostained with a nuclear marker and an epithelial marker to identify epithelial cells as separate from mesenchymal cells. Confocal images are stitched and segmented to identify epithelial and mesenchymal nuclei. Cell-graphs are constructed to model the structural organization of epithelial and mesenchymal cells. Cell-graph metrics are calculated to extract mathematical features that discriminate epithelial vs mesenchymal cells organizations and also distinguish between glands treated with pharmacological inhibitors vs vehicle control. The results indicate that cell-graph features can be used to both describe and predict the developing salivary gland to provide insights into cellular and physical processes driving morphogenesis.