Longitudinal multi-omic changes in the transcriptome and proteome of peripheral blood cells after a 4 Gy total body radiation dose to Rhesus macaques.

BACKGROUND: Non-human primates, such as Rhesus macaques, are a powerful model for studies of the cellular and physiological effects of radiation, development of radiation biodosimetry, and for understanding the impact of radiation on human health. Here, we study the effects of 4 Gy total body irradiation (TBI) at the molecular level out to 28 days and at the cytogenetic level out to 56 days after exposure. We combine the global transcriptomic and proteomic responses in peripheral whole blood to assess the impact of acute TBI exposure at extended times post irradiation. RESULTS: The overall mRNA response in the first week reflects a strong inflammatory reaction, infection response with neutrophil and platelet activation. At 1 week, cell cycle arrest and re-entry processes were enriched among mRNA changes, oncogene-induced senescence and MAPK signaling among the proteome changes. Influenza life cycle and infection pathways initiated earlier in mRNA and are reflected among the proteomic changes during the first week. Transcription factor proteins SRC, TGFß and NFATC2 were immediately induced at 1 day after irradiation with increased transcriptional activity as predicted by mRNA changes persisting up to 1 week. Cell counts revealed a mild / moderate hematopoietic acute radiation syndrome (H-ARS) reaction to irradiation with expected lymphopenia, neutropenia and thrombocytopenia that resolved within 30 days. Measurements of micronuclei per binucleated cell levels in cytokinesis-blocked T-lymphocytes remained high in the range 0.27-0.33 up to 28 days and declined to 0.1 by day 56. CONCLUSIONS: Overall, we show that the TBI 4 Gy dose in NHPs induces many cellular changes that persist up to 1 month after exposure, consistent with damage, death, and repopulation of blood cells.

Comparison of Isolated Lymphocyte and Whole Blood-Based CBMN Assays for Radiation Triage.

Following a mass-casualty nuclear/radiological event, there will be an important need for rapid and accurate estimation of absorbed dose for biological triage. The cytokinesis-block micronucleus (CBMN) assay is an established and validated cytogenetic biomarker used to assess DNA damage in irradiated peripheral blood lymphocytes. Here, we describe an intercomparison experiment between two biodosimetry laboratories, located at Columbia University (CU) and Health Canada (HC) that performed different variants of the human blood CBMN assay to reconstruct dose in human blood, with CU performing the assay on isolated lymphocytes and using semi-automated scoring whereas HC used the more conventional whole blood assay. Although the micronucleus yields varied significantly between the two assays, the predicted doses closely matched up to 4 Gy - the range from which the HC calibration curve was previously established. These results highlight the importance of a robust calibration curve(s) across a wide age range of donors that match the exposure scenario as closely as possible and that will account for differences in methodology between laboratories. We have seen that at low doses, variability in the results may be attributed to variation in the processing while at higher doses the variation is dominated by inter-individual variation in cell proliferation. This interlaboratory collaboration further highlights the usefulness of the CBMN endpoint to accurately reconstruct absorbed dose in human blood after ionizing radiation exposure.

Assessment of Micronuclei Frequency in the Peripheral Blood of Adult and Pediatric Patients Receiving Fractionated Total Body Irradiation.

The cytokinesis-block micronucleus (CBMN) assay is an established method for assessing chromosome damage in human peripheral blood lymphocytes resulting from exposure to genotoxic agents such as ionizing radiation. The objective of this study was to measure cytogenetic DNA damage and hematology parameters in vivo based on MN frequency in peripheral blood lymphocytes (PBLs) from adult and pediatric leukemia patients undergoing hematopoietic stem cell transplantation preceded by total body irradiation (TBI) as part of the conditioning regimen. CBMN assay cultures were prepared from fresh blood samples collected before and at 4 and 24 h after the start of TBI, corresponding to doses of 1.25 Gy and 3.75 Gy, respectively. For both age groups, there was a significant increase in MN yields with increasing dose (p < 0.05) and dose-dependent decrease in the nuclear division index (NDI; p < 0.0001). In the pre-radiotherapy samples, there was a significantly higher NDI measured in the pediatric cohort compared to the adult due to an increase in the percentage of tri- and quadri-nucleated cells scored. Complete blood counts with differential recorded before and after TBI at the 24-h time point showed a rapid increase in neutrophil (p = 0.0001) and decrease in lymphocyte (p = 0.0006) counts, resulting in a highly elevated neutrophil-to-lymphocyte ratio (NLR) of 14.45 ± 1.85 after 3.75 Gy TBI (pre-exposure = 4.62 ± 0.49), indicating a strong systemic inflammatory response. Correlation of the hematological cell subset counts with cytogenetic damage, indicated that only the lymphocyte subset survival fraction (after TBI compared with before TBI) showed a negative correlation with increasing MN frequency from 0 to 1.25 Gy (r = -0.931; p = 0.007). Further, the data presented here indicate that the combination of CBMN assay endpoints (MN frequency and NDI values) and hematology parameters could be used to assess cytogenetic damage and early hematopoietic injury in the peripheral blood of leukemia patients, 24 h after TBI exposure.

Application of the Cytokinesis-Block Micronucleus Assay for High-Dose Exposures Using Imaging Flow Cytometry.

The cytokinesis-block micronucleus assay is a well-established method to assess radiation-induced genetic damage in human cells. This assay has been adapted to imaging flow cytometry (IFC), allowing automated analysis of many cells, and eliminating the need to create microscope slides. Furthermore, to improve the efficiency of assay performance, a small-volume method previously developed was employed. Irradiated human blood samples were cultured, stained, and analyzed by IFC to produce images of the cells. Samples were run using both manual and 96-well plate automated acquisition. Multiple parameter-based image features were collected for each sample, and the results were compared to confirm that these acquisition methods are functionally identical. This paper details the multi-parametric analysis developed and the resulting calibration curves up to 10 Gy. The calibration curves were created using a quadratic random coefficient model with Poisson errors, as well as a logistic discriminant function. The curves were then validated with blinded, irradiated samples, using relative bias and relative mean square error. Overall, the accuracy of the dose estimates was adequate for triage dosimetry (within 1 Gy of the true dose) over 90% of the time for lower doses and about half the time for higher doses, with the lowest success rate between 5 and 6 Gy where the calibration curve reached its peak and there was the smallest change in MN/BNC with dose. This work describes the application of a novel multi-parametric analysis that fits the calibration curves and allows dose estimates up to 10 Gy, which were previously limited to 4 Gy. Furthermore, it demonstrates that the results from samples acquired manually and with the autosampler are functionally similar.

High-throughput measurement of double strand break global repair phenotype in peripheral blood mononuclear cells after long-term cryopreservation.

Peripheral blood mononuclear cells (PBMCs) are a useful model for biochemical assays, particularly for etiological studies. We describe here a method for measuring DNA repair capacity (DRC) in archival cryogenically preserved PBMCs. To model DRC, we measured γ-H2AX repair kinetics in thawed PBMCs after irradiation with 3 Gy gamma rays. Time-dependent fluorescently labeled γ-H2AX levels were measured at five time points from 1 to 20 h, yielding an estimate of global DRC repair kinetics as well as a measure of unrepaired double strand breaks at 20 h. While γ-H2AX levels are traditionally measured by either microscopy or flow-cytometry, we developed a protocol for imaging flow cytometry (IFC) that combines the detailed information of microscopy with the statistical power of flow methods. The visual imaging component of the IFC allows for monitoring aspects such as cellular health and apoptosis as well as fluorescence localization of the γ-H2AX signal, which ensures the power and significance of this technique. Application of a machine-learning based image classification improved flow cytometry fluorescent measurements by identifying apoptotic cells unable to undergo DNA repair. We present here DRC repair parameters from 18 frozen archival PBMCs and 28 fresh blood samples collected from a demographically diverse cohort of women measured in a high-throughput IFC format. This thaw method and assay can be used alone or in conjunction with other assays to measure etiological phenotypes in cryogenic biobanks of PBMCs.

Cytogenetic Damage of Human Lymphocytes in Humanized Mice Exposed to Neutrons and X Rays 24 h After Exposure.

Detonation of an improvised nuclear device highlights the need to understand the risk of mixed radiation exposure as prompt radiation exposure could produce significant neutron and gamma exposures. Although the neutron component may be a relatively small percentage of the total absorbed dose, the large relative biological effectiveness (RBE) can induce larger biological DNA damage and cell killing. The objective of this study was to use a hematopoietically humanized mouse model to measure chromosomal DNA damage in human lymphocytes 24 h after in vivo exposure to neutrons (0.3 Gy) and X rays (1 Gy). The human dicentric and cytokinesis-block micronucleus assays were performed to measure chromosomal aberrations in human lymphocytes in vivo from the blood and spleen, respectively. The mBAND assay based on fluorescent in situ hybridization labeling was used to detect neutron-induced chromosome 1 inversions in the blood lymphocytes of the neutron-irradiated mice. Cytogenetics endpoints, dicentrics and micronuclei showed that there was no significant difference in yields between the 2 irradiation types at the doses tested, indicating that neutron-induced chromosomal DNA damage in vivo was more biologically effective (RBE ∼3.3) compared to X rays. The mBAND assay, which is considered a specific biomarker of high-LET neutron exposure, confirmed the presence of clustered DNA damage in the neutron-irradiated mice but not in the X-irradiated mice, 24 h after exposure.

DNA damage response in peripheral mouse blood leukocytes in vivo after variable, low-dose rate exposure.

Environmental contamination and ingestion of the radionuclide Cesium-137 (137Cs) is a large concern in fallout from a nuclear reactor accident or improvised nuclear device, and highlights the need to develop biological assays for low-dose rate, internal emitter radiation. To mimic low-dose rates attributable to fallout, we have developed a VAriable Dose-rate External 137Cs irradiatoR (VADER), which can provide arbitrarily varying and progressive low-dose rate irradiations in the range of 0.1-1.2 Gy/day, while circumventing the complexities of dealing with radioactively contaminated biomaterials. We investigated the kinetics of mouse peripheral leukocytes DNA damage response in vivo after variable, low-dose rate 137Cs exposure. C57BL/6 mice were placed in the VADER over 7 days with total accumulated dose up to 2.7 Gy. Peripheral blood response including the leukocyte depletion, apoptosis as well as its signal protein p53 and DNA repair biomarker γ-H2AX was measured. The results illustrated that blood leukocyte numbers had significantly dropped by day 7. P53 levels peaked at day 2 (total dose = 0.91 Gy) and then declined; whereas, γ-H2AX fluorescence intensity (MFI) and foci number generally increased with accumulated dose and peaked at day 5 (total dose = 2.08 Gy). ROC curve analysis for γ-H2AX provided a good discrimination of accumulated dose < 2 Gy and ≥ 2 Gy, highlighting the potential of γ-H2AX MFI as a biomarker for dosimetry in a protracted, environmental exposure scenario.

In vitro RABiT measurement of dose rate effects on radiation induction of micronuclei in human peripheral blood lymphocytes.

Developing new methods for radiation biodosimetry has been identified as a high-priority need in case of a radiological accident or nuclear terrorist attacks. A large-scale radiological incident would result in an immediate critical need to assess the radiation doses received by thousands of individuals. Casualties will be exposed to different doses and dose rates due to their geographical position and sheltering conditions, and dose rate is one of the principal factors that determine the biological consequences of a given absorbed dose. In these scenarios, high-throughput platforms are required to identify the biological dose in a large number of exposed individuals for clinical monitoring and medical treatment. The Rapid Automated Biodosimetry Tool (RABiT) is designed to be completely automated from the input of blood sample into the machine to the output of a dose estimate. The primary goal of this paper was to quantify the dose rate effects for RABiT-measured micronuclei in vitro in human lymphocytes. Blood samples from healthy volunteers were exposed in vitro to different doses of X-rays to acute and protracted doses over a period up to 24 h. The acute dose was delivered at ~1.03 Gy/min and the low dose rate exposure at ~0.31 Gy/min. The results showed that the yield of micronuclei decreases with decreasing dose rate starting at 2 Gy, whereas response was indistinguishable from that of acute exposure in the low dose region, up to 0.5 Gy. The results showed a linear-quadratic dose-response relationship for the occurrence of micronuclei for the acute exposure and a linear dose-response relationship for the low dose rate exposure.

The RABiT: high-throughput technology for assessing global DSB repair.

At the Center for High-Throughput Minimally Invasive Radiation Biodosimetry, we have developed a rapid automated biodosimetry tool (RABiT); this is a completely automated, ultra-high-throughput robotically based biodosimetry workstation designed for use following a large-scale radiological event, to perform radiation biodosimetry measurements based on a fingerstick blood sample. High throughput is achieved through purpose built robotics, sample handling in filter-bottomed multi-well plates and innovations in high-speed imaging and analysis. Currently, we are adapting the RABiT technologies for use in laboratory settings, for applications in epidemiological and clinical studies. Our overall goal is to extend the RABiT system to directly measure the kinetics of DNA repair proteins. The design of the kinetic/time-dependent studies is based on repeated, automated sampling of lymphocytes from a central reservoir of cells housed in the RABiT incubator as a function of time after the irradiation challenge. In the present study, we have characterized the DNA repair kinetics of the following repair proteins: γ-H2AX, 53-BP1, ATM kinase, MDC1 at multiple times (0.5, 2, 4, 7 and 24 h) after irradiation with 4 Gy γ rays. In order to provide a consistent dose exposure at time zero, we have developed an automated capillary irradiator to introduce DNA DSBs into fingerstick-size blood samples within the RABiT. To demonstrate the scalability of the laboratory-based RABiT system, we have initiated a population study using γ-H2AX as a biomarker.

γ-H2AX foci are increased in lymphocytes in vivo in young children 1 h after very low-dose X-irradiation: a pilot study.

BACKGROUND: Computed tomography (CT) is an imaging modality involving ionizing radiation. The presence of γ-H2AX foci after low to moderate ionizing radiation exposure has been demonstrated; however it is unknown whether very low ionizing radiation exposure doses from CT exams can induce γ-H2AX formation in vivo in young children. OBJECTIVE: To test whether very low ionizing radiation doses from CT exams can induce lymphocytic γ-H2AX foci (phosphorylated histones used as a marker of DNA damage) formation in vivo in young children. MATERIALS AND METHODS: Parents of participating children signed a consent form. Blood samples from three children (ages 3-21 months) undergoing CT exams involving very low blood ionizing radiation exposure doses (blood doses of 0.22-1.22 mGy) were collected immediately before and 1 h post CT exams. Isolated lymphocytes were quantified for γ-H2AX foci by a technician blinded to the radiation status and dose of the patients. Paired t-tests and regression analyses were performed with significance levels set at P < 0.05. RESULTS: We observed a dose-dependent increase in γ-H2AX foci post-CT exams (P = 0.046) among the three children. Ionizing radiation exposure doses led to a linear increase of foci per cell in post-CT samples (102% between lowest and highest dose). CONCLUSION: We found a significant induction of γ-H2AX foci in lymphocytes from post-CT samples of three very young children. When possible, CT exams should be limited or avoided by possibly applying non-ionizing radiation exposure techniques such as US or MRI.

Interlaboratory comparison of high-throughput protein biomarker assay quantifications for radiation exposure classification.

In the event of a widespread radiological incident, thousands of individuals will require rapid assessment of exposure using validated biodosimetry assays to inform clinical triage. In this scenario, multiple biodosimetry laboratories may be necessary for large-volume sample processing. To meet this need, we have developed a high-throughput assay for the rapid measurement of intracellular protein biomarkers in human peripheral blood samples using an Imaging Flow Cytometry (IFC) platform. The objective of this work was to harmonize and validate the reproducibility of our blood biomarker assay for radiation exposure across three IFC instruments, two located at Columbia University (CU) and the third at Health Canada. The Center for Radiological Research (CRR) at CU served as the central laboratory and reference instrument, where samples were prepared in triplicate, labeled with two radiation responsive leukocyte biomarkers (BAX and phosphor-p53 (Ser37)), and distributed for simultaneous interrogation by each IFC. Initial tests showed that significantly different baseline biomarker measurements were generated on each instrument when using the same acquisition settings, suggesting that harmonization of signal intensities is necessary. Subsequent tests harmonized biomarker measurements after irradiation by modulating laser intensity using two reference materials: unstained samples and standardized rainbow beads. Both methods generated measurements on each instrument without significant differences between the new and references instruments, allowing for the use of one master template to quantify biomarker expression across multiple instruments. Deming regression analyses of 0-5 Gy dose-response curves showed overall good correlation of BAX and p53 values across new and reference instruments. While Bland-Altman analyses indicated low to moderate instrument biases, ROC Curve analyses ultimately show successful discrimination between exposed and unexposed samples on each instrument (AUC values > 0.85).

Development of a human peripheral blood ex vivo model for rapid protein biomarker detection and applications to radiation biodosimetry.

In the event of a widespread radiological incident, thousands of people may be exposed to a wide range of ionizing radiation. In this unfortunate scenario, there will be a need to quickly screen a large number of people to assess the amount of radiation exposure and triage for medical treatment. In our earlier work, we previously identified and validated a panel of radiosensitive protein biomarkers in blood leukocytes, using the humanized-mouse and non-human primate (NHP) models. The objective of this work was to develop a high-throughput imaging flow-cytometry (IFC) based assay for the rapid measurement of protein biomarker expression in human peripheral blood samples irradiated ex vivo. In this assay design, peripheral human blood samples from healthy adult donors were exposed to 0-5 Gy X-irradiation ex vivo and cultured for up to 2 days. Samples were stained with a cocktail of surface antigens (CD66b, CD20, and CD3), fixed and permeabilized, and intracellularly stained for BAX (Bcl-2-associated X) protein, used here as a representative biomarker. Samples were interrogated by IFC, and a uniform analysis template was created to measure biomarker expression in heterogeneous and specific leukocyte subtype populations at each time point. In this human blood ex vivo model, we show that within gated populations of leukocyte subtypes, B-cells are highly radiosensitive with the smallest surviving fraction, followed by T-cells and granulocytes. Dose-dependent biomarker responses were measured in the lymphocytes, B-, and T-cell populations, but not in the granulocytes, with dose-response curves showing increasing fold changes in BAX protein expression up to Day 2 in lymphocyte populations. We present here the successful use of this ex vivo model for the development of radiation dose-response curves of a candidate protein biomarker towards future applications of dose reconstruction and biodosimetry.

High-Throughput γ-H2AX Assay Using Imaging Flow Cytometry.

The γ-H2AX assay is a sensitive and reliable method to evaluate radiation-induced DNA double-strand breaks. The conventional γ-H2AX assay detects individual nuclear foci manually, but is labor-intensive and time-consuming, and hence unsuitable for high-throughput screening in cases of large-scale radiation accidents. We have developed a high-throughput γ-H2AX assay using imaging flow cytometry. This method comprises (1) sample preparation from small volumes of blood in the Matrix™ 96-tube format, (2) automated image acquisition of cells stained with immunofluorescence-labeled γ-H2AX using ImageStream®X, and (3) quantification of γ-H2AX levels and batch processing using the Image Data Exploration and Analysis Software (IDEAS®). This enables the rapid analysis of γ-H2AX levels in several thousand of cells from a small volume of blood with accurate and reliable quantitative measurements for γ-H2AX foci and mean fluorescence levels. This high-throughput γ-H2AX assay could be a useful tool not only for radiation biodosimetry in mass casualty events, but also for large-scale molecular epidemiological studies and individualized radiotherapy.

Improved prediction of breast cancer risk based on phenotypic DNA damage repair capacity in peripheral blood B cells.

Background: Standard Breast Cancer (BC) risk prediction models based only on epidemiologic factors generally have quite poor performance, and there have been a number of risk scores proposed to improve them, such as AI-based mammographic information, polygenic risk scores and pathogenic variants. Even with these additions BC risk prediction performance is still at best moderate. In that decreased DNA repair capacity (DRC) is a major risk factor for development of cancer, we investigated the potential to improve BC risk prediction models by including a measured phenotypic DRC assay. Methods: Using blood samples from the Breast Cancer Family Registry we assessed the performance of phenotypic markers of DRC in 46 matched pairs of individuals, one from each pair with BC (with blood drawn before BC diagnosis) and the other from controls matched by age and time since blood draw. We assessed DRC in thawed cryopreserved peripheral blood mononuclear cells (PBMCs) by measuring γ-H2AX yields (a marker for DNA double-strand breaks) at multiple times from 1 to 20 hrs after a radiation challenge. The studies were performed using surface markers to discriminate between different PBMC subtypes. Results: The parameter Fres, the residual damage signal in PBMC B cells at 20 hrs post challenge, was the strongest predictor of breast cancer with an AUC (Area Under receiver-operator Curve) of 0.89 [95% Confidence Interval: 0.84-0.93] and a BC status prediction accuracy of 0.80. To illustrate the combined use of a phenotypic predictor with standard BC predictors, we combined Fres in B cells with age at blood draw, and found that the combination resulted in significantly greater BC predictive power (AUC of 0.97 [95% CI: 0.94-0.99]), an increase of 13 percentage points over age alone. Conclusions: If replicated in larger studies, these results suggest that inclusion of a fingerstick-based phenotypic DRC blood test has the potential to markedly improve BC risk prediction.

Machine learning approach for quantitative biodosimetry of partial-body or total-body radiation exposures by combining radiation-responsive biomarkers.

During a large-scale radiological event such as an improvised nuclear device detonation, many survivors will be shielded from radiation by environmental objects, and experience only partial-body irradiation (PBI), which has different consequences, compared with total-body irradiation (TBI). In this study, we tested the hypothesis that applying machine learning to a combination of radiation-responsive biomarkers (ACTN1, DDB2, FDXR) and B and T cell counts will quantify and distinguish between PBI and TBI exposures. Adult C57BL/6 mice of both sexes were exposed to 0, 2.0-2.5 or 5.0 Gy of half-body PBI or TBI. The random forest (RF) algorithm trained on ½ of the data reconstructed the radiation dose on the remaining testing portion of the data with mean absolute error of 0.749 Gy and reconstructed the product of dose and exposure status (defined as 1.0 × Dose for TBI and 0.5 × Dose for PBI) with MAE of 0.472 Gy. Among irradiated samples, PBI could be distinguished from TBI: ROC curve AUC = 0.944 (95% CI: 0.844-1.0). Mouse sex did not significantly affect dose reconstruction. These results support the hypothesis that combinations of protein biomarkers and blood cell counts can complement existing methods for biodosimetry of PBI and TBI exposures.

Sex and dose rate effects in automated cytogenetics.

Testing and validation of biodosimetry assays is routinely performed using conventional dose rate irradiation platforms, at a dose rate of approximately 1 Gy/min. In contrast, the exposures from an improvised nuclear device will be delivered over a large range of dose rates with a prompt irradiation component, delivered in less than 1 µs, and a protracted component delivered over hours and days. We present preliminary data from a large demographic study we have undertaken for investigation of age, sex and dose rate effects on dicentric and micronucleus yields. Our data demonstrate reduced dicentric and micronucleus yields at very high dose rates. Additionally, we have seen small differences between males and females, with males having slightly fewer micronuclei and slightly more dicentrics than females, at high doses.

A machine learning method for improving the accuracy of radiation biodosimetry by combining data from the dicentric chromosomes and micronucleus assays.

A large-scale malicious or accidental radiological event can expose vast numbers of people to ionizing radiation. The dicentric chromosome (DCA) and cytokinesis-block micronucleus (CBMN) assays are well-established biodosimetry methods for estimating individual absorbed doses after radiation exposure. Here we used machine learning (ML) to test the hypothesis that combining automated DCA and CBMN assays will improve dose reconstruction accuracy, compared with using either cytogenetic assay alone. We analyzed 1349 blood sample aliquots from 155 donors of different ages (3-69 years) and sexes (49.1% males), ex vivo irradiated with 0-8 Gy at dose rates from 0.08 Gy/day to ≥ 600 Gy/s. We compared the performances of several state-of-the-art ensemble ML methods and found that random forest generated the best results, with R2 for actual vs. reconstructed doses on a testing data subset = 0.845, and mean absolute error = 0.628 Gy. The most important predictor variables were CBMN and DCA frequencies, and age. Removing CBMN or DCA data from the model significantly increased squared errors on testing data (p-values 3.4 × 10-8 and 1.1 × 10-6, respectively). These findings demonstrate the promising potential of combining CBMN and DCA assay data to reconstruct radiation doses in realistic scenarios of heterogeneous populations exposed to a mass-casualty radiological event.

Machine learning methodology for high throughput personalized neutron dose reconstruction in mixed neutron + photon exposures.

We implemented machine learning in the radiation biodosimetry field to quantitatively reconstruct neutron doses in mixed neutron + photon exposures, which are expected in improvised nuclear device detonations. Such individualized reconstructions are crucial for triage and treatment because neutrons are more biologically damaging than photons. We used a high-throughput micronucleus assay with automated scanning/imaging on lymphocytes from human blood ex-vivo irradiated with 44 different combinations of 0-4 Gy neutrons and 0-15 Gy photons (542 blood samples), which include reanalysis of past experiments. We developed several metrics that describe micronuclei/cell probability distributions in binucleated cells, and used them as predictors in random forest (RF) and XGboost machine learning analyses to reconstruct the neutron dose in each sample. The probability of "overfitting" was minimized by training both algorithms with repeated cross-validation on a randomly-selected subset of the data, and measuring performance on the rest. RF achieved the best performance. Mean R2 for actual vs. reconstructed neutron doses over 300 random training/testing splits was 0.869 (range 0.761 to 0.919) and root mean squared error was 0.239 (0.195 to 0.351) Gy. These results demonstrate the promising potential of machine learning to reconstruct the neutron dose component in clinically-relevant complex radiation exposure scenarios.

Gene expression analysis in SV-40 immortalized human corneal epithelial cells cultured with an air-liquid interface.

PURPOSE: To compare the global gene expression profile of stratified epithelia generated in vitro using simian virus 40 (SV40) immortalized human corneal epithelial cells with the previously reported gene expression of normal human corneal epithelia. METHODS: Immortalized cells expanded in submerged culture were grown in an air-liquid interface of liquid permeable collagen-coated filters to foster stratification and differentiation. Stratified epithelia displaying resistances exceeding 300 Ω·cm2 were dissolved in an RNA purification lysis buffer. Purified RNA was used to globally determine gene expression levels using high-density single-channel oligonucleotide microarrays. Raw hybridization readings were converted into relative gene expression levels using Robust Multi-array Average (RMA) algorithm. Expression levels for selected genes were validated by real-time RT-qPCR. The biologic significance of the gene expression profiles was interpreted with the help of several microarray software analysis tools and ad hoc thematical analysis. RESULTS: The stratified cell culture to native epithelial comparison identified over- and under-expression in 22% and 14% of the probed genes, respectively. The larger expression decreases occurred in genes intimately associated with both the stratified epithelial lineage at large such as keratin 14 and the corneal phenotype, such as keratin 12, connexin 43, aldehyde dehydrogenases (ALDHs), and paired box gene 6 (PAX6) and its whole downstream transcriptome. Overexpression related to genes associated with cell cycling stimulation. CONCLUSIONS: The results indicate that the stratified corneal epithelial cell model generated using SV40 immortalized cells may be useful only in certain research applications. Extrapolations of studies with these cells to actual tissue cells should be done with a great deal of caution.