2. Differentially expressed tsRNAs identified in mice subcutaneous adipose tissue
To further explore the differentially expressed tsRNAs in the four groups, we performed pairwise comparisons among the samples, and the identification criteria of differential expression were a fold change (FC) ≥ 1.5 or ≤ 0.67 (p < 0.05). As was shown in Table S4, a total of 37 up-regulated tsRNAs and 18 tsRNAs were identified in HFVDD vs HFVDS (Fig. 1A). 34 up-regulated tsRNAs and 50 down-regulated tsRNAs were identified in HFVDD vs ConVDD (Fig. 1B). 55 up-regulated tsRNAs and 101 down-regulated tsRNAs were identified in HFVDD vs ConVDS (Fig. 1C). 84 up-regulated tsRNAs and 131 down-regulated tsRNAs were identified in HFVDS vs ConVDD (Fig. 1D). 54 up-regulated tsRNAs and 76 down-regulated tsRNAs were identified in HFVDS vs ConVDS (Fig. 1E). 18 up-regulated tsRNAs and 25 down-regulated tsRNAs were identified in ConVDD vs ConVDS (Fig. 1F).
3. Between-group comparisons of the differentially expressed tsRNAs in VD deficiency mice with or without obesity by high-throughput sequencing technology
We listed the top 10 tsRNAs with the most significantly differential expression in each group according to the FCs (Table 1). 5 tsRNAs (tRF5-20-HisGTG-3, tRF5-22-CysGCA-27, tRF5-20-CysGCA-27, 5'tiRNA-36-ValAAC-3, tRF5-18-CysGCA-27) were highly expressed and 10 tsRNAs (mt-5'tiRNA-33-AlaTGC, mt-5'tiRNA-33-SerTGA, mt-5'tiRNA-32-SerTGA and mt-tRF3a-ProTGG, mt-5'tiRNA-38-LeuTAA, mt-5'tiRNA-31-GluTTC, mt-5'tiRNA-36-LeuTAA, mt-tRF3b-TrpTCA, mt-5'tiRNA-37-LeuTAA, mt-5'tiRNA-32-TrpTCA) were low expressed not only in HFVDD vs HFVDS, but aslo in HFVDD vs ConVDS. 5 tsRNAs (tRF3a-GlyGCC-1, 3'tiRNA-50-SerGCT-1, tRF3b-SerGCT-1, tRF3a-ThrAGT-5, 3'tiRNA-41-GlyGCC-6) were highly expressed only in HFVDD vs ConVDS.
Table 1
Part of the significantly differentially expressed tsRNAs in each group (sorted by FC)
Group
|
tsRNA ID
|
type
|
tsRNA-sequence
|
FC
|
HFVDD vs HFVDS (up regulated)
|
|
|
|
|
tRF5-20-HisGTG-3
|
5'tRF
|
GCCGTGATCGTATAGGGGTT
|
3.36
|
|
tRF3a-GlyGCC-1
|
3'tRF
|
TTCCCGGCCCATGCACCA
|
2.80
|
|
tRF5-22-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGGTAGA
|
2.51
|
|
3'tiRNA-50-SerGCT-1
|
3'tiRNA
|
TAATCCATTGTGCTCTGCACGCGTGGGTTCGAATCCCACCTTCGTCGCCA
|
2.48
|
|
tRF3b-SerGCT-1
|
3'tRF
|
TCGAATCCCACCTTCGTCGCCA
|
2.46
|
|
tRF5-20-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGGTA
|
2.46
|
|
tRF3a-ThrAGT-5
|
3'tRF
|
ATCCCAGCAAGGCCTCCA
|
2.28
|
|
5'tiRNA-36-ValAAC-3
|
5'tiRNA
|
GTTTCCGTAGTGTAGTGGTCATCACGCTCGCCTAAC
|
2.23
|
|
3'tiRNA-41-GlyGCC-6
|
3'tiRNA
|
CATGAGGGAGGCCCAGGTTCAATTCCAGGCCCATTGCACCA
|
2.20
|
|
tRF5-18-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGG
|
2.14
|
HFVDD vs HFVDS (down regulated)
|
|
|
|
|
mt-5'tiRNA-33-AlaTGC
|
5'tiRNA
|
GAGGTCTTAGCTTAATTAAAGCAATTGATTTGC
|
0.32
|
|
mt-5'tiRNA-33-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTGA
|
0.39
|
|
mt-5'tiRNA-32-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTG
|
0.50
|
|
mt-tRF3a-ProTGG
|
3'tRF
|
GCTCCTTCTTCTTGACCA
|
0.50
|
|
mt-5'tiRNA-38-LeuTAA
|
5'tiRNA
|
ATTAGGGTGGCAGAGCCAGGAAATTGCGTAAGACTTAA
|
0.53
|
|
mt-5'tiRNA-31-GluTTC
|
5'tiRNA
|
GTTTCTGTAGTTGAATTACAACGATGATTTT
|
0.53
|
|
mt-5'tiRNA-36-LeuTAA
|
5'tiRNA
|
ATTAGGGTGGCAGAGCCAGGAAATTGCGTAAGACTT
|
0.53
|
|
mt-tRF3b-TrpTCA
|
3'tRF
|
CACACAAGTTTAACTTCTGCCA
|
0.55
|
|
mt-5'tiRNA-37-LeuTAA
|
5'tiRNA
|
ATTAGGGTGGCAGAGCCAGGAAATTGCGTAAGACTTA
|
0.55
|
|
mt-5'tiRNA-32-TrpTCA
|
5'tiRNA
|
AGAAGTTTAGGATATACTAGTCCGCGAGCCTT
|
0.55
|
HFVDD vs ConVDD (up regulated)
|
|
|
|
|
mt-tRF5-22-LeuTAA
|
5'tRF
|
ATTAGGGTGGCAGAGCCAGGAA
|
1.83
|
|
tRF5-20-SerCGA-2
|
5'tRF
|
GTCACGGTGGCCGAGTGGTT
|
1.81
|
|
5'tiRNA-34-GlyCCC-1
|
5'tiRNA
|
GCGCCGCTGGTGTAGTGGTATCATGCAAGATTCC
|
1.79
|
|
5'tiRNA-35-AsnGTT-1
|
5'tiRNA
|
GTCTCCGTGGCGCAATCGGTCAGCGCGTTCGGCTG
|
1.78
|
|
5'tiRNA-37-LeuTAG-1
|
5'tiRNA
|
GGTAGCGTGGCCGAGCGGTCTAAGGCGCTGGATTTAG
|
1.76
|
|
tRF5-18-ArgTCT-2
|
5'tRF
|
GGCTCTGTGGCGCAATGG
|
1.76
|
|
tRF5-20-HisGTG-2
|
5'tRF
|
GCCGTGATCGTATAGTGGTT
|
1.75
|
|
5'tiRNA-35-AsnGTT-2
|
5'tiRNA
|
GTCTCTGTGGCGCAATTGGTTAGCGCGTTCGGCTG
|
1.75
|
|
tRF5-22-HisGTG-2
|
5'tRF
|
GCCGTGATCGTATAGTGGTTAG
|
1.75
|
|
tRF5-18-ProAGG-3
|
5'tRF
|
AGCTCATTAGTCTAGGGG
|
1.71
|
HFVDD vs ConVDD (down regulated)
|
|
|
|
|
tRF5-18-iMetCAT-1
|
5'tRF
|
AGCAGAGTGGCGCAGCGG
|
0.32
|
|
tRF5-22-AlaAGC-19
|
5'tRF
|
GTTGGAGATTTAGCTCAGCGGT
|
0.35
|
|
tRF5-18-TrpCCA-3
|
5'tRF
|
GGCCTCGTGGCGCAACGG
|
0.43
|
|
mt-tRF3a-TyrGTA
|
3'tRF
|
ATCCTCTTTTTACCACCA
|
0.46
|
|
mt-tRF3b-TyrGTA
|
3'tRF
|
TTAAATCCTCTTTTTACCACCA
|
0.48
|
|
mt-3'tiRNA-39-HisGTG
|
3'tiRNA
|
GAATCTGACAACAGGAAATAAACCTCCTTATTCACCCCA
|
0.48
|
|
mt-tRF3a-IleGAT
|
3'tRF
|
GCCCTCTTATTTCTACCA
|
0.49
|
|
mt-tRF3b-IleGAT
|
3'tRF
|
TCAAGCCCTCTTATTTCTACCA
|
0.50
|
|
tRF5-18-AlaTGC-5
|
5'tRF
|
GGGGATGTAGCTCAGCGG
|
0.50
|
|
tRF5-18-ValCAC-5
|
5'tRF
|
GTTTTTGTAGTGTAGCGG
|
0.50
|
HFVDD vs ConVDS (up regulated)
|
|
|
|
|
tRF5-20-HisGTG-3
|
5'tRF
|
GCCGTGATCGTATAGGGGTT
|
5.07
|
|
tRF5-22-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGGTAGA
|
2.89
|
|
tRF5-20-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGGTA
|
2.88
|
|
tRF5-20-LysTTT-1
|
5'tRF
|
GCCCGGATAGCTCAGTCGGT
|
2.66
|
|
5'tiRNA-34-ValAAC-3
|
5'tiRNA
|
GTTTCCGTAGTGTAGTGGTCATCACGCTCGCCTA
|
2.65
|
|
tRF5-18-CysGCA-27
|
5'tRF
|
GCGGGTATAGCTCAGGGG
|
2.57
|
|
tRF5-18-CysGCA-17
|
5'tRF
|
TGGGGTATAGCTCAGGGG
|
2.51
|
|
5'tiRNA-35-ValAAC-3
|
5'tiRNA
|
GTTTCCGTAGTGTAGTGGTCATCACGCTCGCCTAA
|
2.34
|
|
tRF5-18-HisGTG-2
|
5'tRF
|
GCCGTGATCGTATAGTGG
|
2.33
|
|
tRF5-22-CysGCA-17
|
5'tiRNA
|
GTTTCCGTAGTGTAGTGGTCATCACGCTCGCCTCA
|
2.25
|
HFVDD vs ConVDS (down regulated)
|
|
|
|
|
mt-5'tiRNA-33-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTGA
|
0.11
|
|
mt-5'tiRNA-32-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTG
|
0.15
|
|
mt-5'tiRNA-33-AlaTGC
|
5'tiRNA
|
GAGGTCTTAGCTTAATTAAAGCAATTGATTTGC
|
0.21
|
|
mt-tRF3a-ProTGG
|
3'tRF
|
GCTCCTTCTTCTTGACCA
|
0.22
|
|
mt-tRF3b-IleGAT
|
3'tRF
|
TCAAGCCCTCTTATTTCTACCA
|
0.25
|
|
mt-tRF3a-AspGTC
|
3'tRF
|
AATCTATATATCTTACCA
|
0.27
|
|
mt-tRF3b-ProTGG
|
3'tRF
|
AGTAGCTCCTTCTTCTTGACCA
|
0.27
|
|
mt-tRF3a-IleGAT
|
3'tRF
|
GCCCTCTTATTTCTACCA
|
0.28
|
|
mt-tRF3a-TyrGTA
|
3'tRF
|
ATCCTCTTTTTACCACCA
|
0.29
|
|
mt-5'tiRNA-30-CysGCA
|
5'tiRNA
|
GGTCTTAAGGTGATATTCATGTCGAATTGC
|
0.29
|
HFVDS vs ConVDD (up regulated)
|
|
|
|
|
5'tiRNA-37-LeuTAG-1
|
5'tiRNA
|
GGTAGCGTGGCCGAGCGGTCTAAGGCGCTGGATTTAG
|
2.25
|
|
5'tiRNA-34-AlaAGC-4
|
5'tiRNA
|
GGGGGATTAGCTCAAATGGTAGAGCGCTCGCTTA
|
2.19
|
|
5'tiRNA-34-AlaAGC-6
|
5'tiRNA
|
GGGGAATTAGCTCAAATGGTAGAGCGCTCGCTTA
|
2.15
|
|
5'tiRNA-33-ProTGG-1
|
5'tiRNA
|
GGCTCGTTGGTCTAGGGGTATGATTCTCGGTTT
|
2.14
|
|
5'tiRNA-33-ProAGG-2
|
5'tiRNA
|
GGCTTGTTGGTCTAGGGGTATGATTCTCACTTA
|
2.09
|
|
5'tiRNA-34-GlyCCC-1
|
5'tiRNA
|
GCGCCGCTGGTGTAGTGGTATCATGCAAGATTCC
|
2.08
|
|
tRF5-31-GlyGCC-1
|
5'tRF
|
GCATGGGTGGTTCAGTGGTAGAATTCTCGCC
|
2.03
|
|
5'tiRNA-33-ProTGG-3
|
5'tiRNA
|
GGTTTGTTGGTCTAGTGGTATGATTCTCAGTTT
|
2.02
|
|
5'tiRNA-33-ProTGG-2
|
5'tiRNA
|
GGCTCGTTGGTCTAGGGGTATGATTCTCGCTTT
|
1.97
|
|
tRF5-18-AlaCGC-5
|
5'tRF
|
TCCCTGGTAGTCTAGTGG
|
1.96
|
HFVDS vs ConVDD (down regulated)
|
|
|
|
|
tRF5-18-iMetCAT-1
|
5'tRF
|
AGCAGAGTGGCGCAGCGG
|
0.28
|
|
tRF3a-MetCAT-4
|
3'tRF
|
GCCTCAGAGAAGGCACCA
|
0.29
|
|
tRF3a-GlyGCC-1
|
3'tRF
|
TTCCCGGCCCATGCACCA
|
0.30
|
|
tRF3a-ThrAGT-5
|
3'tRF
|
ATCCCAGCAAGGCCTCCA
|
0.30
|
|
tRF5-22-AlaAGC-19
|
5'tRF
|
GTTGGAGATTTAGCTCAGCGGT
|
0.31
|
|
tRF5-20-LysCTT-18
|
5'tRF
|
ACCCAGCTAGTTCAGTCGGT
|
0.32
|
|
tRF3a-LysTTT-3
|
3'tRF
|
GTCCCTGTTCAGGCGCCA
|
0.32
|
|
3'tiRNA-41-GlyGCC-6
|
3'tiRNA
|
CATGAGGGAGGCCCAGGTTCAATTCCAGGCCCATTGCACCA
|
0.35
|
|
tRF3a-AlaAGC-1
|
3'tRF
|
TCCCCAGCACCTCCACCA
|
0.36
|
|
tRF3b-MetCAT-4
|
3'tRF
|
TCGAGCCTCAGAGAAGGCACCA
|
0.38
|
HFVDS vs ConVDS (up regulated)
|
|
|
|
|
tRF5-20-AlaAGC-6
|
5'tRF
|
GGGGGTATAGCTCAGTGG
|
2.73
|
|
tRF5-20-AlaAGC-10
|
5'tRF
|
GTAGTCGTGGCCGAGTGGTTAAGG
|
2.43
|
|
5'tiRNA-33-ProTGG-4
|
5'tRF
|
TCCCTGGTGGTCTAGTGGCT
|
2.33
|
|
tRF5-18-AlaCGC-5
|
5'tiRNA
|
GGCTCGTTGGTCTAGGGGTATGATTCTCACTTT
|
2.29
|
|
tRF5-18-AlaAGC-10
|
5'tRF
|
GCATTGGTAGTTCAATGGTAGA
|
2.23
|
|
5'tiRNA-33-ProAGG-2
|
5'tiRNA
|
GACCTCGTGGCACAATGGTAGCACGTCTGACTC
|
2.11
|
|
5'tiRNA-34-ArgACG-1
|
5'tRF
|
GGCTCGTTGGTCTAGGGGTATG
|
2.02
|
|
tRF5-22-ArgCCG-3
|
5'tRF
|
GCCCGGATAGCTCAGTCGGT
|
2.00
|
|
tRF5-20-ArgCCG-3
|
5'tiRNA
|
GGGGATGTAGCTCAGCGGTAGAGCACATGCTTT
|
1.97
|
|
tRF5-20-LysTTT-1
|
5'tiRNA
|
GGCCGGTTAGCTCAGTTGGTTAGAGCATGGTGCTAA
|
1.95
|
HFVDS vs ConVDS (down regulated)
|
|
|
|
|
mt-5'tiRNA-33-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTGA
|
0.28
|
|
mt-5'tiRNA-32-SerTGA
|
5'tiRNA
|
GAGAAAGACATATAGGATATGAGATTGGCTTG
|
0.29
|
|
tRF3a-ThrAGT-5
|
3'tRF
|
ATCCCAGCAAGGCCTCCA
|
0.37
|
|
tRF3a-GlyGCC-1
|
3'tRF
|
TTCCCGGCCCATGCACCA
|
0.39
|
|
mt-tRF3b-IleGAT
|
3'tRF
|
TCAAGCCCTCTTATTTCTACCA
|
0.40
|
|
mt-tRF3a-IleGAT
|
3'tRF
|
GCCCTCTTATTTCTACCA
|
0.44
|
|
mt-tRF3b-ProTGG
|
3'tRF
|
AGTAGCTCCTTCTTCTTGACCA
|
0.44
|
|
mt-tRF3a-ProTGG
|
3'tRF
|
GCTCCTTCTTCTTGACCA
|
0.44
|
|
mt-3'tiRNA-39-HisGTG
|
3'tiRNA
|
GAATCTGACAACAGGAAATAAACCTCCTTATTCACCCCA
|
0.45
|
|
mt-tRF3a-AspGTC
|
3'tRF
|
AATCTATATATCTTACCA
|
0.46
|
ConVDD vs ConVDS (up regulated)
|
|
|
|
|
tRF5-18-iMetCAT-1
|
5'tRF
|
AGCAGAGTGGCGCAGCGG
|
3.41
|
|
tRF5-20-LysCTT-18
|
5'tRF
|
ACCCAGCTAGTTCAGTCGGT
|
2.38
|
|
tRF5-18-TrpCCA-3
|
5'tRF
|
GGCCTCGTGGCGCAACGG
|
2.22
|
|
tRF5-18-AlaTGC-5
|
5'tRF
|
GGGGATGTAGCTCAGCGG
|
1.96
|
|
tRF5-18-ValCAC-5
|
5'tRF
|
GTTTTTGTAGTGTAGCGG
|
1.94
|
|
tRF3b-LeuTAA-5
|
3'tRF
|
TCAAAGCTTGGACGAGTCCCCA
|
1.85
|
|
tRF5-20-LysTTT-2
|
5'tRF
|
GCCTGGATAGCTCAGTCGGT
|
1.79
|
|
3'tiRNA-41-MetCAT-4
|
3'tiRNA
|
TAATCTGAAGGTCCTGAGTTCGAGCCTCAGAGAAGGCACCA
|
1.74
|
|
3'tiRNA-41-LysCTT-7
|
3'tiRNA
|
TAATCCCAGGGTCATGGGTTCGAGCCCCATATTAGGCACCA
|
1.70
|
|
tRF3a-LeuTAA-5
|
3'tRF
|
AGCTTGGACGAGTCCCCA
|
1.68
|
ConVDD vs ConVDS (down regulated)
|
|
|
|
mt-5'tiRNA-30-CysGCA
|
5'tiRNA
|
GGTCTTAAGGTGATATTCATGTCGAATTGC
|
0.38
|
|
5'tiRNA-35-GlyGCC-5
|
5'tiRNA
|
GCATTGGTGGTTCAGTGGTAGAATTCTTGCCTGCC
|
0.39
|
|
mt-5'tiRNA-21-SerGCT
|
5'tiRNA
|
AAGAAAGATTGCAAGAACTGC
|
0.42
|
|
5'tiRNA-33-GlyGCC-1
|
5'tiRNA
|
GCATGGGTGGTTCAGTGGTAGAATTCTCGCCTG
|
0.47
|
|
mt-tRF3b-ProTGG
|
3'tRF
|
AGTAGCTCCTTCTTCTTGACCA
|
0.47
|
|
5'tiRNA-35-GlyGCC-1
|
5'tiRNA
|
GCATGGGTGGTTCAGTGGTAGAATTCTCGCCTGCC
|
0.47
|
|
5'tiRNA-35-GlyGCC-6
|
5'tiRNA
|
GCATGGGTGGTTCAGTGGTAGAATTCTCACCTGCC
|
0.48
|
|
mt-tRF3b-IleGAT
|
3'tRF
|
TCAAGCCCTCTTATTTCTACCA
|
0.51
|
|
5'tiRNA-34-GlyGCC-1
|
5'tiRNA
|
GCATGGGTGGTTCAGTGGTAGAATTCTCGCCTGC
|
0.51
|
|
5'tiRNA-36-GluCTC-6
|
5'tiRNA
|
TCCCTGGCGGCCTAGTGGTTAGGATTCAGTGCTCTC
|
0.51
|
4. Between-group comparisons of the selected differentially expressed tsRNAs in VD deficiency mice with or without obesity verified by RT-PCR
In the next step, in order to validate the differential expression, 2 high candidate tRFs (tRF5-20-HisGTG-3 and tRF5-22-CysGCA-27) and 4 low candidate tRFs (mt-5'tiRNA-33-AlaTGC, mt-5'tiRNA-33-SerTGA, mt-5'tiRNA-32-SerTGA and mt-tRF3a-ProTGG) not only in HFVDD vs HFVDS, but aslo in HFVDD vs ConVDS were measured by RT-PCR. 1 high candidate tRFs (tRF3a-GlyGCC-1) in HFVDD vs ConVDS was measured by RT-PCR. RT-PCR verification further demonstrated that 1 tsRNAs (tRF5-20-HisGTG-3, all P < 0.05) was significantly up-regulated and 1 tsRNA (mt-tRF3a-ProTGG, all P < 0.05) was significantly down-regulated not only in HFVDD vs HFVDS, but aslo in HFVDD vs ConVDS. 1 tsRNAs (tRF5-22-CysGCA-27, P < 0.05) was significantly up-regulated and 3 tsRNA (mt-5'tiRNA-32-SerTGA, mt-5'tiRNA-33-SerTGA and mt-5'tiRNA-33-AlaTGC, all P < 0.05) were significantly down-regulated only in HFVDD vs ConVDS.
4. GO Enrichment, KEGG Pathway Analyses, and Associated Target Genes of the Candidate tsRNAs
The functional annotation of candidate tsRNAs were performed as previously described. The biological processes (BP) terms including cellular process, biological regulation, cellular metabolic process, regulation of biological process, regulation of cellular process, metabolic process, primary metabolic process, organic substance metabolic process, nitrogen compound metabolic process, cellular macromolecule metabolic process were found to highly enriched in the target genes of 3 up-regulated tsRNAs (Fig. 2A). The BP terms including cellular process, biological regulation, regulation of biological process, regulation of cellular process, cellular metabolic process, multicellular organismal process, response to stimulus, metabolic process, primary metabolic process, cell communication were found to highly enriched in the target genes of 4 down-regulated tsRNAs (Fig. 3D). The cellular components (CC) terms including cellular anatomical entity, intracellular, organelle, intracellular organelle, membrane − bounded organelle, intracellular membrane − bounded organelle, cytoplasm, membrane, nucleus, intrinsic component of membrane were found to highly enriched in the target genes of 3 up-regulated tsRNAs (Fig. 3B). The CC terms including cellular anatomical entity, intracellular anatomical structure, organelle, intracellular organelle, membrane, membrane − bounded organelle, intracellular membrane − bounded organelle, cytoplasm, integral component of membrane, intrinsic component of membrane were found to highly enriched in the target genes of 4 down-regulated tsRNAs (Fig. 3E). The molecular functions (MF) terms including binding, protein binding, on binding, heterocyclic compound binding, organic cyclic compound binding, cation binding, metal ion binding, catalytic activity, enzyme binding, identical protein binding were found to highly enriched in the target genes of 3 up-regulated tsRNAs (Fig. 3C). The MF terms including binding, protein binding, ion binding, metal ion binding, cation binding, organic cyclic compound binding, heterocyclic compound binding, catalytic activity, molecular transducer activity, signaling receptor activity were found to highly enriched in the target genes of 4 down-regulated tsRNAs (Fig. 3F).
The top 10 significant pathways of candidate tsRNA were shown in Fig. 4. It was demonstrated that the 3 up-regulated tsRNAs seemed to be most strongly associated with melanoma, breast cancer, chronic myeloid leukemia, FoxO signaling pathway glioma, autophagy-animal, hypertrophic cardiomyopathy, axon guidance, endocrine resistance, PD-L1 expression and PD-1 checkpoint pathway in cancer (Fig. 4A). By contrast, 4 down-regulated tsRNAs seemed to be most strongly associated with olfactory transduction, colorectal cancer, sphingolipid signaling pathway, cushing syndrome, gap junction, wnt signaling pathway, melanogenesis, proteoglycans in cancer, long − term depression, adherens junction(Fig. 4B). .
TargetScan and Miranda were used to predict the target genes of the qPCR-verified DE tsRNAs. There were 403, 330, 4182, 4425, 4324 and 2564 targeted genes of tRF5-20-HisGTG-3, tRF5-22-CysGCA-27, mt-5'tiRNA-33-AlaTGC, mt-5'tiRNA-33-SerTGA, mt-5'tiRNA-32-SerTGA and mt-tRF3a-ProTGG, respectively. The top 20 potential targeted genes of tRF5-20-HisGTG-3, tRF5-22-CysGCA-27, mt-5'tiRNA-33-AlaTGC, mt-5'tiRNA-33-SerTGA, mt-5'tiRNA-32-SerTGA and mt-tRF3a-ProTGG, respectively were listed in Table 2. The co-expression network of the top 20 target genes and differential tsRNA is shown in Fig. 5.
Table 2
The top 20 predicted target genes of tsRNA
tsRNA
|
The top 20 predicted target genes
|
tRF5-20-HisGTG-3
|
Tmem87b, Aco1, Gsto2, Ppp2r5a, Midn, Itsn1, Olfr1157, Ptpn13, Crlf2, Agpat4, Snrk, Sdcbp, Mmab, Agtrap, 2310011J03Rik, Flot2, Ing5, Imp4, 9930012K11Rik, Ubr7
|
tRF5-22-CysGCA-27
|
Akap17b, Chd5, Lonrf2, Grin3a, Ccnj, Ikbip, Rab11b, Abcb9, Prx, Sgsm1, Gpr132, Kdm4a, Ybey, Dnajb6, Rab35, Apc2, Smim24, Syndig1l, Cers3, Tfdp1
|
mt-5_tiRNA-33-AlaTGC
|
Nt5dc3, Cdr1, Dach1, B3galt6, Gm4922, Homer2, Lig3, Zfp59, D2hgdh, Brd8, Gm42669, Ikzf2, Epha4, Vmn1r216, Elf5, Tcl1, Rnf13, Stx19, Slc41a2, Tff3
|
mt-5_tiRNA-33-SerTGA
|
Pou1f1, Olfr2, Sp1, Tnfrsf21, Itgb6, Aqp7, Chd9, Suclg2, Vwa7, Mkx, Cntnap5c, Olfr1151, Sv2c, Zfp560, Kcne1, Stxbp6, Slc35d1, Tox2, Dlgap1, Fn3krp
|
mt-5_tiRNA-32-SerTGA
|
Pou1f1, Olfr2, Tnfrsf21, Itgb6, Sp1, Aqp7, Suclg2, Chd9, Vwa7, Mkx, Cntnap5c, Olfr1151, Kcne1, Stxbp6, Sv2c, Zfp560, Slc35d1, Tox2, Dlgap1, Olfr599
|
mt-tRF3a-ProTGG
|
Cplx2, Emp2, Fam124a, Olfr357, Gcm1, Vmn1r66, Il21r, Pde3b, 1500009C09Rik, Fat3, Atp1a2, Olfr1262, Olfr160, Pigr, Olfr1122, Ms4a12, Slco3a1, Pdzd9, Runx2, Polr3b
|