Portal de Periódicos da CAPES

De Novo Missense Variants in TRAF7 Cause Developmental Delay, Congenital Anomalies, and Dysmorphic Features

2018; Elsevier BV; Volume: 103; Issue: 1 Linguagem: Inglês

10.1016/j.ajhg.2018.06.005

1537-6605

Mari Tokita, Chun‐An Chen, David Chitayat, Ellen F. Macnamara, Jill A. Rosenfeld, Neil A. Hanchard, Andrea M. Lewis, Chester Brown, Ronit Marom, Yunru Shao, Danica Novacic, Lynne A. Wolfe, Colleen E. Wahl, Cynthia J. Tifft, Camilo Toro, Jonathan A. Bernstein, Caitlin L. Hale, Julia Silver, Louanne Hudgins, Amitha Ananth, Andrea Hanson‐Kahn, Shirley Shuster, Pilar Magoulas, Vipulkumar Patel, Wenmiao Zhu, Stella M. Chen, Yanjun Jiang, Pengfei Liu, Christine M. Eng, Dominyka Batkovskyte, Alberto di Ronza, Marco Sardiello, Brendan Lee, Christian P. Schaaf, Yaping Yang, Xia Wang,

NF-κB Signaling Pathways

TRAF7 is a multi-functional protein involved in diverse signaling pathways and cellular processes. The phenotypic consequence of germline TRAF7 variants remains unclear. Here we report missense variants in TRAF7 in seven unrelated individuals referred for clinical exome sequencing. The seven individuals share substantial phenotypic overlap, with developmental delay, congenital heart defects, limb and digital anomalies, and dysmorphic features emerging as key unifying features. The identified variants are de novo in six individuals and comprise four distinct missense changes, including a c.1964G>A (p.Arg655Gln) variant that is recurrent in four individuals. These variants affect evolutionarily conserved amino acids and are located in key functional domains. Gene-specific mutation rate analysis showed that the occurrence of the de novo variants in TRAF7 (p = 2.6 × 10−3) and the recurrent de novo c.1964G>A (p.Arg655Gln) variant (p = 1.9 × 10−8) in our exome cohort was unlikely to have occurred by chance. In vitro analyses of the observed TRAF7 mutations showed reduced ERK1/2 phosphorylation. Our findings suggest that missense mutations in TRAF7 are associated with a multisystem disorder and provide evidence of a role for TRAF7 in human development. TRAF7 is a multi-functional protein involved in diverse signaling pathways and cellular processes. The phenotypic consequence of germline TRAF7 variants remains unclear. Here we report missense variants in TRAF7 in seven unrelated individuals referred for clinical exome sequencing. The seven individuals share substantial phenotypic overlap, with developmental delay, congenital heart defects, limb and digital anomalies, and dysmorphic features emerging as key unifying features. The identified variants are de novo in six individuals and comprise four distinct missense changes, including a c.1964G>A (p.Arg655Gln) variant that is recurrent in four individuals. These variants affect evolutionarily conserved amino acids and are located in key functional domains. Gene-specific mutation rate analysis showed that the occurrence of the de novo variants in TRAF7 (p = 2.6 × 10−3) and the recurrent de novo c.1964G>A (p.Arg655Gln) variant (p = 1.9 × 10−8) in our exome cohort was unlikely to have occurred by chance. In vitro analyses of the observed TRAF7 mutations showed reduced ERK1/2 phosphorylation. Our findings suggest that missense mutations in TRAF7 are associated with a multisystem disorder and provide evidence of a role for TRAF7 in human development. Somatic and germline mutations in proto-oncogenes and tumor suppressor genes are well-known causes of cancer. Many proto-oncogenes and tumor suppressor genes also play a crucial role in human development and as such, germline mutations in these genes can lead to developmental disorders.1Tartaglia M. Mehler E.L. Goldberg R. Zampino G. Brunner H.G. Kremer H. van der Burgt I. Crosby A.H. Ion A. Jeffery S. et al.Mutations in PTPN11, encoding the protein tyrosine phosphatase SHP-2, cause Noonan syndrome.Nat. Genet. 2001; 29: 465-468Crossref PubMed Scopus (1330) Google Scholar, 2Yang W. Wang J. Moore D.C. Liang H. Dooner M. Wu Q. Terek R. Chen Q. Ehrlich M.G. Quesenberry P.J. Neel B.G. Ptpn11 deletion in a novel progenitor causes metachondromatosis by inducing hedgehog signalling.Nature. 2013; 499: 491-495Crossref PubMed Scopus (167) Google Scholar, 3Wang X. Charng W.L. Chen C.A. Rosenfeld J.A. Al Shamsi A. Al-Gazali L. McGuire M. Mew N.A. Arnold G.L. Qu C. et al.Germline mutations in ABL1 cause an autosomal dominant syndrome characterized by congenital heart defects and skeletal malformations.Nat. Genet. 2017; 49: 613-617Crossref PubMed Scopus (31) Google Scholar, 4Hanby A.M. Kelsell D.P. Potts H.W. Gillett C.E. Bishop D.T. Spurr N.K. Barnes D.M. Association between loss of heterozygosity of BRCA1 and BRCA2 and morphological attributes of sporadic breast cancer.Int. J. Cancer. 2000; 88: 204-208Crossref PubMed Scopus (25) Google Scholar, 5Howlett N.G. Taniguchi T. Olson S. Cox B. Waisfisz Q. De Die-Smulders C. Persky N. Grompe M. Joenje H. Pals G. et al.Biallelic inactivation of BRCA2 in Fanconi anemia.Science. 2002; 297: 606-609Crossref PubMed Scopus (963) Google Scholar, 6Snijders Blok L. Madsen E. Juusola J. Gilissen C. Baralle D. Reijnders M.R. Venselaar H. Helsmoortel C. Cho M.T. Hoischen A. et al.DDD StudyMutations in DDX3X are a common cause of unexplained intellectual disability with gender-specific effects on Wnt signaling.Am. J. Hum. Genet. 2015; 97: 343-352Abstract Full Text Full Text PDF PubMed Scopus (157) Google Scholar TRAF7 (Tumor necrosis factor receptor-associated factor 7 [MIM: 606692]) belongs to the multi-functional TRAF family and is involved in multiple biological processes, including ubiquitination and myogenesis.7Zotti T. Vito P. Stilo R. The seventh ring: exploring TRAF7 functions.J. Cell. Physiol. 2012; 227: 1280-1284Crossref PubMed Scopus (59) Google Scholar TRAF7 is also a known mediator of the MAP kinase and NF-κB signaling pathways.7Zotti T. Vito P. Stilo R. The seventh ring: exploring TRAF7 functions.J. Cell. Physiol. 2012; 227: 1280-1284Crossref PubMed Scopus (59) Google Scholar Somatic mutations in TRAF7 have been reported in meningioma and mesothelioma.8Zotti T. Scudiero I. Vito P. Stilo R. The emerging role of TRAF7 in tumor development.J. Cell. Physiol. 2017; 232: 1233-1238Crossref PubMed Scopus (35) Google Scholar Although germline TRAF7 variants have been reported in two individuals with autism spectrum disorder,9Krumm N. Turner T.N. Baker C. Vives L. Mohajeri K. Witherspoon K. Raja A. Coe B.P. Stessman H.A. He Z.X. et al.Excess of rare, inherited truncating mutations in autism.Nat. Genet. 2015; 47: 582-588Crossref PubMed Scopus (365) Google Scholar, 10Neale B.M. Kou Y. Liu L. Ma’ayan A. Samocha K.E. Sabo A. Lin C.F. Stevens C. Wang L.S. Makarov V. et al.Patterns and rates of exonic de novo mutations in autism spectrum disorders.Nature. 2012; 485: 242-245Crossref PubMed Scopus (1282) Google Scholar a clear association between germline defects in TRAF7 and human disease has not been established. Herein we report missense variants in TRAF7 in seven unrelated individuals with overlapping features. This study was performed according to the standards of the ethics committees and the institutional review boards at Baylor College of Medicine, the National Institutes of Health, and Stanford University. Written informed consent was obtained from all study participants. Medical records were reviewed and clinical information was extracted for phenotype analysis. Key clinical features are summarized in Table 1 and are discussed below; comprehensive phenotype information is available in Table S1. Select photographs of the described subjects are presented in Figure 1. Clinical exome sequencing was performed as previously described;11Yang Y. Muzny D.M. Reid J.G. Bainbridge M.N. Willis A. Ward P.A. Braxton A. Beuten J. Xia F. Niu Z. et al.Clinical whole-exome sequencing for the diagnosis of mendelian disorders.N. Engl. J. Med. 2013; 369: 1502-1511Crossref PubMed Scopus (1420) Google Scholar sequencing metrics are summarized in Table S2.Table 1Key Clinical Characteristics of Subjects with TRAF7 VariantsSubject 1Subject 2Subject 3Subject 4Subject 5Subject 6Subject 7Age7.5 years3 weeks43 years28 months1 week20 months8 yearsSexmalemalemalefemalemalefemalemaleVariantp.Arg655Glnp.Arg655Glnp.Arg655Glnp.Arg655Glnp.Thr601Alap.Lys346Glup.Arg371GlyNumber of mutant/total reads by exome sequencing98/20756/107188/39056/120192/41175/13614/105Prenatal ultrasound findingschoroid plexus cyst, two-vessel cordcystic hygroma, cardiac defect, two-vessel cordnonedichorionic/diamniotic twin pregnancy, cardiac defectcardiac defect, sloped forehead, redundant nuchal skin, bell-shaped chestnonecystic hygroma, IUGRPostnatal concernsrespiratory distress, poor feeding, hypotoniarespiratory distresshypersomnia, poor feedingrespiratory distresscardiopulmonary arrest; ECMO complicated by intracranial hemorrhagerespiratory distress with pulmonary edema, jaundiceincreased sleepiness, poor feeding, jaundiceDevelopment Motor delayyesn/ayesyesn/anonormal then regressed Speech delayyesn/ayesyesn/ayesnormal then regressed, severe autismNeurologic features Tonehyper/hypotoniahypertonialeg stiffnessnormalhypotonianormalnormal Seizuresnonoyesnonostaring spells, possible absence seizuresyesDysmorphic features Multiple hair whorlsyesnononoyesnoyes Ptosisyesyesnoyesnonoyes Epicanthal foldsyesyesnoyesnoyesyes Abnormally set or dysplastic earsyesyesyesyesyesyesno Low hairline or excess nuchal skinyesyesnoyesyesnoyes Wide-spaced or inverted nipplesyesyesnoyesyesnoyes Umbilical hernia or diastasis rectiyesyesnoyesnononoStructural findings on neuroimaginghydrocephalus, variant venous anatomy, abnormally shaped corpus callosum and cerebellummildly prominent ventricles and extra cerebral spacessyringohydromyelia of C5-T12cerebral greater than cerebellar atrophynormalpossible sulcal/gyral asymmetry in right frontal/parietal region, mild cortical irregularityprominent anterior and posterior horns of left lateral ventricle, variant venous anatomyEchocardiogram findingssupravalvular pulmonary stenosis, PDA; incomplete right bundle branch blocksevere aortic coarctation, parachute mitral valve, supravalvular mitral ring, PDA, small VSD, hypoplastic left ventriclebicuspid aortic valve, small ASD, PDADORV, pulmonary atresia, hypoplastic left heart, aortopulmonary collaterals, RVHDORV, mitral atresia, univentricular physiologybicuspid aortic valve, moderate aortic valve stenosis, large PDA; prolonged QTnoneDigital anomalies Overlapping digitsyesyesyesnoyesyesno Variant creasesyesyesnoyesyesyesyes Clinodactylynoyesnonoyesnoyes ContracturesyesnoyesnonononoAbbreviations: IUGR, intrauterine growth restriction; n/a, not ascertained; PDA, patent ductus arteriosus; ASD, atrial septal defect; VSD, ventricular septal defect; DORV, double outlet right ventricle; RVH, right ventricular hypertrophy; ECMO, extracorporeal membrane oxygen. Open table in a new tab Abbreviations: IUGR, intrauterine growth restriction; n/a, not ascertained; PDA, patent ductus arteriosus; ASD, atrial septal defect; VSD, ventricular septal defect; DORV, double outlet right ventricle; RVH, right ventricular hypertrophy; ECMO, extracorporeal membrane oxygen. The study cohort consisted of two female and five male subjects ranging in age from 1 week to 43 years. Prenatal histories were notable for antenatal detection of heart anomalies (n = 3), cystic hygroma (n = 2), and two-vessel cord (n = 2). Postnatally, several patients had respiratory distress (n = 4) and poor feeding (n = 3). Congenital heart defects were present in six of seven patients and ranged in type and severity (Tables 1 and S1). Motor and/or speech delay were present to a variable degree in five of the five subjects for whom this outcome could be assessed. Subjects 3 and 7 had documented seizures, and subject 6 had a history of staring spells concerning for absence seizures. Brain imaging revealed variant venous anatomy in two subjects and prominent/enlarged ventricles and/or cerebral atrophy in four subjects (Table S1). Subject 6 was reported to have sulcal/gyral asymmetry of the right frontal parietal region with mild cortical irregularity concerning for possible polymicrogyria. A tethered cord was found in subjects 1, 3, and 4. Sensorineural and mixed hearing loss were reported in subjects 1 and 3, respectively. Cortical blindness was reported in subject 2 and optic atrophy was reported in subjects 2 and 3. Shared dysmorphic features included epicanthal folds (n = 5), ptosis (n = 4), abnormally set or dysplastic ears (n = 6), and widely spaced and/or inverted nipples (n = 5) (Figure 1). A low hairline and/or excess nuchal skin were reported in five subjects, two of whom had a prenatal history of a cystic hygroma. Three subjects were found to have multiple hair whorls. Limb and digital anomalies were present in all seven subjects to a variable degree, including variant palmar/digital crease patterns (n = 6) and overlapping toes (n = 5). Subject 3, the eldest of the cohort, had relatively extensive involvement of the musculoskeletal system including a history of scoliosis with degenerative joint disease, an asymmetric sternum, bilateral avascular necrosis of the hip, tibial malformation, valgus deformity of the ankles, flexion contractures of the hands (Figure 1C), and subluxation of multiple joints in the hands and feet. In summary, the seven individuals described herein share similar clinical features including developmental delay (5/5), congenital heart defects (6/7), limb and digital anomalies (7/7), and dysmorphic facial features (7/7). TRAF7 is a 21-exon gene located at chromosomal region 16p13.3. The longest isoform of the TRAF7 protein (GenBank: NP_115647) is 670 amino acids in length and contains an N-terminal ring finger domain, an adjacent zinc-finger domain, a centrally situated coiled-coil motif, and seven WD40 repeats in the C-terminal domain (Figures 2B and 2C ).7Zotti T. Vito P. Stilo R. The seventh ring: exploring TRAF7 functions.J. Cell. Physiol. 2012; 227: 1280-1284Crossref PubMed Scopus (59) Google Scholar, 12Xu L.G. Li L.Y. Shu H.B. TRAF7 potentiates MEKK3-induced AP1 and CHOP activation and induces apoptosis.J. Biol. Chem. 2004; 279: 17278-17282Abstract Full Text Full Text PDF PubMed Scopus (134) Google Scholar Clinical exome sequencing identified a total of four distinct missense variants in TRAF7 (GenBank: NM_032271.2) in the seven described subjects, including a recurrent c.1964G>A (p.Arg655Gln) variant found in four subjects. All available parental samples were analyzed by exome and/or Sanger sequencing to assess inheritance of the TRAF7 variants. The variants occurred de novo in six subjects (subjects 1–3 and 5–7). For subject 4, a maternal sample was negative for the variant but a paternal sample was not available for testing. The four TRAF7 variants detected in our cohort are absent from the ExAC and gnomAD databases (access date 5/1/2018) and alter highly phylogenetically conserved amino acids (Figure 2A). The in silico variant effect prediction tool MutationTaster predicts all four variants to be disease-causing,13Schwarz J.M. Cooper D.N. Schuelke M. Seelow D. MutationTaster2: mutation prediction for the deep-sequencing age.Nat. Methods. 2014; 11: 361-362Crossref PubMed Scopus (2441) Google Scholar while SIFT and PolyPhen-2 yield consistent results with two exceptions (Table 2).14Adzhubei I. Jordan D.M. Sunyaev S.R. Predicting functional effect of human missense mutations using PolyPhen-2.Curr. Protoc. Hum. Genet. 2013; Chapter 7: Unit7 20PubMed Google Scholar, 15Kumar P. Henikoff S. Ng P.C. Predicting the effects of coding non-synonymous variants on protein function using the SIFT algorithm.Nat. Protoc. 2009; 4: 1073-1081Crossref PubMed Scopus (5027) Google Scholar Subjects 1–4 had a c.1964G>A (p.Arg655Gln) variant in exon 20 that localized to the seventh of seven WD40 domains at the C terminus of TRAF7 (Figures 2B and 2C). The c.1964G>A (p.Arg655Gln) variant has been previously reported as a de novo event in a female with autism spectrum disorder.9Krumm N. Turner T.N. Baker C. Vives L. Mohajeri K. Witherspoon K. Raja A. Coe B.P. Stessman H.A. He Z.X. et al.Excess of rare, inherited truncating mutations in autism.Nat. Genet. 2015; 47: 582-588Crossref PubMed Scopus (365) Google Scholar The recurrent nature of this variant is likely explained by its location within a CpG dinucleotide. Subject 5 had a c.1801A>G (p.Thr601Ala) variant which localized to exon 19 and the sixth WD40 domain. Subject 6 had a c.1036A>G (p.Lys346Glu) variant found in exon 11 affecting an amino acid in the coiled-coiled motif. Finally, subject 7 had a mosaic c.1111C>G (p.Arg371Gly) variant in exon 12 that localized to the coiled-coiled domain. The c.1111C>G (p.Arg371Gly) variant had an allele fraction (AF) of 13% (14 mutant and 105 total reads) by exome sequencing using DNA extracted from the blood of this individual. As a confirmatory study, we used the restriction enzyme, HhaI, to digest the reference allele in both proband and parental samples then Sanger sequenced the digested product. This analysis demonstrated unequivocal presence of the c.1111C>G allele in the proband but not the parental samples, confirming the presence of the mosaic variant in the proband only (Figure S1). We further performed PCR followed by high-throughput sequencing using DNA extracted from the blood of this individual and the parents (Supplemental Material and Methods). We found that the c.1111C>G variant had an AF of 15% at 23,108× coverage in the proband and an AF of 0% at 22,523× and 22,957× coverage in the mother and father, respectively. Notably, the phenotype of subject 7 is not apparently milder than the other six individuals. This may be explained by variability in the levels of mosaicism in critical tissues.Table 2TRAF7 Variants Detected in the CohortSubjects 1–4Subject 5Subject 6Subject 7DNA variantc.1964G>Ac.1801A>Gc.1036A>Gc.1111C>GProtein changep.Arg655Glnp.Thr601Alap.Lys346Glup.Arg371GlyGenomic coordinatechr16:2226351chr16:2226104chr16:2223505chr16:2223813Inheritancede novoaThis variant was de novo in subjects 1–3, and was negative in the mother of subject 4 whose paternal sample was not available.de novode novode novoAA conservationhighhighhighhighMutationTasterdisease-causingdisease-causingdisease-causingdisease-causingSIFTdamagingtolerateddamagingdamagingPolyPhen-2probably damagingprobably damagingbenignpossibly damagingUbiquitination sitenonoyesnoNomenclature is based on transcript GenBank: NM_032271.2. hg19 genomic coordinates are used.a This variant was de novo in subjects 1–3, and was negative in the mother of subject 4 whose paternal sample was not available. Open table in a new tab Nomenclature is based on transcript GenBank: NM_032271.2. hg19 genomic coordinates are used. Using different deletion mutations, several groups have clarified the functional relevance of the recognizable TRAF7 protein domains. The C-terminal WD40 domain has been shown to facilitate TRAF7 interaction with MAP3K3 (MIM: 602539) while TRAF7 homodimerization and subcellular localization appear to rely on the presence of the centrally located coiled-coiled and zinc finger motifs.12Xu L.G. Li L.Y. Shu H.B. TRAF7 potentiates MEKK3-induced AP1 and CHOP activation and induces apoptosis.J. Biol. Chem. 2004; 279: 17278-17282Abstract Full Text Full Text PDF PubMed Scopus (134) Google Scholar, 16Bouwmeester T. Bauch A. Ruffner H. Angrand P.O. Bergamini G. Croughton K. Cruciat C. Eberhard D. Gagneur J. Ghidelli S. et al.A physical and functional map of the human TNF-alpha/NF-kappa B signal transduction pathway.Nat. Cell Biol. 2004; 6: 97-105Crossref PubMed Scopus (869) Google Scholar Of the four TRAF7 mutations detected in our cohort, two are found within the coiled-coiled domain and two in the WD40 repeat region (Figure 2C). Additionally, the p.Lys346Glu variant detected in subject 6 affects an experimentally determined ubiquitination site.17Kim W. Bennett E.J. Huttlin E.L. Guo A. Li J. Possemato A. Sowa M.E. Rad R. Rush J. Comb M.J. et al.Systematic and quantitative assessment of the ubiquitin-modified proteome.Mol. Cell. 2011; 44: 325-340Abstract Full Text Full Text PDF PubMed Scopus (1183) Google Scholar Finally, TRAF7 is predicted to be intolerant to missense variations (z = 3.34).18Lek M. Karczewski K.J. Minikel E.V. Samocha K.E. Banks E. Fennell T. O’Donnell-Luria A.H. Ware J.S. Hill A.J. Cummings B.B. et al.Exome Aggregation ConsortiumAnalysis of protein-coding genetic variation in 60,706 humans.Nature. 2016; 536: 285-291Crossref PubMed Scopus (6602) Google Scholar Thus, the localization of the detected TRAF7 variants to highly conserved amino acids in key functional domains and the predicted intolerance of TRAF7 to missense variation support the notion that missense mutations in TRAF7 may be causative of a human phenotype. We sought to evaluate the likelihood of detecting de novo variants in TRAF7 and the recurrent de novo c.1964G>A change in our exome database by chance. To do so, we performed a gene-specific mutation rate analysis using a previously published approach.19Allen A.S. Berkovic S.F. Cossette P. Delanty N. Dlugos D. Eichler E.E. Epstein M.P. Glauser T. Goldstein D.B. Han Y. et al.Epi4K ConsortiumEpilepsy Phenome/Genome ProjectDe novo mutations in epileptic encephalopathies.Nature. 2013; 501: 217-221Crossref PubMed Scopus (1093) Google Scholar In total there were about 7,700 individuals in our undiagnosed exome cohort at the time of analysis. Focusing on rare variants (population frequency < 1% in gnomAD database) in this cohort, we identified seven individuals with de novo variants in TRAF7 (including subjects 1–3, 5–7, and one additional individual for whom we did not obtain research consent) and three individuals with the recurrent de novo c.1964G>A (p.Arg655Gln) variant (including subjects 1–3 and excluding subject 4 for whom the paternal sample was not available). We used 2.4 × 10−5 as the expected TRAF7 gene-specific rate of de novo missense mutation, which was previously determined based on gene length, local sequence context, and mutation type.20Samocha K.E. Robinson E.B. Sanders S.J. Stevens C. Sabo A. McGrath L.M. Kosmicki J.A. Rehnström K. Mallick S. Kirby A. et al.A framework for the interpretation of de novo mutation in human disease.Nat. Genet. 2014; 46: 944-950Crossref PubMed Scopus (615) Google Scholar After adjusting the p value for 18,892 genes in the Consensus Coding Sequence for genome-wide statistical significance, the likelihood to detect seven individuals with de novo TRAF7 variants and three individuals with the recurrent de novo c.1964G>A (p.Arg655Gln) variant in our cohort was 2.6 × 10−3 and 1.9 × 10−8, respectively (Table S3). These data indicate that our undiagnosed exome cohort is enriched for de novo variants in TRAF7 and renders improbable the detection of seven de novo TRAF7 variants in a cohort this size by chance. Seven TRAF (TNF receptor associated factor) proteins (TRAF1–7) have been characterized in mammals to date and comprise the TRAF family of proteins.7Zotti T. Vito P. Stilo R. The seventh ring: exploring TRAF7 functions.J. Cell. Physiol. 2012; 227: 1280-1284Crossref PubMed Scopus (59) Google Scholar TRAF proteins are intracellular adaptors that bind TNF family receptors and transduce the cellular effects mediated by TNF family ligands.7Zotti T. Vito P. Stilo R. The seventh ring: exploring TRAF7 functions.J. Cell. Physiol. 2012; 227: 1280-1284Crossref PubMed Scopus (59) Google Scholar, 21Dempsey P.W. Doyle S.E. He J.Q. Cheng G. The signaling adaptors and pathways activated by TNF superfamily.Cytokine Growth Factor Rev. 2003; 14: 193-209Crossref PubMed Scopus (401) Google Scholar Mitogen-activated protein kinases (MAPKs) comprise a family of protein kinases that regulate several cellular activities including gene expression, mitosis, movement, metabolism, and programmed death.22Johnson G.L. Lapadat R. Mitogen-activated protein kinase pathways mediated by ERK, JNK, and p38 protein kinases.Science. 2002; 298: 1911-1912Crossref PubMed Scopus (3502) Google Scholar There are three subfamilies of MAPKs in multicellular organisms: ERK1 and ERK2; c-Jun NH2-terminal kinases JNK; and p38.22Johnson G.L. Lapadat R. Mitogen-activated protein kinase pathways mediated by ERK, JNK, and p38 protein kinases.Science. 2002; 298: 1911-1912Crossref PubMed Scopus (3502) Google Scholar Although the phenotypic effects of variants in TRAF7 have not been clearly characterized in humans or in model organisms to our knowledge, a role for TRAF7 in MAPKs pathways has been clearly established. Bouwmeester et al. and Xu et al. showed that TRAF7 interacts by its WD40 domain with MAP3K3 and that coexpression of TRAF7 with MAP3K3 activates the JNK and p38 signaling pathways that affect cell proliferation, differentiation, survival, and migration.12Xu L.G. Li L.Y. Shu H.B. TRAF7 potentiates MEKK3-induced AP1 and CHOP activation and induces apoptosis.J. Biol. Chem. 2004; 279: 17278-17282Abstract Full Text Full Text PDF PubMed Scopus (134) Google Scholar, 16Bouwmeester T. Bauch A. Ruffner H. Angrand P.O. Bergamini G. Croughton K. Cruciat C. Eberhard D. Gagneur J. Ghidelli S. et al.A physical and functional map of the human TNF-alpha/NF-kappa B signal transduction pathway.Nat. Cell Biol. 2004; 6: 97-105Crossref PubMed Scopus (869) Google Scholar, 23Wagner E.F. Nebreda A.R. Signal integration by JNK and p38 MAPK pathways in cancer development.Nat. Rev. Cancer. 2009; 9: 537-549Crossref PubMed Scopus (1845) Google Scholar The ERK1 and ERK2 signaling cascades have not been directly linked to TRAF7, but both are activated by MAP2K1/MAP2K2 (MEK1/2)24Wortzel I. Seger R. The ERK cascade: distinct functions within various subcellular organelles.Genes Cancer. 2011; 2: 195-209Crossref PubMed Scopus (337) Google Scholar which in turn are activated by multiple protein kinases including MAP3K325Ellinger-Ziegelbauer H. Brown K. Kelly K. Siebenlist U. Direct activation of the stress-activated protein kinase (SAPK) and extracellular signal-regulated protein kinase (ERK) pathways by an inducible mitogen-activated protein Kinase/ERK kinase kinase 3 (MEKK) derivative.J. Biol. Chem. 1997; 272: 2668-2674Abstract Full Text Full Text PDF PubMed Scopus (90) Google Scholar, 26Caunt C.J. Sale M.J. Smith P.D. Cook S.J. MEK1 and MEK2 inhibitors and cancer therapy: the long and winding road.Nat. Rev. Cancer. 2015; 15: 577-592Crossref PubMed Scopus (372) Google Scholar—a known TRAF7 binding partner. Collectively, these lines of evidence suggest that the TRAF7 variants detected in our cohort may interfere with TRAF7 regulation of MAPK pathways. To characterize the functional consequence of the TRAF7 variants detected in the described subjects, we performed site-directed mutagenesis to introduce the four observed variants into TRAF7 cDNA, which we then transfected into HEK293T cells (Supplemental Material and Methods). Expression of mutant TRAF7 had no apparent effect on total or phosphorylated p38 and JNK when assessed by immunoblotting (data not shown). Cellular levels of ERK1/2, mediators of a third MAPK cascade, were not affected by mutant TRAF7, but phosphorylation of ERK1/2 was reduced in mutant-transfected cells relative to wild-type (Figure 3A). This reduction was statistically significant with TNFα stimulation for cells expressing the p.Arg655Gln, p.Thr601Ala, or p.Arg371Gly variants (Figures 3D and 3E). Cells expressing the p.Lys346Glu variant demonstrated a consistent but not statistically significant reduction in phosphorylated ERK1/2 (Figures 3B–3E). ERK1/2, also known as MAPK3/MAPK1, act downstream of Ras24Wortzel I. Seger R. The ERK cascade: distinct functions within various subcellular organelles.Genes Cancer. 2011; 2: 195-209Crossref PubMed Scopus (337) Google Scholar and mobilize a signaling cascade involved in cell proliferation, differentiation, and survival.27Frémin C. Saba-El-Leil M.K. Lévesque K. Ang S.L. Meloche S. Functional redundancy of ERK1 and ERK2 MAP kinases during development.Cell Rep. 2015; 12: 913-921Abstract Full Text Full Text PDF PubMed Scopus (57) Google Scholar Functional redundancy between ERK1 and ERK2 has been demonstrated in mice, but ERK1 knockout mice are viable while ERK2 knockout mice die by E11.5 due to defective placental development.28Hatano N. Mori Y. Oh-hora M. Kosugi A. Fujikawa T. Nakai N. Niwa H. Miyazaki J. Hamaoka T. Ogata M. Essential role for ERK2 mitogen-activated protein kinase in placental development.Genes Cells. 2003; 8: 847-856Crossref PubMed Scopus (230) Google Scholar Conditional knockout of ERK2 in the embryo while sparing the placenta results in embryonic survival to E16.5.27Frémin C. Saba-El-Leil M.K. Lévesque K. Ang S.L. Meloche S. Functional redundancy of ERK1 and ERK2 MAP kinases during development.Cell Rep. 2015; 12: 913-921Abstract Full Text Full Text PDF PubMed Scopus (57) Google Scholar Interestingly, ERK2 conditional knockout mice show a phenotype highly consistent with what we have observed in our subjects with TRAF7 variants, including craniofacial abnormalities, cardiovascular malformations, and limb defects.27Frémin C. Saba-El-Leil M.K. Lévesque K. Ang S.L. Meloche S. Functional redundancy of ERK1 and ERK2 MAP kinases during development.Cell Rep. 2015; 12: 913-921Abstract Full Text Full Text PDF PubMed Scopus (57) Google Scholar Additionally, conditional knockout or knockdown of ERK2 in the mouse central nervous system leads to reduced cortical thickness, profound deficits in associative learning, abnormal behaviors related to autism spectrum disorder and anxiety, and impaired long-term memory.29Samuels I.S. Karlo J.C. Faruzzi A.N. Pickering K. Herrup K. Sweatt J.D. Saitta S.C. Landreth G.E. Deletion of ERK2 mitogen-activated protein kinase identifies its key roles in cortical neurogenesis and cognitive function.J. Neurosci. 2008; 28: 6983-6995Crossref PubMed Scopus (206) Google Scholar, 30Satoh Y. Endo S. Nakata T. Kobayashi Y. Yamada K. Ikeda T. Takeuchi A. Hiramoto T. Watanabe Y. Kazama T. ERK2 contributes to the control of social behaviors in mice.J. Neurosci. 2011; 31: 11953-11967Crossref PubMed Scopus (109) Google Scholar, 31Satoh Y. Endo S. Ikeda T. Yamada K. Ito M. Kuroki M. Hiramoto T. Imamura O. Kobayashi Y. Watanabe Y. et al.Extracellular signal-regulated kinase 2 (ERK2) knockdown mice show deficits in long-term memory; ERK2 has a specific function in learning and memory.J. Neurosci. 2007; 27: 10765-10776Crossref PubMed Scopus (138) Google Scholar Thus, the mouse models support that reduced ERK1/2 activation resulting from mutations in TRAF7 represents a plausible mechanism of disease in the described subjects. Further study is needed to illustrate the mechanism of TRAF7 regulation on ERK1/2 and MAPKs pathways. Interestingly, subject 1 had been evaluated for a possible diagnosis of cardio-facio-cutaneous syndrome and subject 2 was suspected of having a Noonan-spectrum condition prior to the detection of the TRAF7 variants by exome sequencing (Table S1). On further review, we noted that subject 3 has features reminiscent of Costello syndrome including a coarse facial appearance, congenital heart defect, and musculoskeletal involvement including scoliosis and contractures. Further studies are needed to illustrate the overlapping and discordant phenotypic features and the possible mechanistic link between RASopathy and TRAF7-related disorder. Our results, along with the studies in human32Tidyman W.E. Rauen K.A. Pathogenetics of the RASopathies.Hum. Mol. Genet. 2016; 25: R123-R132Crossref PubMed Scopus (77) Google Scholar, 33Rauen K.A. The RASopathies.Annu. Rev. Genomics Hum. Genet. 2013; 14: 355-369Crossref PubMed Scopus (544) Google Scholar and mice, suggest that the MAPK pathways need to be tightly regulated during development and that dysregulation can cause developmental disorders. Heterozygous missense variants in TRAF7 have been ascertained as somatically acquired driver mutations in neoplastic tissue. Exome sequencing performed on intraneural perineuromas uncovered TRAF7 variants in the WD40 domain.34Klein C.J. Wu Y. Jentoft M.E. Mer G. Spinner R.J. Dyck P.J. Dyck P.J. Mauermann M.L. Genomic analysis reveals frequent TRAF7 mutations in intraneural perineuriomas.Ann. Neurol. 2017; 81: 316-321Crossref PubMed Scopus (40) Google Scholar Clark et al. sequenced 300 meningiomas and found either missense or in-frame indel mutations in TRAF7 in 25% of samples, predominantly in the WD40 domain.35Clark V.E. Erson-Omay E.Z. Serin A. Yin J. Cotney J. Ozduman K. Avşar T. Li J. Murray P.B. Henegariu O. et al.Genomic analysis of non-NF2 meningiomas reveals mutations in TRAF7, KLF4, AKT1, and SMO.Science. 2013; 339: 1077-1080Crossref PubMed Scopus (562) Google Scholar Goode et al. identified recurrent TRAF7 variants in adenomatoid tumors of the genital tract.36Goode B. Joseph N.M. Stevers M. Van Ziffle J. Onodera C. Talevich E. Grenert J.P. Yeh I. Bastian B.C. Phillips J.J. et al.Adenomatoid tumors of the male and female genital tract are defined by TRAF7 mutations that drive aberrant NF-kB pathway activation.Mod. Pathol. 2018; 31: 660-673Crossref PubMed Scopus (47) Google Scholar Thus, the nature and location of TRAF7 variants detected in tumor tissue mirror the TRAF7 germline events seen in our subjects. Intriguingly, the recurrent TRAF7 variants in somatic cancer, such as c.1562A>G (p.His521Arg) and c.1683C>A (p.Ser561Arg), have not been seen in our cohort or in the gnomAD database, possibly because these variants may be lethal if present as heterozygous constitutional variants. Notably, subject 3, the eldest patient in our cohort, has been diagnosed with a meningioma. Given that somatic changes in TRAF7 have been reported in patients with meningioma and mesothelioma,8Zotti T. Scudiero I. Vito P. Stilo R. The emerging role of TRAF7 in tumor development.J. Cell. Physiol. 2017; 232: 1233-1238Crossref PubMed Scopus (35) Google Scholar additional long-term follow-up of the younger subjects will be necessary to definitively determine whether germline TRAF7 variants predispose to meningioma or other tumor formation. In conclusion, we propose that germline missense variants in TRAF7 are associated with a human phenotype characterized by developmental delay, heart defects, and dysmorphic features. The shared phenotypes among the patients, statistically significant enrichment of de novo TRAF7 variants in our exome cohort, high phylogenetic conservation of the altered amino acids, localization of variants to key functional domains, and reduced ERK1/2 activation resulting from the TRAF7 variants support the putative pathogenicity of these variants. The Department of Molecular and Human Genetics at Baylor College of Medicine derives revenue from molecular genetic testing offered at Baylor Genetics. The authors sincerely thank all the family members for their participation in this study. This work was supported by the departmental fund from the Department of Molecular and Human Genetics at Baylor College of Medicine. This work was supported in part by the Intramural Research Program of the NHGRI and the Common Fund, Office of the Director, NIH. R.M. is supported by the Osteogenesis Imperfecta Foundation Michael Geisman Fellowship. Download .pdf (.46 MB) Help with pdf files Document S1. Supplemental Material and Methods, Figure S1, and Tables S2–S4 Download .xlsx (.02 MB) Help with xlsx files Table S1. Comprehensive Clinical Characteristics of Subjects with TRAF7 Variants Consensus Coding Sequence, https://www.ncbi.nlm.nih.gov/projects/CCDS/CcdsBrowse.cgi?REQUEST=SHOW_STATISTICS#Current_Homo_sapiens_20dbSNP database, https://www.ncbi.nlm.nih.gov/projects/SNP/ExAC Browser, http://exac.broadinstitute.org/GenBank, https://www.ncbi.nlm.nih.gov/genbank/gnomAD Browser, http://gnomad.broadinstitute.org/Mutation Taster, http://www.mutationtaster.org/OMIM, http://www.omim.org/PolyPhen-2, http://genetics.bwh.harvard.edu/pph2/SIFT, http://sift.jcvi.org/UCSC Genome Browser hg19, http://genome.ucsc.edu