The primacy of NF1 loss as the driver of tumorigenesis in neurofibromatosis type 1-associated plexiform neurofibromas.

Neurofibromatosis type 1 (NF1) is a common tumor-predisposition disorder due to germline mutations in the tumor suppressor gene NF1. A virtually pathognomonic finding of NF1 is the plexiform neurofibroma (PN), a benign, likely congenital tumor that arises from bi-allelic inactivation of NF1. PN can undergo transformation to a malignant peripheral nerve sheath tumor, an aggressive soft-tissue sarcoma. To better understand the non-NF1 genetic contributions to PN pathogenesis, we performed whole-exome sequencing, RNASeq profiling and genome-wide copy-number determination for 23 low-passage Schwann cell cultures established from surgical PN material with matching germline DNA. All resected tumors were derived from routine debulking surgeries. None of the tumors were considered at risk for malignant transformation at the time; for example, there was no pain or rapid growth. Deep (~500X) NF1 exon sequencing was also conducted on tumor DNA. Non-NF1 somatic mutation verification was performed using the Ampliseq/IonTorrent platform. We identified 100% of the germline NF1 mutations and found somatic NF1 inactivation in 74% of the PN. One individual with three PNs had different NF1 somatic mutations in each tumor. The median number of somatic mutations per sample, including NF1, was one (range 0-8). NF1 was the only gene that was recurrently somatically inactivated in multiple tumors. Gene Set Enrichment Analysis of transcriptome-wide tumor RNA sequencing identified five significant (FDR<0.01) and seven trending (0.01⩽FDR<0.02) gene sets related to DNA replication, telomere maintenance and elongation, cell cycle progression, signal transduction and cell proliferation. We found no recurrent non-NF1 locus copy-number variation in PN. This is the first multi-sample whole-exome and whole-transcriptome sequencing study of NF1-associated PN. Taken together with concurrent copy-number data, our comprehensive genetic analysis reveals the primacy of NF1 loss as the driver of PN tumorigenesis.

Novel integrative genomic tool for interrogating lithium response in bipolar disorder.

We developed a novel integrative genomic tool called GRANITE (Genetic Regulatory Analysis of Networks Investigational Tool Environment) that can effectively analyze large complex data sets to generate interactive networks. GRANITE is an open-source tool and invaluable resource for a variety of genomic fields. Although our analysis is confined to static expression data, GRANITE has the capability of evaluating time-course data and generating interactive networks that may shed light on acute versus chronic treatment, as well as evaluating dose response and providing insight into mechanisms that underlie therapeutic versus sub-therapeutic doses or toxic doses. As a proof-of-concept study, we investigated lithium (Li) response in bipolar disorder (BD). BD is a severe mood disorder marked by cycles of mania and depression. Li is one of the most commonly prescribed and decidedly effective treatments for many patients (responders), although its mode of action is not yet fully understood, nor is it effective in every patient (non-responders). In an in vitro study, we compared vehicle versus chronic Li treatment in patient-derived lymphoblastoid cells (LCLs) (derived from either responders or non-responders) using both microRNA (miRNA) and messenger RNA gene expression profiling. We present both Li responder and non-responder network visualizations created by our GRANITE analysis in BD. We identified by network visualization that the Let-7 family is consistently downregulated by Li in both groups where this miRNA family has been implicated in neurodegeneration, cell survival and synaptic development. We discuss the potential of this analysis for investigating treatment response and even providing clinicians with a tool for predicting treatment response in their patients, as well as for providing the industry with a tool for identifying network nodes as targets for novel drug discovery.

Pharmacogenomics, ancestry and clinical decision making for global populations.

Pharmacogenomically relevant markers of drug response and adverse drug reactions are known to vary in frequency across populations. We examined minor allele frequencies (MAFs), genetic diversity (FST) and population structure of 1156 genetic variants (including 42 clinically actionable variants) in 212 genes involved in drug absorption, distribution, metabolism and excretion (ADME) in 19 populations (n=1478). There was wide population differentiation in these ADME variants, reflected in the range of mean MAF (ΔMAF) and FST. The largest mean ΔMAF was observed in African ancestry populations (0.10) and the smallest mean ΔMAF in East Asian ancestry populations (0.04). MAFs ranged widely, for example, from 0.93 for single-nucleotide polymorphism (SNP) rs9923231, which influences warfarin dosing to 0.01 for SNP rs3918290 associated with capecitabine metabolism. ADME genetic variants show marked variation between and within continental groupings of populations. Enlarging the scope of pharmacogenomics research to include multiple global populations can improve the evidence base for clinical translation to benefit all peoples.

Genetic screening for von Hippel-Lindau gene mutations in non-syndromic pheochromocytoma: low prevalence and false-positives or misdiagnosis indicate a need for caution.

Genetic testing of tumor susceptibility genes is now recommended in most patients with pheochromocytoma or paraganglioma (PPGL), even in the absence of a syndromic presentation. Once a mutation is diagnosed there is rarely follow-up validation to assess the possibility of misdiagnosis. This study prospectively examined the prevalence of von Hippel-Lindau (VHL) gene mutations among 182 patients with non-syndromic PPGLs. Follow-up in positive cases included comparisons of biochemical and tumor gene expression data in 64 established VHL patients, with confirmatory genetic testing in cases with an atypical presentation. VHL mutations were detected by certified laboratory testing in 3 of the 182 patients with non-syndromic PPGLs. Two of the 3 had an unusual presentation of diffuse peritoneal metastases and substantial increases in plasma metanephrine, the metabolite of epinephrine. Tumor gene expression profiles in these 2 patients also differed markedly from those associated with established VHL syndrome. One patient was diagnosed with a partial deletion by Southern blot analysis and the other with a splice site mutation. Quantitative polymerase chain reaction, multiplex ligation-dependent probe amplification, and comparative genomic hybridization failed to confirm the partial deletion indicated by certified laboratory testing. Analysis of tumor DNA in the other patient with a splice site alteration indicated no loss of heterozygosity or second hit point mutation. In conclusion, VHL germline mutations represent a minor cause of non-syndromic PPGLs and misdiagnoses can occur. Caution should therefore be exercised in interpreting positive genetic test results as the cause of disease in patients with non-syndromic PPGLs.

Chromogranin a expression in phaeochromocytomas associated with von Hippel-Lindau syndrome and multiple endocrine neoplasia type 2.

Chromogranin A (CGA) is a major secretory protein present in the soluble matrix of chromaffin granules of neuroendocrine cells and tumours, such as phaeochromocytomas. CGA has several functions, some of which may be involved in the distinct phenotypic differences of phaeochromocytomas in patients with von Hippel-Lindau (VHL) syndrome compared to multiple endocrine neoplasia type 2 (MEN 2). In this study, we therefore compared tumour and plasma levels of CGA in patients with phaeochromocytoma associated with the two syndromes. We show that phaeochromocytomas from MEN 2 patients express substantially more CGA than tumours from VHL patients at both the mRNA (3-fold greater) and protein (20-fold) level. We further show that relative to increases in plasma catecholamines, patients with phaeochromocytomas associated with MEN 2 have higher plasma concentrations of CGA than those with tumours in VHL syndrome. These data supplement other observations that phaeochromocytomas in VHL compared to MEN 2 patients express lower amounts of catecholamines and other chromaffin granule cargo, such as chromogranin B and neuropeptide Y. Possibly the differences in tumour CGA expression may contribute to differences in secretory vesicle formation and secretion in the two types of tumours. Alternatively the differences in expression in CGA and other secretory constituents may reflect downregulation of the entire regulated secretory pathway in VHL compared to MEN 2 tumours.

Distinct gene expression profiles in norepinephrine- and epinephrine-producing hereditary and sporadic pheochromocytomas: activation of hypoxia-driven angiogenic pathways in von Hippel-Lindau syndrome.

Pheochromocytomas in von Hippel-Lindau (VHL) syndrome produce exclusively norepinephrine, whereas those in multiple endocrine neoplasia type 2 (MEN 2) produce epinephrine. This study examined the pathways activated in VHL-associated pheochromocytomas by comparing gene expression profiles in VHL and MEN 2 tumors in relationship to profiles in sporadic norepinephrine- and epinephrine-producing tumors. Larger and more distinct differences in gene expression among hereditary than sporadic tumors indicated the importance of the underlying mutation to gene expression profiles. Many of the genes over-expressed in VHL compared with MEN 2 tumors were clearly linked to the hypoxia-driven angiogenic pathways that are activated in VHL-associated tumorigenesis. Such genes included those for the glucose transporter, vascular endothelial growth factor (VEGF), placental growth factor, angiopoietin 2, tie-1, VEGF receptor 2 and its coreceptor, neuropilin-1. Other up-regulated genes, such as connective tissue growth factor, cysteine-rich 61, matrix metalloproteinase 1, vascular endothelial cadherin, tenascin C, stanniocalcin 1, and cyclooxygenases 1 and 2 are known to be involved in VEGF-regulated angiogenesis. Shared differences in expression of subsets of genes in norepinephrine- versus epinephrine-producing hereditary and sporadic pheochromocytomas indicated other differences in gene expression that may underlie the biochemical phenotype. Over-expression of the hypoxia-inducible transcription factor, HIF-2alpha, in norepinephrine-predominant sporadic and VHL tumors compared with epinephrine-producing tumors indicates that expression of this gene depends on the noradrenergic biochemical phenotype. The findings fit with the known expression of HIF-2alpha in norepinephrine-producing cells of the sympathetic nervous system and might explain both the development and noradrenergic biochemical phenotype of pheochromocytomas in VHL syndrome.

Angiotensin II AT1 receptor antagonism downregulates stress-related gene expression in brain microvessels from spontaneously hypertensive and normotensive rats.

We studied the effect of treatment with the Angiotensin II AT(1) receptor antagonist candesartan (0.3 mg/kg/day via osmotic minipumps for 4 weeks compared with administration of vehicle) in brain microvessels in adult spontaneously hypertensive rats (SHR) that were vulnerable to stroke and normotensive control rats (WKY). At the dose administered, candesartan normalized blood pressure in SHR without significantly affecting blood pressure in WKY rats. We performed the gene expression analysis in rat brain microvessels using the Affymetrix Gene Chip Expression Analysis Technique. From a total of 8,799 probe array sets analyzed, we found abundant abnormalities in gene expression in SHR. Because stress has been suggested as a precipitant factor in brain ischemia and treatment with AT(1) receptor antagonist candesartan prevents the hormonal and sympathoadrenal reaction to isolation stress and protects from stress-induced gastric ulcers, we focused on the expression of stress-related genes. We found a higher number of probe array sets modified by candesartan treatment in normotensive WKY rats than in hypertensive SHR. AT(1) receptor blockade decreased the transcription levels of the stress-related tyrosine kinase receptor, stathmin, and fibroblast growth receptor genes in WKY and SHR rats. Our results indicate that Angiotensin II and its AT(1) receptors can influence gene expression independently of the effects on blood pressure. In addition, AT(1) receptor regulation of stress-related genes in brain microvessels may explain the proposed association between stress and ischemic disorders of the brain.

Failure of hormone therapy in prostate cancer involves systematic restoration of androgen responsive genes and activation of rapamycin sensitive signaling.

Androgen deprivation therapy for advanced prostate cancer is often effective, but not curative. Molecular pathways mediating the therapeutic response and those contributing to the subsequent hormone-refractory cell growth remain poorly understood. Here, cDNA microarray analysis of human CWR22 prostate cancer xenografts during the course of androgen deprivation therapy revealed distinct global gene expression profiles in primary, regressing and recurrent tumors. Elucidation of the genes involved in the transition between these states implicated specific molecular mechanisms in therapy failure and tumor progression. First, we identified a set of androgen-responsive genes whose expression decreased during the therapy response, but was then systematically restored in the recurrent tumors. In addition, altered expression of genes that encode known targets of rapamycin or that converge on the PI3K/AKT/FRAP pathway was observed in the recurrent tumors. Further suggestion for the involvement of these genes in hormone-refractory prostate cancer came from the observation that cells established from the recurrent xenografts were strongly inhibited in vitro by rapamycin. The results of this functional genomic analysis suggest that the combined effect of re-expression of androgen-responsive genes as well as the activation of rapamycin-sensitive signaling may drive prostate cancer progression, and contribute to the failure of androgen-deprivation therapy.

Gastrointestinal stromal tumors with KIT mutations exhibit a remarkably homogeneous gene expression profile.

Gastrointestinal stromal tumors (GISTs), the most common mesenchymal tumors of the digestive tract, are believed to arise from the interstitial cells of Cajal. GISTs are characterized by mutations in the proto-oncogene KIT that lead to constitutive activation of its tyrosine kinase activity. The tyrosine kinase inhibitor STI571, active against the BCR-ABL fusion protein in chronic myeloid leukemia, was recently shown to be highly effective in GISTs. We used 13,826-element cDNA microarrays to define the expression patterns of 13 KIT mutation-positive GISTs and compared them with the expression profiles of a group of spindle cell tumors from locations outside the gastrointestinal tract. Our results showed a remarkably distinct and uniform expression profile for all of the GISTs. In particular, hierarchical clustering of a subset of 113 cDNAs placed all of the GIST samples into one branch, with a Pearson correlation >0.91. This homogeneity suggests that the molecular pathogenesis of a GIST results from expansion of a clone that has acquired an activating mutation in KIT without the extreme genetic instability found in the common epithelial cancers. The results provide insight into the histogenesis of GIST and the clinical behavior of this therapeutically responsive tumor.

A cyclooxygenase-2 inhibitor (SC-58125) blocks growth of established human colon cancer xenografts.

Selective COX-2 inhibitors reduce adenoma formation and cancer progression in rodent models of colorectal cancer. To assess the therapeutic activity of selective COX-2 inhibitors, we tested the effect of SC-58125 treatment on the growth of human colon carcinoma cells in nude mice. Delaying treatment by 2, 4, or 7 weeks following implantation of the carcinoma cells resulted in a significant inhibition of tumor growth. Furthermore, short-term (48 hours) treatment with SC-58125 was sufficient to attenuate tumor growth for up to 15 days. SC-58125 treatment did not alter the rate at which cells underwent apoptosis, but did result in a delayed progression through the cell cycle at the G(2)/M transition. Accordingly, p34(cdc2) protein levels and activity were decreased following SC-58125 treatment. We conclude that SC-58125 primarily exerts a cytostatic effect in vivo, which is likely to be mediated through inhibition of progression through the G(2)/M phase of the cell cycle.

Two-dimensional DNA electrophoresis identifies novel CpG islands frequently coamplified with MYCN in neuroblastoma.

BACKGROUND: Amplification of the oncogene MYCN in neuroblastoma has been found to correlate with aggressive tumour growth and is used as a predictor of clinical outcome. The MYCN amplicon is known to involve coamplification of extensive DNA regions. Therefore it is possible that other genes are coamplified in this amplicon and that they may play a role in the poor outcome of MYCN amplified tumours. PROCEDURE: We have implemented an approach for the two-dimensional separation of human genomic restriction fragments to detect and isolate as yet unknown amplified sequences in the MYCN amplicon in neuroblastoma. Using this approach we have recently cloned a novel gene referred to as NAG that is frequently coamplified with MYCN in neuroblastoma. RESULTS AND CONCLUSIONS: We report here the identification and cloning of two additional CpG islands that are amplified in neuroblastoma. One contains a sequence that is identical to the first intron of DDX1. The other represents a novel CpG island that is associated with an as yet unidentified gene. We show that the novel CpG island is located in close proximity to the MYCN locus on chromosome 2 and is as frequently coamplified with MYCN in neuroblastoma as NAG and DDX1.

Induction of stearoyl CoA desaturase is associated with high-level induction of emerin RNA.

The genes encoding enzymes involved in fatty acid metabolism are regulated by sterols. Stearoyl CoA desaturase, a key enzyme in the synthesis of unsaturated fatty acyl-CoAs is transcriptionally regulated by fat-free diet and sterols. To identify other genes that are induced in rat liver by fat-free diet we performed an in vivo gene expression profile analysis using DNA microarrays. Here we report that among the genes highly expressed is emerin, an integral protein of the inner nuclear membrane. Mutated or nonexpressed emerin occurs in patients with muscular dystrophy. Sterol regulatory element binding proteins activate the transcription of several sterol regulated genes. To investigate whether sterol regulatory element binding proteins or indirectly cholesterol activates the transcription of stearoyl CoA desaturase and emerin, we cultured Chinese hamster ovary (CHO), either in cholesterol-rich or cholesterol-depleted mediums. We also transiently transfected the cell culture with plasmid encoding sterol regulatory element binding proteins in cholesterol-rich medium. Our data show that cholesterol-supplemented media as well as the transient transfection induced the expression of stearoyl CoA desaturase RNA 3.5- and 7-fold respectively. However, the RNA level of emerin was not altered under these conditions, implying that the parallel induction of emerin is independent of the sterol regulatory element binding regulation pathway.

Elevation of breast carcinoma Nm23-H1 metastasis suppressor gene expression and reduced motility by DNA methylation inhibition.

We hypothesize that elevation of Nm23-H1 expression in micrometastatic breast cancer cells may inhibit their metastatic colonization and further invasion, and induce differentiation, thus resulting in a clinical benefit. The current study investigated the possible contribution of DNA methylation to the regulation of Nm23-H1 expression, based on the observation that two CpG islands are present in its promoter. 5-Aza-2'-deoxycytidine (5-Aza-CdR), a DNA methylation inhibitor, increased the Nm23-H1 expression of 5 of 11 human breast carcinoma cell lines in vitro, including 3 of 3 metastatically competent lines. Increased Nm23-H1 expression was accompanied by a reduction in motility in vitro, with minimal effect on proliferation. Both increased Nm23-H1 expression and decreased motility were observed using low (75 nM) concentrations of 5-Aza-CdR. Array analysis of MDA-MB-231 breast carcinoma cells treated with 5-Aza-CdR confirmed the elevation of nm23-H1 mRNA, whereas relatively few other genes exhibited altered expression. Bisulfite sequencing of the two CpG islands in a panel of cell lines and in 20 infiltrating ductal carcinomas revealed that one island (-3090 bp to -3922 bp) exhibited infrequent differential methylation. The data indicate that DNA methylation inhibitors can directly or indirectly cause both elevation of Nm23-H1 expression and decreased function in one aspect of metastasis, motility.

Monitoring the expression profiles of doxorubicin-induced and doxorubicin-resistant cancer cells by cDNA microarray.

Drug resistance in cancer is a major obstacle to successful chemotherapy. Cancer cells exposed to antitumor drugs may be directly induced to express a subset of genes that could confer resistance, thus allowing some cells to escape killing and form the relapsed resistant tumor. Alternatively, some cancer cells may be expressing an array of genes that could confer intrinsic resistance, and exposure to cytotoxic drugs select for the survival of these cells that form the relapsed tumor. We have used cDNA microarray to monitor the expression profiles of MCF-7 cells that are either transiently treated with doxorubicin or selected for resistance to doxorubicin. Our results showed that transient treatment with doxorubicin altered the expression of a diverse group of genes in a time-dependent manner. A subset of the induced genes was also found to be constitutively overexpressed in cells selected for resistance to doxorubicin. This distinct set of overlapping genes may represent the signature profile of doxorubicin-induced gene expression and resistance in cancer cells. Our studies demonstrate the feasibility of obtaining potential molecular profile or fingerprint of anticancer drugs in cancer cells by cDNA microarray, which might yield further insights into the mechanisms of drug resistance and suggest alternative methods of treatment.

Identification of a larger than 3 Mb deletion including JAG1 in an Alagille syndrome patient with a translocation t(3;20)(q13.3;p12.2).

Alagille syndrome (AGS) is an autosomal dominant, developmental disorder affecting multiple organ systems including liver, heart, vertebrae, eye and face. Recurrent deletions of the 20p12 region led to the localization, and ultimately to the identification of mutations in the Jagged1 gene (JAG1) in AGS patients. A translocation t(3;20)(q13.3;p12.2) in an AGS patient was characterized using fluorescent in situ hybridization (FISH). The involvement of 3q and 20p in this translocation was demonstrated using probes for 3q and 20p. Three overlapping YAC clones, 940D11, 953A2, and 675G11 extending to nearly 4 Mb including the JAG1 were used as probes for FISH analysis to define the translocation breakpoint. The translocated chromosome was found to have a deletion of more than 3 Mb including the entire JAG1 gene. The observation of an accompanying large deletion, revealed by molecular characterization of the t(3;20) translocation, is similar to the only other translocation reported in an AGS patient; a t(2;20) translocation was also found to have a large deletion of the JAG1 region at 20p12.

Development of a prostate cDNA microarray and statistical gene expression analysis package.

A cDNA microarray comprising 5184 different cDNAs spotted onto nylon membrane filters was developed for prostate gene expression studies. The clones used for arraying were identified by cluster analysis of > 35 000 prostate cDNA library-derived expressed sequence tags (ESTs) present in the dbEST database maintained by the National Center for Biotechnology Information. Total RNA from two cell lines, prostate line 8.4 and melanoma line UACC903, was used to make radiolabeled probe for filter hybridizations. The absolute intensity of each individual cDNA spot was determined by phosphorimager scanning and evaluated by a bioinformatics package developed specifically for analysis of cDNA microarray experimentation. Results indicated 89% of the genes showed intensity levels above background in prostate cells compared with only 28% in melanoma cells. Replicate probe preparations yielded results with correlation values ranging from r = 0.90 to 0.93 and coefficient of variation ranging from 16 to 28%. Findings indicate that among others, the keratin 5 and vimentin genes were differentially expressed between these two divergent cell lines. Follow-up northern blot analysis verified these two expression changes, thereby demonstrating the reliability of this system. We report the development of a cDNA microarray system that is sensitive and reliable, demonstrates a low degree of variability, and is capable of determining verifiable gene expression differences between two distinct human cell lines. This system will prove useful for differential gene expression analysis in prostate-derived cells and tissue.

In vivo gene expression profile analysis of human breast cancer progression.

The development and use of molecular-based therapy for breast cancer and other human malignancies will require a detailed molecular genetic analysis of patient tissues. The recent development of laser capture microdissection and high density cDNA arrays now provides a unique opportunity to generate gene expression profiles of cells from various stages of tumor progression as it occurs in the actual neoplastic tissue milieu. We report the combined use of laser capture microdissection and high-throughput cDNA microarrays to monitor in vivo gene expression levels in purified normal, invasive, and metastatic breast cell populations from a single patient. These in vivo gene expression profiles were verified by real-time quantitative PCR and immunohistochemistry. The combined use of laser capture microdissection and cDNA microarray analysis provides a powerful new approach to elucidate the in vivo molecular events surrounding the development and progression of breast cancer and is generally applicable to the study of malignancy.

Hormone therapy failure in human prostate cancer: analysis by complementary DNA and tissue microarrays.

BACKGROUND: The molecular mechanisms underlying the progression of prostate cancer during hormonal therapy have remained poorly understood. In this study, we developed a new strategy for the identification of differentially expressed genes in hormone-refractory human prostate cancer by use of a combination of complementary DNA (cDNA) and tissue microarray technologies. METHODS: Differences in gene expression between hormone-refractory CWR22R prostate cancer xenografts (human prostate cancer transplanted into nude mice) and a xenograft of the parental, hormone-sensitive CWR22 strain were analyzed by use of cDNA microarray technology. To validate the data from cDNA microarrays on clinical prostate cancer specimens, a tissue microarray of specimens from 26 prostates with benign prostatic hyperplasia, 208 primary prostate cancers, and 30 hormone-refractory local recurrences was constructed and used for immunohistochemical detection of protein expression. RESULTS: Among 5184 genes surveyed with cDNA microarray technology, expression of 37 (0.7%) was increased more than twofold in the hormone-refractory CWR22R xenografts compared with the CWR22 xenograft; expression of 135 (2.6%) genes was reduced by more than 50%. The genes encoding insulin-like growth factor-binding protein 2 (IGFBP2) and 27-kd heat-shock protein (HSP27) were among the most consistently overexpressed genes in the CWR22R tumors. Immunohistochemical analysis of tissue microarrays demonstrated high expression of IGFBP2 protein in 100% of the hormone-refractory clinical tumors, in 36% of the primary tumors, and in 0% of the benign prostatic specimens (two-sided P =.0001). Overexpression of HSP27 protein was demonstrated in 31% of the hormone-refractory tumors, in 5% of the primary tumors, and in 0% of the benign prostatic specimens (two-sided P =.0001). CONCLUSIONS: The combination of cDNA and tissue microarray technologies enables rapid identification of genes associated with progression of prostate cancer to the hormone-refractory state and may facilitate analysis of the role of the encoded gene products in the pathogenesis of human prostate cancer.

Identification and cloning of the human homolog (JAG1) of the rat Jagged1 gene from the Alagille syndrome critical region at 20p12.

Notch proteins are a family of closely related transmembrane receptors proven to be instrumental in cell fate decisions. Recently, Notch ligands Delta and Jagged have been identified in Drosophila and rat, respectively. We have isolated the human homolog of the rat Jagged1 gene, JAG1, from a CpG island in a YAC clone covering the Alagille syndrome critical region at chromosome 20p12 (tel-SNAP-D20S186-cen). Alagille syndrome is an autosomal dominant disorder characterized by neonatal jaundice, paucity of intrahepatic bile ducts, and abnormalities of the heart, skeleton, and eyes. The human Jagged1 (JAG1), therefore, appears to be a strong candidate gene for this disease. Here we describe the identification, full-length cDNA cloning, expression patterns, and precise physical location of this gene within the Alagille syndrome critical region.

Mutations in the human Jagged1 gene are responsible for Alagille syndrome.

Alagille syndrome (AGS) is an autosomal-dominant disorder characterized by intrahepatic cholestasis and abnormalities of heart, eye and vertebrae, as well as a characteristic facial appearance. Identification of rare AGS patients with cytogenetic deletions has allowed mapping of the gene of 20p12. We have generated a cloned contig of the critical region and used fluorescent in situ hybridization on cells from patients with submicroscopic deletions to narrow the candidate region to only 250 kb. Within this region we identified JAG1, the human homologue of rat Jagged1, which encodes a ligand for the Notch receptor. Cell-cell Jagged/Notch interactions are known to be critical for determination of cell fates in early development, making this an attractive candidate gene for a developmental disorder in humans. Determining the complete exon-intron structure of JAG1 allowed detailed mutational analysis of DNA samples from non-deletion AGS patients, revealing three frame-shift mutations, two splice donor mutations and one mutation abolishing RNA expression from the altered allele. We conclude that AGS is caused by haploinsufficiency of JAG1.