Phase-specific Expression Research Articles

Division of developmental phases of freshwater leech Whitmania pigra and key genes related to neurogenesis revealed by whole genome and transcriptome analysis

The freshwater leech Whitmania pigra (W. pigra) Whitman (Annelida phylum) is a model organism for neurodevelopmental studies. However, molecular biology research on its embryonic development is still scarce. Here, we described a series of developmental stages of the W. pigra embryos and defined five broad stages of embryogenesis: cleavage stages, blastocyst stage, gastrula stage, organogenesis and refinement, juvenile. We obtained a total of 239.64 Gb transcriptome data of eight representative developmental phases of embryos (from blastocyst stage to maturity), which was then assembled into 21,482 unigenes according to our reference genome sequenced by single-molecule real-time (SMRT) long-read sequencing. We found 3114 genes differentially expressed during the eight phases with phase-specific expression pattern. Using a comprehensive transcriptome dataset, we demonstrated that 57, 49 and 77 DEGs were respectively related to morphogenesis, signal pathways and neurogenesis. 49 DEGs related to signal pathways included 30 wnt genes, 14 notch genes, and 5 hedgehog genes. In particular, we found a cluster consisting of 7 genes related to signal pathways as well as synapses, which were essential for regulating embryonic development. Eight genes cooperatively participated in regulating neurogenesis. Our results reveal the whole picture of W. pigra development mechanism from the perspective of transcriptome and provide new clues for organogenesis and neurodevelopmental studies of Annelida species.

Histochemical analysis of the biphasic properties of formalin pain-induced behavior

The formalin test has been established as a method for evaluating mouse models of pain. Although there have been numerous reports of formalin-pain-induced behavior, few reports of a detailed histochemical analysis of the central nervous system focus on behavioral biphasic properties.To investigate the alternation of spinal neuronal activity with formalin-induced pain, we performed immunofluorescent staining with c-Fos antibodies as neuronal activity markers using acute pain model mice induced by 2% formalin stimulation. As a result, phase-specific expression patterns were observed. In the spinal dorsal horn region, there were many neural activities in the deep region (layers V−VII) in the behavioral first phase and those in the surface region (layers I−III) in the behavioral second phase. Furthermore, we conducted comparative studies using low concentrations (0.25%) of formalin and capsaicin, which did not show distinct behavioral biphasic properties. Neural activity was observed only in the spinal dorsal horn surface region for both stimuli.Our study suggested that the histochemical biphasic nature of formalin-induced pain was attributable to the activity of the deep region of the spinal cord. In the future, treatment strategies focusing on the deep region neuron will lead to the development of effective treatments for allodynia and intractable chronic pain.

Phase-specific transcriptional patterns of the oomycete pathogen Phytophthora sojae unravel genes essential for asexual development and pathogenic processes.

Oomycetes are filamentous microorganisms easily mistaken as fungi but vastly differ in physiology, biochemistry, and genetics. This commonly-held misconception lead to a reduced effectiveness by using conventional fungicides to control oomycetes, thus it demands the identification of novel functional genes as target for precisely design oomycetes-specific microbicide. The present study initially analyzed the available transcriptome data of the model oomycete pathogen, Phytophthora sojae, and constructed an expression matrix of 10,953 genes across the stages of asexual development and host infection. Hierarchical clustering, specificity, and diversity analyses revealed a more pronounced transcriptional plasticity during the stages of asexual development than that in host infection, which drew our attention by particularly focusing on transcripts in asexual development stage to eventually clustered them into 6 phase-specific expression modules. Three of which respectively possessing a serine/threonine phosphatase (PP2C) expressed during the mycelial and sporangium stages, a histidine kinase (HK) expressed during the zoospore and cyst stages, and a bZIP transcription factor (bZIP32) exclusive to the cyst germination stage were selected for down-stream functional validation. In this way, we demonstrated that PP2C, HK, and bZIP32 play significant roles in P. sojae asexual development and virulence. Thus, these findings provide a foundation for further gene functional annotation in oomycetes and crop disease management.

Exploring key developmental phases and phase-specific genes across the entirety of anther development in maize.

Anther development from stamen primordium to pollen dispersal is complex and essential to sexual reproduction. How this highly dynamic and complex developmental process is controlled genetically is not well understood, especially for genes involved in specific key developmental phases. Here we generated RNA sequencing libraries spanning 10 key stages across the entirety of anther development in maize (Zea mays). Global transcriptome analyses revealed distinct phases of cell division and expansion, meiosis, pollen maturation, and mature pollen, for which we detected 50, 245, 42, and 414 phase-specific marker genes, respectively. Phase-specific transcription factor genes were significantly enriched in the phase of meiosis. The phase-specific expression of these marker genes was highly conserved among the maize lines Chang7-2 and W23, indicating they might have important roles in anther development. We explored a desiccation-related protein gene, ZmDRP1, which was exclusively expressed in the tapetum from the tetrad to the uninucleate microspore stage, by generating knockout mutants. Notably, mutants in ZmDRP1 were completely male-sterile, with abnormal Ubisch bodies and defective pollen exine. Our work provides a glimpse into the gene expression dynamics and a valuable resource for exploring the roles of key phase-specific genes that regulate anther development.

MiR-15a and miR-15b modulate natural killer and CD8+T-cell activation and anti-tumor immune response by targeting PD-L1 in neuroblastoma

Neuroblastoma (NB) is an enigmatic and deadliest pediatric cancer to treat. The major obstacles to the effective immunotherapy treatments in NB are defective immune cells and the immune evasion tactics deployed by the tumor cells and the stromal microenvironment. Nervous system development during embryonic and pediatric stages is critically mediated by non-coding RNAs such as micro RNAs (miR). Hence, we explored the role of miRs in anti-tumor immune response via a range of data-driven workflows and in vitro & in vivo experiments. Using the TARGET, NB patient dataset (n=249), we applied the robust bioinformatic workflows incorporating differential expression, co-expression, survival, heatmaps, and box plots. We initially demonstrated the role of miR-15a-5p (miR-15a) and miR-15b-5p (miR-15b) as tumor suppressors, followed by their negative association with stromal cell percentages and a statistically significant negative regulation of T and natural killer (NK) cell signature genes, especially CD274 (PD-L1) in stromal-low patient subsets. The NB phase-specific expression of the miR-15a/miR-15b-PD-L1 axis was further corroborated using the PDX (n=24) dataset. We demonstrated miR-15a/miR-15b mediated degradation of PD-L1 mRNA through its interaction with the 3'-untranslated region and the RNA-induced silencing complex using sequence-specific luciferase activity and Ago2 RNA immunoprecipitation assays. In addition, we established miR-15a/miR-15b induced CD8+T and NK cell activation and cytotoxicity against NB in vitro. Moreover, injection of murine cells expressing miR-15a reduced tumor size, tumor vasculature and enhanced the activation and infiltration of CD8+T and NK cells into the tumors in vivo. We further established that blocking the surface PD-L1 using an anti-PD-L1 antibody rescued miR-15a/miR-15b induced CD8+T and NK cell-mediated anti-tumor responses. These findings demonstrate that miR-15a and miR-15b induce an anti-tumor immune response by targeting PD-L1 in NB.

Comparative Gene Expression Analysis Reveals Similarities and Differences of Chronic Myeloid Leukemia Phases.

Simple SummaryChronic myeloid leukemia (CML) is a blood cancer with a very good long-term prognosis for the majority of patients. Still, some patients relapse or progress under the standard therapy. Therefore, we performed an in-depth computational gene expression analysis to determine similarities and differences between the CML phases at the level of single genes, signaling pathways and gene regulatory networks and further compared treatment-resistant patients to the individual phases to identify associated expression differences. Our study provides several lines of evidence that CML development represents a three rather than a two step process. The identified characteristic gene expression alterations in advanced phases indicate that affected patients could potentially profit from other existing drugs. Moreover, characteristic signaling pathway changes identified for patients for which the standard therapy was inefficient may allow to predict such patients before treatment start and could also provide a basis to develop treatment strategies for them.Chronic myeloid leukemia (CML) is a slowly progressing blood cancer that primarily affects elderly people. Without successful treatment, CML progressively develops from the chronic phase through the accelerated phase to the blast crisis, and ultimately to death. Nowadays, the availability of targeted tyrosine kinase inhibitor (TKI) therapies has led to long-term disease control for the vast majority of patients. Nevertheless, there are still patients that do not respond well enough to TKI therapies and available targeted therapies are also less efficient for patients in accelerated phase or blast crises. Thus, a more detailed characterization of molecular alterations that distinguish the different CML phases is still very important. We performed an in-depth bioinformatics analysis of publicly available gene expression profiles of the three CML phases. Pairwise comparisons revealed many differentially expressed genes that formed a characteristic gene expression signature, which clearly distinguished the three CML phases. Signaling pathway expression patterns were very similar between the three phases but differed strongly in the number of affected genes, which increased with the phase. Still, significant alterations of MAPK, VEGF, PI3K-Akt, adherens junction and cytokine receptor interaction signaling distinguished specific phases. Our study also suggests that one can consider the phase-wise CML development as a three rather than a two-step process. This is in accordance with the phase-specific expression behavior of 24 potential major regulators that we predicted by a network-based approach. Several of these genes are known to be involved in the accumulation of additional mutations, alterations of immune responses, deregulation of signaling pathways or may have an impact on treatment response and survival. Importantly, some of these genes have already been reported in relation to CML (e.g., AURKB, AZU1, HLA-B, HLA-DMB, PF4) and others have been found to play important roles in different leukemias (e.g., CDCA3, RPL18A, PRG3, TLX3). In addition, increased expression of BCL2 in the accelerated and blast phase indicates that venetoclax could be a potential treatment option. Moreover, a characteristic signaling pathway signature with increased expression of cytokine and ECM receptor interaction pathway genes distinguished imatinib-resistant patients from each individual CML phase. Overall, our comparative analysis contributes to an in-depth molecular characterization of similarities and differences of the CML phases and provides hints for the identification of patients that may not profit from an imatinib therapy, which could support the development of additional treatment strategies.

Coordination mechanism of cell and cyanelle division in the glaucophyte alga Cyanophora sudae

In unicellular algae with a single chloroplast, two mechanisms coordinate cell and chloroplast division: the S phase-specific expression of chloroplast division genes and the permission of cell cycle progression from prophase to metaphase by the onset of chloroplast division. This study investigated whether a similar mechanism exists in a unicellular alga with multiple chloroplasts using the glaucophyte alga Cyanophora sudae, which contains four chloroplasts (cyanelles). Cells with eight cyanelles appeared after the S phase arrest with a topoisomerase inhibitor camptothecin, suggesting that the mechanism of S phase-specific expression of cyanelle division genes was conserved in this alga. Inhibition of peptidoglycan synthesis by β-lactam antibiotic ampicillin arrested cells in the S-G2 phase, and inhibition of septum invagination with cephalexin resulted in cells with two nuclei and one cyanelle, despite inhibition of cyanelle division. This indicates that even in the unicellular alga with four chloroplasts, the cell cycle progresses to the M phase following the progression of chloroplast division to a certain division stage. These results suggested that C. sudae has two mechanisms for coordinating cell and cyanelle division, similar to the unicellular algae with a single chloroplast.

Transcriptomic phases of periodontitis lesions using the nonhuman primate model

We used a nonhuman primate model of ligature-induced periodontitis to identify patterns of gingival transcriptomic after changes demarcating phases of periodontitis lesions (initiation, progression, resolution). A total of 18 adult Macaca mulatta (12–22 years) had ligatures placed (premolar, 1st molar teeth) in all 4 quadrants. Gingival tissue samples were obtained (baseline, 2 weeks, 1 and 3 months during periodontitis and at 5 months resolution). Gene expression was analyzed by microarray [Rhesus Gene 1.0 ST Array (Affymetrix)]. Compared to baseline, a large array of genes were significantly altered at initiation (n = 6049), early progression (n = 4893), and late progression (n = 5078) of disease, with the preponderance being up-regulated. Additionally, 1918 genes were altered in expression with disease resolution, skewed towards down-regulation. Assessment of the genes demonstrated specific profiles of epithelial, bone/connective tissue, apoptosis/autophagy, metabolism, regulatory, immune, and inflammatory responses that were related to health, stages of disease, and tissues with resolved lesions. Unique transcriptomic profiles occured during the kinetics of the periodontitis lesion exacerbation and remission. We delineated phase specific gene expression profiles of the disease lesion. Detection of these gene products in gingival crevicular fluid samples from human disease may contribute to a better understanding of the biological dynamics of the disease to improve patient management.

De novo transcriptome assembly from the gonads of a scleractinian coral, Euphyllia ancora: molecular mechanisms underlying scleractinian gametogenesis

BackgroundSexual reproduction of scleractinians has captured the attention of researchers and the general public for decades. Although extensive ecological data has been acquired, underlying molecular and cellular mechanisms remain largely unknown. In this study, to better understand mechanisms underlying gametogenesis, we isolated ovaries and testes at different developmental phases from a gonochoric coral, Euphyllia ancora, and adopted a transcriptomic approach to reveal sex- and phase-specific gene expression profiles. In particular, we explored genes associated with oocyte development and maturation, spermiogenesis, sperm motility / capacitation, and fertilization.Results1.6 billion raw reads were obtained from 24 gonadal samples. De novo assembly of trimmed reads, and elimination of contigs derived from symbiotic dinoflagellates (Symbiodiniaceae) and other organisms yielded a reference E. ancora gonadal transcriptome of 35,802 contigs. Analysis of 4 developmental phases identified 2023 genes that were differentially expressed during oogenesis and 678 during spermatogenesis. In premature/mature ovaries, 631 genes were specifically upregulated, with 538 in mature testes. Upregulated genes included those involved in gametogenesis, gamete maturation, sperm motility / capacitation, and fertilization in other metazoans, including humans. Meanwhile, a large number of genes without homology to sequences in the SWISS-PROT database were also observed among upregulated genes in premature / mature ovaries and mature testes.ConclusionsOur findings show that scleractinian gametogenesis shares many molecular characteristics with that of other metazoans, but it also possesses unique characteristics developed during cnidarian and/or scleractinian evolution. To the best of our knowledge, this study is the first to create a gonadal transcriptome assembly from any scleractinian. This study and associated datasets provide a foundation for future studies regarding gametogenesis and differences between male and female colonies from molecular and cellular perspectives. Furthermore, our transcriptome assembly will be a useful reference for future development of sex-specific and/or stage-specific germ cell markers that can be used in coral aquaculture and ecological studies.

Improved recovery of cell-cycle gene expression in Saccharomyces cerevisiae from regulatory interactions in multiple omics data

BackgroundGene expression is regulated by DNA-binding transcription factors (TFs). Together with their target genes, these factors and their interactions collectively form a gene regulatory network (GRN), which is responsible for producing patterns of transcription, including cyclical processes such as genome replication and cell division. However, identifying how this network regulates the timing of these patterns, including important interactions and regulatory motifs, remains a challenging task.ResultsWe employed four in vivo and in vitro regulatory data sets to investigate the regulatory basis of expression timing and phase-specific patterns cell-cycle expression in Saccharomyces cerevisiae. Specifically, we considered interactions based on direct binding between TF and target gene, indirect effects of TF deletion on gene expression, and computational inference. We found that the source of regulatory information significantly impacts the accuracy and completeness of recovering known cell-cycle expressed genes. The best approach involved combining TF-target and TF-TF interactions features from multiple datasets in a single model. In addition, TFs important to multiple phases of cell-cycle expression also have the greatest impact on individual phases. Important TFs regulating a cell-cycle phase also tend to form modules in the GRN, including two sub-modules composed entirely of unannotated cell-cycle regulators (STE12-TEC1 and RAP1-HAP1-MSN4).ConclusionOur findings illustrate the importance of integrating both multiple omics data and regulatory motifs in order to understand the significance regulatory interactions involved in timing gene expression. This integrated approached allowed us to recover both known cell-cycles interactions and the overall pattern of phase-specific expression across the cell-cycle better than any single data set. Likewise, by looking at regulatory motifs in the form of TF-TF interactions, we identified sets of TFs whose co-regulation of target genes was important for cell-cycle expression, even when regulation by individual TFs was not. Overall, this demonstrates the power of integrating multiple data sets and models of interaction in order to understand the regulatory basis of established biological processes and their associated gene regulatory networks.

Promoter-dependent nuclear RNA degradation ensures cell cycle-specific gene expression

Cell cycle progression depends on phase-specific gene expression. Here we show that the nuclear RNA degradation machinery plays a lead role in promoting cell cycle-dependent gene expression by triggering promoter-dependent co-transcriptional RNA degradation. Single molecule quantification of RNA abundance in different phases of the cell cycle indicates that relative curtailment of gene expression in certain phases is attained even when transcription is not completely inhibited. When nuclear ribonucleases are deleted, transcription of the Saccharomyces cerevisiae G1-specific axial budding gene AXL2 is detected throughout the cell cycle and its phase-specific expression is lost. Promoter replacement abolished cell cycle-dependent RNA degradation and rendered the RNA insensitive to the deletion of nuclear ribonucleases. Together the data reveal a model of gene regulation whereby RNA abundance is controlled by promoter-dependent induction of RNA degradation.

Overlapping Patterns of Gene Expression Between Gametophyte and Sporophyte Phases in the Fern Polypodium amorphum (Polypodiales).

Ferns are unique among land plants in having sporophyte and gametophyte phases that are both free living and fully independent. Here, we examine patterns of sporophytic and gametophytic gene expression in the fern Polypodium amorphum, a member of the homosporous polypod lineage that comprises 80% of extant fern diversity, to assess how expression of a common genome is partitioned between two morphologically, ecologically, and nutritionally independent phases. Using RNA-sequencing, we generated transcriptome profiles for three replicates of paired samples of sporophyte leaf tissue and whole gametophytes to identify genes with significant differences in expression between the two phases. We found a nearly 90% overlap in the identity and expression levels of the genes expressed in both sporophytes and gametophytes, with less than 3% of genes uniquely expressed in either phase. We compare our results to those from similar studies to establish how phase-specific gene expression varies among major land plant lineages. Notably, despite having greater similarity in the identity of gene families shared between P. amorphum and angiosperms, P. amorphum has phase-specific gene expression profiles that are more like bryophytes and lycophytes than seed plants. Our findings suggest that shared patterns of phase-specific gene expression among seed-free plants likely reflect having relatively large, photosynthetic gametophytes (compared to the gametophytes of seed plants that are highly reduced). Phylogenetic analyses were used to further investigate the evolution of phase-specific expression for the phototropin, terpene synthase, and MADS-box gene families.

Identification of candidate pathogenicity determinants of Rhizoctonia solani AG1-IA, which causes sheath blight disease in rice.

Sheath blight disease is one of the predominant diseases of rice and it is caused by the necrotrophic fungal pathogen Rhizoctonia solani. The mechanistic insight about its widespread success as a broad host range pathogen is limited. In this study, we endeavor to identify pathogenicity determinants of R. solani during infection process in rice. Through RNAseq analysis, we identified a total of 65 and 232 R. solani (strain BRS1) genes to be commonly upregulated in three different rice genotypes (PB1, Tetep, and TP309) at establishment and necrotrophic phase, respectively. The induction of genes encoding extracellular protease, ABC transporter, and transcription factors were notable during establishment phase. While during necrotrophic phase, several CAZymes, sugar transporters, cellular metabolism, and protein degradation-related genes were prominently induced. We have also identified few putative secreted effector encoding genes that were upregulated during pathogenesis. The qPCR analysis further validated the phase-specific expression dynamics of some selected putative effectors and pathogenicity-associated genes. Overall, the present study reports identification of key genes and processes that might be crucial for R. solani pathogenesis. The ability to effectively damage host cell wall and survive in hostile plant environment by managing oxidative stress, cytotoxic compounds, etc. is being proposed to be important for pathogenesis of R. solani in rice. The functional characterization of these genes would provide key insights about this important pathosystem and facilitate development of strategies to control this devastating disease.

Transcriptome and metabolite profiling of the infection cycle of Zymoseptoria tritici on wheat reveals a biphasic interaction with plant immunity involving differential pathogen chromosomal contributions and a variation on the hemibiotrophic lifestyle definition.

The hemibiotrophic fungus Zymoseptoria tritici causes Septoria tritici blotch disease of wheat (Triticum aestivum). Pathogen reproduction on wheat occurs without cell penetration, suggesting that dynamic and intimate intercellular communication occurs between fungus and plant throughout the disease cycle. We used deep RNA sequencing and metabolomics to investigate the physiology of plant and pathogen throughout an asexual reproductive cycle of Z. tritici on wheat leaves. Over 3,000 pathogen genes, more than 7,000 wheat genes, and more than 300 metabolites were differentially regulated. Intriguingly, individual fungal chromosomes contributed unequally to the overall gene expression changes. Early transcriptional down-regulation of putative host defense genes was detected in inoculated leaves. There was little evidence for fungal nutrient acquisition from the plant throughout symptomless colonization by Z. tritici, which may instead be utilizing lipid and fatty acid stores for growth. However, the fungus then subsequently manipulated specific plant carbohydrates, including fructan metabolites, during the switch to necrotrophic growth and reproduction. This switch coincided with increased expression of jasmonic acid biosynthesis genes and large-scale activation of other plant defense responses. Fungal genes encoding putative secondary metabolite clusters and secreted effector proteins were identified with distinct infection phase-specific expression patterns, although functional analysis suggested that many have overlapping/redundant functions in virulence. The pathogenic lifestyle of Z. tritici on wheat revealed through this study, involving initial defense suppression by a slow-growing extracellular and nutritionally limited pathogen followed by defense (hyper) activation during reproduction, reveals a subtle modification of the conceptual definition of hemibiotrophic plant infection.

Abstract 3491: Rapid and ultra-sensitive single-cell transcript profiling with droplet digital PCR (ddPCR): Application to cell cycle analysis

Abstract The recognition that tumors and other apparently homogeneous cell populations are often heterogeneous and display stochastic expression of their genes, has motivated the development of methods to profile gene expression in single cells. Among these, RNA-Seq is the most powerful for broad-scale profiling and discovery of genes and pathways influencing biological processes and states of interest. However, given the extent of sample manipulation required to prepare each RNA-Seq library – and the trade-off between sensitivity and sample throughput – it is important to validate discovery findings by orthogonal digital methods with comparable or greater sensitivity and accuracy. Droplet Digital PCR is such a method that can rapidly and cost-effectively replicate (or challenge) RNA-Seq results with minimal sample processing when defined targets are to be verified (Chen et al, 2012). Furthermore, ddPCR is performed in 96-well plates and is well-suited to high throughput studies of focused sets of genes in large numbers of single cells (e.g. hundreds in a day). Here we present a simple and robust workflow for profiling multiplexed, transcript targets in flow-sorted, Jurkat single-cells using Taqman 5’ probe hydrolysis assays. We demonstrate that Bio-Rad's QX100 & QX200 Droplet Digital PCR systems provide absolute counts of transcripts from >100,000 copies to <10 copies per cell, using unpurified and unamplified cDNA synthesized in 30’ in the original wells used for flow-sorting. The results are highly reproducible and the method is currently being used to explore the contribution of the cell-cycle to expression differences among single-cell populations. More than 30 cell-cycle genes have been evaluated by ddPCR in exponentially growing Jurkat cells. Several genes, including CCNB1, CDK1 and DNM2, show an unusually broad range of transcript levels (over 2 logs), suggestive of their known phase-specific cell-cycle roles. Jurkat cells have been FACS sorted for either G0/G1 (“G1”) or G2/M (“G2”) single-cells using DyeCycle Orange to stain for DNA content, and are being assessed for phase-specific expression profiles. Preliminary results indicate that G2 single-cells are enriched for higher CCNB1 and CDK1 transcript levels as expected from their role in promoting early events in mitosis; however, even G2 cells still show expression levels spanning ∼2 logs, suggesting cell-cycle phase heterogeneity and/or bursty transcription within the sorted G2 population. While the average mRNA level of even the house-keeping gene, GAPDH, is several-fold higher in G2 cells as compared to G1 cells, transcripts for CCNB1 and CDK1 average ∼6-fold higher in G2 than in G1. In this work, we demonstrate both the analytical performance and the biological utility of droplet digital PCR for single-cell analysis. This same workflow could be used for fingerprinting and study of sub-populations in tumors. Citation Format: George Karlin-Neumann, Shenglong Wang, Camille Troup, Yann Jouvenot, Eli Hefner. Rapid and ultra-sensitive single-cell transcript profiling with droplet digital PCR (ddPCR): Application to cell cycle analysis. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3491. doi:10.1158/1538-7445.AM2014-3491

A Provisional Gene Regulatory Atlas for Mouse Heart Development

Congenital Heart Disease (CHD) is one of the most common birth defects. Elucidating the molecular mechanisms underlying normal cardiac development is an important step towards early identification of abnormalities during the developmental program and towards the creation of early intervention strategies. We developed a novel computational strategy for leveraging high-content data sets, including a large selection of microarray data associated with mouse cardiac development, mouse genome sequence, ChIP-seq data of selected mouse transcription factors and Y2H data of mouse protein-protein interactions, to infer the active transcriptional regulatory network of mouse cardiac development. We identified phase-specific expression activity for 765 overlapping gene co-expression modules that were defined for obtained cardiac lineage microarray data. For each co-expression module, we identified the phase of cardiac development where gene expression for that module was higher than other phases. Co-expression modules were found to be consistent with biological pathway knowledge in Wikipathways, and met expectations for enrichment of pathways involved in heart lineage development. Over 359,000 transcription factor-target relationships were inferred by analyzing the promoter sequences within each gene module for overrepresentation against the JASPAR database of Transcription Factor Binding Site (TFBS) motifs. The provisional regulatory network will provide a framework of studying the genetic basis of CHD.

Developing a cyclin blueprint as a tool for mapping the cell cycle in GS-NS0

The cell cycle is at the center of growth, productivity, and death of mammalian cell cultures. There exists a need to identify and quantify major landmarks in the cell cycle of industrially relevant mammalian cell lines and its association with productivity; central for designing productivity optimization strategies. Herein, we studied the expression of three cyclins, under both perturbed and unperturbed growth, by flow cytometry in batch cultures of GS-NS0. The perturbed systems involved two different DNA synthesis inhibitors, thymidine and dimethyl sulfoxide (DMSO). This approach enables the establishment of characteristic cyclin profiles, timings, and thresholds. In particular, two G1 class cyclins (D1 and E1), and one G2 cyclin (B1) were investigated. Cyclin B1 showed a clear cell cycle phase-specific expression increasing during G2 phase where it was approximately 40% higher when compared to G1 phase. Similarly, cyclin E1 showed a clear pattern being expressed approximately 10% higher in G1 compared to G2 phase and decreased through S phase. Cyclin D1 expression was fairly invariable throughout the cell cycle phases. The observed patterns provide a blueprint of the cell line's cell cycle, which can be used for the development of biologically accurate and experimentally validated distributed cell cycle models.

Regulation of bcl-2 transcription by estrogen receptor-α and c-Jun in human endometrium

The estrogen-estrogen receptor (ER) signaling pathway plays crucial physiologic roles in not only the control of reproduction, but also in the generation of cancer in the breast and uterus. While some ER target genes have been identified containing the estrogen-responsive element (ERE), others are activated eventually by ER via protein-protein interaction without binding to ERE. In a previous study, we identified that the proliferative phase-specific expression of the bcl-2 gene in glandular cells could be regulated by the binding of c-Jun to its motifs in the promoter. Results from our present study indicate that the menstrual cyclic expression of bcl-2 could be controlled by either direct binding of ERα to ERE in the c-Jun promoter or the interaction of ERα with c-Jun that binds to its motifs in the bcl-2 gene. Intriguingly, the transcriptionally active form of c-Jun phosphorylated at Ser63 was identified binding to its motifs in the bcl-2 gene in a menstrual cyclic non-specific manner. Our study revealed a novel mechanism that transcriptionally regulates the expression of bcl-2 in the normal human endometrium.

Mice with global deletion of mitogen inducible gene 6 display rapid and severe cartilage and subchondral bone damage after ligament and meniscus injury

s / Osteoarthritis and Cartilage 21 (2013) S9–S62 S13 protein subunits was investigated in growth plates of 6 week old mice by immunohistochemistry. Expression of RMRP, its protein subunits and important chondrogenic differentiation markers (Col2a1, Acan, Col10a1, Sox9, Runx2) was also determined by RT-qPCR and immunoblotting in chondrogenically differentiating ATDC5 cells and human bone marrow stem cells (hBMSCs). In addition, RMRP was functionally investigated during chondrogenic differentiation by targeting RMRP expression using a specific siRNA duplex and analysing the outcome on chondrogenic gene expression. Furthermore, chondrogenic induction of dermal fibroblasts was used as a model to test the effect of CHH-associated mutations in RMRP on chondrogenic differentiation of CHH patient-derived fibroblasts. Results: To determine expression of RNase MRP during chondrogenic differentiation in vivo, mouse growth plate sections were immunohistochemically stained for RNaseMRP protein subunits.Without exception, all proteins showed identical growth plate distribution patterns: Resting zone chondrocytes express RNaseMRPproteins,whereas nooronlyweak expression levelswere observed inproliferative chondrocytes. In contrast, expression of RNase MRP proteins was clearly detectable in hypertrophic chondrocytes in the growth plate. This differentiation phase-specific expression of RNase MRP protein subunits was confirmed for Rpp25, Rpp30 and Rpp40 on ATDC5 and hBMSC differentiation models (mRNA); their expressionwas upregulated from day 14 in differentiation onwards, simultaneously with hypertrophic genes. Knockdown of RMRP expression by RNAi resulted in upregulation of RNase MRP mRNA substrates Clb2 and viperin and decreased chondrogenic marker expression (Sox9, Col2a1, Runx2, Col10a1), suggesting a role for RMRP in chondrogenic differentiation. In addition, knockdown of RMRP also resulted in altered Bapx1/Nkx3.2 expression, further substantiating a role in chondrocyte hypertrophy. To verify our results we transdifferentiated CHH patientderived fibroblasts and matched controls to acquire a chondrocyte-like phenotype. Chondrogenic capacity of the CHH patients was clearly disrupted and showed a more than 50% reduced Col10a1 expression as compared to matched control cultures. Conclusions: Our data show that expression of RNase MRP and its subunits is regulated during chondrogenic differentiation, is predominantly associated with chondrocyte hypertrophy and shows functional cross-talk with chondrogenic pathways. These data are the first to show the association of RNase MRP with skeletal development and may provide an explanation for the CHH-associated skeletal phenotype as well as unveiling a potential novel regulator of chondrocyte hypertrophic differentiation in skeletal development. 10 MICE WITH GLOBAL DELETION OF MITOGEN INDUCIBLE GENE 6 DISPLAY RAPID AND SEVERE CARTILAGE AND SUBCHONDRAL BONE DAMAGE AFTER LIGAMENT AND MENISCUS INJURY D.M. Joiner, K.D. Less, E. VanWieren, B.O. Williams. Van Andel Res. Inst., Grand Rapids, MI, USA Purpose: The Epidermal Growth Factor Receptor (EGFR) Signaling pathway plays important roles in osteochondral metabolism and joint disease, which can be triggered by damage to ligaments andmenisci. The goal of this study was to examine the early effect of global deletion of mitogen inducible gene 6 (MIG6), a negative regulator of EGFR, on cartilage and subchondral bone in mice after ligament and meniscus injury. Methods: Global loss of function MIG6 knock-out mice (MIG6-/-) and wild-type (WT) littermates were anesthetized and ligament transections and meniscus removal were performed on the right knee using an established surgical destabilization method. Subchondral bone from surgery knees and untouched contralateral control knees was assessed with x-ray. HE arrows point to bone erosion). Bone damage in surgery knees was greater for MIG6-/animals compared to the control knee and to surgery knees from WT littermates (Figure 2). Ectopic cartilaginous growth pockets were present inMIG6-/mice 4 weeks after ligament and meniscus damage (Figure 1). MIG6-/mice also had increased EGFR expression 4 weeks after surgery (Figure 3) and increased osteoclasts (Figure 4; arrows point to osteoclasts) in the knee joint compared to WT animals. Conclusions: MIG6-/mice displayed rapid and more severe bone and cartilage damage after knee injury compared to WT mice. MIG6 and EGFR could be potential therapeutic targets for preventing cartilage and bone damage and subsequent joint disease development after ligament and meniscus injury.

Dual Effects of IL-1 Overactivity on the Immune System in a Mouse Model of Arthritis due to Deficiency of IL-1 Receptor Antagonist

Previous studies have revealed the significance of cytokine interleukin 1 (IL-1) in the onset and progression of rheumatoid arthritis (RA). The precise molecular mechanisms related to IL-1 underlying RA is still elusive. We conducted a whole genome-wide transcriptomal comparison of wild-type (WT) and arthritis-prone IL-1 receptor antagonist (IL-1rn) deficient BALB/c mice to address this issue. To refine our search efforts, gene expression profiling was also performed on paired wild-type and arthritis-resistant IL-1rn deficient DBA/1 mice as internal controls when identifying causative arthritis candidate genes. Two hundred and fifteen transcripts were found to be dysregulated greater than or equal to 2-fold in the diseased mice. The altered transcriptome in BALB/c mice revealed increased myeloid cell activities and impaired lymphocyte functionality, suggesting dual regulatory effects of IL-1 hyperactivity on immunological changes associated with arthritis development. Phase-specific gene expression changes were identified, such as early increase and late decrease of heat shock protein coding genes. Moreover, common gene expression changes were also observed, especially the upregulation of paired Ig-like receptor A (Pira) in both early and late phases of arthritis. Real-time PCR was performed to validate the expression of Pira and an intervention experiment with a major histocompatibility complex (MHC) class I inhibitor (brefeldin A) was carried out to investigate the role of suppressing Pira activity. We conclude that global pattern changes of common and distinct gene expressions may represent novel opportunities for better control of RA through early diagnosis and development of alternative therapeutic strategies.