A. N. Akimniyazova1
T. K. Niyazova2O. Yu. Yurikova2A. Yu. Pyrkova2,3M. A. Zhanuzakov1
A. T. Ivashchenko3*- 1Higher School of Medicine, Faculty of Medicine and Healthcare, Al-Farabi Kazakh National University, Almaty, Kazakhstan
- 2Department of Biotechnology, Faculty of Biology and Biotechnology, Al-Farabi Kazakh National University, Almaty, Kazakhstan
- 3Center for Bioinformatics and Nanomedicine, Almaty, Kazakhstan
Elucidation of ways to regulate the expression of candidate cancer genes will contribute to the development of methods for cancer diagnosis and therapy. The aim of the present study was to show the role of piRNAs as efficient regulators of mRNA translation of esophageal adenocarcinoma (EAC) candidate genes. We used bioinformatic methods to study the interaction characteristics of up to 200 thousand piRNAs with mRNAs of 38 candidate EAC genes. The piRNAs capable of binding to mRNA of AR, BTG3, CD55, ERBB3, FKBP5, FOXP1, LEP, SEPP1, SMAD4, and TP53 genes with high free energy by the formation of hydrogen bonds between canonical (G-C, A-U) and noncanonical (G-U, A-C) piRNA and mRNA nucleotide pairs were revealed. The organization of piRNA binding sites (BSs) in the mRNA of candidate genes was found to overlap nucleotide sequences to form clusters. Clusters of piRNA BSs were detected in the 5′-untranslated region, coding domain sequence, and 3′-untranslated region of mRNA. Due to the formation of piRNA binding site clusters, compaction of BSs occurs and competition between piRNAs for binding to mRNA of candidate EAC genes occurs. Associations of piRNA and candidate genes were selected for use as markers for the diagnosis of EAC.
Introduction
The biological role of piRNA (PIWI-interacting RNA), despite the long-standing discovery of piRNA (Kim, 2006; Aravin et al., 2007), remains at the stage of speculation. There are no works on the direct interaction of piRNAs with mRNAs of genes. In recent years, there have been reports of the involvement of piRNA in the regulation of many biological processes (Czech et al., 2018; Chen et al., 2021; Zhang et al., 2022a). Since the present work is devoted to the involvement of piRNA in carcinogenesis, here we will review the existing ideas about the role of piRNA in the development of various forms of cancer. It has been shown that piRNAs are abnormally expressed in various forms of cancer and the authors suggest using piRNA in diagnosis and therapy (Liu et al., 2019). In recent years, the role of piRNA in the development of breast cancer has been actively studied. The molecular mechanisms and possible clinical implications of PIWI-interacting RNA have been investigated in this disease (Tan et al., 2019; Qian et al., 2021; Yuan et al., 2021). PIWI-interacting RNAs have been shown to be involved in proliferation, migration, and apoptosis in breast cancer (Kärkkäinen et al., 2021; Liu et al., 2021). The role of piRNA in the pathogenesis and diagnosis of esophageal cancer, colon and stomach cancer, colorectal cancer has been considered (Lin et al., 2019; Halajzadeh et al., 2020; Ameli Mojarad et al., 2022; Jia et al., 2022). The biogenesis and role of piRNA in several cancer localizations at once have been studied (Guo et al., 2020). A number of studies suggest that piRNAs are involved in oncogenesis through their influence on genes, but how this occurs remains unknown (Fonseca et al., 2020; Wang et al., 2022a; Zhang et al., 2022a; Zhang et al., 2022b; Cai et al., 2022; Zhang et al., 2022c; Hanusek et al., 2022). In recent years, the study of piRNAs has continued to receive a great deal of attention. For example, it has been shown that out of 902 piRNAs expressed in somatic cells 527 are synthesized in the brain, of which piR-hsa-92056, piR-hsa-150797, piR-hsa-347751, piR-hsa-1909905, piR-hsa-2476,630, and piR-hsa-2834636 are recommended as biomarkers for Parkinson’s disease (Zhang and Wong, 2022).
Published articles showing the important role of piRNA in gastrointestinal cancer (Lin et al., 2019; Mojarad et al., 2022). The piRNAs dysregulations were reported to promote or suppress the initiation and development of different malignancies, especially gastrointestinal cancers. Recently, suggested the use of piRNAs as potential cancer biomarkers associated with the progression and chemoresistance of gastrointestinal cancer. Hence, this review article aims to focus on the role of piRNAs in gastrointestinal cancer progression, metastasis, and their molecular mechanisms as therapeutic markers for gastrointestinal cancer patients. X. Lin et al., 2022 identified 8,759 piRNAs in the human stomach. Of these, about 50 piRNAs could be differentially expressed compared to controls. It has been shown that piRNA/PIWI protein-mediated DNA methylation regulation mechanisms and methylation changes caused by piRNA/PIWI proteins in different diseases, especially cancers (Jia et al., 2022). Several studies suggest piRNAs as potential cancer biomarkers. Translational studies suggest that piRNAs may regulate key genes and pathways associated with gastric cancer progression, though there is no functional annotation in piRNA (Fonseca et al., 2020). Previously, 6,260 human piRNAs transcriptomes derived from benign and tumor tissues were analyzed, of which 522 piRNAs are expressed in the corresponding tumor tissues (Martinez et al., 2015). The authors suggest that piRNA expression may determine clinical features such as histological subgroups, disease stages, and survival. The piRNAs common to many cancer types may represent a core set of genes that contribute to cancer growth, while piRNAs unique to individual cancer types are likely specific.
The strongest stimulus for the study of the biological role of piRNA was the creation of an extensive database of piRNA properties (Wang et al., 2019; Wang et al., 2021). The publication of changes in piRNA synthesis during ontogenesis has made a significant contribution to the study of piRNA (Zhou et al., 2018). Recently, it was shown that piRNAs can bind to mRNAs of protein-coding genes involved in the development of various diseases (Belkozhayev et al., 2022). This information prompts us to investigate how and which piRNAs can regulate the expression of protein-coding genes. This work aims to show the possibility of the significant influence of piRNAs on the translation process through the interaction of these molecules with mRNAs of candidate genes involved in the development of esophageal adenocarcinoma (EAC).
Materials and methods
The nucleotide sequences of candidate genes of EAC were downloaded from NCBI. The nucleotide sequences of the first 200 thousand piRNAs of the 8,400,000 piRNAs database were taken from Wang et al. (Wang et al., 2021). The piRNA binding sites (BSs) in the mRNAs of genes were predicted using the MirTarget program (Ivashchenko et al., 2014). This program defines the following features of piRNA binding to mRNA: a) the initiation of piRNA binding from the first nucleotide of the mRNAs; c) the schemes of nucleotide interactions between piRNAs and mRNA; d) the free energy ΔG (kJ/mol) of the interaction between piRNA and the mRNA; e) the ratio ΔG/ΔGm (%) is determined for each site (ΔGm equals the free energy of piRNA binding with its fully complementary canonical nucleotide sequence). Only piRNAs whose nucleotides interacted with mRNA through canonical (G-C and A-U) and noncanonical (G-U and A-C) nucleotides with a ΔG/ΔGm value of 90% or more were selected from the calculated data. The 5′-untranslated region (5′UTR), coding domain sequence (CDS), and the 3′-untranslated region (3′UTR) of the mRNAs; c) the schemes of nucleotide interactions between miRNAs, piRNAs and mRNAs d) the free energy of the interaction between miRNAs, piRNAs and the mRNA (ΔG/ΔGm, %); and the ratio ΔG/ΔGm is determined for each site. The MirTarget program finds hydrogen bonds between adenine (A) and uracil (U), guanine (G) and cytosine (C), G and U, and A and C. Regarding the free energy of interactions (∆G), a pair of G and C is equal to 6.37 kJ/mol, a pair of A and U is equal to 4.25 kJ/mol, and a pair of G and U or A and C is equal to 2.12 kJ/mol. The distances between bound A and C (1.04 nm) and G and U (1.02 nm) are similar to those between bound G and C and between bound A and U and are equal to 1.03 nm (Friedman and Honig, 1995; Garg and Heinemann, 2018). The numbers of hydrogen bonds in the G–C, A–U, G–U, and A–C interactions were 3, 2, 1, and 1, respectively. By comparison, MirTarget differs from STarMir Tools, miRWalk programs (Kanoria et al., 2016; Sticht et al., 2018) in terms of finding the BSs of piRNA on the mRNAs in the following: (Kim, 2006): it accounts for the interaction of the piRNAs with mRNA over the entire piRNAs sequence; (Aravin et al., 2007); it considers noncanonical nucleotide pairs; and (Czech et al., 2018) it calculates the free energy of the interaction of the piRNAs with mRNA. When two or more piRNAs are bound with one mRNA, or if the BSs of two different piRNAs coincide in part, the preferred piRNA BS is considered to be the one for which the free energy binding is greater (Leontis et al., 2002; Davis et al., 2005).
Results
Androgen receptor (AR) is one of 38 candidate EAC genes (Bahramian et al., 2022) (Supplementary Table S1). Of the 200 thousand piRNAs we studied, only piR-32860 could bind to the mRNA of the AR gene. A feature of piR-32860 is the presence of two regions of BSs in the CDS mRNA located with overlapping nucleotide sequences, which we called clusters (Supplementary Table S2). The 15 piR-32860 BS starts at 1,287 nt and arranged sequentially through three nucleotides. They form a cluster whose nucleotides (CAG22) encode a polyglutamine of 23 amino acid residues. This triplet arrangement allows three piR-32860 molecules to bind to the mRNA simultaneously. The second cluster of BSs for piR-32860 is located in the CDS mRNA at 2,486 nt (Supplementary Table S2) and encodes for polyglycine. Consequently, several piR-32860 can simultaneously bind to the mRNA of the gene and significantly inhibit AR protein synthesis. The first cluster of BSs with a length of 67 nt is 5.4 times shorter than the total length of 15 piRNAs. The second cluster with a length of 49 nt is 4.4 times shorter than the total length of nine piRNAs. Such compaction of the BSs allows us to reduce significantly the length of the BSs in the coding region so that the rest of the CDS can vary and modify the functional activity of the AR gene.
The EAC suppressor protein BTG3 (Du et al., 2015) is encoded by an mRNA region located from 261 nt to 1,151 nt. The starting BSs of 27 piRNAs are located in the CDS mRNA from 573 nt to 665 nt with overlapping nucleotide sequences of piRNAs (Supplementary Table S3). The BSs of piR-12399, piR-37172, piR-126527, and piR-129369 start at 626 nt, which is approximately in the middle of the first BSs. The BSs cluster length of 27 piRNAs is 118 nt and is 7.6 times shorter than the CDS length and 135 times shorter than the entire mRNA length. Consequently, maintaining the nucleotide conserved nature of the cluster of piRNAs BSs requires a relatively small fraction of the mRNA nucleotide sequence to maintain the dependence of gene expression on 27 piRNAs. The second cluster of BSs of 11 piRNAs is located in the 3′UTR from 1,420 nt to 1,485 nt long by 66 nt and is 72.7 times shorter than the 3′UTR and 236 times shorter than the entire mRNA. This distribution of piRNAs BSs in the mRNA is not random, as it is organized into clusters. The piR-89068 and piR-83600, with a length of 34 nt, bind to mRNA with a ΔG value of −178 kJ/mol and −172 kJ/mol, respectively, indicating a high dependence of gene expression on these piRNAs. The binding of piR-75778 to mRNA occurs by forming almost only canonical nucleotide pairs, as the ΔG/ΔGm value is 99% (Supplementary Table S3). The piRNAs from piR-12399 to piR-129369, i.e., almost the entire interval of 200 thousand piRNAs, bind to the mRNA of the BTG3 gene. The distribution of piRNAs BSs in the form of clusters allows, firstly, compactification of BSs and, secondly, competition for binding to mRNA arises between piRNAs. The piRNA that interacts more strongly with the mRNA or the piRNA that is present in a much higher concentration compared to other piRNAs will predominantly bind to the mRNA.
Eleven piRNA BSs were identified in the mRNA CDS of the CD55 gene, located from 295 nt to 1,617 nt, with an origin from 1,420 nt to 1,458 nt (Supplementary Table S4). That is, the length of the cluster of piRNA BSs is 63 nt, which is 21 times shorter than the length of CDS. Among all the piRNAs, piR-83600 stands out, which binds to the mRNA of the CD55 gene with a free energy of −172 kJ/mol and a ΔG/ΔGm value of 92%. It is very likely that this piRNA can suppress the expression of the candidate CD55 gene associated with the development of EAC (Murao et al., 2016).
Expression of the ERBB3 gene (Yoon et al., 2014) can depend on 26 piRNAs, the BSs of which are located in the 3′UTR from 4,940 nt to 5,097 nt (Supplementary Table S5). These piRNA BSs formed a cluster from 4,940 nt to 5,125 nt, 186 nt long. The 3′UTR was 7.8 times as long as the BSs cluster. A feature of this cluster is the presence of an identical BS for seven piRNAs: piR-5300, piR-5303, piR-5358, piR-6236, piR-7637, piR-65119, piR-96686, each of which interacted with mRNA to approximately the same extent (Supplementary Table S5). The second cluster consisted of the BSs of nine piRNAs (Supplementary Table S5). The cluster length of 35 nt was 41.5 times shorter than the length of the 3′UTR.
The mRNA of the FKBP5 gene (Smith et al., 2016) contained five clusters of piRNA BSs located in 3′UTR the mRNA (Supplementary Table S6). In the first cluster, 20 piRNAs bound from 1,382 nt to 1,422 nt in the length of 72 nt. The second cluster starts from 4,854 nt (piR-123782) to 4,875 nt (piR-136223) and binds 20 piRNAs. This cluster contains five BSs with the same start for piRNA at position 4,862 nt (piR-56610, piR-62204, piR-75328, piR-92776, piR-102326). The ten piRNAs have BSs in the third BSs cluster at 4,940 nt to 4,965 nt. The fourth BSs cluster starts from 6,356 nt to 6,578 nt and contains 88 piRNA BSs. The BSs for 116 piRNAs start at position 7,100 nt to 7,201 nt. The piR-74093, piR-75328, piR-91782, piR-36610, piR-115392, piR-110734, piR-35260, piR-94289, piR-123273, piR-95238. piR-37678, piR-123924, piR-65906, piR-120553, piR-79631, piR-98746, piR-83452, piR-89965 binds with mRNA of FKBP5 gene with ΔG value from −170 kJ/mol to −183 kJ/mol. Of these, piR-65906 binds with a ΔG/ΔGm value of 98%, which means that almost all of it interacts through canonical nucleotide pairs. The piR-45567, piR-59915, piR-56610, piR-91782 bind with a value of ΔG/ΔGm equal to 99%, which also indicates their effectiveness on protein synthesis. Consequently, these piRNAs have the greatest ability to regulate FKBP5 gene expression. The piRNAs from piR-1248 to piR-198479 bind to the mRNA of FKBP5 gene, i.e., practically from the whole range of 200 thousand piRNAs studied. The total length of all clusters is 684 nt and is 9.4 times shorter than the entire 3′UTR region and 10.8 times shorter than the entire mRNA. Nine piRNAs (piR-5300, piR-5303, piR-5358, piR-6236, piR-7637, piR-35905, piR-65119, piR-96686, and piR-101606) from the third cluster bound at one position 7,115 nt (Supplementary Table S6). A group of five piRNAs (piR-45567, piR-93145, piR-96134, piR-123755, and piR-169382) bound at position 7,167 nt.
Expression of the LEP gene (Bain et al., 2014) may also depend on piRNAs, because 162 piRNAs can bind to only 3′UTR mRNA, including piR-23387, which is completely complementary due to canonical nucleotides (Supplementary Table S7). The BSs of 51 piNAs formed a cluster from 3,084 nt to 3,127 nt located in the 3′UTR. The length of the BSs cluster was 41 times less than the length of the 3′UTR. A group of nine piRNAs (piR-5300, piR-5303, piR-5358, piR-6236, piR-7637, piR-35905, piR-65119, piR-96686, and piR-101606) with BSs with 3,101 nt coincides with a group of piRNAs that bind to the mRNA of the FKBP5 gene and includes a group of piRNAs that interact with the mRNA of the ERBB3 gene. The mRNA of the gene contains BSs for three groups each of five piRNAs at positions 3,087 nt, 3,088 nt, and 3,246 nt piR-23387, piR-5299 and piR-1748 interacted with mRNA with ΔG/ΔGm equal to 99% and 100%.
The 46 piRNAs bind to the mRNA of the SEPP1 gene (Takata et al., 2012) from 108 nt to 172 nt (65 nt) at a 5′UTR length of 188 nt, i.e., 2.9 times shorter (Supplementary Table S8). The length of all 46 binding piRNAs is 1,367 nt, and due to the compaction of piRNAs BSs into clusters, competition between piRNAs for binding in the cluster simultaneously arises. At position 123 nt of the 5′UTR of the mRNA of the SEPP1 gene, a group of nine piRNAs was identified, identical to the groups of piRNAs that bind to the mRNA of the ERBB3, FKBP5 and LEP genes.
The 74 piRNAs bound to the 3′UTR of the SMAD4 gene (Weaver et al., 2014). The mRNA and their BSs formed clusters from 4,316 nt to 4,563 nt (Supplementary Table S9). In the cluster, three groups of piRNAs were identified. The first group from 4,335 nt included piR-36976, piR-50827, piR-84920, piR-116971, piR-125671. The second group from 4,337 nt consisted of piR-19014, piR-19076, piR-82088, piR-125598. The third group from 4,401 nt was formed by piR-56309, piR-81517, piR-98951. The cluster was only 247 nt from the 3′UTR length of 6,576 nt, which is 26.6 times the length of the cluster. The piR-82439, piR-84833, piR-100395 interacted with mRNA with a ΔG/ΔGm value of 99%.
There were no clusters of BSs in the mRNA of the TP53 gene (Fisher et al., 2017) for groups with 5 piRNAs or more (Supplementary Table S10). The cluster of 27 piRNAs BSs with a length of 224 nt was 6.3 times less than the 3′UTR length of 1,408 nt. The piR-44059 and piR-83728 bound to the mRNA of the TP53 gene with ΔG kJ/mol equal to -170 kJ/mol and 172 kJ/mol, respectively, i.e., they could strongly inhibit the translation process.
Based on these studies, the associations of piRNA and candidate genes effectively interacting with a ΔG value of -170 kJ/mol or more are proposed as disease markers. Figure 1 shows the interaction patterns of piRNA with mRNA of candidate EAC genes with a ΔG value of -170 kJ/mol or more. No piRNAs with this interaction energy have been identified for the AR gene (Supplementary Table S2). The piR-89068, piR-83600 bind to the CDS mRNA of the BTG3 gene. Only piR-83600 interacts with the CDS mRNA of the CD55 gene. Consequently, piR-83600 can regulate the expression of BTG3 and CD55 genes. These findings show the need to study piRNA and gene associations in diseases specifically, and not just piRNA or gene associations due to the ambiguity of interpreting results separately for piRNA or genes.
FIGURE 1. The schemes of interaction between piRNAs and mRNAs of EAC candidate genes with ΔG values of −170 kJ/mol and more.
In selecting piRNAs as markers of a disease caused by a candidate gene, it is necessary to know how selectively the piRNA alters the expression of that gene. Therefore, it is necessary to establish how this (these) piRNAs can change the expression of other genes. We used 200 thousand piRNAs to identify their interactions with the mRNA of 38 candidate EAC genes (Supplementary Table S1) and showed that only nine genes were targeted by the studied piRNAs. From two piRNAs (AR gene) to 254 piRNAs (FKBP5 gene) that bind to mRNA with a ΔG/ΔGm value of 90% or more could interact with these genes. Increasing the selection criterion for piRNAs (ΔG/ΔGm ratio) from 90% to 95% reduces the number of interacting piRNAs by about three times. In cells, the concentration of piRNAs can vary by several orders of magnitude, then a piRNA interacting is weaker than another piRNA but present in a higher concentration will more effectively suppress the expression of the candidate gene. Therefore, we selected piRNAs with a ΔG/ΔGm value of 90% or higher to account for piRNAs potentially being able to bind to the mRNA of the candidate gene. Figure 2 shows as an example the interaction schemes of piRNAs with the mRNA of candidate EAC genes with ΔG/ΔGm values of more than 97% or higher.
FIGURE 2. The schemes of interaction between piRNAs and mRNAs of EAC candidate genes with ΔG/ΔGm) values greater than 97% or more.
This raises the question of whether the free energy of the interaction of piRNA with mRNA (ΔG) or the degree of complementarity of the binding of piRNA to mRNA (ΔG/ΔGm) is a more important factor in assessing the interaction of piRNA with mRNA. Interestingly, the group of piRNA selected on the principle of free energy greater than -170 kJ/mol and the group of piRNA associating with mRNA with ΔG/ΔGm greater than 90% tends to have different clusters of piRNA BSs, which demonstrates the importance of both criteria for the interaction of piRNA with mRNA.
An important aspect of the choice of associations of piRNAs and candidate genes is that one gene may be involved in the development of several diseases and then that gene, regulated even by one piRNA, may be associated by side effects with another disease. Modern drugs come with annotations about the limitations and side effects of their use. Similarly, piRNAs chosen for diagnosis and therapy should be accompanied by a similar annotation. Therefore, laboratory studies of the proposed piRNAs in combination with candidate genes for these diseases are necessary to the evaluative bioinformatics studies of the involvement of piRNAs in various diseases.
The highest number of piRNAs were associated with a ΔG value of −170 kJ/mol or more with the 3′UTR mRNA of the FKBP5 gene: piR-74093, piR-38669, piR-75328, piR-56610, piR-115392, piR-110734, piR-35260, piR-55529, piR-94289, piR-196859, piR-127715, piR-89432, piR-151720, piR-123273, piR-106904, piR-95238, piR-37678, piR-123924, piR-65906, piR-120553, piR-79631, piR-98746, piR-83452, and piR-89965. Six piRNAs with ΔG values ranging from -170 kJ/mol to -183 kJ/mol bind from 7,160 nt to 7,165 nt in the 3′UTR mRNA of the FKBP5 gene. This indicates that these piRNAs compete with each other and their suppressive effect on FKBP5 protein synthesis is high.
The piR-127715, piR-89432, piR-199381, piR-23387, piR-39121, piR-97487, piR-73272, piR-77917, piR-78890, piR-87066 can bind to 3′UTR mRNA of LEP gene with ΔG value of -170 kJ/mol or more. Note that piR-73272, piR-77917, and piR-78890 bound to the mRNA of the LEP gene with a ΔG value of −183 kJ/mol, −183 kJ/mol, and −189 kJ/mol, respectively, with overlapping nucleotide sequences of BSs. Consequently, these piRNAs cannot bind simultaneously, but their overall effect is high, even if they are present in the cell at different concentrations.
The piR-89432 and piR-151720 binded to a ΔG value of -170 kJ/mol or more in the 5′UTR mRNA of the SEPP1 gene. The piR-125671, piR-82439, piR-125598, piR-98281, piR-94289, piR-83728, and piR-58129 interacted with the 3′UTR mRNA of the SMAD4 gene with a ΔG value of −170 kJ/mol or more. Only piR-83728 and piR-44059 bind to the 3′UTR mRNA of the TP53 gene with a free energy of -172 kJ/mol and −170 kJ/mol, respectively. These results demonstrate the involvement of piR-83728 in the regulation of SMAD4 and TP53 gene expression. This is another example of the need to study piRNA and gene associations rather than individual piRNAs or genes.
We identified piR-12340, piR-48382, piR-48873, piR-55529, piR-59915, piR-96587, piR-99501, and piR-110973 having BSs in mRNA clusters of FKBP5 and SMAD4 genes. Consequently, these piRNAs can simultaneously regulate the expression of FKBP5 and SMAD4 genes.
An even more surprising result was obtained with the ERBB3, FKBP5, LEP, and SEPP1 genes. The mRNAs of these genes contain highly homologous BSs region of 27 piRNAs 54 nt long (Figure 3, (Supplementary Table S11). These data show that the expression of the ERBB3, FKBP5, LEP, and SEPP1 genes can be maintained in a balanced manner by this number of piRNAs. We obtained a similar result when we studied the action of several miRNAs on a number of candidate atherosclerosis genes (Mukushkina et al., 2020).
FIGURE 3. Nucleotide sequences of 27 piRNAs clusters of BSs in the 3′UTR mRNA of the ERBB3, FKBP5, and LEP genes and in the 5′UTR mRNA of the SEPP1 gene.
Discussion
Many human genes are functionally related, so a single gene can exert its effect in several diseases. For example, the FKBP5 gene encodes a Prolyl isomerase that binds to immunosuppressants and interacts functionally with a number of proteins. Changes in FKBP5 gene expression were found in gastric cancer (Kang et al., 2012; Wang et al., 2022b), pancreatic cancer (Hou and Wang, 2012), breast cancer (Xiong et al., 2020), and papillary thyroid carcinoma (Gao et al., 2021). Protein, encoded by the FKBP5 gene, is responsible for stress and metabolic-related disorders, including cancer (Barge et al., 2021). The association of AR gene and FKBP5 gene as independent prognostic indicators for EAC has been shown (Smith et al., 2016). Given these circumstances, the above associations of piRNA and candidate EAC genes recommended for diagnosis of this disease should be considered as probable molecular markers.
Based on these results, several combinations of piRNA and candidate genes can be recommended for the diagnosis. The least expensive way is to use associations of one piRNA and one candidate gene, such as piR-65906 and the FKBP5 gene, or the LEP gene and piR-73272, piR-77917, piR-78890. Another way is to use multiple piRNAs and two candidate genes, such as eight piRNAs and the FKBP5 and SMAD4 genes. Another way to diagnose is to control many piRNAs (piR-5299, piR-5300, piR-5303, piR-5358, piR-6236, piR-7637, piR-19207, piR-34287, piR-65119, piR-89432, piR-92948, piR-94621, piR-95442, piR-96686, piR-100284, piR-101948, piR-108567, piR-127960, piR-151720, and piR-197022) in combination with several candidate genes in EAC (ERBB3, FKBP5, LEP, and SEPP1), which is certainly more appropriate, but more effective because it covers many piRNAs and candidate genes capable of causing EAC.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
AI conceived of the study, drafted the manuscript and gave final approval of the version to be published. MZ conceived of the study, drafted the manuscript. AA conceived of the study and drafted the manuscript and gave final approval of the version to be published. TN conceived of the study and drafted the manuscript. AP software design. OY conceived of the study, drafted the manuscript. All authors made substantial contributions to acquisition of data, to interpretation and modification of the data, were involved in subsequent rounds of revisions, and read and approved the final manuscript.
Funding
The study is supported by Center for Bioinformatics and Nanomedicine, Almaty, Kazakhstan.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1069637/full#supplementary-material
References
Ameli Mojarad, M., Ameli Mojarad, M., and Wang, J. (2022). The function of novel small non-coding RNAs (piRNAs, tRFs) and PIWI protein in colorectal cancer. Cancer Treat. Res. Commun. 31, 100542. doi:10.1016/j.ctarc.2022.100542
Aravin, A. A., Sachidanandam, R., Girard, A., Fejes-Toth, K., and Hannon, G. J. (2007). Developmentally regulated piRNA clusters implicate MILI in transposon control. Science 316 (5825), 744–747. doi:10.1126/science.1142612
Bahramian, S., Sahebi, R., Roohinejad, Z., Delshad, E., Javid, N., Amini, A., et al. (2022). Low expression of LncRNA-CAF attributed to the high expression of HIF1A in esophageal squamous cell carcinoma and gastric cancer patients. Mol. Biol. Rep. 49 (2), 895–905. doi:10.1007/s11033-021-06882-0
Bain, G. H., Collie-Duguid, E., Murray, G. I., Gilbert, F., J., Denison, A., McKiddie, F., et al. (2014). Tumour expression of leptin is associated with chemotherapy resistance and therapy-independent prognosis in gastro-oesophageal adenocarcinomas. Br. J. Cancer 110 (6), 1525–1534. doi:10.1038/bjc.2014.45
Barge, S., Jade, D., Ayyamperumal, S., Manna, P., Borah, J., Nanjan, C. M. J., et al. (2021). Potential inhibitors for FKBP51: An in silico study using virtual screening, molecular docking and molecular dynamics simulation study using virtual screening, molecular docking and molecular dynamics simulation. J. Biomol. Struct. Dyn. 28, 1–13. doi:10.1080/07391102.2021.1994877
Belkozhayev, A., Niyazova, R., Wilson, C., Jainakbayev, N., Pyrkova, A., Ashirbekov, Y., et al. (2022). Bioinformatics analysis of the interaction of miRNAs and piRNAs with human mRNA genes having di- and trinucleotide repeats. Genes 13, 800. doi:10.3390/genes13050800
Cai, A., Hu, Y., Zhou, Z., Qi, Q., Wu, Y., Dong, P., et al. (2022). PIWI-interacting RNAs (piRNAs): Promising applications as emerging biomarkers for digestive system cancer. Front. Mol. Biosci. 27 (9), 848105. doi:10.3389/fmolb.2022.848105
Chen, S., Ben, S., Xin, J., Li, S., Zheng, R., Wang, H., et al. (2021). The biogenesis and biological function of PIWI-interacting RNA in cancer. J. Hematol. Oncol. 14 (1), 93. doi:10.1186/s13045-021-01104-3
Czech, B., Munafò, M., Ciabrelli, F., Eastwood, E. L., Fabry, M. H., Kneuss, E., et al. (2018). piRNA-guided genome defense: From biogenesis to silencing. Annu. Rev. Genet. 52, 131–157. doi:10.1146/annurev-genet-120417-031441
Davis, E., Caiment, F., Tordoir, X., Cavaillé, J., Ferguson-Smith, C., N., Georges, M., et al. (2005). RNAi-mediated allelic trans-interaction at the imprinted Rtl1/Peg11 locus. Curr. Biol. 15, 743–749. doi:10.1016/j.cub.2005.02.060
Du, Y., Liu, P., Zang, W., Wang, Y., Chen, X., Li, M., et al. (2015). BTG3 upregulation induces cell apoptosis and suppresses invasion in esophageal adenocarcinoma. Mol. Cell. Biochem. 404 (1-2), 31–38. doi:10.1007/s11010-015-2363-9
Fisher, O. M., Lord, S. J., Falkenback, D., Clemons, N. J., Eslick, G. D., and Lord, R. V. (2017). The prognostic value of TP53 mutations in oesophageal adenocarcinoma: A systematic review and meta-analysis. Gut 66 (3), 399–410. doi:10.1136/gutjnl-2015-310888
Fonseca, C. G., Azevedo, D., Santos Pinheiro, J., Vidal, A. F., Santos, S., and Ribeiro-Dos-Santos, Â. (2020). piRNAs in gastric cancer: A new approach towards translational research. Int. J. Mol. Sci. 21 (6), 2126. doi:10.3390/ijms21062126
Friedman, R. A., and Honig, B. A. (1995). A free energy analysis of nucleic acid base stacking in aqueous solution. Biophys. J. 69 (4), 1528–1535. doi:10.1016/S0006-3495(95)80023-8
Gao, Z., Yu, F., Jia, H., Ye, Z., and Yao, S. (2021). FK506-binding protein 5 promotes the progression of papillary thyroid carcinoma. J. Int. Med. Res. 49 (4), 3000605211008325. doi:10.1177/03000605211008325
Garg, A., and Heinemann, U. A. (2018). A novel form of RNA double helix based on G·U and C·A+ wobble base pairing. RNA 24, 209–218. doi:10.1261/rna.064048.117
Guo, B., Li, D., Du, L., and Zhu, Xl. (2020). piRNAs: biogenesis and their potential roles in cancer. Cancer Metastasis Rev. 39 (2), 567–575. doi:10.1007/s10555-020-09863-0
Halajzadeh, J., Dana, P. M., Asemi, Z., Mansournia, M. A., and Yousefi, B. (2020). An insight into the roles of piRNAs and PIWI proteins in the diagnosis and pathogenesis of oral.; esophageal.; and gastric cancer. Pathol. Res. Pract. 216, 153112. doi:10.1016/j.prp.2020.153112
Hanusek, K., Poletajew, S., Kryst, P., Piekiełko-Witkowska, A., and Bogusławska, J. (2022). piRNAs and PIWI proteins as diagnostic and prognostic markers of genitourinary cancers. Biomolecules 12 (2), 186. doi:10.3390/biom12020186
Hou, J., and Wang, L. (2012). FKBP5 as a selection biomarker for gemcitabine and Akt inhibitors in treatment of pancreatic cancer. PLoS One 7 (5), e36252. doi:10.1371/journal.pone.0036252
Ivashchenko, A., Berillo, O., Pyrkova, A., Niyazova, R., and Atambayeva, S. (2014). MiR-3960 binding sites with mRNA of human genes. Bioinformation 10, 423–427. doi:10.6026/97320630010423
Jia, D. D., Jiang, H., Zhang, Y. F., Zhang, Y., Qian, L. L., and Zhang, Y. F. (2022). The regulatory function of piRNA/PIWI complex in cancer and other human diseases: The role of DNA methylation. Int. J. Biol. Sci. 18 (8), 3358–3373. doi:10.7150/ijbs.68221
Kang, J. I., Chung, H. C., Jeung, H. C., Kim, S. J., An, S. K., and Namkoong, K. (2012). FKBP5 polymorphisms as vulnerability to anxiety and depression in patients with advanced gastric cancer: A controlled and prospective study. Psychoneuroendocrinology 37 (9), 1569–1576. doi:10.1016/j.psyneuen.2012.02.017
Kanoria, S., Rennie, W., Liu, S., Carmack, C. S., Lu, J., and Ding, Y. (2016). STarMir Tools for prediction of microRNA binding sites. Methods Mol. Biol. 1490, 73–82. doi:10.1007/978-1-4939-6433-8_6
Kärkkäinen, E., Heikkinen, S., Tengström, M., Kosma, V. M., Mannermaa, A., and Hartikainen, J. M. (2021). The debatable presence of PIWI-interacting RNAs in invasive breast cancer. Cancer Med. 10 (11), 3593–3603. doi:10.1002/cam4.3915
Kim, V. N. (2006). Small RNAs just got bigger: Piwi-interacting RNAs (piRNAs) in mammalian testes. Genes Dev. 20, 1993–1997. doi:10.1101/gad.1456106
Leontis, N. B., Stombaugh, J., and Westhof, E. (2002). The non-watson-crick base pairs and their associated isostericity matrices. Nucleic Acids Res. 30, 3497–3531. doi:10.1093/nar/gkf481
Lin, X., Xia, Y., Hu, D., Mao, Q., Yu, Z., Zhang, H., et al. (2019). Transcriptome-wide piRNA profiling in human gastric cancer. Oncol. Rep. 41 (5), 3089–3099. doi:10.3892/or.2019.7073
Liu, T., Wang, J., Sun, L., Li, M., He, X., Jiang, J., et al. (2021). Piwi-interacting RNA-651 promotes cell proliferation and migration and inhibits apoptosis in breast cancer by facilitating DNMT1-mediated PTEN promoter methylation. Cell Cycle 20 (16), 1603–1616. doi:10.1080/15384101.2021.1956090
Liu, Y., Dou, M., Song, X., Dong, Y., Liuk, S., Liu, H., et al. (2019). The emerging role of the piRNA/piwi complex in cancer. Mol. Cancer 18 (1), 123. doi:10.1186/s12943-019-1052-9
Martinez, V. D., Vucic, E. A., Thu, K. L., Hubaux, R., Enfield, K. S., Pikor, L. A., et al. (2015). Unique somatic and malignant expression patterns implicate PIWI-interacting RNAs in cancer-type specific biology. Sci. Rep. 5, 10423. doi:10.1038/srep10423
Mojarad, A. M., Mojarad, A. M., Shojaee, B., and Nazemalhosseini-Mojarad, E. (2022). piRNA: A promising biomarker in early detection of gastrointestinal cancer. Pathol. Res. Pract. 230, 153757. doi:10.1016/j.prp.2021.153757
Mukushkina, D., Aisina, D., Pyrkova, A., Ryskulova, A., Labeit, S., Ivashchenko, A., et al. (2020). In silico prediction of miRNA interactions with candidate atherosclerosis gene mRNAs. Front. genet. 11, 605054. doi:10.3389/fgene.2020.605054
Murao, T., Shiotani, A., Fujita, Y., Yamanaka, Y., Kamada, T., Manabe, N., et al. (2016). Overexpression of CD55 from Barrett's esophagus is associated with esophageal adenocarcinoma risk. J. Gastroenterol. Hepatol. 31 (1), 99–106. doi:10.1111/jgh.13055
Qian, L., Xie, H., Zhang, L., Zhao, Q., Lü, J., and Yu, Z. (2021). Piwi-interacting RNAs: A new class of regulator in human breast cancer. Front. Oncol. 11, 695077. doi:10.3389/fonc.2021.695077
Smith, E., Palethorpe, H. M., Ruszkiewicz, A. R., Edwards, S., Leach, D. A., Underwood, T. J., et al. (2016). Androgen receptor and androgen-responsive gene FKBP5 are independent prognostic indicators for esophageal adenocarcinoma. Dig. Dis. Sci. 61 (2), 433–443. doi:10.1007/s10620-015-3909-0
Sticht, C., De La Torre, C., Parveen, A., and Gretz, N. (2018). miRWalk: An online resource for prediction of microRNA binding sites. PLoS ONE 13, e0206239. doi:10.1371/journal.pone.0206239
Takata, Y., Kristal, A. R., Santella, R. M., King, I. B., Duggan, D. J., Lampe, J. W., et al. (2012). Selenium, selenoenzymes, oxidative stress and risk of neoplastic progression from barrett's esophagus: Results from biomarkers and genetic variants. PLoS One 7 (6), e38612. doi:10.1371/journal.pone.0038612
Tan, L., Mai, D., Zhang, B., Jiang, X., Zhang, J., Bai, R., et al. (2019). PIWI-interacting RNA-36712 restrains breast cancer progression and chemoresistance by interaction with SEPW1 pseudogene SEPW1P RNA. Mol. Cancer 18 (1), 9. doi:10.1186/s12943-019-0940-3
Wang, J., Sh, Y., Zhou, H., Zhang, P., Song, T., Ying, Z., et al. (2021). piRBase: integrating piRNA annotation in all aspects. Nucleic Acids Res. 50, 265–272. doi:10.1093/nar/gkab1012
Wang, J., Zhang, P., Lu, Y., Li, Y., Zheng, Y., Kan, Y., et al. (2019). piRBase: a comprehensive database of piRNA sequences. Nucleic Acids Res. 947, D175–D180. doi:10.1093/nar/gky1043
Wang, X., Gou, L. T., and Liu, M. F. (2022). Noncanonical functions of PIWIL1/piRNAs in animal male germ cells and human diseases. Biol. Reprod. 107 (1), 101–108. doi:10.1093/biolre/ioac073
Wang, Y., Li, T., Yang, L., Zhang, X., Wang, X., Su, X., et al. (2022). Cancer-associated fibroblast-released extracellular vesicles carrying miR-199a-5p induces the progression of gastric cancer through regulation of FKBP5-mediated AKT1/mTORC1 signaling pathway. Cell Cycle, 1–12. doi:10.1080/15384101.2022.2105092
Weaver, J. M. J., Ross-Innes, C. S., Shannon, N., Lynch, A. G., Forshew, T., Barbera, M., et al. (2014). Ordering of mutations in preinvasive disease stages of esophageal carcinogenesis. Nat. Genet. 46 (8), 837–843. doi:10.1038/ng.3013
Xiong, H., Chen, Z., Chen, W., Li, Q., Lin, B., and Jia, Y. (2020). FKBP-related ncRNA-mRNA axis in breast cancer. Genomics 112 (6), 4595–4607. doi:10.1016/j.ygeno.2020.08.017
Yoon, H., H., Sukov, W. R., Shi, Q., Sattler, C, A., Wiktor, A. E., Diasio, R. B., et al. (2014). HER-2/neu gene amplification in relation to expression of HER2 and HER3 proteins in patients with esophageal adenocarcinoma. Cancer 120 (3), 415–424. doi:10.1002/cncr.28435
Yuan, C., Qin, H., Ponnusamy, M., Chen, Y., and Lin, Z. (2021). PIWI‑interacting RNA in cancer: Molecular mechanisms and possible clinical implications (Review). Oncol. Rep. 46 (3), 209. doi:10.3892/or.2021.8160
Zhang, J., Chen, S., and Liu, K. (2022). Structural insights into piRNA biogenesis. Biochim. Biophys. Acta. Gene Regul. Mech. 1865 (2), 194799. doi:10.1016/j.bbagrm.2022.194799
Zhang, L. M., Gao, Q. X., Chen, J., Li, B., Li, M. M., Zheng, L., et al. (2022). A universal catalytic hairpin assembly system for direct plasma biopsy of exosomal PIWI-interacting RNAs and microRNAs. Anal. Chim. Acta 1192, 339382. doi:10.1016/j.aca.2021.339382
Zhang, T., Chen, L., Li, R., Liu, N., Huang, X., and Wong, G. (2022). PIWI-Interacting RNAs in human diseases: Databases and computational models. Brief. Bioinform. 23 (4), bbac217. doi:10.1093/bib/bbac217
Zhang, T., and Wong, G. (2022). Dysregulation of human somatic piRNA expression in Parkinson's disease subtypes and stages. Int. J. Mol. Sci. 23 (5), 2469. doi:10.3390/ijms23052469
Keywords: esophageal adenocarcinoma, piRNA, mRNA, candidate genes, diagnostic
Citation: Akimniyazova AN, Niyazova TK, Yurikova OY, Pyrkova AY, Zhanuzakov MA and Ivashchenko AT (2022) piRNAs may regulate expression of candidate genes of esophageal adenocarcinoma. Front. Genet. 13:1069637. doi: 10.3389/fgene.2022.1069637
Received: 14 October 2022; Accepted: 18 November 2022;
Published: 30 November 2022.
Edited by:
Yuriy L. Orlov, I. M. Sechenov First Moscow State Medical University, RussiaReviewed by:
Stepan Nersisyan, National Research University Higher School of Economics, RussiaMikhail P. Ponomarenko, Institute of Cytology and Genetics (RAS), Russia
Elvira Galieva, Novosibirsk State University, Russia
Copyright © 2022 Akimniyazova, Niyazova, Yurikova, Pyrkova, Zhanuzakov and Ivashchenko. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: A. T. Ivashchenko, YS5pYXZhc2hjaGVua29AZ21haWwuY29t