Whole-genome sequencing identifies EN1 as a determinant of bone density and fracture.

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes.

Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

Compound heterozygous mutations in RIPPLY2 associated with vertebral segmentation defects.

Segmentation defects of the vertebrae (SDV) are caused by aberrant somite formation during embryogenesis and result in irregular formation of the vertebrae and ribs. The Notch signal transduction pathway plays a critical role in somite formation and patterning in model vertebrates. In humans, mutations in several genes involved in the Notch pathway are associated with SDV, with both autosomal recessive (MESP2, DLL3, LFNG, HES7) and autosomal dominant (TBX6) inheritance. However, many individuals with SDV do not carry mutations in these genes. Using whole-exome capture and massive parallel sequencing, we identified compound heterozygous mutations in RIPPLY2 in two brothers with multiple regional SDV, with appropriate familial segregation. One novel mutation (c.A238T:p.Arg80*) introduces a premature stop codon. In transiently transfected C2C12 mouse myoblasts, the RIPPLY2 mutant protein demonstrated impaired transcriptional repression activity compared with wild-type RIPPLY2 despite similar levels of expression. The other mutation (c.240-4T>G), with minor allele frequency <0.002, lies in the highly conserved splice site consensus sequence 5' to the terminal exon. Ripply2 has a well-established role in somitogenesis and vertebral column formation, interacting at both gene and protein levels with SDV-associated Mesp2 and Tbx6. We conclude that compound heterozygous mutations in RIPPLY2 are associated with SDV, a new gene for this condition.

Short-rib polydactyly and Jeune syndromes are caused by mutations in WDR60.

Short-rib polydactyly syndromes (SRPS I-V) are a group of lethal congenital disorders characterized by shortening of the ribs and long bones, polydactyly, and a range of extraskeletal phenotypes. A number of other disorders in this grouping, including Jeune and Ellis-van Creveld syndromes, have an overlapping but generally milder phenotype. Collectively, these short-rib dysplasias (with or without polydactyly) share a common underlying defect in primary cilium function and form a subset of the ciliopathy disease spectrum. By using whole-exome capture and massive parallel sequencing of DNA from an affected Australian individual with SRPS type III, we detected two novel heterozygous mutations in WDR60, a relatively uncharacterized gene. These mutations segregated appropriately in the unaffected parents and another affected family member, confirming compound heterozygosity, and both were predicted to have a damaging effect on the protein. Analysis of an additional 54 skeletal ciliopathy exomes identified compound heterozygous mutations in WDR60 in a Spanish individual with Jeune syndrome of relatively mild presentation. Of note, these two families share one novel WDR60 missense mutation, although haplotype analysis suggested no shared ancestry. We further show that WDR60 localizes at the base of the primary cilium in wild-type human chondrocytes, and analysis of fibroblasts from affected individuals revealed a defect in ciliogenesis and aberrant accumulation of the GLI2 transcription factor at the centrosome or basal body in the absence of an obvious axoneme. These findings show that WDR60 mutations can cause skeletal ciliopathies and suggest a role for WDR60 in ciliogenesis.

Mutations in the gene encoding IFT dynein complex component WDR34 cause Jeune asphyxiating thoracic dystrophy.

Bidirectional (anterograde and retrograde) motor-based intraflagellar transport (IFT) governs cargo transport and delivery processes that are essential for primary cilia growth and maintenance and for hedgehog signaling functions. The IFT dynein-2 motor complex that regulates ciliary retrograde protein transport contains a heavy chain dynein ATPase/motor subunit, DYNC2H1, along with other less well functionally defined subunits. Deficiency of IFT proteins, including DYNC2H1, underlies a spectrum of skeletal ciliopathies. Here, by using exome sequencing and a targeted next-generation sequencing panel, we identified a total of 11 mutations in WDR34 in 9 families with the clinical diagnosis of Jeune syndrome (asphyxiating thoracic dystrophy). WDR34 encodes a WD40 repeat-containing protein orthologous to Chlamydomonas FAP133, a dynein intermediate chain associated with the retrograde intraflagellar transport motor. Three-dimensional protein modeling suggests that the identified mutations all affect residues critical for WDR34 protein-protein interactions. We find that WDR34 concentrates around the centrioles and basal bodies in mammalian cells, also showing axonemal staining. WDR34 coimmunoprecipitates with the dynein-1 light chain DYNLL1 in vitro, and mining of proteomics data suggests that WDR34 could represent a previously unrecognized link between the cytoplasmic dynein-1 and IFT dynein-2 motors. Together, these data show that WDR34 is critical for ciliary functions essential to normal development and survival, most probably as a previously unrecognized component of the mammalian dynein-IFT machinery.

Defects in the IFT-B component IFT172 cause Jeune and Mainzer-Saldino syndromes in humans.

Intraflagellar transport (IFT) depends on two evolutionarily conserved modules, subcomplexes A (IFT-A) and B (IFT-B), to drive ciliary assembly and maintenance. All six IFT-A components and their motor protein, DYNC2H1, have been linked to human skeletal ciliopathies, including asphyxiating thoracic dystrophy (ATD; also known as Jeune syndrome), Sensenbrenner syndrome, and Mainzer-Saldino syndrome (MZSDS). Conversely, the 14 subunits in the IFT-B module, with the exception of IFT80, have unknown roles in human disease. To identify additional IFT-B components defective in ciliopathies, we independently performed different mutation analyses: candidate-based sequencing of all IFT-B-encoding genes in 1,467 individuals with a nephronophthisis-related ciliopathy or whole-exome resequencing in 63 individuals with ATD. We thereby detected biallelic mutations in the IFT-B-encoding gene IFT172 in 12 families. All affected individuals displayed abnormalities of the thorax and/or long bones, as well as renal, hepatic, or retinal involvement, consistent with the diagnosis of ATD or MZSDS. Additionally, cerebellar aplasia or hypoplasia characteristic of Joubert syndrome was present in 2 out of 12 families. Fibroblasts from affected individuals showed disturbed ciliary composition, suggesting alteration of ciliary transport and signaling. Knockdown of ift172 in zebrafish recapitulated the human phenotype and demonstrated a genetic interaction between ift172 and ift80. In summary, we have identified defects in IFT172 as a cause of complex ATD and MZSDS. Our findings link the group of skeletal ciliopathies to an additional IFT-B component, IFT172, similar to what has been shown for IFT-A.

Autosomal dominant spondylocostal dysostosis is caused by mutation in TBX6.

In humans, congenital spinal defects occur with an incidence of 0.5-1 per 1000 live births. One of the most severe syndromes with such defects is spondylocostal dysostosis (SCD). Over the past decade, the genetic basis of several forms of autosomal recessive SCD cases has been solved with the identification of four causative genes (DLL3, MESP2, LFNG and HES7). Autosomal dominant forms of SCD have also been reported, but to date no genetic etiology has been described for these. Here, we have used exome capture and next-generation sequencing to identify a stoploss mutation in TBX6 that segregates with disease in two generations of one family. We show that this mutation has a deleterious effect on the transcriptional activation activity of the TBX6 protein, likely due to haploinsufficiency. In mouse, Tbx6 is essential for the patterning of the vertebral precursor tissues, somites; thus, mutation of TBX6 is likely to be causative of SCD in this family. This is the first identification of the genetic cause of an autosomal dominant form of SCD, and also demonstrates the potential of exome sequencing to identify genetic causes of dominant diseases even in small families with few affected individuals.

Deep-transcriptome and ribonome sequencing redefines the molecular networks of pluripotency and the extracellular space in human embryonic stem cells.

Recent RNA-sequencing studies have shown remarkable complexity in the mammalian transcriptome. The ultimate impact of this complexity on the predicted proteomic output is less well defined. We have undertaken strand-specific RNA sequencing of multiple cellular RNA fractions (>20 Gb) to uncover the transcriptional complexity of human embryonic stem cells (hESCs). We have shown that human embryonic stem (ES) cells display a high degree of transcriptional diversity, with more than half of active genes generating RNAs that differ from conventional gene models. We found evidence that more than 1000 genes express long 5' and/or extended 3'UTRs, which was confirmed by "virtual Northern" analysis. Exhaustive sequencing of the membrane-polysome and cytosolic/untranslated fractions of hESCs was used to identify RNAs encoding peptides destined for secretion and the extracellular space and to demonstrate preferential selection of transcription complexity for translation in vitro. The impact of this newly defined complexity on known gene-centric network models such as the Plurinet and the cell surface signaling machinery in human ES cells revealed a significant expansion of known transcript isoforms at play, many predicting possible alternative functions based on sequence alterations within key functional domains.

Whole exome sequencing is an efficient and sensitive method for detection of germline mutations in patients with phaeochromcytomas and paragangliomas.

BACKGROUND: Genetic testing is recommended when the probability of a disease-associated germline mutation exceeds 10%. Germline mutations are found in approximately 25% of individuals with phaeochromcytoma (PCC) or paraganglioma (PGL); however, genetic heterogeneity for PCC/PGL means many genes may require sequencing. A phenotype-directed iterative approach may limit costs but may also delay diagnosis, and will not detect mutations in genes not previously associated with PCC/PGL. OBJECTIVE: To assess whether whole exome sequencing (WES) was efficient and sensitive for mutation detection in PCC/PGL. METHODS: Whole exome sequencing was performed on blinded samples from eleven individuals with PCC/PGL and known mutations. Illumina TruSeq (Illumina Inc, San Diego, CA, USA) was used for exome capture of seven samples, and NimbleGen SeqCap EZ v3.0 (Roche NimbleGen Inc, Basel, Switzerland) for five samples (one sample was repeated). Massive parallel sequencing was performed on multiplexed samples. Sequencing data were called using Genome Analysis Toolkit and annotated using annovar. Data were assessed for coding variants in RET, NF1, VHL, SDHD, SDHB, SDHC, SDHA, SDHAF2, KIF1B, TMEM127, EGLN1 and MAX. Target capture of five exome capture platforms was compared. RESULTS: Six of seven mutations were detected using Illumina TruSeq exome capture. All five mutations were detected using NimbleGen SeqCap EZ v3.0 platform, including the mutation missed using Illumina TruSeq capture. Target capture for exons in known PCC/PGL genes differs substantially between platforms. Exome sequencing was inexpensive (<$A800 per sample for reagents) and rapid (results <5 weeks from sample reception). CONCLUSION: Whole exome sequencing is sensitive, rapid and efficient for detection of PCC/PGL germline mutations. However, capture platform selection is critical to maximize sensitivity.

A global role for KLF1 in erythropoiesis revealed by ChIP-seq in primary erythroid cells.

KLF1 regulates a diverse suite of genes to direct erythroid cell differentiation from bipotent progenitors. To determine the local cis-regulatory contexts and transcription factor networks in which KLF1 operates, we performed KLF1 ChIP-seq in the mouse. We found at least 945 sites in the genome of E14.5 fetal liver erythroid cells which are occupied by endogenous KLF1. Many of these recovered sites reside in erythroid gene promoters such as Hbb-b1, but the majority are distant to any known gene. Our data suggests KLF1 directly regulates most aspects of terminal erythroid differentiation including production of alpha- and beta-globin protein chains, heme biosynthesis, coordination of proliferation and anti-apoptotic pathways, and construction of the red cell membrane and cytoskeleton by functioning primarily as a transcriptional activator. Additionally, we suggest new mechanisms for KLF1 cooperation with other transcription factors, in particular the erythroid transcription factor GATA1, to maintain homeostasis in the erythroid compartment.

Regulated post-transcriptional RNA cleavage diversifies the eukaryotic transcriptome.

The complexity of the eukaryotic transcriptome is generated by the interplay of transcription initiation, termination, alternative splicing, and other forms of post-transcriptional modification. It was recently shown that RNA transcripts may also undergo cleavage and secondary 5' capping. Here, we show that post-transcriptional cleavage of RNA contributes to the diversification of the transcriptome by generating a range of small RNAs and long coding and noncoding RNAs. Using genome-wide histone modification and RNA polymerase II occupancy data, we confirm that the vast majority of intraexonic CAGE tags are derived from post-transcriptional processing. By comparing exonic CAGE tags to tissue-matched PARE data, we show that the cleavage and subsequent secondary capping is regulated in a developmental-stage- and tissue-specific manner. Furthermore, we find evidence of prevalent RNA cleavage in numerous transcriptomic data sets, including SAGE, cDNA, small RNA libraries, and deep-sequenced size-fractionated pools of RNA. These cleavage products include mRNA variants that retain the potential to be translated into shortened functional protein isoforms. We conclude that post-transcriptional RNA cleavage is a key mechanism that expands the functional repertoire and scope for regulatory control of the eukaryotic transcriptome.

Refining transcriptional programs in kidney development by integration of deep RNA-sequencing and array-based spatial profiling.

BACKGROUND: The developing mouse kidney is currently the best-characterized model of organogenesis at a transcriptional level. Detailed spatial maps have been generated for gene expression profiling combined with systematic in situ screening. These studies, however, fall short of capturing the transcriptional complexity arising from each locus due to the limited scope of microarray-based technology, which is largely based on "gene-centric" models. RESULTS: To address this, the polyadenylated RNA and microRNA transcriptomes of the 15.5 dpc mouse kidney were profiled using strand-specific RNA-sequencing (RNA-Seq) to a depth sufficient to complement spatial maps from pre-existing microarray datasets. The transcriptional complexity of RNAs arising from mouse RefSeq loci was catalogued; including 3568 alternatively spliced transcripts and 532 uncharacterized alternate 3' UTRs. Antisense expressions for 60% of RefSeq genes was also detected including uncharacterized non-coding transcripts overlapping kidney progenitor markers, Six2 and Sall1, and were validated by section in situ hybridization. Analysis of genes known to be involved in kidney development, particularly during mesenchymal-to-epithelial transition, showed an enrichment of non-coding antisense transcripts extended along protein-coding RNAs. CONCLUSION: The resulting resource further refines the transcriptomic cartography of kidney organogenesis by integrating deep RNA sequencing data with locus-based information from previously published expression atlases. The added resolution of RNA-Seq has provided the basis for a transition from classical gene-centric models of kidney development towards more accurate and detailed "transcript-centric" representations, which highlights the extent of transcriptional complexity of genes that direct complex development events.

Stem cell transcriptome profiling via massive-scale mRNA sequencing.

We developed a massive-scale RNA sequencing protocol, short quantitative random RNA libraries or SQRL, to survey the complexity, dynamics and sequence content of transcriptomes in a near-complete fashion. This method generates directional, random-primed, linear cDNA libraries that are optimized for next-generation short-tag sequencing. We surveyed the poly(A)(+) transcriptomes of undifferentiated mouse embryonic stem cells (ESCs) and embryoid bodies (EBs) at an unprecedented depth (10 Gb), using the Applied Biosystems SOLiD technology. These libraries capture the genomic landscape of expression, state-specific expression, single-nucleotide polymorphisms (SNPs), the transcriptional activity of repeat elements, and both known and new alternative splicing events. We investigated the impact of transcriptional complexity on current models of key signaling pathways controlling ESC pluripotency and differentiation, highlighting how SQRL can be used to characterize transcriptome content and dynamics in a quantitative and reproducible manner, and suggesting that our understanding of transcriptional complexity is far from complete.

Kidney transcriptome reveals altered steroid homeostasis in NaS1 sulfate transporter null mice.

Sulfate is essential for human growth and development, and circulating sulfate levels are maintained by the NaS1 sulfate transporter which is expressed in the kidney. Previously, we generated a NaS1-null (Nas1(-/-)) mouse which exhibits hyposulfatemia. In this study, we investigated the kidney transcriptome of Nas1(-/-) mice. We found increased (n=25) and decreased (n=60) mRNA levels of genes with functional roles that include sulfate transport and steroid metabolism. Corticosteroid-binding globulin was the most up-regulated gene (110% increase) in Nas1(-/-) mouse kidney, whereas the sulfate anion transporter-1 (Sat1) was among the most down-regulated genes (>or=50% decrease). These findings led us to investigate the circulating and urinary steroid levels of Nas1(-/-) and Nas1(+/+) mice, which revealed reduced blood levels of corticosterone ( approximately 50% decrease), dehydroepiandrosterone (DHEA, approximately 30% decrease) and DHEA-sulfate ( approximately 40% decrease), and increased urinary corticosterone ( approximately 16-fold increase) and DHEA ( approximately 40% increase) levels in Nas1(-/-) mice. Our data suggest that NaS1 is essential for maintaining a normal metabolic state in the kidney and that loss of NaS1 function leads to reduced circulating steroid levels and increased urinary steroid excretion.

Transcriptional analysis of early lineage commitment in human embryonic stem cells.

BACKGROUND: The mechanisms responsible for the maintenance of pluripotency in human embryonic stem cells, and those that drive their commitment into particular differentiation lineages, are poorly understood. In fact, even our knowledge of the phenotype of hESC is limited, because the immunological and molecular criteria presently used to define this phenotype describe the properties of a heterogeneous population of cells. RESULTS: We used a novel approach combining immunological and transcriptional analysis (immunotranscriptional profiling) to compare gene expression in hESC populations at very early stages of differentiation. Immunotranscriptional profiling enabled us to identify novel markers of stem cells and their differentiated progeny, as well as novel potential regulators of hESC commitment and differentiation. The data show clearly that genes associated with the pluripotent state are downregulated in a coordinated fashion, and that they are co-expressed with lineage specific transcription factors in a continuum during the early stages of stem cell differentiation. CONCLUSION: These findings, that show that maintenance of pluripotency and lineage commitment are dynamic, interactive processes in hESC cultures, have important practical implications for propagation and directed differentiation of these cells, and for the interpretation of mechanistic studies of hESC renewal and commitment. Since embryonic stem cells at defined stages of commitment can be isolated in large numbers by immunological means, they provide a powerful model for studying molecular genetics of stem cell commitment in the embryo.

Dynamic transcription programs during ES cell differentiation towards mesoderm in serum versus serum-freeBMP4 culture.

BACKGROUND: Expression profiling of embryonic stem (ES) cell differentiation in the presence of serum has been performed previously. It remains unclear if transcriptional activation is dependent on complex growth factor mixtures in serum or whether this process is intrinsic to ES cells once the stem cell program has been inactivated. The aims of this study were to determine the transcriptional programs associated with the stem cell state and to characterize mesoderm differentiation between serum and serum-free culture. RESULTS: ES cells were differentiated as embryoid bodies in 10% FBS or serum-free media containing BMP4 (2 ng/ml), and expression profiled using 47 K Illumina(R) Sentrix arrays. Statistical methods were employed to define gene sets characteristic of stem cell, epiblast and primitive streak programs. Although the initial differentiation profile was similar between the two culture conditions, cardiac gene expression was inhibited in serum whereas blood gene expression was enhanced. Also, expression of many members of the Kruppel-like factor (KLF) family of transcription factors changed dramatically during the first few days of differentiation. KLF2 and KLF4 co-localized with OCT4 in a sub-nuclear compartment of ES cells, dynamic changes in KLF-DNA binding activities occurred upon differentiation, and strong bio-informatic evidence for direct regulation of many stem cell genes by KLFs was found. CONCLUSION: Down regulation of stem cell genes and activation of epiblast/primitive streak genes is similar in serum and defined media, but subsequent mesoderm differentiation is strongly influenced by the composition of the media. In addition, KLF family members are likely to be important regulators of many stem cell genes.

Osteonectin downregulates E-cadherin, induces osteopontin and focal adhesion kinase activity stimulating an invasive melanoma phenotype.

Osteonectin is recognised as a marker of metastasis progression in melanoma and has been implicated in the transition from radial to vertical growth phase. A Tetracycline-inducible system was used to regulate Osteonectin protein levels in melanoma cell lines to examine the morphological, biochemical and invasive changes that accompany its altered expression. Assay of protein and phosphorylation changes showed a downregulation of E-cadherin, upregulation of Osteopontin and a corresponding increase in phosphorylation of Focal Adhesion Kinase on Tyr(397) and Tyr(576) concomitant with Osteonectin induction. Melanoma cells overexpressing Osteonectin displayed increased invasive potential, whereas ablation of Osteonectin gene transcription using siRNA suppressed the invasive potential of these cells and resulted in the upregulation of E-cadherin. The recently described interaction of Osteonectin with Integrin Linked Kinase leading to modulation of its activity suggests a mechanism relevant to the loss of E-cadherin and cell adhesion that occurs during melanoma progression. These results indicate a central role for Osteonectin in the regulation of gene expression changes driving the progression of melanoma toward metastasis.

Transcriptional profile reveals altered hepatic lipid and cholesterol metabolism in hyposulfatemic NaS1 null mice.

Sulfate plays an essential role in human growth and development, and its circulating levels are maintained by the renal Na+-SO42- cotransporter, NaS1. We previously generated a NaS1 knockout (Nas1-/-) mouse, an animal model for hyposulfatemia, that exhibits reduced growth and liver abnormalities including hepatomegaly. In this study, we investigated the hepatic gene expression profile of Nas1-/- mice using oligonucleotide microarrays. The mRNA expression levels of 92 genes with known functional roles in metabolism, cell signaling, cell defense, immune response, cell structure, transcription, or protein synthesis were increased (n = 51) or decreased (n = 41) in Nas1-/- mice when compared with Nas1+/+ mice. The most upregulated transcript levels in Nas1-/- mice were found for the sulfotransferase genes, Sult3a1 (approximately 500% increase) and Sult2a2 (100% increase), whereas the metallothionein-1 gene, Mt1, was among the most downregulated genes (70% decrease). Several genes involved in lipid and cholesterol metabolism, including Scd1, Acly, Gpam, Elov16, Acsl5, Mvd, Insig1, and Apoa4, were found to be upregulated (> or = 30% increase) in Nas1-/- mice. In addition, Nas1-/- mice exhibited increased levels of hepatic lipid (approximately 16% increase), serum cholesterol (approximately 20% increase), and low-density lipoprotein (approximately 100% increase) and reduced hepatic glycogen (approximately 50% decrease) levels. In conclusion, these data suggest an altered lipid and cholesterol metabolism in the hyposulfatemic Nas1-/- mouse and provide new insights into the metabolic state of the liver in Nas1-/- mice.

Spatial gene expression in the T-stage mouse metanephros.

The E11.5 mouse metanephros is comprised of a T-stage ureteric epithelial tubule sub-divided into tip and trunk cells surrounded by metanephric mesenchyme (MM). Tip cells are induced to undergo branching morphogenesis by the MM. In contrast, signals within the mesenchyme surrounding the trunk prevent ectopic branching of this region. In order to identify novel genes involved in the molecular regulation of branching morphogenesis we compared the gene expression profiles of isolated tip, trunk and MM cells using Compugen mouse long oligo microarrays. We identified genes enriched in the tip epithelium, sim-1, Arg2, Tacstd1, Crlf-1 and BMP7; genes enriched in the trunk epithelium, Innp1, Itm2b, Mkrn1, SPARC, Emu2 and Gsta3 and genes spatially restricted to the mesenchyme surrounding the trunk, CSPG2 and CV-2, with overlapping and complimentary expression to BMP4, respectively. This study has identified genes spatially expressed in regions of the developing kidney involved in branching morphogenesis, nephrogenesis and the development of the collecting duct system, calyces, renal pelvis and ureter.

Osteonectin/SPARC induction by ectopic beta(3) integrin in human radial growth phase primary melanoma cells.

Expression of the beta(3) integrin subunit in melanoma in situ has been found to correlate with tumor thickness, the ability to invade and metastasize, and poor prognosis. Transition from the radial growth phase (RGP) to the vertical growth phase (VGP) is a critical step in melanoma progression and survival and is distinguished by the expression of beta(3) integrin. The molecular pathways that operate in melanoma cells associated with invasion and metastasis were examined by ectopic induction of the beta(3) integrin subunit in RGP SBcl2 and WM1552C melanoma cells, which converts these cells to a VGP phenotype. We used cDNA representational difference analysis subtractive hybridization between beta(3)-positive and -negative melanoma cells to assess gene expression profile changes accompanying RGP to VGP transition. Fourteen fragments from known genes including osteonectin (also known as SPARC and BM-40) were identified after three rounds of representational difference analysis. Induction of osteonectin was confirmed by Northern and Western blot analysis and immunohistochemistry and correlated in organotypic cultures with the beta(3)-induced progression from RGP to VGP melanoma. Expression of osteonectin was also associated with reduced adhesion to vitronectin, but not to fibronectin. Osteonectin expression was not blocked when melanoma cells were cultured with anti-alpha(v)beta(3) LM609 mAb, mitogen-activated protein kinase, or protein kinase C inhibitors, indicating that other signaling pathway(s) operate through alpha(v)beta(3) integrin during conversion from RGP to VGP.