Expression analysis of G Protein-Coupled Receptors in mouse macrophages

Pruning of the adipocyte peroxisome proliferator-activated receptor γ cistrome by hematopoietic master regulator PU.1

"Master" transcription factors are the gatekeepers of lineage identity. As such, they have been a major focus of efforts to manipulate cell fate for therapeutic purposes. The ETS transcription factor PU.1 has a potent ability to confer macrophage phenotypes on cells already committed to a different lineage, but how it overcomes the presence of other master regulators is not known. The nuclear receptor peroxisome proliferator-activated receptor γ (PPARγ) is the master regulator of the adipose lineage, and its genomic binding pattern in adipocytes is well characterized. Here we show that, when expressed at macrophage levels in mature adipocytes, PU.1 bound a large fraction of its macrophage sites, where it induced chromatin opening and the expression of macrophage target genes. Strikingly, PU.1 markedly reduced the genomic binding of PPARγ without changing its abundance. PU.1 expression repressed genes with nearby adipocyte-specific PPARγ binding sites, while a common macrophage-adipocyte gene expression program was retained. Together, these data reveal unexpected lability within the adipocyte PPARγ cistrome and show that, even in terminally differentiated cells, PU.1 can remodel the cistrome of another master regulator.

Several classical mouse inbred strains, including DBA/2, NOD/Lt, FVB/N, and SJL/J, carry a putative loss-of-function allele of Gpr84

G protein-coupled receptor 84 (GPR84) is a 7-transmembrane protein expressed on myeloid cells that can bind to medium-chain free fatty acids in vitro. Here, we report the discovery of a 2-bp frameshift deletion in the second exon of the Gpr84 gene in several classical mouse inbred strains. This deletion generates a premature stop codon predicted to result in a truncated protein lacking the transmembrane domains 4-7. We sequenced Gpr84 exon 2 from 58 strains representing different groups in the mouse family tree and found that 14 strains are homozygous for the deletion. Some of these strains are DBA/1J, DBA/2J, FVB/NJ, LG/J, MRL/MpJ, NOD/LtJ, and SJL/J. However, the deletion was not found in any of the wild-derived inbred strains analyzed. Haplotype analysis suggested that the deletion originates from a unique mutation event that occurred more than 100 years ago, preceding the development of the first inbred strain (DBA), from a Mus musculus domesticus source. As GPR84 ostensibly plays a role in the biology of myeloid cells, it could be relevant 1) to consider the existence of this Gpr84 nonsense mutation in several mouse strains when choosing a mouse model to study immune processes and 2) to consider reevaluating data obtained using such strains.

Bcl10 links saturated fat overnutrition with hepatocellular NF-kB activation and insulin resistance

Excess serum free fatty acids (FFAs) are fundamental to the pathogenesis of insulin resistance. With high-fat feeding, FFAs activate NF-kB in target tissues, initiating negative crosstalk with insulin signaling. However, the mechanisms underlying FFA-dependent NF-kB activation remain unclear. Here, we demonstrate that the saturated FA, palmitate, requires Bcl10 for NF-kB activation in hepatocytes. Uptake of palmitate, metabolism to diacylglycerol, and subsequent activation of protein kinase C (PKC) appear to mechanistically link palmitate with Bcl10, known as a central component of a signaling complex that, along with CARMA3 and MALT1, activates NF-kB downstream of selected cell surface receptors. Consequently, Bcl10-deficient mice are protected from hepatic NF-kB activation and insulin resistance following brief high-fat diet, suggesting that Bcl10 plays a major role in the metabolic consequences of acute overnutrition. Surprisingly, while CARMA3 also participates in the palmitate response, MALT1 is completely dispensable, thereby revealing an apparent nonclassical role for Bcl10 in NF-kB signaling.

Noncoding variation of the gene for ferritin light chain in hereditary and age-related cataract

PURPOSE

Cataract is a clinically and genetically heterogeneous disorder of the ocular lens and an important cause of visual impairment. The aim of this study was to map and identify the gene underlying autosomal dominant cataract segregating in a four-generation family, determine the lens expression profile of the identified gene, and test for its association with age-related cataract in a case-control cohort.

METHODS

Genomic DNA was prepared from blood leukocytes, and genotyping was performed by means of single-nucleotide polymorphism markers and microsatellite markers. Linkage analyses were performed using the GeneHunter and MLINK programs, and mutation detection was achieved by dideoxy cycle sequencing. Lens expression studies were performed using reverse-transcription polymerase chain reaction (RT-PCR) and in situ hybridization.

RESULTS

Genome-wide linkage analysis with single nucleotide polymorphism markers in the family identified a likely disease-haplotype interval on chromosome 19q (rs888861-[~17Mb]-rs8111640) that encompassed the microsatellite marker D19S879 (logarithm of the odds score [Z]=2.03, recombination distance [θ]=0). Mutation profiling of positional-candidate genes detected a heterozygous, noncoding G-to-T transversion (c.-168G>T) located in the iron response element (IRE) of the gene coding for ferritin light chain (FTL) that cosegregated with cataract in the family. Serum ferritin levels were found to be abnormally elevated (~fourfold), without evidence of iron overload, in an affected family member; this was consistent with a diagnosis of hereditary hyperferritinemia-cataract syndrome. No sequence variations located within the IRE were detected in a cohort of 197 cases with age-related cataract and 102 controls with clear lenses. Expression studies of human FTL, and its mouse counterpart FTL1, in the lens detected RT-PCR amplicons containing full-length protein-coding regions, and strong in situ localization of FTL1 transcripts to the lens equatorial epithelium and peripheral cortex.

CONCLUSIONS

The data are consistent with robust transcription of FTL in the lens, and suggest that whereas variations clustered in the IRE of the FTL gene are directly associated with hereditary hyperferritinemia-cataract syndrome, such IRE variations are unlikely to play a significant role in the genetic etiology of age-related cataract.

Polyunsaturated fatty acids and their metabolites in the pathobiology of schizophrenia

Schizophrenia can be considered as a low-grade systemic inflammatory disease with its origins in the perinatal period. It is likely that genetic, environmental, and nutritional factors interact to induce excess production of pro-inflammatory cytokines that, in turn, damage fetal neurons leading to the adult onset of schizophrenia. Polyunsaturated fatty acids (PUFAs) and their metabolites such as lipoxins, resolvins, protectins, maresins and nitrolipids not only have potent neuroprotective action but also are capable of inhibiting the production of pro-inflammatory cytokines. Decreased formation of PUFAs as a result of low activity of Δ(6) and Δ(5) desaturases can result in an increase in the production of pro-inflammatory cytokines due to the absence of negative control exerted by PUFAs and their anti-inflammatory metabolites that, in turn, may predispose to neuronal damage and development of schizophrenia in adult life. Furthermore, PUFAs are essential for brain growth and development. If this proposal is correct, this implies that perinatal and adult supplementation of PUFAs not only prevents but also helps in the treatment of schizophrenia. Furthermore, synthetic analogs of lipoxins, resolvins, and protectins may be of significant benefit in schizophrenia.

Development, composition, and structural arrangements of the ciliary zonule of the mouse

PURPOSE

Here, we examined the development, composition, and structural organization of the ciliary zonule of the mouse. Fibrillin 1, a large glycoprotein enriched in force-bearing tissues, is a prominent constituent of the mouse zonule. In humans, mutations in the gene for fibrillin 1 (FBN1) underlie Marfan syndrome (MS), a disorder characterized by lens dislocation and other ocular symptoms.

METHODS

Fibrillin expression was analyzed by in situ hybridization. The organization of the zonule was visualized using antibodies to Fbn1, Fbn2, and microfibril-associated glycoprotein-1 (Magp1) in conjunction with 5-ethynyl-2'-deoxyuridine (EdU), an S-phase marker.

RESULTS

Microfibrils, enriched in Fbn2 and Magp1, were prominent components of the temporary vascular tunic of the embryonic lens. Fbn2 expression by nonpigmented ciliary epithelial cells diminished postnatally and there was a concomitant increase in Fbn1 expression, especially in cells located in valleys between the ciliary folds. Zonular fibers projected from the posterior pars plicata to the lens in anterior, equatorial, and posterior groupings. The attachment point of the posterior zonular fibers consisted of a dense meshwork of radially oriented microfibrils that we termed the fibrillar girdle. The fibrillar girdle was located directly above the transition zone, a region of the lens epithelium in which cells commit to terminal differentiation.

CONCLUSIONS

The development and arrangement of the murine ciliary zonule are similar to those of humans, and consequently the mouse eye may be a useful model in which to study ocular complications of MS.

Diacylglycerol kinases: regulated controllers of T cell activation, function, and development

Diacylglycerol kinases (DGKs) are a diverse family of enzymes that catalyze the conversion of diacylglycerol (DAG), a crucial second messenger of receptor-mediated signaling, to phosphatidic acid (PA). Both DAG and PA are bioactive molecules that regulate a wide set of intracellular signaling proteins involved in innate and adaptive immunity. Clear evidence points to a critical role for DGKs in modulating T cell activation, function, and development. More recently, studies have elucidated factors that control DGK function, suggesting an added complexity to how DGKs act during signaling. This review summarizes the available knowledge of the function and regulation of DGK isoforms in signal transduction with a particular focus on T lymphocytes.

Arachidonic acid and lipoxin A4 as possible endogenous anti-diabetic molecules

In both type 1 and type 2 diabetes mellitus, increased production of pro-inflammatory cytokines and reactive oxygen species (ROS) occurs that induce apoptosis of β cells and cause peripheral insulin resistance respectively though the degree of their increased production is higher in type 1 and less in type 2 diabetes mellitus. Despite this, the exact mechanism(s) that lead to increased production of pro-inflammatory cytokines: interleukin-6 (IL-6) and tumor necrosis factor-α (TNF-α) and ROS is not known. Studies showed that plasma concentrations of arachidonic acid (AA) and lipoxin A4 (LXA4) are low in alloxan-induced type 1 diabetes mellitus in experimental animals and patients with type 2 diabetes mellitus. Prior administration of AA, eicosapentaenoic and docosahexaenoic acids (EPA and DHA, respectively) and transgenic animals that produce increased amounts of EPA and DHA acids were protected from chemical-induced diabetes mellitus that was associated with enhanced formation of LXA4 and resolvins, while protectin D1 ameliorated peripheral insulin resistance. AA, LXA4, resolvins and protectins inhibit IL-6 and TNF-α production and suppress ROS generation. Thus, AA and lipoxins, resolvins and protectins may function as endogenous anti-diabetic molecules implying that their administration could be useful in the prevention and management of both types of diabetes mellitus.

Medium-chain fatty acid-sensing receptor, GPR84, is a proinflammatory receptor

G protein-coupled receptor 84 (GPR84) is a putative receptor for medium-chain fatty acids (MCFAs), whose pathophysiological roles have not yet been clarified. Here, we show that GPR84 was activated by MCFAs with the hydroxyl group at the 2- or 3-position more effectively than nonhydroxylated MCFAs. We also identified a surrogate agonist, 6-n-octylaminouracil (6-OAU), for GPR84. These potential ligands and the surrogate agonist, 6-OAU, stimulated [(35)S]GTP binding and accumulated phosphoinositides in a GPR84-dependent manner. The surrogate agonist, 6-OAU, internalized GPR84-EGFP from the cell surface. Both the potential ligands and 6-OAU elicited chemotaxis of human polymorphonuclear leukocytes (PMNs) and macrophages and amplified LPS-stimulated production of the proinflammatory cytokine IL-8 from PMNs and TNFα from macrophages. Furthermore, the intravenous injection of 6-OAU raised the blood CXCL1 level in rats, and the inoculation of 6-OAU into the rat air pouch accumulated PMNs and macrophages in the site. Our results indicate a proinflammatory role of GPR84, suggesting that the receptor may be a novel target to treat chronic low grade inflammation associated-disease.

MICU2, a paralog of MICU1, resides within the mitochondrial uniporter complex to regulate calcium handling

Mitochondrial calcium uptake is present in nearly all vertebrate tissues and is believed to be critical in shaping calcium signaling, regulating ATP synthesis and controlling cell death. Calcium uptake occurs through a channel called the uniporter that resides in the inner mitochondrial membrane. Recently, we used comparative genomics to identify MICU1 and MCU as the key regulatory and putative pore-forming subunits of this channel, respectively. Using bioinformatics, we now report that the human genome encodes two additional paralogs of MICU1, which we call MICU2 and MICU3, each of which likely arose by gene duplication and exhibits distinct patterns of organ expression. We demonstrate that MICU1 and MICU2 are expressed in HeLa and HEK293T cells, and provide multiple lines of biochemical evidence that MCU, MICU1 and MICU2 reside within a complex and cross-stabilize each other's protein expression in a cell-type dependent manner. Using in vivo RNAi technology to silence MICU1, MICU2 or both proteins in mouse liver, we observe an additive impairment in calcium handling without adversely impacting mitochondrial respiration or membrane potential. The results identify MICU2 as a new component of the uniporter complex that may contribute to the tissue-specific regulation of this channel.

Dynamic changes in lipid droplet-associated proteins in the "browning" of white adipose tissues

The morphological and functional differences between lipid droplets (LDs) in brown (BAT) and white (WAT) adipose tissues will largely be determined by their associated proteins. Analysing mRNA expression in mice fat depots we have found that most LD protein genes are expressed at higher levels in BAT, with the greatest differences observed for Cidea and Plin5. Prolonged cold exposure, which induces the appearance of brown-like adipocytes in mice WAT depots, was accompanied with the potentiation of the lipolytic machinery, with changes in ATGL, CGI-58 and G0S2 gene expression. However the major change detected in WAT was the enhancement of Cidea mRNA. Together with the increase in Cidec, it indicates that LD enlargement through LD-LD transference of fat is an important process during WAT browning. To study the dynamics of this phenotypic change, we have applied 4D confocal microscopy in differentiated 3T3-L1 cells under sustained β-adrenergic stimulation. Under these conditions the cells experienced a LD remodelling cycle, with progressive reduction on the LD size by lipolysis, followed by the formation of new LDs, which were subjected to an enlargement process, likely to be CIDE-triggered, until the cell returned to the basal state. This transformation would be triggered by the activation of a thermogenic futile cycle of lipolysis/lipogenesis and could facilitate the molecular mechanism for the unilocular to multilocular transformation during WAT browning. This article is part of a Special Issue entitled Brown and White Fat: From Signaling to Disease.

Spatial coupling of JNK activation to the B cell antigen receptor by tyrosine-phosphorylated ezrin

The ezrin-radixin-moesin proteins regulate B lymphocyte activation via their effect on BCR diffusion and microclustering. This relies on their ability to dynamically tether the plasma membrane with actin filaments that is in turn facilitated by phosphorylation of the conserved threonine residue in the actin-binding domain. In this study, we describe a novel function of ezrin in regulating JNK activation that is mediated by phosphorylation of a tyrosine (Y353) residue that is unconserved with moesin and radixin. BCR, but not CD40, TLR4, or CXCR5 stimulation, induced phosphorylation of ezrin at Y353 in mouse splenic B cells. Ezrin existed in a preformed complex with Syk in unstimulated B cells and underwent Syk-dependent phosphorylation upon anti-IgM stimulation. Y353-phosphorylated ezrin colocalized with the BCR within minutes of stimulation and cotrafficked with the endocytosed BCRs through the early and late endosomes. The T567 residue of ezrin was rephosphorylated in late endosomes and at the plasma membrane at later times of BCR stimulation. Expression of a nonphosphorylatable Y353F mutant of ezrin specifically impaired JNK activation. BCR crosslinking induced the association of Y353-phosphorylated ezrin with JNK and its kinase MAPKK7, as well as spatial colocalization with phosphorylated JNK in the endosomes. The yellow fluorescent protein-tagged Y353F mutant displayed reduced colocalization with the endocytosed BCR as compared with wild-type ezrin-yellow fluorescent protein. Taken together, our data identify a novel role for ezrin as a spatial adaptor that couples JNK signaling components to the BCR signalosome, thus facilitating JNK activation.

Systems-level analysis of genome-wide association data

Genome-wide association studies (GWAS) have emerged as the method of choice for identifying common variants affecting complex disease. In a GWAS, particular attention is placed, for obvious reasons, on single-nucleotide polymorphisms (SNPs) that exceed stringent genome-wide significance thresholds. However, it is expected that many SNPs with only nominal evidence of association (e.g., P < 0.05) truly influence disease. Efforts to extract additional biological information from entire GWAS datasets have primarily focused on pathway-enrichment analyses. However, these methods suffer from a number of limitations and typically fail to lead to testable hypotheses. To evaluate alternative approaches, we performed a systems-level analysis of GWAS data using weighted gene coexpression network analysis. A weighted gene coexpression network was generated for 1918 genes harboring SNPs that displayed nominal evidence of association (P ≤ 0.05) from a GWAS of bone mineral density (BMD) using microarray data on circulating monocytes isolated from individuals with extremely low or high BMD. Thirteen distinct gene modules were identified, each comprising coexpressed and highly interconnected GWAS genes. Through the characterization of module content and topology, we illustrate how network analysis can be used to discover disease-associated subnetworks and characterize novel interactions for genes with a known role in the regulation of BMD. In addition, we provide evidence that network metrics can be used as a prioritizing tool when selecting genes and SNPs for replication studies. Our results highlight the advantages of using systems-level strategies to add value to and inform GWAS.

Development of Physiologically Based Pharmacokinetic Model (PBPK) of BMP2 in Mice

Bone Morphogenetic protein 2 holds great promise for potential applications in the clinic. It is a potent growth factor for the use in the cervical spine surgery (FDA approved 2002) and has been marketed as "Infuse" for treating open tibial shaft fractures (FDA approved 2004). However, its use is limited by several significant side effects that maybe due to its potency and effect on different stem cell populations in the spine. BMP2 is expressed throughout the human body in several tissues and at a very high concentration in the blood. BMP receptors, especially BMP receptor type Ia, is ubiquitously expressed in most tissues. Currently, it is difficult to determine how BMP2 is physiologically distributed in mice or humans and no quantitative models are available. A Physiologically-Based Pharmaco-Kinetic (PBPK) model has been developed to determine steady-state distribution of BMP2 in mice. The multi-compartmental PBPK model represents relevant organ/tissues with physiological accuracy. The organs/tissue compartments chosen were brain, lung, heart, liver, pancreas, kidney, uterus, bone and fat. A blood compartment maintained connectivity among the various organs. Four processes characterized the change in the concentration of the protein in every compartment: blood flow in, blood flow out, protein turnover and receptor binding in the organ. The unique aspects of the model are the determination of elimination using receptor kinetics and generation using protein turnover. The model also predicts steady state concentrations of BMP2 in tissues in mice and may be used for possible scale-up of dosage regimens in humans.

Systems genetic analysis of osteoblast-lineage cells

The osteoblast-lineage consists of cells at various stages of maturation that are essential for skeletal development, growth, and maintenance. Over the past decade, many of the signaling cascades that regulate this lineage have been elucidated; however, little is known of the networks that coordinate, modulate, and transmit these signals. Here, we identify a gene network specific to the osteoblast-lineage through the reconstruction of a bone co-expression network using microarray profiles collected on 96 Hybrid Mouse Diversity Panel (HMDP) inbred strains. Of the 21 modules that comprised the bone network, module 9 (M9) contained genes that were highly correlated with prototypical osteoblast maker genes and were more highly expressed in osteoblasts relative to other bone cells. In addition, the M9 contained many of the key genes that define the osteoblast-lineage, which together suggested that it was specific to this lineage. To use the M9 to identify novel osteoblast genes and highlight its biological relevance, we knocked-down the expression of its two most connected “hub” genes, Maged1 and Pard6g. Their perturbation altered both osteoblast proliferation and differentiation. Furthermore, we demonstrated the mice deficient in Maged1 had decreased bone mineral density (BMD). It was also discovered that a local expression quantitative trait locus (eQTL) regulating the Wnt signaling antagonist Sfrp1 was a key driver of the M9. We also show that the M9 is associated with BMD in the HMDP and is enriched for genes implicated in the regulation of human BMD through genome-wide association studies. In conclusion, we have identified a physiologically relevant gene network and used it to discover novel genes and regulatory mechanisms involved in the function of osteoblast-lineage cells. Our results highlight the power of harnessing natural genetic variation to generate co-expression networks that can be used to gain insight into the function of specific cell-types.

The osteoblast-lineage consists of a range of cells from osteogenic precursors that mature into bone-forming osteoblasts to osteocytes that are entombed in bone. Each cell in the lineage serves a number of distinct and critical roles in the growth and maintenance of the skeleton, as well as many extra-skeletal functions. Over the last decade, many of the major regulatory pathways governing the differentiation and activity of these cells have been discovered. In contrast, little is known regarding the composition or function of gene networks within the lineage. The goal of this study was to increase our understanding of how genes are organized into networks in osteoblasts. Towards this goal, we used microarray gene expression profiles from bone to identify a group of genes that formed a network specific to the osteoblast-lineage. We used the knowledge of this network to identify novel genes that are important for regulating various aspects of osteoblast function. These data improve our understanding of the gene networks operative in cells of the osteoblast-lineage.

The intramembrane protease Sppl2a is required for B cell and DC development and survival via cleavage of the invariant chain

B cell development requires tight regulation to allow for the generation of a diverse repertoire while preventing the development of autoreactive cells. We report, using N-ethyl-N-nitrosourea (ENU)-induced mutagenesis, the identification of a mutant mouse (chompB) with a block in early B cell development. The blockade occurs after the transitional 1 (T1) stage and leads to a decrease in mature B cell subsets and deficits in T cell-dependent antibody responses. Additionally, chompB mice have decreases in myeloid dendritic cells (DCs). The mutation was mapped to the intramembrane protease signal peptide peptidase-like 2a (Sppl2a), a gene not previously implicated in immune cell development. Proteomic analysis identified the invariant chain (CD74) as a key substrate of Sppl2a and suggests that regulated intramembrane proteolysis of CD74 by Sppl2a contributes to B cell and DC survival. Moreover, these data suggest that modulation of Sppl2a may be a useful therapeutic strategy for treatment of B cell dependent autoimmune disorders.

High-throughput identification and characterization of novel, species-selective GPR35 agonists

Drugs targeting the orphan receptor GPR35 have potential therapeutic application in a number of disease areas, including inflammation, metabolic disorders, nociception, and cardiovascular disease. Currently available surrogate GPR35 agonists identified from pharmacologically relevant compound libraries have limited utility due to the likelihood of off-target effects in vitro and in vivo and the variable potency that such ligands exhibit across species. We sought to identify and characterize novel GPR35 agonists to facilitate studies aimed at defining the physiologic role of GPR35. PathHunter β-arrestin recruitment technology was validated as a human GPR35 screening assay, and a high-throughput screen of 100,000 diverse low molecular weight compounds was conducted. Confirmed GPR35 agonists from five distinct chemotypes were selected for detailed characterization using both β-arrestin recruitment and G protein-dependent assays and each of the human, mouse, and rat GPR35 orthologs. These studies identified 4-{(Z)-[(2Z)-2-(2-fluorobenzylidene)-4-oxo-1,3-thiazolidin-5-ylidene]methyl}benzoic acid (compound 1) as the highest potency full agonist of human GPR35 yet described. As with certain other GPR35 agonists, compound 1 was markedly selective for human GPR35, but displayed elements of signal bias between β-arrestin-2 and G protein-dependent assays. Compound 1 also displayed competitive behavior when assessed against the human GPR35 antagonist, ML-145 (2-hydroxy-4-[4-(5Z)-5-[(E)-2-methyl-3-phenylprop-2-enylidene]-4-oxo-2-sulfanylidene-1,3-thiazolidin-3-yl]butanoylamino]benzoic acid). Of the other chemotypes studied, compounds 2 and 3 were selective for the human receptor, but compounds 4 and 5 demonstrated similar activity at human, rat, and mouse GPR35 orthologs. Further characterization of these compounds and related analogs is likely to facilitate a better understanding of GPR35 in health and disease.

The arthritis severity locus Cia5a regulates the expression of inflammatory mediators including Syk pathway genes and proteases in pristane-induced arthritis

Background

Cia5a is a locus on rat chromosome 10 that regulates disease severity and joint damage in two models of rheumatoid arthritis, collagen- and pristane-induced arthritis (PIA). In this study, we aimed to identify cellular and molecular processes regulated by Cia5a using microarray-based gene expression analysis of synovial tissues from MHC identical DA (severe erosive disease) and DA.F344(Cia5a) congenics (mild non-erosive disease) rats.

Results

Synovial tissues from six DA and eight DA.F344(Cia5a) rats were analyzed 21 days after the induction of PIA using the Illumina RatRef-12 BeadChip (21,922 genes) and selected data confirmed with qPCR. There was a significantly increased expression of pro-inflammatory mediators such as Il1b (5-fold), Il18 (3.9-fold), Cxcl1 (10-fold), Cxcl13 (7.5-fold) and Ccl7 (7.9-fold), and proteases like Mmp3 (23-fold), Mmp9 (32-fold), Mmp14 (4.4-fold) and cathepsins in synovial tissues from DA, with reciprocally reduced levels in congenics. mRNA levels of 47 members of the Spleen Tyrosine Kinase (Syk) pathway were significantly increased in DA synovial tissues compared with DA.F344(Cia5a), and included Syk (5.4-fold), Syk-activating receptors and interacting proteins, and genes regulated by Syk such as NFkB, and NAPDH oxidase complex genes. Nuclear receptors (NR) such as Rxrg, Pparg and Rev-erba were increased in the protected congenics, and so was the anti-inflammatory NR-target gene Scd1 (54-fold increase). Tnn (72-fold decrease) was the gene most significantly increased in DA.

Conclusions

Analyses of gene expression in synovial tissues revealed that the arthritis severity locus Cia5a regulates the expression of key mediators of inflammation and joint damage, as well as the expression of members of the Syk pathway. This expression pattern correlates with disease severity and joint damage and along with the gene accounting for Cia5a could become a useful biomarker to identify patients at increased risk for severe and erosive disease. The identification of the gene accounting for Cia5a has the potential to generate a new and important target for therapy and prognosis.

The transcription factor paired-related homeobox 1 (Prrx1) inhibits adipogenesis by activating transforming growth factor-β (TGFβ) signaling

Differentiation of adipocytes from preadipocytes contributes to adipose tissue expansion in obesity. Impaired adipogenesis may underlie the development of metabolic diseases such as insulin resistance and type 2 diabetes. Mechanistically, a well defined transcriptional network coordinates adipocyte differentiation. The family of paired-related homeobox transcription factors, which includes Prrx1a, Prrx1b, and Prrx2, is implicated with regulation of mesenchymal cell fate, including myogenesis and skeletogenesis; however, whether these proteins impact adipogenesis remains to be addressed. In this study, we identify Prrx1a and Prrx1b as negative regulators of adipogenesis. We show that Prrx1a and Prrx1b are down-regulated during adipogenesis in vitro and in vivo. Stable knockdown of Prrx1a/b enhances adipogenesis, with increased expression of peroxisome proliferator-activated receptor-γ, CCAAT/enhancer-binding protein-α and FABP4 and increased secretion of the adipokines adiponectin and chemerin. Although stable low-level expression of Prrx1a, Prrx1b, or Prrx2 does not affect 3T3-L1 adipogenesis, transient overexpression of Prrx1a or Prrx1b inhibits peroxisome proliferator-activated receptor-γ activity. Prrx1 knockdown decreases expression of Tgfb2 and Tgfb3, and inhibition of TGFβ signaling during adipogenesis mimics the effects of Prrx1 knockdown. These data support the hypothesis that endogenous Prrx1 restrains adipogenesis by regulating expression of TGFβ ligands and thereby activating TGFβ signaling. Finally, we find that expression of Prrx1a or Prrx1b in adipose tissue increases during obesity and strongly correlates with Tgfb3 expression in BL6 mice. These observations suggest that increased Prrx1 expression may promote TGFβ activity in adipose tissue and thereby contribute to aberrant adipocyte function during obesity.

Mouse Nkrp1-Clr gene cluster sequence and expression analyses reveal conservation of tissue-specific MHC-independent immunosurveillance

The Nkrp1 (Klrb1)-Clr (Clec2) genes encode a receptor-ligand system utilized by NK cells as an MHC-independent immunosurveillance strategy for innate immune responses. The related Ly49 family of MHC-I receptors displays extreme allelic polymorphism and haplotype plasticity. In contrast, previous BAC-mapping and aCGH studies in the mouse suggest the neighboring and related Nkrp1-Clr cluster is evolutionarily stable. To definitively compare the relative evolutionary rate of Nkrp1-Clr vs. Ly49 gene clusters, the Nkrp1-Clr gene clusters from two Ly49 haplotype-disparate inbred mouse strains, BALB/c and 129S6, were sequenced. Both Nkrp1-Clr gene cluster sequences are highly similar to the C57BL/6 reference sequence, displaying the same gene numbers and order, complete pseudogenes, and gene fragments. The Nkrp1-Clr clusters contain a strikingly dissimilar proportion of repetitive elements compared to the Ly49 clusters, suggesting that certain elements may be partly responsible for the highly disparate Ly49 vs. Nkrp1 evolutionary rate. Focused allelic polymorphisms were found within the Nkrp1b/d (Klrb1b), Nkrp1c (Klrb1c), and Clr-c (Clec2f) genes, suggestive of possible immune selection. Cell-type specific transcription of Nkrp1-Clr genes in a large panel of tissues/organs was determined. Clr-b (Clec2d) and Clr-g (Clec2i) showed wide expression, while other Clr genes showed more tissue-specific expression patterns. In situ hybridization revealed specific expression of various members of the Clr family in leukocytes/hematopoietic cells of immune organs, various tissue-restricted epithelial cells (including intestinal, kidney tubular, lung, and corneal progenitor epithelial cells), as well as myocytes. In summary, the Nkrp1-Clr gene cluster appears to evolve more slowly relative to the related Ly49 cluster, and likely regulates innate immunosurveillance in a tissue-specific manner.

A gene expression atlas of the domestic pig

BACKGROUND

This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering.

RESULTS

The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development.

CONCLUSIONS

As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites http://biogps.org and http://www.macrophages.com/pig-atlas.

A mutation in the FAM36A gene, the human ortholog of COX20, impairs cytochrome c oxidase assembly and is associated with ataxia and muscle hypotonia

The mitochondrial respiratory chain complex IV (cytochrome c oxidase) is a multi-subunit enzyme that transfers electrons from cytochrome c to molecular oxygen, yielding water. Its biogenesis requires concerted expression of mitochondria- and nuclear-encoded subunits and assembly factors. In this report, we describe a homozygous missense mutation in FAM36A from a patient who displays ataxia and muscle hypotonia. The FAM36A gene is a remote, putative ortholog of the fungal complex IV assembly factor COX20. Messenger RNA (mRNA) and protein co-expression analyses support the involvement of FAM36A in complex IV function in mammals. The c.154A>C mutation in the FAM36A gene, a mutation that is absent in sequenced exomes, leads to a reduced activity and lower levels of complex IV and its protein subunits. The FAM36A protein is nearly absent in patient's fibroblasts. Cells affected by the mutation accumulate subassemblies of complex IV that contain COX1 but are almost devoid of COX2 protein. We observe co-purification of FAM36A and COX2 proteins, supporting that the FAM36A defect hampers the early step of complex IV assembly at the incorporation of the COX2 subunit. Lentiviral complementation of patient's fibroblasts with wild-type FAM36A increases the complex IV activity as well as the amount of holocomplex IV and of individual subunits. These results establish the function of the human gene FAM36A/COX20 in complex IV assembly and support a causal role of the gene in complex IV deficiency.

Gene and genome parameters of mammalian liver circadian genes (LCGs)

The mammalian circadian system controls various physiology processes and behavior responses by regulating thousands of circadian genes with rhythmic expressions. In this study, we redefined circadian-regulated genes based on published results in the mouse liver and compared them with other gene groups defined relative to circadian regulations, especially the non-circadian-regulated genes expressed in liver at multiple molecular levels from gene position to protein expression based on integrative analyses of different datasets from the literature. Based on the intra-tissue analysis, the liver circadian genes or LCGs show unique features when compared to other gene groups. First, LCGs in general have less neighboring genes and larger in both genomic and 3'-UTR lengths but shorter in CDS (coding sequence) lengths. Second, LCGs have higher mRNA and protein abundance, higher temporal expression variations, and shorter mRNA half-life. Third, more than 60% of LCGs form major co-expression clusters centered in four temporal windows: dawn, day, dusk, and night. In addition, larger and smaller LCGs are found mainly expressed in the day and night temporal windows, respectively, and we believe that LCGs are well-partitioned into the gene expression regulatory network that takes advantage of gene size, expression constraint, and chromosomal architecture. Based on inter-tissue analysis, more than half of LCGs are ubiquitously expressed in multiple tissues but only show rhythmical expression in one or limited number of tissues. LCGs show at least three-fold lower expression variations across the temporal windows than those among different tissues, and this observation suggests that temporal expression variations regulated by the circadian system is relatively subtle as compared with the tissue expression variations formed during development. Taken together, we suggest that the circadian system selects gene parameters in a cost effective way to improve tissue-specific functions by adapting temporal variations from the environment over evolutionary time scales.

5-hmC in the brain is abundant in synaptic genes and shows differences at the exon-intron boundary

5-hydroxymethylcytosine (5-hmC), a derivative of 5-methylcytosine (5-mC), is abundant in the brain for unknown reasons. Our goal was to characterize the genomic distribution of 5-hmC and 5-mC in human and mouse tissues. We assayed 5-hmC using glucosylation coupled with restriction enzyme digestion, and interrogation on microarrays. We detected 5-hmC enrichment in genes with synapse-related functions in both human and mouse brain. We also identified substantial tissue-specific differential distributions of these DNA modifications at the exon-intron boundary, in both human and mouse. This boundary change was mainly due to 5-hmC in the brain, but due to 5-mC in non-neural contexts. This pattern was replicated in multiple independent datasets and with single molecule sequencing. Moreover, in human frontal cortex, constitutive exons contained higher levels of 5-hmC, relative to alternatively-spliced exons. Our study suggests a novel role for 5-hmC in RNA splicing and synaptic function in the brain.

Comprehensive analysis of genes expressed by rare microchimeric fetal cells in the maternal mouse lung

During pregnancy, cells from each fetus travel into the maternal circulation and organs, resulting in the development of microchimerism. Identification of the cell types in this microchimeric population would permit better understanding of possible mechanisms by which they affect maternal health. However, comprehensive analysis of fetal cells has been hampered by their rarity. In this study, we sought to overcome this obstacle by combining flow cytometry with multidimensional gene expression microarray analysis of fetal cells isolated from the murine maternal lung during late pregnancy. Fetal cells were collected from the lungs of pregnant female mice. cDNA was amplified and hybridized to gene expression microarrays. The resulting fetal cell core transcriptome was interrogated using multiple methods including Ingenuity Pathway Analysis, the BioGPS gene expression database, principal component analysis, the Eurexpress gene expression atlas, and primary literature. Here we report that small numbers of fetal cells can be flow sorted from the maternal lung, facilitating discovery-driven gene expression analysis. We additionally show that gene expression data can provide functional information about fetal cells. Our results suggest that fetal cells in the murine maternal lung are a mixed population, consisting of trophoblasts, mesenchymal stem cells, and cells of the immune system. Detection of trophoblasts and immune cells in the maternal lung may facilitate future mechanistic studies related to the development of immune tolerance and pregnancy-related complications, such as pre-eclampsia. Furthermore, the presence and persistence of mesenchymal stem cells in maternal organs may have implications for long-term postpartum maternal health.

Role of M3 mAChR in in vivo and in vitro models of LPS-induced inflammatory response

OBJECTIVE

We tested the potential role of the mAChR in lipopolysaccharide (LPS)-induced inflammatory response in in vivo and in vitro models and a possible signaling pathway involved in the inflammatory process.

METHODS

Anesthetized mice were challenged with intratracheal LPS to induce acute lung injury. The cytology and histopathology changes, expression of cytokines and pulmonary vascular permeability were used to evaluate the effects of the cholinergic agent. Alveolar macrophage cell line NR8383 was also used to confirm the role of mAChRs and the molecular mechanisms underlying the LPS-induced events.

RESULTS

LPS-induced acute lung injury (ALI) was significantly improved by atropine (a non-selective mAChR antagonist) and 4-DAMP (a M3 mAChR antagonist), as indicated by the diminution of neutrophil infiltration, pulmonary vascular permeability and IL-6 and TNF-α production. LPS-induced TNF-α production from the alveolar macrophage was significantly inhibited by atropine and 4-DAMP, but not pirenzepine (a M1 mAChR antagonist) and methoctramine (a M2 mAChR antagonist). Interestingly, LPS-induced TNF-α production was enhanced by the muscarinic receptor agonist pilocarpine, and treatment with pilocarpine alone was able to trigger TNF-α production from the alveolar macrophage, which was effectively attenuated by 4-DAMP. Western blot analysis showed that LPS-induced degradation of IκBα was strongly blocked by atropine/4-DAMP both in vivo and in vitro, indicating that M3 mAChR was involved in LPS-induced lung inflammation by mediating the NF-κB signaling pathway.

CONCLUSION

Our findings bring the evidence that the blockage of mAChR exerts anti-inflammatory properties, in which the M3 mAChR plays an important role in the LPS-induced lung inflammation.

The BTB-ZF transcription factors

The BTB-ZF (broad-complex, tramtrack and bric-à-brac--zinc finger) proteins are encoded by at least 49 genes in mouse and man and commonly serve as sequence-specific silencers of gene expression. This review will focus on the known physiological functions of mammalian BTB-ZF proteins, which include essential roles in the development of the immune system. We discuss their function in terminally differentiated lymphocytes and the progenitors that give rise to them, their action in hematopoietic malignancy and roles beyond the immune system.

The Raine syndrome protein FAM20C is a Golgi kinase that phosphorylates bio-mineralization proteins

Raine syndrome is caused by mutations in FAM20C, which had been reported to encode a secreted component of bone and teeth. We found that FAM20C encodes a Golgi-localized protein kinase, distantly related to the Golgi-localized kinase Four-jointed. Drosophila also encode a Golgi-localized protein kinase closely related to FAM20C. We show that FAM20C can phosphorylate secreted phosphoproteins, including both Casein and members of the SIBLING protein family, which modulate biomineralization, and we find that FAM20C phosphorylates a biologically active peptide at amino acids essential for inhibition of biomineralization. We also identify autophosphorylation of FAM20C, and characterize parameters of FAM20C’s kinase activity, including its Km, pH and cation dependence, and substrate specificity. The biochemical properties of FAM20C match those of an enzymatic activity known as Golgi casein kinase. Introduction of point mutations identified in Raine syndrome patients into recombinant FAM20C impairs its normal localization and kinase activity. Our results identify FAM20C as a kinase for secreted phosphoproteins and establish a biochemical basis for Raine syndrome.

Identification of RIFL, a novel adipocyte-enriched insulin target gene with a role in lipid metabolism

To identify new genes that are important in fat metabolism, we utilized the Lexicon-Genentech knockout database of genes encoding transmembrane and secreted factors and whole murine genome transcriptional profiling data that we generated for 3T3-L1 in vitro adipogenesis. Cross-referencing null models evidencing metabolic phenotypes with genes induced in adipogenesis led to identification of a new gene, which we named RIFL (refeeding induced fat and liver). RIFL-null mice have serum triglyceride levels approximately one-third of wild type. RIFL transcript is induced >100-fold during 3T3-L1 adipogenesis and is also increased markedly during adipogenesis of murine and human primary preadipocytes. siRNA-mediated knockdown of RIFL during 3T3-L1 adipogenesis results in an ~35% decrease in adipocyte triglyceride content. Murine RIFL transcript is highly enriched in white and brown adipose tissue and liver. Fractionation of WAT reveals that RIFL transcript is exclusive to adipocytes with a lack of expression in stromal-vascular cells. Nutritional and hormonal studies are consistent with a prolipogenic function for RIFL. There is evidence of an approximately eightfold increase in RIFL transcript level in WAT in ob/ob mice compared with wild-type mice. RIFL transcript level in WAT and liver is increased ~80- and 12-fold, respectively, following refeeding of fasted mice. Treatment of 3T3-L1 adipocytes with insulin increases RIFL transcript ≤35-fold, whereas agents that stimulate lipolysis downregulate RIFL. Interestingly, the 198-amino acid RIFL protein is predicted to be secreted and shows ~30% overall conservation with the NH(2)-terminal half of angiopoietin-like 3, a liver-secreted protein that impacts lipid metabolism. In summary, our data suggest that RIFL is an important new regulator of lipid metabolism.

BioGPS and GXD: mouse gene expression data-the benefits and challenges of data integration

Mouse gene expression data are complex and voluminous. To maximize the utility of these data, they must be made readily accessible through databases, and those resources need to place the expression data in the larger biological context. Here we describe two community resources that approach these problems in different but complementary ways: BioGPS and the Mouse Gene Expression Database (GXD). BioGPS connects its large and homogeneous microarray gene expression reference data sets via plugins with a heterogeneous collection of external gene centric resources, thus casting a wide but loose net. GXD acquires different types of expression data from many sources and integrates these data tightly with other types of data in the Mouse Genome Informatics (MGI) resource, with a strong emphasis on consistency checks and manual curation. We describe and contrast the "loose" and "tight" data integration strategies employed by BioGPS and GXD, respectively, and discuss the challenges and benefits of data integration. BioGPS is freely available at http://biogps.org . GXD is freely available through the MGI web site ( www.informatics.jax.org ) or directly at www.informatics.jax.org/expression.shtml .

REST regulates the pool size of the different neural lineages by restricting the generation of neurons and oligodendrocytes from neural stem/progenitor cells

REST is a master repressor of neuronal genes; however, whether it has any role during nervous system development remains largely unknown. Here, we analyzed systematically the role of REST in embryonic stem cells and multipotent neural stem/progenitor (NS/P) cells, including neurogenic and gliogenic NS/P cells derived from embryonic stem (ES) cells or developing mouse embryos. We showed that REST-null ES cells remained pluripotent and generated teratomas consisting of the three germ layers. By contrast, multipotent NS/P cells lacking REST displayed significantly reduced self-renewal capacity owing to reduced cell cycle kinetics and precocious neuronal differentiation. Importantly, although early-born neurogenic NS/P cells that lack REST were capable of differentiating to neurons and glia, the neuronal and oligodendrocytic pools were significantly enlarged and the astrocytic pool was shrunken. However, gliogenic NS/P cells lacking REST were able to generate a normal astrocytic pool size, suggesting that the shrinkage of the astrocytic pool generated from neurogenic NS/P cells lacking REST probably occurs by default. Microarray profiling of early-born NS/P cells lacking REST showed upregulation of neuronal as well as oligodendrocytic genes, specifically those involved in myelination. Furthermore, chromatin immunoprecipitation analyses showed that some of the upregulated oligodendrocytic genes contain an RE1 motif and are direct REST targets. Together, our data support a central role for REST during neural development in promoting NS/P cell self-renewal while restricting the generation and maturation of neurons and oligodendrocytes.

GPR105 ablation prevents inflammation and improves insulin sensitivity in mice with diet-induced obesity

GPR105, a G protein-coupled receptor for UDP-glucose, is highly expressed in several human tissues and participates in the innate immune response. Because inflammation has been implicated as a key initial trigger for type 2 diabetes, we hypothesized that GPR105 (official gene name: P2RY14) might play a role in the initiation of inflammation and insulin resistance in obesity. To this end, we investigated glucose metabolism in GPR105 knockout (KO) and wild-type (WT) mice fed a high-fat diet (HFD). We also examined whether GPR105 regulates macrophage recruitment to liver or adipose tissues by in vivo monocyte tracking and in vitro chemotaxis experiments, followed by transplantation of bone marrow from either KO or WT donors to WT recipients. Our data show that genetic deletion of GPR105 confers protection against HFD-induced insulin resistance, with reduced macrophage infiltration and inflammation in liver, and increased insulin-stimulated Akt phosphorylation in liver, muscle, and adipose tissue. By tracking monocytes from either KO or WT donors, we found that fewer KO monocytes were recruited to the liver of WT recipients. Furthermore, we observed that uridine 5-diphosphoglucose enhanced the in vitro migration of bone marrow-derived macrophages from WT but not KO mice, and that plasma uridine 5-diphosphoglucose levels were significantly higher in obese versus lean mice. Finally, we confirmed that insulin sensitivity improved in HFD mice with a myeloid cell-specific deletion of GPR105. These studies indicate that GPR105 ablation mitigates HFD-induced insulin resistance by inhibiting macrophage recruitment and tissue inflammation. Hence GPR105 provides a novel link between innate immunity and metabolism.

Single-neuron RNA-Seq: technical feasibility and reproducibility

Understanding brain function involves improved knowledge about how the genome specifies such a large diversity of neuronal types. Transcriptome analysis of single neurons has been previously described using gene expression microarrays. Using high-throughput transcriptome sequencing (RNA-Seq), we have developed a method to perform single-neuron RNA-Seq. Following electrophysiology recording from an individual neuron, total RNA was extracted by aspirating the cellular contents into a fine glass electrode tip. The mRNAs were reverse transcribed and amplified to construct a single-neuron cDNA library, and subsequently subjected to high-throughput sequencing. This approach was applied to both individual neurons cultured from embryonic mouse hippocampus, as well as neocortical neurons from live brain slices. We found that the average pairwise Spearman’s rank correlation coefficient of gene expression level expressed as RPKM (reads per kilobase of transcript per million mapped reads) was 0.51 between five cultured neuronal cells, whereas the same measure between three cortical layer 5 neurons in situ was 0.25. The data suggest that there may be greater heterogeneity of the cortical neurons, as compared to neurons in vitro. The results demonstrate the technical feasibility and reproducibility of RNA-Seq in capturing a part of the transcriptome landscape of single neurons, and confirmed that morphologically identical neurons, even from the same region, have distinct gene expression patterns.

Implication of inflammatory macrophages, nuclear receptors, and interferon regulatory factors in increased virulence of pandemic 2009 H1N1 influenza A virus after host adaptation

While pandemic 2009 H1N1 influenza A viruses were responsible for numerous severe infections in humans, these viruses do not typically cause corresponding severe disease in mammalian models. However, the generation of a virulent 2009 H1N1 virus following serial lung passage in mice has allowed for the modeling of human lung pathology in this species. Genetic determinants of mouse-adapted 2009 H1N1 viral pathogenicity have been identified, but the molecular and signaling characteristics of the host response following infection with this adapted virus have not been described. Here we compared the gene expression response following infection of mice with A/CA/04/2009 (CA/04) or the virulent mouse-adapted strain (MA-CA/04). Microarray analysis revealed that increased pathogenicity of MA-CA/04 was associated with the following: (i) an early and sustained inflammatory and interferon response that could be driven in part by interferon regulatory factors (IRFs) and increased NF-κB activation, as well as inhibition of the negative regulator TRIM24, (ii) early and persistent infiltration of immune cells, including inflammatory macrophages, and (iii) the absence of activation of lipid metabolism later in infection, which may be mediated by inhibition of nuclear receptors, including PPARG and HNF1A and -4A, with proinflammatory consequences. Further investigation of these signatures in the host response to other H1N1 viruses of various pathogenicities confirmed their general relevance for virulence of influenza virus and suggested that lung response to MA-CA/04 virus was similar to that following infection with lethal H1N1 r1918 influenza virus. This study links differential activation of IRFs, nuclear receptors, and macrophage infiltration with influenza virulence in vivo.

Pregnancy amelioration of arthritis in SKG mice corresponds with alterations in serum amyloid A3 levels

OBJECTIVES: Pregnancy leads to rheumatoid arthritis remission in humans. The objective of this study was to determine if the SKG mouse could serve as a model for pregnancy-associated inflammatory arthritis amelioration. In addition, the maternal peripheral blood mononuclear cell (PBMC) transcriptome was assessed to define a biomarker associated with remission. METHODS: Cohorts of zymosan-treated pregnant SKG mice and controls were monitored for arthritis progression. Microarray analysis evaluated alterations in gene expression in maternal PBMCs at embryonic day 14.5 (E14.5) between arthritic and pregnancy-remitted mice. A selected target, serum amyloid A3 (SAA3), was further investigated using quantitative reverse transcriptase PCR (qRT-PCR) and an enzyme-linked immunosorbent assay (ELISA). RESULTS: Pregnancy resulted in complete or partial remission in the majority of the zymosan-treated SKG mice. Twenty-seven transcripts were differentially expressed in the PBMCs between arthritic and pregnancy-remitted mice. Expression and plasma SAA3 levels decreased with pregnancy-induced arthritis amelioration and plasma SAA3 levels correlated with arthritis severity. CONCLUSIONS: These results establish the SKG mouse as a model system to study pregnancy-induced amelioration of arthritis. These studies also establish SAA3 as a biomarker of arthritis amelioration in SKG mice. This model can be used to elucidate the molecular and cellular mechanisms underlying the impact of pregnancy on the maternal immune system that results in arthritis amelioration.

Zbtb46 expression distinguishes classical dendritic cells and their committed progenitors from other immune lineages

The zinc finger transcription factor Zbtb46 specifically marks cDCs and their committed precursors and, when overexpressed in BM progenitors, promotes cDC development at the expense of granulocytes.

Distinguishing dendritic cells (DCs) from other cells of the mononuclear phagocyte system is complicated by the shared expression of cell surface markers such as CD11c. In this study, we identified Zbtb46 (BTBD4) as a transcription factor selectively expressed by classical DCs (cDCs) and their committed progenitors but not by plasmacytoid DCs (pDCs), monocytes, macrophages, or other lymphoid or myeloid lineages. Using homologous recombination, we replaced the first coding exon of Zbtb46 with GFP to inactivate the locus while allowing detection of Zbtb46 expression. GFP expression in Zbtb46^gfp/+ mice recapitulated the cDC-specific expression of the native locus, being restricted to cDC precursors (pre-cDCs) and lymphoid organ– and tissue-resident cDCs. GFP⁺ pre-cDCs had restricted developmental potential, generating cDCs but not pDCs, monocytes, or macrophages. Outside the immune system, Zbtb46 was expressed in committed erythroid progenitors and endothelial cell populations. Zbtb46 overexpression in bone marrow progenitor cells inhibited granulocyte potential and promoted cDC development, and although cDCs developed in Zbtb46^gfp/gfp (Zbtb46 deficient) mice, they maintained expression of granulocyte colony-stimulating factor and leukemia inhibitory factor receptors, which are normally down-regulated in cDCs. Thus, Zbtb46 may help enforce cDC identity by restricting responsiveness to non-DC growth factors and may serve as a useful marker to identify rare cDC progenitors and distinguish between cDCs and other mononuclear phagocyte lineages.

Aspirin metabolites are GPR35 agonists

Aspirin is widely used as an anti-inflammatory, anti-platelet, anti-pyretic, and cancer-preventive agent; however, the molecular mode of action is unlikely due entirely to the inhibition of cyclooxygenases. Here, we report the agonist activity of several aspirin metabolites at GPR35, a poorly characterized orphan G protein-coupled receptor. 2,3,5-Trihydroxybenzoic acid, an aspirin catabolite, was found to be the most potent GPR35 agonist among aspirin metabolites. Salicyluric acid, the main metabolite of aspirin, was also active. These results suggest that the GPR35 agonist activity of certain aspirin metabolites may contribute to the clinical features of aspirin.

Conservation and divergence in Toll-like receptor 4-regulated gene expression in primary human versus mouse macrophages

Evolutionary change in gene expression is generally considered to be a major driver of phenotypic differences between species. We investigated innate immune diversification by analyzing interspecies differences in the transcriptional responses of primary human and mouse macrophages to the Toll-like receptor (TLR)-4 agonist lipopolysaccharide (LPS). By using a custom platform permitting cross-species interrogation coupled with deep sequencing of mRNA 5' ends, we identified extensive divergence in LPS-regulated orthologous gene expression between humans and mice (24% of orthologues were identified as "divergently regulated"). We further demonstrate concordant regulation of human-specific LPS target genes in primary pig macrophages. Divergently regulated orthologues were enriched for genes encoding cellular "inputs" such as cell surface receptors (e.g., TLR6, IL-7Rα) and functional "outputs" such as inflammatory cytokines/chemokines (e.g., CCL20, CXCL13). Conversely, intracellular signaling components linking inputs to outputs were typically concordantly regulated. Functional consequences of divergent gene regulation were confirmed by showing LPS pretreatment boosts subsequent TLR6 responses in mouse but not human macrophages, in keeping with mouse-specific TLR6 induction. Divergently regulated genes were associated with a large dynamic range of gene expression, and specific promoter architectural features (TATA box enrichment, CpG island depletion). Surprisingly, regulatory divergence was also associated with enhanced interspecies promoter conservation. Thus, the genes controlled by complex, highly conserved promoters that facilitate dynamic regulation are also the most susceptible to evolutionary change.

Major shifts in the spatio-temporal distribution of lung antioxidant enzymes during influenza pneumonia

With the incessant challenge of exposure to the air we breathe, lung tissue suffers the highest levels of oxygen tension and thus requires robust antioxidant defenses. Furthermore, following injury or infection, lung tissue faces the additional challenge of inflammation-induced reactive oxygen and nitrogen species (ROS/RNS). Little is known about the identity or distribution of lung antioxidant enzymes under normal conditions or during infection-induced inflammation. Using a mouse model of influenza (H1N1 influenza virus A/PR/8/34 [PR8]) in combination with bioinformatics, we identified seven lung-abundant antioxidant enzymes: Glutathione peroxidase 3 (Gpx3), Superoxide dismutase 3 (Sod3), Transferrin (Tf), peroxyredoxin6 (Prdx6), glutathione S-transferase kappa 1 (Gstk1), Catalase (Cat), and Glutathione peroxidase 8 (Gpx8). Interestingly, despite the demand for antioxidants during inflammation, influenza caused depletion in two key antioxidants: Cat and Prdx6. As Cat is highly expressed in Clara cells, virus-induced Clara cell loss contributes to the depletion in Cat. Prdx6 is also reduced due to Clara cell loss, however there is a coincident increase in Prdx6 levels in the alveoli, resulting in only a subtle reduction of Prdx6 overall. Analogously, Gpx3 shifts from the basement membranes underlying the bronchioles and blood vessels to the alveoli, thus maintaining balanced expression. Taken together, these studies identify key lung antioxidants and reveal their distribution among specific cell types. Furthermore, results show that influenza depletes key antioxidants, and that in some cases there is coincident increased expression, consistent with compensatory expression. Given that oxidative stress is known to be a key risk factor during influenza infection, knowledge about the antioxidant repertoire of lungs, and the spatio-temporal distribution of antioxidants, contributes to our understanding of the underlying mechanisms of influenza-induced morbidity and mortality.

A role for epha2 in cell migration and refractive organization of the ocular lens

PURPOSE

The Epha2 receptor is a surprisingly abundant component of the membrane proteome of vertebrate lenses. In humans, genetic studies have linked mutations in EPHA2 to inherited and age-related forms of cataract, but the function of Epha2 in the lens is obscure. To gain insights into the role of Epha2, a comparative analysis of lenses from wild-type and Epha2(-/-) mice was performed.

METHODS

Epha2 distribution was examined using immunocytochemistry and Western blot analysis. Lens optical quality was assessed by laser refractometry. Confocal microscopy was used to analyze cellular phenotypes.

RESULTS

In wild-type lenses, Epha2 was expressed by lens epithelial cells and elongating fibers but was degraded during the later stages of fiber differentiation. Epha2-null lenses retained their transparency, but two key optical parameters, lens shape and internal composition, were compromised in Epha2(-/-) animals. Epha2-null lenses were smaller and more spherical than age-matched wild-type lenses, and laser refractometry revealed a significant decrease in refractive power of the outer cell layers of mutant lenses. In the absence of Epha2, fiber cells deviated from their normal course and terminated at sutures that were no longer centered on the optical axis. Patterning defects were also noted at the level of individual cells. Wild-type fiber cells had hexagonal cross-sectional profiles with membrane protrusions extending from the cell vertices. In contrast, Epha2(-/-) cells had irregular profiles, and protrusions extended from all membrane surfaces.

CONCLUSIONS

These studies indicate that Epha2 is not required for transparency but does play an indispensable role in the cytoarchitecture and refractive quality of the lens.

Role of runt-related transcription factor 3 (RUNX3) in transcription regulation of natural cytotoxicity receptor 1 (NCR1/NKp46), an activating natural killer (NK) cell receptor

Natural cytotoxicity receptor 1 (NCR1), also known as NKp46, is a natural killer (NK) lymphocyte-activating receptor. It is involved in major aspects of NK immune function and shows a high degree of lineage specificity in blood and bone marrow. The nature of its NK-restricted expression is not well understood. In this study, we confirm that human NCR1 NK-specific expression is achieved at the mRNA level. We found two key cis-regulatory elements in the immediate vicinity upstream of the gene. One element acts as an essential promoter, whereas the other acts as a tissue-dependent enhancer/repressor. This latter regulatory element contains a runt related-transcription factor (RUNX) recognition motif that preferentially binds RUNX3. Interfering with RUNX proteins using a dominant negative form results in decreased Ncr1 expression. RUNX3 overexpression had the opposite effect. These findings shed light on the role of RUNX3 in the control of an important NK-activating receptor.

Inflammatory changes in adipose tissue enhance expression of GPR84, a medium-chain fatty acid receptor: TNFα enhances GPR84 expression in adipocytes

In this study we aimed to identify the physiological roles of G protein-coupled receptor 84 (GPR84) in adipose tissue, together with medium-chain fatty acids (MCFAs), the specific ligands for GPR84. In mice, high-fat diet up-regulated GPR84 expression in fat pads. In 3T3-L1 adipocytes, co-culture with a macrophage cell line, RAW264, or TNFα remarkably enhanced GPR84 expression. In the presence of TNFα, MCFAs down-regulated adiponectin mRNA expression in 3T3-L1 adipocytes. Taken together, our results suggest that GPR84 emerges in adipocytes in response to TNFα from infiltrating macrophages and exacerbates the vicious cycle between adiposity and diabesity.

Neuron-specific ELAV/Hu proteins suppress HuR mRNA during neuronal differentiation by alternative polyadenylation

The ubiquitously expressed RNA-binding protein HuR increases the stability and translation of mRNAs encoding growth regulatory proteins that promote proliferation in a variety of cell types. However, the three neuron-specific ELAV/Hu proteins, HuB, HuC and HuD, while binding to the same types of mRNAs, are required instead for neuronal differentiation, and it becomes difficult to reconcile these contrary functions when all four Hu proteins are expressed in the same neuron. HuR mRNA exists as three alternatively polyadenylated variants, a 1.5-kb testes-specific mRNA isoform, a ubiquitous 2.4-kb isoform and a 6.0-kb isoform that we now show is induced during neuronal differentiation and appears to be neuron-specific. This 6.0-kb neuron-specific mRNA isoform is inherently less stable and produces less HuR protein than the ubiquitous 2.4-kb mRNA. Furthermore, we show that neuronal HuB, HuC and HuD, as well as HuR itself, can bind at the 2.4-kb mRNA polyadenylation site, and when overexpressed can affect alternative polyadenylation to generate an extended HuR 3'-UTR that is translationally suppressed. We propose that the regulation of HuR protein expression by alternative polyadenylation allows neurons to post-transcriptionally regulate mRNAs-encoding factors required for proliferation versus differentiation to facilitate neuronal differentiation.

clusterMaker: a multi-algorithm clustering plugin for Cytoscape

Background

In the post-genomic era, the rapid increase in high-throughput data calls for computational tools capable of integrating data of diverse types and facilitating recognition of biologically meaningful patterns within them. For example, protein-protein interaction data sets have been clustered to identify stable complexes, but scientists lack easily accessible tools to facilitate combined analyses of multiple data sets from different types of experiments. Here we present clusterMaker, a Cytoscape plugin that implements several clustering algorithms and provides network, dendrogram, and heat map views of the results. The Cytoscape network is linked to all of the other views, so that a selection in one is immediately reflected in the others. clusterMaker is the first Cytoscape plugin to implement such a wide variety of clustering algorithms and visualizations, including the only implementations of hierarchical clustering, dendrogram plus heat map visualization (tree view), k-means, k-medoid, SCPS, AutoSOME, and native (Java) MCL.

Results

Results are presented in the form of three scenarios of use: analysis of protein expression data using a recently published mouse interactome and a mouse microarray data set of nearly one hundred diverse cell/tissue types; the identification of protein complexes in the yeast Saccharomyces cerevisiae; and the cluster analysis of the vicinal oxygen chelate (VOC) enzyme superfamily. For scenario one, we explore functionally enriched mouse interactomes specific to particular cellular phenotypes and apply fuzzy clustering. For scenario two, we explore the prefoldin complex in detail using both physical and genetic interaction clusters. For scenario three, we explore the possible annotation of a protein as a methylmalonyl-CoA epimerase within the VOC superfamily. Cytoscape session files for all three scenarios are provided in the Additional Files section.

Conclusions

The Cytoscape plugin clusterMaker provides a number of clustering algorithms and visualizations that can be used independently or in combination for analysis and visualization of biological data sets, and for confirming or generating hypotheses about biological function. Several of these visualizations and algorithms are only available to Cytoscape users through the clusterMaker plugin. clusterMaker is available via the Cytoscape plugin manager.

Further analysis of the lens phenotype in Lim2-deficient mice

PURPOSE

Lim2 (MP20) is the second most abundant integral protein of lens fiber cell membranes. A comparative analysis was performed of wild-type and Lim2-deficient (Lim2(Gt/Gt)) mouse lenses, to better define the anatomic and physiologic roles of Lim2.

METHODS

Scanning electron microscopy (SEM) and confocal microscopy were used to assess the contribution of Lim2 to lens tissue architecture. Differentiation-dependent changes in cytoskeletal composition were identified by mass spectrometry and immunoblot analysis. The effects on cell-cell communication were quantified using impedance analysis.

RESULTS

Lim2-null lenses were grossly normal. At the cellular level, however, subtle structural alterations were evident. Confocal microscopy and SEM analysis revealed that cortical Lim2(Gt/Gt) fiber cells lacked the undulating morphology that characterized wild-type fiber cells. On SDS-PAGE analysis the composition of cortical fiber cells from wild-type and Lim2-null lenses appeared similar. However, marked disparities were evident in samples prepared from the lens core of the two genotypes. Several cytoskeletal proteins that were abundant in wild-type core fiber cells were diminished in the cores of Lim2(Gt/Gt) lenses. Electrophysiological measurements indicated a small decrease in the membrane potential of Lim2(Gt/Gt) lenses and a two-fold increase in the effective intracellular resistivity. In the lens core, this may have reflected decreased expression levels of the gap junction protein connexin 46 (Cx46). In contrast, increased resistivity in the outer cell layers of Lim2(Gt/Gt) lenses could not be attributed to decreased connexin expression and may reflect the absence of cell fusions in Lim2(Gt/Gt) lenses.

CONCLUSIONS

Comparative analysis of wild-type and Lim2-deficient lenses has implicated Lim2 in maintenance of cytoskeletal integrity, cell morphology, and intercellular communication.

Identification of quantitative trait loci influencing skeletal architecture in mice: emergence of Cdh11 as a primary candidate gene regulating femoral morphology

Bone strength is influenced by many properties intrinsic to bone, including its mass, geometry, and mineralization. To further advance our understanding of the genetic basis of bone-strength-related traits, we used a large (n = 815), moderately (G(4) ) advanced intercross line (AIL) of mice derived from a high-runner selection line (HR) and the C57BL/6J inbred strain. In total, 16 quantitative trait loci (QTLs) were identified that affected areal bone mineral density (aBMD) and femoral length and width. Four significant (p < .05) and one suggestive (p < .10) QTLs were identified for three aBMD measurements: total body, vertebral, and femoral. A QTL on chromosome (Chr.) 3 influenced all three aBMD measures, whereas the other four QTLs were unique to a single measure. A total of 10 significant and one suggestive QTLs were identified for femoral length (FL) and two measures of femoral width, anteroposterior (AP) and mediolateral (ML). FL QTLs were distinct from loci affecting AP and ML width, and of the 7 AP QTLs, only three affected ML. A QTL on Chr. 8 that explained 7.1% and 4.0% of the variance in AP and ML, respectively, was mapped to a 6-Mb region harboring 12 protein-coding genes. The pattern of haplotype diversity across the QTL region and expression profiles of QTL genes suggested that of the 12, cadherin 11 (Cdh11) was most likely the causal gene. These findings, when combined with existing data from gene knockouts, identify Cdh11 as a strong candidate gene within which genetic variation may affect bone morphology.