Transcriptomic Profiling of Buffalo Spermatozoa Reveals Dysregulation of Functionally Relevant mRNAs in Low-Fertile Bulls

Although, it is known that spermatozoa harbor a variety of RNAs that may influence embryonic development, little is understood about sperm transcriptomic differences in relation to fertility, especially in buffaloes. In the present study, we compared the differences in sperm functional attributes and transcriptomic profile between high- and low-fertile buffalo bulls. Sperm membrane and acrosomal integrity were lower (P < 0.05), while protamine deficiency and lipid peroxidation were higher (P < 0.05) in low- compared to high-fertile bulls. Transcriptomic analysis using mRNA microarray technology detected a total of 51,282 transcripts in buffalo spermatozoa, of which 4,050 transcripts were differentially expressed, and 709 transcripts were found to be significantly dysregulated (P < 0.05 and fold change >1) between high- and low-fertile bulls. Majority of the dysregulated transcripts were related to binding activity, transcription, translation, and metabolic processes with primary localization in the cell nucleus, nucleoplasm, and in cytosol. Pathways related to MAPK signaling, ribosome pathway, and oxidative phosphorylation were dysregulated in low-fertile bull spermatozoa. Using bioinformatics analysis, we observed that several genes related to sperm functional attributes were significantly downregulated in low-fertile bull spermatozoa. Validation of the results of microarray analysis was carried out using real-time qPCR expression analysis of selected genes (YBX1, ORAI3, and TFAP2C). The relative expression of these genes followed the same trend in both the techniques. Collectively, this is the first study to report the transcriptomic profile of buffalo spermatozoa and to demonstrate the dysregulation of functionally relevant transcripts in low-fertile bull spermatozoa. The results of the present study open up new avenues for understanding the etiology for poor fertility in buffalo bulls and to identify fertility biomarkers.

Oral squamous cell carcinoma may originate from bone marrow-derived stem cells

Molecules that demonstrate a clear association with the aggressiveness of oral squamous cell carcinoma (OSCC) have not yet been identified. The current study hypothesized that tumor cells in OSCC have three different origins: Epithelial stem cells, oral tissue stem cells from the salivary gland and bone marrow (BM) stem cells. It was also hypothesized that carcinomas derived from less-differentiated stem cells have a greater malignancy. In the present study, sex chromosome analysis by fluorescence in situ hybridization and/or microdissection PCR was performed in patients with OSCC that developed after hematopoietic stem cell transplantation (HSCT) from the opposite sex. OSCC from 3 male patients among the 6 total transplanted patients were considered to originate from donor-derived BM cells. A total of 2/3 patients had distant metastasis, resulting in a poor prognosis. In a female patient with oral potentially malignant disorder who underwent HSCT, there were 10.7% Y-containing cells in epithelial cells, suggesting that some epithelial cells were from the donor. Subsequently, gene expression patterns in patients with possible BM stem cell-derived OSCC were compared with those in patients with normally developed OSCC by microarray analysis. A total of 3 patients with BM stem cell-derived OSCC exhibited a specific pattern of gene expression. Following cluster analysis by the probes identified on BM stem cell-derived OSCC, 2 patients with normally developed OSCC were included in the cluster of BM stem cell-derived OSCC. If the genes that could discriminate the origin of OSCC were identified, OSCCs were classified into the three aforementioned categories. If diagnosis can be performed based on the origin of the cancer cells, a more specific therapeutic strategy may be implemented to improve prognosis. This would be a paradigm shift in diagnostic and therapeutic strategies for OSCC.

Chromosomal Microarray in Children With Developmental Delay: The Experience of a Tertiary Center in Korea

Background and Objectives: Chromosomal microarray (CMA) is a first-tier genetic test for children with developmental delay (DD), intellectual disability (ID), autism spectrum disorders (ASDs), and multiple congenital anomalies (MCA). In this study, we report our experiences with the use of CMA in Korean children with unexplained DD/ID. Methods: We performed CMA in a cohort of 308 children with DD/ID between January 2010 and September 2020. We also retrospectively reviewed their medical records. The Affymetrix CytoScan 750 K array with an average resolution of 100 kb was used to perform CMA. Results: Comorbid neurodevelopmental disorders were ASD (37 patients; 12.0%), epilepsy (34 patients; 11.0%), and attention deficit hyperactivity disorders (12 patients; 3.9%). The diagnostic yield was 18.5%. Among the 221 copy number variants (CNVs) identified, 70 CNVs (57 patients; 18.5%) were pathogenic. Deletion CNVs were more common among pathogenic CNVs (PCNVs) than in non-PCNVs (P < 0.001). The size difference between PCNVs and non-PCNVs was not significant (P = 0.023). The number of included genes within CNV intervals was significantly higher in PCNVs (average 8.6; 0-347) than in non-PCNVs (average 47.5; 1-386) (P < 0.001). Short stature and hearing difficulty were also more common in the PCNV group than in the non-PCNV group (P = 0.010 and 0.070, respectively). Conclusion: This study provides additional evidence for the usefulness of CMA in genetic testing of children with DD/ID in Korea. The pathogenicity of CNVs correlated with the number of included genes within the CNV interval and deletion type of the CNVs, but not with CNV size.

Upregulation of T Cell Receptor Signaling Pathway Components in Gestational Diabetes Mellitus Patients: Joint Analysis of mRNA and circRNA Expression Profiles

OBJECTIVE

Gestational diabetes mellitus (GDM) is one of the most common complications of pregnancy, and its pathogenesis is still unclear. Studies have shown that circular RNAs (circRNAs) can regulate blood glucose levels by targeting mRNAs, but the role of circRNAs in GDM is still unknown. Therefore, a joint microarray analysis of circRNAs and their target mRNAs in GDM patients and healthy pregnant women was carried out.

METHODS

In this study, microarray analyses of mRNA and circRNA in 6 GDM patients and 6 healthy controls were conducted to identify the differentially expressed mRNA and circRNA in GDM patients, and some of the discovered mRNAs and circRNAs were further validated in additional 56 samples by quantitative realtime PCR (qRT-PCR) and droplet digital PCR (ddPCR).

RESULTS

Gene ontology and pathway analyses showed that the differentially expressed genes were significantly enriched in T cell immune-related pathways. Cross matching of the differentially expressed mRNAs and circRNAs in the top 10 KEGG pathways identified 4 genes (CBLB, ITPR3, NFKBIA, and ICAM1) and 4 corresponding circRNAs (circ-CBLB, circ-ITPR3, circ-NFKBIA, and circ-ICAM1), and these candidates were subsequently verified in larger samples. These differentially expressed circRNAs and their linear transcript mRNAs were all related to the T cell receptor signaling pathway, and PCR results confirmed the initial microarray results. Moreover, circRNA/miRNA/mRNA interactions and circRNA-binding proteins were predicted, and circ-CBLB, circ-ITPR3, and circ-ICAM1 may serve as GDM-related miRNA sponges and regulate the expression of CBLB, ITPR3, NFKBIA, and ICAM1 in cellular immune pathways.

CONCLUSION

Upregulation of T cell receptor signaling pathway components may represent the major pathological mechanism underlying GDM, thus providing a potential approach for the prevention and treatment of GDM.

Aldehyde Dehydrogenase 1-related Genes in Triple-negative Breast Cancer Investigated Using Network Analysis

BACKGROUND/AIM

Aldehyde dehydrogenase 1 (ALDH1) is known as a breast cancer stem cell (CSC) marker. This study aimed to identify genes associated with ALDH1.

MATERIALS AND METHODS

ALDH1-positive and -negative breast cancer cells were isolated using laser capture microdissection from five tissue samples of ALDH1-positive breast cancer patients. Messenger RNA expression levels were compared between ALDH1-positive and -negative cells.

RESULTS

We found 104 differentially expressed genes between ALDH1-positive and -negative cells. Gene ontology and pathway analysis revealed that these genes were correlated with CSC functions and pathways. Network analyses identified 10 genes that were closely associated with ALDH1. We validated these 10 genes utilizing The Cancer Genome Atlas and the Molecular Taxonomy of Breast Cancer International Consortium cohort, and found that they were associated with ALDH1 expression and correlated with Wnt pathway signaling.

CONCLUSION

The 10 genes we identified could be potential targets for CSC therapy of breast cancer.

Genome-wide identification of Chiari malformation type I associated candidate genes and chromosomal variations

Chiari malformation type I (CMI) is a brain malformation that is characterized by herniation of the cerebellum into the spinal canal. Chiari malformation type I is highly heterogeneous; therefore, an accurate explanation of the pathogenesis of the disease is often not possible. Although some studies showed the role of genetics in CMI, the involvement of genetic variations in CMI pathogenesis has not been thoroughly elucidated. Therefore, in the current study we aim to reveal CMI-associated genomic variations in familial cases.Four CMI patients and 7 unaffected healthy members of two distinct families were analyzed. A microarray analysis of the affected and unaffected individuals from two Turkish families with CMI was conducted. Analyses of single nucleotide variations (SNVs) and copy number variations (CNVs) were performed by calculation of B allele frequency (BAF) and log R ratio (LRR) values from whole genome SNV data. Two missense variations, OLFML2A (rs7874348) and SLC4A9 (rs6860077), and a 5'UTR variation of COL4A1 (rs9521687) were significantly associated with CMI. Moreover, 12 SNVs in the intronic regions of FAM155A, NR3C1, TRPC7, ASTN2, and TRAF1 were determined to be associated with CMI. The CNV analysis showed that the 11p15.4 chromosome region is inherited in one of the families. The use of familial studies to explain the molecular pathogenesis of complex diseases such as CMI is crucial. It has been suggested that variations in OLFML2A, SLC4A9, and COL4A1 play a role in CMI molecular pathogenesis. The CNV analysis of individuals in both families revealed a potential chromosomal region, 11p15.4, and risk regions that are associated with CMI.

Preliminary verification of lncRNA ENST00000609755.1 potential ceRNA regulatory network in coronary heart disease

BACKGROUND

This study aims to explore the possible ceRNA regulatory network of lncRNA ENST00000609755.1 in CHD patients based on the population; reveal the possible regulatory mechanism of lncRNA ENST00000609755.1.

METHOD

Microarray analysis were used to identify differentially expressed miRNA, and mRNA profiles between 5 CHD and 5 healthy controls. The lncRNA ENST00000609755.1-miRNA-mRNA ceRNA regulatory network was constructed with lncRNA ENST00000609755.1 as the core based on microarray data and related prediction software (RNAhybird, miRanda, miRWalk 2.0). Furthermore, qRT-PCR was used to verify the expression levels of miRNA and mRNA. t-test and pearson correlation analysis were used to compare the expression differences and correlations of lncRNA, miRNA and mRNA. The receiver operating characteristic (ROC) curve was used to determine the discriminative ability of lncRNA ENST00000609755.1 and its downstream targets.

RESULTS

Totally 25 miRNAs and 953 mRNAs were differentially expressed between CHD and healthy control. The lncRNA ENST00000609755.1- miRNA- mRNA ceRNA regulatory network was constructed (5 miRNA and 58 mRNA). qRT-PCR results suggest that the expression of lncRNA ENST00000609755.1 and ELK1 were up-regulated in CHD group and positively correlated, the expression of miR-150 was down-regulated in CHD, which was negatively correlated with lncRNA ENST00000609755.1 and ELK1. The AUC was 0.777(95%CI, 0.659-0.895) when miRNA-150 and ELK1 was added, which was higher than that of lncRNA ENST00000609755.1 single indicator.

CONCLUSION

LncRNA ENST00000609755.1, miR-150 and ELK1 may have a potential ceRNA regulatory network relationship which could be considered to have a good combined diagnostic value for CHD. Also, preliminarily reveal the possible mechanism of lncRNA ENST00000609755.1 involved in CHD.

Differential expression profile of long noncoding RNAs in chronic HBV infection: New insights into pathogenesis

Increasing studies have revealed that long noncoding RNAs (lncRNAs) might play vital roles in the development and progression of various diseases including viral infectious diseases. However, the expression and biological functions of lncRNAs in chronic hepatitis B virus (HBV) infection remain largely unknown. Therefore, lncRNA microarray was performed to analyze the lncRNAs' and messenger RNAs' (mRNAs) expression profiles in liver tissues from patients with chronic HBV infection. Subsequently, a comprehensive bioinformatics analysis was conducted to investigate the potential functions of the differentially expressed genes. As a result, a total of 203 differentially expressed lncRNAs and 180 mRNAs were identified in chronic HBV infection. The expressions of five differentially expressed lncRNAs were further validated using quantitative real-time polymerase chain reaction. Gene ontology, pathway analysis, and gene set enrichment analysis revealed that differentially expressed lncRNAs might be mainly be involved in cytokine-cytokine receptor interaction and varied biotransformation processes, including fatty acid metabolism, amino acid metabolism, carbon metabolism, and drug metabolism. Additionally, coexpression networks between differentially expressed lncRNAs and mRNAs were constructed to reveal the hub regulator and analyze the functional pathways. This study provided an overview of lncRNA and mRNA expression in liver tissues from patients with chronic HBV infection. These differentially expressed lncRNAs might play crucial roles in the pathogenesis and progression of chronic HBV infection, which deserve further investigation.

A novel 1p33p32.2 deletion involving SCP2, ORC1, and DAB1 genes in a patient with craniofacial dysplasia, short stature, developmental delay, and leukoencephalopathy: A case report

INTRODUCTION

Microdeletion syndromes occur from deletion of 5Mb of a chromosome in approximately 5% of patients with unexplained intellectual disability. Interstitial microdeletions at bands 1p33 and 1p32.2 of the short arm of chromosome 1 are rare and have not been previously reported in relation to disease.

PATIENT CONCERNS

We present a case of a 39-month boy with Pierre Robin sequence, development delay/intellectual disability, growth retardation, short stature, leukoencephalopathy, craniofacial dysplasia, and speech delay. The child was referred to the Child health care department in October 2014 for his delayed language development and aggravated aggression.

DIAGNOSIS

Molecular diagnostic testing with G-band karyotyping was normal but clinical microarray analysis detected a 10 Mb microdeletion at 1p33p32.2.

INTERVENTIONS

The patient received rehabilitation.

OUTCOMES

Three candidate genes were pinpointed to the deleted area, including ORC1, SCP2, and DAB1. Phenotype-genotype analysis suggested that these three genes are likely to be responsible for the main phenotypes observed in the patient, such as microcephaly, growth retardation, short stature, leukoencephalopathy, and development delay/intellectual disability.

CONCLUSIONS

The spectrum of phenotypes this case presented with are likely to be caused by 1p33p32.2 deletion which could represent a new microdeletion syndrome.

Systematic analysis of miRNAs in patients with postmenopausal osteoporosis

MicroRNAs (miRNAs) represent RNA species found in serum. Many miRNAs were observed that were related to osteoporosis and osteopenia. However, expression and function analysis of miRNAs in postmenopausal osteoporosis (PMOP) remain unaddressed. We first compared the miRNA expression of blood samples in postmenopausal women with osteopenia or with osteoporosis via analysis of GSE64433. Bioinformatics analyses were conducted to get the key miRNAs and their functions and pathways. 331 miRNAs were being identified as differentially expressed miRNAs. Among these, 122 miRNA (36.86%) were up-regulated, and the remaining 209 miRNAs (63.14%) were down-regulated. 105 genes were predicted as the targets of these miRNAs. GO enrichment analysis results showed that the miRNAs mainly enriched in DNA binding, ATP binding, gene expression, regulation of the apoptotic process, chromatin binding, and protein kinase binding. KEGG enrichment analysis results demonstrated that the miRNAs mainly enriched in the TGF beta signaling pathway, wnt signaling pathway, JAK-STAT signaling pathway, and androgen receptor signaling pathway. This study identified the abundant differentially expressed miRNAs in the blood samples of postmenopausal women with osteopenia or with osteoporosis. This study may contribute to getting new diagnostic and therapeutic strategies for PMOP.

Integration of gene profile to explore the hub genes of lung adenocarcinoma: A quasi-experimental study

BACKGROUND

Lung cancer is a leading cause of morbidity diseases worldwide, but the key mechanisms of lung cancer remain elusive. This study aims to integrate of GSE 118370 and GSE 32863 profile and identify the key genes and pathway involved in human lung adenocarcinoma.

METHODS

R software (RStudio, Version info: R 3.2.3, Forrester, USA) were utilized to find the differentially expressed genes. All the differentially expressed genes were analyzed by gene ontology, kyoto encyclopedia of genes and genomes. Protein-protein interaction networks were constructed by STRING database and analyzed by Cytohubber and Module. The cancer genome atlas database was used to verification the expression of hub genes. Quantitative reverse transcription-PCR was used to verify the bio-information results.

RESULTS

Sixty-four lung adenocarcinoma and 64 adjacent normal tissues were used for integration analysis. Five hundred ninety-nine co-expression genes were locked. Biological processes mainly enriched in angiogenesis. Cellular component focused on extracellular exosome and molecular function aimed on protein disulfide isomerase activity. Cytohubber analysis showed that GNG11, FPR2, P4HB, PIK3R1, CDC20, ADCY4, TIMP1, IL6, CXC chemokine ligand (CXCL)12, and GAS6 acted as the hub genes during lung adenocarcinoma. Module analysis presented Chemokine signaling pathway was a key pathway. Quantitative reverse transcription-PCR showed that the expression level of GNG11, FPR2, PIK3R1, ADCY4, IL6, CXCL12, and GAS6 were significantly decreased and P4HB, CDC20 and TIMP1 were increased in human adenocarcinoma tissues (P < .05). The cancer genome atlas online analysis showed GNG11 was not associated with survival.

CONCLUSIONS

This study firstly reported GNG11 acting as a hub gene in adenocarcinoma. GNG11 could be used as a biomarker for human adenocarcinoma. Chemokine signaling pathway might play important roles in lung adenocarcinoma.

Differential Gene Expression in Brain and Liver Tissue of Wistar Rats after Rapid Eye Movement Sleep Deprivation

Sleep is essential for the survival of most living beings. Numerous researchers have identified a series of genes that are thought to regulate "sleep-state" or the "deprived state". As sleep has a significant effect on physiology, we believe that lack of total sleep, or particularly rapid eye movement (REM) sleep, for a prolonged period would have a profound impact on various body tissues. Therefore, using the microarray method, we sought to determine which genes and processes are affected in the brain and liver of rats following nine days of REM sleep deprivation. Our findings showed that REM sleep deprivation affected a total of 652 genes in the brain and 426 genes in the liver. Only 23 genes were affected commonly, 10 oppositely, and 13 similarly across brain and liver tissue. Our results suggest that nine-day REM sleep deprivation differentially affects genes and processes in the brain and liver of rats.

Y-Chromosome Genetic Analysis of Modern Polish Population

The study presents a full analysis of the Y-chromosome variability of the modern male Polish population. It is the first study of the Polish population to be conducted with such a large set of data (2,705 individuals), which includes genetic information from inhabitants of all voivodeships, i.e., the first administrative level, in the country and the vast majority of its counties, i.e., the second level. In addition, the available data were divided into clusters corresponding to more natural geographic regions. Genetic analysis included the estimation of F_ST distances, the visualization with the use of multidimensional scaling plots and analysis of molecular variance. Y-chromosome binary haplogroups were classified and visualized with the use of interpolation maps. Results showed that the level of differentiation within Polish population is quite low, but some differences were indicated. It was confirmed that the Polish population is characterized by a high degree of homogeneity, with only slight genetic differences being observed at the regional level. The use of regional clustering as an alternative to counties and voivodeships provided a more detailed view of the genetic structure of the population. Those regional differences identified in the present study highlighted the need for additional division of the population by cultural and ethnic criteria in such studies rather than just by geographical or administrative regionalization.

Circular RNA microarray expression profile and potential function of circ0005875 in clear cell renal cell carcinoma

Background: Circular RNAs (circRNAs), a novel class of endogenous noncoding RNAs, are involved in a variety of diseases, including several types of cancers. We hypothesized that circRNAs are involved in the tumorigenesis and development of clear cell renal cell carcinoma (ccRCC).

Methods: To verify our hypothesis, we explored the circRNA expression profiles in 4 pairs of ccRCC tissues and their adjacent non-carcinoma tissues via microarray analysis. Selected circRNAs were further validated by qPCR. Moreover, hsa_circ_0005875 was selected for further study and the potential clinical values of hsa_circ_0005875 were investigated in 60 pairs of ccRCC tissues and adjacent normal controls. In addition, the role of hsa_circ_0005875 in ccRCC progression were performed using colony formation assay, Transwell assay and Martrigel-Transwell assay respectively. Finally, interactions between the circRNAs and miRNAs were predicted using Arraystar's miRNA target prediction software. Luciferase reporter assays were performed to evaluate the interaction between hsa_circ_0005875 and hsa_miR-145-5p.

Results: The microarray data showed 1988 circRNAs were significantly dysregulated circRNAs, including 1033 upregulated and 955 downregulated ones in the ccRCC tissues. Hsa_circ_0005875 was confirmed to be significantly upregulated in the ccRCC tumor tissues and renal carcinoma cells. Further analysis revealed that hsa_circ_0005875 expression was associated with tumor size, pathological TNM stage, histological differentiation, and lymphatic metastasis. Functional experiments demonstrated that overexpression of hsa_circ_0005875 increased proliferation, migration and invasion abilities. Moreover, bioinformatics analysis and luciferase reporter assays suggest that hsa_circ_0005875 may serve as a ceRNA (competing endogenous RNA) of miR-145-5p to relieve the repressive effect of miR-145-5p on target ZEB2.

Conclusions: These data indicate that hsa_circ_0005875 might play a role in promoting tumor growth and metastasis and be a potential biomarker of ccRCC.

Tocotrienol Rich Fraction Supplementation Modulate Brain Hippocampal Gene Expression in APPswe/PS1dE9 Alzheimer's Disease Mouse Model

Alzheimer’s disease (AD) is a progressive neurodegenerative disorder characterized by loss of memory and other cognitive abilities. AD is associated with aggregation of amyloid-β (Aβ) deposited in the hippocampal brain region. Our previous work has shown that tocotrienol rich fraction (TRF) supplementation was able to attenuate the blood oxidative status, improve behavior, and reduce fibrillary-type Aβ deposition in the hippocampus of an AD mouse model. In the present study, we investigate the effect of 6 months of TRF supplementation on transcriptome profile in the hippocampus of APPswe/PS1dE9 double transgenic mice. TRF supplementation can alleviate AD conditions by modulating several important genes in AD. Moreover, TRF supplementation attenuated the affected biological process and pathways that were upregulated in the AD mouse model. Our findings indicate that TRF supplementation can modulate hippocampal gene expression as well as biological processes that can potentially delay the progression of AD.

Genome-Scale Expression Pattern of Long Non-Coding RNAs in Chinese Uyghur Patients with Parkinson's Disease

BACKGROUND Long non-coding RNAs (lncRNAs) are transcripts thought to regulate gene expression at the post-transcriptional level. Some lncRNAs are associated with Parkinson's disease (PD) and participate in pathological processes of PD. The incidence of PD is relatively high in members of the Uyghur minority living in Xingjiang province of China. This study measured the expression of lncRNAs in the peripheral blood cells of Chinese Uyghur individuals with and without PD and analyzed the possible function of these lncRNAs in the development of PD. MATERIAL AND METHODS Peripheral blood samples were collected from 55 Uyghur patients with PD and 55 healthy volunteers. Total RNA was extracted, and the levels of expression of whole-genome lncRNAs and mRNAs in 10 samples (5 PD and 5 controls) were determined by microarray method. The expression levels of lncRNAs in all 100 subjects were determined by qRT-PCR. The lncRNA expression profiles of PD patients were determined based on lncRNA microarray chip analysis, and differentially expressed lncRNAs were identified. The results of chip analysis were confirmed in a large clinical cohort. RESULTS Comparison of subjects with and without PD identified 32 significantly up-regulated and 18 significantly down-regulated lncRNAs in the PD group. GO analysis showed that mRNAs encoding proteins involved in the regulation of biological processes were differentially expressed, with the inflammatory immune response being the most significantly related pathway. CONCLUSIONS The expression of lncRNAs in peripheral blood differed significantly in PD patients and controls. These differentially expressed lncRNAs may play a role in the development of PD.

Identification and analysis of key lncRNAs in malignant-transformed BEAS-2B cells induced with coal tar pitch by microarray analysis

This study aims to explore the key and differentially expressed long non-coding RNAs (lncRNAs) and elucidates their possible mechanisms in malignant-transformed Human bronchial epithelial (BEAS-2B) cells induced by coal tar pitch extracts (CTPE). BEAS-2B cells were stimulated with 2.4 μg/ml CTPE, then passaged for three times which were named CTPE1 and then passaged until passage 30 (CTPE30). The results showed that cells of CTPE30 appeared abnormal morphology. Furthermore, migration, clonality and proliferation of cells in CTPE group were significantly increased compared with those in control groups. However, the apoptosis of cells in CTPE group was inhibited. A total of 569 differentially expressed mRNAs and 707 differentially expressed lncRNAs were screened out, among which four lncRNAs were validated and were consistent with the microarray results. 32 target genes were screened out by Co-expression network. The study suggests that differentially expressed lncRNAs may play a potential role in lung carcinogenesis.

Differential DNA methylation profile in infants born small-for-gestational-age: association with markers of adiposity and insulin resistance from birth to age 24 months

INTRODUCTION

Prenatal growth restraint followed by rapid postnatal weight gain increases lifelong diabetes risk. Epigenetic dysregulation in critical windows could exert long-term effects on metabolism and confer such risk.

RESEARCH DESIGN AND METHODS

We conducted a genome-wide DNA methylation profiling in peripheral blood from infants born appropriate-for-gestational-age (AGA, n=30) or small-for-gestational-age (SGA, n=21, with postnatal catch-up) at age 12 months, to identify new genes that may predispose to metabolic dysfunction. Candidate genes were validated by bisulfite pyrosequencing in the entire cohort. All infants were followed since birth; cord blood methylation profiling was previously reported. Endocrine-metabolic variables and body composition (dual-energy X-ray absorptiometry) were assessed at birth and at 12 and 24 months.

RESULTS

GPR120 (cg14582356, cg01272400, cg23654127, cg03629447), NKX6.1 (cg22598426, cg07688460, cg17444738, cg12076463, cg10457539), CPT1A (cg14073497, cg00941258, cg12778395) and IGFBP 4 (cg15471812) genes were hypermethylated (GPR120, NKX6.1 were also hypermethylated in cord blood), whereas CHGA (cg13332653, cg15480367, cg05700406), FABP5 (cg00696973, cg10563714, cg16128701), CTRP1 (cg19231170, cg19472078, cg0164309, cg07162665, cg17758081, cg18996910, cg06709009), GAS6 (N/A), ONECUT1 (cg14217069, cg02061705, cg26158897, cg06657050, cg15446043) and SLC2A8 (cg20758474, cg19021975, cg11312566, cg12281690, cg04016166, cg03804985) genes were hypomethylated in SGA infants. These genes were related to β-cell development and function, inflammation, and glucose and lipid metabolism and associated with body mass index, body composition, and markers of insulin resistance at 12 and 24 months.

CONCLUSION

In conclusion, at 12 months, abnormal methylation of GPR120 and NKX6.1 persists and new epigenetic marks further involved in adipogenesis and energy homeostasis arise in SGA infants. These abnormalities may contribute to metabolic dysfunction and diabetes risk later in life.

Identification of potential gene drivers of cutaneous squamous cell carcinoma: Analysis of microarray data

Cutaneous squamous cell carcinoma (cSCC) is a common skin cancer with an increasing incidence. As a pre-cancerous condition, actinic keratosis (AK) has an up to 20% risk of progression to cSCC. This study aims to define the potential genes that associated with genesis and progression of cSCC, thereby further identify critical biomarkers for the prevention, early diagnosis, and effective treatment of cSCC.Two datasets GSE42677 and GSE45216 were downloaded from the GEO. Microarray data analysis was applied to explore the differentially expressed genes (DEGs) between cSCC samples and AK samples. Then functional enrichment analysis, protein-protein interaction (PPI) network, and drug-gene interaction analysis were performed to screen key genes.A total of 711 DEGs, including 238 upregulated genes and 473 downregulated genes, were screened out. DEGs mainly involved in pathways as extracellular matrix (ECM)-receptor interaction, hematopoietic cell lineage, phosphatidylinositol 3-kinase (PI3K-Akt) signaling pathway, and focal adhesion. Candidate genes, including upregulated genes as JUN, filamin A (FLNA), casein kinase 1 delta (CSNK1D), and histone cluster 1 H3 family member f (HIST1H3F), and downregulated genes as androgen receptor (AR), heat shock protein family H member 1 (HSPH1), tropomyosin 1 (TPM1), pyruvate kinase, muscle (PKM), LIM domain and actin binding 1 (LIMA1), and synaptopodin (SYNPO) were screened out. In drug-gene interaction analysis, 13 genes and 44 drugs were identified.This study demonstrates that genes JUN, FLNA, AR, HSPH1, and CSNK1D have the potential to function as targets for diagnosis and treatment of cSCC.

In silico prediction and validation of potential therapeutic genes in pancreatic β-cells associated with type 2 diabetes

Diabetes mellitus is becoming a major health burden worldwide. Pancreatic β-cell death is a characteristic of type 2 diabetes (T2D), but the underlying mechanisms of pancreatic β-cell death remain unknown. Therefore, the aim of the present study was to identify potential targets in the pancreatic islet of T2D. The GSE20966 dataset was obtained from the Gene Expression Omnibus (GEO) database, and differentially expressed genes (DEGs) were identified by using the GEO2R tool. The Gene Ontology terms and Kyoto Encyclopedia of Genes and Genomes Pathway enrichment analysis of DEGs were further assessed using the Database for Annotation, Visualization and Integrated Discovery. Furthermore, protein-protein interaction (PPI) networks were constructed for the up- and downregulated genes using STRING databases and were then visualized with Cytoscape. The body weight, fasting blood glucose (FBG), pancreatic index and biochemistry parameters were measured in db/db mice. Moreover, the morphology of the pancreas was detected by hematoxylin and eosin staining, and hub genes were assessed using reverse transcription-quantitative PCR (RT-qPCR) and western blot analysis. In total, 570 DEGs were screened, including 376 upregulated and 194 downregulated genes, which were associated with 'complement activation, classical pathway', 'proteolysis', 'complement activation' and 'pancreatic secretion pathway'. It was found that the body weight, FBG, alanine aminotransferase, aspartate aminotransferase, total cholesterol, triglycerides, blood urea nitrogen, creatinine, fasting serum insulin, glucagon and low-density lipoprotein cholesterol levels were significantly higher in db/db mice, while high-density lipoprotein cholesterol levels and the pancreatic index were significantly decreased. Furthermore, albumin, interleukin-8, CD44, C-C motif chemokine ligand 2, hepatocyte growth factor, cystic fibrosis transmembrane conductance regulator, histone cluster 1 H2B family member n, mitogen-activated protein kinase 11 and neurotrophic receptor tyrosine kinase 2 were identified as hub genes in PPI network. RT-qPCR and western blotting results demonstrated the same expression trend in hub genes as found by the bioinformatics analysis. Therefore, the present study identified a series of hub genes involved in the progression of pancreatic β-cell, which may help to develop effective therapeutic strategy for T2D.

Expression profile and bioinformatics analysis of circular RNA in intestinal mucosal injury and repair after severe burns

Circular RNA (circRNA) is a novel noncoding RNA that is mostly found in humans and animals. Although the flux of circRNA research has increased in recent years, its precise function is still unclear. Some studies demonstrate that circRNAs can function as microRNA (miRNA) sponges involved in the regulation of competitive endogenous RNAs networks and play a crucial role in many biological processes. Other studies show that circRNAs play multiple biological roles in gastrointestinal diseases. However, the expression characteristics and function of circRNA in intestinal mucosal injury and repair after severe burn have not been reported. This study aims to screen differentially expressed circRNAs in intestinal mucosal injury and repair after severe burns and understand their underlying mechanisms. To test our hypothesis that circRNA may play a role in promoting repair in intestinal mucosa injury after severe burns, we collected the intestinal tissues of three severely burned mice and three pseudo-scalded mice and evaluated the expression of circRNAs via microarray analysis. Quantitative real-time polymerase chain reaction was also used to validate the circRNA microarray data by selecting six based on different multiples, original values, and p values. The host genes of all differentially expressed circRNAs and the downstream target genes of six selected DEcircRNAs were identified by Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway analysis. Meanwhile, we also created a circRNA-miRNA-mRNA network to predict the role and function of circRNAs in intestinal mucosal injury and repair after severe burns.

Differential gene expression in patients with primary mitral valve disease: identifying potential therapeutic targets in the era of precision medicine

Primary (degenerative) mitral valve (MV) disease is a result of structural remodeling due to degenerative and adaptive changes of MV tissue. We hypothesized that in patients with primary MV disease undergoing surgery for severe mitral regurgitation (MR), a distinct genetic expression profile within the MV leaflet tissue could be identified as compared with patients without MV disease. Tissue samples from the MV leaflets of 65 patients undergoing MV surgery for MR due to primary MV disease and 4 control cadavers without MV disease were collected and analyzed. MicroRNA transcripts were hybridized to Illumina HumanHT-12 v4 Beadchips. Ingenuity pathway analyses (IPAs) were conducted to provide biological interpretation. Of the approximately 20 000 genes examined, 4092 (20%) were differentially expressed between patients with primary MV disease and normal controls (false discovery rate<0.05). The differentially expressed genes could be clustered into five regulator effect networks from the Ingenuity Knowledge IPA database with a consistency score of >6. These five networks have been previously implicated in pathophysiological cardiac abnormalities, including inhibited contractility of the heart and fatty acid oxidation as well as activation of apoptosis of smooth muscle cells, cardiac degeneration, and hypertrophy of cardiac cells. MV tissue in patients with primary MV disease demonstrated distinct genetic expression patterns as compared with normal controls. Further studies are necessary to determine whether the molecular pathways identified in this experiment may represent potential therapeutic targets to prevent degeneration of MV tissue leading to severe MR.

Comprehensive analysis of circular RNA expression profiles in cisplatin-resistant non-small cell lung cancer cell lines

Platinum-based drugs such as cisplatin are widely used in combination chemotherapy for non-small cell lung cancer (NSCLC) owing to their high clinical response rate; however, acquired resistance to cisplatin is eventually inevitable. Circular RNAs (circRNAs) are involved in the development of diverse types of cancers, but their connection to cisplatin-resistance in NSCLC has not been studied. In the present study, two cisplatin-resistant NSCLC cell lines (A549/DDP and PC9/DDP) were established by gradually increasing concentrations of cisplatin in the media. The resulting cell lines possessed high resistance to cisplatin and strong proliferation, migration, and colony formation abilities compared to the parental cells. Microarray analysis identified 19,161 circRNAs that were dysregulated in cisplatin-resistant cell lines (fold change abs>2), including 11,915 up-regulated and 7246 down-regulated circRNAs. The expression of the top five up-regulated and down-regulated circRNAs was validated using real-time quantitative polymerase chain reaction. A circRNA-micro RNA (miRNA) network of the top 20 dysregulated circRNAs and their predicted miRNAs was constructed using Cytoscape. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses revealed that the host genes of the identified circRNAs were involved in the regulation of MAP kinase kinase kinase kinase activity, 6-phosphofructo-2-kinase activity, focal adhesion, ErbB signaling, and ECM-receptor interactions, which may contribute to cisplatin resistance in NSCLC. In summary, this is the first report on circRNA profiling in cisplatin-resistant NSCLC cells and it provides new potential targets for the reversal of cisplatin resistance in NSCLC.

Dysregulated circulating SOCS3 and haptoglobin expression associated with stable coronary artery disease and acute coronary syndrome: An integrated study based on bioinformatics analysis and case-control validation

OBJECTIVE

To extensively use blood transcriptome analysis to identify potential diagnostic and therapeutic targets for cardiovascular diseases.

METHODS

Two gene expression datasets (GSE59867 and GSE62646) were downloaded from GEO DataSets to identify altered blood transcriptomes in patients with ST-segment elevation myocardial infarction (STEMI) compared to stable coronary artery disease (CAD). Thereafter, several computational approaches were taken to determine functional roles and regulatory networks of differentially expressed genes (DEGs). Finally, the expression of dysregulated two hub genes-suppressor of cytokine signaling 3 (SOCS3) and haptoglobin (HP)-were validated in a case-control study.

RESULTS

A total of 119 DEGs were identified in the discovery phase, consisting of 71 downregulated genes and 48 upregulated genes; two hub modules consisting of two hub genes-SOCS3 and HP-were identified. In the validation phase, both SOCS3 and HP were significantly downregulated in the stable CAD and acute coronary syndrome (ACS) patients when compared with healthy controls. Meanwhile, HP was significantly upregulated in STEMI patients when compared with stable CAD patients (p=0.041). Logistic regression analysis indicated that: downregulated expression of HP correlated with increased risk of CAD [odds ratio (OR)=0.52, 95% confidence interval (CI)=0.31~0.87, p=0.013]; and downregulated expression of SOCS3 correlated with increased risk of ACS (OR=0.66, 95% CI=0.46~0.94, p=0.023) when age, gender, history of hyperlipidemia, diabetes and hypertension were included as covariates.

CONCLUSION

Future clarification of how SOCS3 and HP influence the pathogenesis of disease may pave the way for the development of novel diagnostic and therapeutic methods.

Investigation of hub genes involved in diabetic nephropathy using biological informatics methods

Background

The aim of this study was to find genes with significantly aberrant expression in diabetic nephropathy (DN) and determine their underlying mechanisms.

Methods

GSE30528 and GSE1009 were obtained by querying the Gene Expression Omnibus (GEO) database. The difference in target gene expression between normal renal tissues and kidney tissues in patients with DN was screened by using the GEO2R tool. Using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database, differentially expressed genes (DEGs) were analysed by Gene Ontology (GO) annotation and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. Then, the protein-protein interactions (PPIs) of DEGs were analyzed by Cytoscape with the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database, and the hub genes in this PPI network were recognized by centrality analysis.

Results

There were 110 genes with significant expression differences between normal and DN tissues. The differences in gene expression involved many functions and expression pathways, such as the formation of the extracellular matrix and the construction of the extracellular domain. The correlation analysis and subgroup analysis of 14 hub genes and the clinical characteristics of DN showed that CTGF, ALB, PDPN, FLT1, IGF1, WT1, GJA1, IGFBP2, FGF9, BMP2, FGF1, BMP7, VEGFA, and TGFBR3 may be involved in the progression of DN.

Conclusions

We confirmed the differentially expressed hub genes and other genes which may be the novel biomarker and target candidates in DN.