[1] So JTY. Treatment of disease with acupuncture. Paradigm Publications. 1987.
[2] Zhang XQ, Hao WS. Holographic meridian scraping therapy. Foreign Languages Press. 2000.
[3] Tsai PS, Lee PH, Wang MY. Demographics, training, and practice patterns of practitioners of folk medicine in Taiwan: a survey of the Taipei Metropolitan Area. J Altern Complement Med, 2008; 14, 1243−8. doi:  10.1089/acm.2008.0316
[4] Nielsen A. Gua Sha: a traditional technique for modern practice. 2nd ed. Churchill Livingstone. 2012.
[5] Tu WZ, Cheng RD, Hu J, et al. Combination treatment with Gua Sha and Blood-letting causes attenuation of systemic inflammation, activated coagulation, tissue ischemia and injury during heatstroke in rats. Chin J Integr Med, 2015; 21, 610−7. doi:  10.1007/s11655-014-1816-4
[6] Yang M, Zhang HY, Yue RZ, et al. Gua Sha attenuates thermal hyperalgesia and decreases proinflammatory cytokine expression in serum in rats with lumbar disc herniation induced by autologous nucleus pulposus. J Tradit Chin Med, 2018; 38, 698−704. doi:  10.1016/S0254-6272(18)30908-7
[7] Wang YW, Xi ZW, Pu B, et al. Gua sha therapy for chronic low back pain: a protocol for systematic review. Medicine (Baltimore), 2020; 99, e20606. doi:  10.1097/MD.0000000000020606
[8] Xu X, Liao HC, Tian YJ, et al. The effect of scraping therapy on perimenopausal syndrome: a synthesis of evidence. Complement Med Res, 2018; 25, 413−6. doi:  10.1159/000490910
[9] Braun M, Schwickert M, Nielsen A, et al. Effectiveness of traditional Chinese "Gua Sha" therapy in patients with chronic neck pain: a randomized controlled trial. Pain Med, 2011; 12, 362−9. doi:  10.1111/j.1526-4637.2011.01053.x
[10] Lauche R, Wübbeling K, Lüdtke R, et al. Randomized controlled pilot study: pain intensity and pressure pain thresholds in patients with neck and low back pain before and after Traditional East Asian "Gua Sha" therapy. Am J Chin Med, 2012; 40, 905−17. doi:  10.1142/S0192415X1250067X
[11] Saha FJ, Brummer G, Lauche R, et al. Gua Sha therapy for chronic low back pain: a randomized controlled trial. Complement Ther Clin Pract, 2019; 34, 64−9. doi:  10.1016/j.ctcp.2018.11.002
[12] Wang XQ, Duan PB, Zheng M, et al. Effect of scraping therapy on chronic low back pain: a systematic review and meta-analysis of randomized controlled trials. J Manipulative Physiol Ther, 2021; 44, 255−70. doi:  10.1016/j.jmpt.2020.08.003
[13] Xie XL, Lu LQ, Zhou XP, et al. Effect of Gua Sha therapy on patients with diabetic peripheral neuropathy: a randomized controlled trial. Complement Ther Clin Pract, 2019; 35, 348−52. doi:  10.1016/j.ctcp.2019.03.018
[14] Ren Q, Yu XY, Liao FJ, et al. Effects of Gua Sha therapy on perimenopausal syndrome: a systematic review and meta-analysis of randomized controlled trials. Complement Ther Clin Pract, 2018; 31, 268−77. doi:  10.1016/j.ctcp.2018.03.012
[15] Yuen JWM, Tsang WWN, Tse SHM, et al. The effects of Gua sha on symptoms and inflammatory biomarkers associated with chronic low back pain: a randomized active-controlled crossover pilot study in elderly. Complement Ther Med, 2017; 32, 25−32. doi:  10.1016/j.ctim.2017.03.010
[16] Chen TT, Liu NH, Liu JX, et al. Gua Sha, a press-stroke treatment of the skin, boosts the immune response to intradermal vaccination. PeerJ, 2016; 4, e2451. doi:  10.7717/peerj.2451
[17] Zhao L, Chen J, Li Y, et al. The Long-term effect of acupuncture for migraine prophylaxis: a randomized clinical trial. JAMA Intern Med, 2017; 177, 508−15. doi:  10.1001/jamainternmed.2016.9378
[18] Zhao L, Li DH, Zheng H, et al. Acupuncture as adjunctive therapy for chronic stable angina: a randomized clinical trial. JAMA Intern Med, 2019; 179, 1388−97. doi:  10.1001/jamainternmed.2019.2407
[19] Chang SC, Kim DH, Jang EY, et al. Acupuncture attenuates alcohol dependence through activation of endorphinergic input to the nucleus accumbens from the arcuate nucleus. Sci Adv, 2019; 5, eaax1342. doi:  10.1126/sciadv.aax1342
[20] He YH, Guo XF, May BH, et al. Clinical evidence for association of acupuncture and acupressure with improved cancer pain: a systematic review and meta-analysis. JAMA Oncol, 2020; 6, 271−8. doi:  10.1001/jamaoncol.2019.5233
[21] Liu SB, Wang ZF, Su YS, et al. A neuroanatomical basis for electroacupuncture to drive the vagal–adrenal axis. Nature, 2021; 598, 641−5. doi:  10.1038/s41586-021-04001-4
[22] Yin HY, Tang Y, Lu SF, et al. Gene Expression Profiles at Moxibustioned Site (ST36): a microarray analysis. Evid Based Complement Alternat Med, 2013; 2013, 890579.
[23] Garcia-Vivas JM, Galaviz-Hernandez C, Fernandez-Retana J, et al. Transcriptomic profiling of adipose tissue in obese women in response to acupuncture catgut embedding therapy with Moxibustion. J Altern Complement Med, 2016; 22, 658−68. doi:  10.1089/acm.2015.0200
[24] Wang Y, Jiang HL, Meng H, et al. Genome-wide transcriptome analysis of hippocampus in rats indicated that TLR/NLR signaling pathway was involved in the pathogenisis of depressive disorder induced by chronic restraint stress. Brain Res Bull, 2017; 134, 195−204. doi:  10.1016/j.brainresbull.2017.07.021
[25] Wang ZQ, Huang Y, Wang D, et al. Genome-wide regulation of acupuncture and moxibustion on ulcerative colitis rats. Evid Based Complement Alternat Med, 2021; 2021, 9945121.
[26] Li J, Peng C, Lai DJ, et al. PET-CT and RNA sequencing reveal novel targets for acupuncture-induced lowering of blood pressure in spontaneously hypertensive rats. Sci Rep, 2021; 11, 10973. doi:  10.1038/s41598-021-90467-1
[27] Patro R, Duggal G, Love MI, et al. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods, 2017; 14, 417−9. doi:  10.1038/nmeth.4197
[28] Howe KL, Achuthan P, Allen J, et al. Ensembl 2021. Nucleic Acids Res, 2021; 49, D884−91. doi:  10.1093/nar/gkaa942
[29] St Laurent G, Shtokalo D, Dong B, et al. VlincRNAs controlled by retroviral elements are a hallmark of pluripotency and cancer. Genome Biol, 2013; 14, R73. doi:  10.1186/gb-2013-14-7-r73
[30] Caron M, St-Onge P, Drouin S, et al. Very long intergenic non-coding RNA transcripts and expression profiles are associated to specific childhood acute lymphoblastic leukemia subtypes. PloS One, 2018; 13, e0207250. doi:  10.1371/journal.pone.0207250
[31] Qi F, Gao F, Cai Y, et al. Complex age- and cancer-related changes in human blood transcriptome—implications for pan-cancer diagnostics. Front Genet, 2021; 12, 746879. doi:  10.3389/fgene.2021.746879
[32] 请参照https://cran.r-project.org/doc/FAQ/R-FAQ.html#Citing-R

R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. 2022
[33]

Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol, 2010; 11, R106. doi:  10.1186/gb-2010-11-10-r106
[34]

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014; 15, 550. doi:  10.1186/s13059-014-0550-8
[35]

Ritchie ME, Phipson B, Wu D, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res, 2015; 43, e47. doi:  10.1093/nar/gkv007
[36]

McInnes L, Healy J, Saul N, et al. UMAP: uniform manifold approximation and projection. J Open Source Softw, 2018; 3, 861. doi:  10.21105/joss.00861
[37]

Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics, 2010; 26, 139−40. doi:  10.1093/bioinformatics/btp616
[38] 请参照https: //cran. r-project. org/web/packages/psych/citation. html 和https://www.scholars.northwestern.edu/en/publications/psych-procedures-for-personality-and-psychological-research

Revelle W. Psych: procedures for psychological, psychometric, and personality research. Evanston, Illinois; 2022
[39]

Gene Ontology Consortium. The gene ontology resource: enriching a GOld mine. Nucleic Acids Res, 2021; 49, D325−34. doi:  10.1093/nar/gkaa1113
[40]

Jassal B, Matthews L, Viteri G, et al. The reactome pathway knowledgebase. Nucleic Acids Res, 2020; 48, D498−503.
[41]

Yu GC, Wang LG, Han YY, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics J Integr Biol, 2012; 16, 284−7. doi:  10.1089/omi.2011.0118
[42]

Kapranov P, St Laurent G, Raz T, et al. The majority of total nuclear-encoded non-ribosomal RNA in a human cell is 'dark matter' un-annotated RNA. BMC Biol, 2010; 8, 149. doi:  10.1186/1741-7007-8-149
[43]

Daca-Roszak P, Zietkiewicz E. Transcriptome variation in human populations and its potential application in forensics. J Appl Genet, 2019; 60, 319−28. doi:  10.1007/s13353-019-00510-1
[44] 请参照https://www.jove.com/t/62528/three-differential-expression-analysis-methods-for-rna-sequencing

Liu SY, Wang ZT, Zhu RH, et al. Three differential expression analysis methods for rna sequencing: limma, EdgeR, DESeq2. J Vis Exp JoVE, 2021; 175, e62528
[45]

Liu X, Zhao JL, Xue LT, et al. A comparison of transcriptome analysis methods with reference genome. BMC Genomics, 2022; 23, 232. doi:  10.1186/s12864-022-08465-0
[46]

Stupnikov A, McInerney CE, Savage KI, et al. Robustness of differential gene expression analysis of RNA-seq. Comput Struct Biotechnol J, 2021; 19, 3470−81. doi:  10.1016/j.csbj.2021.05.040
[47]

Flanagan K, Li WX, Greenblatt EJ, et al. End-to-end pipeline for differential analysis of pausing in ribosome profiling data. STAR Protoc, 2022; 3, 101605. doi:  10.1016/j.xpro.2022.101605
[48]

McDonald JH. Multiple comparisons. In: McDonald JH. Handbook of Biological Statistics. 3rd ed. Sparky House Publishing. 2014, 254–60.
[49]

Berger TC, Vigeland MD, Hjorthaug HS, et al. Neuronal and glial DNA methylation and gene expression changes in early epileptogenesis. PLOS ONE, 2019; 14, e0226575. doi:  10.1371/journal.pone.0226575
[50]

Terme JM, Sesé B, Millán-Ariño L, et al. Histone H1 variants are differentially expressed and incorporated into chromatin during differentiation and reprogramming to pluripotency. J Biol Chem, 2011; 286, 35347−57. doi:  10.1074/jbc.M111.281923
[51]

Harshman SW, Young NL, Parthun MR, et al. H1 histones: current perspectives and challenges. Nucleic Acids Res, 2013; 41, 9593−609. doi:  10.1093/nar/gkt700
[52]

Yang SM, Kim BJ, Toro LN, et al. H1 linker histone promotes epigenetic silencing by regulating both DNA methylation and histone H3 methylation. Proc Natl Acad Sci USA, 2013; 110, 1708−13. doi:  10.1073/pnas.1213266110
[53]

Garg M, Perumalsamy LR, Shivashankar GV, et al. The linker histone H1.2 is an intermediate in the apoptotic response to cytokine deprivation in T-effectors. Int J Cell Biol, 2014; 2014, 674753.
[54]

Kim JM, Kim K, Punj V, et al. Linker histone H1.2 establishes chromatin compaction and gene silencing through recognition of H3K27me3. Sci Rep, 2015; 5, 16714. doi:  10.1038/srep16714
[55]

Kalashnikova AA, Rogge RA, Hansen JC. Linker histone H1 and protein-protein interactions. Biochim Biophys Acta, 2016; 1859, 455−61. doi:  10.1016/j.bbagrm.2015.10.004
[56]

Bauden M, Kristl T, Sasor A, et al. Histone profiling reveals the H1.3 histone variant as a prognostic biomarker for pancreatic ductal adenocarcinoma. BMC Cancer, 2017; 17, 810. doi:  10.1186/s12885-017-3834-z
[57]

Izzo A, Ziegler-Birling C, Hill PWS, et al. Dynamic changes in H1 subtype composition during epigenetic reprogramming. J Cell Biol, 2017; 216, 3017−28. doi:  10.1083/jcb.201611012
[58]

Li ZM, Li YL, Tang M, et al. Destabilization of linker histone H1.2 is essential for ATM activation and DNA damage repair. Cell Res, 2018; 28, 756−70. doi:  10.1038/s41422-018-0048-0
[59]

Sollberger G, Streeck R, Apel F, et al. Linker histone H1.2 and H1.4 affect the neutrophil lineage determination. eLife, 2020; 9, e52563. doi:  10.7554/eLife.52563
[60]

Kamieniarz K, Izzo A, Dundr M, et al. A dual role of linker histone H1.4 Lys 34 acetylation in transcriptional activation. Genes Dev, 2012; 26, 797−802. doi:  10.1101/gad.182014.111
[61]

Weiss T, Hergeth S, Zeissler U, et al. Histone H1 variant-specific lysine methylation by G9a/KMT1C and Glp1/KMT1D. Epigenetics Chromatin, 2010; 3, 7. doi:  10.1186/1756-8935-3-7
[62]

Konishi A, Shimizu S, Hirota J, et al. Involvement of histone H1.2 in apoptosis Induced by DNA double-strand breaks. Cell, 2003; 114, 673−88. doi:  10.1016/S0092-8674(03)00719-0
[63]

Lindenboim L, Borner C, Stein R. Nuclear proteins acting on mitochondria. Biochim Biophys Acta, 2011; 1813, 584−96. doi:  10.1016/j.bbamcr.2010.11.016
[64]

Cascone A, Bruelle C, Lindholm D, et al. Destabilization of the outer and inner mitochondrial membranes by core and linker histones. PLoS ONE, 2012; 7, e35357. doi:  10.1371/journal.pone.0035357
[65]

Medrzycki M, Zhang YZ, McDonald JF, et al. Profiling of linker histone variants in ovarian cancer. Front Biosci, 2012; 17, 396−406. doi:  10.2741/3934
[66]

Duce JA, Smith DP, Blake RE, et al. Linker histone h1 binds to disease associated amyloid-like fibrils. J Mol Biol, 2006; 361, 493−505. doi:  10.1016/j.jmb.2006.06.038
[67]

Okamura H, Yoshida K, Amorim BR, et al. Histone H1.2 is translocated to mitochondria and associates with bak in bleomycin-induced apoptotic cells. J Cell Biochem, 2008; 103, 1488−96. doi:  10.1002/jcb.21537
[68]

Flex E, Martinelli S, Van Dijck A, et al. Aberrant function of the C-terminal tail of HIST1H1E accelerates cellular senescence and causes premature aging. Am J Hum Genet, 2019; 105, 493−508. doi:  10.1016/j.ajhg.2019.07.007
[69]

Millán-Ariño L, Islam ABMMK, Izquierdo-Bouldstridge A, et al. Mapping of six somatic linker histone H1 variants in human breast cancer cells uncovers specific features of H1.2. Nucleic Acids Res, 2014; 42, 4474−93. doi:  10.1093/nar/gku079
[70]

Jenne CN, Urrutia R, Kubes P. Platelets: bridging hemostasis, inflammation, and immunity. Int J Lab Hematol, 2013; 35, 254−61. doi:  10.1111/ijlh.12084
[71]

Gilthorpe JD, Oozeer F, Nash J, et al. Extracellular histone H1 is neurotoxic and drives a pro-inflammatory response in microglia. F1000Res, 2013; 2, 148. doi:  10.12688/f1000research.2-148.v1