An evolutionarily conserved enhancer regulates Bmp4 expression in developing incisor and limb bud.

To elucidate the transcriptional regulation of Bmp4 expression during organogenesis, we used phylogenetic footprinting and transgenic reporter analyses to identify Bmp4 cis-regulatory modules (CRMs). These analyses identified a regulatory region located ∼46 kb upstream of the mouse Bmp4 transcription start site that had previously been shown to direct expression in lateral plate mesoderm. We refined this regulatory region to a 396-bp minimal enhancer, and show that it recapitulates features of endogenous Bmp4 expression in developing mandibular arch ectoderm and incisor epithelium during the initiation-stage of tooth development. In addition, this enhancer directs expression in the apical ectodermal ridge (AER) of the developing limb and in anterior and posterior limb mesenchyme. Transcript profiling of E11.5 mouse incisor dental lamina, together with protein binding microarray (PBM) analyses, allowed identification of a conserved DNA binding motif in the Bmp4 enhancer for Pitx homeoproteins, which are also expressed in the developing mandibular and incisor epithelium. In vitro electrophoretic mobility shift assays (EMSA) and in vivo transgenic reporter mutational analyses revealed that this site supports Pitx binding and that the site is necessary to recapitulate aspects of endogenous Bmp4 expression in developing craniofacial and limb tissues. Finally, Pitx2 chromatin immunoprecipitation (ChIP) demonstrated direct binding of Pitx2 to this Bmp4 enhancer site in a dental epithelial cell line. These results establish a direct molecular regulatory link between Pitx family members and Bmp4 gene expression in developing incisor epithelium.

Delineation of a cellular hierarchy in lung cancer reveals an oncofetal antigen expressed on tumor-initiating cells.

Poorly differentiated tumors in non-small cell lung cancer (NSCLC) have been associated with shorter patient survival and shorter time to recurrence following treatment. Here, we integrate multiple experimental models with clinicopathologic analysis of patient tumors to delineate a cellular hierarchy in NSCLC. We show that the oncofetal protein 5T4 is expressed on tumor-initiating cells and associated with worse clinical outcome in NSCLC. Coexpression of 5T4 and factors involved in the epithelial-to-mesenchymal transition were observed in undifferentiated but not in differentiated tumor cells. Despite heterogeneous expression of 5T4 in NSCLC patient-derived xenografts, treatment with an anti-5T4 antibody-drug conjugate resulted in complete and sustained tumor regression. Thus, the aggressive growth of heterogeneous solid tumors can be blocked by therapeutic agents that target a subpopulation of cells near the top of the cellular hierarchy.

Dynamic changes in the microRNA expression profile reveal multiple regulatory mechanisms in the spinal nerve ligation model of neuropathic pain.

Neuropathic pain resulting from nerve lesions or dysfunction represents one of the most challenging neurological diseases to treat. A better understanding of the molecular mechanisms responsible for causing these maladaptive responses can help develop novel therapeutic strategies and biomarkers for neuropathic pain. We performed a miRNA expression profiling study of dorsal root ganglion (DRG) tissue from rats four weeks post spinal nerve ligation (SNL), a model of neuropathic pain. TaqMan low density arrays identified 63 miRNAs whose level of expression was significantly altered following SNL surgery. Of these, 59 were downregulated and the ipsilateral L4 DRG, not the injured L5 DRG, showed the most significant downregulation suggesting that miRNA changes in the uninjured afferents may underlie the development and maintenance of neuropathic pain. TargetScan was used to predict mRNA targets for these miRNAs and it was found that the transcripts with multiple predicted target sites belong to neurologically important pathways. By employing different bioinformatic approaches we identified neurite remodeling as a significantly regulated biological pathway, and some of these predictions were confirmed by siRNA knockdown for genes that regulate neurite growth in differentiated Neuro2A cells. In vitro validation for predicted target sites in the 3'-UTR of voltage-gated sodium channel Scn11a, alpha 2/delta1 subunit of voltage-dependent Ca-channel, and purinergic receptor P2rx ligand-gated ion channel 4 using luciferase reporter assays showed that identified miRNAs modulated gene expression significantly. Our results suggest the potential for miRNAs to play a direct role in neuropathic pain.

Leukemic IDH1 and IDH2 mutations result in a hypermethylation phenotype, disrupt TET2 function, and impair hematopoietic differentiation.

Cancer-associated IDH mutations are characterized by neomorphic enzyme activity and resultant 2-hydroxyglutarate (2HG) production. Mutational and epigenetic profiling of a large acute myeloid leukemia (AML) patient cohort revealed that IDH1/2-mutant AMLs display global DNA hypermethylation and a specific hypermethylation signature. Furthermore, expression of 2HG-producing IDH alleles in cells induced global DNA hypermethylation. In the AML cohort, IDH1/2 mutations were mutually exclusive with mutations in the α-ketoglutarate-dependent enzyme TET2, and TET2 loss-of-function mutations were associated with similar epigenetic defects as IDH1/2 mutants. Consistent with these genetic and epigenetic data, expression of IDH mutants impaired TET2 catalytic function in cells. Finally, either expression of mutant IDH1/2 or Tet2 depletion impaired hematopoietic differentiation and increased stem/progenitor cell marker expression, suggesting a shared proleukemogenic effect.

An in silico analysis of microRNAs: mining the miRNAome.

Systematic analysis of literature- and experimentally-derived datasets using text mining with ontological enrichment and network modeling revealed global trends in the microRNA (miRNA) interactome. A total of 756 unique miRNAs were resolved from PubMed abstracts and 1165 direct relationships between 270 miRNAs and 581 genes were identified as phrase groups using semantic search techniques. These miRNA:gene interactions were built into a bipartite network (the miRNAome) which displays scale-free degree distribution. Functional classification of miRNA-target genes using PANTHER revealed 189 distinct molecular functions, with significant enrichment of nucleic acid binding, transcription and protein phosphorylation. Pathway analysis revealed a network of 176 miRNAs linked to 368 OMIM disorders via their target genes, which are enriched (p = 0.0047) for disease-associated SNP variations. Reference to a database of drug targets revealed that 24.8% of all published miRNA-targets are targets for drug development programs, while a sub-set (18.2%) are targets for FDA-approved drugs. Consistent with topological analysis of the miRNA-disease network, the most prevalent class of FDA-approved drugs is anti-neoplastic agents against published miRNA-target genes. Linking miRNAs to biological process and diseases reveals distinct co-regulation of phenotypes that could aid in understanding the role miRNA-based gene regulation plays in biological phenomena.

Comparison of human and rat uterine leiomyomata: identification of a dysregulated mammalian target of rapamycin pathway.

Uterine leiomyomata, or fibroids, are benign tumors of the uterine myometrium that significantly affect up to 30% of reproductive-age women. Despite being the primary cause of hysterectomy in the United States, accounting for up to 200,000 procedures annually, the etiology of leiomyoma remains largely unknown. As a basis for understanding leiomyoma pathogenesis and identifying targets for pharmacotherapy, we conducted transcriptional profiling of leiomyoma and unaffected myometrium from humans and Eker rats, the best characterized preclinical model of leiomyomata. A global comparison of mRNA from leiomyoma versus myometrium in human and rat identified a highly significant overlap of dysregulated gene expression in leiomyomata. An unbiased pathway analysis using a method of gene-set enrichment based on the sigPathway algorithm detected the mammalian target of rapamycin (mTOR) pathway as one of the most highly up-regulated pathways in both human and rat tumors. To validate this pathway as a therapeutic target for uterine leiomyomata, preclinical studies were conducted in Eker rats. These rats develop uterine leiomyomata as a consequence of loss of Tsc2 function and up-regulation of mTOR signaling. Inhibition of mTOR in female Eker rats with the rapamycin analogue WAY-129327 for 2 weeks decreased mTOR signaling and cell proliferation in tumors, and treatment for 4 months significantly decreased tumor incidence, multiplicity, and size. These results identify dysregulated mTOR signaling as a component of leiomyoma etiology across species and directly show the dependence of uterine leiomyomata with activated mTOR on this signaling pathway for growth.

Efficacy and pharmacodynamic effects of bosutinib (SKI-606), a Src/Abl inhibitor, in freshly generated human pancreas cancer xenografts.

Recently, Src tyrosine kinase has emerged as an attractive target for anticancer therapy, and Src is overexpressed in pancreatic cancer. The purpose of the study was to investigate the in vivo efficacy and pharmacodynamic effects of bosutinib (SKI-606), a Src/Abl inhibitor, using a panel of human pancreatic tumor xenografts. Surgically resected human pancreatic tumors were implanted into female nude mice and randomized to bosutinib versus control. Src and other pathways were analyzed by Western Blot, IHC, and Affymetrix U133 Plus 2.0 gene arrays. Of 15 patient tumors, 3 patient tumors were found to be sensitive to bosutinib, defined as tumor growth of <45% than that of control tumors. There were no definite differences between sensitive and resistant tumors in the baseline Src kinase pathway protein expression assessed by Western Blot. Caveolin-1 expression, as assessed by reverse transcription-PCR and immunohistochemistry, was frequently higher in sensitive cases. In sensitive tumors, bosutinib resulted in increased apoptosis. Phosphorylation of key signaling molecules downstream of Src, signal transducers and activators of transcription 3, and signal transducers and activators of transcription 3, were significantly inhibited by bosutinib. K-Top Scoring Pairs analysis of gene arrays gave a six-gene classifier that predicted resistance versus sensitivity in six validation cases. These results may aid the clinical development of bosutinib and other Src inhibitors in pancreas cancer.

Molecular analysis of the vaginal response to estrogens in the ovariectomized rat and postmenopausal woman.

BACKGROUND: Vaginal atrophy (VA) is the thinning of the vaginal epithelial lining, typically the result of lowered estrogen levels during menopause. Some of the consequences of VA include increased susceptibility to bacterial infection, pain during sexual intercourse, and vaginal burning or itching. Although estrogen treatment is highly effective, alternative therapies are also desired for women who are not candidates for post-menopausal hormone therapy (HT). The ovariectomized (OVX) rat is widely accepted as an appropriate animal model for many estrogen-dependent responses in humans; however, since reproductive biology can vary significantly between mammalian systems, this study examined how well the OVX rat recapitulates human biology. METHODS: We analyzed 19 vaginal biopsies from human subjects pre and post 3-month 17beta-estradiol treated by expression profiling. Data were compared to transcriptional profiling generated from vaginal samples obtained from ovariectomized rats treated with 17beta-estradiol for 6 hrs, 3 days or 5 days. The level of differential expression between pre- vs. post- estrogen treatment was calculated for each of the human and OVX rat datasets. Probe sets corresponding to orthologous rat and human genes were mapped to each other using NCBI Homologene. RESULTS: A positive correlation was observed between the rat and human responses to estrogen. Genes belonging to several biological pathways and GO categories were similarly differentially expressed in rat and human. A large number of the coordinately regulated biological processes are already known to be involved in human VA, such as inflammation, epithelial development, and EGF pathway activation. CONCLUSION: At the transcriptional level, there is evidence of significant overlap of the effects of estrogen treatment between the OVX rat and human VA samples.

An integrated strategy for analyzing the unique developmental programs of different myoblast subtypes.

An important but largely unmet challenge in understanding the mechanisms that govern the formation of specific organs is to decipher the complex and dynamic genetic programs exhibited by the diversity of cell types within the tissue of interest. Here, we use an integrated genetic, genomic, and computational strategy to comprehensively determine the molecular identities of distinct myoblast subpopulations within the Drosophila embryonic mesoderm at the time that cell fates are initially specified. A compendium of gene expression profiles was generated for primary mesodermal cells purified by flow cytometry from appropriately staged wild-type embryos and from 12 genotypes in which myogenesis was selectively and predictably perturbed. A statistical meta-analysis of these pooled datasets--based on expected trends in gene expression and on the relative contribution of each genotype to the detection of known muscle genes--provisionally assigned hundreds of differentially expressed genes to particular myoblast subtypes. Whole embryo in situ hybridizations were then used to validate the majority of these predictions, thereby enabling true-positive detection rates to be estimated for the microarray data. This combined analysis reveals that myoblasts exhibit much greater gene expression heterogeneity and overall complexity than was previously appreciated. Moreover, it implicates the involvement of large numbers of uncharacterized, differentially expressed genes in myogenic specification and subsequent morphogenesis. These findings also underscore a requirement for considerable regulatory specificity for generating diverse myoblast identities. Finally, to illustrate how the developmental functions of newly identified myoblast genes can be efficiently surveyed, a rapid RNA interference assay that can be scored in living embryos was developed and applied to selected genes. This integrated strategy for examining embryonic gene expression and function provides a substantially expanded framework for further studies of this model developmental system.

Mutant-specific gene programs in the zebrafish.

Hematopoiesis involves the production of stem cells, followed by the orchestrated differentiation of the blood lineages. Genetic screens in zebrafish have identified mutants with defects that disrupt specific stages of hematopoiesis and vasculogenesis, including the cloche, spadetail (tbx16), moonshine (tif1g), bloodless, and vlad tepes (gata1) mutants. To better characterize the blood program, gene expression profiling was carried out in these mutants and in scl-morphants (scl(mo)). Distinct gene clusters were demarcated by stage-specific and mutant-specific gene regulation. These were found to correlate with the transcriptional program of hematopoietic progenitor cells, as well as of the erythroid, myeloid, and vascular lineages. Among these, several novel hematopoietic and vascular genes were detected, for instance, the erythroid transcription factors znfl2 and ncoa4. A specific regulation was found for myeloid genes, as they were more strongly expressed in vlt mutants compared with other erythroid mutants. A unique gene expression pattern of up-regulated isoprenoid synthesis genes was found in cloche and scl(mo), possibly in migrating cells. In conjunction with the high conservation of vertebrate hematopoiesis, the comparison of transcriptional profiles in zebrafish blood mutants represents a versatile and powerful tool to elucidate the genetic regulation of blood and blood vessel development.

Preferred analysis methods for Affymetrix GeneChips revealed by a wholly defined control dataset.

BACKGROUND: As more methods are developed to analyze RNA-profiling data, assessing their performance using control datasets becomes increasingly important. RESULTS: We present a 'spike-in' experiment for Affymetrix GeneChips that provides a defined dataset of 3,860 RNA species, which we use to evaluate analysis options for identifying differentially expressed genes. The experimental design incorporates two novel features. First, to obtain accurate estimates of false-positive and false-negative rates, 100-200 RNAs are spiked in at each fold-change level of interest, ranging from 1.2 to 4-fold. Second, instead of using an uncharacterized background RNA sample, a set of 2,551 RNA species is used as the constant (1x) set, allowing us to know whether any given probe set is truly present or absent. Application of a large number of analysis methods to this dataset reveals clear variation in their ability to identify differentially expressed genes. False-negative and false-positive rates are minimized when the following options are chosen: subtracting nonspecific signal from the PM probe intensities; performing an intensity-dependent normalization at the probe set level; and incorporating a signal intensity-dependent standard deviation in the test statistic. CONCLUSIONS: A best-route combination of analysis methods is presented that allows detection of approximately 70% of true positives before reaching a 10% false-discovery rate. We highlight areas in need of improvement, including better estimate of false-discovery rates and decreased false-negative rates.

Identification of a novel set of genes regulated by a unique liver X receptor-alpha -mediated transcription mechanism.

We have reported previously that liver X receptor-alpha (LXRalpha) can mediate a novel cAMP-dependent increase in renin and c-myc gene transcription by binding as a monomer to a unique regulatory element termed the cAMP-negative response element (CNRE). To determine whether this novel action of LXRalpha has global implications on gene regulation, we employed expression profiling to identify other genes regulated by this unique mechanism. Here we report the existence of a set of known and unknown transcripts regulated in parallel with renin. Querying the Celera Mouse Genome Assembly revealed that a majority of these genes contained the consensus CNRE. We have confirmed the functionality of these CNREs by competition for LXRalpha binding via electrophoretic mobility shift assays (EMSA) and by the use of CNRE decoy molecules documenting the abolishment of the cAMP-mediated gene induction. Taken together, these results demonstrate that the interaction between cAMP-activated LXRalpha and the CNRE enhancer element is responsible for widespread changes in gene expression and identify a set of LXRalpha/cAMP-regulated genes that may have important biological implications.

Gene expression profiling reveals progesterone-mediated cell cycle and immunoregulatory roles of Hoxa-10 in the preimplantation uterus.

Human infertility and recurrent pregnancy loss caused by implantation defects are poorly understood. Hoxa-10-deficient female mice have severe infertility and recurrent pregnancy loss due to defective uterine implantation. Gene expression profiling experiments reveal that Hoxa-10 is an important regulator of two critical events in implantation: stromal cell proliferation and local immunosuppression. At the time of implantation, Hoxa-10 mediates the progesterone-stimulated proliferation of uterine stromal cells. Hoxa-10 mutants express a stromal cell proliferation defect that is accompanied by quantitative or spatial alterations in the expression of two cyclin-dependent kinase inhibitor genes, p57 and p15. Hoxa-10 deficiency also leads to a severe local immunological disturbance, characterized by a polyclonal proliferation of T cells, that occurs in place of the normal progesterone-mediated immunosuppression in the periimplantation uterus.