Gene Signature of Regulatory T Cells Isolated from Children with Selective IgA Deficiency and Common Variable Immunodeficiency

Selective IgA deficiency (SIgAD) is the most common form and common variable immunodeficiency (CVID) is the most symptomatic form of predominant antibody deficiency. Despite differences in the clinical picture, a similar genetic background is suggested. A common feature of both disorders is the occurrence of autoimmune conditions. Regulatory T cells (Tregs) are the major immune cell type that maintains autoimmune tolerance. As the different types of abnormalities of Treg cells have been associated with autoimmune disorders in primary immunodeficiency (PID) patients, in our study we aimed to analyze the gene expression profiles of Treg cells in CVID and SIgAD patients compared to age-matched healthy controls. The transcriptome-wide gene profiling was performed by microarray technology. As a result, we analyzed and visualized gene expression patterns of isolated population of Treg cells. We showed the differences at the gene level between patients with and without autoimmunizations. Our findings suggest that the gene signatures of Treg cells isolated from SIgAD and CVID patients differ from age-matched healthy controls and from each other, presenting transcriptional profiles enriched in innate immune or Th response, respectively. The occurrence of autoimmunity in both types of PID is associated with down-regulation of class I IFNs signaling pathways. In summary, our findings improve our understanding of Treg dysfunctions in patients with common PIDs and associated autoimmunity.

Uncovering lupus nephritis-specific genes and the potential of TNFRSF17-targeted immunotherapy: a high-throughput sequencing study

Introduction

Lupus nephritis (LN) is a severe manifestation of systemic lupus erythematosus (SLE). This study aimed to identify LN specific-genes and potential therapeutic targets.

Methods

We performed high-throughput transcriptome sequencing on peripheral blood mononuclear cells (PBMCs) from LN patients. Healthy individuals and SLE patients without LN were used as controls. To validate the sequencing results, qRT-PCR was performed for 5 upregulated and 5 downregulated genes. Furthermore, the effect of the TNFRSF17-targeting drug IBI379 on patient plasma cells and B cells was evaluated by flow cytometry.

Results

Our analysis identified 1493 and 205 differential genes in the LN group compared to the control and SLE without LN groups respectively, with 70 genes common to both sets, marking them as LN-specific. These LN-specific genes were significantly enriched in the 'regulation of biological quality' GO term and the cell cycle pathway. Notably, several genes including TNFRSF17 were significantly overexpressed in the kidneys of both LN patients and NZB/W mice. TNFRSF17 levels correlated positively with urinary protein levels, and negatively with complement C3 and C4 levels in LN patients. The TNFRSF17-targeting drug IBI379 effectively induced apoptosis in patient plasma cells without significantly affecting B cells.

Discussion

Our findings suggest that TNFRSF17 could serve as a potential therapeutic target for LN. Moreover, IBI379 is presented as a promising treatment option for LN.

Changes in the Dentate Gyrus Gene Expression Profile Induced by Levetiracetam Treatment in Rats with Mesial Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) is one of the most common forms of focal epilepsy. Levetiracetam (LEV) is an antiepileptic drug whose mechanism of action at the genetic level has not been fully described. Therefore, the aim of the present work was to evaluate the relevant gene expression changes in the dentate gyrus (DG) of LEV-treated rats with pilocarpine-induced TLE. Whole-transcriptome microarrays were used to obtain the differential genetic profiles of control (CTRL), epileptic (EPI), and EPI rats treated for one week with LEV (EPI + LEV). Quantitative RT-qPCR was used to evaluate the RNA levels of the genes of interest. According to the results of the EPI vs. CTRL analysis, 685 genes were differentially expressed, 355 of which were underexpressed and 330 of which were overexpressed. According to the analysis of the EPI + LEV vs. EPI groups, 675 genes were differentially expressed, 477 of which were downregulated and 198 of which were upregulated. A total of 94 genes whose expression was altered by epilepsy and modified by LEV were identified. The RT-qPCR confirmed that LEV treatment reversed the increased expression of Hgf mRNA and decreased the expression of the Efcab1, Adam8, Slc24a1, and Serpinb1a genes in the DG. These results indicate that LEV could be involved in nonclassical mechanisms involved in Ca2+ homeostasis and the regulation of the mTOR pathway through Efcab1, Hgf, SLC24a1, Adam8, and Serpinb1a, contributing to reduced hyperexcitability in TLE patients.

Comparative Analysis of Transcriptome Profiles in Patients with Thromboangiitis Obliterans

BACKGROUND

Thromboangiitis obliterans (TAO) causes vascular insufficiency due to chronic inflammation and abrupt thrombosis of the medium and small arteries of the extremities. In our study, we aimed to determine biomarkers for the diagnosis of TAO by evaluating 15 male TAO patients with Shinoya diagnostic criteria and 5 healthy controls who did not have TAO-related symptoms in their family histories.

METHODS

The Clariom D Affymetrix platform was used to conduct microarray analysis on total RNA extracted from whole blood. A total of 477 genes (FC ≤ 5 or >5) common to the fifteen patient and five control samples were selected using comparative microarray analysis; among them, 79 genes were upregulated and 398 genes were downregulated.

RESULTS

According to FC ≤ 10 or >10, in the same TAO patient and control group, 13 genes out of 28 were upregulated, whereas 15 genes were downregulated. The 11 key genes identified according to their mean log2FC values were PLP2, RPL27A, CCL4, FMNL1, EGR1, EIF4A1, RPL9, LAMP2, RNF149, EIF4G2, and DGKZ. The genes were ranked according to their relative expression as follows: FMNL1 > RNF149 > RPL27A > EIF4G2 > EIF4A1 > LAMP2 > EGR1 > PLP2 > DGKZ > RPL9 > CCL4. Using protein-protein interaction network analysis, RPL9, RPL27A, and RPL32 were found to be closely related to EIF4G2 and EIF4A1. The Reactome pathway found pathways linked to 28 genes. These pathways included the immune system, cellular responses to stress, cytokine signaling in the immune system, and signaling by ROBO receptors.

CONCLUSIONS

By figuring out the protein expression levels of the genes that have been found to explain how TAO disease works at the molecular level, it will be possible to figure out how well these chosen transcripts can diagnose and predict the disease.

Cornelia de Lange Syndrome Caused by an Intragenic Heterozygous Deletion in RAD21 Detected through Very-High-Resolution Chromosomal Microarray Analysis

Cornelia de Lange syndrome is a genetic and clinically heterogeneous entity, caused by at least five genes. It is characterized by short stature, gestalt facies, microcephaly, neurodevelopmental disorders, and other anomalies. In this report, we present a 13-year-old female patient with microcephaly, cleft palate, polydactyly, short stature, triangular facies, frontal bossing, a bulbous nose, an overfolded helix, limited pronosupination, and an anomalous uterus. No neurodevelopmental disorders were reported. A chromosomal microarray analysis of 6.5 million markers was performed in the proband and her parents. The results showed a de novo heterozygous microdeletion of exons 9-14 within RAD21, which confirmed the diagnosis of Cornelia de Lange syndrome type 4. Our patient did not show any neurologic phenotype (until the time of diagnosis), although neurodevelopmental disorders are frequently present in patients with Cornelia de Lange syndrome type 4, and despite carrying a deletion that was larger than previously reported. Therefore, unknown genetic modifiers or intrinsic mechanisms of RAD21 variants may exist and should be studied.

Microarray Analysis for Transcriptomic Profiling of Myocardium in Patients with Fatal Myocardial Infarction

Transcriptomic evidence from human myocardium in myocardial infarction (MI) is still not sufficient. Thus, there is a need for studies on human cardiac samples in relation to the clinical data of patients. The purpose of our pilot study was to investigate the transcriptomic profile of myocardium in the infarct zone, in comparison to the remote myocardium, in patients with fatal MI, via microarray analysis. This study included four patients with fatal MI type 1. We selected histologically verified samples from within the infarct area (n = 4) and remote myocardium (n = 4). The whole transcriptome was evaluated using microarray analysis. Differentially expressed genes (DEGs) clustered in the infarct area and in the remote myocardium allowed their differentiation. We identified a total of 1785 DEGs (8.32%) in the infarct area, including 1692 up-regulated (94.79%) and 93 down-regulated (5.21%) genes. The top 10 up-regulated genes were TRAIL, SUCLA2, NAE1, PDCL3, OSBPL5, FCGR2C, SELE, CEP63, ST3GAL3 and C4orf3. In the infarct area, we found up-regulation of seventeen apoptosis-related genes, eleven necroptosis-related, and six necrosis-related genes. Transcriptome profiling of the myocardium in patients with MI remains a relevant area of research for the formation of new scientific hypotheses and a potential way to increase the translational significance of studies into myocardial infarction.

Transcriptome profiling of regulatory T cells from children with transient hypogammaglobulinemia of infancy

Transient hypogammaglobulinemia of infancy (THI) is one of the most common forms of hypogammaglobulinemia in the early childhood. THI is usually associated with chronic, recurrent bacterial and viral infections, life-threatening in some cases, yet its pathogenesis is still largely unknown. As our previous findings indicated the possible role of Treg cells in the pathomechanism of THI, the aim of the current study was to investigate gene expression profile of Treg cells isolated from THI patients. The transcriptome-wide gene profiling was performed using microarray technology on THI patients in two time-points: during (THI-1), and in resolution phase (THI-2) of hypogammaglobulinemia. As a result, a total of 1086 genes were differentially expressed in THI-1 patients, when compared to THI-2 as well as control group. Among them, 931 were up- and 155 downregulated, and part of them encodes genes important for Treg lymphocyte biology and function, i.e. transcription factors/cofactors that regulate FOXP3 expression. Thus, we postulate that Treg cells isolated from THI patients during hypogammaglobulinemia display enhanced suppressor transcriptome signature. Treg expression profile of THI children after normalization of Ig levels largely resembles the results obtained in healthy control group, suggesting THI Treg transcriptome seems to return to that observed in healthy children. Taken together, we suggest that THI pathomechanism is associated not only with transiently elevated Treg cell numbers, but also with their enhanced regulatory/inhibitory functions. These findings expand our knowledge of human Treg cells and may be useful for the future diagnosis or management of THI.

Identification and functional prediction of long non-coding RNA and mRNA related to connective tissue disease-associated interstitial lung diseases

Objective

Recently, the role of long non-coding RNA (lncRNA) in rheumatic immune diseases has attracted widespread attention. However, knowledge of lncRNA in connective tissue disease-associated interstitial lung disease (CTD-ILD) is limited. This study explored the expression profile and possible mechanisms of lncRNA and mRNA in peripheral blood mononuclear cells (PBMCs) of CTD-ILD patients, especially systemic sclerosis (SSc)-ILD and rheumatoid arthritis (RA)-ILD.

Methods

LncRNA microarray analysis identified 240 diferentially expressed lncRNAs and 218 diferentially expressed mRNA in the CTD-ILD group and the connective tissue disease without associated interstitial lung disease (CTD-NILD) group. The bioinformatics analysis of diferential genes has identified several important biological processes and signal pathways, including nuclear factor kappa B (NF-kappa B) signaling pathway, interleukin 17 (IL-17) signaling pathway, B cell receptor signaling pathway. Relative expression levels of five diferentially expressed lncRNAs and one mRNA in 120 SSc and RA patients with or without ILD were detected by quantitative reverse-transcription (PCR).

Results

The ENST00000604692 expression level was significantly higher in the ILD than the without interstitial lung disease (NILD) group; T311354 and arginase-1 were significantly higher in SSc than RA group.

Conclusion

These data suggest that the specific profile of lncRNA in PBMCs of CTD-ILD patients and the potential signal pathways related to the pathogenesis of CTD-ILD, which may provide newfound insights for the diagnosis and treatment of CTD-ILD patients.

Low-Cost High-Throughput Genotyping for Diagnosing Familial Hypercholesterolemia

BACKGROUND

Familial hypercholesterolemia (FH) is a common but underdiagnosed genetic disorder characterized by high low-density lipoprotein cholesterol levels and premature cardiovascular disease. Current sequencing methods to diagnose FH are expensive and time-consuming. In this study, we evaluated the accuracy of a low-cost, high-throughput genotyping array for diagnosing FH.

METHODS

An Illumina Global Screening Array was customized to include probes for 636 variants, previously classified as FH-causing variants. First, its theoretical coverage was assessed in all FH variant carriers diagnosed through next-generation sequencing between 2016 and 2022 in the Netherlands (n=1772). Next, the performance of the array was validated in another sample of FH variant carriers previously identified in the Dutch FH cascade screening program (n=1268).

RESULTS

The theoretical coverage of the array for FH-causing variants was 91.3%. Validation of the array was assessed in a sample of 1268 carriers of whom 1015 carried a variant in LDLR, 250 in APOB, and 3 in PCSK9. The overall sensitivity was 94.7% and increased to 98.2% after excluding participants with variants not included in the array design. Copy number variation analysis yielded a 89.4% sensitivity. In 18 carriers, the array identified a total of 19 additional FH-causing variants. Subsequent DNA analysis confirmed 5 of the additionally identified variants, yielding a false-positive result in 16 subjects (1.3%).

CONCLUSIONS

The FH genotyping array is a promising tool for genetically diagnosing FH at low costs and has the potential to greatly increase accessibility to genetic testing for FH. Continuous customization of the array will further improve its performance.

Data Mining of Microarray Datasets in Translational Neuroscience

Data mining involves the computational analysis of a plethora of publicly available datasets to generate new hypotheses that can be further validated by experiments for the improved understanding of the pathogenesis of neurodegenerative diseases. Although the number of sequencing datasets is on the rise, microarray analysis conducted on diverse biological samples represent a large collection of datasets with multiple web-based programs that enable efficient and convenient data analysis. In this review, we first discuss the selection of biological samples associated with neurological disorders, and the possibility of a combination of datasets, from various types of samples, to conduct an integrated analysis in order to achieve a holistic understanding of the alterations in the examined biological system. We then summarize key approaches and studies that have made use of the data mining of microarray datasets to obtain insights into translational neuroscience applications, including biomarker discovery, therapeutic development, and the elucidation of the pathogenic mechanisms of neurodegenerative diseases. We further discuss the gap to be bridged between microarray and sequencing studies to improve the utilization and combination of different types of datasets, together with experimental validation, for more comprehensive analyses. We conclude by providing future perspectives on integrating multi-omics, to advance precision phenotyping and personalized medicine for neurodegenerative diseases.

Crush Injury and Simulated Flight Effects on Muscle Gene Expression in Female Mice

BACKGROUND

Aeromedical evacuation provides critical care during long-distance transport of injured victims between medical facilities. Often, these victims sustain muscle trauma related to mechanical insults, such as crush. Understanding the effects of flight on injured muscle is important because the aircraft cabin represents an external environment with mild hypoxia-the cabin's altitude is 2,438 m instead of sea level. Because mild hypobaric hypoxia can alter gene expression in normal muscle and affect recovery patterns, it is beneficial to examine whether this type of hypoxia may also alter injury-related genes.

OBJECTIVE

The objective of this study was to verify the hypothesis that differential gene expression occurs in response to mild hypobaric hypoxia exposure in crush-injured muscle during two early recovery (preregeneration stage) time points.

METHODS

Twenty-four female mice were anesthetized, and the right gastrocnemius muscle underwent crush injury. Approximately 24 hours later, mice were exposed to normobaric normoxia or hypobaric hypoxia for 8-9 hours. After 32 or 48 hours of recovery, the mice were euthanized, and the right and left lateral gastrocnemius muscles were collected for microarray and bioinformatics analyses.

RESULTS

The study hypothesis was verified. There were 353 highly upregulated, differentially expressed genes identified in the injured muscle compared to the uninjured muscle. Mid1 was upregulated in both pressure conditions regardless of injury status. There were 52 and 15 differentially expressed genes at 32 and 48 hours postinjury, respectively, in the hypobaric hypoxia-exposed, injured muscle compared to the normobaric normoxia-exposed, injured muscle. The macrophage gene Cd68 correlated with other leukocyte-related genes.

DISCUSSION

These findings expand our understanding of the genetic changes that occur in muscle in response to a crush injury, including those related to the macrophage protein CD68. Nursing interventions addressing adequate functioning after crush muscle injury may need to consider the effects on Cd68 and its closely related genes. In addition, our results suggest a responsiveness of the gene Mid1 to flight-relevant hypobaric hypoxia. Changes in the expression of Mid1 may be appropriate in assessing the long-term health of flight crew members.

A meta-analysis of microarray datasets to identify biological regulatory networks in Alzheimer's disease

Background: Alzheimer's Disease (AD) is an age-related progressive neurodegenerative disorder characterized by mental deterioration, memory deficit, and multiple cognitive abnormalities, with an overall prevalence of ∼2% among industrialized countries. Although a proper diagnosis is not yet available, identification of miRNAs and mRNAs could offer valuable insights into the molecular pathways underlying AD's prognosis. Method: This study aims to utilize microarray bioinformatic analysis to identify potential biomarkers of AD, by analyzing six microarray datasets (GSE4757, GSE5281, GSE16759, GSE28146, GSE12685, and GSE1297) of AD patients, and control groups. Furthermore, this study conducted gene ontology, pathways analysis, and protein-protein interaction network to reveal major pathways linked to probable biological events. The datasets were meta-analyzed using bioinformatics tools, to identify significant differentially expressed genes (DEGs) and hub genes and their targeted miRNAs'. Results: According to the findings, CXCR4, TGFB1, ITGB1, MYH11, and SELE genes were identified as hub genes in this study. The analysis of DEGs using GO (gene ontology) revealed that these genes were significantly enriched in actin cytoskeleton regulation, ECM-receptor interaction, and hypertrophic cardiomyopathy. Eventually, hsa-mir-122-5p, hsa-mir-106a-5p, hsa-mir-27a-3p, hsa-mir16-5p, hsa-mir-145-5p, hsa-mir-12-5p, hsa-mir-128-3p, hsa-mir 3200-3p, hsa-mir-103a-3p, and hsa-mir-9-3p exhibited significant interactions with most of the hub genes. Conclusion: Overall, these genes can be considered as pivotal biomarkers for diagnosing the pathogenesis and molecular functions of AD.

LRP5, SLC6A3, and SOX10 Expression in Conventional Ameloblastoma

Cell proliferation and invasion are characteristic of many tumors, including ameloblastoma, and are important features to target in possible future therapeutic applications.

OBJECTIVE

The objective of this study was the identification of key genes and inhibitory drugs related to the cell proliferation and invasion of ameloblastoma using bioinformatic analysis.

METHODS

The H10KA_07_38 gene profile database was analyzed by Rstudio and ShinyGO Gene Ontology enrichment. String, Cytoscape-MCODE, and Kaplan-Meier plots were generated, which were subsequently validated by RT-qPCR relative expression and immunoexpression analyses. To propose specific inhibitory drugs, a bioinformatic search using Drug Gene Budger and DrugBank was performed.

RESULTS

A total of 204 significantly upregulated genes were identified. Gene ontology enrichment analysis identified four pathways related to cell proliferation and cell invasion. A total of 37 genes were involved in these pathways, and 11 genes showed an MCODE score of ≥0.4; however, only SLC6A3, SOX10, and LRP5 were negatively associated with overall survival (HR = 1.49 (p = 0.0072), HR = 1.55 (p = 0.0018), and HR = 1.38 (p = 0.025), respectively). The RT-qPCR results confirmed the significant differences in expression, with overexpression of >2 for SLC6A3 and SOX10. The immunoexpression analysis indicated positive LRP5 and SLC6A3 expression. The inhibitory drugs bioinformatically obtained for the above three genes were parthenolide and vorinostat.

CONCLUSIONS

We identify LRP5, SLC6A3, and SOX10 as potentially important genes related to cell proliferation and invasion in the pathogenesis of ameloblastomas, along with both parthenolide and vorinostat as inhibitory drugs that could be further investigated for the development of novel therapeutic approaches against ameloblastoma.

Integrated multiple-microarray analysis and mendelian randomization to identify novel targets involved in diabetic nephropathy

Background

Diabetic nephropathy (DN), which is the main cause of renal failure in end-stage renal disease, is becoming a common chronic renal disease worldwide. Mendelian randomization (MR) is a genetic tool that is widely used to minimize confounding and reverse causation when identifying the causal effects of complex traits. In this study, we conducted an integrated multiple microarray analysis and large-scale plasma proteome MR analysis to identify candidate biomarkers and evaluate the causal effects of prospective therapeutic targets in DN.

Methods

Five DN gene expression datasets were selected from the Gene Expression Omnibus. The robust rank aggregation (RRA) method was used to integrate differentially expressed genes (DEGs) of glomerular samples between patients with DN and controls, followed by functional enrichment analysis. Protein quantitative trait loci were incorporated from seven different proteomic genome-wide association studies, and genetic association data on DN were obtained from FinnGen (3676 cases and 283,456 controls) for two-sample MR analysis. External validation and clinical correlation were also conducted.

Results

A total of 82 DEGs (53 upregulated and 29 downregulated) were identified through RRA integrated analysis. The enriched Gene Ontology annotations and Kyoto Encyclopedia of Genes and Genomes pathways of the DEGs were significantly enriched in neutrophil degranulation, neutrophil activation, proteoglycan binding, collagen binding, secretory granule lumen, gluconeogenesis, tricarboxylic acid cycle, and pentose phosphate pathways. MR analysis revealed that the genetically predicted levels of MHC class I polypeptide-related sequence B (MICB), granzyme A (GZMA), cathepsin S (CTSS), chloride intracellular channel protein 5, and ficolin-1 (FCN1) were causally associated with DN risk. Expression validation and clinical correlation analysis showed that MICB, GZMA, FCN1, and insulin-like growth factor 1 may participate in the development of DN, and carbonic anhydrase 2 and lipoprotein lipase may play protective roles in patients with DN.

Conclusion

Our integrated analysis identified novel biomarkers, including MICB and GZMA, which may help further understand the complicated mechanisms of DN and identify new target pathways for intervention.

Early inhalant allergen sensitization at component level: an analysis in atopic Dutch children

Background

Allergic rhinitis is a common respiratory disease in children and sensitization to inhalant allergens plays a significant role in its development. However, limited knowledge exists regarding sensitization profiles of inhalant allergen components in atopic children, particularly in the very young individuals. Understanding these profiles could provide insights into the early development of allergic rhinitis. The objective of this cross-sectional retrospective study was to evaluate the IgE-sensitization profiles to multiple inhalant allergen components and their clinical relevance in Dutch atopic children, with specific focus on children under the age of 4 years.

Methods

A total of 243 atopic children were included in the study and sensitization profiles were analyzed using multiplex microarray analysis (ISAC). Clinical information was obtained from records of a pediatric allergy outpatient clinic between 2011 and 2020. Specific IgE responses to inhalation allergen components from five allergen sources (grass pollen, tree pollen, house dust mite, cat and dog), were examined. The study encompassed children of different age groups and compared those with and without symptoms.

Results

The results demonstrated that sensitization to inhalant allergen components was present in 92% of the cohort. Sensitization was already evident at a young age (87%), including infancy, with a rapid increase in prevalence after 1 year of age. House dust mite emerged as the most predominant sensitizing allergen in early childhood, followed by tree pollen in later years. Sensitization patterns were similar between symptomatic and asymptomatic children, although symptomatic children exhibited higher frequencies and values. The sensitization profiles in very young children were comparable to those of children across all age groups.

Conclusion

These findings highlight the presence of sensitization to inhalant allergen components and the early onset of allergic rhinitis before the age of 4, including infancy, in Dutch atopic children. Notable allergen molecules in Dutch atopic children under the age of 4 years include Bet v 1, Fel d 1, Der f 1, Der p 1, Der p 10 and Phl p 4, with house dust mite sensitization being the most common among Dutch infants. Moreover, the prevalence of sensitization to inhalant allergens in this Dutch cohort surpassed that of general European populations, emphasizing the importance of early assessment and management of allergic rhinitis in young atopic children.

Integrated High-Throughput Bioinformatics (Microarray, RNA-Seq, and RNA Interaction) and qRT-PCR Investigation of BMPR1B Axis as a Potential Diagnostic Biomarker of Isfahan Breast Cancer

Background

According to the bioinformatics analyses and previous studies, bone morphogenetic protein receptor type 1B (BMPR1B) dysregulation could remarkably affect breast cancer (BC) status as a potential biomarker and tumor suppressor. Therefore, the analysis of the expression level of BMPR1B and other relevant biological factors such as microRNAs, long non-coding RNAs, downstream proteins in the relevant signaling pathways, and finding the accurate biological mechanism of BMPR1B could be helpful for a better understanding of BC pathogenicity and discovering the new treatment methods and drugs.

Materials and Methods

R Studio software (4.0.2) was used for microarray data analyses. GSE31448 dataset was downloaded by GEOquery package and analyzed by limma package. STRING and miRWalk online databases and Cytoscape software were used for interaction analyses. Quantitative measurement of BMPR1B expression level was performed by qRT-PCR experiment.

Result

Microarray and real-time PCR analysis revealed that BMPR1B has a significant downregulation in the transforming growth factor (TGF)-beta and bone morphogenic protein (BMP) signaling pathways in BC samples. BMPR1B is a potential diagnostic biomarker, regulated by hsa-miR-181a-5p. Also, BMPR1B regulates the function of BMP2, BMP6, SMAD4, SMAD5, and SMAD6 proteins.

Discussion

BMPR1B have a significant role in the development of BC by regulating the potential proteins' function, playing the diagnostic biomarker role, and regulation of TGF-beta and BMP signaling pathways. The high amount of BMPR1B protein helps in increasing the survival rate of the patients.

Correlation of DNA methylation and lymph node metastasis in papillary thyroid carcinoma

BACKGROUND

Papillary thyroid carcinoma (PTC) is the most common type of thyroid cancer with a primarily good prognosis, and its 10-year survival rate is over 90%. However, PTC is prone to early lymph node metastasis.

METHODS

Thyroid cancer tissues from PTC patients with lymphatic metastasis and normal tissues were collected for DNA methylation analysis. Different methylation sites, different methylation regions, gene-enriched pathways, and protein-protein interactions (PPIs) were analyzed.

RESULTS

There were 1004 differentially methylated sites in the PTC group versus the control group; these involved 479 hypermethylated sites in 415 related genes, 525 hypomethylated sites in 482 related genes, 64 differentially methylated regions located in the CpG island region, 34 differentially methylated genes closely related to thyroid cancer, and 17 genes with differentially methylated genes in the DNA promoter region.

CONCLUSION

NDRG4 hypermethylation and FOXO3, ZEB2, and CDK6 hypomethylation were associated with PTC lymph node metastasis.

Differential expression profile of microRNAs in the lung tissues of coal workers with pneumoconiosis and patients with silicosis

The purpose of this study was to characterize the microRNA (miRNA) profile of the lung tissues from coal workers' pneumoconiosis (CWP) and silicosis and to analyze the changes in downstream genes, biological processes, and signaling pathways based on the differently expressed miRNAs. Lung tissues from three CWP patients, eight silicosis patients, and four healthy controls were collected and analyzed for their miRNA profiles using Affymetrix® GeneChip® miRNA Arrays. Differentially expressed miRNAs (DEMs) were identified between the different groups. The miRanda and TargetScan databases were used to predict the putative target genes, and volcano and heat maps were drawn. Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment analyses were then performed to screen the DEMs-associated biological process and signaling pathways, respectively. Further identification with a comprehensive literature research involving particle exposure, fibrosis, inflammation and lung cancer were used to further screen DEMs of CWP and silicosis. Microarray data showed that 375 and 88 miRNAs were differentially expressed in CWP and silicosis lung tissues compared with healthy lung tissues, while 34 miRNAs were differentially expressed in CWP compared with silicosis lung tissues. The GO and KEGG pathway analyses showed that, the target genes were mainly enriched in the TGF-β, MAPK, p53 and other signal pathways. These results provided insight into the miRNA-related underlying mechanisms of CWP and silicosis, and they provided new clues for miRNAs as biomarkers for the diagnosis and differential diagnosis of these two diseases.

Classification and survival prediction in early-stage cirrhosis by gene expression profiling

Liver cirrhosis has been increasingly diagnosed at an early stage owing to the non-invasive diagnostic techniques. However, it is difficult to identify patients at high risk of disease progression. Screening cirrhotic patients with poor prognosis who are most in need of surveillance is still challenging. Gene expression data GSE15654 and GSE14520 were downloaded for performing unsupervised clustering analysis. The prognostic differences between the different clusters were explored by Cox regression. Integrative analysis of gene expression signature, immune cell enrichments and clinical characterization was performed for different clusters. Two distinctive subclasses were identified in HCV-related GSE15654, and Kaplan-Meier analysis indicated that subtype 2 had lower survival rates than subtype 1 (p = 0.0399). Further analysis revealed subtype 2 had a higher density of follicular T helper cells, resting natural killer cells and M0, M2 macrophages while subtype 1 with a higher fraction of naive B cells, memory B cells, resting memory CD 4 T cells, activated natural killer cells and monocytes. 226 differentially expressed genes were identified between the two subtypes, and Reactome analysis showed the mainly enriched pathways were biological oxidations and fatty acid metabolism. Five hub genes (AKT1, RPS16, CDC42, CCND1 and PCBP2) and three significant modules were extracted from the PPI network. The results were validated in HBV-related GSE14520 cohort. We identified two subtypes of patients with different prognosis for hepatitis C-related early-stage liver cirrhosis. Bioinformatics analysis of the gene expression and immune cell profile may provide fresh insight into understanding the prognosis difference.

SNPs Sets in Codifying Genes for Xenobiotics-Processing Enzymes Are Associated with COPD Secondary to Biomass-Burning Smoke

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide; the main risk factors associated with the suffering are tobacco smoking (TS) and chronic exposure to biomass-burning smoke (BBS). Different biological pathways have been associated with COPD, especially xenobiotic or drug metabolism enzymes. This research aims to identify single nucleotide polymorphisms (SNPs) profiles associated with COPD from two expositional sources: tobacco smoking and BBS. One thousand-five hundred Mexican mestizo subjects were included in the study and divided into those exposed to biomass-burning smoke and smokers. Genome-wide exome genotyping was carried out using Infinium Exome-24 kit arrays v. 1.2. Data quality control was conducted using PLINK 1.07. For clinical and demographic data analysis, Rstudio was used. Eight SNPs were found associated with COPD secondary to TS and seven SNPs were conserved when data were analyzed by genotype. When haplotype analyses were carried out, five blocks were predicted. In COPD secondary to BBS, 24 SNPs in MGST3 and CYP family genes were associated. Seven blocks of haplotypes were associated with COPD-BBS. SNPs in the ARNT2 and CYP46A1 genes are associated with COPD secondary to TS, while in the BBS comparison, SNPs in CYP2C8, CYP2C9, MGST3, and MGST1 genes were associated with increased COPD risk.

Disruption of the anterior commissure in Olig2 deficient mice

In the present study, we examined neural circuit formation in the forebrain of the Olig2 knockout (Olig2-KO) mouse model and found disruption of the anterior commissure at the late foetal stage. Axon bundles of the anterior commissure encountered the wall of the third ventricle and ceased axonal extension. L1-CAM immunohistochemistry showed that Olig2-KO mice lose decussation formation in the basal forebrain. DiI tracing revealed that the thin bundles of the anterior commissure axons crossed the midline but ceased further extension into the deep part of the contralateral side. Furthermore, some fractions of DiI-labelled axons were oriented dorsolaterally, which was not observed in the control mouse forebrain. The rostral part of the third ventricle was much wider in the Olig2-KO mice than in wild-type mice, which likely resulted in the delay of midline fusion and subsequent delay and malformation of the anterior commissure. We analysed gene expression alterations in the Olig2-KO mice using a public database and found multiple genes, which are related to axon guidance and epithelial-mesenchymal transition, showing subtle expression changes. These results suggest that Olig2 is essential for anterior commissure formation, likely by regulating multiple biological processes.

Determination of BRCAness Phenotype in Breast Tumors for the Appointment of Neoadjuvant Chemotherapy Based on Platinum and Taxanes

The concept of BRCAness was developed because of similarities between sporadic and hereditary breast cancer. BRCAness defines the pathogenesis and treatment sensitivity of many types of cancer, as well as the presence of a defect in the homologous recombination repair of tumor cells simulating the loss of BRCA1 or BRCA2, as in the presence of germline mutations. The question of treatment effectiveness for BRCA-like tumors is controversial and open. Thus, the aim of this work was to study the effectiveness of neoadjuvant chemotherapy (NAC) in BRCA-deficient breast cancer patients without germline mutations. The study involved 130 patients with breast cancer in stages IIA-IIIB. The treatment regimen included neoadjuvant chemotherapy, surgery, and adjuvant chemotherapy. The materials used were tumor samples from before and after chemotherapy. DNA and RNA were isolated from the tumor material. RNA was used to assess the expression level of BRCA1, while DNA was used for methyl-sensitive PCR. A microarray analysis was performed on high-density DNA chips from an Affymetrix CytoScan^TM HD Array to assess DNA copy number aberration (CNA status) and loss of heterozygosity. A statistical analysis was performed using the Statistica 8.0 application package. It was noted that the existence of copy number aberrations in genes was statistically significantly associated with tumor treatment response and disease prognosis. Patients with partial regression had a statistically significantly higher amount of deletion than patients without an objective response (5/25 patients; 16%), as shown in the general sample of patients (52.9% versus 27.1%, respectively) at p = 0.0001 and in patients treated with anthracycline-containing regimen (p = 0.0001). In addition, it was shown that patients with BRCA1 deletion had higher rates of metastatic-free survival (log rank test, p = 0.009). BRCAness patients had a higher rate of 5-year metastatic survival, but not of treatment efficacy. The prospective study showed the positive effect of assessing the BRCAness phenotype of a tumor before treatment and of prescribing personalized NAC regimens. The objective response rate was statistically significantly more often observed in the group of patients with personalized chemotherapy (85.0% (34/40 patients) versus 62.3% (56/90 patients); p = 0.007). Despite the controversial effectiveness of BRCA-like tumor treatment, our data showed high predictive and prognostic significance of the BRCAness phenotype for the personalization of platinum and taxane regimens.

Pathway expression analysis

This paper introduces a pathway expression framework as an approach for constructing derived biomarkers. The pathway expression framework incorporates the biological connections of genes leading to a biologically relevant model. Using this framework, we distinguish between shedding subjects post-infection and all subjects pre-infection in human blood transcriptomic samples challenged with various respiratory viruses: H1N1, H3N2, HRV (Human Rhinoviruses), and RSV (Respiratory Syncytial Virus). Additionally, pathway expression data is used for selecting discriminatory pathways from these experiments. The classification results and selected pathways are benchmarked against standard gene expression based classification and pathway ranking methodologies. We find that using the pathway expression data along with selected pathways, which have minimal overlap with high ranking pathways found by traditional methods, improves classification rates across experiments.

Conventional Ameloblastoma. A Case Report with Microarray and Bioinformatic Analysis

Ameloblastoma is a rare benign epithelial odontogenic neoplasm, but with great clinical implications, as despite its benignity and slow growth, most cases are locally aggressive with a significant recurrence rate. Histological, cellular, or molecular analyses of its pathogenesis have confirmed the complexity of this neoplasm. We present the case of a 20-year-old patient with a suggestive clinical and radiographic diagnosis of ameloblastoma. An incisional biopsy was obtained confirming the diagnosis of conventional ameloblastoma. Left hemimandibulectomy and plate reconstruction were performed. Histopathological analysis of the surgical specimen confirmed the conventional ameloblastoma with a plexiform pattern and significant areas of cystic degeneration and amyloid-like-like deposits. Additionally, a microarray was carried out with bioinformatic analysis for the enrichment, protein interaction, and determination of eight hub genes (CRP, BCHE, APP, AKT1, AGT, ACTC1, ADAM10, and APOA2) related to their pathogenesis.

Circular RNAs expression profiles and bioinformatics analysis in bronchopulmonary dysplasia

BACKGROUND

Bronchopulmonary dysplasia (BPD) has long been considered the most challenging chronic lung disease for neonatologists and researchers due to its complex pathological mechanisms and difficulty in prediction. Growing evidence indicates that BPD is associated with the dysregulation of circular RNAs (circRNAs). Therefore, we aimed to explore the expression profiles of circRNAs and investigate the underlying molecular network associated with BPD.

METHODS

Peripheral blood was collected from very-low-birth-weight (VLBW) infants at 5-8 days of life to extract PBMCs. Microarray analysis and qRT-PCR tests were performed to determine the differentially expressed circRNAs (DEcircRNAs) between BPD and non-BPD VLBW infants. Simultaneous analysis of GSE32472 was conducted to obtain differentially expressed mRNAs (DEmRNA) from BPD infants. The miRNAs were predicted by DEcircRNAs and DEmRNAs of upregulated, respectively, and then screened for overlapping ones. GO and KEGG analysis was performed following construction of the competing endogenous RNA regulatory network (ceRNA) for further investigation.

RESULTS

A total of 65 circRNAs (52 upregulated and 13 downregulated) were identified as DEcircRNAs between the two groups (FC >2.0 and p.adj <0.05). As a result, the ceRNA network was constructed based on three upregulated DEcircRNAs validated by qRT-PCR (hsa_circ_0007054, hsa_circ_0057950, and hsa_circ_0120151). Bioinformatics analysis indicated these DEcircRNAs participated in response to stimulus, IL-1 receptor activation, neutrophil activation, and metabolic pathways.

CONCLUSIONS

In VLBW infants with a high risk for developing BPD, the circRNA expression profiles in PBMCs were significantly altered in the early post-birth period, suggesting immune dysregulation caused by infection and inflammatory response already existed.

Morinda officinalis polysaccharide enable suppression of osteoclastic differentiation by exosomes derived from rat mesenchymal stem cells

CONTEXT

Morinda officinalis F.C. How. (MO) (Rubiaceae) can strengthen bone function.

OBJECTIVE

To examine the functional mechanism and effect of MO polysaccharides (MOPs) in rats with glucocorticoid-induced osteoporosis (GIOP).

MATERIALS AND METHODS

Rats with GIOP were treated with 5, 15 or 45 mL/kg of MOP [n = 15 for each dose, intraperitoneal (i.p.) injection every other day for 8 weeks]. The body weight of rats and histomorphology of bone tissues were examined. Bone marrow mesenchymal stem cells (BMSCs)-derived exosomes (Exo) were collected and identified. Bone marrow-derived macrophages (BMMs) were induced to differentiate into osteoclasts and treated with BMSC-Exo for in vitro studies.

RESULTS

MOP reduced the body weight (5, 15, or 45 mg/kg MOP vs. phosphate-buffered saline: 8%, 15% and 25%, p < 0.01), elevated the bone volume to tissue volume (BV/TV), mean trabecular thickness (Tb.Th), mean trabecular number (Tb.N) and mean connectivity density (Conn.D) (40-86%, p < 0.01), decreased the mean trabecular separation/spacing (Tb.Sp) (22-37%, p < 0.01), increased the cortical bone continuity (35-90%, p < 0.01) and elevated RUNX family transcription factor 2 and RANK levels (5-12%, p < 0.01), but suppressed matrix metallopeptidase 9 and cathepsin K levels (9-20%, p < 0.01) in femur tissues. BMSC-Exo from MOP-treated rats (MOP-Exo) suppressed osteoclastic differentiation and proliferation of BMMs. The downregulation of microRNA-101-3p (miR-101-3p) or the upregulation of prostaglandin-endoperoxide synthase 2 (PTGS2) blocked the functions of MOP-Exo.

DISCUSSION AND CONCLUSIONS

MOP inhibits osteoclastic differentiation and could potentially be used for osteoporosis management. This suppression may be enhanced by the upregulation of miR-101-3p or the inhibition of PTGS2.

Microarray Analysis Reveals Changes in tRNA-Derived Small RNAs (tsRNAs) Expression in Mice with Septic Cardiomyopathy

Background: tRNA-derived small RNAs (tsRNAs) as a novel non-coding RNA have been studied in many cardiovascular diseases, but the relationship between tsRNAs and septic cardiomyopathy has not been investigated. We sought to analyze changes of the expression profile of tsRNAs in septic cardiomyopathy and reveal an important role for tsRNAs. Methods: We constructed a sepsis model by cecal ligation and puncture (CLP) in mice, and microarray analysis was used to find differentially expressed tsRNAs. Quantitative real-time PCR was used to verify the expression of tsRNAs and the interference effect of angiogenin (ANG), a key nuclease producing tsRNAs. Bioinformatics analysis was used to predict target genes and functions. CCK-8 and LDH release assays were used to detect cell viability and cell death. Results: A total of 158 tsRNAs were screened, of which 101 were up-regulated and 57 were down-regulated. A total of 8 tsRNAs were verified by qPCR, which was consistent with microarray results. Gene Ontology (GO) annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses suggest that these tsRNAs may be associated with the Wnt signaling pathway and participate in cellular process. The expression of tsRNAs decreased after the interference of the key nuclease ANG, while CCK-8 suggested a corresponding decrease in cell viability and an increase in the release of LDH (cell death), indicating that tsRNAs can protect cardiomyocytes during the development of septic cardiomyopathy, reduced cardiomyocyte death. Conclusions: A total of 158 tsRNAs changed significantly in septic cardiomyopathy, and these tsRNAs may play a protective role in the development of septic cardiomyopathy.

[Corrigendum] Role of 4‑aminobutyrate aminotransferase (ABAT) and the lncRNA co‑expression network in the development of myelodysplastic syndrome

Subsequently to the publication of the above article, and a Corrigendum that has already been published with the intention of showing corrected versions of Figs. 1 and 8 (DOI: 10.3892/or.2022.8348; published online on June 14, 2022), the authors have belatedly realized that the revisions made to Fig. 8 necessitated changes that should have been introduced into Fig. 9, although these were not attended to in the first corrigendum. Essentially, Fig. 8 was revised as the cell apoptosis and cell proliferation assays therein were poorly presented, which made the interpretation of the data difficult; Fig. 9 showed the fractions of apoptotic cells in the SKM‑1 and THP‑1 cell lines with lncENST00000444102 overexpression as this pertained to Fig. 8. A revised version of Fig. 9, presenting the analysis of the data shown in the revised version of Fig. 8, is shown opposite. In addition to the revision of Fig. 9, the sentence starting on p. 517, left‑hand column, line 12 ["The flow cytometric apoptosis assay revealed that lncENST00000444102 overexpression promoted tumor cells to undergo apoptosis compared to control cells (P<0.001, Fig. 9)"] should be replaced with the following text, to reflect the change in the level of statistical significance: 'The flow cytometric apoptosis assay revealed that lncENST00000444102 overexpression promoted tumor cells to undergo apoptosis compared to control cells (P<0.01, Fig. 9)". Note that the revisions made to Figs. 8 and 9 in this paper have not had a major impact on the reported results, and do not affect the overall conclusions reported in the study. All the authors agree to the publication of this corrigendum. The authors are grateful to the Editor of Oncology Reports for allowing them the opportunity to publish this additional Corrigendum; furthermore, they apologize for any inconvenience caused to the readership of the Journal. [Oncology Reports 42: 509‑520, 2019; DOI: 10.3892/or.2019.7175].

Whole-genome identification and construction of the lncRNA-mRNA co-expression network in patients with actinic keratosis

BACKGROUND

Actinic keratosis (AK) is a common premalignant lesion induced by chronic exposure to ultraviolet radiation and may develop into invasive cutaneous squamous carcinoma (cSCC). The identification of specific biomarkers in AK are still unclear. Long non-coding RNAs (lncRNAs), as transcripts of more than 200 nucleotides, significantly involving in multiple biologic processes, especially in the development of tumors.

METHODS

In our study, we obtained data from RNA-sequencing analysis using two AK lesion tissues and three normal cutaneous tissues to comparatively analyze the differentially expressed (DE) lncRNAs and messenger RNAs (mRNAs). Firstly, we used microarray analyses to identify DE lncRNAs and DE mRNAs. Secondly, we performed Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis to analyze the primary function and find out significant pathways of these DE mRNA and lncRNAs. Finally, we used the top ten DE lncRNAs to construct a lncRNA-mRNA co-expression network.

RESULTS

Our results showed that there were a total of 2,097 DE lncRNAs and 2,043 DE mRNAs identified. GO and KEGG analysis and the lncRNA-mRNA co-expression network (using the top 10 DE lncRNAs comprises 130 specific co-expressed mRNAs to construct) indicated that lncRNA uc011fnr.2 may negatively regulate SCIMP and Toll-like receptor 4 (TLR4) and play an important role in Janus kinase-signal transducer and activator of transcription 3 (JAK-STAT3) signaling pathway of AK.

CONCLUSIONS

lncRNA uc011fnr.2 may play an important role in JAK-STAT3 signaling pathway of AK by modulating SCIMP, TLR4 and IL-6. Further research is required to validate the value of lncRNA uc011fnr.2 in the progression of AK.

Serum microRNAs as non-invasive diagnostic biomarkers for intrahepatic cholestasis of pregnancy

OBJECTIVES

Intrahepatic cholestasis of pregnancy (IHCP) causes itching, preterm birth, and stillbirth. However, there is no accurate diagnostic method for IHCP. Currently, circulating microRNAs (miRNAs) have become candidate biomarkers for the diagnosis of multiple diseases. Here, we investigated the diagnostic value of miRNAs in IHCP and aimed to predict the molecular mechanism of IHCP pathogenesis.

METHODS

We analyzed differentially expressed miRNAs in both women with IHCP and normal pregnant women. The selected candidate miRNAs were validated in 46 IHCP cases and 46 normal pregnant subjects, and we constructed receiver operator characteristic curves of miRNAs. Pearson correlations between levels of total bile acid (TBA) and differentially expressed miRNAs were also calculated. In addition, we clustered functionally significant biological pathways using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses.

RESULTS

The expression levels of 13 miRNAs were remarkably upregulated while the other 35 miRNAs were significantly downregulated, in women with IHCP (P≤0.05) when compared with healthy pregnant women. The areas under the curves of miRNA-7706, miRNA-877-3p, and miRNA-128-3p were higher than 0.90, indicating more reliable diagnosis of IHCP. The Pearson analysis showed that the levels of these miRNAs were positively correlated to TBA level. Additionally, the results of bioinformatics analysis revealed that the differentially expressed miRNAs mainly influenced fatty acid biosynthesis, the endoplasmic reticulum ubiquitin ligase complex, and the p53, and mammalian target of rapamycin (mTOR) and AMP-activated protein kinase (AMPK) signaling pathways.

CONCLUSION

The panel of three-miRNAs (miRNA-7706, miRNA-877-3p, and miRNA-128-3p) may be a useful noninvasive diagnostic biomarker of IHCP.

A comprehensive bioinformatic analysis revealed novel MicroRNA biomarkers of Parkinson's disease

Parkinson's disease (PD) is categorized as a neurodegenerative disorder. Different studies have focused on the role of microRNAs (miRNAs) on PD progression. Due to its complexity in initiation and progression, a considerable requirement has arisen to identify novel miRNA biomarkers in a noninvasive manner. In silico analysis has been used to select differentially expressed miRNAs (DE-miRNAs) and key pathways in this disease. In this manner, several data sets of different neurodegenerative diseases have been analyzed to purify the findings of the present study. Totally, 15 DE miRNAs showed significant changes compared to healthy controls and other neurodegenerative diseases. Then, the targets of the miRNAs were predicted through miRTarBase and TargetScan databases. Besides, enrichment analysis was implemented for predicted target genes. Most of the target genes were enriched in the TRAIL signaling pathway, Regulation of nucleobase, nucleoside, nucleotide and nucleic acid metabolism, protein serine/threonine kinase activity, and Cytoplasm. Moreover, a protein-protein interaction network was constructed to find the most key DE miRNAs and targets in this disease. The results of the present study may help researchers shed light on the discovery of novel biomarkers for PD.

Scalable Additive Construction of Arrayed Microstructures with Encoded Properties for Bioimaging

Microarrays are essential components of analytical instruments. The elements of microarrays may be imbued with additional functionalities and encodings using composite materials and structures, but traditional microfabrication methods present substantial barriers to fabrication, design, and scalability. In this work, a tool-free technique was reported to additively batch-construct micromolded, composite, and arrayed microstructures. The method required only a compatible carrier fluid to deposit a material onto a substrate with some topography. Permutations of this basic fabrication approach were leveraged to gain control over the volumes and positions of deposited materials within the microstructures. As a proof of concept, cell micro-carrier arrays were constructed to demonstrate a range of designs, compositions, functionalities, and applications for composite microstructures. This approach is envisioned to enable the fabrication of complex composite biological and synthetic microelements for biosensing, cellular analysis, and biochemical screening.

SNAI1 is upregulated during muscle regeneration and represses FGF21 and ATF3 expression by directly binding their promoters

During skeletal myogenesis, the zinc-finger transcription factors SNAI1 and SNAI2, are expressed in proliferating myoblasts and regulate the transition to terminally differentiated myotubes while repressing pro-differentiation genes. Here, we demonstrate that SNAI1 is upregulated in vivo during the early phase of muscle regeneration induced by bupivacaine injury. Using shRNA-mediated gene silencing in C2C12 myoblasts and whole-transcriptome microarray analysis, we identified a collection of genes belonging to the endoplasmic reticulum (ER) stress pathway whose expression, induced by myogenic differentiation, was upregulated in absence of SNAI1. Among these, key ER stress genes, such as Atf3, Ddit3/Chop, Hspa5/Bip, and Fgf21, a myokine involved in muscle differentiation, were strongly upregulated. Furthermore, by promoter mutant analysis and Chromatin immune precipitation assay, we demonstrated that SNAI1 represses Fgf21 and Atf3 in proliferating myoblasts by directly binding to multiple E boxes in their respective promoter regions. Together, these data describe a new regulatory mechanism of myogenic differentiation involving the direct repressive action of SNAI1 on ER stress and Fgf21 expression, ultimately contributing to maintaining the proliferative and undifferentiated state of myoblasts.

Association of Ligamentum Flavum Hypertrophy with Adolescent Idiopathic Scoliosis Progression-Comparative Microarray Gene Expression Analysis

The role of the ligamentum flavum (LF) in the pathogenesis of adolescent idiopathic scoliosis (AIS) is not well understood. Using magnetic resonance imaging (MRI), we investigated the degrees of LF hypertrophy in 18 patients without scoliosis and on the convex and concave sides of the apex of the curvature in 22 patients with AIS. Next, gene expression was compared among neutral vertebral LF and LF on the convex and concave sides of the apex of the curvature in patients with AIS. Histological and microarray analyses of the LF were compared among neutral vertebrae (control) and the LF on the apex of the curvatures. The mean area of LF in the without scoliosis, apical concave, and convex with scoliosis groups was 10.5, 13.5, and 20.3 mm2, respectively. There were significant differences among the three groups (p < 0.05). Histological analysis showed that the ratio of fibers (Collagen/Elastic) was significantly increased on the convex side compared to the concave side (p < 0.05). Microarray analysis showed that ERC2 and MAFB showed significantly increased gene expression on the convex side compared with those of the concave side and the neutral vertebral LF cells. These genes were significantly associated with increased expression of collagen by LF cells (p < 0.05). LF hypertrophy was identified in scoliosis patients, and the convex side was significantly more hypertrophic than that of the concave side. ERC2 and MAFB genes were associated with LF hypertrophy in patients with AIS. These phenomena are likely to be associated with the progression of scoliosis.

Identification of Small Regions of Overlap from Copy Number Variable Regions in Patients with Hypospadias

Hypospadias is a common form of congenital atypical sex development that is often associated with other congenital comorbidities. Many genes have been associated with the condition, most commonly single sequence variations. Further investigations of recurrent and overlapping copy number variations (CNVs) have resulted in the identification of genes and chromosome regions associated with various conditions, including differences of sex development (DSD). In this retrospective study, we investigated the DECIPHER database, as well as an internal institutional database, to identify small recurrent CNVs among individuals with isolated and syndromic hypospadias. We further investigated these overlapping recurrent CNVs to identify 75 smallest regions of overlap (SROs) on 18 chromosomes. Some of the genes within these SROs may be considered potential candidate genes for the etiology of hypospadias and, occasionally, additional comorbid phenotypes. This study also investigates for the first time additional common phenotypes among individuals with hypospadias and overlapping CNVs. This study provides data that may aid genetic counseling and management of individuals with hypospadias, as well as improve understanding of its underlying genetic etiology and human genital development overall.

Contribution of FBLN5 to Unstable Plaques in Carotid Atherosclerosis via mir128 and mir532-3p Based on Bioinformatics Prediction and Validation

FBLN5, a member of the short fibulins in the fibulin family of extracellular matrix/matricellular proteins, is involved in interactions with components of the basement membrane and extracellular matrix proteins. It plays key roles in endothelial tissues in many vascular diseases. In this study, the relationship between FBLN5 and carotid atherosclerotic plaque stability as well as the regulatory roles of miRNAs were evaluated. Differential gene expression analyses and weighted gene co-expression network analysis (WGCNA) based on the GSE163154 dataset (including 16 samples without intraplaque hemorrhage and 27 samples with intraplaque hemorrhage) in GEO revealed that FBLN5 is related to plaque stability and is the most significantly differentially expressed gene. LASSO regression was used to evaluate genes obtained from the intersection of differentially expressed genes and clinically significant modules identified by WGCNA. A prediction model based on eight genes, including FBLN5, was constructed and showed an accuracy of 0.951 based on an ROC analysis. Low FBLN5 expression in plaque tissues was confirmed by immunohistochemistry and western blotting. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses showed that FBLN5 acted mainly by the maintenance of the cellular matrix and reactive oxygen species production. miRNAs upstream of these eight predictive genes, including FBLN5, were identified and used to construct a network diagram. These results revealed that hsa-mir-128 and hsa-mir-532–3p were upstream regulatory factors of FBLN5, as verified by PCR assays of human plaque tissues demonstrating that both miRNAs were significantly up-regulated. Therefore, FBLN5 may play an important role in carotid atherosclerosis via hsa-mir-128 and hsa-mir-532–3p as well as become an essential target for treatment.

Investigation of hub genes involved in Turner syndrome using biological informatics methods

BACKGROUND

This study aimed to explore candidate genes and their potential interaction mechanism critical to the pathophysiology of Turner syndrome by using the Gene Expression Omnibus database.

METHODS

GSE58435 data set was obtained by querying the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were screened using R and subsequently annotated by Gene Ontology. Functional enrichment analysis was performed based on the Kyoto Encyclopedia of Genes and Genomes database for annotation, visualization, and integrated discovery. A protein-protein interaction network of different genes was constructed based on the STRING database, in which hub genes were explored through Cytoscape software. The expression of the hub genes was verified by analyzing the gene expression in the GSE46687 data set.

RESULTS

A total of 733 differential genes were identified. These differentially expressed genes were significantly enriched in nucleoplasm and nucleus. Their molecular function was concentrated on DNA binding and transcription, coronary artery, and adipose tissue development. According to the annotation of Kyoto Encyclopedia of Genes and Genomes, the identified DEGs were mainly enriched in inflammatory mediator regulation of TRP channels, osteoclast differentiation. A total of 10 hub genes (HIST1H2BA, TRIM71, HIST1H2BB, HIST1H4D, TNF, TP53BP1, CDCA8, EGF, HMG20B, and BCL9) were identified from the constructed protein-protein interaction network. These genes were discovered to be highly expressed in osteoclasts, ovaries, digestive tract, blood, and lymphatic tissues through the online application of human protein atlas.

CONCLUSION

In this study, 733 DEGs and 10 hub genes were identified. They would be new candidate targets in Turner syndrome.

Long non-coding RNA screening and identification of potential biomarkers for type 2 diabetes

BACKGROUND

To investigate new lncRNAs as molecular markers of T2D.

METHODS

We used microarrays to identify differentially expressed lncRNAs and mRNAs from five patients with T2D and paired controls. Through bioinformatics analysis, qRT-PCR validation, ELISA, and receiver operating characteristic (ROC) curve analysis of 100 patients with T2D and 100 controls to evaluate the correlation between lncRNAs and T2D, and whether lncRNAs could be used in the diagnosis of T2D patients.

RESULTS

We identified 68 and 74 differentially expressed lncRNAs and mRNAs, respectively. The top five upregulated lncRNAs are ENST00000381108.3, ENST00000515544.1, ENST00000539543.1, ENST00000508174.1, and ENST00000564527.1, and the top five downregulated lncRNAs are TCONS_00017539, ENST00000430816.1, ENST00000533203.1, ENST00000609522.1, and ENST00000417079.1. The top five upregulated mRNAs are Q59H50, CYP27A1, DNASE1L3, GRIP2, and lnc-TMEM18-12, and the top five downregulated mRNAs are GSTM4, PODN, GLYATL2, ZNF772, and CLTC. Examination of lncRNA-mRNA interaction pairs indicated that the target gene of lncRNA XR_108954.2 is E2F2. Multiple linear regression analysis showed that XR_108954.2 (r = 0.387, p < 0.01) and E2F2 (r = 0.368, p < 0.01) expression levels were positively correlated with glucose metabolism indicators. Moreover, E2F2 was positively correlated with lipid metabolism indicators (r = 0.333, p < 0.05). The area under the ROC curve was 0.704 (95% CI: 0.578-0.830, p = 0.05) for lncRNA XR_108954.2 and 0.653 (95% CI: 0.516-0.790, p = 0.035) for E2F2.

CONCLUSIONS

This transcriptome analysis explored the aberrantly expressed lncRNAs and identified E2F2 and lncRNA XR_108954.2 as potential biomarkers for patients with T2D.

Differential Transcriptome Profiling Unveils Novel Deregulated Gene Signatures Involved in Pathogenesis of Alzheimer's Disease

Alzheimer’s disease (AD) is a neurodegenerative disorder that is characterized by a progressive loss of cognitive functions at a higher level than normal aging. Although the apolipoprotein (APOE) gene is a major risk factor in developing AD, other genes have also been reported to be linked with complex phenotypes. Therefore, this genome-wide expression study explored differentially expressed genes as possible novel biomarkers involved in AD. The mRNA expression dataset, GSE28146, containing 15 sample data composed of 7 AD cases from the hippocampus region with age-matched control (n = 8, >80 years), was analyzed. Using “affy” R-package, mRNA expression was calculated, while pathway enrichment analysis was performed to determine related biological processes. Of 58 differentially expressed genes, 44 downregulated and 14 upregulated genes were found to be significantly (p < 0.001) altered. The pathway enrichment analysis revealed two altered genes, i.e., dynein light chain 1 (DYNLL1) and kalirin (KLRN), associated with AD in the elderly population. The majority of genes were associated with retrograde endocannabinoid as well as vascular endothelial growth factors affecting the complex phenotypes. The DYNLL1 and KLRN genes may be involved with AD and Huntington’s disease (HD) phenotypes and represent a common genetic basis of these diseases. However, the hallmark of AD is dementia, while the classic motor sign of HD includes chorea. Our data warrant further investigation to identify the role of these genes in disease pathogenesis.

Identification of Significant Genes and Pathways Associated with Tenascin-C in Cancer Progression by Bioinformatics Analysis

Background

Tenascin-C (TNC) is a large glycoprotein of the extracellular matrix which associated with poor clinical outcomes in several malignancies. TNC over-expression is repeatedly observed in several cancer tissues and promotes several processes in tumor progression. Until quite recently, more needs to be known about the potential mechanisms of TNC as a key player in cancer progression and metastasis.

Materials and Methods

In the present study, we performed a bioinformatics analysis of breast and colorectal cancer expression microarray data to survey TNC role and function with holistic view. Gene expression profiles were analyzed to identify differentially expressed genes (DEGs) between normal samples and cancer biopsy samples. The protein-protein interaction (PPI) networks of the DEGs with CluePedia plugin of Cytoscape software were constructed. Furthermore, after PPI network construction, gene-regulatory networks analysis was performed to predict long noncoding RNAs and microRNAs associated with TNC and cluster analysis was performed. Using the Clue gene ontology (GO) plugin of Cytoscape software, the GO and pathway enrichment analysis were performed.

Results

PPI and DEGs-miRNA-lncRNA regulatory networks showed TNC is a significant node in a huge network, and one of the main gene with high centrality parameters. Furthermore, from the regulatory level perspective, TNC could be significantly impressed by miR-335-5p. GO analysis results showed that TNC was significantly enriched in cancer-related biological processes.

Conclusions

It is important to identify the TNC underlying molecular mechanisms in cancer progression, which may be clinically useful for tumor-targeting strategies. Bioinformatics analysis provides an insight into the significant roles that TNC plays in cancer progression scenarios.

Identification of hub genes and pathways in bladder cancer using bioinformatics analysis

Bladder cancer (BC) is the most common malignant tumor of urinary tract system. The aim of this study was to investigate the genetic signatures of bladder cancer (BC) and identify its potential molecular mechanisms. The gene expression profiles of GSE3167 (50 samples, including 41BC and 9 non-cancerous urothelial cells) was downloaded from the GEO database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) were performed to identify enriched pathways, and a protein-protein interaction (PPI) network was used to identify hub genes and for module analysis. Moreover, we conducted expression and survival analyses to screen and validate hub genes. In total, 1528 DEGs were identified in bladder cancer (BC), including 1212 up-regulated genes and 316 down-regulated genes. Up-regulated differentially expressed genes (DEGs) were significantly enriched in negative regulation of macromolecule metabolic process, macromolecule catabolic process, proteolysis and regulation of cell death, while the down-regulated differentially expressed genes (DEGs) were mainly involved in cell surface receptor linked signal transduction, ion transport, cell-cell signaling and defense response. The top 10 hub genes with the highest degrees were selected from the PPI network. These genes included HSP90AA1, MYH11, MYL9, CNN1, ACTC1, RAN, ENO1, HNRNPC, ACTG2 and YWHAZ. From sub-networks, we found these genes were involved in the proteasome, pathways in cancer and cell cycle. Hence, the identified DEGs and hub genes may be beneficial to elucidate the mechanisms underlying BC.

Pseudo-pheochromocytoma due to obstructive sleep apnea: a case report

SUMMARY

Obstructive sleep apnea (OSA) is a condition of intermittent nocturnal upper airway obstruction. OSA increases sympathetic drive which may result in clinical and biochemical features suggestive of pheochromocytoma. We present the case of a 65-year-old male with a 2.9-cm left adrenal incidentaloma on CT, hypertension, symptoms of headache, anxiety and diaphoresis, and persistently elevated 24-h urine norepinephrine (initially 818 nmol/day (89-470)) and normetanephrine (initially 11.2 µmol/day (0.6-2.7)). He was started on prazosin and underwent left adrenalectomy. Pathology revealed an adrenal corticoadenoma with no evidence of pheochromocytoma. Over the next 2 years, urine norepinephrine and normetanephrine remained significantly elevated with no MIBG avid disease. Years later, he was diagnosed with severe OSA and treated with continuous positive airway pressure. Urine testing done once OSA was well controlled revealed complete normalization of urine norepinephrine and normetanephrine with substantial symptom improvement. It was concluded that the patient never had a pheochromocytoma but rather an adrenal adenoma with biochemistry and symptoms suggestive of pheochromocytoma due to untreated severe OSA. Pseudo-pheochromocytoma is a rare presentation of OSA and should be considered on the differential of elevated urine catecholamines and metanephrines in the right clinical setting.

LEARNING POINTS

Obstructive sleep apnea (OSA) is a common condition among adults. OSA may rarely present as pseudo-pheochromocytoma with symptoms of pallor, palpitations, perspiration, headache, or anxiety. OSA should be considered on the differential of elevated urine catecholamines and metanephrines, especially in patients with negative metaiodobenzylguanidine (MIBG) scan results.

Characterization of a Caffeic Acid 8-O-Methyltransferase from Citrus and Its Function in Nobiletin Biosynthesis

Nobiletin (3',4',5,6,7,8-hexamethoxyflavone) is a polymethoxylated flavonoid specifically accumulated in citrus fruit with numerous beneficial effects to human health. In this study, a novel O-methyltransferase (CitOMT2) was isolated from three citrus varieties, Ponkan mandarin (Citrus reticulata Blanco), Nou 6 ("King mandarin" × "Mukaku-kishu"), and Satsuma mandarin (Citrus unshiu Marc.), and its functions were characterized in vitro. The gene expression results showed that CitOMT2 was highly expressed in the two nobiletin abundant varieties of Ponkan mandarin and Nou 6. However, the expression level of CitOMT2 was low in the flavedo of Satsuma mandarin, in which only a small amount of nobiletin was accumulated. Functional analysis suggested that CitOMT2 was a caffeic acid 8-O-methyltransferase, and it catalyzed the O-methylation of 7,8-dihydroxyflavone at 8-OH. As the methylation of flavone at 8-OH was required for nobiletin biosynthesis, the results presented in this study suggested that CitOMT2 was a key gene regulating nobiletin accumulation in citrus fruit.

Expression profiles of long noncoding and messenger RNAs in epicardial adipose tissue derived from patients with coronary atherosclerosis

BACKGROUND

Given its close anatomical location to the heart and its endocrine properties, attention on epicardial adipose tissue (EAT) has increased.

OBJECTIVE

This study investigated the expression profiles of long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) in EAT derived from patients with coronary artery disease (CAD).

METHODS

EAT samples from 8 CAD and 8 non-CAD patients were obtained during open-heart surgery. The expression of lncRNAs and mRNAs in each EAT sample was investigated using microarray analysis and further verified using reverse transcription-quantitative polymerase chain reaction.

RESULTS

Overall, 1,093 differentially expressed mRNAs and 2,282 differentially expressed lncRNAs were identified in EAT from CAD vs non-CAD patients. Analysis using Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes showed that these differentially expressed genes were mainly enriched in various inflammatory, immune, and metabolic processes. They were also involved in osteoclast differentiation, B cell receptor and adipocytokine signaling, and insulin resistance pathways. Additionally, lncRNA-mRNA and lncRNA-target pathway networks were built to identify potential core genes (e.g. Lnc-CCDC68-2:1, AC010148.1, NONHSAT104810) involved in atherosclerosis pathogenesis.

CONCLUSION

In summary, lncRNA and mRNA profiles in EAT were markedly different between CAD and non-CAD patients. Our study identifies several potential key genes and pathways that may participate in atherosclerosis development.

Differential Expression of Long Non-Coding RNAs and Their Role in Rodent Neuropathic Pain Models

Neuropathic pain, which is accompanied by an unpleasant sensation, affects the patient’s quality of life severely. Considering the complexity of the neuropathic pain, there are huge unmet medical needs for it while current effective therapeutics remain far from satisfactory. Accordingly, exploration of mechanisms of neuropathic pain could provide new therapeutic insights. While numerous researches have pointed out the contribution of sensory neuron-immune cell interactions, other mechanisms of action, such as long non-coding RNAs (lncRNAs), also could contribute to the neuropathic pain observed in vivo. LncRNAs have more than 200 nucleotides and were originally considered as transcriptional byproducts. However, recent studies have suggested that lncRNAs played a significant role in gene regulation and disease pathogenesis. A substantial number of long non-coding RNAs were expressed differentially in neuropathic pain models. Besides, therapies targeting specific lncRNAs can significantly ameliorate the development of neuropathic pain, which reveals the contribution of lncRNAs in the generation and maintenance of neuropathic pain and provides a new therapeutic strategy. The primary purpose of this review is to introduce recent studies of lncRNAs on different neuropathic pain models.

Expression profiling of lncRNAs and mRNAs in placental site trophoblastic tumor (PSTT) by microarray

As a rare type of gestational trophoblastic disease, placental site trophoblastic tumor (PSTT) is originated from intermediate trophoblast cells. Long noncoding RNAs (lncRNAs) regulate numerous biological process. However, the role of lncRNAs in PSTT remains poorly understood. In the present study, expression levels of lncRNAs and mRNAs in four human PSTT tissues and four normal placental villi were investigated. The results of microarray were validated by the reverse transcription and quantitative real-time polymerase reaction (RT-qPCR) and immunohistochemistry analyses. Furthermore, GO and KEGG pathway analyses were performed to identify the underlying biological processes and signaling pathways of aberrantly expressed lncRNAs and mRNAs. We also conducted the coding-non-coding gene co-expression (CNC) network to explore the interaction of altered lncRNAs and mRNAs. In total, we identified 1247 up-regulated lncRNAs and 1013 down-regulated lncRNAs as well as 828 up-regulated mRNAs and 1393 down-regulated mRNAs in PSTT tissues compared to normal villi (fold change ≥ 2.0, p < 0.05). GO analysis showed that mitochondrion was the most significantly down-regulated GO term, and immune response was the most significantly up-regulated term. A CNC network profile based on six confirmed lncRNAs (NONHSAT114519, NR_103711, NONHSAT003875, NONHSAT136587, NONHSAT134431, NONHSAT102500) as well as 354 mRNAs was composed of 497 edges. GO and KEGG analyses indicated that interacted mRNAs were enriched in the signal-recognition particle (SRP)-dependent cotranslational protein targeting to membrane and Ribosome pathway. It contributes to expand the understanding of the aberrant lncRNAs and mRNAs profiles of PSTT, which may be helpful for the exploration of new diagnosis and treatment of PSTT.

Analysis of Gene Expression Profiles in the Liver of Rats With Intrauterine Growth Retardation

BACKGROUND

Intrauterine growth restriction (IUGR) is highly associated with fetal as well as neonatal morbidity, mortality, and an increased risk metabolic disease development later in life. The mechanism involved in the increased risk has not been established. We compared differentially expressed genes between the liver of appropriate for gestational age (AGA) and IUGR rat models and identified their effects on molecular pathways involved in the metabolic syndrome.

METHODS

We extracted RNA from the liver of IUGR and AGA rats and profiled gene expression by microarray analysis. GO function and KEGG pathway enrichment analyses were conducted using the Search Tool for the Retrieval of Interacting Genes database. Then, the Cytoscape software was used to visualize regulatory interaction networks of IUGR-related genes. The results were further verified via quantitative reverse transcriptase PCR analysis.

RESULTS

In this study, 815 genes were found to be markedly differentially expressed (fold-change >1.5, p < 0.05) between IUGR and AGA, with 347 genes elevated and 468 suppressed in IUGR, relative to AGA. Enrichment and protein-protein interaction network analyses of target genes revealed that core genes including Ppargc1a, Prkaa2, Slc2a1, Rxrg, and Gcgr, and pathways, including the PPAR signaling pathway and FoxO signaling pathway, had a potential association with metabolic syndrome development in IUGR. We also confirmed that at the mRNA level, five genes involved in glycometabolism were differentially expressed between IUGR and AGA.

CONCLUSION

Our findings elucidate on differential gene expression profiles in IUGR and AGA. Moreover, they elucidate on the pathogenesis of IUGR-associated metabolic syndromes. The suggested candidates are potential biomarkers and eventually intended to treat them appropriately.

Role of hsa_circ_0066966 in proliferation and migration of hepatitis B virus-related liver cancer cells

A large proportion of liver cancer cases is caused by hepatitis B virus (HBV) infection. In recent years, an increasing number of reports have indicated that circular RNAs (circRNAs) exert regulatory effects in cancer development, whereas the role of circRNAs in HBV-positive liver cancer requires further investigation. In the present study, abnormally expressed circRNAs were identified in HBV-positive liver cancer cells through microarray analysis. A total of 1,493 differentially expressed circRNAs [absolute fold-change (FC) ≥2] in HBV-positive liver cancer cells were detected, of which 171 were upregulated and 1,322 were downregulated. Subsequently, Gene Ontology enrichment analysis indicated that the genes of dysregulated circRNAs were mainly involved in regulating Sertoli cell differentiation and development, as well as telomeric DNA binding. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that most of these genes were enriched in cancer-related signaling pathways, including the MAPK and Hippo signaling pathways. Next, the expression levels of the top-10 dysregulated circRNAs were verified in HBV-positive liver cancer cells through reverse transcription-quantitative PCR. Among them, hsa_circ_0066966 had the highest absolute Log₂FC value and was abnormally increased in HBV-positive liver cancer cells. Functional experiments further verified that knockdown of hsa_circ_0066966 had a significant inhibitory effect on the proliferation and migration of HBV-positive liver cancer cells. By contrast, overexpression of hsa_circ_0066966 in HBV-negative liver cancer cells resulted in the opposite effect. In conclusion, in the present study, comprehensive circRNA profiling in HBV-positive liver cancer cells indicated that hsa_circ_0066966 may regulate the progression of HBV-positive liver cancer.

Screening and Identification of Muscle-Specific Candidate Genes via Mouse Microarray Data Analysis

Muscle tissue is involved with every stage of life activities and has roles in biological processes. For example, the blood circulation system needs the heart muscle to transport blood to all parts, and the movement cannot be separated from the participation of skeletal muscle. However, the process of muscle development and the regulatory mechanisms of muscle development are not clear at present. In this study, we used bioinformatics techniques to identify differentially expressed genes specifically expressed in multiple muscle tissues of mice as potential candidate genes for studying the regulatory mechanisms of muscle development. Mouse tissue microarray data from 18 tissue samples was selected from the GEO database for analysis. Muscle tissue as the treatment group, and the other 17 tissues as the control group. Genes expressed in the muscle tissue were different to those in the other 17 tissues and identified 272 differential genes with highly specific expression in muscle tissue, including 260 up-regulated genes and 12 down regulated genes. is the genes were associated with the myofibril, contractile fibers, and sarcomere, cytoskeletal protein binding, and actin binding. KEGG pathway analysis showed that the differentially expressed genes in muscle tissue were mainly concentrated in pathways for AMPK signaling, cGMP PKG signaling calcium signaling, glycolysis, and, arginine and proline metabolism. A PPI protein interaction network was constructed for the selected differential genes, and the MCODE module used for modular analysis. Five modules with Score > 3.0 are selected. Then the Cytoscape software was used to analyze the tissue specificity of differential genes, and the genes with high degree scores collected, and some common genes selected for quantitative PCR verification. The conclusion is that we have screened the differentially expressed gene set specific to mouse muscle to provide potential candidate genes for the study of the important mechanisms of muscle development.