LIM Kinase Regulation of Cytoskeletal Dynamics is Required for Salivary Gland Branching Morphogenesis

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

Systemic Response to Microgravity: Utilizing GeneLab Datasets to Identify Molecular Targets for Future Hypotheses-Driven Spaceflight Studies

Biological risks associated with microgravity are a major concern for long-term space travel. Although determination of risk has been a focus for NASA research, data examining systemic (i.e., multi- or pan-tissue) responses to space flight are sparse. To perform our analysis, we utilized the NASA GeneLab database which is a publicly available repository containing a wide array of omics results from experiments conducted with: i) with different flight conditions (space shuttle (STS) missions vs. International Space Station (ISS); ii) a variety of tissues; and 3) assays that measure epigenetic, transcriptional, and protein expression changes. Meta-analysis of the transcriptomic data from 7 different murine and rat data sets, examining tissues such as liver, kidney, adrenal gland, thymus, mammary gland, skin, and skeletal muscle (soleus, extensor digitorum longus, tibialis anterior, quadriceps, and gastrocnemius) revealed for the first time, the existence of potential master regulators coordinating systemic responses to microgravity in rodents. We identified p53, TGF(beta)1 and immune related pathways as the highly prevalent pan-tissue signaling pathways that are affected by microgravity. Some variability in the degree of change in their expression across species, strain and time of flight was also observed. Interestingly, while certain skeletal muscle (gastrocnemius and soleus) exhibited an overall down-regulation of these genes, some other muscle types such as the extensor digitorum longus, tibialis anterior and quadriceps, showed an up-regulated expression, indicative of potential compensatory mechanisms to prevent microgravity-induced atrophy. Key genes isolated by unbiased systems analyses displayed a major overlap between tissue types and flight conditions and established TGF(beta)1 to be the most connected gene across all data sets. Finally, a set of microgravity responsive miRNA signature was identified and based on their predicted functional state and subsequent impact on health, a theoretical health risk score was calculated. The genes and miRNAs identified from our analyses can be targeted for future research involving efficient countermeasure design. Our study thus exemplifies the utility of GeneLab data repository to aid in the process of performing novel hypothesis based spaceflight research aimed at elucidating the global impact of environmental stressors at multiple biological scales

Cell-Based Multi-Parametric Model of Cleft Progression during Submandibular Salivary Gland Branching Morphogenesis

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

GeneLab: Omics Database for Spaceflight Experiments

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)

Systemic Microgravity Response: Utilizing GeneLab to Develop Hypotheses for Spaceflight Risks

Biological risks associated with microgravity are a major concern for long-term space travel. Although determination of risk has been a focus for NASA research, data examining systemic (i.e., multi- or pan-tissue) responses to space flight are sparse. To perform our analysis, we utilized the NASA GeneLab database which is a publicly available repository containing a wide array of omics results from experiments conducted with: i) with different flight conditions (space shuttle (STS) missions vs. International Space Station (ISS); ii) a variety of tissues; and 3) assays that measure epigenetic, transcriptional, and protein expression changes. Meta-analysis of the transcriptomic data from 7 different murine and rat data sets, examining tissues such as liver, kidney, adrenal gland, thymus, mammary gland, skin, and skeletal muscle (soleus, extensor digitorum longus, tibialis anterior, quadriceps, and gastrocnemius) revealed for the first time, the existence of potential master regulators coordinating systemic responses to microgravity in rodents. We identified p53, TGF1 and immune related pathways as the highly prevalent pan-tissue signaling pathways that are affected by microgravity. Some variability in the degree of change in their expression across species, strain and time of flight was also observed. Interestingly, while certain skeletal muscle (gastrocnemius and soleus) exhibited an overall down-regulation of these genes, some other muscle types such as the extensor digitorum longus, tibialis anterior and quadriceps, showed an up-regulated expression, indicative of potential compensatory mechanisms to prevent microgravity-induced atrophy. Key genes isolated by unbiased systems analyses displayed a major overlap between tissue types and flight conditions and established TGF1 to be the most connected gene across all data sets. Finally, a set of microgravity responsive miRNA signature was identified and based on their predicted functional state and subsequent impact on health, a theoretical health risk score was calculated. The genes and miRNAs identified from our analyses can be targeted for future research involving efficient countermeasure design. Our study thus exemplifies the utility of GeneLab data repository to aid in the process of performing novel hypothesis based spaceflight research aimed at elucidating the global impact of environmental stressors at multiple biological scales

Genetic Correlation with the DNA Repair Assay in Mice Exposed to High-LET

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

Variability in Galactic Cosmic Radiation- Induced DNA Damage Response in Inbred Mice Is Modulated by Genetics

In radiation biology, the ability to predict cancer risk associated with exposure to low doses of high-LET (Linear Energy Transfer) ionizing radiation remains a challenge. Epidemiological methods lack the sensitivity and power to provide detailed risk estimates for cancer and ignore individual sensitivity. We have hypothesized that DNA repair capacity is the primary factor differentiating peoples radiation sensitivity. We previously showed in immortalized human cell lines that characterizing the dose and time dependence of p53-binding protein 1 (53BP1) foci formation in the nucleus following X-rays exposure is sufficient to predict DNA repair response to any other LET in the same cell line. We now tested this hypothesis across a population of mice with different genetic background. Fibroblast cells were extracted and cultivated from 76 individual mice from 15 different strains and exposed to HZE (high (H) atomic number (Z) and energy (E) galactic cosmic ray particles) particles and X-rays. Individual radiation sensitivities were investigated by high throughput measurement of DNA repair kinetics that evaluated 53bp1 foci numbers as a surrogate for DNA double-strand breaks at various times post-irradiation. Instead of just counting foci which can be hard to distinguish for high-LET or high doses, we also took into account the track structure of high-LET particles to compute the remaining number of unrepaired tracks as a function of time post-irradiation. As expected, the percentage of unrepaired track over a 48 hours follow-up period increased with LET. In addition, repair rate was modulated by genetics, with animals from the same strain showing small variance while large rate differences were observed between strains. Radiation strain sensitivity ranking was estimated based on repair rates from exposure to each LET evaluated in this work, and ranking for high-LET correlated better with ranking from high dose of X-ray, not low dose. At the in-vivo level, drops in T-cells and B-cells number measured 24 hours after 0.1 Gy (Gray) X-ray exposure, correlated with slower DNA repair kinetic in fibroblast cells of the same strains of mice. At the genomic level, mouse genome wide association (GWA) analysis identified seven significant genetic loci on chromosomes 2, 3, 7, 10, 11, 13 and 19 with different significance depending on the LET. Interestingly, for the two highest LET, a common locus on Chromosome 10 was identified with high enrichment for DNA repair associated genes.Overall, this work suggests that repair kinetics of primary skin fibroblasts is a good surrogate marker for in-vivo radiation sensitivities in other tissues and that this response is modulated by genetics. Our study also confirms that DNA repair kinetics following high doses of X-ray can be used to predict radiation sensitivity to high-LET

Mouse genome-wide association studies and systems genetics uncover the genetic architecture associated with hepatic pharmacokinetic and pharmacodynamic properties of a constrained ethyl antisense oligonucleotide targeting Malat1

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 ciseQTL 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

Multiscale Feature Analysis of Salivary Gland Branching Morphogenesis

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

Genetic Modification and Recombination of Salivary Gland Organ Cultures

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 cavity1-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 vivo5. 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-1116,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 glands8. 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