Next Article in Journal
Impact of Long-Term Pyriproxyfen Exposure on the Genetic Structure and Diversity of Aedes aegypti and Aedes albopictus in Manaus, Amazonas, Brazil
Previous Article in Journal
Molecular Cloning, Tissue Distribution and Antiviral Immune Response of Duck Src
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Systematic Review

Non-Thermal Effects of Terahertz Radiation on Gene Expression: Systematic Review and Meta-Analysis

1
School of Life Science and Technology, Xi’an Jiaotong University (XJTU), Xi’an 710049, China
2
Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi’an Jiaotong University, Xi’an 710049, China
*
Authors to whom correspondence should be addressed.
Genes 2024, 15(8), 1045; https://doi.org/10.3390/genes15081045
Submission received: 22 July 2024 / Revised: 5 August 2024 / Accepted: 7 August 2024 / Published: 8 August 2024
(This article belongs to the Section Genes & Environments)

Abstract

:
With the advancement of terahertz technology, unveiling the mysteries of terahertz has had a profound impact on the field of biomedicine. However, the lack of systematic comparisons for gene expression signatures may diminish the effectiveness and efficiency of identifying common mechanisms underlying terahertz effects across diverse research findings. We performed a comprehensive review and meta-analysis to compile patterns of gene expression profiles associated with THz radiation. Thorough bibliographic reviews were conducted, utilizing the PubMed, Embase, Web of Science, and ProQuest databases to extract references from published articles. Raw CEL files were obtained from Gene Expression Omnibus and preprocessed using Bioconductor packages. This systematic review (Registration No. CDR42024502937) resulted in a detailed analysis of 13 studies (14 papers). There are several possible mechanisms and pathways through which THz radiation could cause biological changes. While the established gene expression results are largely associated with immune response and inflammatory markers, other genes demonstrated transcriptional outcomes that may unravel unknown functions. The enrichment of genes primarily found networks associated with broader stress responses. Altogether, the findings showed that THz can induce a distinct transcriptomic profile that is not associated with a microthermal cellular response. However, it is impossible to pinpoint a single gene or family of genes that would accurately and reliably justify the patterns of gene expression response under THz exposure.

1. Introduction

Terahertz electromagnetic waves are a gap bridging the electronics and the photonics regions of the electromagnetic spectrum [1,2]. In the last two decades, researchers have realized that terahertz (THz) technology can be used for real-life applications. By definition, THz radiation has photon energy in the milli-electron volt (meV). Hence, it is physically impossible for THz radiation to ionize components of biomacromolecules such as DNA and proteins. The vibration and rotation energy levels of many biological macromolecules such as DNA, RNA, and lipids are in the THz band [3,4]. Frequencies from 0.10 to 2.53 THz are the most investigated in the published studies. This could be explained by the fact that sources and detectors are more accommodating in these ranges [5]. Frequency becomes less penetrative as its value increases. In addition, a peak absorption is observed with biological tissues at 0.1 THz due to interactions between THz waves and water molecules increasing transcriptional activity [6].
The history of THz bioeffects studies has witnessed disinterest, rebirth, and promise. For a long time, this field has remained underdeveloped and inadequately explored because of a scarcity of appropriate sources and detection methods [7,8]. The advancement of manufacturing and the increased accessibility of THz-related instruments have propelled research on THz bioeffects and the field is now gaining momentum. One crucial aspect that was initiated early on and holds significant importance is the impact of terahertz radiation on gene expression.
More than 150 articles have been published on THz medical applications on PubMed, yet only around 50 address the effects on living cells and tissue [9]. The environmental exposures to THz have increased in the last two decades, and this trend may continue in the coming years.
Traditionally, the effects following the radiation of the magnitude of THz on biological molecules and the understanding of their mechanisms have not been of real concern [10,11]. On the one hand, the results of investigations are still pending as to whether exposure to THz should be prevented to avoid biological and potential health risks associated with the radiation energy emitted by widespread devices [11]. On the other hand, the question whether they should be reduced by limiting the deployment of these technologies until adequate safety standards or margins for longer periods are reached by consensus is still being debated [12]. Therefore, answers to queries related to the biological effects of terahertz radiation exposure on organisms are worth exploring.
The choice of the range of THz electromagnetic frequency is supported by the idea that this frequency increases water molecule rotation and interacts with biomolecules such as proteins. However, those claims are difficult to verify experimentally. Output power is the amount of energy emitted and absorbed in the biological material, leading to an increase in the thermal effects. The thermal effect following THz pulses absorption is, to an extent, connected to power density [13,14]. Computational and experimental models both showed the damage threshold to be 7.16 W/cm2. Uncontrolled output power and excessive thermal exposure cause structural changes and harmful thermal effects [15].
Alteration in gene expression is one of the first steps towards changes in intracellular signaling pathways and directly modifies cell dynamics and plasticity. It is argued that transcription activity regulators and macromolecular changes can be triggered by THz and lead to gene expression shifts. THz irradiation action on the latter is thought to be catalytic and linked to several interdependent factors such as transcription factors, promoters, and the pathways involved [16].
While non-ionizing, exposure to THz frequencies may hinder oxidative stress repair mechanisms, increasing the risk of cancer through the impairment of cellular components, including DNA [17]. In the present literature, reports on the efficacy of THz radiation as a biologic modulator exhibit inconsistency on various biological levels—some studies reported strong effects to subtle differences at the DNA level and some correlations between THz and biological effects [18,19,20], where a few studies documented small changes in the differential response of genes; some suggest its promising therapeutic potential because of its ability to boost immune cells [21,22]. However, the mechanisms underlying these changes remain elusive and the data used in the models related to the subject are very preliminary.
The systematization and standardization of the results is mainly impeded by fundamental differences in the experiment protocols.
The objectives of this review were as follows: (1) conduct an updated systematic review of the available literature on the relationships between terahertz radiation exposure and gene expression; (2) summarize key parameters in the experimental studies of terahertz irradiation; (3) highlight groups of genes and the pathways that are enriched. The outcome of this review could contribute to broadening the evidence of the terahertz radiation’s effect on gene expression and downstream consequences.

2. Materials and Methods

2.1. Protocol and Registration

The protocol of this research is available in PROSPERO (CDR42024502937). This systematic review adhered to the international standard called Preferred Reporting Items for Systematic Meta-Analyses (PRISMA) [23].

2.2. Eligibility Criteria

We included original articles up until 28 February 2024, including animals, cell, and in vitro models—the studies involved investigated terahertz irradiation followed by a gene expression profile assessment between irradiated and control groups. The reporting gene expression studies must have used high-throughput RNA sequencing or a microarray for quantification. The pre-specified exclusion criteria were as follows: (1) no relevant outcome measure reported; (2) gene expression measured using qPCR, (3) not an original full research paper (reviews papers, conference proceedings, books, editorials, comments, research letters, no detailed methodology described, as well as manuscripts without an English full text available). Finally, only articles published in English were included.

2.3. Information Sources

Electronic database (PubMed, Embase, Web of Science, and ProQuest) searches were conducted between 2 January 2024 and 1 March 2024.

2.4. Search

The full research strategies used within these databases are presented in Supplementary File S1 of this manuscript. The following components were included in the search string or other subject heading terms. The search string was conducted by M.N.D. The desired references were obtained by combining the main keywords and their synonyms. For example: (1) “terahertz radiation”, “terahertz wave”, “terahertz pulse”, and (2) “gene expression”, “gene expression profiling”, and “transcriptome” were some of the compound keywords.

2.5. Study Selection

Search results across databases were merged using Endnote. After duplicates were removed, two researchers (M.N.D, S.S) separately screened titles and abstracts of retrieved records for full-text selection. In the case of a disagreement on which article’s full text should be screened, consensus was reached by discussion. The full-text screening constituted the last step of the selection process. Thereafter, in pairs, (M.N.D, S.S) we independently screened full-text articles for inclusion. Again, in the case of a disagreement or uncertainties, consensus was reached based on inclusion and exclusion criteria or by discussion. A third party was consulted if necessary.

2.6. Data Collection Process

In this step, reviewers independently recorded relevant information from studies. A customized Excel sheet was used to extract the following: (1) author name; (2) year of publication; (3) country; (4) methodological information such as (a) animals/cell type models, (b) frequency range, (c) average power, (d) peak power, (e) pulse rate, repetition rate; (5) duration of exposure; (6) time before RNA extraction; (7) method of genetic analysis; (8) fold change; (9) number of differentially expressed genes; (10) pathways analysis method; and (11) biological outcomes. The GEO [24] search was conducted using (THz) OR (Terahertz). The results were further filtered using Entry type = Series, then Organism = mus musculus, and finally Study type = Expression profiling by array. Raw CEL files were downloaded using the GEOquery package (version 2.66.0) [25]. The full search string used was as follows: (Terahertz[All Fields] OR THz[All Fields]) AND “Mus musculus”[porgn] AND (“gse”[Filter] AND “Expression profiling by array”[Filter]).

2.7. Bioinformatical Analysis

All the following steps were performed in R (version 4.2.3), from Bioconductor (version 3.16) [26]. The raw expression files of the Affymetrix and Agilent files were normalized using the fRMA function in the fRMA package (version 1.50.0) [27] and quantile normalization, respectively, in R. Probe IDs were annotated to match the corresponding gene symbol. The limma R (version 3.54.2) package [28] was used to screen for DEGs separately between the control and the exposed group. DEGs were selected to have a logFC (log fold change) > 1 and a p-value < 0.05. Over-representation was performed with clusterProfiler (version 4.6.2) [29] using C5 canonical pathways from MsigDB (version 7.5.1) [30]. This was considered significative if the adjusted p.value was less than 0.05. The R-script from pre-processing to DEGs screening and enrichment analysis can be found in Supplementary File S1.

3. Results

3.1. Study Selection

After duplicated articles were removed, a total of 625 articles were retrieved, of which 64 were from PubMed, 61 from Embase, 66 from Scopus, 191 from Web of Science Core Collection, and 243 from ProQuest. Three-hundred and eighty (380) studies were found to be eligible for the review. After titles, the abstracts and keywords were analyzed; 65 articles were sought for full-text screening. When the final eligibility was manually screened, 14 articles were found to meet the inclusion criteria. Figure 1 shows the Prisma flowchart describing the process of the systematic literature search to identify all eligible studies. Thereafter, each of these papers were thoroughly and autonomously read by one of the reviewers to identify and extract meaningful information. The datasets comprising five microarrays were the following: SuperSeries [GSE41106] (GSE41083, GSE41084, GSE41085), GSE44671, and GSE178729; these were extracted from the Gene Expression Omnibus (GEO) using the GEOquery package [25].

3.2. Study Characteristics

The studies considered and their respective parameters are summarized in Table 1.
The countries that are a part of the affiliations in the studies are China (four articles), Canada (three articles), USA (three articles), Russia (two articles), and South Korea (two articles).
The analysis of the selected primary studies showed that they were conducted on eight different cellular models and one bacterial species (E. coli); the in vitro studies corresponded to C57BL/6J and BALB/c mice. It is essential to provide as many details as possible in any exposure study about the source of irradiation. We have listed the studies and summarized all the reported parameters in the studies.
As seen in Table 1, the parameters of the THz fluctuations in the irradiation sources in the experimental settings exhibit a manifest variability in their characteristics: frequency, intensity (average and peak power density), and operation regime (pulsed and CW). All experiments were carried out under monitored thermal conditions to ensure that temperature was not related to the observed effects.
The frequency of radiation in our studies varied from 0.1 THz [31] to 10 THz [10]. Exposure time varied substantially between 10 and 12 h. The shortest time exposure was 10 min [32] and the longest was 12 h exposure [16]. In most of the studies, the biological materials were harvested at one selected time point for mRNA quantification, from 10 min [33] to 48 h [34]. In one study, mRNA was extracted 28 days after the repeated treatment of THz exposure [22].
Table 1. Study characteristics of the included studies.
Table 1. Study characteristics of the included studies.
StudyCountryBiological ModelSample (n)THz Exposure ParametersRNA Extraction Time Post-Exposure
f (THz)Pulse RatePulse DurationI Average (mW/cm2)I Peak (mW/cm2)T (°C)Δ (°C) IncreaseExposure Duration
In vitro studies
Bock et al. [10]USAMouse stem cellsDuplicates~101 kHz35 fs1~3021.61NC (0.28)9 hNA
Alexandrov et al. [16]USAMouse stem cellsTriplicates2.521 kHz35 fs~0226–27NC12 hImmediat.
Tivota et al. [32]CanadaArtificial human skin tissueQuadruplicates0.2–2.51 kHzNA573321<0.710 min30 min
Bogomazova et al. [35]RussiaHuman embryonic stem cellsDuplicates2.3CWNA110NA24~11 h2 h
Echchgadda et al. [36]USAHuman jurkat T lymphocytes cell lineTriplicates2.52CWNA532NA37~640 min4 h
Hough et al. [37,38]Canada3D human skin tissueQuadruplicates0.6>1 kHz50 fs7473.837<110 min30 min
Zhao et al. [39]ChinaRPE cells, HCE cells, Müller cellsQuadruplicates0.71 kHz0.5 ps~12~37<0.26 h15 h
Shang et al. [31]ChinaRat primary hippocampal neuronsTriplicates0.1NANA33NA37NC20 min2 h
Peltek et al. [33]RussiaE. coliDuplicates1.25–3.755.6–22.440–100 ps0.140.837215 min10 min
Zhao et al. [40]ChinaMouse cortical neurons and oligodendrocytes precursor cellsQuadruplicates3.11 kHzNA70 μW/cm2NA35NC15 mn (x3)/day–3 h/dayNA
Cheon et al. [34]South KoreaMelanoma cellsTriplicates1.61 kHz35 fs750NA38NA30 min0–48 h
Animal model studies
Kim et al. [8]South KoreaC57BL/6J and BALB/cTriplicates2.501 kHz3100.32 μW/cm2NA22NC1 h24 h
Xu et al. [22]ChinaFemale C57BL/65–50.14NANA90NARTNC10 min28 days
Notes. f = frequency; I = power density (mW/cm2); I average = average incident irradiance; I peak = peak power density; T = ambient temperature; Immediat. = immediately; h = hour; min= minutes; Δ = increase on temperature during experiment; NC = no temperature change; NA = not reported; RT = room temperature.
All studies evaluated gene expression changes but differed widely in technologies for estimating gene expression changes and statistical cut-offs. All included studies measured gene expression adequately in accordance with the standard and indicated that the specific changes noticed may primarily depend on the general responsiveness of the biological material towards THz irradiation. Microarray was revealed to be the most popular method of investigating THz bioeffect.
Nine out of thirteen studies used gene expression microarrays to assess gene expression while four studies employed next-generation RNA-sequencing techniques with the Illumina platform. All studies resulted in two distinct expression patterns consisting of upregulation and downregulation. The gene expression technology used differentially expressed genes and databases for enrichment analysis and pathway analysis results; the biological outcomes for each study reported in the reviewed studies are synthesized in Table 2.
This search in GEO outputted seven results. After manual screening, only five studies were eligible (Table 3). One was a super-series of the three studies already available (GSE41106), and the other was composed of only one sample (GSE23888); therefore, these were not considered.
We performed a functional enrichment of the selected differentially expressed genes. The unique and shared DEGs are presented in an upset plot (Figure 2). We then constructed a network of the significantly enriched pathway from the five datasets retrieved in GEO, using Limma and ClusterProfiler. The results are presented in Figure 3. Among pathways that were the most enriched were “Oxalate Membrane Transporter Activity”, “Palmar Hyperhidrosis”, “Positive Regulation of Fatty Acid Biosynthetic Process”, “Positive Regulation of Acute Inflammatory Response”, and “Positive Regulation of Fatty Acid Metabolic Process” in response to THz stress response.

4. Discussion

This review systematically summarized an updated literature review of the effects of THz on gene expression profiling studies in different groups.
Only 13 studies met our inclusion criteria. Our knowledge of how terahertz-based wave irradiation can affect live tissues at a cellular level is scarce. Interest in THz technologies continues to grow but fewer studies have investigated the biological effects of terahertz-based wave irradiation on gene expression. As of now, no reviews summarizing the present evidence to draw more acute conclusions have been conducted on studies where changes at the transcript (mRNA) level were assessed as the outcome.

4.1. Risk of Bias Assessment

The case–control study was the only study design used to explore the link between THz and gene expression changes, and therefore, no element of the study design seemed to be the origin of the various results. Thus, it showed a low risk of bias. At least two replicates were used in each study, and they all specified the number of replicates, thereby representing a low risk of bias. As for the quality of exposure, all studies reported the exposure duration, with some missing details in the other parameters, indicating a high risk of bias. With regard to temperature, changes were reportedly monitored in all studies, while one study, Cheon et al. [34], did not report temperature change, meaning a low risk of bias. RNA extraction time post-exposure was reported in all but two studies: Bock et al. [10] and Zhao et al. [40]. Overall, the assessment of the risk of bias revealed that four studies provided all the parameters, while the remaining nine studies had at least one parameter missing.

4.2. Summary of Evidence

Although pathway analyses were performed for most studies, their different pathway nomenclature made comparisons throughout studies that are not comparable one to another. For completeness and to ensure fair comparisons, the mainly reported differentially expressed genes and enriched pathways are presented by topic in terms of the main findings on cellular types.
The skin was the most used model [8,22,32,34,37,38]. It is likely that more people will have to be subjected to electromagnetic THz waves in the near future, as their use has been envisioned to be fast-growing and will likely find wider applications in categories such as patients undergoing THz diagnostics, the security and the military fields, material science, art conservation workers, and communication networks. This raises concerns about dermal tissue that is continuously radiated with THz bandwidth and the adverse ocular effects. The skin is the principal recipient organ of terahertz applied to the body [41]. The depth and width penetration of THz radiation into the skin is estimated to be 0.1 to 0.3 mm and depends on the frequency of the wavelength. Hence, probable hazards may involve the skin’s superficial cellular layers [4]. Moreover, the skin is considered a good model since it plays the role of buffer between the external and internal environment; therefore, the electromagnetic energy is highly absorbed by the skin [18]. However, terahertz absorption by the skin is a complex process. In this perspective, Titova et al. [32] tested THz effects in dermatologic conditions and inflammatory skin diseases by radiating artificial human skin tissue. The keratinization formation of the cornified envelope genes was downregulated. Their findings set grounds for the possible use of high-THz pulses for therapeutic applications by targeting expression patterns of inflammatory-regulating genes, as they found a low expression of immunologic proteins (SERPINS). In the same line of work, Hough et al. [37,38] did a series of experiments on 3D skin tissue exposed to prolonged trains of THz. Both studies showed that terahertz pulses dysregulated pathways involved in human cancer initiation, maintenance, and progression in skin cells. More recently, Cheon et al. [34], in a 48 h irradiation period, assessed DNA methylation at 1.6 THz. They stimulated the global magnitude of DNA methylation by 16% in melanoma cancerous samples. The authors found a significantly decreased mRNA expression of the JUN/FOS transcription factor and inflammation-related gene (CXCL8). The findings also revealed that THz modulates transcriptional changes with potential applications in cancer, since abnormalities in transcriptional activities are known to initiate cancer. The result is of broader health importance due to the methylation that occurred in irradiated cells. Kim et al. [8] proposed a mechanism by which THz activates wound response in skin but the initial signaling trigger remains unclear. Xu et al. [22] explored the role THz might play in driving the repair of spinal injury. THz stimulation across control and treated groups elucidated molecular mechanisms promoting the rehabilitation of SCI (Spinal Cord Injury) after 28 days of daily THz irradiation treatment for 15 min.
Concerning immune cells, Echchgadda [36] investigated changes that take place immediately in jurkat cells exposed to THz compared to thermally heated cells. The results indicated that the response induced by THz could not be perfectly mimicked by heated cells. These studies warrant further investigations to support new insights into the potential applicability of THz in clinical settings.
Concerning corneal cells, Zhao et al. [39] found that the expression of numerous genes, known as protooncogenes, which are involved in normal cell growth and anaplastic lymphoma kinase (ALK), was enhanced, and human centromere-associated protein E (CENPE) expression was the only under-expressed gene in human corneal epithelial cells.
Out of the 13 studies, 3 studies [10,16,35] covered THz effects on stem cells. Stem cells are unspecialized cells capable of infinite renewal. They are able to intrinsically differentiate into any specialized cell of an organism. Developmental potency is diminished in each step, which means perturbations of the hallmarks of stem cells may cause developmental disorders and tissue depletion [42].
Bock et al. [10] was the earliest study to explore gene expression using microarrays and the first shred of evidence of THz radiation on stem cells. These studies explored to what extent THz radiation affects the differentiation capabilities of mouse MSCs. Both investigations [10,16] resulted in expression profile changes in genes involved in differentiation. This brings about skeptical questions of whether THz radiation modifies self-renewal ability and the narrow differentiation capabilities of stem cells. They found increased expression transcriptional regulator activity induced by THz radiation through promoting the differentiation of stem cells into adipocytes. This early research indicated the non-thermal biological effects of THz and theorized that the variation in specific genes could be explained by the DNA breathing model, which put forward the claim that the responses of genes to THz radiation are prone to the promoter structure, therefore enabling transcription as a result of denaturation bubbles upon THz exposure [43]. Furthermore, it was hypothesized that THz radiation could be a biological modifier tool for the differentiation of specialized cells.
Concerning neuronal cells, two studies were found. Zhao et al. [40] showed that Cytotoxicity effects were not observed in response to THz exposure. A downstream analysis of the pathways results appeared to be in line with the morphological changes. Shang et al. [31] reported an increased expression of Aqp5 and Mmp3; meanwhile, they observed downregulation in Itpr1, Atp2b4, Apof, and Erf in neonatal rat neuron cells. The transcriptional activity of AP-1 was significantly diminished to binding sites, implying a possible direct modification of the interaction between DNA and TFs.
Peltek et al. [33] conducted the only study on bacteria. They showed that THz pulses induced a nonthermal response in gene expression. Electromagnetism leads to a stress response mediated by “pili adhesion”, “cell aggregation”, and “assembly of the septal ring”. The study revealed one characteristic of THz radiation for the targeted activation or repression of genes in bacteria.
Altogether, using the advantages of having datasets from GEO, we obtained, overall, 24 samples. A total of 415 differentially expressed genes were screened and put together. Of which there were 21 from GSE41083, 98 from GSE41084, only 4 from GSE41085, 259 from GSE44671, and finally 33 from GSE178729. The meta-analysis focused on consistent changes in every single dataset. All unique differentially expressed genes were included in the study. This method was able to capture both shared and individual signals that would otherwise be missed in single datasets. The meta-analysis demonstrated that there is a high heterogeneity change during exposure. The heterogeneity may arise from various factors, including dose dependency, differences in types of biological materials being tested, and variations in the parameters of the terahertz wave. Nevertheless, there was no prominent consistent change in gene expression linked to heat stress. Considering all of this together, we conclude that several frequencies and power densities are capable of promoting non-thermal changes in the gene expression of different cell types and in regulating pathways.

4.3. Exposure Parameters

Two main factors are to be considered when assessing the biological effects of electromagnetic THz radiation; one is the irradiation parameters and the second is the nature of the biological material. Although there are some consistent and highly mixed patterns of effects shown by the authors, the diversity in exposure parameters and biological models made it hard to identify and single out one or several biomarkers. In light of this, it can be said that the degree of THz absorption and the biological response to the stress with regard to THz exposure depends on various factors, including the frequency, power density, and exposure time, as well as the biological system properties and cell types, which implies the heterogeneity in differential responses. Moreover, the biological impacts could be attributed to modifications in the cell molecular structure in THz-exposed samples [44,45].
These observations raise two problems of great importance in the field of terahertz exposure: the reproducibility and replicability of the results and the impact of the biological model. The lack of uniformity among protocols renders increasingly difficult challenges to replicating studies individually. The findings of terahertz experiments are difficult to explain and compare because of the variety of models, frequencies, exposure systems, power densities, and exposure durations. Additionally, attention is often drawn to the variations in frequencies and exposure conditions to compare differences between THz studies, and as a consequence, the influence of the biological material is, in many instances, minimized.
Moreover, it is hypothesized that exposure in terms of incident power density on a culture medium has a limited impact on the cells, which may have been far away from where the THz energy is absorbed. Or, in the case where it reaches the cells, it might not be ubiquitously distributed. Thus, better methods of irradiation techniques are necessary for high-quality THz biological effects studies.

4.4. Gene Expression Technology

Microarray and recently RNA-Seq have become powerful tools to solve a variety of biological problems in study designs such as exposed versus unexposed groups, with the biological goal of detecting possible statistically differential expression [46,47].
RNA-seq appears to perform better in gene expression profiling compared to microarray. Log2 fold change (FC) and false discovery rate (FDR) are often used to screen differentially expressed genes (DEGs). DEGs results tables consist of all genes that have an FDR of less than 0.05 as hard thresholds and are statistically associated with a particular condition. It is noteworthy that a decrease in FC results in an increase in the number of DEGs [48]. Differences in thresholding applied to fold change will give different sets of DEGs as a consequence of different processes, signaling pathways, and functions involved as a result of biological experiments [49]. The assumption is that genes that react the most have the highest fold change; therefore, expression changes by a factor of 1.5 or 2 are much more interesting than genes that do not really change. The above methods have drawbacks. Trivial phenotypical changes can be caused by a trivial variation in expression. Variability in FC is beneficial in reducing or increasing gene expression data, the power of analysis, and background noise. Fold change increase leads towards a decrease in the number of DEGs. Genes can be significantly expressed with low fold changes. The mathematical and biological meaning resulting from global differential gene expression is not always the same. It is good to identify DEGs but it is also crucial to put it into a biological framework. The biological context of fold change is likely linked to both the genes and the experimental context [50].
In addition, pathways analyses are derived only from about 1% of the original data collected versus 99% discarded with their biological meanings. An increasing number of researchers apply gene ontology or pathway enrichment analysis to better fit the study purpose. Gene expression and transcriptomic technologies produce a tremendous amount of quantitative data, usually an extended list of transcripts that are overwhelming to comprehend. Roughly 150 to 600 genes are hard to decipher and overwhelming. In addition, some genes are promoters and transcription factors that control the activity of other interesting sets of genes. Genes tend to be team players in networks in order to exhibit a distinct biological function. A way to deal with this issue is the use of bioinformatics tools to establish and investigate which of the genes share a pathway or network and, therefore, work towards a unique function. Pathway analysis helps to derive meaning from these high-throughput experiments with different features, validate or confirm expected results, generate new hypotheses, and design subsequent experiments.
GO (Gene Ontology) [51], KEGG (Kyoto Encyclopedia for Genes and Genomes) [52], IPA (Ingenuity Pathway Analysis) [53], and Reactome [54] provide pathways analysis approaches such as over-representation analysis (ORA) and GSEA (Gene Set Enrichment Analysis). One problem encountered with ORA and GSEA is not only the major differences between these curated databases in their terminology but also the differences in the data underlying pathway definitions. Each pathway is defined by a set of genes that may not exactly be the same and render pathway descriptions slightly differently across databases. Terms may overlap across databases or be unique to one database. Whereas GO is more specialized towards cellular function, KEGG is oriented towards metabolism pathways, while Reactome has a lot more disease-related pathways than GO and KEGG.
Furthermore, in gene interactions, the promoter directions are not considered, thereby producing false results in certain conditions. Moreover, gene expression studies using an observational design between non-exposed and exposed groups only provide a brief picture or snapshot of the biological process at a given time.
Time points might not be optimal for capturing gene expression alterations in the biological model used. The limited range of time points in RNA collection produces a weaker response or may provide false-negative results [55]. Longer studies and follow-ups may permit a better dissection of the biological variations. While it is impossible to determine the extent of which the platform type may have influenced the results, lots of the variability reported in the studies might be a direct result of the usage of different gene expression assessment technologies and pipelines for biological quantification and the applied statistical methods.

4.5. Limitations

One of the strengths of this study resides in its wide eligibility criteria, providing considerable heterogeneity among the included studies and an exhaustive search of multiple databases. We employed PRISMA guidelines and searched multiple databases. As in most systematic reviews, our review is limited by the heterogeneity of the outcome definitions. The studies differed with regard to parameters that they have used. A second limitation of this review is that the heterogeneity of the studies did not allow a strong meta-analysis to offer a larger picture of the differentially expressed gene pathways due to differences in the study methodology. Despite these limitations, we were still able to synthesize the data by listing the reported gene expression patterns. An additional limitation was excluding studies published not in English.

5. Conclusions

To conclude, this systematic review is the first that has investigated terahertz radiation studies using microarray or RNA-Seq. We undertook a review to assess the degree to which studies that had focused on gene expression led to overall conclusions. The review also highlighted the methodological variations across studies. A standardization of the methods for reporting results could enhance our understanding of THz biological effects and indicate which family genes are the most likely to exert their influence on THz stress response. There are several limitations of the research that need to be improved in coming studies, among which are better-designed longitudinal studies, mainly with larger samples and more than one exploratory method for more accurate observations. Our investigation consolidates and advances our understanding of the gene expression underlying THz exposure.
Finally, in the future, additional and more comprehensive studies should focus on the epigenetic underpinnings of the range of responses or evidence. Moreover, current breakthroughs in the field of bioinformatics such as metabolomics profiling and array-based methylation profiling could fill the voids and provide a more complete picture of the complex connection between THz and gene expression.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15081045/s1, Table S1. Gene set enrichment analysis (C5 Gene Ontology) of differentially expressed genes (p.adjust < 0.05).

Author Contributions

Conceptualization, M.N.D.; methodology, M.N.D.; formal analysis, M.N.D.; data curation, M.N.D.; writing—original draft preparation, M.N.D.; writing—review and editing; M.N.D., S.S., Q.Z. and S.Z.; supervision, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China and China Academy of Engineering Physics (U2330108).

Institutional Review Board Statement

This research was registered with Prospero (Registration No: CDR42024502937).

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Nagatsuma, T. Terahertz technologies: Present and future. IEICE Electron. Express 2011, 8, 1127–1142. [Google Scholar] [CrossRef]
  2. Sun, L.; Zhao, L.; Peng, R.-Y. Research progress in the effects of terahertz waves on biomacromolecules. Mil. Med. Res. 2021, 8, 28. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, M.; Liu, J.; Liang, W.; Lu, B.; Fan, P.; Song, Y.; Wang, M.; Wu, Y.; Cai, X. Recent advances and research progress on microsystems and bioeffects of terahertz neuromodulation. Microsyst. Nanoeng. 2023, 9, 143. [Google Scholar] [CrossRef] [PubMed]
  4. Romanenko, S.; Begley, R.; Harvey, A.R.; Hool, L.; Wallace, V.P. The interaction between electromagnetic fields at megahertz, gigahertz and terahertz frequencies with cells, tissues and organisms: Risks and potential. J. R. Soc. Interface 2017, 14, 20170585. [Google Scholar] [CrossRef] [PubMed]
  5. Hintzsche, H.; Jastrow, C.; Kleine-Ostmann, T.; Kärst, U.; Schrader, T.; Stopper, H. Terahertz electromagnetic fields (0.106 THz) do not induce manifest genomic damage in vitro. PLoS ONE 2012, 7, e46397. [Google Scholar] [CrossRef] [PubMed]
  6. Cherkasova, O.P.; Serdyukov, D.S.; Nemova, E.F.; Ratushnyak, A.S.; Kucheryavenko, A.S.; Dolganova, I.N.; Xu, G.; Skorobogatiy, M.; Reshetov, I.V.; Timashev, P.S.; et al. Cellular effects of terahertz waves. J. Biomed. Opt. 2021, 26, 090902. [Google Scholar] [CrossRef] [PubMed]
  7. Dhillon, S.S.; Vitiello, M.S.; Linfield, E.H.; Davies, A.G.; Hoffmann, M.C.; Booske, J.; Paoloni, C.; Gensch, M.; Weightman, P.; Williams, G.P.; et al. The 2017 terahertz science and technology roadmap. J. Phys. Appl. Phys. 2017, 50, 043001. [Google Scholar] [CrossRef]
  8. Kim, K.-T.; Park, J.; Jo, S.J.; Jung, S.; Kwon, O.S.; Gallerano, G.P.; Park, W.-Y.; Park, G.-S. High-power femtosecond-terahertz pulse induces a wound response in mouse skin. Sci. Rep. 2013, 3, 2296. [Google Scholar] [CrossRef]
  9. Ayyalasomayajula, S.; Willows, S.; Hough, C.; Hegmann, F.; Kulka, M. Terahertz radiation may not alter rat basophilic leukemia cell membrane permeability to propidium iodide. Spectrum 2021, 1–11. [Google Scholar] [CrossRef]
  10. Bock, J.; Fukuyo, Y.; Kang, S.; Phipps, M.L.; Alexandrov, L.B.; Rasmussen, K.Ø.; Bishop, A.R.; Rosen, E.D.; Martinez, J.S.; Chen, H.-T.; et al. Mammalian Stem Cells Reprogramming in Response to Terahertz Radiation. PLoS ONE 2010, 5, e15806. [Google Scholar] [CrossRef]
  11. Zeni, O.; Gallerano, G.P.; Perrotta, A.; Romanò, M.; Sannino, A.; Sarti, M.; D’Arienzo, M.; Doria, A.; Giovenale, E.; Lai, A.; et al. Cytogenetic Observations in Human Peripheral Blood Leukocytes following In Vitro Exposure to Thz Radiation: A Pilot Study. Health Phys. 2007, 92, 349. [Google Scholar] [CrossRef] [PubMed]
  12. Korenstein-Ilan, A.; Barbul, A.; Hasin, P.; Eliran, A.; Gover, A.; Korenstein, R. Terahertz Radiation Increases Genomic Instability in Human Lymphocytes. Radiat. Res. 2008, 170, 224–234. [Google Scholar] [CrossRef] [PubMed]
  13. Li, K.; Sasaki, K.; Watanabe, S.; Shirai, H. Relationship between power density and surface temperature elevation for human skin exposure to electromagnetic waves with oblique incidence angle from 6 GHz to 1 THz. Phys. Med. Biol. 2019, 64, 065016. [Google Scholar] [CrossRef] [PubMed]
  14. Wilmink, G.J.; Grundt, J.E. Invited Review Article: Current State of Research on Biological Effects of Terahertz Radiation. J. Infrared Millim. Terahertz Waves 2011, 32, 1074–1122. [Google Scholar] [CrossRef]
  15. Dalzell, D.R.; McQuade, J.; Vincelette, R.; Ibey, B.; Payne, J.; Thomas, R.; Roach, W.P.; Roth, C.L.; Wilmink, G.J. Damage thresholds for terahertz radiation. In Proceedings of the Optical Interactions with Tissues and Cells XXI, San Francisco, CA, USA, 25–27 January 2010; Jansen, E.D., Thomas, R.J., Eds.; SPIE: Bellingham, WA, USA, 2010; p. 75620M. [Google Scholar] [CrossRef]
  16. Alexandrov, B.S.; Phipps, M.L.; Alexandrov, L.B.; Booshehri, L.G.; Erat, A.; Zabolotny, J.; Mielke, C.H.; Chen, H.-T.; Rodriguez, G.; Rasmussen, K.Ø.; et al. Specificity and heterogeneity of terahertz radiation effect on gene expression in mouse mesenchymal stem cells. Sci. Rep. 2013, 3, 1184. [Google Scholar] [CrossRef] [PubMed]
  17. Havas, M. When theory and observation collide: Can non-ionizing radiation cause cancer? Environ. Pollut. 2017, 221, 501–505. [Google Scholar] [CrossRef] [PubMed]
  18. Amicis, A.D.; Sanctis, S.D.; Cristofaro, S.D.; Franchini, V.; Lista, F.; Regalbuto, E.; Giovenale, E.; Gallerano, G.P.; Nenzi, P.; Bei, R.; et al. Biological effects of in vitro THz radiation exposure in human foetal fibroblasts. Mutat. Res. Toxicol. Environ. Mutagen. 2015, 793, 150–160. [Google Scholar] [CrossRef] [PubMed]
  19. Sitnikov, D.S.; Ilina, I.V.; Revkova, V.A.; Rodionov, S.A.; Gurova, S.A.; Shatalova, R.O.; Kovalev, A.V.; Ovchinnikov, A.V.; Chefonov, O.V.; Konoplyannikov, M.A.; et al. Effects of high intensity non-ionizing terahertz radiation on human skin fibroblasts. Biomed. Opt. Express 2021, 12, 7122–7138. [Google Scholar] [CrossRef] [PubMed]
  20. Amini, T.; Jahangiri, F.; Ameri, Z.; Hemmatian, M.A. A Review of Feasible Applications of THz Waves in Medical Diagnostics and Treatments. J. Lasers Med. Sci. 2021, 12, e92. [Google Scholar] [CrossRef]
  21. Titova, L.V.; Ayesheshim, A.K.; Golubov, A.; Fogen, D.; Rodriguez-Juarez, R.; Hegmann, F.A.; Kovalchuk, O. Intense THz pulses cause H2AX phosphorylation and activate DNA damage response in human skin tissue. Biomed. Opt. Express 2013, 4, 559. [Google Scholar] [CrossRef]
  22. Xu, D.; Liu, R.; Li, B.; Zhang, H.; Yang, Y. Effect and mechanism of terahertz irradiation in repairing spinal cord injury in mice. Gene 2023, 860, 147218. [Google Scholar] [CrossRef]
  23. Page, M.J.; McKenzie, J.E.; Bossuyt, P.M.; Boutron, I.; Hoffmann, T.C.; Mulrow, C.D.; Shamseer, L.; Tetzlaff, J.M.; Akl, E.A.; Brennan, S.E.; et al. The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. BMJ 2021, 372, n71. [Google Scholar] [CrossRef]
  24. Edgar, R. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30, 207–210. [Google Scholar] [CrossRef] [PubMed]
  25. Davis, S.; Meltzer, P.S. GEOquery: A bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics 2007, 23, 1846–1847. [Google Scholar] [CrossRef]
  26. Gentleman, R.C.; Carey, V.J.; Bates, D.M.; Bolstad, B.; Dettling, M.; Dudoit, S.; Ellis, B.; Gautier, L.; Ge, Y.; Gentry, J.; et al. Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004, 5, R80. [Google Scholar] [CrossRef]
  27. McCall, M.N.; Jaffee, H.A.; Irizarry, R.A. fRMA ST: Frozen robust multiarray analysis for Affymetrix Exon and Gene ST arrays. Bioinformatics 2012, 28, 3153–3154. [Google Scholar] [CrossRef]
  28. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  29. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. OMICS J. Integr. Biol. 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  30. Liberzon, A.; Birger, C.; Thorvaldsdóttir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef]
  31. Shang, S.; Wu, X.; Zhang, Q.; Zhao, J.; Hu, E.; Wang, L.; Lu, X. 0.1 THz exposure affects primary hippocampus neuron gene expression via alternating transcription factor binding. Biomed. Opt. Express 2021, 12, 3729. [Google Scholar] [CrossRef]
  32. Titova, L.V.; Ayesheshim, A.K.; Golubov, A.; Rodriguez-Juarez, R.; Woycicki, R.; Hegmann, F.A.; Kovalchuk, O. Intense THz pulses down-regulate genes associated with skin cancer and psoriasis: A new therapeutic avenue? Sci. Rep. 2013, 3, 2363. [Google Scholar] [CrossRef] [PubMed]
  33. Peltek, S.; Meshcheryakova, I.; Kiseleva, E.; Oshchepkov, D.; Rozanov, A.; Serdyukov, D.; Demidov, E.; Vasiliev, G.; Vinokurov, N.; Bryanskaya, A.; et al. E. coli aggregation and impaired cell division after terahertz irradiation. Sci. Rep. 2021, 11, 20464. [Google Scholar] [CrossRef] [PubMed]
  34. Cheon, H.; Hur, J.K.; Hwang, W.; Yang, H.-J.; Son, J.-H. Epigenetic modification of gene expression in cancer cells by terahertz demethylation. Sci. Rep. 2023, 13, 4930. [Google Scholar] [CrossRef]
  35. Bogomazova, A.N.; Vassina, E.M.; Goryachkovskaya, T.N.; Popik, V.M.; Sokolov, A.S.; Kolchanov, N.A.; Lagarkova, M.A.; Kiselev, S.L.; Peltek, S.E. No DNA damage response and negligible genome-wide transcriptional changes in human embryonic stem cells exposed to terahertz radiation. Sci. Rep. 2015, 5, 7749. [Google Scholar] [CrossRef] [PubMed]
  36. Echchgadda, I.; Grundt, J.E.; Cerna, C.Z.; Roth, C.C.; Payne, J.A.; Ibey, B.L.; Wilmink, G.J. Terahertz Radiation: A Non-contact Tool for the Selective Stimulation of Biological Responses in Human Cells. IEEE Trans. Terahertz Sci. Technol. 2016, 6, 54–68. [Google Scholar] [CrossRef]
  37. Hough, C.M.; Purschke, D.N.; Huang, C.; Titova, L.V.; Kovalchuk, O.; Warkentin, B.J.; Hegmann, F.A. Topology-Based Prediction of Pathway Dysregulation Induced by Intense Terahertz Pulses in Human Skin Tissue Models. J. Infrared Millim. Terahertz Waves 2018, 39, 887–898. [Google Scholar] [CrossRef]
  38. Hough, C.M.; Purschke, D.N.; Huang, C.; Titova, L.V.; Kovalchuk, O.V.; Warkentin, B.J.; Hegmann, F.A. Intense terahertz pulses inhibit Ras signaling and other cancer-associated signaling pathways in human skin tissue models. J. Phys. Photonics 2021, 3, 034004. [Google Scholar] [CrossRef]
  39. Zhao, J.-W.; He, M.-X.; Dong, L.-J.; Li, S.-X.; Liu, L.-Y.; Bu, S.-C.; Ouyang, C.-M.; Wang, P.-F.; Sun, L.-L. Effect of terahertz pulse on gene expression in human eye cells. Chin. Phys. B 2019, 28, 048703. [Google Scholar] [CrossRef]
  40. Zhao, X.; Zhang, M.; Liu, Y.; Liu, H.; Ren, K.; Xue, Q.; Zhang, H.; Zhi, N.; Wang, W.; Wu, S. Terahertz exposure enhances neuronal synaptic transmission and oligodendrocyte differentiation in vitro. iScience 2021, 24, 103485. [Google Scholar] [CrossRef]
  41. Bucur-Portase, R.-C. Introduction to the Biological Effects of Terahertz Radiation. In Trends in Terahertz Technology; Meziani, Y.M., Velázquez-Pérez, J.E., Eds.; IntechOpen: London, UK, 2023. [Google Scholar] [CrossRef]
  42. Zakrzewski, W.; Dobrzyński, M.; Szymonowicz, M.; Rybak, Z. Stem cells: Past, present, and future. Stem Cell Res. Ther. 2019, 10, 68. [Google Scholar] [CrossRef]
  43. Alexandrov, B.S.; Rasmussen, K.Ø.; Bishop, A.R.; Usheva, A.; Alexandrov, L.B.; Chong, S.; Dagon, Y.; Booshehri, L.G.; Mielke, C.H.; Phipps, M.L.; et al. Non-thermal effects of terahertz radiation on gene expression in mouse stem cells. Biomed. Opt. Express 2011, 2, 2679. [Google Scholar] [CrossRef] [PubMed]
  44. Nikitkina, A.I.; Bikmulina, P.Y.; Gafarova, E.R.; Kosheleva, N.V.; Efremov, Y.M.; Bezrukov, E.A.; Butnaru, D.V.; Dolganova, I.N.; Chernomyrdin, N.V.; Cherkasova, O.P.; et al. Terahertz radiation and the skin: A review. J. Biomed. Opt. 2021, 26, 043005. [Google Scholar] [CrossRef] [PubMed]
  45. Vilagosh, Z.; Lajevardipour, A.; Wood, A.W. Computational phantom study of frozen melanoma imaging at 0.45 terahertz. Bioelectromagnetics 2019, 40, 118–127. [Google Scholar] [CrossRef] [PubMed]
  46. Corchete, L.A.; Rojas, E.A.; Alonso-López, D.; De Las Rivas, J.; Gutiérrez, N.C.; Burguillo, F.J. Systematic comparison and assessment of RNA-seq procedures for gene expression quantitative analysis. Sci. Rep. 2020, 10, 19737. [Google Scholar] [CrossRef]
  47. Mathur, R.; Rotroff, D.; Ma, J.; Shojaie, A.; Motsinger-Reif, A. Gene set analysis methods: A systematic comparison. BioData Min. 2018, 11, 8. [Google Scholar] [CrossRef] [PubMed]
  48. Zhao, S.; Fung-Leung, W.-P.; Bittner, A.; Ngo, K.; Liu, X. Comparison of RNA-Seq and microarray in transcriptome profiling of activated T cells. PLoS ONE 2014, 9, e78644. [Google Scholar] [CrossRef]
  49. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef]
  50. McCarthy, D.J.; Smyth, G.K. Testing significance relative to a fold-change threshold is a TREAT. Bioinformatics 2009, 25, 765–771. [Google Scholar] [CrossRef] [PubMed]
  51. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef]
  52. Kanehisa, M.; Goto, S.; Sato, Y.; Kawashima, M.; Furumichi, M.; Tanabe, M. Data, information, knowledge and principle: Back to metabolism in KEGG. Nucleic Acids Res. 2014, 42, D199–D205. [Google Scholar] [CrossRef]
  53. Knowledge and Databases, Bioinforma. Softw. QIAGEN Digit. Insights (n.d.). Available online: https://digitalinsights.qiagen.com/knowledge-and-databases/ (accessed on 14 June 2024).
  54. Jassal, B.; Matthews, L.; Viteri, G.; Gong, C.; Lorente, P.; Fabregat, A.; Sidiropoulos, K.; Cook, J.; Gillespie, M.; Haw, R.; et al. The reactome pathway knowledgebase. Nucleic Acids Res. 2019, 48, D498–D503. [Google Scholar] [CrossRef] [PubMed]
  55. Leszczynski, D.; De Pomerai, D.; Koczan, D.; Stoll, D.; Franke, H.; Albar, J.P. Five years later: The current status of the use of proteomics and transcriptomics in EMF research. Proteomics 2012, 12, 2493–2509. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study selection following Prisma flowchart.
Figure 1. Study selection following Prisma flowchart.
Genes 15 01045 g001
Figure 2. UpSet plot visualizing the distribution and overlap between sets of unique or shared differentially expressed genes and their intersection among microarray datasets (Alexandrov et al. [16]; Kim et al. [8]; Zhao et al. [40]). The horizontal bars depict the number of DEGs. The vertical bars display common elements in the sets, represented with dots under each bar.
Figure 2. UpSet plot visualizing the distribution and overlap between sets of unique or shared differentially expressed genes and their intersection among microarray datasets (Alexandrov et al. [16]; Kim et al. [8]; Zhao et al. [40]). The horizontal bars depict the number of DEGs. The vertical bars display common elements in the sets, represented with dots under each bar.
Genes 15 01045 g002
Figure 3. Enrichment analysis: (a) significant genes enriched from MsigDB (c5 canonical pathways); (b) significant genes in set/total genes in set (nk/k); (c) Cnet plot of significant genes enriched.
Figure 3. Enrichment analysis: (a) significant genes enriched from MsigDB (c5 canonical pathways); (b) significant genes in set/total genes in set (nk/k); (c) Cnet plot of significant genes enriched.
Genes 15 01045 g003aGenes 15 01045 g003b
Table 2. Gene expression profile approach and study outcomes.
Table 2. Gene expression profile approach and study outcomes.
StudyGene Expression Profile MethodLog 2 Fold ChangeNumbers of DEGsPathwaysBiological Outcomes
Bock et al. [10]MicroarrayNA2204NAGene related to the following: Adiponectin ↑; GLUT 4 ↑; FABP4 ↑; PPARG ↑;
PPARA ↓; LXR ↓; Cyclin B ↓; TBP ↓
Alexandrov et al. [16]Microarray020GOGene related to the following: Adiponectin ↑; FABP4 ↑; PPARG ↑;
Nfe2l2 ↓; Gem-; Slco4a1 ↓
Gem ↓; Slco4a1 ↓; Arg2 ↓; Ifitm6 ↓; Hspa1a ↓; Hspa1b ↓; Hipk1 ↓; Slc6a12 ↓; 4921509J17Rik ↓; Hspa4l ↓; Rhpn2 ↓; Akr1b3 ↓; Flrt3 ↓; Tet2+; Ctgf ↑; Skil ↓; Fgd3 ↓; Nfe2l2 ↑; Dclre1b ↓; Fn1↓
Kim et al. [8]Microarray1.5149IPAGene related to the following: inflammation (TGF- β ↑); wound healing response (Bmp2 ↑, Cd44 ↑, Thbs1 ↑, Serpine1 ↑, Krt6a ↑, Lep ↑, Sprrb1 ↑)
Tivota et al. [32]Microarray0.6442 Gene related to the following: psoriasis (S100 family: S100A11 ↓, S100A15 ↓, S100A12 ↓; Small Proline Rich Region: SPRR1B ↓, SPRR2A ↓, SPRR3 ↓, SPRR2B ↓, SPRR2C ↓, LCE3D ↓, involucrin (IVL) ↓; Serine protease inhibitor: SerpinB3 ↓, SerpinB4 ↓, SerpinB7 ↓, SerpinB13 ↓); atopic dermatitis ↓; epidermal hyperplasia ↓; dermatitis ↓; cancer (DEFB103A ↓, DEFB4 ↓, DEFB1 ↓, LOC728454 ↓, CD24 ↓); apoptotic signaling pathways (S100A ↓, S100P ↓)
Echchgadda et al. [36]Microarray2531IPAGene related to the following: cytoplasm (ARSB ↓, CYP19A1 ↓, CYP2D6+, CYP2E1 ↓, CYP4X1 ↑, GOT1L1 ↑, MAOB ↑, PCYT1B ↑, PLD1 ↑, TAT ↑, ATP2A1 ↑, BHMT ↑, CAD ↓, CTPS2 ↑, SARDH ↑, ACTA1 ↑, ACTA2 ↑, MYH4 ↓, RHOB ↑, RHOJ ↑, WASL ↑, STON2 ↑, MAP2K5 ↑, PDE1A ↑, PIK3C2A ↑, PLCB4 ↓, PLD1 ↑, RAB11B ↑, RAB7A ↑, RGS7 ↓, SOS2 ↓); plasma membrane (CD44 ↓, ATP2B3 ↑, CDH19 ↑, DSP ↑, ITGB1 ↑, ITGB3 ↓, ADRA1D ↑, DRD1 ↓, FGFR3 ↑, GABBR2 ↑, GNA11 ↑, NPR3 ↑, NTRK2 ↑, PTGER3 ↑,); cell shape, adhesion and cytoskeleton organization; immune response (HLA-DPA1 ↓, HLA-DQA1 ↑, IGHG1 ↑, IL4 ↑); transcription (ATF2 ↑, APBB1 ↑, JUN ↑, TNL2 ↑, RGS12 ↑); extracellular space (CANT1 ↓, FGF10 ↓, FGF7 ↓, PLAU ↑, ALB ↓, APOE ↑)
Echchgadda et al. [36]Microarray2531IPAGene related to the following: cytoplasm (ARSB ↓, CYP19A1 ↓, CYP2D6+, CYP2E1 ↓, CYP4X1 ↑, GOT1L1 ↑, MAOB ↑, PCYT1B ↑, PLD1 ↑, TAT ↑, ATP2A1 ↑, BHMT ↑, CAD ↓, CTPS2 ↑, SARDH ↑, ACTA1 ↑, ACTA2 ↑, MYH4 ↓, RHOB ↑, RHOJ ↑, WASL ↑, STON2 ↑, MAP2K5 ↑, PDE1A ↑, PIK3C2A ↑, PLCB4 ↓, PLD1 ↑, RAB11B ↑, RAB7A ↑, RGS7 ↓, SOS2 ↓); plasma membrane (CD44 ↓, ATP2B3 ↑, CDH19 ↑, DSP ↑, ITGB1 ↑, ITGB3 ↓, ADRA1D ↑, DRD1 ↓, FGFR3 ↑, GABBR2 ↑, GNA11 ↑, NPR3 ↑, NTRK2 ↑, PTGER3 ↑,); cell shape, adhesion and cytoskeleton organization; immune response (HLA-DPA1 ↓, HLA-DQA1 ↑, IGHG1 ↑, IL4 ↑); transcription (ATF2 ↑, APBB1 ↑, JUN ↑, TNL2 ↑, RGS12 ↑); extracellular space (CANT1 ↓, FGF10 ↓, FGF7 ↓, PLAU ↑, ALB ↓, APOE ↑)
Bogomazova et al. [35]Microarray1.573GOGene related to the following: mitochondrial ribosome (MRPL34 ↑, MRPL43 ↑, MRPL55 ↑, MRPS24 ↑)
Hough et al. [37,38]Microarray1.51681GO and KEGGGene related to the following: Ras–Raf–Mek–Erk Cascade (KRAS ↓, RAS ↓, MAPK3 ↓); apoptosis (AKT ↑); pro-inflammatory AND cytokine–cytokine receptor interaction (CXCL16 ↓, BMP7 ↓, CCL8 ↑, CXCL5 ↑, CCL20 ↑, Il-6 ↑, IL-24 ↑); glioma pathway (CXCL16 ↓, BMP7 ↓, CCL8 ↑, CXCL5 ↑, CCL20 ↑, Il-6 ↑, IL-24 ↑, HGF ↑, RAP1B ↑, FGF2 ↑, FGF7 ↑, RAC2 ↓, RRAS ↓, EFNA1 ↓, PGF ↓, calmodulin family: CALML5 ↓, CALML3 ↓, CALML1 ↓, CAMK families ↓), calmodulin (CaM ↓), calmodulin-like (CaLM ↓), and calmodulin-dependent kinase (CaMK) families, GNAI3 ↓, EDNRB ↑, PPP3CC ↑.
Zhao et al. [39]RNA-Seq1.5HCE+
Müller+
RPE = 766
GOGene related to the following: cell growth, differentiation and migration (RET ↑, ALK ↑, ROS1 ↑); fibroblast growth factor (FGFR2 ↑); stable spindle microtubule (CENPE ↓); MUC16 ↓; COL1A2 ↑
Peltek et al. [33]RNA-Seq2 (up)
1 (down)
741KEGGGene related to the following: cell aggregation (tdcABCDEFGR, and matA-F genes); suppression of cell motility (gene yjjQ); suppression of cell division (dicABCF, FtsZ, and minCDE); adhesin synthesis (sfmACDHF); cell envelope stabilization (yjbEFGH, gfcA)|pilus ↑; molecular functions of fimbrial porins ↑; organization and assembli of pili ↑; cell adhesion ↑; cytoplasmic genes ↓; genes of the respiratory chain ↓; ribosomes ↓; processes of aerobic respiration ↓; translation ↓
Shang et al. [31]RNA-Seq1165GO and KEGGGene related to the following: ribosome (small subunit: Rps18, Rps28; large subunit: Rp136, Rp136A, Rp135; mitoribosomal: Mrps18c, Mrps15, Mrps11; Mrp114; Mrp118, Mrp158, Mrp121) ↑; calcium signaling pathway; axon guidance; insulin signaling pathway; GTPase
MF: Rab GTPase binding; guanyl-nucleotide exchange factor activity; GTPase binding; small GTPase binding; Ras GTPase binding
Aquaporin (Aqp5) ↑; matrix metallopeptidase (Mmp3) ↑; transcriptional factor (Erf) ↓; alipoprotein (Apof) ↓; plasma membrane calcium ATPase (Atp2b4) ↓; intracellular receptor for inositol (Itpr1)↓
Zhao et al. [40]Microarray1.51249GO and KEGGGene related to the following: neuron projection ↑; synapse organization ↑; dendritic spine ↑; actin binding ↓; cytoskeleton organization ↓; response to stimulus ↓
Ppp1r2 ↑, Camk2b ↑, Fmr1 ↑, Nlgn1 ↑, Nrxn1 ↑, Rapgef ↑, Brwd3 ↓, Crx ↓, Olfr18 ↓, Svil ↓, Tas2r1 ↓, Tnf ↓
Cheon et al. [34]Microarray2 KEGGGene related to the following: cancer, apoptosis pathways, and inflammation related-gene (transcription factor for AP-1: JUN, FOS, CXCL8 ↓)
Xu et al. [22]RNA-Seq1137GO and KEGGGene related to the following: inflammation; immune regulation; IL-17 signaling pathway (alipoprotein C-I ↑; alipoprotein C-IV ↑; immunoglobulin lambda constant 1 ↑; immunoglobulin heavy constant γ 2B ↑)
Notes. ↑ upregulation; ↓ downregulation.
Table 3. Dataset collected and included in the meta-analysis.
Table 3. Dataset collected and included in the meta-analysis.
DatasetAuthorYearGSE IDPlatformCell TypeTotal Sample
Dataset_1Alexandrov et al. [16]2013GSE41083GPL1261 [Mouse430_2] AffymetrixMouse stem cells6
Dataset_2Alexandrov et al. [16]2013GSE41084GPL1261 [Mouse430_2] AffymetrixMouse stem cells6
Dataset_3Alexandrov et al. [16]2013GSE41085GPL1261 [Mouse430_2] AffymetrixMouse stem cells6
Dataset_4Kim et al. [8]2013GSE44671GPL6246 [MoGene-1_0-st] AffymetrixC57BL/6J and BALB/c skin cell6
Dataset_5Zhao et al. [40]2021GSE178729GPL21163 Array Agilent-074809Mouse cortical neurons and oligodendrocytes precursor cells6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dione, M.N.; Shang, S.; Zhang, Q.; Zhao, S.; Lu, X. Non-Thermal Effects of Terahertz Radiation on Gene Expression: Systematic Review and Meta-Analysis. Genes 2024, 15, 1045. https://doi.org/10.3390/genes15081045

AMA Style

Dione MN, Shang S, Zhang Q, Zhao S, Lu X. Non-Thermal Effects of Terahertz Radiation on Gene Expression: Systematic Review and Meta-Analysis. Genes. 2024; 15(8):1045. https://doi.org/10.3390/genes15081045

Chicago/Turabian Style

Dione, Mactar Ndiaga, Sen Shang, Qi Zhang, Sicheng Zhao, and Xiaoyun Lu. 2024. "Non-Thermal Effects of Terahertz Radiation on Gene Expression: Systematic Review and Meta-Analysis" Genes 15, no. 8: 1045. https://doi.org/10.3390/genes15081045

APA Style

Dione, M. N., Shang, S., Zhang, Q., Zhao, S., & Lu, X. (2024). Non-Thermal Effects of Terahertz Radiation on Gene Expression: Systematic Review and Meta-Analysis. Genes, 15(8), 1045. https://doi.org/10.3390/genes15081045

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop