Comprehensive Multi-omics Analysis Reveals Mitochondrial Stress as a Central Biological Hub for Spaceflight Impact.

Spaceflight is known to impose changes on human physiology with unknown molecular etiologies. To reveal these causes, we used a multi-omics, systems biology analytical approach using biomedical profiles from fifty-nine astronauts and data from NASA's GeneLab derived from hundreds of samples flown in space to determine transcriptomic, proteomic, metabolomic, and epigenetic responses to spaceflight. Overall pathway analyses on the multi-omics datasets showed significant enrichment for mitochondrial processes, as well as innate immunity, chronic inflammation, cell cycle, circadian rhythm, and olfactory functions. Importantly, NASA's Twin Study provided a platform to confirm several of our principal findings. Evidence of altered mitochondrial function and DNA damage was also found in the urine and blood metabolic data compiled from the astronaut cohort and NASA Twin Study data, indicating mitochondrial stress as a consistent phenotype of spaceflight.

Genetic association and characterization of FSTL5 in isolated clubfoot.

Talipes equinovarus (clubfoot, TEV) is a congenital rotational foot deformity occurring in 1 per 1000 births with increased prevalence in males compared with females. The genetic etiology of isolated clubfoot (iTEV) remains unclear. Using a genome-wide association study, we identified a locus within FSTL5, encoding follistatin-like 5, significantly associated with iTEV. FSTL5 is an uncharacterized gene whose potential role in embryonic and postnatal development was previously unstudied. Utilizing multiple model systems, we found that Fstl5 was expressed during later stages of embryonic hindlimb development, and, in mice, expression was restricted to the condensing cartilage anlage destined to form the limb skeleton. In the postnatal growth plate, Fstl5 was specifically expressed in prehypertrophic chondrocytes. As Fstl5 knockout rats displayed no gross malformations, we engineered a conditional transgenic mouse line (Fstl5LSL) to overexpress Fstl5 in skeletal osteochondroprogenitors. We observed that hindlimbs were slightly shorter and that bone mineral density was reduced in adult male, but not female, Prrx1-cre;Fstl5LSL mice compared with control. No overt clubfoot-like deformity was observed in Prrx1-cre;Fstl5LSL mice, suggesting FSTL5 may function in other cell types to contribute to iTEV pathogenesis. Interrogating published mouse embryonic single-cell expression data showed that Fstl5 was expressed in cell lineage subclusters whose transcriptomes were associated with neural system development. Moreover, our results suggest that lineage-specific expression of the Fstl genes correlates with their divergent roles as modulators of transforming growth factor beta and bone morphogenetic protein signaling. Results from this study associate FSTL5 with iTEV and suggest a potential sexually dimorphic role for Fstl5 in vivo.

Diagnostic yield and clinical impact of exome sequencing in early-onset scoliosis (EOS).

BACKGROUND: Early-onset scoliosis (EOS), defined by an onset age of scoliosis less than 10 years, conveys significant health risk to affected children. Identification of the molecular aetiology underlying patients with EOS could provide valuable information for both clinical management and prenatal screening. METHODS: In this study, we consecutively recruited a cohort of 447 Chinese patients with operative EOS. We performed exome sequencing (ES) screening on these individuals and their available family members (totaling 670 subjects). Another cohort of 13 patients with idiopathic early-onset scoliosis (IEOS) from the USA who underwent ES was also recruited. RESULTS: After ES data processing and variant interpretation, we detected molecular diagnostic variants in 92 out of 447 (20.6%) Chinese patients with EOS, including 8 patients with molecular confirmation of their clinical diagnosis and 84 patients with molecular diagnoses of previously unrecognised diseases underlying scoliosis. One out of 13 patients with IEOS from the US cohort was molecularly diagnosed. The age at presentation, the number of organ systems involved and the Cobb angle were the three top features predictive of a molecular diagnosis. CONCLUSION: ES enabled the molecular diagnosis/classification of patients with EOS. Specific clinical features/feature pairs are able to indicate the likelihood of gaining a molecular diagnosis through ES.

Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci.

Adolescent idiopathic scoliosis (AIS) is the most common musculoskeletal disorder of childhood development. The genetic architecture of AIS is complex, and the great majority of risk factors are undiscovered. To identify new AIS susceptibility loci, we conducted the first genome-wide meta-analysis of AIS genome-wide association studies, including 7956 cases and 88 459 controls from 3 ancestral groups. Three novel loci that surpassed genome-wide significance were uncovered in intragenic regions of the CDH13 (P-value_rs4513093 = 1.7E-15), ABO (P-value_ rs687621 = 7.3E-10) and SOX6 (P-value_rs1455114 = 2.98E-08) genes. Restricting the analysis to females improved the associations at multiple loci, most notably with variants within CDH13 despite the reduction in sample size. Genome-wide gene-functional enrichment analysis identified significant perturbation of pathways involving cartilage and connective tissue development. Expression of both SOX6 and CDH13 was detected in cartilage chondrocytes and chromatin immunoprecipitation sequencing experiments in that tissue revealed multiple HeK27ac-positive peaks overlapping associated loci. Our results further define the genetic architecture of AIS and highlight the importance of vertebral cartilage development in its pathogenesis.

Distal chromosome 16p11.2 duplications containing SH2B1 in patients with scoliosis.

INTRODUCTION: Adolescent idiopathic scoliosis (AIS) is a common musculoskeletal disorder with strong evidence for a genetic contribution. CNVs play an important role in congenital scoliosis, but their role in idiopathic scoliosis has been largely unexplored. METHODS: Exome sequence data from 1197 AIS cases and 1664 in-house controls was analysed using coverage data to identify rare CNVs. CNV calls were filtered to include only highly confident CNVs with >10 average reads per region and mean log-ratio of coverage consistent with single-copy duplication or deletion. The frequency of 55 common recurrent CNVs was determined and correlated with clinical characteristics. RESULTS: Distal chromosome 16p11.2 microduplications containing the gene SH2B1 were found in 0.7% of AIS cases (8/1197). We replicated this finding in two additional AIS cohorts (8/1097 and 2/433), resulting in 0.7% (18/2727) of all AIS cases harbouring a chromosome 16p11.2 microduplication, compared with 0.06% of local controls (1/1664) and 0.04% of published controls (8/19584) (p=2.28×10-11, OR=16.15). Furthermore, examination of electronic health records of 92 455 patients from the Geisinger health system showed scoliosis in 30% (20/66) patients with chromosome 16p11.2 microduplications containing SH2B1 compared with 7.6% (10/132) of controls (p=5.6×10-4, OR=3.9). CONCLUSIONS: Recurrent distal chromosome 16p11.2 duplications explain nearly 1% of AIS. Distal chromosome 16p11.2 duplications may contribute to scoliosis pathogenesis by directly impairing growth or by altering expression of nearby genes, such as TBX6. Individuals with distal chromosome 16p11.2 microduplications should be screened for scoliosis to facilitate early treatment.

Bone Marrow Endothelial Cells Influence Function and Phenotype of Hematopoietic Stem and Progenitor Cells after Mixed Neutron/Gamma Radiation.

The bone marrow (BM) microenvironment plays a crucial role in the maintenance and regeneration of hematopoietic stem (HSC) and progenitor cells (HSPC). In particular, the vascular niche is responsible for regulating HSC maintenance, differentiation, and migration of cells in and out of the BM. Damage to this niche upon exposure to ionizing radiation, whether accidental or as a result of therapy, can contribute to delays in HSC recovery and/or function. The ability of BM derived-endothelial cells (BMEC) to alter and/or protect HSPC after exposure to ionizing radiation was investigated. Our data show that exposure of BMEC to ionizing radiation resulted in alterations in Akt signaling, increased expression of PARP-1, IL6, and MCP-1, and decreased expression of MMP1 and MMP9. In addition, global analysis of gene expression of HSC and BMEC in response to mixed neutron/gamma field (MF) radiation identified 60 genes whose expression was altered after radiation in both cell types, suggesting that a subset of genes is commonly affected by this type of radiation. Focused gene analysis by RT-PCR revealed two categories of BMEC alterations: (a) a subset of genes whose expression was altered in response to radiation, with no additional effect observed during coculture with HSPC, and (b) a subset of genes upregulated in response to radiation, and altered when cocultured with HSPC. Coculture of BMEC with CD34+ HSPC induced HSPC proliferation, and improved BM function after MF radiation. Nonirradiated HSPC exhibited reduced CD34 expression over time, but when irradiated, they maintained higher CD34 expression. Nonirradiated HSPC cocultured with nonirradiated BMEC expressed lower levels of CD34 expression compared to nonirradiated alone. These data characterize the role of each cell type in response to MF radiation and demonstrate the interdependence of each cell's response to ionizing radiation. The identified genes modulated by radiation and coculture provide guidance for future experiments to test hypotheses concerning specific factors mediating the beneficial effects of BMEC on HSPC. This information will prove useful in the search for medical countermeasures to radiation-induced hematopoietic injury.

Proteomic and phosphoproteomic characterization of cardiovascular tissues after long term exposure to simulated space radiation.

Introduction: It may take decades to develop cardiovascular dysfunction following exposure to high doses of ionizing radiation from medical therapy or from nuclear accidents. Since astronauts may be exposed continually to a complex space radiation environment unlike that experienced on Earth, it is unresolved whether there is a risk to cardiovascular health during long-term space exploration missions. Previously, we have described that mice exposed to a single dose of simplified Galactic Cosmic Ray (GCR5-ion) develop cardiovascular dysfunction by 12 months post-radiation. Methods: To investigate the biological basis of this dysfunction, here we performed a quantitative mass spectrometry-based proteomics analysis of heart tissue (proteome and phosphoproteome) and plasma (proteome only) from these mice at 8 months post-radiation. Results: Differentially expressed proteins (DEPs) for irradiated versus sham irradiated samples (fold-change ≥1.2 and an adjusted p-value of ≤0.05) were identified for each proteomics data set. For the heart proteome, there were 87 significant DEPs (11 upregulated and 76 downregulated); for the heart phosphoproteome, there were 60 significant differentially phosphorylated peptides (17 upregulated and 43 downregulated); and for the plasma proteome, there was only one upregulated protein. A Gene Set Enrichment Analysis (GSEA) technique that assesses canonical pathways from BIOCARTA, KEGG, PID, REACTOME, and WikiPathways revealed significant perturbation in pathways in each data set. For the heart proteome, 166 pathways were significantly altered (36 upregulated and 130 downregulated); for the plasma proteome, there were 73 pathways significantly altered (25 upregulated and 48 downregulated); and for the phosphoproteome, there were 223 pathways significantly affected at 0.1 adjusted p-value cutoff. Pathways related to inflammation were the most highly perturbed in the heart and plasma. In line with sustained inflammation, neutrophil extracellular traps (NETs) were demonstrated to be increased in GCR5-ion irradiated hearts at 12-month post irradiation. NETs play a fundamental role in combating bacterial pathogens, modulating inflammatory responses, inflicting damage on healthy tissues, and escalating vascular thrombosis. Discussion: These findings suggest that a single exposure to GCR5-ion results in long-lasting changes in the proteome and that these proteomic changes can potentiate acute and chronic health issues for astronauts, such as what we have previously described with late cardiac dysfunction in these mice.

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than fivefold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here, we sought to define the roles of PAX1 and newly identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); p=7.07E-11, OR = 1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wild-type. By genetic targeting we found that wild-type Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, the latter suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2 or tamoxifen treatment significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a PAX1-COL11a1-MMP3 signaling axis in spinal chondrocytes.

Adolescent idiopathic scoliosis (AIS) is a twisting deformity of the spine that occurs during periods of rapid growth in children worldwide. Children with severe cases of AIS require surgery to stop it from getting worse, presenting a significant financial burden to health systems and families. Although AIS is known to cluster in families, its genetic causes and its inheritance pattern have remained elusive. Additionally, AIS is known to be more prevalent in females, a bias that has not been explained. Advances in techniques to study the genetics underlying diseases have revealed that certain variations that increase the risk of AIS affect cartilage and connective tissue. In humans, one such variation is near a gene called Pax1, and it is female-specific. The extracellular matrix is a network of proteins and other molecules in the space between cells that help connect tissues together, and it is particularly important in cartilage and other connective tissues. One of the main components of the extracellular matrix is collagen. Yu, Kanshour, Ushiki et al. hypothesized that changes in the extracellular matrix could affect the cartilage and connective tissues of the spine, leading to AIS. To show this, the scientists screened over 100,000 individuals and found that AIS is associated with variants in two genes coding for extracellular matrix proteins. One of these variants was found in a gene called Col11a1, which codes for one of the proteins that makes up collagen. To understand the relationship between Pax1 and Col11a1, Yu, Kanshour, Ushiki et al. genetically modified mice so that they would lack the Pax1 gene. In these mice, the activation of Col11a1 was reduced in the mouse spine. They also found that the form of Col11a1 associated with AIS could not suppress the activation of a gene called Mmp3 in mouse cartilage cells as effectively as unmutated Col11a1. Going one step further, the researchers found that lowering the levels of an estrogen receptor altered the activation patterns of Pax1, Col11a1, and Mmp3 in mouse cartilage cells. These findings suggest a possible mechanism for AIS, particularly in females. The findings of Yu, Kanshour, Ushiki et al. highlight that cartilage cells in the spine are particularly relevant in AIS. The results also point to specific molecules within the extracellular matrix as important for maintaining proper alignment in the spine when children are growing rapidly. This information may guide future therapies aimed at maintaining healthy spinal cells in adolescent children, particularly girls.

Computational approaches for discovery of common immunomodulators in fungal infections: towards broad-spectrum immunotherapeutic interventions.

BACKGROUND: Fungi are the second most abundant type of human pathogens. Invasive fungal pathogens are leading causes of life-threatening infections in clinical settings. Toxicity to the host and drug-resistance are two major deleterious issues associated with existing antifungal agents. Increasing a host's tolerance and/or immunity to fungal pathogens has potential to alleviate these problems. A host's tolerance may be improved by modulating the immune system such that it responds more rapidly and robustly in all facets, ranging from the recognition of pathogens to their clearance from the host. An understanding of biological processes and genes that are perturbed during attempted fungal exposure, colonization, and/or invasion will help guide the identification of endogenous immunomodulators and/or small molecules that activate host-immune responses such as specialized adjuvants. RESULTS: In this study, we present computational techniques and approaches using publicly available transcriptional data sets, to predict immunomodulators that may act against multiple fungal pathogens. Our study analyzed data sets derived from host cells exposed to five fungal pathogens, namely, Alternaria alternata, Aspergillus fumigatus, Candida albicans, Pneumocystis jirovecii, and Stachybotrys chartarum. We observed statistically significant associations between host responses to A. fumigatus and C. albicans. Our analysis identified biological processes that were consistently perturbed by these two pathogens. These processes contained both immune response-inducing genes such as MALT1, SERPINE1, ICAM1, and IL8, and immune response-repressing genes such as DUSP8, DUSP6, and SPRED2. We hypothesize that these genes belong to a pool of common immunomodulators that can potentially be activated or suppressed (agonized or antagonized) in order to render the host more tolerant to infections caused by A. fumigatus and C. albicans. CONCLUSIONS: Our computational approaches and methodologies described here can now be applied to newly generated or expanded data sets for further elucidation of additional drug targets. Moreover, identified immunomodulators may be used to generate experimentally testable hypotheses that could help in the discovery of broad-spectrum immunotherapeutic interventions. All of our results are available at the following supplementary website: http://bioinformatics.cs.vt.edu/~murali/supplements/2013-kidane-bmc.

Staphylococcus aureus isolates from children with clinically differentiated osteomyelitis exhibit distinct transcriptomic signatures.

There is substantial genomic heterogeneity among Staphylococcus aureus isolates of children with acute hematogenous osteomyelitis (AHO) but transcriptional behavior of clinically differentiated strains has not been previously described. This study evaluates transcriptional activity of S. aureus isolates of children with AHO that may regulate metabolism, biosynthesis, or virulence during bacterial growth and pathogenesis. In vitro growth kinetics were compared between three S. aureus clinical isolates from children with AHO who had mild, moderate, and severe illness. Total RNA sequencing was performed for each isolate at six separate time points throughout the logarithmic phase of growth. The NASA RNA-Sequencing Consensus Pipeline was used to identify differentially expressed genes allowing for 54 comparisons between the three isolates during growth. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment pathways were used to evaluate transcriptional variation in metabolism, biosynthesis pathways and virulence potential of the isolates. The S. aureus isolates demonstrated differing growth kinetics under standardized conditions with the mild isolate having higher optical densities with earlier and higher peak rates of growth than that of the other isolates (p<0.001). Enrichment pathway analysis established distinct transcriptional signatures according to both sampling time and clinical severity. Moderate and severe isolates demonstrated pathways of bacterial invasion, S. aureus infection, quorum sensing and two component systems. In comparison, the mild strain favored biosynthesis and metabolism. These findings suggest that transcriptional regulation during the growth of S. aureus may impact the pathogenetic mechanisms involved in the progression of severity of illness in childhood osteomyelitis. The clinical isolates studied demonstrated a tradeoff between growth and virulence. Further investigation is needed to evaluate these transcriptional pathways in an animal model or during active clinical infections of children with AHO.

Molecular Dissection of Somatic Skeletal Disease in Neurofibromatosis Type 1.

Neurofibromatosis type 1 (NF1) is a tumor predisposition syndrome caused by heterozygous NF1 gene mutations. Patients with NF1 present with pleiotropic somatic secondary manifestations, including development of bone pseudarthrosis after fracture. Somatic NF1 gene mutations were reproducibly identified in patient-derived pseudarthrosis specimens, suggesting a local mosaic cell population including somatic pathologic cells. The somatic cellular pathogenesis of NF1 pseudarthroses remains unclear, though defects in osteogenesis have been posited. Here, we applied time-series single-cell RNA-sequencing (scRNA-seq) to patient-matched control and pseudarthrosis-derived primary bone stromal cells (BSCs). We show that osteogenic specification to an osteoblast progenitor cell population was evident for control bone-derived cells and haploinsufficient pseudarthrosis-derived cells. Similar results were observed for somatic patient fracture-derived NF1-/- cells; however, expression of genetic pathways associated with skeletal mineralization were significantly reduced in NF1-/- cells compared with fracture-derived NF1+/- cells. In mice, we show that Nf1 expressed in bone marrow osteoprogenitors is required for the maintenance of the adult skeleton. Results from our study implicate impaired Clec11a-Itga11-Wnt signaling in the pathogenesis of NF1-associated skeletal disease. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Association of genetic variation in COL11A1 with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is a common and progressive spinal deformity in children that exhibits striking sexual dimorphism, with girls at more than five-fold greater risk of severe disease compared to boys. Despite its medical impact, the molecular mechanisms that drive AIS are largely unknown. We previously defined a female-specific AIS genetic risk locus in an enhancer near the PAX1 gene. Here we sought to define the roles of PAX1 and newly-identified AIS-associated genes in the developmental mechanism of AIS. In a genetic study of 10,519 individuals with AIS and 93,238 unaffected controls, significant association was identified with a variant in COL11A1 encoding collagen (α1) XI (rs3753841; NM_080629.2_c.4004C>T; p.(Pro1335Leu); P=7.07e-11, OR=1.118). Using CRISPR mutagenesis we generated Pax1 knockout mice (Pax1-/-). In postnatal spines we found that PAX1 and collagen (α1) XI protein both localize within the intervertebral disc (IVD)-vertebral junction region encompassing the growth plate, with less collagen (α1) XI detected in Pax1-/- spines compared to wildtype. By genetic targeting we found that wildtype Col11a1 expression in costal chondrocytes suppresses expression of Pax1 and of Mmp3, encoding the matrix metalloproteinase 3 enzyme implicated in matrix remodeling. However, this suppression was abrogated in the presence of the AIS-associated COL11A1P1335L mutant. Further, we found that either knockdown of the estrogen receptor gene Esr2, or tamoxifen treatment, significantly altered Col11a1 and Mmp3 expression in chondrocytes. We propose a new molecular model of AIS pathogenesis wherein genetic variation and estrogen signaling increase disease susceptibility by altering a Pax1-Col11a1-Mmp3 signaling axis in spinal chondrocytes.

Core mitochondrial genes are down-regulated during SARS-CoV-2 infection of rodent and human hosts.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral proteins bind to host mitochondrial proteins, likely inhibiting oxidative phosphorylation (OXPHOS) and stimulating glycolysis. We analyzed mitochondrial gene expression in nasopharyngeal and autopsy tissues from patients with coronavirus disease 2019 (COVID-19). In nasopharyngeal samples with declining viral titers, the virus blocked the transcription of a subset of nuclear DNA (nDNA)-encoded mitochondrial OXPHOS genes, induced the expression of microRNA 2392, activated HIF-1α to induce glycolysis, and activated host immune defenses including the integrated stress response. In autopsy tissues from patients with COVID-19, SARS-CoV-2 was no longer present, and mitochondrial gene transcription had recovered in the lungs. However, nDNA mitochondrial gene expression remained suppressed in autopsy tissue from the heart and, to a lesser extent, kidney, and liver, whereas mitochondrial DNA transcription was induced and host-immune defense pathways were activated. During early SARS-CoV-2 infection of hamsters with peak lung viral load, mitochondrial gene expression in the lung was minimally perturbed but was down-regulated in the cerebellum and up-regulated in the striatum even though no SARS-CoV-2 was detected in the brain. During the mid-phase SARS-CoV-2 infection of mice, mitochondrial gene expression was starting to recover in mouse lungs. These data suggest that when the viral titer first peaks, there is a systemic host response followed by viral suppression of mitochondrial gene transcription and induction of glycolysis leading to the deployment of antiviral immune defenses. Even when the virus was cleared and lung mitochondrial function had recovered, mitochondrial function in the heart, kidney, liver, and lymph nodes remained impaired, potentially leading to severe COVID-19 pathology.

TARGETED DOWN REGULATION OF CORE MITOCHONDRIAL GENES DURING SARS-COV-2 INFECTION.

Defects in mitochondrial oxidative phosphorylation (OXPHOS) have been reported in COVID-19 patients, but the timing and organs affected vary among reports. Here, we reveal the dynamics of COVID-19 through transcription profiles in nasopharyngeal and autopsy samples from patients and infected rodent models. While mitochondrial bioenergetics is repressed in the viral nasopharyngeal portal of entry, it is up regulated in autopsy lung tissues from deceased patients. In most disease stages and organs, discrete OXPHOS functions are blocked by the virus, and this is countered by the host broadly up regulating unblocked OXPHOS functions. No such rebound is seen in autopsy heart, results in severe repression of genes across all OXPHOS modules. Hence, targeted enhancement of mitochondrial gene expression may mitigate the pathogenesis of COVID-19.

NASA GeneLab RNA-seq consensus pipeline: standardized processing of short-read RNA-seq data.

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.

The cartilage matrisome in adolescent idiopathic scoliosis.

The human spinal column is a dynamic, segmented, bony, and cartilaginous structure that protects the neurologic system and simultaneously provides balance and flexibility. Children with developmental disorders that affect the patterning or shape of the spine can be at risk of neurologic and other physiologic dysfunctions. The most common developmental disorder of the spine is scoliosis, a lateral deformity in the shape of the spinal column. Scoliosis may be part of the clinical spectrum that is observed in many developmental disorders, but typically presents as an isolated symptom in otherwise healthy adolescent children. Adolescent idiopathic scoliosis (AIS) has defied understanding in part due to its genetic complexity. Breakthroughs have come from recent genome-wide association studies (GWAS) and next generation sequencing (NGS) of human AIS cohorts, as well as investigations of animal models. These studies have identified genetic associations with determinants of cartilage biogenesis and development of the intervertebral disc (IVD). Current evidence suggests that a fraction of AIS cases may arise from variation in factors involved in the structural integrity and homeostasis of the cartilaginous extracellular matrix (ECM). Here, we review the development of the spine and spinal cartilages, the composition of the cartilage ECM, the so-called "matrisome" and its functions, and the players involved in the genetic architecture of AIS. We also propose a molecular model by which the cartilage matrisome of the IVD contributes to AIS susceptibility.

Modulation of nucleosomal DNA accessibility via charge-altering post-translational modifications in histone core.

BACKGROUND: Controlled modulation of nucleosomal DNA accessibility via post-translational modifications (PTM) is a critical component to many cellular functions. Charge-altering PTMs in the globular histone core-including acetylation, phosphorylation, crotonylation, propionylation, butyrylation, formylation, and citrullination-can alter the strong electrostatic interactions between the oppositely charged nucleosomal DNA and the histone proteins and thus modulate accessibility of the nucleosomal DNA, affecting processes that depend on access to the genetic information, such as transcription. However, direct experimental investigation of the effects of these PTMs is very difficult. Theoretical models can rationalize existing observations, suggest working hypotheses for future experiments, and provide a unifying framework for connecting PTMs with the observed effects. RESULTS: A physics-based framework is proposed that predicts the effect of charge-altering PTMs in the histone core, quantitatively for several types of lysine charge-neutralizing PTMs including acetylation, and qualitatively for all phosphorylations, on the nucleosome stability and subsequent changes in DNA accessibility, making a connection to resulting biological phenotypes. The framework takes into account multiple partially assembled states of the nucleosome at the atomic resolution. The framework is validated against experimentally known nucleosome stability changes due to the acetylation of specific lysines, and their effect on transcription. The predicted effect of charge-altering PTMs on DNA accessibility can vary dramatically, from virtually none to a strong, region-dependent increase in accessibility of the nucleosomal DNA; in some cases, e.g., H4K44, H2AK75, and H2BK57, the effect is significantly stronger than that of the extensively studied acetylation sites such H3K56, H3K115 or H3K122. Proximity to the DNA is suggestive of the strength of the PTM effect, but there are many exceptions. For the vast majority of charge-altering PTMs, the predicted increase in the DNA accessibility should be large enough to result in a measurable modulation of transcription. However, a few possible PTMs, such as acetylation of H4K77, counterintuitively decrease the DNA accessibility, suggestive of the repressed chromatin. A structural explanation for the phenomenon is provided. For the majority of charge-altering PTMs, the effect on DNA accessibility is simply additive (noncooperative), but there are exceptions, e.g., simultaneous acetylation of H4K79 and H3K122, where the combined effect is amplified. The amplification is a direct consequence of the nucleosome-DNA complex having more than two structural states. The effect of individual PTMs is classified based on changes in the accessibility of various regions throughout the nucleosomal DNA. The PTM's resulting imprint on the DNA accessibility, "PTMprint," is used to predict effects of many yet unexplored PTMs. For example, acetylation of H4K44 yields a PTMprint similar to the PTMprint of H3K56, and thus acetylation of H4K44 is predicted to lead to a wide range of strong biological effects. CONCLUSION: Charge-altering post-translational modifications in the relatively unexplored globular histone core may provide a precision mechanism for controlling accessibility to the nucleosomal DNA.

NASA GeneLab Project: Bridging Space Radiation Omics with Ground Studies.

Accurate assessment of risks of long-term space missions is critical for human space exploration. It is essential to have a detailed understanding of the biological effects on humans living and working in deep space. Ionizing radiation from galactic cosmic rays (GCR) is a major health risk factor for astronauts on extended missions outside the protective effects of the Earth's magnetic field. Currently, there are gaps in our knowledge of the health risks associated with chronic low-dose, low-dose-rate ionizing radiation, specifically ions associated with high (H) atomic number (Z) and energy (E). The NASA GeneLab project ( https://genelab.nasa.gov/ ) aims to provide a detailed library of omics datasets associated with biological samples exposed to HZE. The GeneLab Data System (GLDS) includes datasets from both spaceflight and ground-based studies, a majority of which involve exposure to ionizing radiation. In addition to detailed information on radiation exposure for ground-based studies, GeneLab is adding detailed, curated dosimetry information for spaceflight experiments. GeneLab is the first comprehensive omics database for space-related research from which an investigator can generate hypotheses to direct future experiments, utilizing both ground and space biological radiation data. The GLDS is continually expanding as omics-related data are generated by the space life sciences community. Here we provide a brief summary of the space radiation-related data available at GeneLab.

Cellular responses and gene expression profile changes due to bleomycin-induced DNA damage in human fibroblasts in space.

Living organisms in space are constantly exposed to radiation, toxic chemicals or reactive oxygen species generated due to increased levels of environmental and psychological stresses. Understanding the impact of spaceflight factors, microgravity in particular, on cellular responses to DNA damage is essential for assessing the radiation risk for astronauts and the mutation rate in microorganisms. In a study conducted on the International Space Station, confluent human fibroblasts in culture were treated with bleomycin for three hours in the true microgravity environment. The degree of DNA damage was quantified by immunofluorescence staining for γ-H2AX, which is manifested in three types of staining patterns. Although similar percentages of these types of patterns were found between flight and ground cells, there was a slight shift in the distribution of foci counts in the flown cells with countable numbers of γ-H2AX foci. Comparison of the cells in confluent and in exponential growth conditions indicated that the proliferation rate between flight and the ground may be responsible for such a shift. We also performed a microarray analysis of gene expressions in response to bleomycin treatment. A qualitative comparison of the responsive pathways between the flown and ground cells showed similar responses with the p53 network being the top upstream regulator. The microarray data was confirmed with a PCR array analysis containing a set of genes involved in DNA damage signaling; with BBC3, CDKN1A, PCNA and PPM1D being significantly upregulated in both flight and ground cells after bleomycin treatment. Our results suggest that whether microgravity affects DNA damage response in space can be dependent on the cell type and cell growth condition.

Microgravity induces proteomics changes involved in endoplasmic reticulum stress and mitochondrial protection.

On Earth, biological systems have evolved in response to environmental stressors, interactions dictated by physical forces that include gravity. The absence of gravity is an extreme stressor and the impact of its absence on biological systems is ill-defined. Astronauts who have spent extended time under conditions of minimal gravity (microgravity) experience an array of biological alterations, including perturbations in cardiovascular function. We hypothesized that physiological perturbations in cardiac function in microgravity may be a consequence of alterations in molecular and organellar dynamics within the cellular milieu of cardiomyocytes. We used a combination of mass spectrometry-based approaches to compare the relative abundance and turnover rates of 848 and 196 proteins, respectively, in rat neonatal cardiomyocytes exposed to simulated microgravity or normal gravity. Gene functional enrichment analysis of these data suggested that the protein content and function of the mitochondria, ribosomes, and endoplasmic reticulum were differentially modulated in microgravity. We confirmed experimentally that in microgravity protein synthesis was decreased while apoptosis, cell viability, and protein degradation were largely unaffected. These data support our conclusion that in microgravity cardiomyocytes attempt to maintain mitochondrial homeostasis at the expense of protein synthesis. The overall response to this stress may culminate in cardiac muscle atrophy.