Inspiration4 astronaut sample collection
The exposure to radiation of the Inspiration4 mission crew was 2.2 mGy-Eq based on NASA Q quality factors and 50% FAX ModelEffective dose (Lens: 2.6 mGy, Skin: 2.7 mGy, Heart: 1.8 mGy, Central Nervous System: 1.3 mGy). Blood samples were collected before (L-92, L-44, and L-3 days) and after (R+1, R+45, and R+82 days) spaceflight. TGFβ pathway, Wnt pathway, FGF pathway, PDGFα pathway, RAAS pathway, SMAD3 pathway, Angiotensin II pathway, Fibronectin pathway, Specific fibrosis genes suggested in the paper has been used for identifying DEGs from I4 PBMCs and subpopulations. Seurat was used to normalize RNA count data and calculate average expression of each gene. The Seurat FindMarker function was used to calculate I4 DEGs (P < 0.05). Average expression of those genes was used to generate heatmap and heatmaps are generated by heatmap R package. For pathway enrichment analysis, fgsea package was used with I4 DEGs and selected pathways related to TGFb, Wnt, FGF, PDGFa, RAAS, SMAD3, Angiotensin II, and Fibronectin. All human studies were done with ethical approvals with established and approved IRBs at Weill Cornell Medicine.
JAXA Cell-Free Epigenome (CFE) Study transcriptomic data
Plasma CFE RNA differential expression data and study protocols were shared through NASA’s GeneLab platform with accession number: OSD-530 47. Briefly, blood samples were collected from 6 astronauts before, during, and after the spaceflight on the ISS. Total RNA was purified from plasma samples and processed for RNA-seq analysis. Mean expression values were obtained from normalized read counts of 6 astronauts for each time point. Heatmaps were made using the custom made fibrotic genes and pathways in Table S1 on the normalized values per time point using R package pheatmap version 1.0.12. Functional cluster analysis on the TGFβ pathway was done by selecting individual gene clusters for functional analysis using the hierarchical clustering of the target gene set. Using the default analysis settings, we used the multi-query feature of g:Profiler 63 to identify functional enrichments within and across the clusters.
Murine simulated space environment experiments
15-week +/- 3-day old, C57Bl/6J wild-type (wt) female mice were purchased from Jackson Laboratories and housed at Brookhaven National Laboratory (BNL, Upton, NY). Upon arrival to BNL, mice were quarantined and acclimated to a standard 12:12 h light:dark cycle, with controlled temperature/humidity for 1-week prior to cage acclimation. Food and water were given ad libitum, and standard bedding was changed once per week. Mice were cage acclimated (n=10 mice per group; 2 mice per cage to maintain social interaction) 3-days prior to HU, followed by 14-days either normally loaded (NL) or hindlimb unloaded (HU, see details below). Irradiation was administered on day 13 and blood tissues were collected 24-hours post-irradiation and post-euthanasia by CO2 overdose, followed by cervical dislocation. Blood was collected via the abdominal IVC in EDTA-coated tubes (0.5 M) and plasma was separated by centrifugation at 1,700Xg for 15 minutes. Plasma was collected and frozen at -80°C for storage. A 100 μl aliquot of cellular fraction was frozen and stored at -80°C for RNA analyses, while the remaining cellular fraction was lysed with 1x RBC lysis buffer (Thermo Fisher Scientific) followed by flow cytometric preparation and analyses, as described below. All organs were flash frozen and collected at the time of dissection. Specifically, the heart was cut in half vertically with the left half of the heart flash frozen and the right half the heart fixed in 4% paraformaldehyde. All tissues were stored at either -80°C or 4°C. Body weight tracking was performed on days -3, 0, 7 and 14.
Hindlimb unloading (HU) was performed using the adjusted Morey-Holton method for social housing (n=2 per cage) 52. Briefly, mice were suspended from the tail using non-invasive traction tape attached to an adjustable pulley mounted on the top of a standard rat cage. The adjustable nature of the device allows the user to position the animal in a head-down position (approximately 30° to the horizon) once attached. The crossbar height, lateral crossbar position, and chain length are adjustable. The pulley is free to move along the crossbar that spans the length of the top of the cage enabling the mouse to navigate freely on one half of the cage using its forepaws and to interact with the second mouse in the cage, without enabling ambulation of the limbs. Individual food and water sources are provided to each mouse and replaced daily. Nestlets were provided to each mouse to encourage nesting. Control mice were housed in identical cages with normal ambulation.
On day 13, mice were transported from the Brookhaven Laboratory Animal Facility at BNL to the NASA Space Radiation Laboratory (NSRL) facility by animal care staff and were transferred to individual HU boxes. The following doses of irradiation were administered; simplified GCR sim (0.5 Gy), SPE (1 Gy), and Sham control (0 Gy). Mice were positioned in the plateau region of the Bragg curve and irradiated at room temperature in individual HU boxes to enable whole body irradiation while maintaining hindlimb suspension. The NSRL physics staff performed the dosimetry. To simulate GCR, we used the simplified GCR simulation of ions, energy, and doses determined by a NASA consensus formula that consists of 5 ions: protons at 1000 MeV, 28Si at 600 MeV/n, 4He at 250 MeV/n, 16O at 350 MeV/n, 56Fe at 600 MeV/n, and protons at 250 MeV. This dose of radiation is equivalent to what an astronaut is predicted to receive in deep space during a Mars mission, though it is modeled as a single exposure over 25 minutes instead of the actual chronic exposures over 1.5 years. Further, GCR sim is a mixture of high and low LET ions in a ratio of 15% to 85%, respectively. To simulate SPE, we used a total dose of 1 Gy protons with energy ranges from 50 MeV to 150 MeV. For all irradiation treatments a 60x60 beam was utilized at NSRL. For radiation dose equivalence and biological reference, we used 5 Gy gamma irradiation with the cesium resource available at BNL, in the absence of HU due to resource limitations. Sham controls were treated similarly to GCR/SPE irradiated mice, including HU cage boxes and beam line (without irradiation) for the same duration as GCR simulation, i.e. 25 minutes.
Antagomir countermeasure treatment in the mice
Inhibitors, also referred to as antagomirs, to the following three specific miRNAs were purchased and designed from AUM biotech (Philadelphia, PA, USA): miR-16-5p, miR-125b-5p, and let-7a-5p. Specifically, the antagomirs were identified as AUMantagomir oligos. Each antagomir was given at a dose of 10 mg/kg combined for a total of 30 mg/kg total AUMantagomir dose by intraperitoneal (IP) injection every 3 days starting from 1 day prior to HU setup until 1 day prior irradiation.
Antibodies
The following antibodies were used in the analyses: Recombinant Anti-TGF beta 1 antibody [EPR18163] (1:1000 in dilution, abcam, Cambridge, UK), Anti-Smad3 antibody (1:1000 in dilution, abcam, Cambridge, UK), and Anti-Collagen I antibody (1:1000 in dilution, abcam, Cambridge, UK). Vinculin (E1E9V) XP® (1:2000 in dilution, Cell Signaling Technology, MA, USA), and GAPDH (D16H11) XP® (1:8000 in dilution, Cell Signaling Technology, MA, USA) were used as housekeeping proteins. Anti-rabbit IgG, HRP-linked (1:4000 in dilution, Cell Signaling Technology, MA, USA) served as a secondary antibody for immunoblot analyses.
Histological analysis
For histological analysis, formalin-fixed paraffin-embedded (FFPE) mouse cardiac tissue was cut into 6 µm thick serial sections and dried overnight. The next day, Masson’s trichrome staining was performed to visualize the extent of fibrosis in cardiac tissue. In brief, tissue sections were deparaffinized by ROTI®Histol (Carl Roth, Germany) and rehydrated followed by the incubation with Bouin’s fixative (Carl Roth, Germany) for ≥1 h at 60°C. Thereafter the following steps were performed to achieve the appropriate staining: washing in dH2O; incubating in Weigert’s iron hematoxylin (Carl Roth, Germany); differentiation in 0.5% hydrogen chloride (HCl)/70% EtOH solution; bluing in tap water; incubation in Biebrich scarlet solution; washing in dH2O; incubation in phosphomolybdic/phosphotungstic acid, followed by aniline blue; washing in dH2O; incubating in 1% acetic acid, and a final washing step. Then the slices were dehydrated by an ethanol row finishing with 2 times Xylene (Carl Roth, Germany). Finally, the stained sections were mounted with Entellan® (Merck Millipore, MA, US) and dried overnight.
Microscopical analyses of the stained sections were performed using an inverse microscope Axiovert 200 (Axio Observer, Zeiss, Germany) and an EC Plan-NEOFLUAR 20x/0.5 objective (Zeiss, Germany). For quantitative analysis, depending on the sample size, 15 pictures were taken by Axiocam MRc3 (Zeiss, Germany) and AxioVision SE64 Rel. 4.8 (Zeiss, Germany) from each tissue sample to generate a mostly randomized collection of sections to be analyzed. To differentiate and quantify the fibrotic (blue) area in relation to the pure cardiac muscle areas (cardiomyocytes stained in red) per image, the software ImageJ (NIH, LOCI, US) was used. Therefore, a color threshold was set as follows: Hue 100-225, Saturation 0-255, and Brightness 0-180. Furthermore, to adequately assess the amount of red-colored myocardial tissue, a color deconvolution using Masson’s Trichrome Vector was performed for each picture. The threshold was then set between 0-90 or 0-150, depending on the intensity of the staining. Within each sample stained, the same threshold was used. Finally, the relative fibrotic area (RFA) was calculated by dividing the pixels of the blue area by the pixels' total area of the tissue section, expressed as a percentage.
Immunoblotting
Shock frozen mouse cardiac specimen was lysed by RIPA lysis buffer system (150 mM NaCl, 50 mM Tris-HCl, 1% Nonidet P-40, 0.5% sodium deoxycholate, 0.1% SDS) with phenylmethylsulfonyl fluoride (PMSF) and protease inhibitor. Total protein concentration was determined using BCA assay (Thermo Fisher Scientific, NH, USA). Equal quantities of protein lysates were separated by SDS-PAGE on a 4-15% Mini-PROTEAN TGX Gel (Bio-Rad Laboratories, Austria) with 4x Laemmli buffer (Bio-Rad Laboratories, Austria) serving as the loading buffer. The separated proteins were further transferred to polyvinylidene fluoride (PVDF) Immun-Blot PVDF membranes (Bio-Rad Laboratories, Austria) by tank blotting. Membranes were blocked with 5% non-fat dry milk powder in TRIS-buffered saline with Tween-20 (TBS-T) for 1 h and were incubated overnight at 4°C, as indicated above. The next day, the membranes were incubated with a HRP-linked secondary antibody for 1 h at room temperature. The chemiluminescent signal was detected using a Raytest Stella 8005 (Elysia-Raytest, Germany) and the camera control software XStella Version 1.00.011 (Elysia-Raytest, Germany). In order to normalize the proteins of interest, the housekeeping proteins Vinculin or GAPDH were further incubated overnight at 4°C and the bands were detected the next day as indicated.
Data analysis and normalization were performed in Image Lab 6.0.1 software (Bio-Rad, Vienna, Austria) by generating a multichannel image. Thereafter, the lanes and bands were detected manually, and the normalization channel (housekeeping proteins) and reference lane (reference sample) selected. Furthermore, the background subtraction was set to a disk size of 70.0 mm. The normalized quantities of the samples were calculated automatically and are given as the adjusted volume corrected by the normalization factor, named “Norm. Vol. (Int)”: Thereafter, the relative protein volume (RPV) for each sample was calculated in relation to the reference sample, which was termed as 100%. Consequently, the RPV of the proteins of interest may exceed 100% as the RVP is in relation to the 100% band of the reference sample.
Statistical analysis
The statistical analysis was performed using GraphPad PRISM 9 software (GraphPad-Software, San Diego, CA, USA). The data were tested for Gaussian approximation by Shapiro-Wilk test. Grouped analyses of the RFA and Western Blot results were performed using two-way ANOVA with Tukey’s multiple comparisons test. If the bar charts included normally distributed data, the figures are given as mean and standard deviation (SD); not normally distributed data is given as median and interquartile range (IQR; 𝐼𝑄𝑅=𝑄0.75−𝑄0.25). If the figures include both, normally distributed, as well as not normally distributed values, all datasets in the graphs are given as median and IQR. Correlation analyses were also performed according to the individual Gaussian approximation of the datasets. The nonparametric data was analyzed using Spearman correlation. In the figures of the correlation analyses, the simple linear regression lines are depicted. The alpha threshold and the confidence level were set as p≤0.05 for all analyses and the results are termed as statistically significant.