According to microecological theory, the ecological balance of the nervous system, endocrine system, metabolic system, microbial system and immune system is the essence of life. The thyroid is an important endocrine gland that regulates multiple metabolic pathways by producing thyroid hormones and has a great impact on the regulation of the internal environment. Several investigations have indicated that the gut microbiota may indirectly participate in disease development and might become a promising biomarker for immune-related adverse events. Studies have found that intestinal flora are related to various thyroid diseases, and the “thyroid-gut axis” theory has been proposed accordingly[4,5]. The mechanism by which the thyroid-gut axis maintains homeostasis under external stimulation deserves further study.
As one of the environmental stressors affecting people worldwide, low-dose ionizing radiation (LDR) exposure is usually from medical radiotherapy, radiation-contaminated areas and human space exploration[6,7]. The studies previously reported have indicated that the biological responses of LDR exposure on epidemiology, animals, and cells are paradoxical due to discrepancies in the dose rate. The role of LDR in inducing apoptosis, stimulating antioxidant capacity and repairing DNA damage has been well established. Conversely, LDR exposure from radiotherapy for the treatment of autoimmune diseases or carcinoma is associated with an increasing risk of cardiovascular disease or secondary carcinoma. For the thyroid, the levels of thyroid stimulating hormone (TSH), free triiodothyronine (fT3) and free thyroxine (fT4) have been reported to be affected by LDR exposure. In addition, changes between external irradiation and the gut microbiota have been preliminarily confirmed. Thus, we want to know how the thyroid-gut axis responds to radiation exposure and what roles the thyroid and gut microbiota play in an individual’s response to radiation injury.
The objective of our research was to investigate changes in thyroid hormone levels and gut microbiota composition in LDR-exposed mice by enzyme-linked immunosorbent assay (ELISA) and 16S rRNA gene sequencing, respectively. Our research revealed the abundances of microbial species in the gut ecosystem following dose-modulated LDR exposure. These investigations predicted possible LDR-induced gut microbiota disorders based on variations in microbial species abundance. Our research aims to explore the influence of LDR on homeostasis from the “thyroid-gut axis” perspective for the first time and may provide a concrete conceptual and analytical basis for the additional study into the long-term impact of LDR on human health.
MATERIALS AND METHODS
Mice and samples
Thirty 6-wk-old female wild-type C57BL/6J mice were divided equally into five groups after one week of adaptive feeding at the Science and Innovation Center of Shanghai Tenth People’s Hospital (ID Number: SHDSYY-2022-4381). The low-dose irradiation conditions were as follows: Dose rate, 0.05 Gy/min; 1 min/time; and once every other day. Irradiation was completed by the Irradiation Center of the Institute of Radiation Medicine, Tongji University (160 keV, 1 mA). The two groups were the test group and control group. And the test group was further divided into four groups according to the radiation dose, namely, the 0.05 Gy group, 0.10 Gy group, 0.15 Gy group, and 0.20 Gy group. The general condition of the mice was observed during the experiment.
Assays of thyroid function and acquisition of thyroid tissue
Peripheral blood was collected from the mouse tail vein to measure thyroid hormone levels when the mice completed the planned dose of radiation exposure within 24 h. Serum levels of TSH, fT4, and fT3 were assayed by an ELISA kit (Rayto RT-6100). Histological sections of mouse thyroid tissue were prepared by paraffin embedding and stained with hematoxylin and eosin (HE).
DNA extraction and 16S rRNA gene sequencing
Stool samples were stored in a sterile container and snap-frozen with liquid nitrogen following collection at -80 °C (Sarstedt, 80.734.311, Germany). Using the E.Z.N.A. Soil Kit for DNA extraction (Omega, Bio-Tek, Norcross, GA, United States) and using a Nanodrop 2000 to measure the concentration of bacterial DNA (Thermo Scientific, Wilmington, United States). All the fecal samples were stored at -20 °C until ready for sequencing. Then, the V3-V4 region of the bacterial 16S rRNA was amplified by polymerase chain reaction with barcode index primers using TransStart FastPfu DNA Polymerase (TransGen BioTech, AP221-02, Beijing, China) on an ABI GeneAmp 9700 polymerase chain reaction system (United States). Amplicons were then purified by gel extraction (AxyPrep DNA Gel Extraction Kit, Axygen Biosciences, United States) and quantified using QuantiFluor-ST (Promega, United States). The adapter sequences were added to the purified amplicons using a TruSeqTM DNA Sample Prep Kit (Illumina, San Diego, United States), and paired-end sequencing was performed using an Illumina MiSeq instrument (Illumina, San Diego, United States) (Figure 1).
Figure 1 Graphical abstract.
Schematic showing sample workflow from sequencing to analysis. TSH: Thyroid stimulating hormone; ELISA: Enzyme-linked immunosorbent assay; ROC: Receiver operating characteristic; KEGG: Kyoto Encyclopedia of Genes and Genomes; LEfSe: Linear discriminant analysis effect size.
The paired-end readings obtained by MiSeq were first merged by overlapping sequences. Filtering and quality control allow the sequences to be optimized. Operational taxonomic units (OTUs) were clustered using a 97% similarity cutoff with UPARSE (version 10.0). Cohort was randomly generated for a cross-sectional study that included 24 mice exposed to radiation and 6 controls. The alpha diversity was analyzed in terms of the species richness, including CHAO, ACE, SHANNON, and SIMPSON indices. Beta diversity was used to calculate structural variation by principal coordinates analysis (PCoA) and nonmetric multidimensional scaling (NMDS) and was mapped from principal coordinate analysis. Functional pathway analysis was conducted in accordance with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database based on our microbial information.
Data that fit the normal distribution are expressed as the mean ± SD, and nonnormally distributed data are expressed as the median (interquartile). The t test or rank-sum test was used for comparing the thyroid hormone levels between the control group (n = 6) and test group (n = 24). Significant differences in the diversity index and microbial abundance between groups was determined by the Mann-Whitney U test. The microbial taxa with differences between the two groups were identified by linear discriminant analysis (LDA) effect size (LEfSe), with an LDA threshold value of > 2.0. The correlation between microbiota abundance and thyroid function was analyzed by Spearman correlation coefficient. The area under the curve (AUC) reflected the diagnostic accuracy of the bacterial indicators based on the receiver operating characteristic (ROC) curve. Unless otherwise stated, all analyses were performed using SPSS Statistics (Version 20.0). Statistical significance was set at P < 0.05 for analyses.
LDR caused changes in TSH levels in the mice
To understand changes in thyroid function in the irradiated mice, the levels of thyroid hormones in peripheral blood were determined by ELISA, including TSH, fT3, and fT4. As shown in Figure 2A, the TSH level of mice in the test mice (1.487 ± 0.347 mU/L) was obviously lower than the control mice (1.799 ± 0.289 mU/L), t = 2.188, P = 0.037). The level of TSH was the lowest in the 0.05 Gy group (1.262 ± 0.254 mU/L), which was the lowest among the four test subgroups, and in the 0.20 Gy group, it was the highest (1.664 ± 0.264 mU/L) (Figure 2A and B). The change in thyroid hormone levels increased with increasing radiation dose and was consistent with TSH regulation in the test mice. The average levels of fT3 and fT4 in the test mice [386.560 (323.973, 442.594) ng/mL and 6.211 (5.018, 8.536) ng/mL, respectively] were higher than those in the control groups [372.664 (260.701, 384.492) ng/mL and 4.428 (3.175, 6.781) ng/mL, respectively]; nevertheless, the difference was not statistically confirmed between the two groups (Figure 2C and E). Moreover, the changes in fT3 and fT4 Levels from the different dose groups were consistent with the regulation of TSH (Figure 2D and F). Figure 3 shows thyroid tissue from the different groups of mice, and there were no obvious structural abnormalities.
Figure 2 The level of thyroid hormone in mice exposed to radiation.
A: Box and whisker plots of the levels of thyroid stimulating hormone (TSH) in the control (n = 6) and test (n = 24) groups. The TSH level of the mice in the test group was lower than that in the control mice (t = 2.188, P = 0.037); B: Histogram of the levels of TSH from the control and test groups. Statistically significant differences were not confirmed between the two groups; C: Box and whisker plots of the level of fT3 in the control (n = 6) and test groups (n = 24); D: Histogram of the level of fT3 in the control and test subgroups; E: Box and whisker plots of the level of fT4 in the control (n = 6) and test groups (n = 24); F: Histogram of the level of fT4 in the control and test subgroups. TSH: Thyroid stimulating hormone. aP < 0.05.
Figure 3 Thyroid tissue from the mice exposed to radiation (HE, × 400).
A-D: Thyroid tissue of radiation-exposed mice from the 0.05 Gy group, 0.10 Gy group, 0.15 Gy group, and 0.20 Gy group. TSH: Thyroid stimulating hormone.
Overview of the 16S rRNA gene sequencing data
An overview of the thyroid hormone levels and taxonomic data of all mice is visualized in Figure 4. High-throughput sequencing results of the 16S rRNA V3-V4 hypervariable regions obtained from stool samples revealed that the tested taxa were abundant, including 11 phyla, 50 families, 110 genera, and 424 OTUs. Rank-abundance curves according to OTU levels indicated that the species were distributed evenly and richly (Figure 5A). Figure 5B indicates that the curves flattened as increasing the sample size, demonstrating that the sequencing data were adequate and appropriately capturing microbial diversity.
Figure 4 Taxonomic features of the thyroid hormone level and gut microbiota in the control mice (n = 6) and radiation-exposed group (n = 24).
A-D: The thyroid hormone level for all subjects (A), with an overview of the 5 most abundant phyla (B), 10 most abundant families (C) and 15 most abundant genera (D) in fecal samples identified by 16S rRNA gene sequencing. The taxa levels were determined by summing the operational taxonomic units and are presented by relative abundance.
Figure 5 Bioinformatic analysis of 16S rRNA gene sequencing data in the exploratory cohort at the operational taxonomic units level.
A and B: Rank-abundance curves and rarefaction curves (SOB index) of the gut microbiota; C-F: Box and whisker plots of the alpha diversity indices for richness (Chao index and Ace index) and diversity (Shannon index and Simpson index) of the bacterial communities in control mice (n = 6) and the test mice (n = 24); G-J: Box and whisker plots of the alpha diversity indices for richness (Chao index and Ace index) and diversity (Shannon index and Simpson index) of the bacterial communities among the test subgroups (0.05 Gy, 0.10 Gy, 0.15 Gy, and 0.20 Gy; n = 6 in each group); K and L: Principal coordinates analysis (PCoA) and nonmetric multidimensional scaling (NMDS) analyses indicated that the microbial diversity differed nonsignificantly between the control and test groups based on the unweighted UniFrac distance (PCoA, R = -0.065, P = -0.715; NMDS, R = -0.065, P = 0.695); M and N: PCoA and NMDS analyses indicated that the microbial diversity differed nonsignificantly among the test subgroups based on the unweighted UniFrac distance (PCoA, R = 0.004, P = 0.24; NMDS, R = 0.215, P = 0.695).
LDR does not cause gut microbiota diversity changes in irradiated mice
According to the outcome of alpha diversity analysis, the results showed that the test group had no differences in the diversity and abundances of certain microbes compared with the control group. The Chao and Ace indices, reflecting richness, indicated no considerable difference was found between the two groups (P > 0.05) (Figure 5C and D). Moreover, the Shannon and Simpson indices, which represent alpha diversity, indicated no meaningful difference between the two groups (P > 0.05, Figure 5E and F). Figure 5G-J shows the Chao, Ace, Shannon and Simpson indices among the different subgroups from the test group.
PCoA and NMDS analyses of the unweighted UniFrac distance were conducted on evaluation the overall microbiota diversity. PCoA indicated that there were no statistically significant variation in diversity of the microbiota between the test and control mice (R = -0.065, P = 0.715, Figure 5K). The NMDS results demonstrated no differences between the two groups (R = -0.065, P = 0.695, Figure 5L). Figure 5M and N show the results of the PCoA and NMDS analyses in the four test subgroups. These findings suggest that LDR in this study does not lead to altered microbiota diversity.
Changes in the composition of the gut microbiota in irradiated mice
The taxon compositions were further determined in the two groups. Eleven phyla, 50 families, and 110 genera were present in the samples from both groups. Figure 6A shows that Firmicutes, Bacteroidetes, Patescibacteria, and Deferribacteres were the four dominant phyla, with total abundances accounting for 97.83% ± 1.23% and 97.68% ± 0.75% of the microbiota of the test and control groups, respectively. The two dominant phyla were Firmicutes and Bacteroidetes in both the test and control group. At the phylum level, no meaningful variation in the relative abundances of the microorganisms was observed between the two groups (Figure 6B). We found that the Firmicutes-to-Bacteroidetes (F/B) ratio in the test group (median = 1.623 ± 0.679) was higher than controls (median = 1.266 ± 0.345), the differences were not confirmed by statistics. In addition, the microbiota was further analyzed according to the distribution at the family and genus levels.
Figure 6 The gut microbiota composition in the control and test groups.
A: Bar plots show the gut microbiota composition at the phylum level; B: Box and whisker plots of the Firmicutes/Bacteroidetes ratio in the control (n = 6) and test groups (n = 24), and there was no difference between the groups; C: Bar plots show the gut microbiota composition at the family level; D: Bar plots show the gut microbiota composition at the genus level.
LEfSe analysis results indicated that the control group had relatively high abundance in two families and four genera (LDA > 2, P < 0.05, Figure 7A). Muribaculaceae, Lachnospiraceae, Ruminococcaceae, Lactobacillaceae and Bacteroidaceae were the five dominant families, with total abundances accounting for 83.04% ± 4.59% and 80.68% ± 6.93% of the microbiomes of the test group and control group, respectively (Figure 7A). At the family level, the test group had significantly lower levels of Moraxellaceae and Enterobacteriaceae (all P < 0.05, Figure 7B). We also found that the abundance of Lachnospiraceae in the test group (20.05% ± 10.79%) was significantly higher than that in the controls (15.08% ± 10.02%), which increased in a dose-dependent manner following LDR exposure (Figure 7B and C).
Figure 7 Altered gut microbiota composition in the radiation-exposed group compared with the control group.
A: Identification of the significantly different species in the radiation-exposed group and controls at the family and genus levels as determined by linear discriminant analysis (LDA) effect size (LEfSe). Only taxa meeting a significant linear discriminant analysis threshold value of > 2.0 and P < 0.05 are shown; B: Comparison of Moraxellaceae, Enterobacteriaceae, and Lachnospiraceae abundances between the test group and control group at the family level; C: Histogram of the abundance of Lachnospiraceae in the control and test subgroups; D: Comparison of the differential bacterial abundances between the test group and controls at the genus level, including uncultured_bacterium_g_Acinetobacter, uncultured_bacterium_o_Mollicutes_RF39, uncultured_bacterium_g_Citrobacter, and uncultured_bacterium_g_Lactococcus. aP < 0.05.
Uncultured_bacterium_f_Muribaculaceae, Lactobacillus, uncultured_bacterium_f_Lachnospiraceae, Bacteroides, and uncultured_bacterium_g_Ruminococcaceae_UCG-014 were the five dominant genera, with total abundances accounting for 59.87% ± 11.92% and 67.05% ± 10.38% of the microbiomes of the test group and control group, respectively (Figure 6C and D). At the genus level, the test group had significantly lower levels of uncultured_bacterium_g_Acinetobacter, uncultured_bacterium_o_ Mollicutes_RF39, uncultured_bacterium_g_Citrobacter, and uncultured_bacterium_g_Lactococcus (all P < 0.05, Figure 7D).
ROC curve analysis reveals the diagnostic value of specific intestinal microbiota
Based on the altered gut microbiota composition between the test and control groups, we found that six kinds of bacteria might serve as diagnostic biomarkers. At the family level, Moraxellaceae and Enterobacteriaceae are the two different bacteria between the two groups, and we selected them to serve as diagnostic biomarkers. ROC curves were generated to evaluate the diagnostic accuracy of these two families [AUC value of 0.819 (95% confidence interval: 0.640-0.999)] (Figure 8A).
Figure 8 The diagnostic accuracy of the bacteria based on the linear discriminant analysis effect size results.
ROC: Receiver operating characteristic.
At the genus level, we found four different bacteria and selected them to serve as diagnostic biomarkers, including uncultured_bacterium_g_Acinetobacter, uncultured_bacterium_o_Mollicutes_RF39, uncultured_bacterium_g_Citrobacter, and uncultured_bacterium_g_Lactococcus. ROC curves were generated to evaluate the diagnostic accuracy of these four genera [AUC value of 0.882 (confidence interval: 0.728-1.000)]. At the same time, the accuracy of the prognostic variables of two families and four genera was also analyzed by ROC curves [AUC value of 0.854 (confidence interval: 0.673-1.000)] (Figure 8).
The test groups included the 0.05 Gy group, 0.10 Gy group, 0.15 Gy group, and 0.20 Gy group according to the total radiation dose the mice received. On the basis of these results, we investigated the following relationships among thyroid function tests, absorbed dose, and levels of the gut microorganisms: f_Moraxellaceae, f_Enterobacteriaceae, uncultured_bacterium_g_Acinetobacter, uncultured_ bacterium_o_Mollicutes_RF39, uncultured_bacterium_g_Citrobacter, and uncultured_bacterium_g_Lactococcus. We made an attempt to demarcate the functional influence of these flora according to the KEGG database to determine which biological pathways they may participate. Spearman correlation coefficient testing showed that Moraxellaceae (r = -0.789, P = 7.75E-7), Enterobacteriaceae (r=-0.662, P = 2.92E-4), uncultured_bacterium_g_Acinetobacter (r = -0.789, P = 7.75E-7), uncultured_bacterium_o_Mollicutes_RF39 (r = -0.365, P = 0.0106), and uncultured_bacterium_g_Lactococcus (r = -0.667, P = 4.77E-3) levels were significantly negatively correlated with the absorbed dose (Figure 9).
Figure 9 The relationship between the dose and predictive bacteria as well as thyroid function tests.
Positive correlations are represented in red, negative correlations are colored in blue. aP < 0.05; bP < 0.01; cP < 0.001.
The OTU abundance table was standardized to acquire information about the prediction of metabolic function in the KEGG pathway (Figure 10). According to the effect size, the top five pathways included carbohydrate metabolism, nucleotide metabolism, glycan biosynthesis and metabolism, xenobiotics biodegradation and metabolism, and metabolism of other amino acids. The effect size of the genetic information processing function is second only to metabolic function, which includes translation, replication and repair, folding, sorting and degradation, and transcription.
Figure 10 The histogram shows the Kyoto Encyclopedia of Genes and Genomes categories that differed significantly between the test group and controls.
A: Metabolic function; B: Processing function.
As an important organ of the endocrine system, the thyroid participates in homeostatic regulation. Recent studies demonstrated that the gut microbiota, a sophisticated and metabolically active community, is closely related to homeostasis in the host. Although LDR exposure has been demonstrated to induce contradictory effects previously, i.e., beneficial or harmful, the impacts of LDR exposure on the “thyroid-gut axis” have yet to be clarified. Therefore, we analyzed the association among LDR exposure, thyroid function and the gut microbiota in mice in this research. Our group indicated that as the LDR radiation dose increased, subclinical thyroid hyperthyroidism appeared in mice, and the observed species and OTU abundances showed no obvious changes in the radiated groups compared with the controls. It was thought that the low doses of radiation applied in this paper have an influence on the thyroid function of mice but are unable to alter the richness of the intestinal flora.
The thyroid is a radiation-sensitive organ, and internal radionuclide contamination external radiation and exposure (such as therapeutic radiation) may lead to thyroid dysfunction[14,15]. Our research found that the TSH level in the test mice was lower than the controls. Double-strand breaks can lead to cytotoxic, mutagenic and carcinogenic effects if unrepaired or misrepaired, although they only constitute a minority of DNA lesions induced by LDR, which will affect cell function[16,17]. Ishizaki et al demonstrated a very small increase in γH2AX-positive foci in a human diploid cell line after exposure to LDR, which indicates the presence of DNA double-strand breaks. Some studies have also suggested that radiation-induced thyroid dysfunction is induced by damage to small thyroid vessels and to the gland capsule, the mechanisms of which include direct thyroid cell injury from radiation damage. Because of the low thyroid dose in patients treated with total body irradiation, the incidence of Graves’ disease accompanying this management was also reported in a total body irradiation series. Strikingly, our results indicate pronounced modulation of the TSH level in mice exposed to 0.05 Gy and 0.10 Gy, while those exposed to a much higher dose (0.20 Gy) consistently clustered with controls in terms of thyroid hormone levels. We speculate that there is an intricate dose-dependent response in thyroid function to ionizing radiation, with enhanced sensitivity for the lowest doses employed here. Therefore, thyroid function and microbial communities might be directly affected by LDR exposure.
Firmicutes and Bacteroidetes were the two dominant microbiota in these two groups. We found that the Firmicutes proportion in the irradiated group was higher than in the controls, while the Bacteroidetes proportion was lower in the irradiated group. One study about environmental radionuclides on the gut microbiota from Lavrinienko came to similar conclusions that the F/B ratio was increased almost 2-fold compared with an unirradiated group. Previous investigation have also revealed that patients with metabolic diseases such as obesity, diabetes, and nonalcoholic fatty liver disease have a reduced F/B ratio compared to healthy individuals[22-24]. Recently, our group demonstrated that thyroid diseases, including Graves’ disease and thyroid carcinoma, can cause changes in the F/B ratio[25,26]. Therefore, we speculate that the composition and function of microbial communities might be affected by LDR exposure directly or indirectly through changes in thyroid function.
Compared with the controls, the abundances of Moraxellaceae and Enterobacteriaceae were decreased in the test group at the family level. Obvious differences in the levels of four species were observed between the two groups at the genus level, and uncultured_bacterium_g_Acinetobacter, uncultured_ bacterium_o_Mollicutes_RF39, uncultured_bacterium_g_Citrobacter and uncultured_bacterium_g_Lactococcus were reduced in the test group. On the one hand, the decreased abundance of the bacterium might be directly caused by radiation damage. On the other hand, changes in microbial abundance inhibit the development of gastrointestinal mucositis and promote intestinal repair mechanisms. We hypothesize that metabolites from bacteria impair intestinal barrier function. By regulating bacterial abundance and reducing the release of proinflammatory factors, the intestinal barrier could repair the damage from LDR exposure. Another significant finding in this research was that the abundance of Lachnospiraceae increased in a dose-dependent manner following LDR exposure, which was in accordance with the performance outcome of the group of mice that have recovered from high-dose radiation to live a normal life span.
The synthesis of short-chain fatty acids (SCFAs) by fermenting polysaccharides in food has made members of the Lachnospira family well known. SCFAs are essential substrates for the maintenance of the intestinal epithelium, regulation of the immune system and inflammatory response[30,31]. Broad inhibition of systemic inflammation, through reduction the levels of proinflammatory cytokines or through induction of anti-inflammatory cytokine interleukin-10, mediated by SCFAs may also potentially provide radioprotection. Elevated abundances of Lachnospiraceae have been shown to be connected with the repair of gastrointestinal after exposure to radiation.
Additionally, we employed KEGG pathway analysis to characterize changes in the gut microbiota associated with the stress from LDR exposure. KEGG pathway analysis based on OTU abundance indicated that biological metabolism was predicted to influence functional activities, including carbohydrate metabolism, glycan biosynthesis and metabolism, and nucleotide metabolism. The gut microbiota composition and metabolic pathways were influenced by LDR exposure. Moreover, processing functions were second only to metabolic functions according to the effect size, including translation, replication and repair. As the most significant influence of radiation, base modifications are a common type of DNA lesion. KEGG pathway analysis may suggest a relationship between radiation damage and functional repair. Nevertheless, our observations provide circumstantial the proof of host-microbiome interactions; more investigations are required to uncover the crosstalk within LDR-responsive alterations and host-derived targets or signals.
Our research did not consider the following issues, and these limitations will also become the focus of follow-up research. First, the effects of exposure dose and recovery time on the thyroid and gut microbiota in the mice were not considered. Second, this was a descriptive study in nature and lacked evidence regarding the underlying mechanisms of the responses to LDR exposure, which could be used to determine the alterations in metabolites from the gut microbiota. Finally, our experiment was a single-center study that lacked external validation. Despite these limitations, our research is one of the few to indicate the effects of LDR on homeostasis from the “thyroid-gut axis” perspective. In summary, our study constitutes a baseline analysis that can not only be used to further study the effects of LDR on human health but also provide potential targets for preventing or predicting radiation damage.
The thyroid-gut axis has a great influence on the maintenance of human health.
We know very little about the effects of low-dose ionizing radiation (LDR) on thyroid hormone levels and gut microbiota composition.
To investigate the potential effects by administering low-dose X-ray radiation to male C57BL/6J mice.
Peripheral blood was collected for enzyme-linked immunosorbent assay (ELISA), and stool samples were taken for 16S ribosomal RNA (rRNA) gene sequencing after irradiation.
We found that LDR caused changes in thyroid stimulating hormone (TSH) levels in the irradiated mice, suggesting a dose-dependent response in thyroid function to ionizing radiation. No changes in the diversity and richness of the gut microbiota were observed in the LDR-exposed group in comparison to the controls. The abundance of Moraxellaceae and Enterobacteriaceae decreased in the LDR-exposed groups compared with the controls, and the Lachnospiraceae abundance increased in a dose-dependent manner in the radiated groups. And the abundances of uncultured_bacterium_g_Acinetobacter, uncultured_bacterium_o_Mollicutes_RF39, uncultured_bacterium_g_Citrobacter, and uncul-tured_ bacterium_g_Lactococcus decreased in the radiated groups at the genus level, which showed a correlation with radiation exposure and diagnostic efficacy. Analysis of functional metabolic pathways revealed that biological metabolism was predicted to have an effect on functional activities, such as nucleotide metabolism, carbohydrate metabolism, and glycan biosynthesis and metabolism. Furthermore, KEGG pathway annotation also suggested that changes in the gut microbiota were related to processing functions, including translation, replication and repair.
LDR can change thyroid function and the gut microbiota, and changes in the abundances of bacteria are correlated with the radiation dose.
We administered low-dose X-ray radiation to male C57BL/6J mice to investigate the potential effects. Peripheral blood was collected for ELISA, and stool samples were taken for 16S rRNA gene sequencing after irradiation.