Transcriptome Analysis of gene in Mouse M-1 cells revealing the functional mechanism of Astragali radix extract

Astragali Radix (AR),the dried root of legumes, belongs to the Qi- invigorating herbs in traditional Chinese medicine and plays an important role in the treatment of many diseases. In order to understand the mechanism of action of AR extract, we used AR extract to treat M-1, mouse kidney cells, and used transcriptome sequencing technology to detect the genomic transcription level of the cells under the action of AR at different concentrations and times. Result The results showed that after a low concentration of AR treatments on the cells, the expression of genes related to cell growth and cellular immune response changed signicantly, among which multiple genes are related to mitochondrial function, while high concentrations of AR affected the expression of histones and disease-related genes. It showed that the low concentration of AR extract can achieve the effect of invigorating Qi by regulating the function of mitochondria. In addition, several important genes and pathways were identied as potential targets of AR activation.

levels [2]. Flavonoids are mainly ingested as glycosides, which have exhibited many potential bene cial activities in animal research and human trials [3]. APS is a type of water-soluble heteropolysaccharide with biological activity, which has the characteristics of low toxic and side effects, no residue and intolerance. The biological functions of APS include regulating immune function, anti-aging, anti-tumor [4], lowering blood sugar, lowering blood fat, anti-brosis, antibacterial, radiation protection and antiviral effects [5]. AR extract was the isolated effective ingredients from AR dry root, and its biological effects was as effective as dry AR root.
In the literature, the existing high-throughput sequencing studies related to AR are either the analysis of the therapeutic effects of the main components of AR on disease models [6,7], or the analysis of the effects of AR extracts on plants development [8]. Except for the transcriptome analysis of compound medicines containing AR [9], there is no transcriptome analysis research on the effect of AR extract to illustrate the healing effect of AR. For this study, we applied the extract of AR, a traditional Chinese medicine for invigorating Qi, to mouse M-1 cells, and used modern transcriptomics methods to screen differently expression genes during cell culture to study the effects of traditional Chinese medicine on cells. It is a bridge between traditional medicine and modern medical research, which helps us to understand the functional mechanism of the Qi-invigorating herb.

Result
Transcriptome analysis of the effect of low-dose AR on M-1 cell In order to reveal the molecular mechanism of the effects of AR on cells after different periods of exposure time, the low-dose AR extract (50 µg/ml ) was added to the cell culture medium, and the cell pellets were collected at two time points of 4 h and 24 h. Use high-throughput sequencing to remove predicted genes (gene names start with GM and then coded with numbers), and compare gene expression changes between the control group and treatment group at the same time point. The results showed that 439 Differently expressed genes (DEGs) (including 233 up-regulated genes and 206 down-regulated genes) and 475 (including 268 up-regulated genes and 207 down-regulated genes) were screened out between 4 h and 24 h after AR treatment (|fold change|>2 and p < 0.05) (Fig. 1, A,B). The result showed that the AR extract made a greater impact on cell growth after AR treating the cells.
So as to understand the mechanism of AR more clearly, we chose to analyze the 64 common differentially expression genes(Co-DEGs) (Fig. 1C). Gene Ontology (GO) analysis classi es and enriches differential genes from three different dimensions: Biological Process(BP), Cellular Component(CC), and Molecular Function (MF), thereby to re ect the dynamic change process between different groups. It was found that there were more enriched genes in biological regulation, metabolic process and cell development in BP term. The cell part, organelle part, membrane and other enriched genes in CC term were more abundant, and the binding, catalytic activity, molecular function regulation and other enriched genes in MF term were more enriched. The most enriched BP terms included response to UV-A, synthesis and regulation of MHC class II biosynthetic processes (Tlr4, Spi1), response to toxic substances, immune response (Spi1, C3, Rnf26, Il2, Tlr4) and protein modi cation. The number of enriched genes related to protein modi cation was the largest among them. Receptor complex (Adrb3, Tlr4, Itga10) and cilium were the most enriched CC terms. The most enriched MF terms were growth factor activity (Cers1, Tgfa, Il2), carbohydrate transmembrane transporter activity (Slc2a4, Aqp3) ( Fig. 1, D-G). It showed that the AR extract had a greater regulatory effect on cell growth and immune response-related genes after AR acting on cells.
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed to explore the most signi cantly enriched DEG pathway (Fig. 1, H). Among the top 20 KEGG pathways, the most enriched were Chagas disease (American trypanosomiasis), NOD-like receptor signaling pathway (Nek7, Antxr1, Tlr4), PI3K-AKT signaling pathway (Il2, Tgfa, Tlr4, Itga10 ), HTLV-I infection(Spi1, Il2, Rasl2-9), and pathways in cancer. NOD-like receptor signaling pathway played a key role in pathogen recognition and innate immune response; PI3K-AKT signaling pathway was involved in the regulation of many basic cellular processes, indicating that the effects of AR extract on cells changed cellular immunerelated genes.
Transcriptome analysis of the effect of high-dose AR on M-1 cell The effects of high-dose AR extract on M-1 cells were studied, in which the concentration of AR was increased to 5 times the original (250 µg/ml). As in the previous experiment, the cells were cultured for 4 hours with AR extract, and then the cells were harvested. By comparing the gene expression changes between the control group and treatment group, DEGs were analyzed to identify the effects of highconcentration AR extract on the cells. The experimental results showed that a total of 64 DEGs were detected, of which 35 genes were up-regulated and 29 genes were down-regulated (|fold change|>2 and p < 0.05) (Fig. 2, A). The result revealed that the high-concentration AR extract had an signi cant effect on cell growth after treating the cells.
In order to clarify the effective mechanism of the high concentration of AR extract on the growth of cells, the 64 DEGs that were screened were further analyzed. GO analysis showed that BP was mainly concentrated in DNA, nucleosome assembly, negative regulation of immune system process and T cell migration. CC terms were mainly enriched in the intrinsic component of plasma membrane, integral component of plasma membrane, and intrinsic component of membrane. The MF terms were enriched in gated channel activity, and protein heterodimerization activity (Fig. 2, B-D). It showed that the high concentration of AR extract mainly acted on basic life processes such as DNA assembly and material transmembrane transportation. Changes in these processes were more likely to change the cell growth state.
KEGG pathway enrichment analysis of differential genes was performed, and the results showed that among the top 20 signaling pathways, the most enriched KEGG pathways were related to diseases such as alcohol abuse, systemic lupus erythematosus, renin secretion, viral carcinogenesis, taste transmission, and Shigellosis (Fig. 2, E). It was also found that the genes enriched in the pathways of alcoholism and Systemic lupus erythematosus were all histones ( Table 1). As it was known, histones was involved in the assembly of nucleosomes [10], and histone modi cations were the main type of epigenetic modi cation. Histones were modi ed by a large number of enzymes and related to a variety of cancers [11,12]. The results showed that AR extract affected the gene expression of histones genes and different concentration of AR had different effect effects, which indicated that AR extract affected the assembly of nucleosomes in multiple manners, which caused the different clinic functions of AR. We performed Venn analysis on the DEGs in the same time period of high dose and low dose, and found a total of 17 Co-DEGs (Table 2). 8 co-DEGs (2810468N07Rik, Il27ra, Dnah5, AC123954.1, Pcdhgc4, Gp1bb, Clcn1, Ccl20) were found to have opposite expression levels in the two groups. The genes expression of high dose AR among these co-DEGs were positively correlated with the occurrence of diseases, which demonstrated that the KEGG pathways altered by high-concentration AR extract were related to the occurrence of diseases, and had side effects on cells growth.  Discussion AR is a traditional Chinese medicine for invigorating Qi that is widely used in clinical Chinese medicine, and its clinical e cacy has been universally recognized by modern medicine [13]. Qi is a relatively abstract concept in traditional Chinese medicine. The theory of traditional Chinese medicine believes that Qi is the driving force of human life activities and one of the basic substances that maintain the normal life activities of the human body from different aspects. Studies have shown that AR extract has complex components, including Astragaloside, Astragalus avonoids, APS and amino acids, which play a role in immune regulation, multi-organ protection, hypoglycemia, anti-virus, and anti-tumor [13,14]. As so far, the cellular biological mechanism of AR curative effect is not clear.
In order to understand the AR effects on cells through modern technology, we used AR extract to treat cells for different times and concentrations. With the help of transcriptomics technology, we analyzed the changes in cell gene transcription levels and reveal the signi cance of AR in cell life activities. The results showed that treatment of cells with a low concentration of AR extract can screen out proli c DEGs, and the co-DEGs of different treatment times were mainly enriched in cell growth, metabolic processes, immune response and other related pathways, which were consistent with the multiple clinical function of AR in Chinese medicine [13]. After a high dose of AR acts on cells, its signal pathways are mainly concentrated on disease-related pathways such as systemic lupus erythematosus, viral carcinogenesis, and Shigellosis. The results well explained that the high dose of AR extract may cause some side-effect in healthy people.
The growth and development of cells is usually a complex process involving various regulatory molecules. Cell growth depends on the comprehensive effects of biochemical and molecular processes that synthesize its various components inside the cell. In our previous experimental results, it was found that low dose of AR increased the content of ATP and GSH in cells [15]. ATP is the carrier of energy released by cells during respiration and directly provides energy for various life activities. Mitochondrial homeostasis is important for maintaining mitochondrial function and cell stability [16]. Low concentration of AR acted on the differential genes screened by cells, and multiple genes become potential target genes, such as ceramide synthase 1(Cers1), aquaporin 3 (Aqp3), Interleukin-2 (Il2) and so on (Result 3.1). Cers1 is a transduction stress signal that regulates lipid metabolism from the endoplasmic reticulum to the outer mitochondrial membrane [17]. The transmembrane transport of H 2 O 2 requires the participation of Aqp. H 2 O 2 crosses the cytoplasmic membrane from outside the cell into the cytoplasm with the assistance of Aqp3 [18]. Aqp3 can also inhibit cell apoptosis by reducing the genes expression related Wnt/GSK-3β/β-catenin pathway. Toll-like receptor 4(Tlr4) participates in a variety of autoimmune diseases by acting as a bridge between innate immunity and acquired immunity. The activation of Tlr4 induces mitochondrial dynamic imbalance and promotes the occurrence of diseases [19]. NIMA-related kinase 7 (Nek7) is a multifunctional kinase that affects centrosome replication, mitochondrial regulation, intracellular protein transport, DNA repair and mitotic spindle assembly [20], and participates in NF-κB signaling. (Il2) is a pro-in ammatory cytokine that can regulate cell growth, proliferation, and apoptosis. After Il2 stimulation, the quality of mitochondria changes, or the production of ROS in cells increases [21,22]. In addition, there are some differential gene functions related to mitochondria, such as the perredoxin (Prx) family that controls the level of H 2 O 2 in the cell and regulates signal transduction [23]. Glutathione S-transferase (GST) family members 7 (Gstm7) exerts the detoxi cation effect of electrophilic compounds by binding to reduced glutathione [24]. The transcription of E3 ligase proteins Fbxo15 and Ring nger protein 26 (Rnf26) is involved in the immune regulation function of cells [25,26]. These genes are directly and indirectly related to the function of mitochondria. And it can be seen from Fig. 1, E that in the BP of co-DEGs, protein modi cation related genes had the most enrichment. It showed that low concentrations of AR may play a role in changing the function of mitochondria by modifying the protein after expression of mitochondrial-related genes, thereby improving cellular immune function and cell growth.
The decrease of mitochondrial membrane potential can cause a series of biochemical changes in the mitochondrial membrane, which is an early event in the process of cell apoptosis [27]. High concentration of AR reduces the production of ATP and lowers the membrane potential of mitochondria [15]. Excessive mitochondrial ssion can cause mitochondrial ROS overload and inhibit the production of mitochondrial ATP [28]. As it showed in Table 1, the high concentration of AR extract produced drug toxicity for cell growth and caused cellular immune response. Among the co-DEGs shared with the low-concentration AR extract, the expression levels of multiple genes showed opposite results ( Table 2). The expression patterns of these genes under the action of high concentrations of AR were mostly related to the occurrence of diseases. Studies have found that human neutrophils treated with Interleukin-27 (Il27) swelled the cells, ruptured mitochondria and condensed the nucleus [29]. Chemokine (C-C motif) ligand 20(CCL20) is expressed in various human tissues and immune cells, and involved in local immune cell recruitment [30]. CCL20 is signi cantly up-regulated in a variety of cancer tissues. CCL20 produced by liver cancer cells recruits CCR6 + CD5 + B cells and induces blood vessels generated to promote tumor growth [31]. In our experiments, it was also found that Ccl20 expression was highly expressed under the action of high concentrations of AR. The mutation of Clcn1 is related to myotonia [32]. GPIb-IX-V complex is expressed on megakaryocytes and platelets, which acts as a receptor for von Willebrand factor and mediates platelet-subendothelial interactions, and is a major participant in hemostasis. The glycoprotein Ib (GPIb) -IX-V complex contains four different transmembrane subunits: GPIbα, GPIbβ, GPIX and GPV. Mutations in GPIbβ caused GPIb-IX-V defects [33]. Some experiments have found that high concentrations of AR can shorten the lifespan of bees and also have an inhibitory effect on the growth of cancer cells [34,35]. It was obvious that high concentrations of drugs had adverse effects on cells and may promote the occurrence of diseases. In addition, the high dose of AR extract also changes the expression of multiple histones, which can cause some genes to change their expression patterns at the DNA level. This also shows that the high concentration of AR extract may have drug toxicity on cells, resulting in changes in gene expression patterns related to diseases.
Based on the RNA-Seq analysis, we gained further understanding of the molecular mechanism of AR, and also learned the preliminary functional mechanism for the side effects of the high drug concentration AR.
Several genes were identi ed as potential targets of AR, which were expected to become treatment strategies for diseases caused by abnormal expression of the them. Our future research will focus on animal model to verify these hypotheses.

Conclusion
Using RNA-Seq, the DEGs related to the effects of the AR extract on the cells were identi ed. The GO and KEGG function enrichment analysis of DEGs showed that a low concentration of AR extract was bene cial to cell growth and enhance the expression of genes related to cellular immune response. Many genes are related to mitochondrial function, and high concentrations of AR may cause abnormal expression of multiple disease-related genes in cells. Here, for the rst time, we use AR extract to directly act on healthy cells, and explain the role of AR from the perspective of mitochondria. This results have a guiding role in accurately understanding the pharmacological effects of AR and its molecular mechanism, which is helpful for application of AR extract in disease prevention and treatment.

Extraction of Astragali Radix
Astragalus root of 50 grams (g) was sliced to pieces and powdered. The powder was soaked in 500 ml distilled water for 1 hour (h) and then boiled at 100℃ for 30 minutes. The mixed liquid was ltered to get the supernatant, and the insoluble matter was again mixed with 400 ml distilled water and boiled at 100 ℃ for 20 minutes. The supernatant was obtained by the same way and mixed with the former supernatant. The mixture was concentrated by rotary evaporation apparatus to obtain approximate 7.2 g of extracted product.

Cells culture and treatment
Mouse cell line M1 (ATCC, CRL-2038) was grown in Dulbecco's modi ed Eagle's medium (DMEM) supplemented with 10 % fetal bovine serum (FBS) and 1% penicillin-streptomycin. Cells were cultured at 37°C in a 5% CO 2 incubator. During the experimental period, cells were divided into three groups, which were low dose, high dose and control group (CK), respectively. Low dose extraction of Astragalus (50 µg/ml) was added to the supernatant of cells and incubated for 4 h and 24 h respectively, and then the supernatant was discarded. High dose extraction of Astragalus (250 µg/ml) was added to the supernatant of cells and incubated for 4 h, followed by removing of AR supernatant. Cells with no AR treatment were used as controls, which were processed with the same methods except that same volume of PBS was added instead of AR extraction. The cells were trypsinized, washed with PBS, and the cell pellet was collected and washed with PBS again. The experiment was conducted in three biological replicates.

RNA extraction and construction of cDNA library
Total RNA of M1 cells treated with extraction of Astragalus or mock treated was extracted and isolated using TaKaRa MiniBEST Universal RNA Extraction Kit (9767) according to the manufacturer's instructions. The quantity and purity of RNA was assessed by Agarose gel electrophoresis.
RNA was then sent to Sangon Biotech (Shanghai) Co., Ltd for transcriptome analysis. mRNA was isolated from total RNA and then reverse transcribed into cDNA to construct the cDNA library. Sequences were initial analyzed using Illumina HiSeqTM system, Shanghai Sangon Biotech, Shanghai, China.

Library quality identi cation and sequence comparison analysis
The quality of raw data, sequences obtained from the Illumina HiSeqTM system, were further assessed by FastQC software. Nucleic acid sequences containing low quality base, unknown base, etc. were excluded. The remaining sequences after quality control were compared with the reference genome of Mus musculus. Expression abundance of genes was calculated in terms of transcripts per million (TPM

Consent for publication
Not applicable.

Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Competing interests
The authors have no competing interest to report.

ENN Research Fund
Authors' contributions BW designed the speci c plan of the experiment and the modi cation of the article; QF performed cell culture experiment and RNA extraction; CZ was a major contributor in analysising data and writing the manuscript; YC wrote the rst draft of the article; ML extracted astragalus Radix; YC and LW was responsible for the drug concentration exploration before the experiment; JT was assisted in data analysis; ZC was responsible for the reference review; YZ assisted in RNA extraction; XB provided conceptual contributions, overall design and guidance. All authors read and approved the nal manuscript.