The composition of the intestinal microbiota influences the development of inflammatory disorders. However, associating inflammatory diseases with specific microbial members of the microbiota is challenging, because clinically detectable inflammation and its treatment can alter the microbiota's composition. Immunologic checkpoint blockade with ipilimumab, a monoclonal antibody that blocks cytotoxic T-lymphocyte-associated antigen-4 (CTLA-4) signalling, is associated with new-onset, immune-mediated colitis. Here we conduct a prospective study of patients with metastatic melanoma undergoing ipilimumab treatment and correlate the pre-inflammation faecal microbiota and microbiome composition with subsequent colitis development. We demonstrate that increased representation of bacteria belonging to the Bacteroidetes phylum is correlated with resistance to the development of checkpoint-blockade-induced colitis. Furthermore, a paucity of genetic pathways involved in polyamine transport and B vitamin biosynthesis is associated with an increased risk of colitis. Identification of these biomarkers may enable interventions to reduce the risk of inflammatory complications following cancer immunotherapy.
A subset of cancer patients treated with immune checkpoint blockade develops colitis. Here the authors show that lower abundance of Bacteroidetes and vitamin B biosynthetic modules in fecal samples of melanoma patients can predict their susceptibility to colitis following anti-CTLA-4 treatment.
There are conflicting reports on the impact of soy on breast carcinogenesis. This study examines the effects of soy supplementation on breast cancer-related genes and pathways.
Women (n = 140) with early-stage breast cancer were randomly assigned to soy protein supplementation (n = 70) or placebo (n = 70) for 7 to 30 days, from diagnosis until surgery. Adherence was determined by plasma isoflavones: genistein and daidzein. Gene expression changes were evaluated by NanoString in pre- and posttreatment tumor tissue. Genome-wide expression analysis was performed on posttreatment tissue. Proliferation (Ki67) and apoptosis (Cas3) were assessed by immunohistochemistry.
Plasma isoflavones rose in the soy group (two-sided Wilcoxon rank-sum test, P < .001) and did not change in the placebo group. In paired analysis of pre- and posttreatment samples, 21 genes (out of 202) showed altered expression (two-sided Student’s t-test, P < .05). Several genes including FANCC and UGT2A1 revealed different magnitude and direction of expression changes between the two groups (two-sided Student’s t-test, P < .05). A high-genistein signature consisting of 126 differentially expressed genes was identified from microarray analysis of tumors. This signature was characterized by overexpression (>2-fold) of cell cycle transcripts, including those that promote cell proliferation, such as FGFR2, E2F5, BUB1, CCNB2, MYBL2, CDK1, and CDC20 (P < .01). Soy intake did not result in statistically significant changes in Ki67 or Cas3.
Gene expression associated with soy intake and high plasma genistein defines a signature characterized by overexpression of FGFR2 and genes that drive cell cycle and proliferation pathways. These findings raise the concerns that in a subset of women soy could adversely affect gene expression in breast cancer.
Epigenetic regulation of key transcriptional programs is a critical mechanism that controls hematopoietic development, and, thus, aberrant expression patterns or mutations in epigenetic regulators occur frequently in hematologic malignancies. We demonstrate that the Polycomb protein L3MBTL1, which is monoallelically deleted in 20q- myeloid malignancies, represses the ability of stem cells to drive hematopoietic-specific transcriptional programs by regulating the expression of SMAD5 and impairing its recruitment to target regulatory regions. Indeed, knockdown of L3MBTL1 promotes the development of hematopoiesis and impairs neural cell fate in human pluripotent stem cells. We also found a role for L3MBTL1 in regulating SMAD5 target gene expression in mature hematopoietic cell populations, thereby affecting erythroid differentiation. Taken together, we have identified epigenetic priming of hematopoietic-specific transcriptional networks, which may assist in the development of therapeutic approaches for patients with anemia.
•L3MBTL1 is a chromatin-binding protein that represses SMAD5 expression•Lack of L3MBTL1 primes the hematopoietic development of pluripotent stem cells•L3MBTL1 regulates erythroid differentiation
Nimer, Perna, and colleagues report that the Polycomb protein L3MBTL1, deleted in 20q- myeloid malignancies, plays a regulatory role during the differentiation of human iPS cells. Upon lentiviral knockdown of L3MBTL1, iPS cells stay poised for hematopoietic development and upregulate BMP/SMAD target gene expression, leading to impaired neural development.
The mutually exclusive pattern of the major driver oncogenes in lung cancer suggests that other mutually exclusive oncogenes exist. We performed a systematic search for tyrosine kinase (TK) fusions by screening all TKs for aberrantly high RNA expression levels of the 3′ kinase domain (KD) exons relative to more 5′ exons.
We studied 69 patients (including 5 never smokers and 64 current or former smokers) with lung adenocarcinoma negative for all major mutations in KRAS, EGFR, BRAF, MEK1, and HER2, and for ALK fusions (termed “pan-negative”). A NanoString-based assay was designed to query the transcripts of 90 TKs at two points: 5′ to the KD and within the KD or 3′ to it. Tumor RNAs were hybridized to the NanoString probes and analyzed for outlier 3′ to 5′ expression ratios. Presumed novel fusion events were studied by rapid amplification of cDNA ends (RACE) and confirmatory RT-PCR and FISH.
We identified 1 case each of aberrant 3′ to 5′ ratios in ROS1 and RET. RACE isolated a GOPC-ROS1 (FIG-ROS1) fusion in the former and a KIF5B-RET fusion in the latter, both confirmed by RT-PCR. The RET rearrangement was also confirmed by FISH. The KIF5B-RET patient was one of only 5 never smokers in this cohort.
The KIF5B-RET fusion defines an additional subset of lung cancer with a potentially targetable driver oncogene enriched in never smokers with “pan-negative” lung adenocarcinomas. We also report for the first time in lung cancer the GOPC-ROS1 fusion previously characterized in glioma.
lung cancer; kinase; gene fusion; RET; ROS1; ALK
To investigate the relationship between LDH-A expression, lactate concentration, cell metabolism and metastases in murine 4T1 breast tumors.
Inhibition of LDH-A expression and protein levels were achieved in a metastatic breast cancer cell line (4T1) using shRNA technology. The relationship between tumor LDH-A protein levels and lactate concentration (measured by magnetic resonance spectroscopic imaging-MRSI) and metastases was assessed.
LDH-A knockdown cells (KD9) showed a significant reduction in LDH-A protein and LDH activity, less acid production, decreased transwell migration and invasion, lower proliferation, reduced glucose utilization and glycolysis and increase in oxygen consumption, ROS and cellular ATP levels, compared to control (NC) cells cultured in 25 mM glucose. In vivo studies showed lower lactate levels in KD9, KD5, KD317 tumors than in NC or 4T1 wild-type tumors (p<0.01), and a linear relationship between tumor LDH-A protein expression and lactate concentration. Metastases were delayed and primary tumor growth rate decreased.
We show for the first time that LDH-A knockdown inhibited the formation of metastases, and was accompanied by in vivo changes in tumor cell metabolism. Lactate MRSI can be used as a surrogate to monitor targeted inhibition of LDH-A in a pre-clinical setting and provides a non-invasive imaging strategy to monitor LDH-A targeted therapy. This imaging strategy can be translated to the clinic to identify and monitor patients who are at high risk of developing metastatic disease.
Breast cancer; Metastases; Lactate; LDH-A; MRI; MRS
Biomarkers based on detecting prostate cancer-specific transcripts are associated with inferior outcomes, but their validation in a clinical context is lacking.
To determine whether detecting prostate cancer enhanced transcripts in whole blood using an analytically valid assay has prognostic significance relative to circulating tumor cell (CTC) enumeration.
Design, Setting, and Participants
The predictive value for overall survival of the detection in whole blood by reverse transcription real-time polymerase chain reaction (RT-PCR) of KLK3, KLK2, HOXB13, GRHL2, and FOXA1 was studied in 97 men with metastatic castration-resistant prostate cancer (mCRPC).
2.5ml of blood was collected in PAXgene tubes for total RNA extraction and 7.5 ml for CTC enumeration from patients with progressive mCRPC.
Outcome Measurements and Statistical Analysis
Prostate cancer enriched genes were detected using a sensitive RT-PCR assay in whole blood from patients with mCRPC. Analytical validity of the assay was established in a clinical laboratory environment. The frequency of detecting transcripts was compared to CTC enumeration using CellSearch® in an independent data set and survival associations were explored by concordance probability estimate (CPE).
Results and Limitations
Two or more genes were detected by PCR in 53% (51 of 97, 95% CI 43–63%) of patients, and unfavorable CTC counts (≥5cells) were seen in 46% (45 of 97, 95% CI 36–56%). Importantly, transcripts were detectable in 11 of 52 patients with favorable CTC counts (21%, 95% CI 8–35%). Transcript detection predicted overall survival in a proportional hazards model. Significantly, the predictive accuracy of RT-PCR detection in combination with CTC enumeration had a CPE of 0.752 (SE=0.038), although limited by the number of patients.
This validated RT-PCR assay detecting prostate-specific RNA in whole blood is prognostic for survival, and may assess patient risk complimentary with CellSearch CTC enumeration. Its clinical utility is being prospectively explored.
biomarker; circulating tumor cells; prostate cancer; prostate-specific markers
To profile the subgingival oral microbiota abundance and diversity in never-treated, new-onset rheumatoid arthritis (NORA) patients.
Periodontal disease (PD) status, clinical activity and sociodemographic factors were determined in patients with NORA, chronic RA (CRA) and healthy subjects. Massively parallel pyrosequencing was used to compare the composition of subgingival microbiota and establish correlations between presence/abundance of bacteria and disease phenotypes. Anti-P. gingivalis antibodies were tested to assess prior exposure.
The more advanced forms of periodontitis are already present at disease onset in NORA patients. The subgingival microbiota of NORA is distinct from controls. In most cases, however, these differences can be attributed to PD severity and are not inherent to RA. The presence and abundance of P. gingivalis is directly associated with PD severity as well, is not unique to RA, and does not correlate with anti-citrullinated peptide antibody (ACPA) titers. Overall exposure to P. gingivalis is similar in RA and controls, observed in 78.4% and 83.3%, respectively. Anaeroglobus geminatus correlated with ACPA/RF presence. Prevotella and Leptotrichia species are the only characteristic taxa in the NORA group irrespective of PD status.
NORA patients exhibit a high prevalence of PD at disease onset, despite their young age and paucity of smoking history. The subgingival microbiota of NORA patients is similar to CRA and healthy subjects of comparable PD severity. Although colonization with P. gingivalis correlates with PD severity, overall exposure is similar among groups. The role of A. geminatus and Prevotella/Leptotrichia species in this process merits further study.
A large number of computational methods have been developed for analyzing differential gene expression in RNA-seq data. We describe a comprehensive evaluation of common methods using the SEQC benchmark dataset and ENCODE data. We consider a number of key features, including normalization, accuracy of differential expression detection and differential expression analysis when one condition has no detectable expression. We find significant differences among the methods, but note that array-based methods adapted to RNA-seq data perform comparably to methods designed for RNA-seq. Our results demonstrate that increasing the number of replicate samples significantly improves detection power over increased sequencing depth.
Angiosarcomas (AS) represent a heterogeneous group of malignant vascular tumors that may occur spontaneously as primary tumors or secondarily after radiation therapy or in the context of chronic lymphedema. Most secondary AS have been associated with MYC oncogene amplification, while the role of MYC abnormalities in primary AS is not well defined. Twenty-two primary and secondary AS were analyzed by array-comparative genomic hybridization (aCGH) and by deep sequencing of small RNA libraries. By aCGH and subsequently confirmed by FISH, MYC amplification was identified in three of six primary tumors and in eight of 12 secondary AS. We have also found MAML1 as a new potential oncogene in MYC-amplified AS. Significant up-regulation of the miR17-92 cluster was observed in MYC-amplified AS compared to AS lacking MYC amplification and the control group (other vascular tumors, non-vascular sarcomas). Moreover, MYC-amplified AS were associated with a significantly lower expression of thrombospondin-1 (THBS1) than AS without MYC amplification or controls. Altogether, our study implicates MYC amplification not only in the pathogenesis of secondary AS but also in a subset of primary AS. Thus, MYC amplification may play a crucial role in the angiogenic phenotype of AS through up-regulation of the miR-17-92 cluster, which subsequently downregulates THBS1, a potent endogenous inhibitor of angiogenesis.
Cancer gene fusions that encode a chimeric protein are often characterized by an intragenic discontinuity in the RNA expression levels of the exons that are 5′ or 3′ to the fusion point in one or both of the fusion partners due to differences in the levels of activation of their respective promoters. Based on this, we developed an unbiased, genome-wide bioinformatic screen for gene fusions using Affymetrix Exon array expression data. Using a training set of 46 samples with different known gene fusions, we developed a data analysis pipeline, the “Fusion Score (FS) model”, to score and rank genes for intragenic changes in expression. In a separate discovery set of 41 tumor samples with possible unknown gene fusions, the FS model generated a list of 552 candidate genes. The transcription factor gene NCOA2 was one of the candidates identified in a mesenchymal chondrosarcoma. A novel HEY1-NCOA2 fusion was identified by 5′ RACE, representing an in-frame fusion of HEY1 exon 4 to NCOA2 exon 13. RT-PCR or FISH evidence of this HEY1-NCOA2 fusion was present in all additional mesenchymal chondrosarcomas tested with a definitive histologic diagnosis and adequate material for analysis (n=9) but was absent in 15 samples of other subtypes of chondrosarcomas. We also identified a NUP107-LGR5 fusion in a dedifferentiated liposarcoma but analysis of 17 additional samples did not confirm it as a recurrent event in this sarcoma type. The novel HEY1-NCOA2 fusion appears to be the defining and diagnostic gene fusion in mesenchymal chondrosarcomas.
Differences between TLR-deficient mouse colonies occur from extended husbandry in isolation that are communicated to offspring by maternal transmission.
The intestinal microbiota contributes to the development of the immune system, and conversely, the immune system influences the composition of the microbiota. Toll-like receptors (TLRs) in the gut recognize bacterial ligands. Although TLR signaling represents a major arm of the innate immune system, the extent to which TLRs influence the composition of the intestinal microbiota remains unclear. We performed deep 16S ribosomal RNA sequencing to characterize the complex bacterial populations inhabiting the ileum and cecum of TLR- and MyD88-deficient mice. The microbiota of MyD88- and TLR-deficient mouse colonies differed markedly, with each colony harboring distinct and distinguishable bacterial populations in the small and large intestine. Comparison of MyD88-, TLR2-, TLR4-, TLR5-, and TLR9-deficient mice and their respective wild-type (WT) littermates demonstrated that the impact of TLR deficiency on the composition of the intestinal microbiota is minimal under homeostatic conditions and after recovery from antibiotic treatment. Thus, differences between TLR-deficient mouse colonies reflected long-term divergence of the microbiota after extended husbandry in isolation from each other. Long-term breeding of isolated mouse colonies results in changes of the intestinal microbiota that are communicated to offspring by maternal transmission, which account for marked compositional differences between WT and mutant mouse strains.
GVHD is associated with significant shifts in the composition of the intestinal microbiota in human and mouse models; manipulating the microbiota can alter the severity of GVHD in mice.
Despite a growing understanding of the link between intestinal inflammation and resident gut microbes, longitudinal studies of human flora before initial onset of intestinal inflammation have not been reported. Here, we demonstrate in murine and human recipients of allogeneic bone marrow transplantation (BMT) that intestinal inflammation secondary to graft-versus-host disease (GVHD) is associated with major shifts in the composition of the intestinal microbiota. The microbiota, in turn, can modulate the severity of intestinal inflammation. In mouse models of GVHD, we observed loss of overall diversity and expansion of Lactobacillales and loss of Clostridiales. Eliminating Lactobacillales from the flora of mice before BMT aggravated GVHD, whereas reintroducing the predominant species of Lactobacillus mediated significant protection against GVHD. We then characterized gut flora of patients during onset of intestinal inflammation caused by GVHD and found patterns mirroring those in mice. We also identified increased microbial chaos early after allogeneic BMT as a potential risk factor for subsequent GVHD. Together, these data demonstrate regulation of flora by intestinal inflammation and suggest that flora manipulation may reduce intestinal inflammation and improve outcomes for allogeneic BMT recipients.
Recurrent mutations in isocitrate dehydrogenase 1 (IDH1) and IDH2 have been identified in gliomas, acute myeloid leukaemias (AML) and chondrosarcomas, and share a novel enzymatic property of producing 2-hydroxyglutarate (2HG) from α-ketoglutarate1-6. Here we report that 2HG-producing IDH mutants can prevent the histone demethylation that is required for lineage-specific progenitor cells to differentiate into terminally differentiated cells. In tumour samples from glioma patients, IDH mutations were associated with a distinct gene expression profile enriched for genes expressed in neural progenitor cells, and this was associated with increased histone methylation. To test whether the ability of IDH mutants to promote histone methylation contributes to a block in cell differentiation in non-transformed cells, we tested the effect of neomorphic IDH mutants on adipocyte differentiation in vitro. Introduction of either mutant IDH or cell-permeable 2HG was associated with repression of the inducible expression of lineage-specific differentiation genes and a block to differentiation. This correlated with a significant increase in repressive histone methylation marks without observable changes in promoter DNA methylation. Gliomas were found to have elevated levels of similar histone repressive marks. Stable transfection of a 2HG-producing mutant IDH into immortalized astrocytes resulted in progressive accumulation of histone methylation. Of the marks examined, increased H3K9 methylation reproducibly preceded a rise in DNA methylation as cells were passaged in culture. Furthermore, we found that the 2HG-inhibitable H3K9 demethylase KDM4C was induced during adipocyte differentiation, and that RNA-interference suppression of KDM4C was sufficient to block differentiation. Together these data demonstrate that 2HG can inhibit histone demethylation and that inhibition of histone demethylation can be sufficient to block the differentiation of non-transformed cells.
Liposarcoma remains the most common mesenchymal cancer, with a mortality rate of 60% among patients with this disease. To address the present lack of therapeutic options, we embarked upon a study of microRNA (miRNA) expression alterations associated with liposarcomagenesis with the goal of exploiting differentially expressed miRNAs and the gene products they regulate as potential therapeutic targets. MicroRNA expression was profiled in samples of normal adipose tissue, well-differentiated liposarcoma, and dedifferentiated liposarcoma by both deep sequencing of small RNA libraries and hybridization-based Agilent microarrays. The expression profiles discriminated liposarcoma from normal adipose tissue and well-differentiated from dedifferentiated disease. We defined over 40 miRNAs that were dysregulated in dedifferentiated liposarcomas in both the sequencing and the microarray analysis. The upregulated miRNAs included two cancer-associated species (miR-21, miR-26a), and the downregulated miRNAs included two species that were highly abundant in adipose tissue (miR-143, miR-145). Restoring miR-143 expression in dedifferentiated liposarcoma cells inhibited proliferation, induced apoptosis, and decreased expression of BCL2, TOP2A, PRC1, and PLK1. The downregulation of PRC1 and its docking partner PLK1 suggests that miR-143 inhibits cytokinesis in these cells. In support of this idea, treatment with a PLK1 inhibitor potently induced G2/M growth arrest and apoptosis in liposarcoma cells. Taken together, our findings suggest that miR-143 re-expression vectors or selective agents directed at miR-143 or its targets may have therapeutic value in dedifferentiated liposarcoma.
Gastric cancer may be subdivided into three distinct subtypes –proximal, diffuse, and distal gastric cancer– based on histopathologic and anatomic criteria. Each subtype is associated with unique epidemiology. Our aim is to test the hypothesis that these distinct gastric cancer subtypes may also be distinguished by gene expression analysis.
Patients with localized gastric adenocarcinoma being screened for a phase II preoperative clinical trial (NCI 5917) underwent endoscopic biopsy for fresh tumor procurement. 4–6 targeted biopsies of the primary tumor were obtained. Macrodissection was performed to ensure >80% carcinoma in the sample. HG-U133A GeneChip (Affymetrix) was used for cDNA expression analysis, and all arrays were processed and analyzed using the Bioconductor R-package.
Between November 2003 and January 2006, 57 patients were screened to identify 36 patients with localized gastric cancer who had adequate RNA for expression analysis. Using supervised analysis, we built a classifier to distinguish the three gastric cancer subtypes, successfully classifying each into tightly grouped clusters. Leave-one-out cross validation error was 0.14, suggesting that >85% of samples were classified correctly. Gene set analysis with the False Discovery Rate set at 0.25 identified several pathways that were differentially regulated when comparing each gastric cancer subtype to adjacent normal stomach.
Subtypes of gastric cancer that have epidemiologic and histologic distinction are also distinguished by gene expression data. These preliminary data suggest a new classification of gastric cancer with implications for improving our understanding of disease biology and identification of unique molecular drivers for each gastric cancer subtype.
Gastric Cancer; cDNA Expression; classification; pathways
Uveal melanoma is the most common intraocular cancer. There are no effective therapies for metastatic disease. Mutations in GNAQ, the gene encoding an alpha subunit of heterotrimeric G proteins, are found in 40% of uveal melanomas.
We sequenced exon 5 of GNAQ and GNA11, a paralogue of GNAQ, in 713 melanocytic neoplasms of different types (186 uveal melanomas, 139 blue nevi, 106 other nevi, and 282 other melanomas). We sequenced exon 4 of GNAQ and GNA11 in 453 of these samples and in all coding exons of GNAQ and GNA11 in 97 uveal melanomas and 45 blue nevi.
We found somatic mutations in exon 5 (affecting Q209) and in exon 4 (affecting R183) in both GNA11 and GNAQ, in a mutually exclusive pattern. Mutations affecting Q209 in GNA11 were present in 7% of blue nevi, 32% of primary uveal melanomas, and 57% of uveal melanoma metastases. In contrast, we observed Q209 mutations in GNAQ in 55% of blue nevi, 45% of uveal melanomas, and 22% of uveal melanoma metastases. Mutations affecting R183 in either GNAQ or GNA11 were less prevalent (2% of blue nevi and 6% of uveal melanomas) than the Q209 mutations. Mutations in GNA11 induced spontaneously metastasizing tumors in a mouse model and activated the mitogen-activated protein kinase pathway.
Of the uveal melanomas we analyzed, 83% had somatic mutations in GNAQ or GNA11. Constitutive activation of the pathway involving these two genes appears to be a major contributor to the development of uveal melanoma. (Funded by the National Institutes of Health and others.)
Prior to gastrulation in the mouse, all endodermal cells arise from the primitive
endoderm of the blastocyst stage embryo. Primitive endoderm and its derivatives
are generally referred to as extra-embryonic endoderm (ExEn) because the
majority of these cells contribute to extra-embryonic lineages encompassing the
visceral endoderm (VE) and the parietal endoderm (PE). During gastrulation, the
definitive endoderm (DE) forms by ingression of cells from the epiblast. The DE
comprises most of the cells of the gut and its accessory organs. Despite their
different origins and fates, there is a surprising amount of overlap in marker
expression between the ExEn and DE, making it difficult to distinguish between
these cell types by marker analysis. This is significant for two main reasons.
First, because endodermal organs, such as the liver and pancreas, play important
physiological roles in adult animals, much experimental effort has been directed
in recent years toward the establishment of protocols for the efficient
derivation of endodermal cell types in vitro. Conversely,
factors secreted by the VE play pivotal roles that cannot be attributed to the
DE in early axis formation, heart formation and the patterning of the anterior
nervous system. Thus, efforts in both of these areas have been hampered by a
lack of markers that clearly distinguish between ExEn and DE. To further
understand the ExEn we have undertaken a comparative analysis of three ExEn-like
cell lines (END2, PYS2 and XEN). PYS2 cells are derived from embryonal
carcinomas (EC) of 129 strain mice and have been characterized as parietal
endoderm-like , END2 cells are derived from P19 ECs and
described as visceral endoderm-like, while XEN cells are derived from blastocyst
stage embryos and are described as primitive endoderm-like. Our analysis
suggests that none of these cell lines represent a bona fide
single in vivo lineage. Both PYS2 and XEN cells represent mixed
populations expressing markers for several ExEn lineages. Conversely END2 cells,
which were previously characterized as VE-like, fail to express many markers
that are widely expressed in the VE, but instead express markers for only a
subset of the VE, the anterior visceral endoderm. In addition END2 cells also
express markers for the PE. We extended these observations with microarray
analysis which was used to probe and refine previously published data sets of
genes proposed to distinguish between DE and VE. Finally, genome-wide pathway
analysis revealed that SMAD-independent TGFbeta signaling through a TAK1/p38/JNK
or TAK1/NLK pathway may represent one mode of intracellular signaling shared by
all three of these lines, and suggests that factors downstream of these pathways
may mediate some functions of the ExEn. These studies represent the first step
in the development of XEN cells as a powerful molecular genetic tool to study
the endodermal signals that mediate the important developmental functions of the
extra-embryonic endoderm. Our data refine our current knowledge of markers that
distinguish various subtypes of endoderm. In addition, pathway analysis suggests
that the ExEn may mediate some of its functions through a non-classical MAP
Kinase signaling pathway downstream of TAK1.
Owing to the dynamic nature of the transcriptome, gene expression profiling is a promising tool for discovery of disease-related genes and biological pathways. In the present study, we examined gene expression in whole blood of 12 patients with CAD (coronary artery disease) and 12 healthy control subjects. Furthermore, ten patients with CAD underwent whole-blood gene expression analysis before and after the completion of a cardiac rehabilitation programme following surgical coronary revascularization. mRNA and miRNA (microRNA) were isolated for expression profiling. Gene expression analysis identified 365 differentially expressed genes in patients with CAD compared with healthy controls (175 up- and 190 down-regulated in CAD), and 645 in CAD rehabilitation patients (196 up- and 449 down-regulated post-rehabilitation). Biological pathway analysis identified a number of canonical pathways, including oxidative phosphorylation and mitochondrial function, as being significantly and consistently modulated across the groups. Analysis of miRNA expression revealed a number of differentially expressed miRNAs, including hsa-miR-140-3p (control compared with CAD, P=0.017), hsa-miR-182 (control compared with CAD, P=0.093), hsa-miR-92a and hsa-miR-92b (post- compared with pre-exercise, P<0.01). Global analysis of predicted miRNA targets found significantly reduced expression of genes with target regions compared with those without: hsa-miR-140-3p (P=0.002), hsa-miR-182 (P=0.001), hsa-miR-92a and hsa-miR-92b (P=2.2×10−16). In conclusion, using whole blood as a ‘surrogate tissue’ in patients with CAD, we have identified differentially expressed miRNAs, differentially regulated genes and modulated pathways which warrant further investigation in the setting of cardiovascular function. This approach may represent a novel non-invasive strategy to unravel potentially modifiable pathways and possible therapeutic targets in cardiovascular disease.
coronary artery disease (CAD); gene expression; microRNA (miRNA); mitochondrion; oxidative phosphorylation; rehabilitation programme; ATP5I, ATP synthase, H+ transporting, mitochondrial F0 complex, subunit E; ATP5L, ATP synthase, H+ transporting, mitochondrial F0 complex, subunit G; CABG, coronary artery bypass graft; CAD, coronary artery disease; CASP3, caspase 3; COX7C, cytochrome c oxidase subunit VIIc, CRP, C-reactive protein; FDR, False Discovery Rate; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; miRNA, microRNA; NDUFA1, NADH dehydrogenase (ubiquinone) 1α subcomplex 1; NDUFB3, NADH dehydrogenase (ubiquinone) 1β subcomplex 3; qRT-PCR, quantitative real-time PCR; ROS, reactive oxygen species; UQCRQ, ubiquinol-cytochrome c reductase, complex III subunit VII; UTR, untranslated region
In many approaches to the inference and modeling of regulatory interactions using microarray data, the expression of the gene coding for the transcription factor is considered to be an accurate surrogate for the true activity of the protein it produces. There are many instances where this is inaccurate due to post-translational modifications of the transcription factor protein. Inference of the activity of the transcription factor from the expression of its targets has predominantly involved linear models that do not reflect the nonlinear nature of transcription. We extend a recent approach to inferring the transcription factor activity based on nonlinear Michaelis-Menten kinetics of transcription from maximum likelihood to fully Bayesian inference and give an example of how the model can be further developed.
We present results on synthetic and real microarray data. Additionally, we illustrate how gene and replicate specific delays can be incorporated into the model.
We demonstrate that full Bayesian inference is appropriate in this application and has several benefits over the maximum likelihood approach, especially when the volume of data is limited. We also show the benefits of using a non-linear model over a linear model, particularly in the case of repression.
Angiosarcomas (ASs) represent a heterogeneous group of malignant vascular tumors that may occur spontaneously as primary tumors or secondarily after radiation therapy or in the context of chronic lymphedema. Most secondary ASs have been associated with MYC oncogene amplification, whereas the role of MYC abnormalities in primary AS is not well defined. Twenty-two primary and secondary ASs were analyzed by array-comparative genomic hybridization (aCGH) and by deep sequencing of small RNA libraries. By aCGH and subsequently confirmed by fluorescence in situ hybridization, MYC amplification was identified in three out of six primary tumors and in 8 out of 12 secondary AS. We have also found MAML1 as a new potential oncogene in MYC-amplified AS. Significant upregulation of the miR-17-92 cluster was observed in MYC-amplified AS compared to AS lacking MYC amplification and the control group (other vascular tumors, nonvascular sarcomas). Moreover, MYC-amplified ASs were associated with a significantly lower expression of thrombospondin-1 (THBS1) than AS without MYC amplification or controls. Altogether, our study implicates MYC amplification not only in the pathogenesis of secondary AS but also in a subset of primary AS. Thus, MYC amplification may play a crucial role in the angiogenic phenotype of AS through upregulation of the miR-17-92 cluster, which subsequently downregulates THBS1, a potent endogenous inhibitor of angiogenesis. © 2012 Wiley Periodicals, Inc.