Search

Osteocytes are master regulators of the skeleton. We mapped the transcriptome of osteocytes from different skeletal sites, across age and sexes in mice to reveal genes and molecular programs that control this complex cellular-network. We define an osteocyte transcriptome signature of 1239 genes that distinguishes osteocytes from other cells. 77% have no previously known role in the skeleton and are enriched for genes regulating neuronal network formation, suggesting this programme is important in osteocyte communication. We evaluated 19 skeletal parameters in 733 knockout mouse lines and reveal 26 osteocyte transcriptome signature genes that control bone structure and function. We showed osteocyte transcriptome signature genes are enriched for human orthologs that cause monogenic skeletal disorders (P = 2.4 × 10 -22 ) and are associated with the polygenic diseases osteoporosis (P = 1.8 × 10 -13 ) and osteoarthritis (P = 1.6 × 10 -7 ). Thus, we reveal the molecular landscape that regulates osteocyte network formation and function and establish the importance of osteocytes in human skeletal disease.

Osteoclasts are large multinucleated bone-resorbing cells formed by the fusion of monocyte/macrophage-derived precursors that are thought to undergo apoptosis once resorption is complete. Here, by intravital imaging, we reveal that RANKL-stimulated osteoclasts have an alternative cell fate in which they fission into daughter cells called osteomorphs. Inhibiting RANKL blocked this cellular recycling and resulted in osteomorph accumulation. Single-cell RNA sequencing showed that osteomorphs are transcriptionally distinct from osteoclasts and macrophages and express a number of non-canonical osteoclast genes that are associated with structural and functional bone phenotypes when deleted in mice. Furthermore, genetic variation in human orthologs of osteomorph genes causes monogenic skeletal disorders and associates with bone mineral density, a polygenetic skeletal trait. Thus, osteoclasts recycle via osteomorphs, a cell type involved in the regulation of bone resorption that may be targeted for the treatment of skeletal diseases., Competing Interests: Declaration of interests F.H.E. and S.S. are employees and shareholders of BioVinc LLC, who has licensed patents for fluorescent probe compositions used in this work. P.I.C. has grant funding from Amgen. Other authors have no competing financial interests., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

Osteoarthritis causes pain and functional disability for over 500 million people worldwide. To develop disease-stratifying tools and modifying therapies, we need a better understanding of the molecular basis of the disease in relevant tissue and cell types. Here, we study primary cartilage and synovium from 115 patients with osteoarthritis to construct a deep molecular signature map of the disease. By integrating genetics with transcriptomics and proteomics, we discover molecular trait loci in each tissue type and omics level, identify likely effector genes for osteoarthritis-associated genetic signals and highlight high-value targets for drug development and repurposing. These findings provide insights into disease aetiopathology, and offer translational opportunities in response to the global clinical challenge of osteoarthritis.

Osteoarthritis causes debilitating pain and disability, resulting in a considerable socioeconomic burden, yet no drugs are available that prevent disease onset or progression. Here, we develop, validate and use rapid-throughput imaging techniques to identify abnormal joint phenotypes in randomly selected mutant mice generated by the International Knockout Mouse Consortium. We identify 14 genes with functional involvement in osteoarthritis pathogenesis, including the homeobox gene Pitx1, and functionally characterize 6 candidate human osteoarthritis genes in mouse models. We demonstrate sensitivity of the methods by identifying age-related degenerative joint damage in wild-type mice. Finally, we phenotype previously generated mutant mice with an osteoarthritis-associated polymorphism in the Dio2 gene by CRISPR/Cas9 genome editing and demonstrate a protective role in disease onset with public health implications. We hope this expanding resource of mutant mice will accelerate functional gene discovery in osteoarthritis and offer drug discovery opportunities for this common, incapacitating chronic disease.

Context: The X-linked immunoglobulin superfamily, member 1 (IGSF1), gene is highly expressed in the hypothalamus and in pituitary cells of the POU1F1 lineage. Human loss-of-function mutations in IGSF1 cause central hypothyroidism, hypoprolactinemia, and macroorchidism. Additionally, most affected adults exhibit higher than average IGF-1 levels and anecdotal reports describe acromegaloid features in older subjects. However, somatotrope function has not yet been formally evaluated in this condition., Objective: We aimed to evaluate the role of IGSF1 in human and murine somatotrope function., Patients, Design, and Setting: We evaluated 21 adult males harboring hemizygous IGSF1 loss-of-function mutations for features of GH excess, in an academic clinical setting., Main Outcome Measures: We compared biochemical and tissue markers of GH excess in patients and controls, including 24-hour GH profile studies in 7 patients. Parallel studies were undertaken in male Igsf1-deficient mice and wild-type littermates., Results: IGSF1-deficient adult male patients demonstrated acromegaloid facial features with increased head circumference as well as increased finger soft-tissue thickness. Median serum IGF-1 concentrations were elevated, and 24-hour GH profile studies confirmed 2- to 3-fold increased median basal, pulsatile, and total GH secretion. Male Igsf1-deficient mice also demonstrated features of GH excess with increased lean mass, organ size, and skeletal dimensions and elevated mean circulating IGF-1 and pituitary GH levels., Conclusions: We demonstrate somatotrope neurosecretory hyperfunction in IGSF1-deficient humans and mice. These observations define a hitherto uncharacterized role for IGSF1 in somatotropes and indicate that patients with IGSF1 mutations should be evaluated for long-term consequences of increased GH exposure., (© Endocrine Society 2019.)

Objective: Bone loss in anorexia nervosa and following bariatric surgery is associated with an elevated circulating concentration of the gastrointestinal, anorexigenic hormone, peptide YY (PYY). Selective deletion of the PYY receptor Y1R in osteoblasts or Y2R in the hypothalamus results in high bone mass, but deletion of PYY in mice has resulted in conflicting skeletal phenotypes leading to uncertainty regarding its role in the regulation of bone mass. As PYY analogs are under development for treatment of obesity, we aimed to clarify the relationship between PYY and bone mass., Methods: The skeletal phenotype of Pyy knockout (KO) mice was investigated during growth (postnatal day P14) and adulthood (P70 and P186) using X-ray microradiography, micro-CT, back-scattered electron scanning electron microscopy (BSE-SEM), histomorphometry and biomechanical testing., Results: Bones from juvenile and Pyy KO mice were longer (P < 0.001), with decreased bone mineral content (P < 0.001). Whereas, bones from adult Pyy KO mice had increased bone mineral content (P < 0.05) with increased mineralisation of both cortical (P < 0.001) and trabecular (P < 0.001) compartments. Long bones from adult Pyy KO mice were stronger (maximum load P < 0.001), with increased stiffness (P < 0.01) and toughness (P < 0.05) compared to wild-type (WT) control mice despite increased cortical vascularity and porosity (P < 0.001). The increased bone mass and strength in Pyy KO mice resulted from increases in trabecular (P < 0.01) and cortical bone formation (P < 0.05)., Conclusions: These findings demonstrate that PYY acts as a negative regulator of osteoblastic bone formation, implicating increased PYY levels in the pathogenesis of bone loss during anorexia or following bariatric surgery., (Copyright © 2019 The Authors. Published by Elsevier Inc. All rights reserved.)

Osteoporosis is characterized by low bone mineral density (BMD) and fragility fracture and affects over 200 million people worldwide. Bone quality describes the material properties that contribute to strength independently of BMD, and its quantitative analysis is a major priority in osteoporosis research. Tissue mineralization is a fundamental process requiring calcium and phosphate transporters. Here we identify impaired bone quality and strength in Slc20a2 -/- mice lacking the phosphate transporter SLC20A2. Juveniles had abnormal endochondral and intramembranous ossification, decreased mineral accrual, and short stature. Adults exhibited only small reductions in bone mass and mineralization but a profound impairment of bone strength. Bone quality was severely impaired in Slc20a2 -/- mice: yield load (-2.3 SD), maximum load (-1.7 SD), and stiffness (-2.7 SD) were all below values predicted from their bone mineral content as determined in a cohort of 320 wild-type controls. These studies identify Slc20a2 as a physiological regulator of tissue mineralization and highlight its critical role in the determination of bone quality and strength. © 2019 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals Inc., (© 2019 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals Inc.)

Background: Thyroid hormones act in bone and cartilage via thyroid hormone receptor alpha (TRα). In the absence of triiodothyronine (T3), TRα interacts with co-repressors, including nuclear receptor co-repressor-1 (NCoR1), which recruit histone deacetylases (HDACs) and mediate transcriptional repression. Dominant-negative mutations of TRα cause resistance to thyroid hormone alpha (RTHα; OMIM 614450), characterized by excessive repression of T3 target genes leading to delayed skeletal development, growth retardation, and bone dysplasia. Treatment with thyroxine has been of limited benefit, even in mildly affected individuals, and there is a need for new therapeutic strategies. It was hypothesized that (i) the skeletal manifestations of RTHα are mediated by the persistent TRα/NCoR1/HDAC repressor complex containing mutant TRα, and (ii) treatment with the HDAC inhibitor suberoylanilide hydroxamic acid (SAHA) would ameliorate these manifestations. Methods: The skeletal phenotypes of (i) Thra1 PV/+ mice, a well characterized model of RTHα; (ii) Ncor1 ΔID/ΔID mice, which express an NCoR1 mutant that fails to interact with TRα; and (iii) Thra1 PV/+ Ncor1 ΔID/ΔID double-mutant adult mice were determined. Wild-type, Thra1 PV/+ , Ncor1 ΔID/ΔID , and Thra1 PV/+ Ncor1 ΔID/ΔID double-mutant mice were also treated with SAHA to determine whether HDAC inhibition results in amelioration of skeletal abnormalities. Results: Thra1 PV/+ mice had a severe skeletal dysplasia, characterized by short stature, abnormal bone morphology, and increased bone mineral content. Despite normal bone length, Ncor1 ΔID/ΔID mice displayed increased cortical bone mass, mineralization, and strength. Thra1 PV/+ Ncor1 ΔID/ΔID double-mutant mice displayed only a small improvement of skeletal abnormalities compared to Thra1 PV/+ mice. Treatment with SAHA to inhibit histone deacetylation had no beneficial or detrimental effects on bone structure, mineralization, or strength in wild-type or mutant mice. Conclusions: These studies indicate treatment with SAHA is unlikely to improve the skeletal manifestations of RTHα. Nevertheless, the findings (i) confirm that TRα1 has a critical role in the regulation of skeletal development and adult bone mass, (ii) suggest a physiological role for alternative co-repressors that interact with TR in skeletal cells, and (iii) demonstrate a novel role for NCoR1 in the regulation of adult bone mass and strength.

In the version of this article initially published, in Fig. 5a, the data in the right column of 'DAAM2 gRNA1' were incorrectly plotted as circles indicating 'untreated' rather than as squares indicating 'treated'. The error has been corrected in the HTML and PDF versions of the article.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10 -75 ) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

Levothyroxine (LT4) is a form of thyroid hormone used to treat hypothyroidism. In the brain, T4 is converted to the active form T3 by type 2 deiodinase (D2). Thus, it is intriguing that carriers of the Thr92Ala polymorphism in the D2 gene (DIO2) exhibit clinical improvement when liothyronine (LT3) is added to LT4 therapy. Here, we report that D2 is a cargo protein in ER Golgi intermediary compartment (ERGIC) vesicles, recycling between ER and Golgi. The Thr92-to-Ala substitution (Ala92-D2) caused ER stress and activated the unfolded protein response (UPR). Ala92-D2 accumulated in the trans-Golgi and generated less T3, which was restored by eliminating ER stress with the chemical chaperone 4-phenyl butyric acid (4-PBA). An Ala92-Dio2 polymorphism-carrying mouse exhibited UPR and hypothyroidism in distinct brain areas. The mouse refrained from physical activity, slept more, and required additional time to memorize objects. Enhancing T3 signaling in the brain with LT3 improved cognition, whereas restoring proteostasis with 4-PBA eliminated the Ala92-Dio2 phenotype. In contrast, primary hypothyroidism intensified the Ala92-Dio2 phenotype, with only partial response to LT4 therapy. Disruption of cellular proteostasis and reduced Ala92-D2 activity may explain the failure of LT4 therapy in carriers of Thr92Ala-DIO2.

This chapter describes the use of point projection digital microradiography for rapid imaging and quantitation of bone mineral content in mice.

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

The shapes of homologous skeletal elements in the vertebrate forelimb and hindlimb are distinct, with each element exquisitely adapted to their divergent functions. Many of the signals and signalling pathways responsible for patterning the developing limb bud are common to both forelimb and hindlimb. How disparate morphologies are generated from common signalling inputs during limb development remains poorly understood. We show that, similar to what has been shown in the chick, characteristic differences in mouse forelimb and hindlimb cartilage morphology are maintained when chondrogenesis proceeds in vitro away from the endogenous limb bud environment. Chondrogenic nodules that form in high-density micromass cultures derived from forelimb and hindlimb buds are consistently different in size and shape. We described analytical tools we have developed to quantify these differences in nodule morphology and demonstrate that characteristic hindlimb nodule morphology is lost in the absence of the hindlimb-restricted limb modifier gene Pitx1. Furthermore, we show that ectopic expression of Pitx1 in the forelimb is sufficient to generate nodule patterns characteristic of the hindlimb. We also demonstrate that hindlimb cells are less adhesive to the tissue culture substrate and, within the limb environment, to the extracellular matrix and to each other. These results reveal autonomously programmed differences in forelimb and hindlimb cartilage precursors of the limb skeleton are controlled, at least in part, by Pitx1 and suggest this has an important role in generating distinct limb-type morphologies. Our results demonstrate that the micromass culture system is ideally suited to study cues governing morphogenesis of limb skeletal elements in a simple and experimentally tractable in vitro system that reflects in vivo potential.

The forelimbs and hindlimbs of vertebrates are bilaterally symmetric. The mechanisms that ensure symmetric limb formation are unknown but they can be disrupted in disease. In Holt-Oram Syndrome (HOS), caused by mutations in TBX5, affected individuals have left-biased upper/forelimb defects. We demonstrate a role for the transcription factor Tbx5 in ensuring the symmetric formation of the left and right forelimb. In our mouse model, bilateral hypomorphic levels of Tbx5 produces asymmetric forelimb defects that are consistently more severe in the left limb than the right, phenocopying the left-biased limb defects seen in HOS patients. In Tbx hypomorphic mutants maintained on an INV mutant background, with situs inversus, the laterality of defects is reversed. Our data demonstrate an early, inherent asymmetry in the left and right limb-forming regions and that threshold levels of Tbx5 are required to overcome this asymmetry to ensure symmetric forelimb formation., Competing Interests: The authors have declared that no competing interests exist.

Mutations in components of the intraflagellar transport (IFT) machinery required for assembly and function of the primary cilium cause a subset of human ciliopathies characterized primarily by skeletal dysplasia. Recently, mutations in the IFT-A gene IFT144 have been described in patients with Sensenbrenner and Jeune syndromes, which are associated with short ribs and limbs, polydactyly and craniofacial defects. Here, we describe an N-ethyl-N-nitrosourea-derived mouse mutant with a hypomorphic missense mutation in the Ift144 gene. The mutant twinkle-toes (Ift144(twt)) phenocopies a number of the skeletal and craniofacial anomalies seen in patients with human skeletal ciliopathies. Like other IFT-A mouse mutants, Ift144 mutant embryos display a generalized ligand-independent expansion of hedgehog (Hh) signalling, in spite of defective ciliogenesis and an attenuation of the ability of mutant cells to respond to upstream stimulation of the pathway. This enhanced Hh signalling is consistent with cleft palate and polydactyly phenotypes in the Ift144(twt) mutant, although extensive rib branching, fusion and truncation phenotypes correlate with defects in early somite patterning and may reflect contributions from multiple signalling pathways. Analysis of embryos harbouring a second allele of Ift144 which represents a functional null, revealed a dose-dependent effect on limb outgrowth consistent with the short-limb phenotypes characteristic of these ciliopathies. This allelic series of mouse mutants provides a unique opportunity to uncover the underlying mechanistic basis of this intriguing subset of ciliopathies.

The bones of the vertebrate limb form by the process of endochondral ossification, whereby limb mesenchyme condenses to form an intermediate cartilage scaffold that is then replaced by bone. Although Indian hedgehog (IHH) is known to control hypertophic differentiation of chondrocytes during this process, the role of hedgehog signaling in the earlier stages of chondrogenesis is less clear. We have conditionally inactivated the hedgehog receptor Ptc1 in undifferentiated limb mesenchyme of the mouse limb using Prx1-Cre, thus inducing constitutively active ligand-independent hedgehog signaling. In addition to major patterning defects, we observed a marked disruption to the cartilage elements in the limbs of Prx1-Cre:Ptc1(c/c) embryos. Using an in vitro micromass culture system we show that this defect lies downstream of mesenchymal cell condensation and likely upstream of chondrocyte differentiation. Despite early increases in levels of chondrogenic genes, soon after mesenchymal condensation the stromal layer of Prx1-Cre:Ptc1(c/c)-derived micromass cultures is characterized by a loss of cell integrity, which is associated with increased cell death and a striking decrease in Alcian blue staining cartilage nodules. Furthermore, inhibition of the hedgehog pathway activation using cyclopamine was sufficient to essentially overcome this chondrogenic defect in both micromass and ex vivo explant assays of Prx1-Cre:Ptc1(c/c) limbs. These data demonstrate for the first time the inhibitory effect of cell autonomously activated hedgehog signaling on chondrogenesis, and stress the importance of PTC1 in maintaining strict control of signaling levels during this phase of skeletal development.

The limb has long been considered a paradigm for organogenesis because of its simplicity and ease of manipulation. However, it has become increasingly clear that the processes required to produce a perfectly formed limb involve complex molecular interactions across all three axes of limb development. Old models have evolved with acquisition of molecular knowledge, and in more recent times mathematical modeling approaches have been invoked to explain the precise spatio-temporal regulation of gene networks that coordinate limb patterning and outgrowth. This review focuses on recent advances in our understanding of vertebrate limb development, highlighting the signaling interactions required to lay down the pattern on which the processes of differentiation will act to ultimately produce the final limb., (Copyright 2010 Elsevier Inc. All rights reserved.)

Pitrm1 is a zinc metalloendopeptidase that has been implicated in Alzheimer's disease and mitochondrial peptide degradation, but to date no major role in embryonic development has been documented. In a screen for genes regulated by hedgehog signaling in the mouse limb, we showed that expression of Pitrm1 is upregulated in response to loss of the Gli3 transcription factor. Here we confirm spatial changes in Pitrm1 expression in the Gli3 mutant mouse limb and examine Pitrm1 expression in Shh null and Ptch1 conditional deletion mouse mutants. In wild-type mice, Pitrm1 is expressed in a number of developing tissues known to be patterned by Sonic hedgehog, including the limbs, face, cortex, hippocampus, cerebellum, tectum, sub-mandibular gland, lung, genital tubercle, hair follicles, and the enamel knot of the teeth. Additionally, Pitrm1 is expressed in Pax3-expressing myoblast progenitors in the limb, the dermomyotome, and developing muscles of the face and torso., ((c) 2009 Wiley-Liss, Inc.)

The vertebrate hedgehog receptor patched 1 (Ptc1) is crucial for negative regulation of the sonic hedgehog (Shh) pathway during anterior-posterior patterning of the limb. We have conditionally inactivated Ptc1 in the mesenchyme of the mouse limb using Prx1-Cre. This results in constitutive activation of hedgehog (Hh) signalling during the early stages of limb budding. Our data suggest that variations in the timing and efficiency of Cre-mediated excision result in differential forelimb and hindlimb phenotypes. Hindlimbs display polydactyly (gain of digits) and a molecular profile similar to the Gli3 mutant extra-toes. Strikingly, forelimbs are predominantly oligodactylous (displaying a loss of digits), with a symmetrical, mirror-image molecular profile that is consistent with re-specification of the anterior forelimb to a posterior identity. Our data suggest that this is related to very early inactivation of Ptc1 in the forelimb perturbing the gene regulatory networks responsible for both the pre-patterning and the subsequent patterning stages of limb development. These results establish the importance of the downstream consequences of Hh pathway repression, and identify Ptc1 as a key player in limb patterning even prior to the onset of Shh expression.

The NET/Nlz family of zinc finger transcription factors contribute to aspects of developmental growth and patterning across evolutionarily diverse species. To date, however, these molecules remain largely uncharacterized in mouse and chick. We previously reported that limb bud expression of Zfp503, the mouse orthologue of zebrafish nlz2/znf503, is dependent on Gli3. Here, we show that Zfp503/Znf503 is expressed in a restricted pattern during mouse and chick embryogenesis, with particularly dynamic expression in the developing limbs, face, somites, and brain. We also add to our previous data on Gli3 regulation by showing that the anterior domain of Zfp503 expression in the mouse limb is responsive to genetic and nongenetic manipulation of hedgehog signaling. Finally, we demonstrate that posterior expression of Znf503 in the chick limb is responsive to bone morphogenetic protein (BMP) signaling, indicating that Zfp503/Znf503 may act at the nexus of multiple signaling pathways in development., ((c) 2008 Wiley-Liss, Inc.)

Background: Some years ago we established an N-ethyl-N-nitrosourea screen for modifiers of transgene variegation in the mouse and a preliminary description of the first six mutant lines, named MommeD1-D6, has been published. We have reported the underlying genes in three cases: MommeD1 is a mutation in SMC hinge domain containing 1 (Smchd1), a novel modifier of epigenetic gene silencing; MommeD2 is a mutation in DNA methyltransferase 1 (Dnmt1); and MommeD4 is a mutation in Smarca 5 (Snf2h), a known chromatin remodeler. The identification of Dnmt1 and Smarca5 attest to the effectiveness of the screen design., Results: We have now extended the screen and have identified four new modifiers, MommeD7-D10. Here we show that all ten MommeDs link to unique sites in the genome, that homozygosity for the mutations is associated with severe developmental abnormalities and that heterozygosity results in phenotypic abnormalities and reduced reproductive fitness in some cases. In addition, we have now identified the underlying genes for MommeD5 and MommeD10. MommeD5 is a mutation in Hdac1, which encodes histone deacetylase 1, and MommeD10 is a mutation in Baz1b (also known as Williams syndrome transcription factor), which encodes a transcription factor containing a PHD-type zinc finger and a bromodomain. We show that reduction in the level of Baz1b in the mouse results in craniofacial features reminiscent of Williams syndrome., Conclusions: These results demonstrate the importance of dosage-dependent epigenetic reprogramming in the development of the embryo and the power of the screen to provide mouse models to study this process.

Craniofacial anomalies are a common feature of human congenital dysmorphology syndromes, suggesting that genes expressed in the developing face are likely to play a wider role in embryonic development. To facilitate the identification of genes involved in embryogenesis, we previously constructed an enriched cDNA library by subtracting adult mouse liver cDNA from that of embryonic day (E)10.5 mouse pharyngeal arch cDNA. From this library, 273 unique clones were sequenced and known proteins binned into functional categories in order to assess enrichment of the library (1). We have now selected 31 novel and poorly characterised genes from this library and present bioinformatic analysis to predict proteins encoded by these genes, and to detect evolutionary conservation. Of these genes 61% (19/31) showed restricted expression in the developing embryo, and a subset of these was chosen for further in silico characterisation as well as experimental determination of subcellular localisation based on transient transfection of predicted full-length coding sequences into mammalian cell lines. Where a human orthologue of these genes was detected, chromosomal localisation was determined relative to known loci for human congenital disease.