Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture

Purpose: Limbal stem cell deficiency is a challenging clinical problem and the current treatment involves replenishing the depleted limbal stem cell (LSC) pool by either limbal tissue transplantation or use of cultivated limbal epithelial cells (LEC). Our experience of cultivating the LEC on denuded human amniotic membrane using a feeder cell free method, led to identification of mesenchymal cells of limbus (MC-L), which showed phenotypic resemblance to bone marrow derived mesenchymal stem cells (MSC-BM). To understand the transcriptional profile of these cells, microarray experiments were carried out.Methods: RNA was isolated from cultured LEC, MC-L and MSC-BM and microarray experiments were carried out by using Agilent chip (4x44 k). The microarray data was validated by using Realtime and semiquntitative reverse transcription polymerase chain reaction. Results: The microarray analysis revealed specific gene signature of LEC and MC-L, and also their complementary role related to cytokine and growth factor profile, thus supporting the nurturing roles of the MC-L. We have also observed similar and differential gene expression between MC-L and MSC-BM.Conclusions: This study represents the first extensive gene expression analysis of limbal explant culture derived epithelial and mesenchymal cells and as such reveals new insight into the biology, ontogeny, and in vivo function of these cells.

Fourteen gene GBM prognostic signature identifies association of immune response pathway and mesenchymal subtype with high risk group

Background: Recent research on glioblastoma (GBM) has focused on deducing gene signatures predicting prognosis. The present study evaluated the mRNA expression of selected genes and correlated with outcome to arrive at a prognostic gene signature. Methods: Patients with GBM (n = 123) were prospectively recruited, treated with a uniform protocol and followed up. Expression of 175 genes in GBM tissue was determined using qRT-PCR. A supervised principal component analysis followed by derivation of gene signature was performed. Independent validation of the signature was done using TCGA data. Gene Ontology and KEGG pathway analysis was carried out among patients from TCGA cohort. Results: A 14 gene signature was identified that predicted outcome in GBM. A weighted gene (WG) score was found to be an independent predictor of survival in multivariate analysis in the present cohort (HR = 2.507; B = 0.919; p < 0.001) and in TCGA cohort. Risk stratification by standardized WG score classified patients into low and high risk predicting survival both in our cohort (p = <0.001) and TCGA cohort (p = 0.001). Pathway analysis using the most differentially regulated genes (n = 76) between the low and high risk groups revealed association of activated inflammatory/immune response pathways and mesenchymal subtype in the high risk group. Conclusion: We have identified a 14 gene expression signature that can predict survival in GBM patients. A network analysis revealed activation of inflammatory response pathway specifically in high risk group. These findings may have implications in understanding of gliomagenesis, development of targeted therapies and selection of high risk cancer patients for alternate adjuvant therapies.

A 16-Gene Signature Distinguishes Anaplastic Astrocytoma from Glioblastoma

Anaplastic astrocytoma (AA; Grade III) and glioblastoma (GBM; Grade IV) are diffusely infiltrating tumors and are called malignant astrocytomas. The treatment regimen and prognosis are distinctly different between anaplastic astrocytoma and glioblastoma patients. Although histopathology based current grading system is well accepted and largely reproducible, intratumoral histologic variations often lead to difficulties in classification of malignant astrocytoma samples. In order to obtain a more robust molecular classifier, we analysed RT-qPCR expression data of 175 differentially regulated genes across astrocytoma using Prediction Analysis of Microarrays (PAM) and found the most discriminatory 16-gene expression signature for the classification of anaplastic astrocytoma and glioblastoma. The 16-gene signature obtained in the training set was validated in the test set with diagnostic accuracy of 89%. Additionally, validation of the 16-gene signature in multiple independent cohorts revealed that the signature predicted anaplastic astrocytoma and glioblastoma samples with accuracy rates of 99%, 88%, and 92% in TCGA, GSE1993 and GSE4422 datasets, respectively. The protein-protein interaction network and pathway analysis suggested that the 16-genes of the signature identified epithelial-mesenchymal transition (EMT) pathway as the most differentially regulated pathway in glioblastoma compared to anaplastic astrocytoma. In addition to identifying 16 gene classification signature, we also demonstrated that genes involved in epithelial-mesenchymal transition may play an important role in distinguishing glioblastoma from anaplastic astrocytoma.

The potential role of a late gene expression factor, lef2, from Bombyx mori nuclear polyhedrosis virus in very late gene transcription and DNA replication

Several late gene expression factors (Lefs) have been implicated in fostering high levels of transcription from the very late gene promoters of polyhedrin and p10 from baculoviruses. We cloned and characterized from Bombyx mori nuclear polyhedrosis virus a late gene expression factor (Bmlef2) that encodes a 209-amino-acid protein harboring a Cys-rich C-terminal domain. The temporal transcription profiles of lef2 revealed a 1.2-kb transcript in both delayed early and late periods after virus infection. Transcription start site mapping identified the presence of an aphidicolin-sensitive late transcript arising from a TAAG motif located at -352 nucleotides and an aphidicolin-insensitive early transcript originating from a TTGT motif located 35 nucleotides downstream to a TATA box at -312 nucleotides, with respect to the +1 ATG of lef2. BmLef2 trans-activated very late gene expression from both polyhedrin and p10 promoters in transient expression assays. Internal deletion of the Cys-rich domain from the C-terminal region abolished the transcriptional activation. Inactivation of Lef2 synthesis by antisense lef2 transcripts drastically reduced the very late gene transcription but showed little effect on the expression from immediate early promoter. Decrease in viral DNA synthesis and a reduction in virus titer were observed only when antisense lef2 was expressed under the immediate early (ie-1) promoter. Furthermore, the antisense experiments suggested that lef2 plays a direct role in very late gene transcription.

Role of TATATAA element in the regulation of tRNA(1)(Gly) gene expression in Bombyx mori is position dependent

Transcription of tRNA genes by RNA polymerase III is controlled by the internal conserved sequences within the coding region and the immediate upstream flanking sequences. A highly transcribed copy of glycyl tRNA gene tRNA1(Gly)-1 from Bombyx mori is down regulated by sequences located much farther upstream in the region -150 to -300 nucleotides (nt), with respect to the +1 nt of tRNA. The negative regulatory effect has been narrowed down to a sequence motif 'TATATAA', a perfect consensus recognised by the TATA binding protein, TBP. This sequence element, when brought closer to the transcription start point, on the other hand, exerts a positive effect by promoting transcription of the gene devoid of other cis regulatory elements. The identity of the nuclear protein interacting with this 'TATATAA' element to TBP has been established by antibody and mutagenesis studies. The 'TATATAA' element thus influences the transcription of tRNA genes positively or negatively in a position-dependent manner either by recruitment or sequestration of TBP from the transcription machinery.

Genome-Wide Gene Expression Analysis Reveals a Dynamic Interplay between Luteotropic and Luteolytic Factors in the Regulation of Corpus Luteum Function in the Bonnet Monkey (Macaca radiata)

Although LH is essential for survival and function of the corpus luteum (CL) in higher primates, luteolysis occurs during nonfertile cycles without a discernible decrease in circulating LH levels. Using genome-wide expression analysis, several experiments were performed to examine the processes of luteolysis and rescue of luteal function in monkeys. Induced luteolysis with GnRH receptor antagonist (Cetrorelix) resulted in differential regulation of 3949 genes, whereas replacement with exogenous LH (Cetrorelix plus LH) led to regulation of 4434 genes (1563 down-regulation and 2871 up-regulation). A model system for prostaglandin (PG) F-2 alpha-induced luteolysis in the monkey was standardized and demonstrated that PGF(2 alpha) regulated expression of 2290 genes in the CL. Analysis of the LH-regulated luteal transcriptome revealed that 120 genes were regulated in an antagonistic fashion by PGF(2 alpha). Based on the microarray data, 25 genes were selected for validation by real-time RT-PCR analysis, and expression of these genes was also examined in the CL throughout the luteal phase and from monkeys treated with human chorionic gonadotropin (hCG) to mimic early pregnancy. The results indicated changes in expression of genes favorable to PGF(2 alpha) action during the late to very late luteal phase, and expressions of many of these genes were regulated in an opposite manner by exogenous hCG treatment. Collectively, the findings suggest that curtailment of expression of downstream LH-target genes possibly through PGF(2 alpha) action on the CL is among the mechanisms underlying cross talk between the luteotropic and luteolytic signaling pathways that result in the cessation of luteal function, but hCG is likely to abrogate the PGF(2 alpha)-responsive gene expression changes resulting in luteal rescue crucial for the maintenance of early pregnancy. (Endocrinology 150: 1473-1484, 2009)

Standardization and validation of an induced ovulation model system in buffalo cows: Characterization of gene expression changes in the periovulatory follicle

In bovines characterization of biochemical and molecular determinants of the dominant follicle before and during different time intervals after gonadotrophin surge requires precise identification of the dominant follicle from a follicular wave. The objectives of the present study were to standardize an experimental model in buffalo cows for accurately identifying the dominant follicle of the first wave of follicular growth and characterize changes in follicular fluid hormone concentrations as well as expression patterns of various genes associated with the process of ovulation. From the day of estrus (day 0), animals were subjected to blood sampling and ultrasonography for monitoring circulating progesterone levels and follicular growth. On day 7 of the cycle, animals were administered a PGF2α analogue (Tiaprost Trometamol, 750 μg i.m.) followed by an injection of hCG (2000 IU i.m.) 36 h later. Circulating progesterone levels progressively increased from day 1 of the cycle to 2.26 ± 0.17 ng/ml on day 7 of the cycle, but declined significantly after PGF2α injection. A progressive increase in the size of the dominant follicle was observed by ultrasonography. The follicular fluid estradiol and progesterone concentrations in the dominant follicle were 600 ± 16.7 and 38 ± 7.6 ng/ml, respectively, before hCG injection and the concentration of estradiol decreased to 125.8 ± 25.26 ng/ml, but concentration of progesterone increased to 195 ± 24.6 ng/ml, 24 h post-hCG injection. Inh-α and Cyp19A1 expressions in granulosa cells were maximal in the dominant follicle and declined in response to hCG treatment. Progesterone receptor, oxytocin and cycloxygenase-2 expressions in granulosa cells, regarded as markers of ovulation, were maximal at 24 h post-hCG. The expressions of genes belonging to the super family of proteases were also examined; Cathepsin L expression decreased, while ADAMTS 3 and 5 expressions increased 24 h post-hCG treatment. The results of the current study indicate that sequential treatments of PGF2α and hCG during early estrous cycle in the buffalo cow leads to follicular growth that culminates in ovulation. The model system reported in the present study would be valuable for examining temporo-spatial changes in the periovulatory follicle immediately before and after the onset of gonadotrophin surge.

DNA structural variability as a factor in gene expression and evolution

Redundant DNA can buffer sequence dependent structural deviations from an ideal double helix. Buffering serves a mechanistic function by reducing extraneous conformational effects which could interfere with readout or which would impose energetic constraints on evolution. It also serves an evolutionary function by allowing for gradual variations in conformation-dependent regulation of gene expression. Such gradualism is critical for the rate of evolution. The buffer structure concept provides a new interpretation for repetitive DNA and for exons and introns.

Regulation of cytochrome P-450 gene expression. Studies with a cloned probe

A cDNA clone for cytochrome P-450e, a phenobarbitone-inducible species in rat liver, has been isolated and characterized. With the use of this cloned DNA, an attempt has been initiated to elucidate the factors regulating the cytochrome P-450 gene expression. Inhibitors of heme synthesis such as cobalt chloride and 3-amino-1,2,4-triazole block the induction of cytochrome P-450e by phenobarbitone at the level of transcription. This is evident from the decrease in the rate of synthesis of cytochrome P-450e, a decrease in the levels of specific translatable messenger RNA, a decrease in the specific cytoplasmic and nuclear messenger RNA contents, and nuclear transcription of cytochrome P-450e gene, as revealed by hybridization to the cloned probe, under these conditions. It is proposed that heme is a general regulator of cytochrome P-450 gene expression at the level of transcription, whereas the drug or its metabolite would impart the specificity needed for the induction of a particular species.

A comparative study of tbe regulation of cytochrome P-450 and glutathione transferase gene expression in rat liver

A cDNA clone for the Ya subunit of glutathione transferase from rat liver was constructed in E.coli. The clone hybridized to Ya and Yc subunit messenger RNAs. On the basis of experiments involving cell-free translation and hybridization to the cloned probe, it was shown that prototype inducers of cytochrome P-450 such as phenobarbitone and 3-methylcholanthrene as well as inhibitors such as CoCl2 and 3-amino-l,2,4-triazole enhanced the glutathione transferase (Ya+Yc) messenger RNA contents in rat liver. A comparative study with the induction of cytochrome P-450 (b+e) by phenobarbitone revealed that the drug manifested a striking increase in the nuclear pre-messenger RNAs for the cytochrome at 12 hr, but did not significantly affect the same in the case of glutathione transferase (Ya+Yc). 3-Amino-l, 2,4-tnazole and CoCl- blocked the phenobarbitone mediated increase in cytochrome P-450 (b+e) nuclear pre-messenger RNAs. These compounds did not significantly affect the glutathione transferase (Ya+Yc) nuclear pre-messenger RNA levels. The polysomal, poly (A)- containing messenger RNAs for cytochrome P-450 (b+e) increased by 12–15 fold after phenobarbitone administration, reached a maximum around 16hr and then decreased sharply. In comparison, the increase in the case of glutathione transferase (Ya+Yc) mesenger RNAs was sluggish and steady and a value of 3–4 fold was reached around 24 hr. Run-off transcription rates for cytochrome P-450 (b+e) increased by nearly 15 fold in 4 hr after phenobarbitone administration, whereas the increase for glutathione transferase (Ya+Yc) was only 2.0 fold. At 12 hr after the drug administration, the glutathione transferase (Ya+Yc) transcription rates were near normal. Administration of 3-amino-l,2,4-triazole and CoCl2 blocked the phenobarbitone-mediated increase in the transcription of cytochrome P-450 (b+e) messenger RNAs. These compounds at best had only marginal effects on the transcription of glutathione transferase (Ya+Yc) messenger RNAs. The half-life of cytochrome P-450 (b+e) messenger RNA was estimated to be 3–4 hr, whereas that for glutathione transferase (Ya+Yc) was found to be 8-9 hr. Administration of phenobarbitone enhanced the half-life of glutathione transferase (Ya+Yc) messenger RNA by nearly two fold. It is suggested that while transcription activation may play a primary role in the induction of cytochrome P-450 (b+e), the induction of glutathione transferase (Ya+Yc) may essentially involve stabilization of the messenger RNAs.

Positional Preferences of Polypurine/Polypyrimidine Tracts in Saccharomyces cerevisiae Genome: Implications for cis Regulation of Gene Expression

The complete genome of the baker's yeast S. cerevisiae was analyzed for the presence of polypurine/polypyrimidine (poly[pu/py]) repeats and their occurrences were classified on the basis of their location within and outside open reading frames (ORFs). The analysis reveals that such sequence motifs are present abundantly both in coding as well as noncoding regions. Clear positional preferences are seen when these tracts occur in noncoding regions. These motifs appear to occur predominantly at a unit nucleosomal length both upstream and downstream of ORFs. Moreover, there is a biased distribution of polypurines in the coding strands when these motifs occur within open reading frames. The significance of the biased distribution is discussed with reference to the occurrence of these motifs in other known mRNA sequences and expressed sequence tags. A model for cis regulation of gene expression is proposed based on the ability of these motifs to form an intermolecular triple helix structure when present within the coding region and/or to modulate nucleosome positioning via enhanced histone affinity when present outside coding regions.

Distress call-induced gene expression in the brain of the Indian short-nosed fruit bat, Cynopterus sphinx

Individuals in distress emit audible vocalizations to either warn or inform conspecifics. The Indian short-nosed fruit bat, Cynopterus sphinx, emits distress calls soon after becoming entangled in mist nets, which appear to attract conspecifics. Phase I of these distress calls is longer and louder, and includes a secondary peak, compared to phase II. Activity-dependent expression of egr-1 was examined in free-ranging C. sphinx following the emissions and responses to a distress call. We found that the level of expression of egr-1 was higher in bats that emitted a distress call, in adults that responded, and in pups than in silent bats. Up-regulated cDNA was amplified to identify the target gene (TOE1) of the protein Egr-1. The observed expression pattern Toe1 was similar to that of egr-1. These findings suggest that the neuronal activity related to recognition of a distress call and an auditory feedback mechanism induces the expression of Egr-1. Co-expression of egr-1 with Toe1 may play a role in initial triggering of the genetic mechanism that could be involved in the consolidation or stabilization of distress call memories.

Gene expression profiling during Forskolin induced differentiation of BeWo cells by differential display RT-PCR

The differentiation of cytotrophoblasts into syncytiotrophoblasts in the placenta has been employed as a model to investigate stage specific expression as well as regulation of genes during this process. While the cytotrophoblasts are highly invasive and proliferative with relatively less capacity to synthesize pregnancy related proteins, the multinucleated syncytiotrophoblasts are non-proliferative and non-invasive. However, syncytiotrophoblasts are the site of synthesis of a variety of protein, peptide and steroid hormones as well as several growth factors. Both the freshly isolated cytotrophoblasts from human placenta as well as the BeWo cell, a choriocarcinoma cell line model which retain several characteristic of cytotrophoblasts has been employed by us to study regulation of differentiation. In the present study, we have employed the differential display RT-PCR analysis (DD-RT-PCR) to evaluate gene expression changes during Forskolin induced in vitro differentiation of BeWo cells. We have identified several genes which are differentially expressed during differentiation and the differential expression of 10 transcripts was confirmed by Northern blot analysis. Based on the identity of the transcripts an attempt has been made to relate the known function of the gene products, to changes observed during differentiation. Of the several transcripts, one of the transcripts, namely Secretory Leukocyte Protease Inhibitor (SLPI) which is known to have multiple functions was found to increase 15-fold in the syntiotrophoblast.

Stochastic simulations of the origins and implications of long-tailed distributions in gene expression

Gene expression noise results in protein number distributions ranging from long-tailed to Gaussian. We show how long-tailed distributions arise from a stochastic model of the constituent chemical reactions and suggest that, in conjunction with cooperative switches, they lead to more sensitive selection of a subpopulation of cells with high protein number than is possible with Gaussian distributions. Single-cell-tracking experiments are presented to validate some of the assumptions of the stochastic simulations. We also examine the effect of DNA looping on the shape of protein distributions. We further show that when switches are incorporated in the regulation of a gene via a feedback loop, the distributions can become bimodal. This might explain the bimodal distribution of certain morphogens during early embryogenesis.

Regulation of cytochrome P-450b/e gene expression by a heme- and phenobarbitone-modulated transcription factor

The cloned DNA fragment of the cytochrome P-450b/e gene containing the upstream region from position -179 through part of the first exon is faithfully transcribed in freeze-thawed rat liver nuclei. Phenobarbitone treatment of the animal strikingly increases this transcription, and the increase is blocked by cycloheximide (protein synthesis inhibitor) or CoCl2 (heme biosynthetic inhibitor) treatment of animals. This picture correlates very well with the reported cytochrome P-450b/e mRNA levels in vivo and run-on transcription rates in vitro under these conditions. The upstream region (from position -179) was assessed for protein binding with nuclear extracts by nitrocellulose filter binding, gel retardation, DNase I treatment ("footprinting"), and Western blot analysis. Phenobarbitone treatment dramatically increases protein binding to the upstream region, an increase once again blocked by cycloheximide or CoCl2 treatments. Addition of heme in vitro to heme-deficient nuclei and nuclear extracts restores the induced levels of transcription and protein binding to the upstream fragment, respectively. Thus, drug-mediated synthesis and heme-modulated binding of a transcription factor(s) appear involved in the transcriptional activation of the cytochrome P-450b/e genes, and an 85-kDa protein may be a major factor in this regard.