Global comparisons of gene expression profiles between species provide significant insight into gene regulation, evolutionary processes and disease mechanisms. In this work, we describe a flexible and intuitive approach for global expression profiling of closely related species, using high-density exon arrays designed for a single reference genome. The high-density probe coverage of exon arrays allows us to select identical sets of perfect-match probes to measure expression levels of orthologous genes. This eliminates a serious confounding factor in probe affinity effects of species-specific microarray probes, and enables direct comparisons of estimated expression indexes across species. Using a newly designed Affymetrix exon array, with eight probes per exon for approximately 315,000 exons in the human genome, we conducted expression profiling in corresponding tissues from humans, chimpanzees and rhesus macaques. Quantitative real-time PCR analysis of differentially expressed candidate genes is highly concordant with microarray data, yielding a validation rate of 21/22 for human versus chimpanzee differences, and 11/11 for human versus rhesus differences. This method has the potential to greatly facilitate biomedical and evolutionary studies of gene expression in nonhuman primates and can be easily extended to expression array design and comparative analysis of other animals and plants.

In our previous study, it was well defined that IGFBP7 was an important tumor suppressor gene in colorectal cancer (CRC). We aimed to uncover the downstream molecules responsible for IGFBP7's behaviour in this study.

To discover new rheumatoid arthritis (RA) risk loci, we systematically examined 370 SNPs from 179 independent loci with P < 0.001 in a published meta-analysis of RA genome-wide association studies (GWAS) of 3,393 cases and 12,462 controls. We used Gene Relationships Across Implicated Loci (GRAIL), a computational method that applies statistical text mining to PubMed abstracts, to score these 179 loci for functional relationships to genes in 16 established RA disease loci. We identified 22 loci with a significant degree of functional connectivity. We genotyped 22 representative SNPs in an independent set of 7,957 cases and 11,958 matched controls. Three were convincingly validated: CD2-CD58 (rs11586238, P = 1 x 10(-6) replication, P = 1 x 10(-9) overall), CD28 (rs1980422, P = 5 x 10(-6) replication, P = 1 x 10(-9) overall) and PRDM1 (rs548234, P = 1 x 10(-5) replication, P = 2 x 10(-8) overall). An additional four were replicated (P < 0.0023): TAGAP (rs394581, P = 0.0002 replication, P = 4 x 10(-7) overall), PTPRC (rs10919563, P = 0.0003 replication, P = 7 x 10(-7) overall), TRAF6-RAG1 (rs540386, P = 0.0008 replication, P = 4 x 10(-6) overall) and FCGR2A (rs12746613, P = 0.0022 replication, P = 2 x 10(-5) overall). Many of these loci are also associated to other immunologic diseases.

Glycine receptors (GlyRs) play important roles in regulating hippocampal neural network activity and spinal nociception. Here we show that, in cultured rat hippocampal (HIP) and spinal dorsal horn (SDH) neurons, 17-beta-estradiol (E2) rapidly and reversibly reduced the peak amplitude of whole-cell glycine-activated currents (IGly). In outside-out membrane patches from HIP neurons devoid of nuclei, E2 similarly inhibited IGly, suggesting a non-genomic characteristic. Moreover, the E2 effect on IGly persisted in the presence of the calcium chelator BAPTA, the protein kinase inhibitor staurosporine, the classical ER (i.e. ERalpha and ERbeta) antagonist tamoxifen, or the G-protein modulators, favoring a direct action of E2 on GlyRs. In HEK293 cells expressing various combinations of GlyR subunits, E2 only affected the IGly in cells expressing alpha2, alpha2beta or alpha3beta subunits, suggesting that either alpha2-containing or alpha3beta-GlyRs mediate the E2 effect observed in neurons. Furthermore, E2 inhibited the GlyR-mediated tonic current in pyramidal neurons of HIP CA1 region, where abundant GlyR alpha2 subunit is expressed. We suggest that the neuronal GlyR is a novel molecular target of E2 which directly inhibits the function of GlyRs in the HIP and SDH regions. This finding may shed new light on premenstrual dysphoric disorder and the gender differences in pain sensation at the CNS level.

Nudix hydrolases (Nd) is a widespread superfamily, which is found in all classes of organism, hydrolyse a wide range of organic pyrophosphates and has a 'housecleaning' function. The previous study showed that Trichinella spiralis Nd (TsNd) bound to intestinal epithelial cells (IECs), and the vaccination of mice with T7 phage-displayed TsNd polypeptides produced protective immunity. The aim of this study was to clone, express and identify the full-length TsNd and to investigate its immune protection against T. spiralis infection.

The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.

The ability to quantify cellular heterogeneity is a major advantage of single-cell technologies. However, statistical methods often treat cellular heterogeneity as a nuisance. We present a novel method to characterize differences in expression in the presence of distinct expression states within and among biological conditions. We demonstrate that this framework can detect differential expression patterns under a wide range of settings. Compared to existing approaches, this method has higher power to detect subtle differences in gene expression distributions that are more complex than a mean shift, and can characterize those differences. The freely available R package scDD implements the approach.

The orchestration of brain function requires complex gene regulatory networks that are modulated, in part, by microRNAs (miRNAs). These noncoding RNAs associate with argonaute (Ago) proteins in order to direct posttranscriptional gene suppression via base pairing with target transcripts. In order to better understand how miRNAs contribute to human-specialized brain processes and neurological phenotypes, identifying their targets is of paramount importance. Here, we address the latter by profiling Ago2:RNA interactions using HITS-CLIP to generate a transcriptome-wide map of miRNA binding sites in human brain. We uncovered ∼ 7,000 stringent Ago2 binding sites that are highly enriched for conserved sequences corresponding to abundant brain miRNAs. This interactome points to functional miRNA:target pairs across >3,000 genes and represents a valuable resource for accelerating our understanding of miRNA functions in brain. We demonstrate the utility of this map for exploring clinically relevant miRNA binding sites that may facilitate the translation of genetic studies of complex neuropsychiatric diseases into therapeutics.

Mesenchymal stem cells (MSCs) are multipotent progenitors, which give rise to several lineages, including bone, cartilage and fat. Epidermal growth factor (EGF) stimulates cell growth, proliferation and differentiation. EGF acts by binding with high affinity to epidermal growth factor receptor (EGFR) on the cell surface and stimulating the intrinsic protein tyrosine kinase activity of its receptor, which initiates a signal transduction cascade causing a variety of biochemical changes within the cell and regulating cell proliferation and differentiation. We have identified BMP9 as one of the most osteogenic BMPs in MSCs. In this study, we investigate if EGF signalling cross-talks with BMP9 and regulates BMP9-induced osteogenic differentiation. We find that EGF potentiates BMP9-induced early and late osteogenic markers of MSCs in vitro, which can be effectively blunted by EGFR inhibitors Gefitinib and Erlotinib or receptor tyrosine kinase inhibitors AG-1478 and AG-494 in a dose- and time-dependent manner. Furthermore, EGF significantly augments BMP9-induced bone formation in the cultured mouse foetal limb explants. In vivo stem cell implantation experiment reveals that exogenous expression of EGF in MSCs can effectively potentiate BMP9-induced ectopic bone formation, yielding larger and more mature bone masses. Interestingly, we find that, while EGF can induce BMP9 expression in MSCs, EGFR expression is directly up-regulated by BMP9 through Smad1/5/8 signalling pathway. Thus, the cross-talk between EGF and BMP9 signalling pathways in MSCs may underline their important roles in regulating osteogenic differentiation. Harnessing the synergy between BMP9 and EGF should be beneficial for enhancing osteogenesis in regenerative medicine.

Although it has been known for many years that T. spiralis muscle larvae (ML) can not invade intestinal epithelial cells unless they are exposed to the intestinal milieu and activated into intestinal infective larvae (IIL), which genes in IIL are involved in the process of invasion is still unknown. In this study, suppression subtractive hybridization (SSH) was performed to identify differentially expressed genes between IIL and ML. SSH library was constructed using cDNA generated from IIL as the 'tester'. About 110 positive clones were randomly selected from the library and sequenced, of which 33 T. spiralis genes were identified. Thirty encoded proteins were annotated according to Gene Ontology Annotation in terms of molecular function, biological process, and cellular localization. Out of 30 annotated proteins, 16 proteins (53.3%) had binding activity and 12 proteins (40.0%) had catalytic activity. The results of real-time PCR showed that the expression of nine genes (Ts7, Ndr family protein; Ts8, serine/threonine-protein kinase polo; Ts11, proteasome subunit beta type-7; Ts17, nudix hydrolase; Ts19, ovochymase-1; Ts22, fibronectin type III domain protein; Ts23, muscle cell intermediate filament protein OV71; Ts26, neutral and basic amino acid transport protein rBAT and Ts33, FACT complex subunit SPT16) from 33 T. spiralis genes in IIL were up-regulated compared with that of ML. The present study provide a group of the potential invasion-related candidate genes and will be helpful for further studies of mechanisms by which T. spiralis infective larvae recognize and invade the intestinal epithelial cells.

Combinatorial interplay among transcription factors (TFs) is an important mechanism by which transcriptional regulatory specificity is achieved. However, despite the increasing number of TFs for which either binding specificities or genome-wide occupancy data are known, knowledge about cooperativity between TFs remains limited. To address this, we developed a computational framework for predicting genome-wide co-binding between TFs (CCAT, Combinatorial Code Analysis Tool), and applied it to Drosophila melanogaster to uncover cooperativity among TFs during embryo development. Using publicly available TF binding specificity data and DNaseI chromatin accessibility data, we first predicted genome-wide binding sites for 324 TFs across five stages of D. melanogaster embryo development. We then applied CCAT in each of these developmental stages, and identified from 19 to 58 pairs of TFs in each stage whose predicted binding sites are significantly co-localized. We found that nearby binding sites for pairs of TFs predicted to cooperate were enriched in regions bound in relevant ChIP experiments, and were more evolutionarily conserved than other pairs. Further, we found that TFs tend to be co-localized with other TFs in a dynamic manner across developmental stages. All generated data as well as source code for our front-to-end pipeline are available at http://cat.princeton.edu.

Human pluripotent stem cells (hPSCs) have been differentiated to astroglia, but the utilization of hPSC-derived astroglia as cell therapy for neurological diseases has not been well studied. Astroglia are heterogeneous, and not all astroglia are equivalent in promoting neural repair. A prerequisite for cell therapy is to derive defined cell populations with superior therapeutic effects. Here we use an Olig2-GFP human embryonic stem cell (hESC) reporter to demonstrate that hESC-derived Olig2(+) progenitors generate a subtype of previously uncharacterized astroglia (Olig2PC-Astros). These Olig2PC-Astros differ substantially from astroglia differentiated from Olig2-negative hESC-derived neural progenitor cells (NPC-Astros), particularly in their neuroprotective properties. When grafted into brains subjected to global ischaemia, Olig2PC-Astros exhibit superior neuroprotective effects and improved behavioural outcome compared to NPC-Astros. Thus, this new paradigm of human astroglial differentiation is useful for studying the heterogeneity of human astroglia, and the unique Olig2PC-Astros may constitute a new cell therapy for treating cerebral ischaemia and other neurological diseases.

Oscillatory gene expression is fundamental to development, but technologies for monitoring expression oscillations are limited. We have developed a statistical approach called Oscope to identify and characterize the transcriptional dynamics of oscillating genes in single-cell RNA-seq data from an unsynchronized cell population. Applying Oscope to a number of data sets, we demonstrated its utility and also identified a potential artifact in the Fluidigm C1 platform.

Previous studies have shown that vitreous stimulates degradation of the tumor suppressor protein p53 and that knockdown of phosphatidylinositol 5-phosphate 4-kinases (PI5P4Kα and -β) abrogates proliferation of p53-deficient cells. The purpose of this study was to determine whether vitreous stimulated expression of PI5P4Kα and -β and whether suppression of PI5P4Kα and -β would inhibit vitreous-induced cellular responses and experimental proliferative vitreoretinopathy (PVR).

Astrocytes, once considered passive support cells, are increasingly appreciated as dynamic regulators of neuronal development and function, in part via secreted factors. The extent to which they similarly regulate oligodendrocytes or proliferation and differentiation of oligodendrocyte progenitor cells (OPCs) is less understood. Here, we generated astrocytes from human pluripotent stem cells (hiPSC-Astros) and demonstrated that immature astrocytes, as opposed to mature ones, promote oligodendrogenesis in vitro. In the PVL mouse model of neonatal hypoxic/ischemic encephalopathy, associated with cerebral palsy in humans, transplanted immature hiPSC-Astros promoted myelinogenesis and behavioral outcome. We further identified TIMP-1 as a selectively upregulated component secreted from immature hiPSC-Astros. Accordingly, in the rat PVL model, intranasal administration of conditioned medium from immature hiPSC-Astros promoted oligodendrocyte maturation in a TIMP-1-dependent manner. Our findings suggest stage-specific developmental interactions between astroglia and oligodendroglia and have important therapeutic implications for promoting myelinogenesis.

Rheumatoid arthritis (RA) affects millions world-wide. While anti-TNF treatment is widely used to reduce disease progression, treatment fails in ∼one-third of patients. No biomarker currently exists that identifies non-responders before treatment. A rigorous community-based assessment of the utility of SNP data for predicting anti-TNF treatment efficacy in RA patients was performed in the context of a DREAM Challenge (http://www.synapse.org/RA_Challenge). An open challenge framework enabled the comparative evaluation of predictions developed by 73 research groups using the most comprehensive available data and covering a wide range of state-of-the-art modelling methodologies. Despite a significant genetic heritability estimate of treatment non-response trait (h(2)=0.18, P value=0.02), no significant genetic contribution to prediction accuracy is observed. Results formally confirm the expectations of the rheumatology community that SNP information does not significantly improve predictive performance relative to standard clinical traits, thereby justifying a refocusing of future efforts on collection of other data.

To identify the biomarkers in early period of acute myocardial infarction (AMI) in rat serum and reveal the effective mechanism of a Traditional Chinese Medicine (TCM) named Shexiang Baoxin Pill (SBP).

In a broad range of organisms, Piwi-interacting RNAs (piRNAs) have emerged as core components of a surveillance system that protects the genome by silencing transposable and repetitive elements. A vast proportion of piRNAs is produced from discrete genomic loci, termed piRNA clusters, which are generally embedded in heterochromatic regions. The molecular mechanisms and the factors that govern their expression are largely unknown. Here, we show that Cutoff (Cuff), a Drosophila protein related to the yeast transcription termination factor Rai1, is essential for piRNA production in germline tissues. Cuff accumulates at centromeric/pericentromeric positions in germ-cell nuclei and strongly colocalizes with the major heterochromatic domains. Remarkably, we show that Cuff is enriched at the dual-strand piRNA cluster 1/42AB and is likely to be involved in regulation of transcript levels of similar loci dispersed in the genome. Consistent with this observation, Cuff physically interacts with the Heterochromatin Protein 1 (HP1) variant Rhino (Rhi). Our results unveil a link between Cuff activity, heterochromatin assembly and piRNA cluster expression, which is critical for stem-cell and germ-cell development in Drosophila.

Helicobacter pylori infection is associated with several extragastric conditions including dyslipidemia and metabolic syndrome. This study aimed to investigate additional metabolic parameters associated with H. pylori infection in a Chinese population.

Human pluripotent stem cells (hPSCs) are susceptible to numerical and structural chromosomal alterations during long-term culture. We show that mitotic errors occur frequently in hPSCs and that prometaphase arrest leads to very rapid apoptosis in undifferentiated but not in differentiated cells. hPSCs express high levels of proapoptotic protein NOXA in undifferentiated state. Knocking out NOXA by CRISPR or upregulation of the anti-apoptosis gene BCL-XL significantly reduced mitotic cell death, allowing the survival of aneuploid cells and the formation of teratomas significantly larger than their wild-type parental hPSCs. These results indicate that the normally low threshold of apoptosis in hPSCs can safeguard their genome integrity by clearing cells undergoing abnormal division. The amplification of BCL2L1 on chromosome 20q11.21, a frequent mutation in hPSCs, although not directly oncogenic, reduces the sensitivity of hPSCs to damage caused by erroneous mitosis and increases the risk of gaining aneuploidy.