A DNA microarray survey of gene expression in normal human tissues

A Review of Computational Approach for S-system-based Modeling of Gene Regulatory Network

Inference of gene regulatory network (GRN) from time series microarray data remains as a fascinating task for computer science researchers to understand the complex biological process that occurred inside a cell. Among the different popular models to infer GRN, S-system is considered as one of the promising non-linear mathematical tools to model the dynamics of gene expressions, as well as to infer the GRN. S-system is based on biochemical system theory and power law formalism. By observing the value of kinetic parameters of S-system model, it is possible to extract the regulatory relationships among genes. In this review, several existing intelligent methods that were already proposed for inference of S-system-based GRN are explained. It is observed that finding out the most suitable and efficient optimization technique for the accurate inference of all kinds of networks, i.e., in-silico, in-vivo, etc., with less computational complexity is still an open research problem to all. This paper may help the beginners or researchers who want to continue their research in the field of computational biology and bioinformatics.

Crosstalk between SOX Genes and Long Non-Coding RNAs in Glioblastoma

Glioblastoma (GBM) continues to be the most devastating primary brain malignancy. Despite significant advancements in understanding basic GBM biology and enormous efforts in developing new therapeutic approaches, the prognosis for most GBM patients remains poor with a median survival time of 15 months. Recently, the interplay between the SOX (SRY-related HMG-box) genes and lncRNAs (long non-coding RNAs) has become the focus of GBM research. Both classes of molecules have an aberrant expression in GBM and play essential roles in tumor initiation, progression, therapy resistance, and recurrence. In GBM, SOX and lncRNAs crosstalk through numerous functional axes, some of which are part of the complex transcriptional and epigenetic regulatory mechanisms. This review provides a systematic summary of current literature data on the complex interplay between SOX genes and lncRNAs and represents an effort to underscore the effects of SOX/lncRNA crosstalk on the malignant properties of GBM cells. Furthermore, we highlight the significance of this crosstalk in searching for new biomarkers and therapeutic approaches in GBM treatment.

Ionizing Radiation and Estrogen Affecting Growth Factor Genes in an Experimental Breast Cancer Model

Genes associated with growth factors were previously analyzed in a radiation- and estrogen-induced experimental breast cancer model. Such in vitro experimental breast cancer model was developed by exposure of the immortalized human breast epithelial cell line, MCF-10F, to low doses of high linear energy transfer (LET) α particle radiation (150 keV/μm) and subsequent growth in the presence or absence of 17β-estradiol. The MCF-10F cell line was analyzed in different stages of transformation after being irradiated with either a single 60 cGy dose or 60/60 cGy doses of alpha particles. In the present report, the profiling of differentially expressed genes associated with growth factors was analyzed in their relationship with clinical parameters. Thus, the results indicated that Fibroblast growth factor2 gene expression levels were higher in cells transformed by radiation or in the presence of ionizing radiation; whereas the fibroblast growth factor-binding protein 1gene expression was higher in the tumor cell line derived from this model. Such expressions were coincident with higher values in normal than malignant tissues and with estrogen receptor (ER) negative samples for both gene types. The results also showed that transforming growth factor alpha gene expression was higher in the tumor cell line than the tumorigenic A5 and the transformed A3 cell line, whereas the transforming growth factor beta receptor 3 gene expression was higher in A3 and A5 than in Tumor2 cell lines and the untreated controls and the E cell lines. Such gene expression was accompanied by results indicating negative and positive receptors for transforming growth factor alpha and the transforming growth factor beta receptor 3, respectively. Such expressions were low in malignant tissues when compared with benign ones. Furthermore, Fibroblast growth factor2, the fibroblast growth factor-binding protein 1, transforming growth factor alpha, the transforming growth factor beta receptor 3, and the insulin growth factor receptor gene expressions were found to be present in all BRCA patients that are BRCA-Basal, BRCA-LumA, and BRCA-LumB, except in BRCA-Her2 patients. The results also indicated that the insulin growth factor receptor gene expression was higher in the tumor cell line Tumor2 than in Alpha3 cells transformed by ionizing radiation only; then, the insulin growth factor receptor was higher in the A5 than E cell line. The insulin growth factor receptor gene expression was higher in breast cancer than in normal tissues in breast cancer patients. Furthermore, Fibroblast growth factor2, the fibroblast growth factor-binding protein 1, transforming growth factor alpha, the transforming growth factor beta receptor 3, and the insulin growth factor receptor gene expression levels were in stages 3 and 4 of breast cancer patients. It can be concluded that, by using gene technology and molecular information, it is possible to improve therapy and reduce the side effects of therapeutic radiation use. Knowing the different genes involved in breast cancer will make possible the improvement of clinical chemotherapy.

Cell Adhesion Molecules Affected by Ionizing Radiation and Estrogen in an Experimental Breast Cancer Model

Cancer develops in a multi-step process where environmental carcinogenic exposure is a primary etiological component, and where cell-cell communication governs the biological activities of tissues. Identifying the molecular genes that regulate this process is essential to targeting metastatic breast cancer. Ionizing radiation can modify and damage DNA, RNA, and cell membrane components such as lipids and proteins by direct ionization. Comparing differential gene expression can help to determine the effect of radiation and estrogens on cell adhesion. An in vitro experimental breast cancer model was developed by exposure of the immortalized human breast epithelial cell line MCF-10F to low doses of high linear energy transfer α particle radiation and subsequent growth in the presence of 17β-estradiol. The MCF-10F cell line was analyzed in different stages of transformation that showed gradual phenotypic changes including altered morphology, increase in cell proliferation relative to the control, anchorage-independent growth, and invasive capability before becoming tumorigenic in nude mice. This model was used to determine genes associated with cell adhesion and communication such as E-cadherin, the desmocollin 3, the gap junction protein alpha 1, the Integrin alpha 6, the Integrin beta 6, the Keratin 14, Keratin 16, Keratin 17, Keratin 6B, and the laminin beta 3. Results indicated that most genes had greater expression in the tumorigenic cell line Tumor2 derived from the athymic animal than the Alpha3, a non-tumorigenic cell line exposed only to radiation, indicating that altered expression levels of adhesion molecules depended on estrogen. There is a significant need for experimental model systems that facilitate the study of cell plasticity to assess the importance of estrogens in modulating the biology of cancer cells.

Transformation to ischaemia tolerance of frog brain function corresponds to dynamic changes in mRNA co-expression across metabolic pathways

Neural activity is costly and requires continuous ATP from aerobic metabolism. Brainstem motor function of American bullfrogs normally collapses after minutes of ischaemia, but following hibernation, it becomes ischaemia-tolerant, generating output for up to 2 h without oxygen or glucose delivery. Transforming the brainstem to function during ischaemia involves a switch to anaerobic glycolysis and brain glycogen. We hypothesized that improving neural performance during ischaemia involves a transcriptional program for glycogen and glucose metabolism. Here we measured mRNA copy number of genes along the path from glycogen metabolism to lactate production using real-time quantitative PCR. The expression of individual genes did not reflect enhanced glucose metabolism. However, the number of co-expressed gene pairs increased early into hibernation, and by the end, most genes involved in glycogen metabolism, glucose transport and glycolysis exhibited striking linear co-expression. By contrast, co-expression of genes in the Krebs cycle and electron transport chain decreased throughout hibernation. Our results uncover reorganization of the metabolic transcriptional network associated with a shift to ischaemia tolerance in brain function. We conclude that modifying gene co-expression may be a critical step in synchronizing storage and use of glucose to achieve ischaemia tolerance in active neural circuits.

Tear Lipocalin and Lipocalin-Interacting Membrane Receptor

Tear lipocalin is a primate protein that was recognized as a lipocalin from the homology of the primary sequence. The protein is most concentrated in tears and produced by lacrimal glands. Tear lipocalin is also produced in the tongue, pituitary, prostate, and the tracheobronchial tree. Tear lipocalin has been assigned a multitude of functions. The functions of tear lipocalin are inexorably linked to structural characteristics that are often shared by the lipocalin family. These characteristics result in the binding and or transport of a wide range of small hydrophobic molecules. The cavity of tear lipocalin is formed by eight strands (A-H) that are arranged in a β-barrel and are joined by loops between the β-strands. Recently, studies of the solution structure of tear lipocalin have unveiled new structural features such as cation-π interactions, which are extant throughout the lipocalin family. Lipocalin has many unique features that affect ligand specificity. These include a capacious and a flexible cavity with mobile and short overhanging loops. Specific features that confer promiscuity for ligand binding in tear lipocalin will be analyzed. The functions of tear lipocalin include the following: antimicrobial activities, scavenger of toxic and tear disruptive compounds, endonuclease activity, and inhibition of cysteine proteases. In addition, tear lipocalin binds and may modulate lipids in the tears. Such actions support roles as an acceptor for phospholipid transfer protein, heteropolymer formation to alter viscosity, and tear surface interactions. The promiscuous lipid-binding properties of tear lipocalin have created opportunities for its use as a drug carrier. Mutant analogs have been created to bind other molecules such as vascular endothelial growth factor for medicinal use. Tear lipocalin has been touted as a useful biomarker for several diseases including breast cancer, chronic obstructive pulmonary disease, diabetic retinopathy, and keratoconus. The functional possibilities of tear lipocalin dramatically expanded when a putative receptor, lipocalin-interacting membrane receptor was identified. However, opposing studies claim that lipocalin-interacting membrane receptor is not specific for lipocalin. A recent study even suggests a different function for the membrane protein. This controversy will be reviewed in light of gene expression data, which suggest that tear lipocalin has a different tissue distribution than the putative receptor. But the data show lipocalin-interacting membrane receptor is expressed on ocular surface epithelium and that a receptor function here would be rational.

Hepcidin Downregulation Correlates With Disease Aggressiveness And Immune Infiltration in Liver Cancers

Background

Hepcidin is a polypeptide hormone mainly produced by hepatocytes to modulate systemic iron balance. A drastic downregulation of the hepcidin gene was found in liver cancers. However, there is a paucity of information about the clinical significance of hepcidin gene downregulation in liver cancers.

Methods

Hepcidin expression profiles were assessed using multiple public datasets via several bioinformatics platforms. Clinical and pathological information was utilized to stratify patients for comparison. Patient survival outcomes were evaluated using the Kaplan-Meier plotter, a meta-analysis tool. Tumor immune infiltration was analyzed using the single sample gene set enrichment analysis (ssGSEA) approach on the Cancer Genome Atlas (TCGA) dataset. Hepcidin antagonist Fursultiamine was used to treat liver cancer HepG2 and Huh7 cells together with Sorafenib.

Results

Hepcidin gene was predominantly expressed in benign liver tissues but drastically decreased in liver cancer tissues. Hepcidin reduction in liver cancers correlated with risk factors like non-alcoholic fatty liver disease (NAFLD) and liver fibrosis, as well as cancer grade and tumor stage. Hepcidin downregulation was associated with a rapid cancer progression and worse disease-specific survival, especially in patients of the White race without alcohol consumption history. Hepcidin expression in liver cancer tissues positively correlated with the bone morphogenetic protein-6 (BPM6)/interleukin-6 (IL6) cytokines and cytotoxic immune infiltration. Blocking hepcidin action with its antagonist Fursultiamine moderately reduced Sorafenib-induced apoptotic cell death in HepG2 and Huh7 cells.

Conclusion

Hepcidin downregulation in liver cancers correlated with liver cancer risk factors, cancer aggressiveness, cytotoxic immune cell infiltration, and patient survival outcomes. BMP6/IL6 pathway insufficiency is a potential cause of hepcidin downregulation in liver cancers.

Heterogeneity of G protein activation by the calcium-sensing receptor

The calcium-sensing receptor (CaSR) is a G protein-coupled receptor that plays a fundamental role in extracellular calcium (Ca2+e) homeostasis by regulating parathyroid hormone release and urinary calcium excretion. The CaSR has been described to activate all four G protein subfamilies (Gαq/11, Gαi/o, Gα12/13, Gαs), and mutations in the receptor that cause hyper/hypocalcaemia, have been described to bias receptor signalling. However, many of these studies are based on measurements of second messengers or gene transcription that occurs many steps downstream of receptor activation and can represent convergence points of several signalling pathways. Therefore, to assess CaSR-mediated G protein activation directly, we took advantage of a recently described NanoBiT G protein dissociation assay system. Our studies, performed in HEK293 cells stably expressing CaSR, demonstrate that Ca2+e stimulation activates all Gαq/11 family and several Gαi/o family proteins, although Gαz was not activated. CaSR stimulated dissociation of Gα12/13 and Gαs from Gβ-subunits, but this occurred at a slower rate than that of other Gα-subunits. Investigation of cDNA expression of G proteins in three tissues abundantly expressing CaSR, the parathyroids, kidneys and pancreas, showed Gα11, Gαz, Gαi1 and Gα13 genes were highly expressed in parathyroid tissue, indicating CaSR most likely activates Gα11 and Gαi1 in parathyroids. In kidney and pancreas, the majority of G proteins were similarly expressed, suggesting CaSR may activate multiple G proteins in these cells. Thus, these studies validate a single assay system that can be used to robustly assess CaSR variants and biased signalling and could be utilised in the development of new pharmacological compounds targeting CaSR.

The Orphan Nuclear Receptor Gene NR0B2 Is a Favorite Prognosis Factor Modulated by Multiple Cellular Signal Pathways in Human Liver Cancers

BACKGROUND

Liver cancer is a leading cause of cancer death worldwide, and novel prognostic factor is needed for early detection and therapeutic responsiveness monitoring. The orphan nuclear receptor NR0B2 was reported to suppress liver cancer development in a mouse model, and its expression levels were reduced in liver cancer tissues and cell lines due to hypermethylation within its promoter region. However, it is not clear if NR0B2 expression is associated with cancer survival or disease progression and how NR0B2 gene expression is regulated at the molecular level.

METHODS

Multiple cancer databases were utilized to explore NR0B2 gene expression profiles crossing a variety of human cancers, including liver cancers, on several publicly assessable bioinformatics platforms. NR0B2 gene expression with or without kinase inhibitor treatment was analyzed using the qPCR technique, and NR0B2 protein expression was assessed in western blot assays. Two human hepatocellular carcinoma cell lines HepG2 and Huh7, were used in these experiments. NR0B2 gene activation was evaluated using NR0B2 promoter-driven luciferase reporter assays.

RESULTS

NR0B2 gene is predominantly expressed in liver tissue crossing human major organs or tissues, but it is significantly downregulated in liver cancers. NR0B2 expression is mostly downregulated in most common cancers but also upregulated in a few intestinal cancers. NR0B2 gene expression significantly correlated with patient overall survival status in multiple human malignancies, including lung, kidney, breast, urinary bladder, thyroid, colon, and head-neck cancers, as well as liposarcoma and B-cell lymphoma. In liver cancer patients, higher NR0B2 expression is associated with favorite relapse-free and progression-free survival, especially in Asian male patients with viral infection history. In addition, NR0B2 expression negatively correlated with immune infiltration and PIK3CA and PIK3CG gene expression in liver cancer tissues. In HepG2 and Huh7 cells, NR0B2 expression at the transcription level was drastically reduced after MAPK inhibition but was significantly enhanced after PI3K inhibition.

CONCLUSION

NR0B2 gene expression is altered mainly in most human malignancies and significantly reduced in liver cancers. NR0B2 is a prognosis factor for patient survival in liver cancers. MAPK and PI3K oppositely modulate NR0B2 expression, and NR0B2 gene upregulation might serve as a therapeutic responsiveness factor in anti-PI3K therapy for liver cancer.

Over-expression of RALYL suppresses the progression of ovarian clear cell carcinoma through inhibiting MAPK and CDH1 signaling pathways

Background: The molecular mechanism in the progression of ovarian clear cell carcinoma (OCCC) remains unclear. Objective: This study aimed to investigate the potential function of RAYLY in OCCC. Methods: To validate RAYLY expression, immunohistochemistry, quantitative real-time PCR and western blotting were performed in OCCC tissues and the cell lines of OCCC and epithelial ovarian carcinoma (EOC). Subsequently, the biological effects of RALYL were evaluated through colony formation, and cell proliferation, migration and invasion assays. Finally, RNA-sequencing and gene set enrichment analysis (GSEA) were conducted to explore potential mechanism of RALYL in OCCC. Results: In our study, RALYL was significantly down-regulated in a majority of OCCC tissues compared to adjacent non-tumorous tissues, and OCCC cells had a lower expression level of RALYL than that of EOC cells. OCCC patients with high RALYL expression had a better pathological stage and prognosis. In vitro, over-expression of RALYL inhibited cell proliferation, migration and invasion in OCCC. GSEA analysis and western blot indicated an enrichment of MAPK and CDH1 signaling pathways in OCCC cells without RALYL over-expression. Conclusions: RALYL played an important role in the progression of OCCC, and might serve as a potential prognostic biomarker and novel therapeutic target for OCCC.

eBreCaP: extreme learning-based model for breast cancer survival prediction

Breast cancer is the second leading cause of death in the world. Breast cancer research is focused towards its early prediction, diagnosis, and prognosis. Breast cancer can be predicted on omics profiles, clinical tests, and pathological images. The omics profiles comprise of genomic, proteomic, and transcriptomic profiles that are available as high-dimensional datasets. Survival prediction is carried out on omics data to predict early the onset of disease, relapse, reoccurrence of diseases, and biomarker identification. The early prediction of breast cancer is desired for the effective treatment of patients as delay can aggravate the staging of cancer. In this study, extreme learning machine (ELM) based model for breast cancer survival prediction named eBreCaP is proposed. It integrates the genomic (gene expression, copy number alteration, DNA methylation, protein expression) and pathological image datasets; and trains them using an ensemble of ELM with the six best-chosen models suitable to be applied on integrated data. eBreCaP has been evaluated on nine performance parameters, namely sensitivity, specificity, precision, accuracy, Matthews correlation coefficient, area under curve, area under precision-recall, hazard ratio, and concordance Index. eBreCaP has achieved an accuracy of 85% for early breast cancer survival prediction using the ensemble of ELM with gradient boosting.

Selecting Appropriate Reference Genes for Quantitative Real-Time Polymerase Chain Reaction Studies in Isolated and Cultured Ocular Surface Epithelia

The introduction of tissue engineering has allowed scientists to push the boundaries and treat seriously damaged ocular surface epithelia. They have managed to do this through the development of biological substitutes that restore, maintain or improve tissue function. To ensure the generation of a therapeutically safe and effective graft, knowledge on the transcriptional profile of native and cultured ocular surface epithelia is of undeniable value. Gene expression studies are, however, only as reliable as their proper selection of internal reaction controls or reference genes. In this study, we determined the expression stability of a number of reference genes: 18s rRNA, ACTB, ATP5B, CyC1, EIF4A2, GAPDH, RPL13A, SDHA, TOP1, UBC, and YWHAZ in primary isolates as well as in ex vivo cultured ocular surface epithelia explants (day 0 and/or day 14). Expression stability of the reference genes was assessed with both the geNorm and NormFinder software that use a pairwise comparison and a model-based approach, respectively. Our results extend the general recommendation of using multiple reference genes for normalization purposes to our model systems and provide an overview of several references genes that are likely to be stable in similar culture protocols.

S100P is a molecular determinant of E-cadherin function in gastric cancer

Background

E-cadherin has been awarded a key role in the aetiology of both sporadic and hereditary forms of gastric cancer. In this study, we aimed to identify molecular interactors that influence the expression and function of E-cadherin associated to cancer.

Methods

A data mining approach was used to predict stomach-specific candidate genes, uncovering S100P as a key candidate. The role of S100P was evaluated through in vitro functional assays and its expression was studied in a gastric cancer tissue microarray (TMA).

Results

S100P was found to contribute to a cancer pathway dependent on the context of E-cadherin function. In particular, we demonstrated that S100P acts as an E-cadherin positive regulator in a wild-type E-cadherin context, and its inhibition results in decreased E-cadherin expression and function. In contrast, S100P is likely to be a pro-survival factor in gastric cancer cells with loss of functional E-cadherin, contributing to an oncogenic molecular program. Moreover, expression analysis in a gastric cancer TMA revealed that S100P expression impacts negatively among patients bearing Ecad⁻ tumours, despite not being significantly associated with overall survival on its own.

Conclusions

We propose that S100P has a dual role in gastric cancer, acting as an oncogenic factor in the context of E-cadherin loss and as a tumour suppressor in a functional E-cadherin setting. The discovery of antagonist effects of S100P in different E-cadherin contexts will aid in the stratification of gastric cancer patients who may benefit from S100P-targeted therapies.

Graphical abstract

Insights into the Signal Transduction Pathways of Mouse Lung Type II Cells Revealed by Transcription Factor Profiling in the Transcriptome

Alveolar type II cells constitute a small fraction of the total lung cell mass. However, they play an important role in many cellular processes including trans-differentiation into type I cells as well as repair of lung injury in response to toxic chemicals and respiratory pathogens. Transcription factors are the regulatory proteins dynamically modulating DNA structure and gene expression. Transcription factor profiling in microarray datasets revealed that several members of AP1, ATF, NF-kB, and C/EBP families involved in diverse responses were expressed in mouse lung type II cells. A transcriptional factor signature consisting of Cebpa, Srebf1, Stat3, Klf5, and Elf3 was identified in lung type II cells, Sox9+ pluripotent lung stem cells as well as in mouse lung development. Identification of the transcription factor profile in mouse lung type II cells will serve as a useful resource and facilitate the integrated analysis of signal transduction pathways and specific gene targets in a variety of physiological conditions.

Prognostic and Clinicopathological Significance of SERTAD1 in Various Types of Cancer Risk: A Systematic Review and Retrospective Analysis

SERTAD/TRIP-Br genes are considered as a key nuclear transcriptional player in diverse mechanisms of cell including carcinogenesis. The Oncomine™-Online Platform was used for differential expression and biological insights. Kaplan-Meier survival estimated by KM-plotter/cBioPortal/PrognoScan with 95% CI. SERTAD1 was found significantly elevated levels in most of tumor samples. Kaplan-Meier Plotter results distinctly showed the SERTAD1 over-expression significantly reduced median overall-survival (OS) of patients in liver (n = 364/Logrank-test p = 0.0015), ovarian (n = 655/Logrank-test p = 0.00011) and gastric (n = 631/Logrank-test p = 0.1866). Increased level of SERTAD1 has a significantly higher survival rate in the initial time period, but after 100 months slightly reduced OS (n = 26/Logrank-test p = 0.34) and RFS in HER2 positive breast cancer patients. In meta-analysis, cancer patients with higher SERTAD1 mRNA fold resulted worse overall survival than those with lower SERTAD1 levels. Heterogeneity was observed in the fixed effect model analysis DFS [Tau² = 0.0.073, Q (df = 4) = 15.536 (p = 0.004), I² = 74.253], DSS [Tau² = 1.015, Q (df = 2) = 33.214, (p = 0.000), I² = 93.973], RFS [Tau² = 0.492, Q (df = 7) = 71.133 (p = 0.000), I² = 90.159] (Figure 5). OS [Tau² = 0.480, Q (df = 17) = 222.344 (p = 0.000), I² = 92.354]. Lastly, SERTAD1 involved in several signaling cascades through interaction and correlation with many candidate factors as well as miRNAs. This meta-analysis demonstrates a robust evidence of an association between higher or lower SERTAD1, alteration and without alteration of SERTAD1 in cancers in terms of survival and cancer invasiveness.

The Multiple Levels of Mitonuclear Coregulation

Together, the nuclear and mitochondrial genomes encode the oxidative phosphorylation (OXPHOS) complexes that reside in the mitochondrial inner membrane and enable aerobic life. Mitochondria maintain their own genome that is expressed and regulated by factors distinct from their nuclear counterparts. For optimal function, the cell must ensure proper stoichiometric production of OXPHOS subunits by coordinating two physically separated and evolutionarily distinct gene expression systems. Here, we review our current understanding of mitonuclear coregulation primarily at the levels of transcription and translation. Additionally, we discuss other levels of coregulation that may exist but remain largely unexplored, including mRNA modification and stability and posttranslational protein degradation.

Human primitive brain displays negative mitochondrial-nuclear expression correlation of respiratory genes

Oxidative phosphorylation (OXPHOS), a fundamental energy source in all human tissues, requires interactions between mitochondrial (mtDNA)- and nuclear (nDNA)-encoded protein subunits. Although such interactions are fundamental to OXPHOS, bi-genomic coregulation is poorly understood. To address this question, we analyzed ∼8500 RNA-seq experiments from 48 human body sites. Despite well-known variation in mitochondrial activity, quantity, and morphology, we found overall positive mtDNA-nDNA OXPHOS genes' co-expression across human tissues. Nevertheless, negative mtDNA-nDNA gene expression correlation was identified in the hypothalamus, basal ganglia, and amygdala (subcortical brain regions, collectively termed the "primitive" brain). Single-cell RNA-seq analysis of mouse and human brains revealed that this phenomenon is evolutionarily conserved, and both are influenced by brain cell types (involving excitatory/inhibitory neurons and nonneuronal cells) and by their spatial brain location. As the "primitive" brain is highly oxidative, we hypothesized that such negative mtDNA-nDNA co-expression likely controls for the high mtDNA transcript levels, which enforce tight OXPHOS regulation, rather than rewiring toward glycolysis. Accordingly, we found "primitive" brain-specific up-regulation of lactate dehydrogenase B (LDHB), which associates with high OXPHOS activity, at the expense of LDHA, which promotes glycolysis. Analyses of co-expression, DNase-seq, and ChIP-seq experiments revealed candidate RNA-binding proteins and CEBPB as the best regulatory candidates to explain these phenomena. Finally, cross-tissue expression analysis unearthed tissue-dependent splice variants and OXPHOS subunit paralogs and allowed revising the list of canonical OXPHOS transcripts. Taken together, our analysis provides a comprehensive view of mito-nuclear gene co-expression across human tissues and provides overall insights into the bi-genomic regulation of mitochondrial activities.

A novel approach for human whole transcriptome analysis based on absolute gene expression of microarray data

Background

In spite of the emergence of RNA sequencing (RNA-seq), microarrays remain in widespread use for gene expression analysis in the clinic. There are over 767,000 RNA microarrays from human samples in public repositories, which are an invaluable resource for biomedical research and personalized medicine. The absolute gene expression analysis allows the transcriptome profiling of all expressed genes under a specific biological condition without the need of a reference sample. However, the background fluorescence represents a challenge to determine the absolute gene expression in microarrays. Given that the Y chromosome is absent in female subjects, we used it as a new approach for absolute gene expression analysis in which the fluorescence of the Y chromosome genes of female subjects was used as the background fluorescence for all the probes in the microarray. This fluorescence was used to establish an absolute gene expression threshold, allowing the differentiation between expressed and non-expressed genes in microarrays.

Methods

We extracted the RNA from 16 children leukocyte samples (nine males and seven females, ages 6-10 years). An Affymetrix Gene Chip Human Gene 1.0 ST Array was carried out for each sample and the fluorescence of 124 genes of the Y chromosome was used to calculate the absolute gene expression threshold. After that, several expressed and non-expressed genes according to our absolute gene expression threshold were compared against the expression obtained using real-time quantitative polymerase chain reaction (RT-qPCR).

Results

From the 124 genes of the Y chromosome, three genes (DDX3Y, TXLNG2P and EIF1AY) that displayed significant differences between sexes were used to calculate the absolute gene expression threshold. Using this threshold, we selected 13 expressed and non-expressed genes and confirmed their expression level by RT-qPCR. Then, we selected the top 5% most expressed genes and found that several KEGG pathways were significantly enriched. Interestingly, these pathways were related to the typical functions of leukocytes cells, such as antigen processing and presentation and natural killer cell mediated cytotoxicity. We also applied this method to obtain the absolute gene expression threshold in already published microarray data of liver cells, where the top 5% expressed genes showed an enrichment of typical KEGG pathways for liver cells. Our results suggest that the three selected genes of the Y chromosome can be used to calculate an absolute gene expression threshold, allowing a transcriptome profiling of microarray data without the need of an additional reference experiment.

Discussion

Our approach based on the establishment of a threshold for absolute gene expression analysis will allow a new way to analyze thousands of microarrays from public databases. This allows the study of different human diseases without the need of having additional samples for relative expression experiments.

Development of unique cytotoxic chimeric antigen receptors based on human scFv targeting B7H6

As a stress-inducible natural killer (NK) cell ligand, B7H6 plays a role in innate tumor immunosurveillance and is a fairly tumor selective marker expressed on a variety of solid and hematologic cancer cells. Here, we describe the isolation and characterization of a new family of single chain fragment variable (scFv) molecules targeting the human B7H6 ligand. Through directed evolution of a yeast surface displayed non-immune human-derived scFv library, eight candidates comprising a single family of clones differing by up to four amino acid mutations and exhibiting nM avidities for soluble B7H6-Ig were isolated. A representative clone re-formatted as an scFv-CH1-Fc molecule demonstrated specific binding to both B7H6-Ig and native membrane-bound B7H6 on tumor cell lines with a binding avidity comparable to the previously characterized B7H6-targeting antibody, TZ47. Furthermore, these clones recognized an epitope distinct from that of TZ47 and the natural NK cell ligand NKp30, and demonstrated specific activity against B7H6-expressing tumor cells when expressed as a chimeric antigen receptor (CAR) in T cells.

Diverse type 2 diabetes genetic risk factors functionally converge in a phenotype-focused gene network

Type 2 Diabetes (T2D) constitutes a global health burden. Efforts to uncover predisposing genetic variation have been considerable, yet detailed knowledge of the underlying pathogenesis remains poor. Here, we constructed a T2D phenotypic-linkage network (T2D-PLN), by integrating diverse gene functional information that highlight genes, which when disrupted in mice, elicit similar T2D-relevant phenotypes. Sensitising the network to T2D-relevant phenotypes enabled significant functional convergence to be detected between genes implicated in monogenic or syndromic diabetes and genes lying within genomic regions associated with T2D common risk. We extended these analyses to a recent multiethnic T2D case-control exome of 12,940 individuals that found no evidence of T2D risk association for rare frequency variants outside of previously known T2D risk loci. Examining associations involving protein-truncating variants (PTV), most at low population frequencies, the T2D-PLN was able to identify a convergent set of biological pathways that were perturbed within four of five independent T2D case/control ethnic sets of 2000 to 5000 exomes each. These same pathways were found to be over-represented among both known monogenic or syndromic diabetes genes and genes within T2D-associated common risk loci. Our study demonstrates convergent biology amongst variants representing different classes of T2D genetic risk. Although convergence was observed at the pathway level, few of the contributing genes were found in common between different cohorts or variant classes, most notably between the exome variant sets which suggests that future rare variant studies may be better focusing their power onto a single population of recent common ancestry.

Cancer cell line specific co-factors modulate the FOXM1 cistrome

ChIP-seq has been commonly applied to identify genomic occupation of transcription factors (TFs) in a context-specific manner. It is generally assumed that a TF should have similar binding patterns in cells from the same or closely related tissues. Surprisingly, this assumption has not been carefully examined. To this end, we systematically compared the genomic binding of the cell cycle regulator FOXM1 in eight cell lines from seven different human tissues at binding signal, peaks and target genes levels. We found that FOXM1 binding in ER-positive breast cancer cell line MCF-7 are distinct comparing to those in not only other non-breast cell lines, but also MDA-MB-231, ER-negative breast cancer cell line. However, binding sites in MDA-MB-231 and non-breast cell lines were highly consistent. The recruitment of estrogen receptor alpha (ERα) caused the unique FOXM1 binding patterns in MCF-7. Moreover, the activity of FOXM1 in MCF-7 reflects the regulatory functions of ERα, while in MDA-MB-231 and non-breast cell lines, FOXM1 activities regulate cell proliferation. Our results suggest that tissue similarity, in some specific contexts, does not hold precedence over TF-cofactors interactions in determining transcriptional states and that the genomic binding of a TF can be dramatically affected by a particular co-factor under certain conditions.

miR-629-3p may serve as a novel biomarker and potential therapeutic target for lung metastases of triple-negative breast cancer

Background

Different breast cancer subtypes show distinct tropisms for sites of metastasis. Notably, the lung is the most common site for the first distant recurrence in triple-negative breast cancer (TNBC). The identification of novel biomarkers for lung metastasis is of great importance to improving the outcome of TNBC. In this study, we sought to identify a microRNA (miRNA)-based biomarker and therapeutic target for lung metastasis of TNBC.

Methods

A total of 669 patients without de novo stage IV TNBC were recruited for this study. miRNA profiling was conducted in the discovery cohort. Diagnostic accuracy and prognostic values of candidate miRNAs were evaluated in the training and validation cohorts, respectively. The biological functions of candidate miRNAs, as well as potential targets, were further evaluated through bioinformatic analysis as well as by performing in vitro and in vivo assays.

Results

In the discovery set, we found that miR-629-3p was specifically upregulated in both metastatic foci (fold change 144.16, P < 0.0001) and primary tumors (fold change 74.37, P = 0.004) in patients with lung metastases. In the training set, the ROC curve showed that miR-629-3p yielded high diagnostic accuracy in discriminating patients with lung metastasis from patients without recurrence (AUC 0.865, 95% CI 0.800–0.930, P < 0.0001). Although miR-629-3p predicted poor overall survival and disease-free survival in the validation set, it failed to show significance after multivariate analysis. Notably, logistic regression analyses confirmed that miR-629-3p was an independent risk factor for lung metastasis (OR 4.1, 95% CI 2.5–6.6, P < 0.001). Inhibition of miR-629-3p drastically attenuated the viability and migration of TNBC cells, and it markedly suppressed lung metastasis in vivo. Furthermore, we identified the leukemia inhibitory factor receptor (LIFR), a well-known metastatic suppressive gene, to be a direct target of miR-629-3p.

Conclusions

miR-629-3p may serve as a novel biomarker and potential therapeutic target for lung metastases of TNBC mediated via LIFR.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0865-y) contains supplementary material, which is available to authorized users.

A novel FOXA1/ESR1 interacting pathway: A study of Oncomine™ breast cancer microarrays

Forkhead box protein A1 (FOXA1) is essential for the growth and differentiation of breast epithelium, and has a favorable outcome in breast cancer (BC). Elevated FOXA1 expression in BC also facilitates hormone responsiveness in estrogen receptor (ESR)-positive BC. However, the interaction between these two pathways is not fully understood. FOXA1 and GATA binding protein 3 (GATA3) along with ESR1 expression are responsible for maintaining a luminal phenotype, thus suggesting the existence of a strong association between them. The present study utilized the Oncomine™ microarray database to identify FOXA1:ESR1 and FOXA1:ESR1:GATA3 co-expression co-regulated genes. Oncomine™ analysis revealed 115 and 79 overlapping genes clusters in FOXA1:ESR1 and FOXA1:ESR1:GATA3 microarrays, respectively. Five ESR1 direct target genes [trefoil factor 1 (TFF1/PS2), B-cell lymphoma 2 (BCL2), seven in absentia homolog 2 (SIAH2), cellular myeloblastosis viral oncogene homolog (CMYB) and progesterone receptor (PGR)] were detected in the co-expression clusters. To further investigate the role of FOXA1 in ESR1-positive cells, MCF7 cells were transfected with a FOXA1 expression plasmid, and it was observed that the direct target genes of ESR1 (PS2, BCL2, SIAH2 and PGR) were significantly regulated upon transfection. Analysis of one of these target genes, PS2, revealed the presence of two FOXA1 binding sites in the vicinity of the estrogen response element (ERE), which was confirmed by binding assays. Under estrogen stimulation, FOXA1 protein was recruited to the FOXA1 site and could also bind to the ERE site (although in minimal amounts) in the PS2 promoter. Co-transfection of FOXA1/ESR1 expression plasmids demonstrated a significantly regulation of the target genes identified in the FOXA1/ESR1 multi-arrays compared with only FOXA1 transfection, which was suggestive of a synergistic effect of ESR1 and FOXA1 on the target genes. In summary, the present study identified novel FOXA1, ESR1 and GATA3 co-expressed genes that may be involved in breast tumorigenesis.

Allele-specific splicing effects on DKKL1 and ZNF419 transcripts in HeLa cells

Allele-specific splicing is the production of different RNA isoforms from different alleles of a gene. Altered splicing patterns such as exon skipping can have a dramatic effect on the final protein product yet have traditionally proven difficult to predict. We investigated the splicing effects of a set of nine single nucleotide polymorphisms (SNPs) which are predicted to have a direct impact on mRNA splicing, each in a different gene. Predictions were based on SNP location relative to splice junctions and intronic/exonic splicing elements, combined with an analysis of splice isoform expression data from public sources. Of the nine genes tested, six SNPs led to direct impacts on mRNA splicing as determined by the splicing reporter minigene assay and RT-PCR in human HeLa cells, of which four were allele-specific effects. These included previously unreported alternative splicing patterns in the genes ZNF419 and DKKL1. Notably, the SNP in ZNF419, a transcription factor, leads to the deletion of a DNA-binding domain from the protein and is associated with an expression QTL, while the SNP in DKKL1 leads to shortened transcripts predicted to produce a truncated protein. We conclude that the impact of SNP mutations on mRNA splicing, and its biological relevance, can be predicted by integrating SNP position with available data on relative isoform abundance in human cell lines.

A high-throughput fluorimetric microarray with enhanced fluorescence and suppressed "coffee-ring" effects for the detection of calcium ions in blood

A rapid, ultrasensitive, and high-throughput fluorimetric microarray method has been developed using hydrophobic pattern as the microarray substrate and 3-aminopropyltriethoxysilane-coupled carboxylic acid calcium (APS-CCA) as the fluorescent probes for sensing Ca2+ ions in blood. The hydrophobic pattern of the developed Ca2+ analysis microarray could largely suppress the "coffee-ring" effects to facilitate the better distribution density of testing microspots toward the high-throughput detections, and especially prevent the cross-contamination of the multiple samples between adjacent microspots. Moreover, the use of APS matrix could endow the CCA probe the enhanced environmental stability and fluorescence intensity, which is about 2.3-fold higher than that of free CCA. The interactions between APS-CCA and Ca2+ ions were systematically characterized by UV-vis and fluorescence measurements including microscopy imaging. It was demonstrated that the fluorimetric microarray could display the strong capacity of specifically sensing Ca2+ ions with the minimal interferences from blood backgrounds. Such an APS-CCA-based fluorimetric microarray can allow for the analysis of Ca2+ ions down to 0.0050 mM in blood, promising a highly sensitive and selective detection candidate for Ca2+ ions to be applied in the clinical laboratory.

Systematic Analysis of Drug Targets Confirms Expression in Disease-Relevant Tissues

It is commonly assumed that drug targets are expressed in tissues relevant to their indicated diseases, even under normal conditions. While multiple anecdotal cases support this hypothesis, a comprehensive study has not been performed to verify it. We conducted a systematic analysis to assess gene and protein expression for all targets of marketed and phase III drugs across a diverse collection of normal human tissues. For 87% of gene-disease pairs, the target is expressed in a disease-affected tissue under healthy conditions. This result validates the importance of confirming expression of a novel drug target in an appropriate tissue for each disease indication and strengthens previous findings showing that targets of efficacious drugs should be expressed in relevant tissues under normal conditions. Further characterization of the remaining 13% of gene-disease pairs revealed that most genes are expressed in a different tissue linked to another disease. Our analysis demonstrates the value of extensive tissue specific expression resources.both in terms of tissue and cell diversity as well as techniques used to measure gene expression.

A 3-methylcrotonyl-CoA carboxylase deficient human skin fibroblast transcriptome reveals underlying mitochondrial dysfunction and oxidative stress

Isolated 3-methylcrotonyl-CoA carboxylase (MCC) deficiency is an autosomal recessive inherited metabolic disease of leucine catabolism with a highly variable phenotype. Apart from extensive mutation analyses of the MCCC1 and MCCC2 genes encoding 3-methylcrotonyl-CoA carboxylase (EC 6.4.1.4), molecular data on MCC deficiency gene expression studies in human tissues is lacking. For IEMs, unbiased '-omics' approaches are starting to reveal the secondary cellular responses to defects in biochemical pathways. Here we present the first whole genome expression profile of immortalized cultured skin fibroblast cells of two clinically affected MCC deficient patients and two healthy individuals generated using Affymetrix(®)HuExST1.0 arrays. There were 16191 significantly differentially expressed transcript IDs of which 3591 were well annotated and present in the predefined knowledge database of Ingenuity Pathway Analysis software used for downstream functional analyses. The most noticeable feature of this MCCA deficient skin fibroblast transcriptome was the typical genetic hallmark of mitochondrial dysfunction, decreased antioxidant response and disruption of energy homeostasis, which was confirmed by mitochondrial functional analyses. The MCC deficient transcriptome seems to predict oxidative stress that could alter the complex secondary cellular response that involve genes of the glycolysis, the TCA cycle, OXPHOS, gluconeogenesis, β-oxidation and the branched-chain fatty acid metabolism. An important emerging insight from this human MCCA transcriptome in combination with previous reports is that chronic exposure to the primary and secondary metabolites of MCC deficiency and the resulting oxidative stress might impact adversely on the quality of life and energy levels, irrespective of whether MCC deficient individuals are clinically affected or asymptomatic.

Supplementation of essential fatty acids to Holstein calves during late uterine life and first month of life alters hepatic fatty acid profile and gene expression

Linoleic acid is an essential dietary fatty acid (FA). However, how the supplementation of linoleic acid during uterine and early life may modify the FA profile and transcriptome regulation of the liver, and performance of preweaned dairy calves is unknown. Our objective was to evaluate the effect of supplementation of essential FA to Holstein calves during late uterine and early life on their hepatic FA profile and global gene expression at 30 d of age. During the last 8 wk of pregnancy, Holstein cattle (n=96) were fed either no fat supplement (control), a saturated FA supplement enriched with C18:0, or an unsaturated FA supplement enriched with linoleic acid. Male calves (n=40) born from these dams were fed a milk replacer (MR) with either low (LLA) or high linoleic acid (HLA) concentration as the sole feedstuff during the first 30 d. Liver biopsy was performed at 30 d of age, and microarray analysis was performed on 18 liver samples. Total concentration of FA in liver were greater in calves fed LLA compared with those fed HLA MR (8.2 vs. 7.1%), but plasma concentrations of total FA did not differ due to MR diets. The FA profiles of plasma and liver of calves were affected differently by the prepartum diets. Specifically, the FA profile in liver was affected moderately by the feeding of fat prepartum, but the profiles did not differ due to the type of FA fed prepartum. The type of MR fed during the first 30 d of life had major effects on both plasma and liver FA profiles, resembling the type of fat fed. Plasma and liver of calves fed LLA MR had greater percentage of medium-chain FA (C12:0 and C14:0), whereas plasma and liver from calves fed HLA MR had greater percentages of linoleic and α-linolenic acids. Dams fed fat or a specific type of FA modified the expression of some genes in liver of calves, particularly those genes involved in biological functions and pathways related to upregulation of lipid metabolism and downregulation of inflammatory responses. Feeding HLA instead of LLA MR modified the expression of hepatic genes, including genes predicted to decrease infections and to increase lipid utilization and protein synthesis. Research evaluating the effect of FA supplementation during uterine and neonatal life on the future productivity of the neonate is warranted.

T-LAK Cell-Originated Protein Kinase (TOPK) as a Prognostic Factor and a Potential Therapeutic Target in Ovarian Cancer

BACKGROUND

We aimed to clarify the clinical significance of TOPK (T-lymphokine-activated killer cell-originated protein kinase) expression in ovarian cancer and evaluate the possible effect of TOPK inhibitors, OTS514 and OTS964, on ovarian cancer cells.

METHODS

TOPK expression was examined by immunohistochemistry using 163 samples with epithelial ovarian cancer (EOC). TOPK protein level and FOXM1 transcriptional level in ovarian cancer cell lines were examined by Western blot and RT-PCR, respectively. Half-maximum inhibitory concentration (IC₅₀) values against TOPK inhibitors were examined by the MTT assay. Using the peritoneal dissemination model of ES-2 ovarian cancer cells, we examined the in vivo efficacy of OTS514. In addition, the cytotoxic effect of OTS514 and OTS964 on 31 patient-derived primary ovarian cancer cells was examined.

RESULTS

TOPK was expressed very highly in 84 (52%) of 163 EOC tissues, and high TOPK expression was significantly associated with poor progression-free survival and overall survival in early-stage cases of EOC (P = 0.008 and 0.006, respectively). Both OTS514 and OTS964 showed significant growth-inhibitory effect on ovarian cancer cell lines with IC₅₀ values of 3.0 to 46 nmol/L and 14 to 110 nmol/L, respectively. TOPK protein and transcriptional levels of FOXM1 were reduced by TOPK inhibitor treatment. Oral administration of OTS514 significantly elongated overall survival in the ES-2 abdominal dissemination xenograft model, compared with vehicle control (P < 0.001). Two drugs showed strong growth-inhibitory effect on primary ovarian cancer cells regardless of tumor sites or histological subtypes.

CONCLUSIONS

Our results demonstrated the clinical significance of high TOPK expression and potential of TOPK inhibitors to treat ovarian cancer. Clin Cancer Res; 22(24); 6110-7. ©2016 AACR.

Ly6E/K Signaling to TGFβ Promotes Breast Cancer Progression, Immune Escape, and Drug Resistance

Stem cell antigen Sca-1 is implicated in murine cancer stem cell biology and breast cancer models, but the role of its human homologs Ly6K and Ly6E in breast cancer are not established. Here we report increased expression of Ly6K/E in human breast cancer specimens correlates with poor overall survival, with an additional specific role for Ly6E in poor therapeutic outcomes. Increased expression of Ly6K/E also correlated with increased expression of the immune checkpoint molecules PDL1 and CTLA4, increased tumor-infiltrating T regulatory cells, and decreased natural killer (NK) cell activation. Mechanistically, Ly6K/E was required for TGFβ signaling and proliferation in breast cancer cells, where they contributed to phosphorylation of Smad1/5 and Smad2/3. Furthermore, Ly6K/E promoted cytokine-induced PDL1 expression and activation and binding of NK cells to cancer cells. Finally, we found that Ly6K/E promoted drug resistance and facilitated immune escape in this setting. Overall, our results establish a pivotal role for a Ly6K/E signaling axis involving TGFβ in breast cancer pathophysiology and drug response, and highlight this signaling axis as a compelling realm for therapeutic invention. Cancer Res; 76(11); 3376-86. ©2016 AACR.

Detection bias in microarray and sequencing transcriptomic analysis identified by housekeeping genes

This work includes the original data used to discover the gene ontology bias in transcriptomic analysis conducted by microarray and high throughput sequencing (Zhang et al., 2015) [1]. In the analysis, housekeeping genes were used to examine the differential detection ability by microarray and sequencing because these genes are probably the most reliably detected. The genes included here were compiled from 15 human housekeeping gene studies. The provided tables here comprise of detailed chromosomal location, detection breadth, normalized expression level, exon count, total exon length, and total intron length of each concerned gene and their related transcripts. We hope this information can help researchers better understand the differences in gene ontology-bias we discussed (Zhang et al., 2015) [1] and can encourage further improvement on these two technology platforms.

A common molecular signature in ASD gene expression: following Root 66 to autism

Several gene expression experiments on autism spectrum disorders have been conducted using both blood and brain tissue. Individually, these studies have advanced our understanding of the molecular systems involved in the molecular pathology of autism and have formed the bases of ongoing work to build autism biomarkers. In this study, we conducted an integrated systems biology analysis of 9 independent gene expression experiments covering 657 autism, 9 mental retardation and developmental delay and 566 control samples to determine if a common signature exists and to test whether regulatory patterns in the brain relevant to autism can also be detected in blood. We constructed a matrix of differentially expressed genes from these experiments and used a Jaccard coefficient to create a gene-based phylogeny, validated by bootstrap. As expected, experiments and tissue types clustered together with high statistical confidence. However, we discovered a statistically significant subgrouping of 3 blood and 2 brain data sets from 3 different experiments rooted by a highly correlated regulatory pattern of 66 genes. This Root 66 appeared to be non-random and of potential etiologic relevance to autism, given their enriched roles in neurological processes key for normal brain growth and function, learning and memory, neurodegeneration, social behavior and cognition. Our results suggest that there is a detectable autism signature in the blood that may be a molecular echo of autism-related dysregulation in the brain.

Transcriptome analysis of complex I-deficient patients reveals distinct expression programs for subunits and assembly factors of the oxidative phosphorylation system

Background

Transcriptional control of mitochondrial metabolism is essential for cellular function. A better understanding of this process will aid the elucidation of mitochondrial disorders, in particular of the many genetically unsolved cases of oxidative phosphorylation (OXPHOS) deficiency. Yet, to date only few studies have investigated nuclear gene regulation in the context of OXPHOS deficiency. In this study we performed RNA sequencing of two control and two complex I-deficient patient cell lines cultured in the presence of compounds that perturb mitochondrial metabolism: chloramphenicol, AICAR, or resveratrol. We combined this with a comprehensive analysis of mitochondrial and nuclear gene expression patterns, co-expression calculations and transcription factor binding sites.

Results

Our analyses show that subsets of mitochondrial OXPHOS genes respond opposingly to chloramphenicol and AICAR, whereas the response of nuclear OXPHOS genes is less consistent between cell lines and treatments. Across all samples nuclear OXPHOS genes have a significantly higher co-expression with each other than with other genes, including those encoding mitochondrial proteins. We found no evidence for complex-specific mRNA expression regulation: subunits of different OXPHOS complexes are similarly (co-)expressed and regulated by a common set of transcription factors. However, we did observe significant differences between the expression of nuclear genes for OXPHOS subunits versus assembly factors, suggesting divergent transcription programs. Furthermore, complex I co-expression calculations identified 684 genes with a likely role in OXPHOS biogenesis and function. Analysis of evolutionarily conserved transcription factor binding sites in the promoters of these genes revealed almost all known OXPHOS regulators (including GABP, NRF1/2, SP1, YY1, E-box factors) and a set of novel candidates (ELK1, KLF7, SP4, EHF, ZNF143, and TEL2).

Conclusions

OXPHOS genes share an expression program distinct from other genes encoding mitochondrial proteins, indicative of targeted nuclear regulation of a mitochondrial sub-process. Within the subset of OXPHOS genes we established a difference in expression between mitochondrial and nuclear genes, and between nuclear genes encoding subunits and assembly factors. Most transcription regulators of genes that co-express with complex I are well-established factors for OXPHOS biogenesis. For the remaining six factors we here suggest for the first time a link with transcription regulation in OXPHOS deficiency.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1883-8) contains supplementary material, which is available to authorized users.

Identification of putative drug targets for human sperm-egg interaction defect using protein network approach

Background

Sperm-egg interaction defect is a significant cause of in-vitro fertilization failure for infertile cases. Numerous molecular interactions in the form of protein-protein interactions mediate the sperm-egg membrane interaction process. Recent studies have demonstrated that in addition to experimental techniques, computational methods, namely protein interaction network approach, can address protein-protein interactions between human sperm and egg. Up to now, no drugs have been detected to treat sperm-egg interaction disorder, and the initial step in drug discovery research is finding out essential proteins or drug targets for a biological process. The main purpose of this study is to identify putative drug targets for human sperm-egg interaction deficiency and consider if the detected essential proteins are targets for any known drugs using protein-protein interaction network and ingenuity pathway analysis.

Results

We have created human sperm-egg protein interaction networks with high confidence, including 106 nodes and 415 interactions. Through topological analysis of the network with calculation of some metrics, such as connectivity and betweenness centrality, we have identified 13 essential proteins as putative drug targets. The potential drug targets are from integrins, fibronectins, epidermal growth factor receptors, collagens and tetraspanins protein families. We evaluated these targets by ingenuity pathway analysis, and the known drugs for the targets have been detected, and the possible effective role of the drugs on sperm-egg interaction defect has been considered. These results showed that the drugs ocriplasmin (Jetrea©), gefitinib (Iressa©), erlotinib hydrochloride (Tarceva©), clingitide, cetuximab (Erbitux©) and panitumumab (Vectibix©) are possible candidates for efficacy testing for the treatment of sperm-egg interaction deficiency. Further experimental validation can be carried out to confirm these results.

Conclusion

We have identified the first potential list of drug targets for human sperm-egg interaction defect using the protein interaction network approach. The essential proteins or potential drug targets were found using topological analysis of the protein network. These putative targets are promising for further experimental validation. These study results, if validated, may develop drug discovery techniques for sperm-egg interaction defect and also improve assisted reproductive technologies to avoid in-vitro fertilization failure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0186-7) contains supplementary material, which is available to authorized users.

Comprehensive comparison of large-scale tissue expression datasets

For tissues to carry out their functions, they rely on the right proteins to be present. Several high-throughput technologies have been used to map out which proteins are expressed in which tissues; however, the data have not previously been systematically compared and integrated. We present a comprehensive evaluation of tissue expression data from a variety of experimental techniques and show that these agree surprisingly well with each other and with results from literature curation and text mining. We further found that most datasets support the assumed but not demonstrated distinction between tissue-specific and ubiquitous expression. By developing comparable confidence scores for all types of evidence, we show that it is possible to improve both quality and coverage by combining the datasets. To facilitate use and visualization of our work, we have developed the TISSUES resource (http://tissues.jensenlab.org), which makes all the scored and integrated data available through a single user-friendly web interface.

Effects of Slit3 silencing on the invasive ability of lung carcinoma A549 cells

Slit proteins function as chemorepellents in axon guidance and neuronal migration by binding to cognate Robo receptors. The Slit/Robo signaling pathway is also involved in the regulation of tumor cell metastasis. However, whether the Slit/Robo signaling pathway exerts prometastatic or antimetastasis functions remains controversial. To date, most of the research on Slit/Robo has focused on Slit2, and the effects of Slit3 on metastasis remain largely unknown. Based on the Oncomine database, overall expression of Slit3 is low in tumor tissues compared to its level in normal tissues. The underlying mechanism for slit3 silencing in tumor tissues is likely related to hypermethylation of the slit3 promoter. However, lung carcinomas appear to be an exception. Several studies have reported that the frequency of Slit3 methylation in lung cancers is far lower than the frequency of Slit2. In the present study, high Slit3 expression at the mRNA level, yet not at the protein level, was detected in lung adenocarcinoma A549 cells. The function of Slit3 in tumor migration and invasion was examined by silencing of Slit3 expression in A549 cells. Silencing of Slit3 promoted proliferation, migration and invasion of A549 cells and induced epithelial-mesenchymal transition by downregulation of E-cadherin and upregulation of vimentin. The inhibitory effects of Slit3 on tumor migration and invasion are likely related to matrix metalloproteinases (MMPs). Silencing of Slit3 in the A549 cells enhanced MMP2 and MMP9 expression. These results indicate that Slit3 is a potential tumor suppressor in lung adenocarcinoma.

Do housekeeping genes exist?

The searching of human housekeeping (HK) genes has been a long quest since the emergence of transcriptomics, and is instrumental for us to understand the structure of genome and the fundamentals of biological processes. The resolved genes are frequently used in evolution studies and as normalization standards in quantitative gene-expression analysis. Within the past 20 years, more than a dozen HK-gene studies have been conducted, yet none of them sampled human tissues completely. We believe an integration of these results will help remove false positive genes owing to the inadequate sampling. Surprisingly, we only find one common gene across 15 examined HK-gene datasets comprising 187 different tissue and cell types. Our subsequent analyses suggest that it might not be appropriate to rigidly define HK genes as expressed in all tissue types that have diverse developmental, physiological, and pathological states. It might be beneficial to use more robustly identified HK functions for filtering criteria, in which the representing genes can be a subset of genome. These genes are not necessarily the same, and perhaps need not to be the same, everywhere in our body.

Submicroscopic deletions at 13q32.1 cause congenital microcoria

Congenital microcoria (MCOR) is a rare autosomal-dominant disorder characterized by inability of the iris to dilate owing to absence of dilator pupillae muscle. So far, a dozen MCOR-affected families have been reported worldwide. By using whole-genome oligonucleotide array CGH, we have identified deletions at 13q32.1 segregating with MCOR in six families originating from France, Japan, and Mexico. Breakpoint sequence analyses showed nonrecurrent deletions in 5/6 families. The deletions varied from 35 kbp to 80 kbp in size, but invariably encompassed or interrupted only two genes: TGDS encoding the TDP-glucose 4,6-dehydratase and GPR180 encoding the G protein-coupled receptor 180, also known as intimal thickness-related receptor (ITR). Unlike TGDS which has no known function in muscle cells, GPR180 is involved in the regulation of smooth muscle cell growth. The identification of a null GPR180 mutation segregating over two generations with iridocorneal angle dysgenesis, which can be regarded as a MCOR endophenotype, is consistent with the view that deletions of this gene, with or without the loss of elements regulating the expression of neighboring genes, are the cause of MCOR.

Regulated genes in mesenchymal stem cells and gastric cancer

AIM: To investigate the genes regulated in mesenchymal stem cells (MSCs) and diffuse-type gastric cancer (GC), gene expression was analyzed.

METHODS: Gene expression of MSCs and diffuse-type GC cells were analyzed by microarray. Genes related to stem cells, cancer and the epithelial-mesenchymal transition (EMT) were extracted from human gene lists using Gene Ontology and reference information. Gene panels were generated, and messenger RNA gene expression in MSCs and diffuse-type GC cells was analyzed. Cluster analysis was performed using the NCSS software.

RESULTS: The gene expression of regulator of G-protein signaling 1 (RGS1) was up-regulated in diffuse-type GC cells compared with MSCs. A panel of stem-cell related genes and genes involved in cancer or the EMT were examined. Stem-cell related genes, such as growth arrest-specific 6, musashi RNA-binding protein 2 and hairy and enhancer of split 1 (Drosophila), NOTCH family genes and Notch ligands, such as delta-like 1 (Drosophila) and Jagged 2, were regulated.

CONCLUSION: Expression of RGS1 is up-regulated, and genes related to stem cells and NOTCH signaling are altered in diffuse-type GC compared with MSCs.

The human gastrointestinal tract-specific transcriptome and proteome as defined by RNA sequencing and antibody-based profiling

BACKGROUND

The gastrointestinal tract (GIT) is subdivided into different anatomical organs with many shared functions and characteristics, but also distinct differences. We have combined a genome-wide transcriptomics analysis with immunohistochemistry-based protein profiling to describe the gene and protein expression patterns that define the human GIT.

METHODS

RNA sequencing data derived from stomach, duodenum, jejunum/ileum and colon specimens were compared to gene expression levels in 23 other normal human tissues analysed with the same method. Protein profiling based on immunohistochemistry and tissue microarrays was used to sub-localize the corresponding proteins with GIT-specific expression into sub-cellular compartments and cell types.

RESULTS

Approximately 75% of all human protein-coding genes were expressed in at least one of the GIT tissues. Only 51 genes showed enriched expression in either one of the GIT tissues and an additional 83 genes were enriched in two or more GIT tissues. The list of GIT-enriched genes with validated protein expression patterns included various well-known but also previously uncharacterised or poorly studied genes. For instance, the colon-enriched expression of NXPE family member 1 (NXPE1) was established, while NLR family, pyrin domain-containing 6 (NLRP6) expression was primarily found in the human small intestine.

CONCLUSIONS

We have applied a genome-wide analysis based on transcriptomics and antibody-based protein profiling to identify genes that are expressed in a specific manner within the human GIT. These genes and proteins constitute important starting points for an improved understanding of the normal function and the different states of disease associated with the GIT.

Key regulators in prostate cancer identified by co-expression module analysis

Background

Prostate cancer (PrCa) is the most commonly diagnosed cancer in men in the world. Despite the fact that a large number of its genes have been investigated, its etiology remains poorly understood. Furthermore, most PrCa candidate genes have not been rigorously replicated, and the methods by which they biologically function in PrCa remain largely unknown.

Results

Aiming to identify key players in the complex prostate cancer system, we reconstructed PrCa co-expressed modules within functional gene sets defined by the Gene Ontology (GO) annotation (biological process, GO_BP). We primarily identified 118 GO_BP terms that were well-preserved between two independent gene expression datasets and a consequent 55 conserved co-expression modules within them. Five modules were then found to be significantly enriched with PrCa candidate genes collected from expression Quantitative Trait Loci (eQTL), somatic copy number alteration (SCNA), somatic mutation data, or prognostic analyses. Specifically, two transcription factors (TFs) (NFAT and SP1) and three microRNAs (hsa-miR-19a, hsa-miR-15a, and hsa-miR-200b) regulating these five candidate modules were found to be critical to the development of PrCa.

Conclusions

Collectively, our results indicated that genes with similar functions may play important roles in disease through co-expression, and modules with different functions could be regulated by similar genetic components, such as TFs and microRNAs, in a synergistic manner.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1015) contains supplementary material, which is available to authorized users.

Comparison of Metabolic Network between Muscle and Intramuscular Adipose Tissues in Hanwoo Beef Cattle Using a Systems Biology Approach

The interrelationship between muscle and adipose tissues plays a major role in determining the quality of carcass traits. The objective of this study was to compare metabolic differences between muscle and intramuscular adipose (IMA) tissues in the longissimus dorsi (LD) of Hanwoo (Bos taurus coreanae) using the RNA-seq technology and a systems biology approach. The LD sections between the 6th and 7th ribs were removed from nine (each of three cows, steers, and bulls) Hanwoo beef cattle (carcass weight of 430.2 ± 40.66 kg) immediately after slaughter. The total mRNA from muscle, IMA, and subcutaneous adipose and omental adipose tissues were isolated and sequenced. The reads that passed quality control were mapped onto the bovine reference genome (build bosTau6), and differentially expressed genes across tissues were identified. The KEGG pathway enrichment tests revealed the opposite direction of metabolic regulation between muscle and IMA. Metabolic gene network analysis clearly indicated that oxidative metabolism was upregulated in muscle and downregulated in IMA. Interestingly, pathways for regulating cell adhesion, structure, and integrity and chemokine signaling pathway were upregulated in IMA and downregulated in muscle. It is thus inferred that IMA may play an important role in the regulation of development and structure of the LD tissues and muscle/adipose communication.

Benzo pyrene-induced DNA adducts and gene expression profiles in target and non-target organs for carcinogenesis in mice

BACKGROUND

Gene expression changes induced by carcinogens may identify differences in molecular function between target and non-target organs. Target organs for benzo[a]pyrene (BaP) carcinogenicity in mice (lung, spleen and forestomach) and three non-target organs (liver, colon and glandular stomach) were investigated for DNA adducts by 32P-postlabelling, for gene expression changes by cDNA microarray and for miRNA expression changes by miRNA microarray after exposure of animals to BaP.

RESULTS

BaP-DNA adduct formation occurred in all six organs at levels that did not distinguish between target and non-target. cDNA microarray analysis showed a variety of genes modulated significantly by BaP in the six organs and the overall gene expression patterns were tissue specific. Gene ontology analysis also revealed that BaP-induced bioactivities were tissue specific; eight genes (Tubb5, Fos, Cdh1, Cyp1a1, Apc, Myc, Ctnnb1 and Cav) showed significant expression difference between three target and three non-target organs. Additionally, several gene expression changes, such as in Trp53 activation and Stat3 activity suggested some similarities in molecular mechanisms in two target organs (lung and spleen), which were not found in the other four organs. Changes in miRNA expression were generally tissue specific, involving, in total, 21/54 miRNAs significantly up- or down-regulated.

CONCLUSIONS

Altogether, these findings showed that DNA adduct levels and early gene expression changes did not fully distinguish target from non-target organs. However, mechanisms related to early changes in p53, Stat3 and Wnt/β-catenin pathways may play roles in defining BaP organotropism.

Unbiased functional clustering of gene variants with a phenotypic-linkage network

Groupwise functional analysis of gene variants is becoming standard in next-generation sequencing studies. As the function of many genes is unknown and their classification to pathways is scant, functional associations between genes are often inferred from large-scale omics data. Such data types—including protein–protein interactions and gene co-expression networks—are used to examine the interrelations of the implicated genes. Statistical significance is assessed by comparing the interconnectedness of the mutated genes with that of random gene sets. However, interconnectedness can be affected by confounding bias, potentially resulting in false positive findings. We show that genes implicated through de novo sequence variants are biased in their coding-sequence length and longer genes tend to cluster together, which leads to exaggerated p-values in functional studies; we present here an integrative method that addresses these bias. To discern molecular pathways relevant to complex disease, we have inferred functional associations between human genes from diverse data types and assessed them with a novel phenotype-based method. Examining the functional association between de novo gene variants, we control for the heretofore unexplored confounding bias in coding-sequence length. We test different data types and networks and find that the disease-associated genes cluster more significantly in an integrated phenotypic-linkage network than in other gene networks. We present a tool of superior power to identify functional associations among genes mutated in the same disease even after accounting for significant sequencing study bias and demonstrate the suitability of this method to functionally cluster variant genes underlying polygenic disorders.

Plenty of gene variants have been associated with a disease, yet most of the heritability, along with the molecular basis, of common diseases remains unexplained. However, it is widely thought that the products of genes whose mutations are implicated in the same disease function together in the same biological pathways and it is the disruption of these pathways that underlies the disease. Such pathways are not well defined and their identification could help elucidate disease mechanisms. Consequently, groupwise functional analyses of gene variants to identify common disease-relevant pathways are becoming standard in next-generation sequencing studies, but we find that these analyses are confounded by coding-sequence length bias. We control for these bias and describe a phenotype-based approach which outperforms other methods in discerning functional associations among the disease-associated genes. We also demonstrate the suitability of this method to functionally dissect the gene variants underlying a complex disorder, the identified functional clusters offering insight into disease mechanisms.

Current advances in systems and integrative biology

Systems biology has gained a tremendous amount of interest in the last few years. This is partly due to the realization that traditional approaches focusing only on a few molecules at a time cannot describe the impact of aberrant or modulated molecular environments across a whole system. Furthermore, a hypothesis-driven study aims to prove or disprove its postulations, whereas a hypothesis-free systems approach can yield an unbiased and novel testable hypothesis as an end-result. This latter approach foregoes assumptions which predict how a biological system should react to an altered microenvironment within a cellular context, across a tissue or impacting on distant organs. Additionally, re-use of existing data by systematic data mining and re-stratification, one of the cornerstones of integrative systems biology, is also gaining attention. While tremendous efforts using a systems methodology have already yielded excellent results, it is apparent that a lack of suitable analytic tools and purpose-built databases poses a major bottleneck in applying a systematic workflow. This review addresses the current approaches used in systems analysis and obstacles often encountered in large-scale data analysis and integration which tend to go unnoticed, but have a direct impact on the final outcome of a systems approach. Its wide applicability, ranging from basic research, disease descriptors, pharmacological studies, to personalized medicine, makes this emerging approach well suited to address biological and medical questions where conventional methods are not ideal.

Comparative transcriptomic analysis by RNA-seq to discern differential expression of genes in liver and muscle tissues of adult Berkshire and Jeju Native Pig

RNA-seq is being rapidly adopted for the profiling of the transcriptomes in different areas of biology, especially in the studies related to gene regulation. The discovery of differentially expressed genes (DEGs) between adult animals of Jeju Native Pig (JNP) and Berkshire breeds of Sus scrofa, is of particular interest for the current study. For the better understanding of the gene expression profiles of the liver and longissimus dorsi muscle, DEGs were identified via RNA-seq. Sequence reads were obtained from Illumina HiSeq2000 and mapped to the pig reference genome (Sscrofa10.2) using Tophat2. We identified 169 and 39 DEGs in the liver and muscle of JNP respectively, by comparison with Berkshire breed. Out of all identified genes, 41 genes in the liver and 9 genes in the muscle have given significant expression. Gene ontology (GO) terms of developmental process and KEGG pathway analysis showed that metabolic, immune response and protein binding were commonly enriched pathways in the two tissues. Further the heat map analysis by ArrayStar has shown the different levels of expression in JNP with respect to the Berkshire breed. The validation through real time PCR and western blotting also confirmed the differential expression of genes in both breeds. Genes pertaining to metabolic process and inflammatory and immune system are more enriched in Berkshire breed. This comparative transcriptome analysis of two tissues suggests a subset of novel marker genes which expressed differently between the JNP and Berkshire.