Dorsal Striatum Transcriptome Profile Profound Shift in Repeated Aggression Mouse Model Converged to Networks of 12 Transcription Factors after Fighting Deprivation

Both aggressive and aggression-deprived (AD) species represent pathologic cases intensely addressed in psychiatry and substance abuse disciplines. Previously, we reported that AD mice displayed a higher aggressive behavior score than the aggressive group, implying the manifestation of a withdrawal effect. We employed an animal model of chronic social conflicts, curated in our lab for more than 30 years. In the study, we pursued the task of evaluating key events in the dorsal striatum transcriptome of aggression experienced mice and AD species compared to controls using RNA-Seq profiling. Aggressive species were subjected to repeated social conflict encounters (fights) with regular positive (winners) experience in the course of 20 consecutive days (A20 group). This led to a profoundly shifted transcriptome expression profile relative to the control group, outlined by more than 1000 differentially expressed genes (DEGs). RNA-Seq cluster analysis revealed that elevated cyclic AMP (cAMP) signaling cascade and associated genes comprising 170 differentially expressed genes (DEGs) in aggressive (A20) species were accompanied by a downturn in the majority of other metabolic/signaling gene networks (839 DEGs) via the activation of transcriptional repressor DEGs. Fourteen days of a consecutive fighting deprivation period (AD group) featured the basic restoration of the normal (control) transcriptome expression profile yielding only 62 DEGs against the control. Notably, we observed a network of 12 coordinated DEG Transcription Factor (TF) activators from 62 DEGs in total that were distinctly altered in AD compared to control group, underlining the distinct transcription programs featuring AD group, partly retained from the aggressive encounters and not restored to normal in 14 days. We found circadian clock TFs among them, reported previously as a withdrawal effect factor. We conclude that the aggressive phenotype selection with positive reward effect (winning) manifests an addiction model featuring a distinct opioid-related withdrawal effect in AD group. Along with reporting profound transcriptome alteration in A20 group and gaining some insight on its specifics, we outline specific TF activator gene networks associated with transcriptional repression in affected species compared to controls, outlining Nr1d1 as a primary candidate, thus offering putative therapeutic targets in opioid-induced withdrawal treatment.

Download Full-text

RNA-Seq Analysis Identifies Differentially Expressed Genes in Subcutaneous Adipose Tissue in Qaidaford Cattle, Cattle-Yak, and Angus Cattle

Animals ◽

10.3390/ani9121077 ◽

2019 ◽

Vol 9 (12) ◽

pp. 1077 ◽

Cited By ~ 2

Author(s):

Chengchuang Song ◽

Yongzhen Huang ◽

Zhaoxin Yang ◽

Yulin Ma ◽

Buren Chaogetu ◽

...

Keyword(s):

Adipose Tissue ◽

Meat Quality ◽

Differentially Expressed Genes ◽

Transcriptome Analysis ◽

Cell Metabolism ◽

Differentially Expressed ◽

Receptor Interaction ◽

Control Group ◽

Rna Seq ◽

Angus Cattle

In the beef industry, fat tissue is closely related to meat quality. In this study, high-throughput RNA sequencing was utilized for adipose tissue transcriptome analysis between cattle-yak, Qaidamford cattle, and Angus cattle. The screening and identification of differentially expressed genes (DEGs) between different breeds of cattle would facilitate cattle breeding. Compared to Angus cattle adipose tissue, a total of 4167 DEGs were identified in cattle-yak adipose tissue and 3269 DEGs were identified in Qaidamford cattle adipose tissue. Considering cattle-yak as a control group, 154 DEGs were identified in Qaidamford cattle adipose tissue. GO analysis indicated the significant enrichment of some DEGs related to lipid metabolism. The KEGG pathway database was also used to map DEGs and revealed that most annotated genes were involved in ECM-receptor interaction and the PI3K-Akt signal pathway, which are closely related to cell metabolism. Eight selected DEGs related to adipose tissue development or metabolism were verified by RT-qPCR, indicating the reliability of the RNA-seq data. The results of this comparative transcriptome analysis of adipose tissue and screening DEGs suggest several candidates for further investigations of meat quality in different cattle breeds.

Download Full-text

Effect of hyperoside on cervical cancer cells and transcriptome analysis of differentially expressed genes

Cancer Cell International ◽

10.1186/s12935-019-0953-4 ◽

2019 ◽

Vol 19 (1) ◽

Cited By ~ 4

Author(s):

Weikang Guo ◽

Hui Yu ◽

Lu Zhang ◽

Xiuwei Chen ◽

Yunduo Liu ◽

...

Keyword(s):

Cervical Cancer ◽

Cell Proliferation ◽

Differentially Expressed Genes ◽

Interaction Network ◽

Differentially Expressed ◽

Control Group ◽

Rna Seq ◽

Network Construction ◽

Target Network ◽

Key Genes

Abstract Background Hyperoside (Hy) is a plant-derived quercetin 3-d-galactoside that exhibits inhibitory activities on various tumor types. The objective of the current study was to explore Hy effects on cervical cancer cell proliferation, and to perform a transcriptome analysis of differentially expressed genes. Methods Cervical cancer HeLa and C-33A cells were cultured and the effect of Hy treatment was determined using the Cell Counting Kit-8 (CCK-8) assay. After calculating the IC50 of Hy in HeLa and C-33A cells, the more sensitive to Hy treatment cell type was selected for RNA-Seq. Differentially expressed genes (DEGs) were identified by comparing gene expression between the Hy and control groups. Candidate genes were determined through DEG analysis, protein interaction network (PPI) construction, PPI module analysis, transcription factor (TF) prediction, TF-target network construction, and survival analysis. Finally, the key candidate genes were verified by RT-qPCR and western blot. Results Hy inhibited HeLa and C33A cell proliferation in a dose- and time-dependent manner, as determined by the CCK-8 assay. Treatment of C-33A cells with 2 mM Hy was selected for the subsequent experiments. Compared with the control group, 754 upregulated and 509 downregulated genes were identified after RNA-Seq. After functional enrichment, 74 gene ontology biological processes and 43 Kyoto Encyclopedia of Genes and Genomes pathways were obtained. According to the protein interaction network (PPI), PPI module analysis, TF-target network construction, and survival analysis, the key genes MYC, CNKN1A, PAX2, TFRC, ACOX2, UNC5B, APBA1, PRKACA, PEAR1, COL12A1, CACNA1G, PEAR1, and CCNA2 were detected. RT-qPCR was performed on the key genes, and Western blot was used to verify C-MYC and TFRC. C-MYC and TFRC expressions were lower and higher than the corresponding values in the control group, respectively, in accordance with the results from the RNA-Seq analysis. Conclusion Hy inhibited HeLa and C-33A cell proliferation through C-MYC gene expression reduction in C-33A cells and TFRC regulation. The results of the current study provide a theoretical basis for Hy treatment of cervical cancer.

Download Full-text

Study on the Effect of Macrophages on Vascular Endothelium in Mice With Different TCM Syndromes of Dyslipidemia and its Biological Basis Based on RNA-Seq Technology

Frontiers in Pharmacology ◽

10.3389/fphar.2021.665635 ◽

2021 ◽

Vol 12 ◽

Author(s):

Jing Chen ◽

Chao Ye ◽

Zheng Yang ◽

Tieshan Wang ◽

Bing Xu ◽

...

Keyword(s):

Quality Assessment ◽

Differentially Expressed Genes ◽

Normal Control ◽

Normal Control Group ◽

Differentially Expressed ◽

Control Group ◽

Biological Processes ◽

Rna Seq ◽

Pathogenic Factors ◽

Protective Processes

Background: “Treating the same disease with different methods” is a Traditional Chinese medicine (TCM) therapeutic concept suggesting that, while patients may be diagnosed with the same disease, they may also have different syndromes that require distinct drug administrations.Objective: This study aimed to identify the differentially expressed genes and related biological processes in dyslipidemia in relation to phlegm–dampness retention (PDR) syndrome and spleen and kidney Yang deficiency (SKYD) syndrome using transcriptomic analysis.Methods: Ten ApoE−/− mice were used for the establishment of dyslipidemic disease–syndrome models via multifactor-hybrid modeling, with five in the PDR group and five in the SKYD group. Additionally, five C57BL/6J mice were employed as a normal control group. Test model-quality aortic endothelial macrophages in mice were screened using flow cytometry. Transcriptomic analysis was performed for macrophages using RNA-Seq.Results: A quality assessment of the disease–syndrome model showed that levels of lipids significantly increased in the PDR and SKYD groups, compared to the normal control group, p < 0.05. Applying, in addition, hematoxylin and eosin staining of aorta, the disease model was also successfully established. A quality assessment of the syndrome models showed that mice in the PDR group presented with typical manifestations of PDR syndrome, and mice in the SKYD group had related manifestations of SKYD syndrome, indicating that the syndrome models were successfully constructed as well. After comparing the differentially expressed gene expressions in macrophages of the dyslipidemic mice with different syndromes, 4,142 genes were identified with statistical significance, p < 0.05. Gene ontology analysis for the differentially expressed genes showed that the biological process of difference between the PDR group and the SKYD group included both adverse and protective processes.Conclusion: The differentially expressed genes between PDR syndrome and SKYD syndrome indicate different biological mechanisms between the onsets of the two syndromes. They have distinctive biological processes, including adverse and protective processes that correspond to the invasion of pathogenic factors into the body and the fight of healthy Qi against pathogenic factors, respectively, according to TCM theory. Our results provide biological evidence for the TCM principle of “treating the same disease with different treatments.”

Download Full-text

Complement C3a Not C4a Activates Osteoclasts By Regulating the PI3K/PDK1/SGK3 Pathway in Patients with Multiple Myeloma

Blood ◽

10.1182/blood-2019-127022 ◽

2019 ◽

Vol 134 (Supplement_1) ◽

pp. 4416-4416

Author(s):

Rong Fu ◽

Fengjuan Jiang ◽

Hui Liu ◽

Zhaoyun Liu ◽

Siyang Yan ◽

...

Keyword(s):

Differentially Expressed Genes ◽

Bone Disease ◽

Cathepsin K ◽

Differentially Expressed ◽

Control Group ◽

Rna Seq ◽

Relative Expression ◽

Absorption Area ◽

Osteoclast Resorption

Objective: Myeloma bone disease (MBD) is the most common complication of multiple myeloma (MM). We found that the serum levels of C3/C4 in MM patients were significantly positive correlated with the severity of bone disease in our previously study. Thus, we conduct detailed studies to explore the effect and potential mechanism of C3a/C4a on osteoclasts in patients with MM. Methods: By adding C3a/C4a to culture system of osteoclasts induced from mononuclear cells in vitro, and the expression of related gene, number and function of osteoclasts were detected. RNA-Seq analysis was used to detect the differentially expressed genes on osteoclasts between the complement group and control group, and the possible signal pathways were analyzed. Quantitative real-time PCR (qRT-PCR), Western blot and pathway inhibitors were used for further validation. R esults: In vitro, the osteoclasts area per view induced by C3a 1μg/ml (52.794±13.027 %) and 10μg/ml (51.797±12.464 %) was significantly increased than control group (0μg/ml) (33.668±8.427 %) (P<0.001; P<0.001), the mRNA relative expression of osteoclasts related genes OSCAR/TRAP/RANKL/Cathepsin K induced by C3a 1μg/ml (median 5.041; 3.726; 1.638; 4.752) and 10μg/ml (median 5.140; 3.702; 2.250; 5.172) was significantly increased than control group (0μg/ml) (median 3.137; 2.004; 0.573; 2.257) (1μg/ml P=0.001; P=0.003; P<0.001; P=0.008; 1μg/ml P<0.001; P=0.019; P<0.001; P=0.002), and the absorption area of osteoclast resorption pit per view induced by C3a 1μg/ml (51.464±11.983 %) and 10μg/ml (50.219±12.067 %) was also significantly increased than control group (0μg/ml) (33.845±8.331 %) (P<0.001; P<0.001) in NDMM patients. There was no difference among the osteoclasts area, relative expression of osteoclasts related genes and absorption area of osteoclast resorption pit between C4a (1μg/ml and 10μg/ml) group and control group (0μg/ml). RNA-Seq was performed on total RNA of osteoclasts induced by C3a in 1μg/ml group and 0μg/ml group of 4 patients with NDMM. There were 184 differentially expressed genes that were detected by RNA-Seq analysis. KEGG Pathway enrichment bubble chart shows C3a may through Phosphoinositide 3-kinase (PI3K) signaling pathways (including PI3K-Akt pathway and AKT-independent signaling pathway) promotes the proliferation of osteoclast. Upregulated differentially expressed genes in this pathway among at least 3 patients with sequencing were validated by qRT-PCR and Western Blot. It was found that the relative expression level of Phosphoinositide dependent kinase-1 (PDK1) / Serum and glucocorticoid inducible protein kinases (SGK3) genes (median 2.078; 4.428) in C3a group (1μg/ml) was significantly higher than control group (0μg/ml) (median 1.336; 1.714) (P<0.001; P=0.001). The relative grayscale levels of PDK1/ P-SGK3 protein (1.785±0.323; 2.190±0.274) in C3a group (1μg/ml) was significantly stronger than control group (0μg/ml) (0.8653±0.588; 0.176±0.152) (P=0.034; P<0.001). Under the action of C3a in patients with NDMM, osteoclasts area per view in SGK inhibitor (EMD638683) 1μM group (39.244±9.089 %) and 10μM group (39.299±9.587 %) significantly reduced than control group (0μM) (54.884±12.837 %) (P<0.001; P<0.001), the relative expression of osteoclast related genes OSCAR/RANKL/TRAP/Cathepsin K in EMD638683 1μM group (median 0.869; 1.097; 0.902; 1.328) and 10μM group (median 0.703; 1.391; 0.843; 1.418) significantly decreased than control group (0μM) (median 2.270; 3.024; 2.208; 3.237) (1μM P=0.015; P=0.002; P=0.003; P=0.015; 10μM P=0.012; P=0.006; P<0.001; P=0.017), and the absorption area per view of osteoclast resorption pit in EMD638683 1μM (35.383±7.794 %) group and 10μM group (32.886±8.993 %) significantly reduced than control group (0μM) (49.358±11.856 %) (P < 0.001; P < 0.001). Conclusions:ComplementC3a activates osteoclasts by regulating the PI3K/PDK1/SGK3 pathway in patients with MM, thus leading to the occurrence of bone diseases. SGK inhibitor has a significant inhibitory effect on osteoclasts in patients with MM. This study provide important evidences for the search for new therapeutic targets and strategies for myeloma bone disease patients. Disclosures No relevant conflicts of interest to declare.

Download Full-text

Transcriptome analysis of fasting caecotrophy on hepatic lipid metabolism in New Zealand rabbits

10.1101/493957 ◽

2018 ◽

Author(s):

yadong wang ◽

Huifen Xu ◽

Guirong Sun ◽

Mingming Xue ◽

Shuaijie Sun ◽

...

Keyword(s):

Growth Rate ◽

Lipid Metabolism ◽

Growth Performance ◽

Differentially Expressed Genes ◽

Molecular Mechanisms ◽

Effective Means ◽

Differentially Expressed ◽

Control Group ◽

Rna Seq ◽

Final Body Weight

Background: Rabbit produce two kinds of feces: hard and soft feces, and they have a preference for consuming the latter. Although this habit of rabbits has been reported for many years, little is known on whether this behavior will impact growth performance and metabolism. The RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of the present study was to investigate the effects of fasting caecotrophy on growth performance and lipid metabolism in rabbits. Results: Our results indicated that, compared with the control group, the final body weight, weight gain, liver weight, specific growth rate and feed conversion ratio were all decreased in the experimental group (P<0.05). Oil red staining of the liver tissue indicated that fasting caecotrophy resulted in decrease of lipid droplet accumulation. RNA sequencing (RNA-seq) analysis revealed a total of 301.2 million raw reads approximately 45.06 Gb of high-quality clean data. The data were mapped to the rabbit genome (http://www.ensembl.org/Oryctolagus_cuniculus). After a five step filtering process, 14964 genes were identified, including 444 differentially expressed genes (P<0.05, foldchange≥1). Especially for remarkable changes of genes related to lipid metabolism, RT-PCR further validated the remarkable decrease of these genes in fasting caecotrophy group, including CYP7A1, PPARG, ABCA1, ABCB1, ABCG1, GPAM, SREBP, etc. KEGG annotation of the differentially expressed genes indicated that the main pathways affected were retinol metabolism, pentose and glucuronide interactions, starch and sucrose metabolism, fatty acid degradation, steroid hormone biosynthesis. Conclusion: In conclusion, the present study revealed that banning caecotrophy reduced growth rate and altered lipid metabolism, our results laid instructive basis for rabbit feeding and production. These data also provides a reference for studying the effects of soft feces on other small herbivores.

Download Full-text

Exon-level estimates improve the detection of differentially expressed genes in RNA-seq studies

RNA Biology ◽

10.1080/15476286.2020.1868151 ◽

2021 ◽

pp. 1-8

Author(s):

Arfa Mehmood ◽

Asta Laiho ◽

Laura L. Elo

Keyword(s):

Differentially Expressed Genes ◽

Differentially Expressed ◽

Rna Seq ◽

Exon Level

Download Full-text

PSVII-2 Differentially expressed genes and their biological function in skeletal muscle of calves born from cows with or without protein supplementation during mid-gestation

Journal of Animal Science ◽

10.1093/jas/skaa054.293 ◽

2020 ◽

Vol 98 (Supplement_3) ◽

pp. 165-166

Author(s):

Elisa B Carvalho ◽

Letícia P Sanglard ◽

Karolina B Nascimento ◽

Javier M Meneses ◽

Daniel R Casagrande ◽

...

Keyword(s):

Skeletal Muscle ◽

Differentially Expressed Genes ◽

Muscle Development ◽

Negative Binomial ◽

Muscle Protein ◽

Basal Diet ◽

Protein Supplementation ◽

Differentially Expressed ◽

Control Group ◽

Q Value

Abstract Gestating cows have an increased nutrient demand to meet the needs of developing the fetus and the mid-gestation is a critical period for the fetal skeletal muscle development. The aim of this study was to evaluate the skeletal muscle transcriptome in the progeny as a function of the maternal protein nutrition during mid-gestation. Eleven Tabapuã cows and their male calves were used in this study. In the first third of gestation (0 to 100 days of gestation; dg), all cows were kept on pasture. From 100 to 200 dg, the control group (CTRL; 7 animals) received a basal diet achieving 5.5% crude protein (CP), whereas the supplemented group (SUPPL; 4 animals) received a basal diet plus protein supplementation (40% CP). After 200 dg, all animals received the same diet. Weaning was performed at 205 ± 7.5 days of age and animals were kept on pasture until reaching 240 days of age, when they were transferred to a feedlot. Muscle samples were collected at 260 days of age and RNA was extracted for RNA-seq analysis. Gene expression data was analyzed with a negative binomial model to identify (q-value ≤ 0.05) differentially expressed genes (DEG) between treatments. A total of 716 DEG were identified (289 DEG up-regulated and 427 down-regulated in SUPPL group; q-value ≤ 0.05). From the 10 most significant down-regulated DEG in the SUPPL group, two genes associated with apoptotic process were identified: MAPK8IP1 and GRINA, with log2 Fold-Changes (log2FC) of 1.04 and 0.49, respectively. From the 10 most significant up-regulated DEG in the SUPPL group, mTOR was identified, with log2FC=0.31. This is a well-known gene involved in muscle protein synthesis. In conclusion, maternal protein supplementation during mid-gestation affects the expression of genes related to energy metabolism and muscle development, which can lead to long-term impacts on production efficiency.

Download Full-text

HIV-1 Infection Transcriptomics: Meta-Analysis of CD4+ T Cells Gene Expression Profiles

Viruses ◽

10.3390/v13020244 ◽

2021 ◽

Vol 13 (2) ◽

pp. 244 ◽

Cited By ~ 1

Author(s):

Antonio Victor Campos Coelho ◽

Rossella Gratton ◽

João Paulo Britto de Melo ◽

José Leandro Andrade-Santos ◽

Rafael Lima Guimarães ◽

...

Keyword(s):

Gene Expression ◽

T Cells ◽

Differentially Expressed Genes ◽

Cd4 T Cells ◽

Expression Profiles ◽

Meta Analysis ◽

Differentially Expressed ◽

Rna Seq ◽

Infected Cells ◽

Hiv 1

HIV-1 infection elicits a complex dynamic of the expression various host genes. High throughput sequencing added an expressive amount of information regarding HIV-1 infections and pathogenesis. RNA sequencing (RNA-Seq) is currently the tool of choice to investigate gene expression in a several range of experimental setting. This study aims at performing a meta-analysis of RNA-Seq expression profiles in samples of HIV-1 infected CD4+ T cells compared to uninfected cells to assess consistently differentially expressed genes in the context of HIV-1 infection. We selected two studies (22 samples: 15 experimentally infected and 7 mock-infected). We found 208 differentially expressed genes in infected cells when compared to uninfected/mock-infected cells. This result had moderate overlap when compared to previous studies of HIV-1 infection transcriptomics, but we identified 64 genes already known to interact with HIV-1 according to the HIV-1 Human Interaction Database. A gene ontology (GO) analysis revealed enrichment of several pathways involved in immune response, cell adhesion, cell migration, inflammation, apoptosis, Wnt, Notch and ERK/MAPK signaling.

Download Full-text

Phytophthora infestans Sporangia Produced in Culture and on Tomato Leaflet Lesions Show Marked Differences in Indirect Germination Rates, Aggressiveness, and Global Transcription Profiles

Molecular Plant-Microbe Interactions ◽

10.1094/mpmi-09-18-0255-ta ◽

2019 ◽

Vol 32 (5) ◽

pp. 515-526 ◽

Cited By ~ 5

Author(s):

William E. Fry ◽

Sean P. Patev ◽

Kevin L. Myers ◽

Kan Bao ◽

Zhangjun Fei

Keyword(s):

Phytophthora Infestans ◽

Differentially Expressed Genes ◽

Differentially Expressed ◽

Rna Seq ◽

Moist Chamber ◽

Field Isolates ◽

Rxlr Effectors ◽

Transcription Profiles ◽

Independent Field ◽

Pure Cultures

Sporangia of Phytophthora infestans from pure cultures on agar plates are typically used in lab studies, whereas sporangia from leaflet lesions drive natural infections and epidemics. Multiple assays were performed to determine if sporangia from these two sources are equivalent. Sporangia from plate cultures showed much lower rates of indirect germination and produced much less disease in field and moist-chamber tests. This difference in aggressiveness was observed whether the sporangia had been previously incubated at 4°C (to induce indirect germination) or at 21°C (to prevent indirect germination). Furthermore, lesions caused by sporangia from plates produced much less sporulation. RNA-Seq analysis revealed that thousands of the >17,000 P. infestans genes with a RPKM (reads per kilobase of exon model per million mapped reads) >1 were differentially expressed in sporangia obtained from plate cultures of two independent field isolates compared with sporangia of those isolates from leaflet lesions. Among the significant differentially expressed genes (DEGs), putative RxLR effectors were overrepresented, with almost half of the 355 effectors with RPKM >1 being up- or downregulated. DEGs of both isolates include nine flagellar-associated genes, and all were down-regulated in plate sporangia. Ten elicitin genes were also detected as DEGs in both isolates, and nine (including INF1) were up-regulated in plate sporangia. These results corroborate previous observations that sporangia produced from plates and leaflets sometimes yield different experimental results and suggest hypotheses for potential mechanisms. We caution that use of plate sporangia in assays may not always produce results reflective of natural infections and epidemics.

Download Full-text

Bioinformatics-based Screening of Key Genes for Saponin Metabolism in Quinoa

10.21203/rs.3.rs-139481/v1 ◽

2021 ◽

Author(s):

Chengang Guo ◽

Zhimin wei ◽

Wei Lyu ◽

Yanlou Geng

Keyword(s):

Differentially Expressed Genes ◽

Principal Component ◽

Enrichment Analysis ◽

Hierarchical Cluster ◽

Gene Set Enrichment Analysis ◽

Differentially Expressed ◽

Differential Analysis ◽

Rna Seq ◽

Protein Coding ◽

Key Genes

Abstract Quinoa saponins have complex, diverse and evident physiologic activities. However, the key regulatory genes for quinoa saponin metabolism are not yet well studied. The purpose of this study was to explore genes closely related to quinoa saponin metabolism. In this study, the significantly differentially expressed genes in yellow quinoa were firstly screened based on RNA-seq technology. Then, the key genes for saponin metabolism were selected by gene set enrichment analysis (GSEA) and principal component analysis (PCA) statistical methods. Finally, the specificity of the key genes was verified by hierarchical clustering. The results of differential analysis showed that 1654 differentially expressed genes were achieved after pseudogenes deletion. Therein, there were 142 long non-coding genes and 1512 protein-coding genes. Based on GSEA analysis, 116 key candidate genes were found to be significantly correlated with quinoa saponin metabolism. Through PCA dimension reduction analysis, 57 key genes were finally obtained. Hierarchical cluster analysis further demonstrated that these key genes can clearly separate the four groups of samples. The present results could provide references for the breeding of sweet quinoa and would be helpful for the rational utilization of quinoa saponins.

Download Full-text