Microarray analyses of peripheral blood cells identifies unique gene expression signature in psoriatic arthritis

Psoriatic arthritis (PsA) is a chronic and erosive form of arthritis of unknown cause. We aimed to characterize the PsA phenotype using gene expression profiling and comparing it with healthy control subjects and patients rheumatoid arthritis (RA). Peripheral blood cells (PBCs) of 19 patients with active PsA and 19 age- and sex-matched control subjects were used in the analyses of PsA, with blood samples collected in PaxGene tubes. A significant alteration in the pattern of expression of 313 genes was noted in the PBCs of PsA patients on Affymetrix U133A arrays: 257 genes were expressed at reduced levels in PsA, and 56 genes were expressed at increased levels, compared with controls. Downregulated genes tended to cluster to certain chromosomal regions, including those containing the psoriasis susceptibility loci PSORS1 and PSORS2. Among the genes with the most significantly reduced expression were those involved in downregulation or suppression of innate and acquired immune responses, such as SIGIRR, STAT3, SHP1, IKBKB, IL-11RA, and TCF7, suggesting inappropriate control that favors proin-flammatory responses. Several members of the MAPK signaling pathway and tumor suppressor genes showed reduced expression. Three proinflammatory genes--S100A8, S100A12, and thioredoxin--showed increased expression. Logistic regression and recursive partitioning analysis determined that one gene, nucleoporin 62 kDa, could correctly classify all controls and 94.7% of the PsA patients. Using a dataset of 48 RA samples for comparison, the combination of two genes, MAP3K3 followed by CACNA1S, was enough to correctly classify all RA and PsA patients. Thus, PBC gene expression profiling identified a gene expression signature that differentiated PsA from RA, and PsA from controls. Several novel genes were differentially expressed in PsA and may prove to be diagnostic biomarkers or serve as new targets for the development of therapies.     

Interval-valued analysis for discriminative gene selection and tissue sample classification using microarray data

An important application of gene expression data is to classify samples in a variety of diagnostic fields. However, high dimensionality and a small number of noisy samples pose significant challenges to existing classification methods. Focused on the problems of overfitting and sensitivity to noise of the dataset in the classification of microarray data, we propose an interval-valued analysis method based on a rough set technique to select discriminative genes and to use these genes to classify tissue samples of microarray data. We first select a small subset of genes based on interval-valued rough set by considering the preference-ordered domains of the gene expression data, and then classify test samples into certain classes with a term of similar degree. Experiments show that the proposed method is able to reach high prediction accuracies with a small number of selected genes and its performance is robust to noise.     

MicroRNAs hsa-miR-99b,hsa-miR-330, hsa-miR-126 and hsa-miR-30c: Potential Diagnostic Biomarkers in Natural Killer (NK) Cells of Patients with Chronic Fatigue Syndrome (CFS)/ Myalgic Encephalomyelitis (ME)

Background

Chronic Fatigue Syndrome (CFS/ME) is a complex multisystem disease of unknown aetiology which causes debilitating symptoms in up to 1% of the global population. Although a large cohort of genes have been shown to exhibit altered expression in CFS/ME patients, it is currently unknown whether microRNA (miRNA) molecules which regulate gene translation contribute to disease pathogenesis. We hypothesized that changes in microRNA expression in patient leukocytes contribute to CFS/ME pathology, and may therefore represent useful diagnostic biomarkers that can be detected in the peripheral blood of CFS/ME patients.

Methods

miRNA expression in peripheral blood mononuclear cells (PBMC) from CFS/ME patients and healthy controls was analysed using the Ambion Bioarray V1. miRNA demonstrating differential expression were validated by qRT-PCR and then replicated in fractionated blood leukocyte subsets from an independent patient cohort. The CFS/ME associated miRNA identified by these experiments were then transfected into primary NK cells and gene expression analyses conducted to identify their gene targets.

Results

Microarray analysis identified differential expression of 34 miRNA, all of which were up-regulated. Four of the 34 miRNA had confirmed expression changes by qRT-PCR. Fractionating PBMC samples by cell type from an independent patient cohort identified changes in miRNA expression in NK-cells, B-cells and monocytes with the most significant abnormalities occurring in NK cells. Transfecting primary NK cells with hsa-miR-99b or hsa-miR-330-3p, resulted in gene expression changes consistent with NK cell activation but diminished cytotoxicity, suggesting that defective NK cell function contributes to CFS/ME pathology.

Conclusion

This study demonstrates altered microRNA expression in the peripheral blood mononuclear cells of CFS/ME patients, which are potential diagnostic biomarkers. The greatest degree of miRNA deregulation was identified in NK cells with targets consistent with cellular activation and altered effector function.     

基于ceRNA网络识别登革热诊断标志物

通过比较登革热患者和健康人群转录组数据,识别差异基因,构建失调ceRNA网络,筛选关键基因富集分析,解析潜在生物学功能,助力登革热诊断标志物的研究。从GEO数据库下载登革热外周血芯片数据,识别差异基因并进行富集分析。结合miRNA-mRNA互作数据,利用超几何算法和皮尔森相关性计算方法识别登革热失调ceRNA互作对,使用Cytoscape软件可视化ceRNA网络与模块挖掘,对网络模块进行功能富集及外部数据验证表达模式。筛选出251个差异基因,发现其富集在细胞周期等生物学通路中。经外部数据验证,网络模块基因的表达趋势与训练集数据大致相同,表明模块基因在登革热疾病中的潜在诊断效能。本研究可为确定有效的疾病诊断分子标志物提供思路。     

EEG spectral coherence data distinguish chronic fatigue syndrome patients from healthy controls and depressed patients-A case control study

Background  

Previous studies suggest central nervous system involvement in chronic fatigue syndrome (CFS), yet there are no established diagnostic criteria. CFS may be difficult to differentiate from clinical depression. The study's objective was to determine if spectral coherence, a computational derivative of spectral analysis of the electroencephalogram (EEG), could distinguish patients with CFS from healthy control subjects and not erroneously classify depressed patients as having CFS.     

Analysis of developmental gene conservation in the Actinomycetales using DNA/DNA microarray comparisons

Based on available genome sequences, Actinomycetales show significant gene synteny across a wide range of species and genera. In addition, many genera show varying degrees of complex morphological development. Using the presence of gene synteny as a basis, it is clear that an analysis of gene conservation across the Streptomyces and various other Actinomycetales will provide information on both the importance of genes and gene clusters and the evolution of morphogenesis in these bacteria. Genome sequencing, although becoming cheaper, is still relatively expensive for comparing large numbers of strains. Thus, a heterologous DNA/DNA microarray hybridization dataset based on a Streptomyces coelicolor microarray allows a cheaper and greater depth of analysis of gene conservation. This study, using both bioinformatical and microarray approaches, was able to classify genes previously identified as involved in morphogenesis in Streptomyces into various subgroups in terms of conservation across species and genera. This will allow the targeting of genes for further study based on their importance at the species level and at higher evolutionary levels.     

Reporter Gene Silencing in Targeted Mouse Mutants Is Associated with Promoter CpG Island Methylation

Targeted mutations in mouse disrupt local chromatin structure and may lead to unanticipated local effects. We evaluated targeted gene promoter silencing in a group of six mutants carrying the tm1a Knockout Mouse Project allele containing both a LacZ reporter gene driven by the native promoter and a neo selection cassette. Messenger RNA levels of the reporter gene and targeted gene were assessed by qRT-PCR, and methylation of the promoter CpG islands and LacZ coding sequence were evaluated by sequencing of bisulfite-treated DNA. Mutants were stratified by LacZ staining into presumed Silenced and Expressed reporter genes. Silenced mutants had reduced relative quantities LacZ mRNA and greater CpG Island methylation compared with the Expressed mutant group. Within the silenced group, LacZ coding sequence methylation was significantly and positively correlated with CpG Island methylation, while promoter CpG methylation was only weakly correlated with LacZ gene mRNA. The results support the conclusion that there is promoter silencing in a subset of mutants carrying the tm1a allele. The features of targeted genes which promote local silencing when targeted remain unknown.     

Predicting hepatitis B virus-positive metastatic hepatocellular carcinomas using gene expression profiling and supervised machine learning

Hepatocellular carcinoma (HCC) is one of the most common and aggressive human malignancies. Its high mortality rate is mainly a result of intra-hepatic metastases. We analyzed the expression profiles of HCC samples without or with intra-hepatic metastases. Using a supervised machine-learning algorithm, we generated for the first time a molecular signature that can classify metastatic HCC patients and identified genes that were relevant to metastasis and patient survival. We found that the gene expression signature of primary HCCs with accompanying metastasis was very similar to that of their corresponding metastases, implying that genes favoring metastasis progression were initiated in the primary tumors. Osteopontin, which was identified as a lead gene in the signature, was over-expressed in metastatic HCC; an osteopontin-specific antibody effectively blocked HCC cell invasion in vitro and inhibited pulmonary metastasis of HCC cells in nude mice. Thus, osteopontin acts as both a diagnostic marker and a potential therapeutic target for metastatic HCC.     

DMD and IL1RAPL1: two large adjacent genes localized within a common fragile site (FRAXC) have reduced expression in cultured brain tumors

Common fragile sites (CFSs) are large regions of profound genomic instability found in all individuals. Spanning the center of the two most frequently expressed CFS regions, FRA3B (3p14.3) and FRA16D (16q23.2), are the 1.5 Mb FHIT gene and the 1.0 Mb WWOX gene. These genes are frequently deleted and/or altered in many different cancers. Both FHIT and WWOX have been demonstrated to function as tumor suppressors, both in vitro and in vivo. A number of other large CFS genes have been identified and are also frequently inactivated in multiple cancers. Based on these data, several additional very large genes were tested to determine if they were derived from within CFS regions, but DCC and RAD51L1 were not. However, the 2.0 Mb DMD gene and its immediately distal neighbor, the 1.8 Mb IL1RAPL1 gene are CFS genes contained within the FRAXC CFS region (Xp21.2-->p21.1). They are abundantly expressed in normal brain but were dramatically underexpressed in every brain tumor cell line and xenograft (derived from an intracranial model of glioblastoma multiforme) examined. We studied the expression of eleven other large CFS genes in the same panel of brain tumor cell lines and xenografts and found reduced expression of multiple large CFS genes in these samples. In this report we show that there is selective loss of specific large CFS genes in different cancers that does not appear to be mediated by the relative instability within different CFS regions. Further, the inactivation of multiple large CFS genes in xenografts and brain tumor cell lines may help to explain why this type of cancer is highly aggressive and associated with a poor clinical outcome.     

Integrated analysis of DNA methylation and gene expression profiles identified S100A9 as a potential biomarker in ulcerative colitis

Ulcerative colitis (UC) is a prevalent relapsing-remitting inflammatory bowel disease whose pathogenetic mechanisms remain elusive. In the present study, colonic biopsies samples from three UC patients treated in the Traditional Chinese Medicine Hospital and three healthy controls were obtained. The genome-wide mRNA and lncRNA expression of the samples were profiled through Agilent gene expression microarray. Moreover, the genome-wide DNA methylation dataset of normal and UC colon tissues was also downloaded from GEO for a collaborative analysis. Differential expression of lncRNA (DELs) and mRNAs (DEMs) in UC samples compared with healthy samples were identified by using limma Bioconductor package. Differentially methylated promoters (DMPs) in UC samples compared with controls were obtained through comparing the average methylation level of CpGs located at promoters by using t-test. Functional enrichment analysis was performed by the DAVID. STRING database was applied to the construction of gene functional interaction network. As a result, 2090 DEMs and 1242 DELs were screened out in UC samples that were closely associated with processes related to complement and coagulation cascades, osteoclast differentiation vaccinia, and hemorrhagic diseases. A total of 90 DEMs and 72 DELs were retained for the construction of functional network for the promoters of their corresponding genes were identified as DMPs. S100A9, HECW2, SOD3 and HIX0114733 showed high interaction degrees in the functional network, and expression of S100A9 was confirmed to be significantly elevated in colon tissues of UC patients compared with that of controls by qRT-PCR that was consistent with gene microarray analysis. These indicate that S100A9 could potentially be used as predictive biomarkers in UC.     

Extracting three-way gene interactions from microarray data

MOTIVATION: It is an important and difficult task to extract gene network information from high-throughput genomic data. A common approach is to cluster genes using pairwise correlation as a distance metric. However, pairwise correlation is clearly too simplistic to describe the complex relationships among real genes since co-expression relationships are often restricted to a specific set of biological conditions/processes. In this study, we described a three-way gene interaction model that captures the dynamic nature of co-expression relationship between a gene pair through the introduction of a controller gene. RESULTS: We surveyed 0.4 billion possible three-way interactions among 1000 genes in a microarray dataset containing 678 human cancer samples. To test the reproducibility and statistical significance of our results, we randomly split the samples into a training set and a testing set. We found that the gene triplets with the strongest interactions (i.e. with the smallest P-values from appropriate statistical tests) in the training set also had the strongest interactions in the testing set. A distinctive pattern of three-way interaction emerged from these gene triplets: depending on the third gene being expressed or not, the remaining two genes can be either co-expressed or mutually exclusive (i.e. expression of either one of them would repress the other). Such three-way interactions can exist without apparent pairwise correlations. The identified three-way interactions may constitute candidates for further experimentation using techniques such as RNA interference, so that novel gene network or pathways could be identified.     

Five candidate biomarkers associated with the diagnosis and prognosis of cervical cancer

Purpose: Cervical cancer (CC) is one of the most general gynecological malignancies and is associated with high morbidity and mortality. We aimed to select candidate genes related to the diagnosis and prognosis of CC.Methods: The mRNA expression profile datasets were downloaded. We also downloaded RNA-sequencing gene expression data and related clinical materials from TCGA, which included 307 CC samples and 3 normal samples. Differentially expressed genes (DEGs) were obtained by R software. GO function analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of DEGs were performed in the DAVID dataset. Using machine learning, the optimal diagnostic mRNA biomarkers for CC were identified. We used qRT-PCR and Human Protein Atlas (HPA) database to exhibit the differences in gene and protein levels of candidate genes.Results: A total of 313 DEGs were screened from the microarray expression profile datasets. DNA methyltransferase 1 (DNMT1), Chromatin Assembly Factor 1, subunit B (CHAF1B), Chromatin Assembly Factor 1, subunit A (CHAF1A), MCM2, CDKN2A were identified as optimal diagnostic mRNA biomarkers for CC. Additionally, the GEPIA database showed that the DNMT1, CHAF1B, CHAF1A, MCM2 and CDKN2A were associated with the poor survival of CC patients. HPA database and qRT-PCR confirmed that these genes were highly expressed in CC tissues.Conclusion: The present study identified five DEmRNAs, including DNMT1, CHAF1B, CHAF1A, MCM2 and Kinetochore-related protein 1 (KNTC1), as potential diagnostic and prognostic biomarkers of CC.     

Identification of key biomarkers for thyroid cancer by integrative gene expression profiles

Thyroid cancer is a frequently diagnosed malignancy and the incidence has been increased rapidly in recent years. Despite the favorable prognosis of most thyroid cancer patients, advanced patients with metastasis and recurrence still have poor prognosis. Therefore, the molecular mechanisms of progression and targeted biomarkers were investigated for developing effective targets for treating thyroid cancer. Eight chip datasets from the gene expression omnibus database were selected and the inSilicoDb and inSilicoMerging R/Bioconductor packages were used to integrate and normalize them across platforms. After merging the eight gene expression omnibus datasets, we obtained one dataset that contained the expression profiles of 319 samples (188 tumor samples plus 131 normal thyroid tissue samples). After screening, we identified 594 significantly differentially expressed genes (277 up-regulated genes plus 317 down-regulated genes) between the tumor and normal tissue samples. The differentially expressed genes exhibited enrichment in multiple signaling pathways, such as p53 signaling. By building a protein–protein interaction network and module analysis, we confirmed seven hub genes, and they were all differentially expressed at all the clinical stages of thyroid cancer. A diagnostic seven-gene signature was established using a logistic regression model with the area under the receiver operating characteristic curve (AUC) of 0.967. Seven robust candidate biomarkers predictive of thyroid cancer were identified, and the obtained seven-gene signature may serve as a useful marker for thyroid cancer diagnosis and prognosis.     

Cluster-Rasch models for microarray gene expression data

Background

We propose two different formulations of the Rasch statistical models to the problem of relating gene expression profiles to the phenotypes. One formulation allows us to investigate whether a cluster of genes with similar expression profiles is related to the observed phenotypes; this model can also be used for future prediction. The other formulation provides an alternative way of identifying genes that are over- or underexpressed from their expression levels in tissue or cell samples of a given tissue or cell type.

Results

We illustrate the methods on available datasets of a classification of acute leukemias and of 60 cancer cell lines. For tumor classification, the results are comparable to those previously obtained. For the cancer cell lines dataset, we found four clusters of genes that are related to drug response for many of the 90 drugs that we considered. In addition, for each type of cell line, we identified genes that are over- or underexpressed relative to other genes.

Conclusions

The cluster-Rasch model provides a probabilistic model for describing gene expression patterns across samples and can be used to relate gene expression profiles to phenotypes.     

Reliable gene signatures for microarray classification: assessment of stability and performance

MOTIVATION: Two important questions for the analysis of gene expression measurements from different sample classes are (1) how to classify samples and (2) how to identify meaningful gene signatures (ranked gene lists) exhibiting the differences between classes and sample subsets. Solutions to both questions have immediate biological and biomedical applications. To achieve optimal classification performance, a suitable combination of classifier and gene selection method needs to be specifically selected for a given dataset. The selected gene signatures can be unstable and the resulting classification accuracy unreliable, particularly when considering different subsets of samples. Both unstable gene signatures and overestimated classification accuracy can impair biological conclusions. METHODS: We address these two issues by repeatedly evaluating the classification performance of all models, i.e. pairwise combinations of various gene selection and classification methods, for random subsets of arrays (sampling). A model score is used to select the most appropriate model for the given dataset. Consensus gene signatures are constructed by extracting those genes frequently selected over many samplings. Sampling additionally permits measurement of the stability of the classification performance for each model, which serves as a measure of model reliability. RESULTS: We analyzed a large gene expression dataset with 78 measurements of four different cartilage sample classes. Classifiers trained on subsets of measurements frequently produce models with highly variable performance. Our approach provides reliable classification performance estimates via sampling. In addition to reliable classification performance, we determined stable consensus signatures (i.e. gene lists) for sample classes. Manual literature screening showed that these genes are highly relevant to our gene expression experiment with osteoarthritic cartilage. We compared our approach to others based on a publicly available dataset on breast cancer. AVAILABILITY: R package at http://www.bio.ifi.lmu.de/~davis/edaprakt     

Peptidylpropyl isomerase B (PPIB): a suitable reference gene for mRNA quantification in peripheral whole blood

Quantitative real-time RT-PCR is a very powerful technique for measuring gene expression at the mRNA level. In order to compare mRNA expression in different experimental or clinical conditions, expression of a target gene has to be normalized to an appropriate internal standard, which is generally a housekeeping gene. In our study, we have tested several housekeeping genes in peripheral whole blood of healthy volunteers and patients suffering from inflammatory diseases. A first analysis of 91 samples illustrated that the mRNA expression of peptidylpropyl isomerase B (PPIB) encoding for cyclophilin B protein, is more stable than beta actin and glyceraldehyde-3-phosphate dehydrogenase, which are both commonly selected as internal standard. Among the three genes tested, beta actin displayed the highest inter-sample variation of expression. The constancy of PPIB mRNA expression was further confirmed by 214 additional samples. In conclusion, we showed that PPIB, in contrast to beta actin and glyceraldehyde-3-phosphate dehydrogenase, is a suitable housekeeping gene in human peripheral blood.