ORE identifies extreme expression effects enriched for rare variants.

MOTIVATION: Non-coding rare variants (RVs) may contribute to Mendelian disorders but have been challenging to study due to small sample sizes, genetic heterogeneity and uncertainty about relevant non-coding features. Previous studies identified RVs associated with expression outliers, but varying outlier definitions were employed and no comprehensive open-source software was developed. RESULTS: We developed Outlier-RV Enrichment (ORE) to identify biologically-meaningful non-coding RVs. We implemented ORE combining whole-genome sequencing and cardiac RNAseq from congenital heart defect patients from the Pediatric Cardiac Genomics Consortium and deceased adults from Genotype-Tissue Expression. Use of rank-based outliers maximized sensitivity while a most extreme outlier approach maximized specificity. Rarer variants had stronger associations, suggesting they are under negative selective pressure and providing a basis for investigating their contribution to Mendelian disorders. AVAILABILITY AND IMPLEMENTATION: ORE, source code, and documentation are available at https://pypi.python.org/pypi/ore under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Noonan syndrome: clinical aspects and molecular pathogenesis.

Noonan syndrome (NS) is a relatively common, clinically variable and genetically heterogeneous developmental disorder characterized by postnatally reduced growth, distinctive facial dysmorphism, cardiac defects and variable cognitive deficits. Other associated features include ectodermal and skeletal defects, cryptorchidism, lymphatic dysplasias, bleeding tendency, and, rarely, predisposition to hematologic malignancies during childhood. NS is caused by mutations in the PTPN11, SOS1, KRAS, RAF1, BRAF and MEK1 (MAP2K1) genes, accounting for approximately 70% of affected individuals. SHP2 (encoded by PTPN11), SOS1, BRAF, RAF1 and MEK1 positively contribute to RAS-MAPK signaling, and possess complex autoinhibitory mechanisms that are impaired by mutations. Similarly, reduced GTPase activity or increased guanine nucleotide release underlie the aberrant signal flow through the MAPK cascade promoted by most KRAS mutations. More recently, a single missense mutation in SHOC2, which encodes a cytoplasmic scaffold positively controlling RAF1 activation, has been discovered to cause a closely related phenotype previously termed Noonan-like syndrome with loose anagen hair. This mutation promotes aberrantly acquired N-myristoylation of the protein, resulting in its constitutive targeting to the plasma membrane and dysregulated function. PTPN11, BRAF and RAF1 mutations also account for approximately 95% of LEOPARD syndrome, a condition which resembles NS phenotypically but is characterized by multiple lentigines dispersed throughout the body, café-au-lait spots, and a higher prevalence of electrocardiographic conduction abnormalities, obstructive cardiomyopathy and sensorineural hearing deficits. These recent discoveries demonstrate that the substantial phenotypic variation characterizing NS and related conditions can be ascribed, in part, to the gene mutated and even the specific molecular lesion involved.

PTPN11 analysis for the prenatal diagnosis of Noonan syndrome in fetuses with abnormal ultrasound findings.

Noonan syndrome (NS) is an autosomal dominant disorder characterized by short stature, congenital heart defects and distinctive facies. The disorder is genetically heterogeneous with approximately 50% of patients having PTPN11 mutations. Prenatally, the diagnosis of NS has been suspected following certain ultrasound findings, such as cystic hygroma, increased nuchal translucency (NT) and hydrops fetalis. Studies of fetuses with cystic hygroma have suggested an NS prevalence of 1-3%. A retrospective review was performed to assess the utility of PTPN11 testing based on prenatal sonographic findings (n = 134). The most commonly reported indications for testing were increased NT and cystic hygroma. Analysis showed heterozygous missense mutations in 12 fetuses, corresponding to a positive test rate of 9%. PTPN11 mutations were identified in 16% and 2% of fetuses with cystic hygroma and increased NT, respectively. Among fetuses with isolated cystic hygroma, PTPN11 mutation prevalence was 11%. The mutations observed in the three fetuses with hydrops fetalis had previously been reported as somatic cancer mutations. Prenatal PTPN11 testing has diagnostic and possible prognostic properties that can aid in risk assessment and genetic counseling. As NS is genetically heterogeneous, negative PTPN11 testing cannot exclude the diagnosis and further study is warranted regarding the other NS genes.

Thiamine pyrophosphate: an essential cofactor for the alpha-oxidation in mammals--implications for thiamine deficiencies?

The identification of 2-hydroxyphytanoyl-CoA lyase (2-HPCL), a thiamine pyrophosphate (TPP)-dependent peroxisomal enzyme involved in the alpha-oxidation of phytanic acid and of 2-hydroxy straight chain fatty acids, pointed towards a role of TPP in these processes. Until then, TPP had not been implicated in mammalian peroxisomal metabolism. The effect of thiamine deficiency on 2-HPCL and alpha-oxidation has not been studied, nor have possible adverse effects of deficient alpha-oxidation been considered in the pathogenesis of diseases associated with thiamine shortage, such as thiamine-responsive megaloblastic anemia (TRMA). Experiments with cultured cells and animal models showed that alpha-oxidation is controlled by the thiamine status of the cell/tissue/organism, and suggested that some pathological consequences of thiamine starvation could be related to impaired alpha-oxidation. Whereas accumulation of phytanic acid and/or 2-hydroxyfatty acids or their alpha-oxidation intermediates in TRMA patients given a normal supply of thiamine is unlikely, this may not be true when malnourished.

PTPN11 mutations play a minor role in isolated congenital heart disease.

PTPN11 missense mutations cause approximately 50% of Noonan syndrome, an autosomal dominant disorder presenting with various congenital heart defects, most commonly valvar pulmonary stenosis, and hypertrophic cardiomyopathy. Atrioventricular septal defects and coarctation of the aorta occur in 15% and 9%, respectively. The aim of this study was to determine if PTPN11 mutations exist in non-syndromic patients with these two relevant forms of congenital heart disease. The 15 coding PTPN11 exons and their intron boundaries from subjects with atrioventricular septal defects (n = 24) and coarctation of the aorta (n = 157) were analyzed using denaturing high performance liquid chromatography and sequenced if abnormal. One subject with an atrioventricular septal defect but no other known medical problems had a c.127C > T transition in exon 2, predicting a p.L43F substitution. This mutation affected the phosphotyrosine-binding region in the N-terminal src homology 2 domain and was close to a Noonan syndrome mutation (p.T42A). An otherwise healthy patient with aortic coarctation had a silent c.540C > T change in exon 5 corresponding to p.D180D. Our study showed that PTPN11 mutations are rarely found in two isolated forms of congenital heart disease that commonly occur in Noonan syndrome. The p.L43F mutation belongs to a rare class of PTPN11 mutations altering the phosphotyrosine-binding region. These mutations are not predicted to alter the autoinhibition of the PTPN11 protein product, SHP-2, which is the mechanism for the vast majority of mutations causing Noonan syndrome. Future studies will be directed towards understanding these rare phosphotyrosine binding region mutants.

Cathepsin K deficiency in pycnodysostosis results in accumulation of non-digested phagocytosed collagen in fibroblasts.

The rare osteosclerotic disease, pycnodysostosis, is characterized by decreased osteoclastic bone collagen degradation due to the absence of active cathepsin K. Although this enzyme is primarily expressed by osteoclasts, there is increasing evidence that it may also be present in other cells, including fibroblasts. Since fibroblasts are known to degrade collagen intracellularly following phagocytosis, we analyzed various soft connective tissues (periosteum, perichondrium, tendon, and synovial membrane) from a 13-week-old human fetus with pycnodysostosis for changes in this collagen digestion pathway. In addition, the same tissues from cathepsin K-deficient and control mice were analyzed. Microscopic examination of the human fetal tissues showed that cross-banded collagen fibrils had accumulated in lysosomal vacuoles of fibroblasts. Morphometric analysis of periosteal fibroblasts revealed that the volume density of collagen-containing vacuoles was 18 times higher than in fibroblasts of control patients. A similar accumulation was seen in periosteal fibroblasts of three children with pycnodysostosis. In contrast to the findings in humans, an accumulation of internalized collagen was not apparent in fibroblasts of mice with cathepsin K deficiency. Our observations indicate that the intracellular digestion of phagocytosed collagen by fibroblasts is inhibited in humans with pycnodysostosis, but probably not in the mouse model mimicking this disease. The data strongly suggest that cathepsin K is a crucial protease for this process in human fibroblasts. Murine fibroblasts may have other proteolytic activities that are expressed constitutively or up regulated in response to a deficiency of cathepsin K. This may explain why cathepsin K-deficient mice lack the dysostotic features that are prominent in patients with pycnodysostosis.

Exclusion of PTPN11 mutations in Costello syndrome: further evidence for distinct genetic etiologies for Noonan, cardio-facio-cutaneous and Costello syndromes.

Costello syndrome (CS) is a rare, multiple congenital anomaly syndrome with characteristic dysmorphic features, cardiac anomalies and a tendency to develop certain cancers. Phenotypically there is some overlap with other genetic disorders, notably cardio-facio-cutaneous (CFC) syndrome and Noonan syndrome (NS), suggesting that these syndromes may be allelic. We recently identified PTPN11, which encodes the non-receptor protein tyrosine phosphatase, SHP-2, as a major NS disease gene. In this report, we screened a cohort of 27 patients, with the clinical diagnosis of CS, for PTPN11 mutations using denaturing high performance liquid chromatography analysis. No mutations of the PTPN11 gene were found in the CS patients. Common polymorphisms in introns 6 and 7 and exon 8 were identified in four individuals. With our previous exclusion of PTPN11 mutations in CFC syndrome, these data suggest distinct genetic etiologies for Noonan, CFC and Costello syndromes.

Mutations in PTPN11, encoding the protein tyrosine phosphatase SHP-2, cause Noonan syndrome.

Noonan syndrome (MIM 163950) is an autosomal dominant disorder characterized by dysmorphic facial features, proportionate short stature and heart disease (most commonly pulmonic stenosis and hypertrophic cardiomyopathy). Webbed neck, chest deformity, cryptorchidism, mental retardation and bleeding diatheses also are frequently associated with this disease. This syndrome is relatively common, with an estimated incidence of 1 in 1,000-2,500 live births. It has been mapped to a 5-cM region (NS1) [corrected] on chromosome 12q24.1, and genetic heterogeneity has also been documented. Here we show that missense mutations in PTPN11 (MIM 176876)-a gene encoding the nonreceptor protein tyrosine phosphatase SHP-2, which contains two Src homology 2 (SH2) domains-cause Noonan syndrome and account for more than 50% of the cases that we examined. All PTPN11 missense mutations cluster in interacting portions of the amino N-SH2 domain and the phosphotyrosine phosphatase domains, which are involved in switching the protein between its inactive and active conformations. An energetics-based structural analysis of two N-SH2 mutants indicates that in these mutants there may be a significant shift of the equilibrium favoring the active conformation. This implies that they are gain-of-function changes and that the pathogenesis of Noonan syndrome arises from excessive SHP-2 activity.

Novel TFAP2B mutations that cause Char syndrome provide a genotype-phenotype correlation.

To elucidate further the role, in normal development and in disease pathogenesis, of TFAP2B, a transcription factor expressed in neuroectoderm, we studied eight patients with Char syndrome and their families. Four novel mutations were identified, three residing in the basic domain, which is responsible for DNA binding, and a fourth affecting a conserved PY motif in the transactivation domain. Functional analyses of the four mutants disclosed that two, R225C and R225S, failed to bind target sequence in vitro and that all four had dominant negative effects when expressed in eukaryotic cells. Our present findings, combined with data about two previously identified TFAP2B mutations, show that dominant negative effects consistently appear to be involved in the etiology of Char syndrome. Affected individuals in the family with the PY motif mutation, P62R, had a high prevalence of patent ductus arteriosus but had only mild abnormalities of facial features and no apparent hand anomalies, a phenotype different from that associated with the five basic domain mutations. This genotype-phenotype correlation supports the existence of TFAP2 coactivators that have tissue specificity and are important for ductal development but less critical for craniofacial and limb development.

Cloning and characterization of a novel mouse AP-2 transcription factor, AP-2delta, with unique DNA binding and transactivation properties.

AP-2 transcription factors are sequence-specific DNA-binding proteins expressed in neural crest and other tissues during mammalian development. Three mammalian genes, AP-2alpha, AP-2beta, and AP-2gamma, have been reported previously. A partial predicted AP-2 gene was identified in tandem with AP-2beta on human chromosome 6p12-p21.1. The orthologous mouse gene, which we named Ap-2delta, was identified from a fetal mouse head cDNA library. Northern analysis revealed two transcripts in embryonic and newborn mouse brain, with markedly higher steady-state levels in the former. The predicted Ap-2delta protein comprised 452 amino acids and was highly similar to other AP-2 proteins across the DNA-binding and dimerization domains. Ap-2delta formed homodimers and heterodimers in vitro, bound an optimized AP-2 consensus DNA sequence, and transactivated gene expression in eukaryotic cells. Ap-2delta dimers bound poorly to an AP-2 binding sequence from the human metallothionein IIa promoter in vitro, revealing a sequence specificity not previously observed among other AP-2 proteins. The PY motif and critical residues in the transactivation domain, which are highly conserved in the AP-2 family and believed necessary for transactivation, were divergent in Ap-2delta. The unique protein sequence and functional features of Ap-2delta suggest mechanisms, besides tissue-specific AP-2 gene expression, for specific control of target gene activation.

Slc19a2: cloning and characterization of the murine thiamin transporter cDNA and genomic sequence, the orthologue of the human TRMA gene.

Recently, our group and others cloned the TRMA disease gene, SLC19A2, which encodes a thiamin transporter. Here, we report the cloning and characterization of the full-length cDNA and genomic sequences of mouse Slc19a2. The Slc19a2 cDNA contained a 1494-bp open-reading frame, and had 5'- and 3'-untranslated regions of 189 and 1857 bp, respectively. A putative GC-rich, TATA-less promoter was identified in genomic sequence directly upstream of the identified 5' end. The Slc19a2 gene spanned 16.3 kb and was organized into six exons, a gene structure conserved with the human orthologue. The predicted Slc19a2 protein, like SLC19A2, was predicted to have 12 transmembrane domains and shared a number of other conserved sequence motifs with the human orthologue, including one potential N-glycosylation site (N(63)) and several potential phosphorylation sites. Comparison of the Slc19a2 amino acid sequence with those of the other known SLC19A solute carriers highlighted interesting patterns of conservation and divergence in various domains, allowing insight into potential structure-function relationships. The identification of the mouse Slc19a2 cDNA and genomic sequences will facilitate the generation of an animal model of TRMA, permitting future studies of disease pathogenesis.

Genetic basis of syndromes associated with congenital heart disease.

Numerous syndromes affecting patients have phenotypes that include congenital heart defects (CHDs). These disorders have fascinated physicians for many years, raising questions about how seemingly disparate aspects of human development can be perturbed together in striking, but consistent, ways. Paralleling the major advances in human genetics during recent decades, we have come to understand that some of these syndromes arise from gross defects in chromosomal number, some from subtler alterations in genomic regions, and still others from point mutations in specific genes. These disorders, largely mendelian in nature, have provided researchers with the wherewithal to discover disease genes underlying CHD. Although some of these medical conditions are relatively rare, their solution has often provided insights that could be applied toward understanding the basis of nonsyndromic CHD. In this review, recent progress toward uncovering the molecular basis of several forms of syndromic CHD is discussed.

Impact of the reduced folate carrier on the accumulation of active thiamin metabolites in murine leukemia cells.

The thiamin transporter encoded by SLC19A2 and the reduced folate carrier (RFC1) share 40% homology at the protein level, but the thiamin transporter does not mediate transport of folates. By using murine leukemia cell lines that express no, normal, or high levels of RFC1, we demonstrate that RFC1 does not mediate thiamin influx. However, high level RFC1 expression substantially reduced accumulation of the active thiamin coenzyme, thiamin pyrophosphate (TPP). This decreased level of TPP, synthesized intracellularly from imported thiamin, resulted from RFC1-mediated efflux of TPP. This conclusion was supported by the following observations. (i) Efflux of intracellular TPP was increased in cells with high expression of RFC1. (ii) Methotrexate inhibits TPP influx. (iii) TPP competitively inhibits methotrexate influx. (iv) Loading cells, which overexpress RFC1 to high levels of methotrexate to inhibit competitively RFC1-mediated TPP efflux, augment TPP accumulation. (v) There was an inverse correlation between thiamin accumulation and RFC1 activity in cells grown at a physiological concentration of thiamin. The modulation of thiamin accumulation by RFC1 in murine leukemia cells suggests that this carrier may play a role in thiamin homeostasis and could serve as a modifying factor in thiamin nutritional deficiency as well as when the high affinity thiamin transporter is mutated.

Char syndrome: an additional family with polythelia, a new finding.

This report describes a father and daughter with Char syndrome, a rare autosomal dominant disorder. Both affected individuals had typical face, strabismus, and foot anomalies. The girl also had a patent ductus arteriosus. In addition, both patients had polythelia (supernumerary nipples), a finding not described before in the Char syndrome.

Genetics of peanut allergy: a twin study.

BACKGROUND: The role of genetics in the etiology of peanut allergy is unknown. For complex genetic traits, twin studies can provide information on the relative contribution of genetic factors to a disease, as the relative confounding effects of environmental factors are markedly decreased. OBJECTIVE: This study was performed to search for evidence that genetic factors influence peanut allergy by comparing the concordance rate for this allergy among monozygotic and dizygotic twins. METHODS: Twin pairs with at least one member with peanut allergy were ascertained through the Food Allergy Network by advertisements in the organization's newsletters and Web site. Individuals with peanut allergy or parental surrogates were interviewed by telephone. A full atopic history was obtained, and peanut allergy and zygosity were determined using previously validated questionnaires. Heritability of peanut allergy was determined using univariate genetic model fitting by maximum likelihood with the Mx statistical modeling software package. RESULTS: Seventy-five twin pairs were recruited. Seventeen pairs were excluded because of unconvincing peanut allergy histories (9 pairs, including 4 of uncertain zygosity) or because one twin had reportedly never ingested peanut (8 pairs). The median age of the 58 remaining twin pairs was 5 years (range 1 to 58 years). Seventy individuals had peanut allergy. In addition to convincing histories of peanut allergy, 52 (74%) had been tested (skin prick testing with or without radioallergosorbent assay) and all had positive reactions to peanut. Twenty-nine of the 70 had experienced >1 reaction to peanut; 29 of 70 had multisystem reactions. Among the monozygotic pairs (n = 14), 9 were concordant for peanut allergy (pairwise concordance, 64.3%) and among dizygotic pairs (n = 44), 3 were concordant for peanut allergy (pairwise concordance, 6.8%; chi(2) = 21.38, P <.0001). Heritability of peanut allergy was estimated at 81. 6% (95% confidence interval 41.6% to 99.7%) with model fitting using a population prevalence of peanut allergy of 0.4%. CONCLUSIONS: The significantly higher concordance rate of peanut allergy among monozygotic twins suggests strongly that there is a significant genetic influence on peanut allergy.

Mutations in TFAP2B cause Char syndrome, a familial form of patent ductus arteriosus.

Char syndrome is an autosomal dominant trait characterized by patent ductus arteriosus, facial dysmorphism and hand anomalies. Using a positional candidacy strategy, we mapped TFAP2B, encoding a transcription factor expressed in neural crest cells, to the Char syndrome critical region and identified missense mutations altering conserved residues in two affected families. Mutant TFAP2B proteins dimerized properly in vitro, but showed abnormal binding to TFAP2 target sequence. Dimerization of both mutants with normal TFAP2B adversely affected transactivation, demonstrating a dominant-negative mechanism. Our work shows that TFAP2B has a role in ductal, facial and limb development and suggests that Char syndrome results from derangement of neural-crest-cell derivatives.

Gaucher disease: the origins of the Ashkenazi Jewish N370S and 84GG acid beta-glucosidase mutations.

Type 1 Gaucher disease (GD), a non-neuronopathic lysosomal storage disorder, results from the deficient activity of acid beta-glucosidase (GBA). Type 1 disease is panethnic but is more prevalent in individuals of Ashkenazi Jewish (AJ) descent. Of the causative GBA mutations, N370S is particularly frequent in the AJ population, (q approximately .03), whereas the 84GG insertion (q approximately .003) occurs exclusively in the Ashkenazim. To investigate the genetic history of these mutations in the AJ population, short tandem repeat (STR) markers were used to map a 9.3-cM region containing the GBA locus and to genotype 261 AJ N370S chromosomes, 60 European non-Jewish N370S chromosomes, and 62 AJ 84GG chromosomes. A highly conserved haplotype at four markers flanking GBA (PKLR, D1S1595, D1S2721, and D1S2777) was observed on both the AJ chromosomes and the non-Jewish N370S chromosomes, suggesting the occurrence of a founder common to both populations. Of note, the presence of different divergent haplotypes suggested the occurrence of de novo, recurrent N370S mutations. In contrast, a different conserved haplotype at these markers was identified on the 84GG chromosomes, which was unique to the AJ population. On the basis of the linkage disequilibrium (LD) delta values, the non-Jewish European N370S chromosomes had greater haplotype diversity and less LD at the markers flanking the conserved haplotype than did the AJ N370S chromosomes. This finding is consistent with the presence of the N370S mutation in the non-Jewish European population prior to the founding of the AJ population. Coalescence analyses for the N370S and 84GG mutations estimated similar coalescence times, of 48 and 55.5 generations ago, respectively. The results of these studies are consistent with a significant bottleneck occurring in the AJ population during the first millennium, when the population became established in Europe.

Recent advances in the understanding of genetic causes of congenital heart defects.

The clinical approach to children with congenital heart defects (CHD) has been revolutionized during the past four decades by developments in diagnostics and therapeutics. In contrast, a profound understanding of the causes of the majority of CHD has only begun to emerge within the past few years. Prior epidemiological studies suggested that Mendelian disorders constituted a very small percentage of CHD and that polygenic inheritance was responsible for the majority of cases. Recent discoveries, largely achieved with molecular genetic studies, have provided new insights into the genetic basis of heart malformations. These studies have shown that CHD caused by single gene or single locus defects is more common than had been suspected. In addition, a higher percentage of heart malformations occur in the context of familial disease than was evident previously. In this review, molecular genetic studies of specific heart lesions and syndromes with CHD are reviewed. Progress on the Human Genome Project has accelerated identification of genes for Mendelian traits with heart defects, and it is anticipated that disease genes for most single gene traits will be known within a few years. Future challenges include utilizing this emerging genetic information to improve diagnosis and treatment of children with CHD, and harnessing the power of genomics to analyze isolated heart defects with complex inheritance patterns.

Malignant fibrous histiocytoma: inherited and sporadic forms have loss of heterozygosity at chromosome bands 9p21-22-evidence for a common genetic defect.

Hereditary cancers represent a unique opportunity to investigate the genetic etiology of their more common sporadic forms. We recently established genetic linkage for the rare autosomal-dominant bone dysplasia/cancer syndrome, diaphyseal medullary stenosis with malignant fibrous histiocytoma (DMS-MFH), to a 3-cM region on chromosome bands 9p21-22. This hereditary cancer syndrome is characterized by bone infarctions, cortical growth abnormalities, pathologic fractures, and painful debilitation. Most notably, 35% of affected individuals develop bone MFH, a sarcoma that, in its sporadic form, accounts for 6% of all bone cancers. To determine whether the hereditary and sporadic forms of bone MFH are genetically linked, we performed loss of heterozygosity (LOH) studies of the DMS-MFH critical region. In addition to the hereditary specimen, 71% (5/7) of informative sporadic bone MFH specimens displayed LOH for markers within that same region. Definition of the minimal region of LOH overlap effectively limited the DMS-MFH gene to a 2-cM region between markers D9S736 and D9S171. In summary, these studies suggest that a common genetic etiology underlies the autosomal-dominant and sporadic forms of this sarcoma and provide the basis for identifying the putative MFH tumor suppressor gene. Genes Chromosomes Cancer 27:191-195, 2000.