Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 07 December 2023
Sec. RNA

Impact of MICA 3′UTR allelic variability on miRNA binding prediction, a bioinformatic approach

  • 1Faculty of Medicine, Immunology Program, Institute of Biomedical Sciences (ICBM), Universidad de Chile, Santiago, Chile
  • 2School of Biochemical Engineering, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile
  • 3Faculty of Medicine, Clinical and Molecular Pharmacology Program, Institute of Biomedical Sciences (ICBM), Universidad de Chile, Santiago, Chile
  • 4Gastroenterology and Hepatology Department, University Medical Center Groningen, Groningen, Netherlands

MicroRNAs (miRNAs) are small non-coding RNAs that participate as powerful genetic regulators. MiRNAs can interfere with cellular processes by interacting with a broad spectrum of target genes under physiological and pathological states, including cancer development and progression. Major histocompatibility complex major histocompatibility complex class I-related chain A (MICA) belongs to a family of proteins that bind the natural-killer group 2, member D (NKG2D) receptor on Natural Killer cells and other cytotoxic lymphocytes. MICA plays a crucial role in the host’s innate immune response to several disease settings, including cancer. MICA harbors various single nucleotide polymorphisms (SNPs) located in its 3′-untranslated region (3′UTR), a characteristic that increases the complexity of MICA regulation, favoring its post-transcriptional modulation by miRNAs under physiological and pathological conditions. Here, we conducted an in-depth analysis of MICA 3′UTR sequences according to each MICA allele described to date using NCBI database. We also systematically evaluated interactions between miRNAs and their putative targets on MICA 3′UTR containing SNPs using in silico analysis. Our in silico results showed that MICA SNPs rs9266829, rs 1880, and rs9266825, located in the target sequence of miRNAs hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-93, hsa-miR-1207.5p, and hsa-miR-711 could modify the binding free energy between −8.62 and −18.14 kcal/mol, which may affect the regulation of MICA expression. We believe that our results may provide a starting point for further exploration of miRNA regulatory effects depending on MICA allelic variability; they may also be a guide to conduct miRNA in silico analysis for other highly polymorphic genes.

1 Introduction

MicroRNAs (miRNAs) are small non-coding RNAs containing 20–22 nucleotides (nts) that can significantly change gene expression through post-transcriptional regulation by driving target mRNA deadenylation, degradation and translational repression (Huntzinger and Izaurralde, 2011; Shah and Shah, 2020). MiRNAs can affect the expression of one or multiple genes under physiological and pathological states, and dysregulation mechanisms affecting miRNAs can potentially alter biological pathways involved in cancer development and progression. In mammalians, the biogenesis of miRNAs is a highly conserved process that involves the action of type II and III polymerases (Borchert et al., 2006), resulting in the transcription, from a DNA sequence, of a single strand of primary miRNA (pri-miRNA) (∼1,000 nts). Pri-miRNA is processed by a complex formed by the RNA-binding protein DiGeorge Syndrome Critical Region 8 (DGCR8) and a ribonuclease III enzyme, DROSHA, giving rise to a pre-miRNA transcript (∼60 nts), which is further processed by the RNAse III endonuclease Dicer, resulting in a mature miRNA duplex. Argonaute 2 (AGO2), a protein with catalytic activity, associates with miRNA, and this complex generates the miRNA-induced silencing complex (miRISC), which drives the interactions between the miRNAs and their messenger RNAs (mRNAs) targets (OBrien et al., 2018). This interaction depends on the particular configuration between miRNA:mRNA complementarity, which includes the characteristic annealing in the “seed sequence region” (six to eight base pairs from the miRNA 5‘end), and also non-canonical pairing (including nucleotides in miRNA 3′end) conferring stability to the RNA interaction and promoting the effect of miRISC (Brennecke et al., 2005; Wang, 2014). These recognition patterns can be damped by slight changes in the target nucleotide sequence, such as DNA single nucleotide polymorphisms (SNPs) associated with the 3′untranslated region (3′UTR) of a given mRNA, and potentially modify the effect of functional miRNA (Cai et al., 2009). Thus, SNPs in the 3′UTR may abolish or create novel miRNA binding sites (Manikandan and Munirajan, 2014), and polymorphic genes could be more prone to this type of plasticity.

The major histocompatibility complex (MHC) class I-related chain A (MICA) is a stress-induced ligand that plays an important role in the immune response to tumors. It binds to natural-killer group 2, member D (NKG2D) receptor on Natural Killer (NK) cells and other cytotoxic lymphocytes, which triggers the activation of effector cell cytolytic activity and becomes the signal for interferon-gamma (IFN-y) production, thus activating the host immune response mainly against intracellular pathogens and tumor cells (Mah and Cooper, 2016; Fuertes et al., 2021). MICA is a highly polymorphic gene located in chromosome 6 at 6p21.3; it contains six exons: exon 1 encodes the leader peptide, exons 2-4 encode the three extracellular domains α1, α2, and α3, respectively, exon 5 encodes the transmembrane region, and exon 6 encodes the hydrophobic carboxy-terminal cytoplasmic tail or it can be part of 3′UTR in the case of MICA allele * 008 (Vitiani et al., 1998; Li et al., 1999). Currently, more than two hundred MICA full-sequences have been described, according to the IPD-IMGT/HLA Database (IPD-IMGT/HLA, 2023). MICA*008 is the most prevalent allele in the human population, followed by MICA*002; together, they comprise almost 50% of the MICA alleles described to date (The Allele Frequency Net Database, 2022). Some of them have been associated with increased risks of developing cancer, including gastric (MICA*009/049) (Toledo-Stuardo et al., 2021), and colorectal cancer (MICA*012:01) (Ding et al., 2020). Therefore, MICA allelic distribution may contribute to disease risks in specific human populations (Lucas et al., 2008; Wenda et al., 2013; Wang et al., 2016; Yamakawa et al., 2017; Klussmeier et al., 2020).

MiRNAs that have a direct effect on the MICA transcript have been associated with degradation and reduction of MICA expression on the tumor cell membrane, which may affect the cytotoxic response mediated by NK cells (Xie et al., 2014). Hence, miRNAs address a fine regulation of transcript, regulating MICA levels mainly in stress conditions (Stern-Ginossar et al., 2008). Post-transcriptional regulation of MICA by miRNAs has been reported in different types of cancer and clinically associated with relevant tumor characteristics and disease progression. Some of these miRNAs include hsa-miR-25, hsa-miR-93, hsa-miR-106a, hsa-miR-106b and hsa-miR-373, which were detected in hepatocellular carcinoma cell lines (Kishikawa et al., 2013; Wu et al., 2014; Yang et al., 2015); hsa-miR-146b-5p, which has been related with primary tumor-derived papillary thyroid carcinoma cells (Al-Abdallah et al., 2020); hsa-miR-519a-3p, detected in breast cancer (Breunig et al., 2017), hsa-miR-302c and hsa-miR-520c, present in multiple cancer cell lines (Kasumi-1, K-562, MDA-MB-231, MCF-7, and HEK-293T) (Min et al., 2013). In contrast, a miRNA that upregulates MICA expression has also been reported, such as miR-125b, whose upregulation leads to increased MICA levels in human myeloma cells, affecting the immune response against tumors (Abruzzese et al., 2016). The presence of SNPs located in the 3′UTR region of the MICA gene with very strong linkage disequilibrium (LD) allowed the classification of the MICA 3′UTR into seven types, UTR1-7 (Luo et al., 2013). These SNPs can affect the strength of miRNAs interaction. MICA gene harbors more than 300 SNPs (MICA ENSG00000204520, variant table - Ensembl release 109 (Gene: MICA ENSG00000233051, 2023; Martin et al., 2023)) that may permit a gain or loss of allele-dependent miRNA target sequence recognition, which gives plasticity to miRNA regulation mechanisms, thus affecting protein expression. Such versatility of miRNA targets conferred by genetic variants of MICA gene has not been completely approached.

Here, we studied miRNAs with potential novel binding sites to SNPs located in MICA 3′UTR by in silico analysis; we also discuss their possible implications in the regulation of MICA expression in cancer, emphasizing the relevance of MICA alleles in this post-transcriptional modulation. We hypothesize that the miRNAs evaluated herein could determine MICA expression depending on the allele. We believe our data contribute to the knowledge of how MICA alleles variability can influence miRNA binding, and present rational in silico evidence for future MICA expression studies.

2 Methods

2.1 MICA 3′UTRs identification and analysis

The most frequent MICA alleles worldwide were selected for analysis: MICA*001, *002, *004, *007, *008, *009, *010, *011, *012, *015, *017, *018, *019, *027, *049, *068 (http://www.allelefrequencies.net). NCBI reference sequence used for analysis was NM_001177519.3. The MICA 3′UTR mRNA polymorphic sequences were reviewed on IPD-IMGT/HLA Database, and the initial sequence and length of 3′ UTR regions were analyzed by using the IMGTplatform (IPD-IMGT/HLA, 2023; Lefranc et al., 2003)).

2.1.1 Classification of MICA 3′UTR types

In previous studies, the classification of MICA 3′UTR types was based on the presence of 9 SNPs; rs9266826 (G/A), rs9266827 (A/G), rs9266828 (G/C), rs9266929 (A/G), rs140390705 (C/T), rs 1880 (C/A), rs113015830 (A/T), rs1131904 (A/G), and rs9266831 (G/A) were grouped into UTR1-UTR7 (Luo et al., 2013). However, in the present study, we introduced a new 3′UTR classification that considers: 1) the introduction of two SNPs, rs9266825(C/A) and rs 1882 (A/G), as they are located in regions with LD with the other 9 SNPs described (from position 1,000 to 1,327 (NM_001177519.3)); 2) exclusion of UTR5 and UTR6, which harbor the SNPs rs140390705 allele T and rs113015830 allele A, respectively, as they have a lower global frequency (0.0006 and 0.0044, respectively), and they have not been identified as part of MICA alleles reported in the IPD-IMGT/HLA Database (Table 1); 3) in consideration of a guanine insertion in the transmembrane region at position 952 located in exon 5 (MICA *008), and a deletion of guanine at position 892 in exon 4 (MICA *015, and *017), both of which change the reading frame, resulting in premature stop codons TAA or TGA at positions +997 and +1,035, respectively. The rest of the alleles maintain TAG stop codon at position +1,154. Thus, MICA*008 3′UTR is longer than the 3′UTR from MICA *015, and *017 MICA alleles and the other MICA alleles described to date. MICA alleles reported in the IPD-IMGT/HLA Database are included in Table 1.

TABLE 1
www.frontiersin.org

TABLE 1. Diversity of 3′UTR of the MICA gene and novel classification.

2.2 Recognition of miRNAs binding sites to SNPs located in MICA 3′UTR

To estimate the miRNAs with potential binding sites to different MICA 3′UTR alleles, we performed a primary search for miRNAs with described binding capacity to all MICA 3′UTRs, according to the SNPs harbored in this region, using miRbase search tool (https://www.mirbase.org/). We also used the miRBD target prediction database to search for additional miRNAs with targets on SNP-containing MICA 3′UTR (Chen and Wang, 2020).

2.3 Bioinformatic analysis of energy and type of interaction between miRNAs and MICA 3′UTR

The strength of the interaction between the selected miRNAs and their target sequence (Gibbs free energy: ΔG) was analyzed using IntaRNA web server (Bernhart et al., 2008; Chitsaz et al., 2009; Mann et al., 2017). The minimal free energy (MFE) for miRNA selection, to cover possible miRNA::RNA interactions, was lower than −8.00 kcal/mol. Only miRNAs interaction with canonical seed sequence features, according to TargetScan, was considered (8-mer, 7mer-A1, 7mer-m8) (Agarwal et al., 2015). We also analyzed previously reported miRNAs with MICA 3′UTR binding sequences not containing SNPs using the same search tools. The general strategy for miRNA search and data analysis is illustrated in Figure 1.

FIGURE 1
www.frontiersin.org

FIGURE 1. Workflow of the methodological approach for target prediction of miRNAs with interaction sites in SNPs located in the 3′UTR of MICA. Initially, we defined the 3′UTR of MICA alleles to further search for miRNA target sequences containing SNPs. MICA transcript nucleotide sequence NM_001177519.3 was used for this purpose, and the different 3′UTRs were analyzed using the mirbase miRNA search tools [39–41] (Table 1). We also used the miRBD target prediction database to look for additional miRNAs with targets on SNPs in MICA 3′UTR [42]. Next, we analyzed the strength of the ΔG between the selected miRNAs and their targets; such analysis was carried out using IntaRNA web server [43–45]. We selected those miRNAs with the highest probability of effect (minimal energy profile - ΔG: < −8 kcal/mol); the type of canonical miRNA interaction (8-mer, 7mer-A1, 7mer-m8) was then analyzed using TargetScan [46]. The miRNAs with MICA targets that have been previously reported were further analyzed according to their target sequence with SNPs in MICA, the ΔG (IntaRNA), and the type of miRNA interaction (TargetScan).

2.4 Evaluation of the impact of SNPs on miRNA target sites

To evaluate the impact of SNPs in MICA alleles on miRNA interaction, we used the difference of ΔG between SNP alleles, as previously reported by Manikandan & Kannan Munirajan (Manikandan and Munirajan, 2014).

ΔMFE=MFE of miRNA[mRNA duplexancestral allele]-MFE of miRNA[mRNA duplexderived allele]

In addition, to complement the identification of gain/loss of miRNA binding sites by SNPs, we used the bioinformatic tool miRNASNP (Liu et al., 2021).

2.5 Folding of 3′UTR secondary structure analysis

To predict the probability of 3′UTR RNA folding, as well as its thermodynamic parameters, we used the RNAfold WebServer from ViennaRNA WebServices (Lorenz et al., 2011) and the UNAfold tool, respectively.

2.6 Phylogenetic tree analysis of the MICA 3′UTR types

Coding DNA HLA sequences from IPD-IMGT/HLA (https://www.ebi.ac.uk/ipd/imgt/hla/) were aligned using Clustal Omega EMBL-EBI Webserver applying the Neighbour-joining tree without distance corrections analysis (Madeira et al., 2022).

3 Results

3.1 Analysis of MICA 3′UTR according to alleles and SNPs

To gain a comprehensive understanding of the 3′UTR regions of the most frequent MICA alleles, we conducted an analysis to define the specific characteristics of each allele, taking into account the presence of stop codons and SNPs. The MICA*015 and *017 alleles lack guanine at position 892, the last nt in exon 4. These changes in the genetic sequence result in a modification of the reading frame, leading to the generation of premature stop codons: TAA (position +997) for MICA*015 and TGA (position +1,035) for MICA*017 (Wang et al., 2016) (Figure 2, 3). The rest of the alleles maintain TAG stop codon (position +1,154). In turn, the MICA*008 allele includes a guanine insertion in the transmembrane region at position 952 located in exon 5. Thus, MICA*008 3′UTR harbors a longer sequence, with 327 nts (position +1,000 to +1,327), compared to the 3′UTR from alleles *015 and *017, which present 289 nts (position +1,038 to +1,327); the remaining MICA alleles contained 170 nts (position +1,157 to +1,327). MICA 3′UTR of *015 and *017 alleles, as well as the remaining MICA 3′UTR sequences, comprised nucleotides from positions 1,038 to 1,327 and 1,157 to 1,327, respectively, using as reference the NCBI Reference Sequence NM_001177519.3. Figure 2 summarizes the lengths of the 3′UTR sequences according to MICA alleles.

FIGURE 2
www.frontiersin.org

FIGURE 2. Schematic representation of the gene structure of main MICA alleles studied. MICA gene is composed of six exons; the length of MICA alleles sequences depends on deletions or insertion in exon 5. The MICA∗008 allele includes a guanine insertion in the transmembrane region at position 952 located in exon 5. The MICA *015 and *017 alleles lack guanine at position 892 in exon 4. These modifications change the reading frame, resulting in the generation of premature stop codons TAA or TGA at position +997 and +1,035 of these alleles, respectively. The rest of the alleles maintain TAG stop codon at position +1,154. MICA*008 3′UTR is longer than the 3′UTR from alleles *015 and *017, as well as the rest of MICA alleles analyzed. MICA alleles showed variations in the number of repetitions of the GCT triplet within exon 5, known as short tandem repetitions (STRs) (Table inserted). These triplets encode alanine in the protein, leading to different lengths of this region across MICA alleles.

FIGURE 3
www.frontiersin.org

FIGURE 3. Binding sites of miRNA on different MICA 3′UTRs. The different types of 3′UTRs of MICA alleles show a combination of nine polymorphisms (highlighted in the light grey zone) and aligned in each case (white zone). The presence of polymorphisms leads to differences in the binding of miRNAs. The miRNAs and the binding sites sequences in the MICA 3′UTRs are also detailed (dark grey zone). (A) Hsa-miR-1184 shows a preferential binding site within UTR1a, as this region is a coding sequence (CDS) in the other alleles. Furthermore, hsa-miR-711 and hsa-miR335-5p have a binding site preferentially within UTR1a and UTR2c. (B) hsa-miR-1207-5p shows two different binding sites. Hsa-miR-17-5p, hsa-miR106a-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-106b, and hsa-miR-93 also show multiple binding sites that involve polymorphisms (bold nucleotide). The hsa-miR-520d-3p shows a binding site to a non-polymorphic region. (C) hsa-miR-7156-3p also shows a binding site to a non-polymorphic region, binding to all types of MICA 3′UTRs. gDNA: genomic sequence of DNA.

Considering the variations in the MICA 3′UTR sequences, including differences in the length and the presence of SNPs, we adopted a novel categorization system. We subdivided the 3′UTRs into distinct categories: UTR1a (*008), UTR1b (*004, *009, *049), UTR1c (*010, *019, *027), UTR2a (*002), UTR2b (*007), UTR2c (*015, *017, *068), UTR4 (*001, *012, *018), and UTR7 (*011). UTR3, UTR5, and UTR6, as classified by Luo and cols (Luo et al., 2013), were excluded from our analysis due to the presence of a deletion (in UTR3) or because of their low global allele frequencies for the SNPs rs140390705 (for UTR5) and rs113015830 (for UTR6). The new classification aligns with the clade grouping in the phylogenetic tree analysis (Supplementary Figure S1) based on sequence similarity. This novel UTR classification reflects that the UTR2a, UTR2b, UTR4, and UTR7 belong to the LI lineage, while UTR1a, UTR1b, and UTR1c originate from the LII lineage (Chen and Gyllensten, 2014). Additionally, the 3′UTR RNAs have different conformational ensembles and shift folding dynamics. For instance, UTR1a, UTR1b/c, UTR2a/b, UTR2c, UTR4, and UTR7 have a ΔG of folding at 37°C of −111.40, −47.80, −52.20, −96.30, −48.60, and −48.0 kcal/mol, respectively (Supplementary Table S2). This indicates that the presence of SNPs changes the thermodynamics of the secondary structure folding of 3′UTR RNAs, which may subsequently shift miRNA binding sites.

3.2 Impact of SNPs within the MICA 3′UTRs on the gain and loss of binding of miRNAs

In our analysis, we examined the potential impact of SNPs within MICA 3′UTR on miRNA binding. Using the bioinformatic tool miRNASNP, we identified that the presence of SNPs can result in the gain or loss of binding sites for specific miRNAs, as outlined in Table 2. For instance, the presence of rs9266826 (G/A) generates binding sites to 6 miRNAs with a favorable ΔG value, such as hsa-miR-22-3p (−17.91 kcal/mol) and hsa-miR-4267 (−18.48 kcal/mol); conversely, 5 miRNAs, including hsa-miR-6078 (−21.18 kcal/mol), loose their binding sites to MICA 3′UTR in the presence of the same SNP. On the other hand, the presence of SNP rs9266829 (A/G) does not confer any miRNA binding gain but generates the loss of binding of hsa-miR-636 (−14.36 kcal/mol), hsa-miR-758-3p (−11.31 kcal/mol), and hsa-miR-581 (−13.83 kcal/mol). These findings acquire relevance when considering the dysregulation of miRNAs in cancer, as gaining or losing their binding to MICA transcripts could potentially affect the expression of the protein.

TABLE 2
www.frontiersin.org

TABLE 2. Gain and loss of binding of miRNAs for presence of SNPs in the MICA 3′UTRs.

3.3 Prediction of main putative miRNAs targets on different types of MICA 3′UTRs

We considered different types of MICA 3′UTR to predict potential miRNA targets with a high probability of binding to MICA 3′UTR. These included the specific combinations of SNPs described on Table 1 and Figure 3. Our analysis identified 12 miRNAs with potential binding sites to MICA 3′UTR variants, which exhibited diverse strength of interaction (ΔG ranged from −8.00 to −18.14 kcal/mol). The miRNA that meet these characteristics are hsa-miR-1184, hsa-miR-711, hsa-miR-335-5p, hsa-miR-1207-5p, hsa-miR-20b-5p, hsa-miR-93, hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-520d-3p, hsa-miR-106b, and hsa-miR-7156-3p, according to the nt position to which they bind (Supplementary Table S1). These miRNA fulfilled our predefined criteria, as they bind to their target seed sequences from canonical binding sites with Gibbs free energy values lower than −8.00 kcal/mol. Notable, some of these miRNAs, such as hsa-miR-1207-5p, hsa-miR-17-5p, hsa-miR-106a-5p, hsa-miR-106b, hsa-miR-20a-5p, and hsa-miR-93, showed multiple binding sites on MICA 3′UTR bearing SNPs, thus increasing the possibility of post-transcriptional regulation of MICA. Figure 3 shows a compilation of the DNA sequences encoding the 3′UTR for multiple MICA alleles, highlighting the presence of nine SNPs of interest within this region, along with the corresponding sequences of miRNAs and their seed sequences found in this region. We observed that hsa-miR-1184 exhibited a preferential binding site to the 3′UTR of MICA*008 (MFE = −13.9 kcal/mol due to the presence of a premature codon stop. However, other MICA alleles did not show the same binding capacity, since this site corresponds to their coding sequences. Additionally, a binding site for hsa-miR-711 found on the 3′UTRs of MICA*015 and MICA*017 (UTR2c), which could be affected by the presence of the rs9266825 SNPs, did not affect the 3′UTRs of other MICA alleles. These findings highlight the binding of miRNAs to specific MICA alleles, as different alleles exhibit distinct sequences.

The accessibility of these miRNAs to 3′UTR sequences depends on the unpaired or paired nts involved in the binding in the seed sequence, which are generated from secondary RNA structures of each 3′UTR type, as depicted in Figure 4. For example, UTR1a and UTR1b, as well as UTR2, UTR4, and UTR7 exhibit different RNA folding patterns, which determines the availability of binding sites to miRNA. In the case of hsa-miR-1207-5p, this miRNA could access both, 3′UTR1a and 3′UTR1b, since the nt in position 2 in its seed sequence is unpaired. Conversely, hsa-520d-3p could access 3′UTR1b but not 3′UTR1a, whereas hsa-miR-335-5p showed the opposite pattern of accessibility.

FIGURE 4
www.frontiersin.org

FIGURE 4. Accessibility of miRNAs to 3′UTR RNA fold of UTR1a (MICA*008), UTR1b (MICA*004, *009, *049), and UTR2c (MICA*015, *017, *). The base pair probability of 3′UTR RNA folding is indicated in the colored scale, where a lesser probability is highlighted in blue and a greater one in red. These structures presented folding capacity with the most favorable thermodynamic energy (Supplementary Table S2). The start nucleotide involved in the interaction of miRNAs with each UTR is indicated with arrows. (A) The folding structure of UTR1a is the largest and in consequence the folding is different to other 3′UTR types. The start nucleotide that binds to hsa-miR-1184, hsa-miR- 1207-5p, and hsa-miR-335-5p with UTR1a show that this latter has an unpaired nucleotide in its RNA folding structure, indicating the greater accessibility of these miRNAs to RNA. (B) The folding of UTR1b is similar to the UTR2a, UTR2b, UTR4, and UTR7 (not shown). This structure presents unpaired nucleotides that favor the binding of hsa-miR-1207-5p, hsa-miR-520d-3p, and a group of miRNAs (hsa-miR-17-5p, hsa-miR-106a-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, and hsa-miR-106b). (C) The UTR2c folding has few nucleotides with a high probability of base pairs; however, this structure was the best prediction given by RNAfold for this sequence. The hsa-miR-520d-3p may be the only one with accessibility to this structure, since the start nucleotide involved in the binding is unpaired.

3.4 The presence of SNPs in MICA 3′UTR within the binding site of miRNAs has a direct impact on their interaction strength

The presence of SNPs within the 3′UTR region can disrupt the binding site for certain miRNAs, which affects the strength of their interaction. Among the analyzed SNPs, rs9266829 (A/G) showed the highest potential as a target for miR-106a-5p, miR-106b, miR-17-5p, miR-20a-5p, miR-20b-5p, and miR-93, as shown in Table 3. Furthermore, the strength of interaction between these miRNAs and their respective target sites exhibited allele specificity, with UTR1a and UTR2a presenting higher ΔMFE values, ranging from 1 to 3 (Table 3). The difference in the binding energy between miR-20b-5p and UTR1a or UTR7 is the highest (ΔMFE = 2.9). Notably, miR-106b displayed potential binding exclusively to UTR1 (rs9266829A) in all its variations, including UTR1a, UTR1b, and UTR1c, while the presence of rs9266829G resulted in the loss of binding site for this miRNA. Intriguingly, the secondary structure of these MICA 3′UTRs revealed an unpaired nucleotide in the sequence involved in hsa-miR-106b binding (Figure 4), suggesting enhanced accessibility of the miRNA to its targets.

TABLE 3
www.frontiersin.org

TABLE 3. Impact of SNPs on MICA alleles on miRNA interaction.

The binding site of hsa-miR-711, located within the rs9266825 variant, displayed binding energy only with the UTR1a and UTR2c alleles, resulting in MFE = −18.14 and −16.44, respectively). Although certain miRNA binding sites did not contain any SNPs, the evaluation of their interaction strength with different MICA 3′UTR regions indicated that variations in the binding energy could impact miRNA binding capacity.

3.5 SNPs outside the binding site of miRNAs indirectly influence their interaction strength with MICA 3′UTR

We identified an indirect influence of SNPs in the binding energy of hsa-miR-520d-3p, which exhibits variability depending on the type of 3′UTR. Even though hsa-miR-520d-3p does not bind to UTR1a, the accessibility of this miRNA to the secondary structure of the majority of the 3′UTRs is favorable due to the presence of unpaired nucleotides in the 3′UTR sequence, such as in the case of UTR1b, UTR2c, UTR4 and UTR7. Nevertheless, in the UTR1a, the nucleotide involved in the binding is already paired, thus affecting the binding potential of this miRNA. This observation suggests that the presence of SNPs harbored in the 3′UTR sequences could indirectly influence the binding sites of specific miRNAs due to changes in their RNA folding.

4 Discussion

The 3′UTR regions of MICA alleles display considerable diversity, resulting in variations in the binding energy to their potential miRNA. Consequently, these differences in binding probability can lead to diverse regulatory effects on MICA expression, particularly in the context of miRNA dysregulation in cancer and inflammatory processes. Notably, SNPs located within putative miRNA binding sites in the MICA 3′UTR have been shown to exert significant regulatory influence by affecting miRNA binding, either enhancing or reducing this phenomenon. In our study, we investigated the collective impact of SNPs within the MICA 3′UTR region, specifically focusing on the distinct 3′UTR types observed across different MICA alleles. Through bioinformatic approaches, we aimed to improve the prediction of miRNAs binding potential that may exert functional effects, thereby enhancing our understanding of MICA protein regulation.

The presence of SNPs in the 3′UTR region on genes involved in cancer has revealed that disruptions of the regulatory role of miRNAs can contribute to the development of cancer and autoimmune diseases (Yu et al., 2020; Tchacrome et al., 2022) SNPs in this region have been proposed as biomarkers of prognosis, prevention (Jeon et al., 2021), susceptibility and chemotherapy response (Cipollini et al., 2014), or even associated with cancer risk (Yu et al., 2020). The presence of multiple SNPs within this region explains, in part, the observed variability in cancer-related gene expression among patients. For instance, in non-small cell lung cancer (NSCLC), the rs9266825 and rs1882 MICA variants, which were included in our bioinformatic analysis, have shown a significant association with patient overall survival, suggesting these alleles a potential prognostic marker for NSCLC clinical outcomes (Xu et al., 2014). Additionally, MICA alleles have been associated with higher susceptibility to inflammatory diseases, mainly autoimmune diseases, such as ankylosing spondylitis (MICA *007) (Zhou et al., 2014), Behçet’s syndrome (MICA*009 and/or *019) (Hughes et al., 2005; Munoz-Saa et al., 2006), ulcerative colitis (MICA *007), and inflammatory bowel disease with type 2 peripheral arthropathy (MICA*008) (Orchard et al., 2001). As a whole MICA alleles may contribute to the pathophysiology and the different of miRNAs binding sites offered by these multiple alleles may impact on MICA regulation.

Hence, the presence of SNPs within MICA 3′UTR may prevent or generate binding sites to miRNA, resulting in the gain or loss of miRNA regulators of MICA gene expression impacting on cellular performance in the tumoral microenvironment, induction of damage, or tissue rejection. In the case of gaining or improved interaction strength to a miRNA binding site, these SNPs may promote MICA transcript degradation, resulting in reduced expression of the protein on the cell membrane. In cancer, reduced MICA expression diminishes the recognition of NK cells via the NKG2D activating receptor, thereby avoiding the antitumoral immune response (Xing and Ferrari de Andrade, 2020). As an example of this possible outcome, our analysis showed that UTR4 have increased ΔG versus UTR1a/b/c for miR-20b-5p (Table 3), a miRNA that is upregulated in hepatocellular cancer (Zhang et al., 2014), suggesting that the presence of UTR4 could enhance down-expression of MICA by miR-20b-5p in this type of cancer; however, function assays are necessary to elucidate this MICA regulatory mechanism. On the other hand, the loss of binding sites for miRNAs may disrupt the control exerted by particular miRNA on MICA transcription. This disruption could lead to enhanced MICA expression, consequently inducing the activation of NK cells via the NKG2D receptor. This scenario has been found in patients with autoimmune disease, such as systemic lupus erythematosus (Spada et al., 2015); the magnitude of this response may rely on the specific MICA variant, which may possibly result in increased tissue damage.

In the tumor environment, certain MICA variants that express methionine in the 129-amino acid (MICA-129 Met), such as the allelic variants *001, *002, *007, *011, *012, *015, *017, and *018 have been reported to exhibit a higher affinity for NKG2D receptor and potentially lead to a strong cytotoxic response (Isernhagen et al., 2016). However, MICA expression should be viewed with caution, because of cleavage susceptibility by metalloproteases expressed by tumor cells (Waldhauer et al., 2008; Chitadze et al., 2013). This cleavage results in the release of a soluble form of MICA, which has the capacity to downregulate the NKG2D receptor, thereby reducing the cytotoxic response of NK cells (Groh et al., 2002). Furthermore, upon translation, MICA*008 variant is also released via exosomes, which exhibits a potent downregulation effect on the NKG2D receptor (Ashiru et al., 2013). Consequently, the pathological effects of miRNAs depend on several factors, including MICA3′UTR alleles, miRNA expression profile within each cell type, and the stage of tumor development.

To provide a comprehensive assessment of the SNPs present in each type of MICA 3′UTR on the potential impact on miRNA binding, we grouped all MICA variants described to date into clusters to redefine the types of 3′UTR. This novel classification not only enhances our understanding of their features but also enables the application of our methodology to identify potential targets for miRNA binding to other genes. It is interesting to note the different thermodynamics of the secondary structure of the 3′UTRs, whose characteristics may be related to the presence of SNPs and to the length of their sequences. These changes in the secondary structure and the thermodynamic stability of 3ÚTR could implicate functional consequences, such as changes in protein expression (Chen et al., 2008; Ceolin et al., 2016).

Alternative polyadenylation and alternative splicing are mechanisms that contribute to the 3′UTR diversity, promoting extension or shortening of this region, and influence mRNA regulation by gain or loss of miRNAs binding sites (Han et al., 2018). These alternate 3′UTR isoforms can impact transcriptional regulation, affecting stability, translation and RNA localization (Tian and Manley, 2013). Although our study did not analyze alternative polyadenylation patterns in MICA alleles, future analysis should take these patterns into account, particularly considering the type of cell, as they are often tissue-specific and have potential to influence on biological cellular processes (MacDonald, 2019; Arora et al., 2022).

Herein we propose new and specific miRNAs with a high probability of binding to different 3′ UTR types on MICA alleles applying free-access bioinformatic tools. We identified differential binding sites and binding energies between miRNA and 3′UTR types. We found that hsa-miR-1184 and hsa-miR-335-5p may regulate mainly UTR1a, while hsa-miR-711 may regulate UTR1a and UTR2c; in turn, hsa-miR-106b may regulate all subtypes of UTR1. Also, hsa-miR-1207-5p, hsa-miR-20b-5p, hsa-miR-17-5p, hsa-miR-520d-3p, and hsa-miR-7156-3p showed binding potential to different types of UTRs depending on the energy binding strength. It would be interesting to perform in vitro studies using these miRNAs to analyze their effect on MICA expression, especially in cancer, as MICA plays a crucial role in immune surveillance against this disease. Additionally, we also have identified several miRNAs (hsa-miR-20a-5p, hsa-miR-93, hsa-miR-106a-5p, and hsa-mir-106b) that have been previously reported to be upregulated and act as down-regulators of MICA expression in various types of cancer (Stern-Ginossar et al., 2008; Kishikawa et al., 2013; Wu et al., 2014; Xie et al., 2014; Yang et al., 2015; Jafarzadeh-Samani et al., 2017; Shen et al., 2017; Tang et al., 2019; Shekari et al., 2020; Awad et al., 2021), which supports our prediction approach.

Although the binding site of hsa-miR-1184 is not polymorphic, this miRNA may regulate the expression of MICA in patients carrying UTR1a, whereas its hybridization could be less efficient in other UTRs, since its binding region corresponds to its coding sequences. In addition, hsa-miR-1184 plays a significant role in various types of cancers and its impact on MICA expression should be addressed. In colorectal cancer, the overexpression of this miRNA inhibits cell proliferation and promotes apoptosis by circular RNA hsa_circ_0128,846 (Chen S. et al., 2020), whereas in triple-negative breast cancer (TNBC) subtype tissue and cell lines it is linked to enhanced proliferation, migration and invasion by circular RNA lar hsa_circ_0000732 (Chen et al., 2022). Furthermore, the downregulation of this miRNA promotes tumorigenesis (Wang et al., 2020a) in colorectal cancer through IGFBP2 and ITGA3 genes, and in bladder cancer, it is associated with accelerated tumor progression (Jiang P. et al., 2020; Yang et al., 2020), hsa-miR-1184 has been proposed as a potential serum biomarker to distinguish between benign prostatic hyperplasia (BPH) and prostate cancer (Knyazev et al., 2016). Thus, hsa-miR-1184 is an interesting target for further studies on its regulatory role in cancer-related processes.

Hsa-miR-711, which may regulate MICA expression on patients carrying UTR1a and UTR2c, has been associated with osteosarcoma-associated proliferation, migration, and invasion; it is regulated by circular RNA circ_0008792 (Chen L. et al., 2020). Conversely, in gastric cancer, the expression of hsa-miR-711 is low and its overexpression downregulates adhesion molecules, inhibiting the epithelial-mesenchymal transition (EMT) (Xiao et al., 2018). Similarly, this miRNA has been found in extracellular vesicles (EV), which inhibits the expression of CD44 in bone marrow mesenchymal stem cells in chronic myelogenous leukemia (Jiang Y-H. et al., 2020). The presence of miRNAs in extracellular vesicles suggest that MICA expression on tumor cells can be regulated by miRNA carried into EVs.

Herein we also highlight hsa-miR-1207-5 miRNA, which displayed binding sites for different 3′UTR: two binding sites for UTR1a UTR1b, and UTR1c and one binding site for other MICA UTRs. In both cases, this miRNA could have a favorable accessibility to the secondary structure of 3′UTR. Also, hsa-miR-1207-5p could bind with high strength to UTR1 and UTR4; in turn, it may bind to UTR2 and UTR7 in another region with lower interaction force, thus giving this miRNA a versatile effect depending on the genetic component and the tumor microenvironment. The involvement of hsa-miR-1207-5p in various cancer-related processes has been extensively documented. The evidence shows it is associated with tumor suppressive activity in several types of cancer (Wang and Wu, 2017; Cui et al., 2019; Yan et al., 2020); it has also been regarded as a negative regulator of the expression of several genes involved in epithelial-to-mesenchymal transition (EMT) (Dang et al., 2016; Qin et al., 2016). In colorectal cancer, this miRNA suppresses cell proliferation, migration and invasion through targeting FMNL2, a downstream molecule of a circular RNA related as an oncogene in this type of cancer (Yan et al., 2020); its downregulation has also been associated with shorter overall survival rate in patients (Wang and Wu, 2017). Recently, hsa-miR-1207-5p has been proposed as a biomarker for diagnosis and prognosis of colorectal cancer (Wang et al., 2020b). In lung cancer, it participates in macrophage function control and facilitates the release of proinflammatory chemokines throughout its target colony stimulating factor 1 (CDF1); it can also inhibit STAT3 and AKT signaling (Dang et al., 2016). In gastric cancer, hsa-miR-1207-5p exerts its effects on telomerase reverse transcriptase mRNA (Chen et al., 2014), and promotes proliferation by a long noncoding RNA BC032469 (Lu et al., 2016). In nasopharyngeal carcinoma, hsa-miR-1207 has been shown to be inhibited by non-coding RNA Pvt1 or long non-coding RNA 319, which may act as a sponge to miRNAs, inducing cellular proliferation via PI3K/AKT (Xiao et al., 2018), or carcinogenesis via KLF12 signaling pathways (Nunziato et al., 2019), respectively. In prostate cancer, hsa-miR-1207-5p is downregulated and associates with the androgen receptor via FNDC1, a protein that contains a conserved protein domain of fibronectin and it is upregulated in this type of cancer (Das et al., 2016). In contrast, in breast cancer, this miRNA is upregulated and promotes tumor cell proliferation through STAT6 regulation (Yan et al., 2017). Taken together, these findings suggest the importance of uncovering the targets of miRNA, as well as deciphering their regulatory potential over gene expression, including MICA. The identification of the expression levels of miRNAs in each type of cancer, the presence of SNPs in miRNAs (Ryan et al., 2010), and their effect on MICA expression could give us a major comprehension of their biological effects in this context. Such molecular features are worth exploring using in vitro studies in order to move beyond the limitations of bioinformatic tools alone (Rojo Arias and Busskamp, 2019).

The effect of regulation of gene expression by miRNAs binding to 3′UTR has been widely studied, so has their regulation through coding sequences (CDS), although with less pronounced effects (Friedrich et al., 2020). This is in agreement with a number of reports analyzing CDS-located miRNA binding sites and may be due to the accessibility of the target site in the context of the translation process (Gu et al., 2009). In this work, we focused on studying solely the regulation of MICA 3′UTR by miRNA targets; however, it will be valuable to use bioinformatic tools that include the coding sequences of the genes (Liu et al., 2013). The integration of CDS into bioinformatic tools may improve the efficiency and accuracy of binding prediction relative to algorithm trained with 3′UTR information (Bertolazzi et al., 2020; Lin et al., 2023).

While our initial bioinformatics analysis focused on elucidating the potential interactions between miRNAs and MICA 3′UTR allelic sequences, we acknowledge that this computational approach does not directly unveil the concrete impacts of the selected miRNAs. In this sense, we recognize the imperative need to bridge this gap and carry out experimental approaches to substantiate and delineate the actual effects of miRNAs in the context of our work.

In conclusion, MICA post-transcriptional regulation by miRNAs should be carefully studied, considering the diversity of MICA 3′UTR types within specific MICA alleles. The presence of SNPs in this region could promote gain or loss of binding sites or, alternatively, significantly impact the binding energy; enhancement or detriment of binding may also have a huge importance in the context of cancer and other diseases. Our in silico approach sets the stage for future investigations in cancer and inflammatory models, where empirical experiments will be pivotal to establish conclusive insights, thus advancing our understanding of the intricate relationship between miRNAs and MICA allelic variability.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

KT-S: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing–original draft, Writing–review and editing. CR: Formal Analysis, Investigation, Visualization, Writing–original draft, Writing–review and editing, Methodology, Validation. IC: Formal Analysis, Investigation, Visualization, Writing–original draft, Writing–review and editing. ST: Formal Analysis, Investigation, Visualization, Writing–original draft, Writing–review and editing. YL: Formal Analysis, Investigation, Visualization, Writing–original draft, Writing–review and editing. CA: Formal Analysis, Investigation, Visualization, Writing–original draft, Writing–review and editing. KD-C: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Validation, Visualization, Writing–original draft, Writing–review and editing, Supervision. MM: Conceptualization, Data curation, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing–original draft, Writing–review and editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by National Agency of Research and Development ANID/FONDECYT GRANTS 1221031 (MCM), 3210367 (KD-C), 3230454 (KT-S), 3220181 (YL); ANID Scholarship 1221031 (IC), 21231772 (ST); ANID/FONDEF IDeA ID20I10106 (MCM) and Anillo Regular de Ciencia y/o Tecnología 2021 ACT210068 (CA).

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.2023.1273296/full#supplementary-material

References

Abruzzese, M. P., Bilotta, M. T., Fionda, C., Zingoni, A., Soriani, A., Vulpis, E., et al. (2016). Inhibition of bromodomain and extra-terminal (BET) proteins increases NKG2D ligand MICA expression and sensitivity to NK cell-mediated cytotoxicity in multiple myeloma cells: role of cMYC-IRF4-miR-125b interplay. J. Hematol. Oncol. 9, 134. doi:10.1186/s13045-016-0362-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Agarwal, V., Bell, G. W., Nam, J.-W., and Bartel, D. P. (2015). Predicting effective microRNA target sites in mammalian mRNAs. eLife 4, e05005. doi:10.7554/eLife.05005

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Abdallah, A., Jahanbani, I., Mehdawi, H., Ali, R. H., Al-Brahim, N., Mojiminiyi, O., et al. (2020). Down-regulation of the human major histocompatibility complex class I chain-related gene A (MICA) and its receptor is mediated by microRNA-146b-5p and is a potential mechanism of immunoediting in papillary thyroid carcinoma. Exp. Mol. Pathology 113, 104379. doi:10.1016/j.yexmp.2020.104379

PubMed Abstract | CrossRef Full Text | Google Scholar

Arora, A., Goering, R., Lo, H. Y. G., Lo, J., Moffatt, C., and Taliaferro, J. M. (2022). The role of alternative polyadenylation in the regulation of subcellular RNA localization. Front. Genet. 12, 818668. doi:10.3389/fgene.2021.818668

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashiru, O., López-Cobo, S., Fernández-Messina, L., Pontes-Quero, S., Pandolfi, R., Reyburn, H. T., et al. (2013). A GPI anchor explains the unique biological features of the common NKG2D-ligand allele MICA*008. Biochem. J. 454, 295–302. doi:10.1042/BJ20130194

PubMed Abstract | CrossRef Full Text | Google Scholar

Awad, A. R., Youness, R. A., Ibrahim, M., Motaal, A. A., El-Askary, H. I., ReemA, A., et al. (2021). An acetylated derivative of vitexin halts MDA-MB-231 cellular progression and improves its immunogenic profile through tuning miR- 20a-MICA/B axis. Nat. Prod. Res. 35, 3126–3130. doi:10.1080/14786419.2019.1686372

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernhart, S. H., Hofacker, I. L., Will, S., Gruber, A. R., and Stadler, P. F. (2008). RNAalifold: improved consensus structure prediction for RNA alignments. BMC Bioinforma. 9, 474. doi:10.1186/1471-2105-9-474

PubMed Abstract | CrossRef Full Text | Google Scholar

Bertolazzi, G., Benos, P. V., Tumminello, M., and Coronnello, C. (2020). An improvement of ComiR algorithm for microRNA target prediction by exploiting coding region sequences of mRNAs. BMC Bioinforma. 21, 201. doi:10.1186/s12859-020-3519-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Borchert, G. M., Lanier, W., and Davidson, B. L. (2006). RNA polymerase III transcribes human microRNAs. Nat. Struct. Mol. Biol. 13, 1097–1101. doi:10.1038/nsmb1167

PubMed Abstract | CrossRef Full Text | Google Scholar

Brennecke, J., Stark, A., Russell, R. B., and Cohen, S. M. (2005). Principles of MicroRNA–target recognition. PLOS Biol. 3, e85. doi:10.1371/journal.pbio.0030085

PubMed Abstract | CrossRef Full Text | Google Scholar

Breunig, C., Pahl, J., Küblbeck, M., Miller, M., Antonelli, D., Erdem, N., et al. (2017). MicroRNA-519a-3p mediates apoptosis resistance in breast cancer cells and their escape from recognition by natural killer cells. Cell Death Dis. 8, e2973. doi:10.1038/cddis.2017.364

PubMed Abstract | CrossRef Full Text | Google Scholar

Cai, Y., Yu, X., Hu, S., and Yu, J. (2009). A brief review on the mechanisms of miRNA regulation. Genomics Proteomics Bioinforma. 7, 147–154. doi:10.1016/S1672-0229(08)60044-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ceolin, L., Romitti, M., Siqueira, D. R., Ferreira, C. V., Scapineli, J. O., Assis-Brazil, B., et al. (2016). Effect of 3′UTR RET variants on RET mRNA secondary structure and disease presentation in medullary thyroid carcinoma. PLOS ONE 11, e0147840. doi:10.1371/journal.pone.0147840

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D., and Gyllensten, U. (2014). MICA polymorphism: biology and importance in cancer. Carcinogenesis 35, 2633–2642. doi:10.1093/carcin/bgu215

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, D., Wang, M., Zhang, H., Zhou, S., and Luo, C. (2022). Estrogen receptor β2 (ERβ2)-mediated upregulation of hsa_circ_0000732 promotes tumor progression via sponging microRNA-1184 in triple-negative breast cancer (TNBC). Inflamm. Res. 71, 255–266. doi:10.1007/s00011-021-01536-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Lü, M.-H., Zhang, D., Hao, N.-B., Fan, Y.-H., Wu, Y.-Y., et al. (2014). miR-1207-5p and miR-1266 suppress gastric cancer growth and invasion by targeting telomerase reverse transcriptase. Cell Death Dis. 5, e1034. doi:10.1038/cddis.2013.553

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, L., Shan, Y., Zhang, H., Wang, H., and Chen, Y. (2020b). Up-regulation of Hsa_circ_0008792 inhibits osteosarcoma cell invasion and migration and promotes apoptosis by regulating hsa-miR-711/ZFP1. Onco Targets Ther. 13, 2173–2181. doi:10.2147/OTT.S239256

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, R. W., Bemis, L. T., Amato, C. M., Myint, H., Tran, H., Birks, D. K., et al. (2008). Truncation in CCND1 mRNA alters miR-16-1 regulation in mantle cell lymphoma. Blood 112, 822–829. doi:10.1182/blood-2008-03-142182

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Wang, Y., Xu, M., Zhang, L., Su, Y., Wang, B., et al. (2020a). miR-1184 regulates the proliferation and apoptosis of colon cancer cells via targeting CSNK2A1. Mol. Cell. Probes 53, 101625. doi:10.1016/j.mcp.2020.101625

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., and Wang, X. (2020). miRDB: an online database for prediction of functional microRNA targets. Nucleic Acids Res. 48, D127–D131. doi:10.1093/nar/gkz757

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitadze, G., Lettau, M., Bhat, J., Wesch, D., Steinle, A., Fürst, D., et al. (2013). Shedding of endogenous MHC class I-related chain molecules A and B from different human tumor entities: heterogeneous involvement of the “a disintegrin and metalloproteases” 10 and 17. Int. J. Cancer 133, 1557–1566. doi:10.1002/ijc.28174

PubMed Abstract | CrossRef Full Text | Google Scholar

Chitsaz, H., Salari, R., Sahinalp, S. C., and Backofen, R. (2009). A partition function algorithm for interacting nucleic acid strands. Bioinformatics 25, i365–i373. doi:10.1093/bioinformatics/btp212

PubMed Abstract | CrossRef Full Text | Google Scholar

Cipollini, M., Landi, S., and Gemignani, F. (2014). MicroRNA binding site polymorphisms as biomarkers in cancer management and research. PGPM 7, 173–191. doi:10.2147/PGPM.S61693

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, M., Chang, Y., Fang, Q.-G., Du, W., Wu, J.-F., Wang, J.-H., et al. (2019). Non-coding RNA Pvt1 promotes cancer stem cell–like traits in nasopharyngeal cancer via inhibiting miR-1207. Pathol. Oncol. Res. 25, 1411–1422. doi:10.1007/s12253-018-0453-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Dang, W., Qin, Z., Fan, S., Wen, Q., Lu, Y., Wang, J., et al. (2016). miR-1207-5p suppresses lung cancer growth and metastasis by targeting CSF1. Oncotarget 7, 32421–32432. doi:10.18632/oncotarget.8718

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, D. K., Naidoo, M., Ilboudo, A., Park, J. Y., Ali, T., Krampis, K., et al. (2016). miR-1207-3p regulates the androgen receptor in prostate cancer via FNDC1/fibronectin. Exp. Cell Res. 348, 190–200. doi:10.1016/j.yexcr.2016.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, W., Ma, Y., Zhu, W., Pu, W., Zhang, J., Qian, F., et al. (2020). MICA ∗012:01 allele facilitates the metastasis of KRAS-mutant colorectal cancer. Front. Genet. 11, 511. doi:10.3389/fgene.2020.00511

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedrich, M., Vaxevanis, C. K., Biehl, K., Mueller, A., and Seliger, B. (2020). Targeting the coding sequence: opposing roles in regulating classical and non-classical MHC class I molecules by miR-16 and miR-744. J. Immunother. Cancer 8, e000396. doi:10.1136/jitc-2019-000396

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuertes, M. B., Domaica, C. I., and Zwirner, N. W. (2021). Leveraging NKG2D ligands in immuno-oncology. Front. Immunol. 12, 713158. doi:10.3389/fimmu.2021.713158

PubMed Abstract | CrossRef Full Text | Google Scholar

Gene: MICA (ENSG00000233051) (2023). Variant table - Homo_sapiens - Ensembl genome browser 110. Available at: https://www.ensembl.org/Homo_sapiens/Gene/Variation_Gene/Table?db=core;g=ENSG00000233051;r=HSCHR6_MHC_COX_CTG1:2880151-2896610 (Accessed July 1, 2023).

Google Scholar

Groh, V., Wu, J., Yee, C., and Spies, T. (2002). Tumour-derived soluble MIC ligands impair expression of NKG2D and T-cell activation. Nature 419, 734–738. doi:10.1038/nature01112

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, S., Jin, L., Zhang, F., Sarnow, P., and Kay, M. A. (2009). Biological basis for restriction of microRNA targets to the 3′ untranslated region in mammalian mRNAs. Nat. Struct. Mol. Biol. 16, 144–150. doi:10.1038/nsmb.1552

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, S., Kim, D., Shivakumar, M., Lee, Y.-J., Garg, T., Miller, J. E., et al. (2018). The effects of alternative splicing on miRNA binding sites in bladder cancer. PLOS ONE 13, e0190708. doi:10.1371/journal.pone.0190708

PubMed Abstract | CrossRef Full Text | Google Scholar

Hughes, E. H., Collins, R. W. M., Kondeatis, E., Wallace, G. R., Graham, E. M., Vaughan, R. W., et al. (2005). Associations of major histocompatibility complex class I chain-related molecule polymorphisms with Behcet’s disease in Caucasian patients. Tissue Antigens 66, 195–199. doi:10.1111/j.1399-0039.2005.00465.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huntzinger, E., and Izaurralde, E. (2011). Gene silencing by microRNAs: contributions of translational repression and mRNA decay. Nat. Rev. Genet. 12, 99–110. doi:10.1038/nrg2936

PubMed Abstract | CrossRef Full Text | Google Scholar

IPD-IMGT/HLA (2023). Database. Available at: https://www.ebi.ac.uk/ipd/imgt/hla/(Accessed July 1, 2023).

Google Scholar

Isernhagen, A., Malzahn, D., Bickeböller, H., and Dressel, R. (2016). Impact of the MICA-129Met/val dimorphism on nkg2d-mediated biological functions and disease risks. Front. Immunol. 7, 588. doi:10.3389/fimmu.2016.00588

PubMed Abstract | CrossRef Full Text | Google Scholar

Jafarzadeh-Samani, Z., Sohrabi, S., Shirmohammadi, K., Effatpanah, H., Yadegarazari, R., and Saidijam, M. (2017). Evaluation of miR-22 and miR-20a as diagnostic biomarkers for gastric cancer. Chin. Clin. Oncol. 6, 16. doi:10.21037/cco.2017.03.01

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeon, Y.-J., Cho, S.-H., Kim, E.-J., Ryu, C.-S., Park, H.-S., Kim, J.-W., et al. (2021). 3′-UTR polymorphisms in thymidylate synthase with colorectal cancer prevalence and prognosis. J. Personalized Med. 11, 537. doi:10.3390/jpm11060537

CrossRef Full Text | Google Scholar

Jiang, P., Zhu, Y., Xu, Z., Chen, X., and Xie, L. (2020a). Interference with circBC048201 inhibits the proliferation, migration, and invasion of bladder cancer cells through the miR-1184/ITGA3 axis. Mol. Cell Biochem. 474, 83–94. doi:10.1007/s11010-020-03835-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y.-H., Liu, J., Lin, J., Li, S.-Q., Xu, Y.-M., Min, Q.-H., et al. (2020b). K562 cell-derived exosomes suppress the adhesive function of bone marrow mesenchymal stem cells via delivery of miR-711. Biochem. Biophys. Res. Commun. 521, 584–589. doi:10.1016/j.bbrc.2019.10.096

PubMed Abstract | CrossRef Full Text | Google Scholar

Kishikawa, T., Otsuka, M., Yoshikawa, T., Ohno, M., Takata, A., Shibata, C., et al. (2013). Regulation of the expression of the liver cancer susceptibility gene MICA by microRNAs. Sci. Rep. 3, 2739. doi:10.1038/srep02739

PubMed Abstract | CrossRef Full Text | Google Scholar

Klussmeier, A., Massalski, C., Putke, K., Schäfer, G., Sauter, J., Schefzyk, D., et al. (2020). High-throughput MICA/B genotyping of over two million samples: Workflow and allele frequencies. Front. Immunol. 11, 314. doi:10.3389/fimmu.2020.00314

PubMed Abstract | CrossRef Full Text | Google Scholar

Knyazev, E. N., Fomicheva, K. A., Mikhailenko, D. S., Nyushko, K. M., Samatov, T. R., Alekseev, B. Y., et al. (2016). Plasma levels of hsa-miR-619-5p and hsa-miR-1184 differ in prostatic benign hyperplasia and cancer. Bull. Exp. Biol. Med. 161, 108–111. doi:10.1007/s10517-016-3357-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lefranc, M.-P., Pommié, C., Ruiz, M., Giudicelli, V., Foulquier, E., Truong, L., et al. (2003). IMGT unique numbering for immunoglobulin and T cell receptor variable domains and Ig superfamily V-like domains. Dev. Comp. Immunol. 27, 55–77. doi:10.1016/s0145-305x(02)00039-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P., Willie, S. T., Bauer, S., Morris, D. L., Spies, T., and Strong, R. K. (1999). Crystal structure of the MHC class I homolog MIC-A, a gammadelta T cell ligand. Immunity 10, 577–584. doi:10.1016/s1074-7613(00)80057-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, B. C., Katneni, U., Jankowska, K. I., Meyer, D., and Kimchi-Sarfaty, C. (2023). In silico methods for predicting functional synonymous variants. Genome Biol. 24, 126. doi:10.1186/s13059-023-02966-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Mallick, B., Long, D., Rennie, W. A., Wolenc, A., Carmack, C. S., et al. (2013). CLIP-based prediction of mammalian microRNA binding sites. Nucleic Acids Res. 41, e138. doi:10.1093/nar/gkt435

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C.-J., Fu, X., Xia, M., Zhang, Q., Gu, Z., and Guo, A.-Y. (2021). miRNASNP-v3: a comprehensive database for SNPs and disease-related variations in miRNAs and miRNA targets. Nucleic Acids Res. 49, D1276–D1281. doi:10.1093/nar/gkaa783

PubMed Abstract | CrossRef Full Text | Google Scholar

Lorenz, R., Bernhart, S. H., Höner zu Siederdissen, C., Tafer, H., Flamm, C., Stadler, P. F., et al. (2011). ViennaRNA package 2.0. Algorithms Mol. Biol. 6, 26. doi:10.1186/1748-7188-6-26

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, M.-H., Tang, B., Zeng, S., Hu, C.-J., Xie, R., Wu, Y.-Y., et al. (2016). Long noncoding RNA BC032469, a novel competing endogenous RNA, upregulates hTERT expression by sponging miR-1207-5p and promotes proliferation in gastric cancer. Oncogene 35, 3524–3534. doi:10.1038/onc.2015.413

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, D., Campillo, J. A., López-Hernández, R., Martínez-García, P., López-Sánchez, M., Botella, C., et al. (2008). Allelic diversity of MICA gene and MICA/HLA-B haplotypic variation in a population of the Murcia region in southeastern Spain. Hum. Immunol. 69, 655–660. doi:10.1016/j.humimm.2008.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, J., Tian, W., Liu, X. X., Yu, J., Li, L., and Pan, F. (2013). The genetic architecture of 3’untranslated region of the MICA gene: polymorphisms and haplotypes. Hum. Immunol. 74, 1321–1325. doi:10.1016/j.humimm.2013.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

MacDonald, C. C. (2019). Tissue-specific mechanisms of alternative polyadenylation: testis, brain, and beyond (2018 update). WIREs RNA 10, e1526. doi:10.1002/wrna.1526

PubMed Abstract | CrossRef Full Text | Google Scholar

Madeira, F., Pearce, M., Tivey, A. R. N., Basutkar, P., Lee, J., Edbali, O., et al. (2022). Search and sequence analysis tools services from EMBL-EBI in 2022. Nucleic Acids Res. 50, W276–W279. doi:10.1093/nar/gkac240

PubMed Abstract | CrossRef Full Text | Google Scholar

Mah, A. Y., and Cooper, M. A. (2016). Metabolic regulation of natural killer cell IFN-γ production. CRI 36, 131–147. doi:10.1615/CritRevImmunol.2016017387

CrossRef Full Text | Google Scholar

Manikandan, M., and Munirajan, A. K. (2014). Single nucleotide polymorphisms in MicroRNA binding sites of oncogenes: implications in cancer and pharmacogenomics. OMICS A J. Integr. Biol. 18, 142–154. doi:10.1089/omi.2013.0098

PubMed Abstract | CrossRef Full Text | Google Scholar

Mann, M., Wright, P. R., and Backofen, R. (2017). IntaRNA 2.0: enhanced and customizable prediction of RNA–RNA interactions. Nucleic Acids Res. 45, W435-W439–W439. doi:10.1093/nar/gkx279

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, F. J., Amode, M. R., Aneja, A., Austine-Orimoloye, O., Azov, A. G., Barnes, I., et al. (2023). Ensembl 2023. Nucleic Acids Res. 51, D933–D941. doi:10.1093/nar/gkac958

PubMed Abstract | CrossRef Full Text | Google Scholar

Min, D., Lv, X.-b., Wang, X., Zhang, B., Meng, W., Yu, F., et al. (2013). Downregulation of miR-302c and miR-520c by 1,25(OH)2D3 treatment enhances the susceptibility of tumour cells to natural killer cell-mediated cytotoxicity. Br. J. Cancer 109, 723–730. doi:10.1038/bjc.2013.337

PubMed Abstract | CrossRef Full Text | Google Scholar

Munoz-Saa, I., Cambra, A., Pallarés, L., Espinosa, G., Juan, A., Pujalte, F., et al. (2006). Allelic diversity and affinity variants of MICA are imbalanced in Spanish patients with Behçet’s disease. Scand. J. Immunol. 64, 77–82. doi:10.1111/j.1365-3083.2006.01780.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nunziato, M., Esposito, M. V., Starnone, F., Diroma, M. A., Calabrese, A., Del Monaco, V., et al. (2019). A multi-gene panel beyond BRCA1/BRCA2 to identify new breast cancer-predisposing mutations by a picodroplet PCR followed by a next-generation sequencing strategy: a pilot study. Anal. Chim. Acta 1046, 154–162. doi:10.1016/j.aca.2018.09.032

PubMed Abstract | CrossRef Full Text | Google Scholar

OBrien, J., Hayder, H., Zayed, Y., and Peng, C. (2018). Overview of MicroRNA biogenesis, mechanisms of actions, and circulation. Front. Endocrinol. 9, 402. doi:10.3389/fendo.2018.00402

CrossRef Full Text | Google Scholar

Orchard, T. R., Dhar, A., Simmons, J. D., Vaughan, R., Welsh, K. I., and Jewell, D. P. (2001). MHC class I chain-like gene A (MICA) and its associations with inflammatory bowel disease and peripheral arthropathy. Clin. Exp. Immunol. 126, 437–440. doi:10.1046/j.1365-2249.2001.01662.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Qin, Z., He, W., Tang, J., Ye, Q., Dang, W., Lu, Y., et al. (2016). MicroRNAs provide feedback regulation of epithelial-mesenchymal transition induced by growth factors. J. Cell. Physiology 231, 120–129. doi:10.1002/jcp.25060

PubMed Abstract | CrossRef Full Text | Google Scholar

Rojo Arias, J. E., and Busskamp, V. (2019). Challenges in microRNAs’ targetome prediction and validation. Neural Regen. Res. 14, 1672–1677. doi:10.4103/1673-5374.257514

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryan, B. M., Robles, A. I., and Harris, C. C. (2010). Genetic variation in microRNA networks: the implications for cancer research. Nat. Rev. Cancer 10, 389–402. doi:10.1038/nrc2867

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, V., and Shah, J. (2020). Recent trends in targeting miRNAs for cancer therapy. J. Pharm. Pharmacol. 72, 1732–1749. doi:10.1111/jphp.13351

PubMed Abstract | CrossRef Full Text | Google Scholar

Shekari, N., Javadian, M., Ghaffari, S., Baradaran, B., Darabi, M., and Kazemi, T. (2020). DHA abolishes the detrimental effect of docetaxel on downregulation of the MICA via decreasing the expression level of MicroRNA-20a in gastric cancer. J. Gastrointest. Canc 51, 545–551. doi:10.1007/s12029-019-00280-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Shen, J., Pan, J., Du, C., Si, W., Yao, M., Xu, L., et al. (2017). Silencing NKG2D ligand-targeting miRNAs enhances natural killer cell-mediated cytotoxicity in breast cancer. Cell Death Dis. 8, e2740. doi:10.1038/cddis.2017.158

PubMed Abstract | CrossRef Full Text | Google Scholar

Spada, R., Rojas, J. M., Pérez-Yagüe, S., Mulens, V., Cannata-Ortiz, P., Bragado, R., et al. (2015). NKG2D ligand overexpression in lupus nephritis correlates with increased NK cell activity and differentiation in kidneys but not in the periphery. J. Leukoc. Biol. 97, 583–598. doi:10.1189/jlb.4A0714-326R

PubMed Abstract | CrossRef Full Text | Google Scholar

Stern-Ginossar, N., Gur, C., Biton, M., Horwitz, E., Elboim, M., Stanietsky, N., et al. (2008). Human microRNAs regulate stress-induced immune responses mediated by the receptor NKG2D. Nat. Immunol. 9, 1065–1073. doi:10.1038/ni.1642

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, S., Fu, H., Xu, Q., and Zhou, Y. (2019). miR-20a regulates sensitivity of colorectal cancer cells to NK cells by targeting MICA. Biosci. Rep. 39, BSR20180695. doi:10.1042/BSR20180695

PubMed Abstract | CrossRef Full Text | Google Scholar

Tchacrome, I., Zhu, Q., Saleh, M. A., and Zou, Y. (2022). Diseases association with the polymorphic major histocompatibility complex class I related chain a: MICA gene. Transpl. Immunol. 75, 101665. doi:10.1016/j.trim.2022.101665

PubMed Abstract | CrossRef Full Text | Google Scholar

The Allele Frequency Net Database (2022). Allele, haplotype and genotype frequencies in Worldwide Populations. Available at: http://www.allelefrequencies.net/ (Accessed October 31, 2022).

Google Scholar

Tian, B., and Manley, J. L. (2013). Alternative cleavage and polyadenylation: the long and short of it. Trends Biochem. Sci. 38, 312–320. doi:10.1016/j.tibs.2013.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Toledo-Stuardo, K., Ribeiro, C. H., Canals, A., Morales, M., Gárate, V., Rodríguez-Siza, J., et al. (2021). Major histocompatibility complex class I-related chain A (MICA) allelic variants associate with susceptibility and prognosis of gastric cancer. Front. Immunol. 12, 645528. doi:10.3389/fimmu.2021.645528

PubMed Abstract | CrossRef Full Text | Google Scholar

Vitiani, L. r., Potolicchio, I., D’Amato, M., Baricordi, O. r., and Sorrentino, R. (1998). MICA exon 5 microsatellite typing by DNA heteroduplex analysis: a new polymorphism in the transmembrane region. Tissue Antigens 51, 309–311. doi:10.1111/j.1399-0039.1998.tb03108.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Waldhauer, I., Goehlsdorf, D., Gieseke, F., Weinschenk, T., Wittenbrink, M., Ludwig, A., et al. (2008). Tumor-associated MICA is shed by ADAM proteases. Cancer Res. 68, 6368–6376. doi:10.1158/0008-5472.CAN-07-6768

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, W. Y., Tian, W., Zhu, F. m., Liu, X. X., Li, L. X., and Wang, F. (2016). MICA, MICB polymorphisms and linkage disequilibrium with HLA-B in a Chinese Mongolian population. Scand. J. Immunol. 83, 456–462. doi:10.1111/sji.12437

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. (2014). Composition of seed sequence is a major determinant of microRNA targeting patterns. Bioinformatics 30, 1377–1383. doi:10.1093/bioinformatics/btu045

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Chen, Y., Liu, W., Liu, T., and Sun, D. (2020a). Hsa_circ_0128846 promotes tumorigenesis of colorectal cancer by sponging hsa-miR-1184 and releasing AJUBA and inactivating Hippo/YAP signalling. J. Cell. Mol. Med. 24, 9908–9924. doi:10.1111/jcmm.15590

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Li, L., Xiao, W., and Xu, Q. (2020b). Plasma microRNA-1207-5p as a potential biomarker for diagnosis and prognosis of colorectal cancer. Clin. Lab. 66. doi:10.7754/Clin.Lab.2020.191269

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., and Wu, X. (2017). The role of MicroRNA-1207-5p in colorectal cancer. Clin. Lab. 63, 1875–1882. doi:10.7754/Clin.Lab.2017.170625

PubMed Abstract | CrossRef Full Text | Google Scholar

Wenda, S., Faé, I., Sanchez-Mazas, A., Nunes, J. M., Mayr, W. R., and Fischer, G. F. (2013). The distribution of MICA alleles in an Austrian population: evidence for increasing polymorphism. Hum. Immunol. 74, 1295–1299. doi:10.1016/j.humimm.2013.06.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Zhang, X.-J., Shi, K.-Q., Chen, Y.-P., Ren, Y.-F., Song, Y.-J., et al. (2014). Hepatitis B surface antigen inhibits MICA and MICB expression via induction of cellular miRNAs in hepatocellular carcinoma cells. Carcinogenesis 35, 155–163. doi:10.1093/carcin/bgt268

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiao, W.-S., Li, D.-F., Tang, Y.-P., Chen, Y.-Z., Deng, W.-B., Chen, J., et al. (2018). Inhibition of epithelial-mesenchymal transition in gastric cancer cells by miR-711-mediated downregulation of CD44 expression. Oncol. Rep. 40, 2844–2853. doi:10.3892/or.2018.6681

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, J., Liu, M., Li, Y., Nie, Y., Mi, Q., and Zhao, S. (2014). Ovarian tumor-associated microRNA-20a decreases natural killer cell cytotoxicity by downregulating MICA/B expression. Cell Mol. Immunol. 11, 495–502. doi:10.1038/cmi.2014.30

PubMed Abstract | CrossRef Full Text | Google Scholar

Xing, S., and Ferrari de Andrade, L. (2020). NKG2D and MICA/B shedding: a ‘tag game’ between NK cells and malignant cells. Clin. Transl. Immunol. 9, e1230. doi:10.1002/cti2.1230

CrossRef Full Text | Google Scholar

Xu, J., Tian, S., Yin, Z., Wu, S., Liu, L., Qian, Y., et al. (2014). MicroRNA-binding site SNPs in deregulated genes are associated with clinical outcome of non-small cell lung cancer. Lung Cancer 85, 442–448. doi:10.1016/j.lungcan.2014.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamakawa, R. H., Saito, P. K., Gelmini, G. F., Silva, J. S. da, Bicalho, M. da G., and Borelli, S. D. (2017). MICA diversity and linkage disequilibrium with HLA-B alleles in renal-transplant candidates in southern Brazil. PLOS ONE 12, e0176072. doi:10.1371/journal.pone.0176072

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, C., Chen, Y., Kong, W., Fu, L., Liu, Y., Yao, Q., et al. (2017). PVT1-derived miR-1207-5p promotes breast cancer cell growth by targeting STAT6. Cancer Sci. 108, 868–876. doi:10.1111/cas.13212

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, Y., Su, M., and Qin, B. (2020). CircHIPK3 promotes colorectal cancer cells proliferation and metastasis via modulating of miR-1207-5p/FMNL2 signal. Biochem. Biophys. Res. Commun. 524, 839–846. doi:10.1016/j.bbrc.2020.01.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Qian, H., Fang, Z., Xu, A., Zhao, S., Liu, B., et al. (2020). Silencing circular RNA VANGL1 inhibits progression of bladder cancer by regulating miR-1184/IGFBP2 axis. Cancer Med. 9, 700–710. doi:10.1002/cam4.2650

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, H., Lan, P., Hou, Z., Guan, Y., Zhang, J., Xu, W., et al. (2015). Histone deacetylase inhibitor SAHA epigenetically regulates miR-17-92 cluster and MCM7 to upregulate MICA expression in hepatoma. Br. J. Cancer 112, 112–121. doi:10.1038/bjc.2014.547

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, A., Zhang, J., Mei, Y., Zhong, H., Chen, S., and Song, Q. (2020). Correlation between single nucleotide polymorphisms of an miRNA binding site in the 3’UTR of PTEN and risk of cervical cancer among the han Chinese. Genet. Test. Mol. Biomarkers 24, 381–389. doi:10.1089/gtmb.2019.0269

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Xu, Z., Zhou, X., Zhang, H., Yang, N., Wu, Y., et al. (2014). Loss of expression of MHC class I-related chain A (MICA) is a frequent event and predicts poor survival in patients with hepatocellular carcinoma. Int. J. Clin. Exp. Pathol. 7, 3123–3131.

PubMed Abstract | Google Scholar

Zhou, X., Wang, J., Zou, H., Ward, M. M., Weisman, M. H., Espitia, M. G., et al. (2014). MICA, a gene contributing strong susceptibility to ankylosing spondylitis. Ann. Rheum. Dis. 73, 1552–1557. doi:10.1136/annrheumdis-2013-203352

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: MICA, SNP, miRNA, post-transcriptional regulation, allelic variation, 3′UTR

Citation: Toledo-Stuardo K, Ribeiro CH, Campos I, Tello S, Latorre Y, Altamirano C, Dubois-Camacho K and Molina MC (2023) Impact of MICA 3′UTR allelic variability on miRNA binding prediction, a bioinformatic approach. Front. Genet. 14:1273296. doi: 10.3389/fgene.2023.1273296

Received: 05 August 2023; Accepted: 13 November 2023;
Published: 07 December 2023.

Edited by:

Nicholas Delihas, Stony Brook Medicine, United States

Reviewed by:

Eric Ramirez-Salazar, Wistar Institute, United States
Rana A. Youness, German International University in Cairo (GIU), Egypt

Copyright © 2023 Toledo-Stuardo, Ribeiro, Campos, Tello, Latorre, Altamirano, Dubois-Camacho and Molina. 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: Maria Carmen Molina, mcmolinas@uchile.cl; Karen Dubois-Camacho, kdubois@med.uchile.cl

These authors share senior authorship

Disclaimer: 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.