Exome Aggregation Consortium Research Articles

Assessment of PARP4 as a candidate breast cancer susceptibility gene.

PARP4 has been proposed as a candidate breast cancer susceptibility gene. However, its function and involvement in breast carcinogenesis is unclear. We sought to determine the variant frequency of PARP4 in BRCA-negative women referred for genetic testing from Singapore and to perform functional analyses of PARP4. Next-generation sequencing of PARP4 was conducted for 198 BRCA-negative cases from Singapore. Three independent case-control association analyses of PARP4 were performed for (1) our Singaporean cohort, (2) three dbGaP datasets, and (3) cases from TCGA, with controls from the Exome Aggregation Consortium (ExAC). PARP4 knockout cells were generated utilizing the CRISPR-Cas9 approach in MDA-MB-231 (breast cancer) and MCF10A (normal breast) cell lines, and colony formation, cell proliferation, and migration assays carried out. Candidate variants in PARP4 were identified in 5.5% (11/198) of our Singapore cohort. Case-control association studies for our cases and the dbGaP datasets showed no significant association. However, a significant association was observed for PARP4 variants when comparing 988 breast cancer cases from the TCGA provisional data and 53,105 controls from ExAC (ALL) (OR 0.249, 95% CI 0.139-0.414, P = 2.86 × 10-11). PARP4 knockout did not affect the clonogenicity, proliferation rate, and migration of normal breast cells, but appeared to decrease the proliferation rate and clonogenicity of breast cancer cells. Taken together, our results do not support that PARP4 functions as a cancer susceptibility gene. This study highlights the importance of performing functional analyses for candidate cancer predisposition genes.

MetaDome: Pathogenicity analysis of genetic variants through aggregation of homologous human protein domains.

The growing availability of human genetic variation has given rise to novel methods of measuring genetic tolerance that better interpret variants of unknown significance. We recently developed a concept based on protein domain homology in the human genome to improve variant interpretation. For this purpose, we mapped population variation from the Exome Aggregation Consortium (ExAC) and pathogenic mutations from the Human Gene Mutation Database (HGMD) onto Pfam protein domains. The aggregation of these variation data across homologous domains into meta‐domains allowed us to generate amino acid resolution of genetic intolerance profiles for human protein domains. Here, we developed MetaDome, a fast and easy‐to‐use web server that visualizes meta‐domain information and gene‐wide profiles of genetic tolerance. We updated the underlying data of MetaDome to contain information from 56,319 human transcripts, 71,419 protein domains, 12,164,292 genetic variants from gnomAD, and 34,076 pathogenic mutations from ClinVar. MetaDome allows researchers to easily investigate their variants of interest for the presence or absence of variation at corresponding positions within homologous domains. We illustrate the added value of MetaDome by an example that highlights how it may help in the interpretation of variants of unknown significance. The MetaDome web server is freely accessible at https://stuart.radboudumc.nl/metadome.

Germline loss-of-function variants in the BARD1 gene are associated with early-onset familial breast cancer but not ovarian cancer

BackgroundThe role of the BARD1 gene in breast cancer (BC) and ovarian cancer (OC) predisposition remains elusive, as published case-control investigations have revealed controversial results. We aimed to assess the role of deleterious BARD1 germline variants in BC/OC predisposition in a sample of 4920 BRCA1/2-negative female BC/OC index patients of the German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC).MethodsA total of 4469 female index patients with BC, 451 index patients with OC, and 2767 geographically matched female control individuals were screened for loss-of-function (LoF) mutations and potentially damaging rare missense variants in BARD1. All patients met the inclusion criteria of the GC-HBOC for germline testing and reported at least one relative with BC or OC. Additional control datasets (Exome Aggregation Consortium, ExAC; Fabulous Ladies Over Seventy, FLOSSIES) were included for the calculation of odds ratios (ORs).ResultsWe identified LoF variants in 23 of 4469 BC index patients (0.51%) and in 36 of 37,265 control individuals (0.10%), resulting in an OR of 5.35 (95% confidence interval [CI] = 3.17–9.04; P < 0.00001). BARD1-mutated BC index patients showed a significantly younger mean age at first diagnosis (AAD; 42.3 years, range 24–60 years) compared with the overall study sample (48.6 years, range 17–92 years; P = 0.00347). In the subgroup of BC index patients with an AAD < 40 years, an OR of 12.04 (95% CI = 5.78–25.08; P < 0.00001) was observed. An OR of 7.43 (95% CI = 4.26–12.98; P < 0.00001) was observed when stratified for an AAD < 50 years. LoF variants in BARD1 were not significantly associated with BC in the subgroup of index patients with an AAD ≥ 50 years (OR = 2.29; 95% CI = 0.82–6.45; P = 0.11217). Overall, rare and predicted damaging BARD1 missense variants were significantly more prevalent in BC index patients compared with control individuals (OR = 2.15; 95% CI = 1.26–3.67; P = 0.00723). Neither LoF variants nor predicted damaging rare missense variants in BARD1 were identified in 451 familial index patients with OC.ConclusionsDue to the significant association of germline LoF variants in BARD1 with early-onset BC, we suggest that intensified BC surveillance programs should be offered to women carrying pathogenic BARD1 gene variants.

Mutation in NPPA causes atrial fibrillation by activating inflammation and cardiac fibrosis in a knock‐in rat model

Atrial fibrillation (AF) affects >30 million individuals worldwide. However, no genetic mutation from human patients with AF has been linked to inflammation. Here, we show that AF-associated human variant p.Ile138Thr in natriuretic peptide A (NPPA) encoding the atrial natriuretic peptide (ANP) causes inflammation, fibroblast activation, atrial fibrosis, and AF in knock-in (KI) rats. Variant p.Ile138Thr inhibits the interaction between ANP and its receptor natriuretic peptide receptor A and reduces intracellular cGMP levels. RNA sequencing and follow-up analyses showed that mutant ANP (mANP) activates multiple innate immunity pathways, including TNF-α, NF-κB, and IL-1β signaling. mANP induces differentiation of cardiac fibroblasts (CFs) to myofibroblasts and promotes CF proliferation and fibrosis. These results suggest that NPPA variant p.Ile138Thr causes AF by activating TNF-α, NF-κB, and IL-1β signaling, inflammation, and fibrosis. Multiple computational programs suggest that p.Ile138Thr is damaging or deleterious. Based on the 2015 American College of Medical Genetics and Genomics Standards and Guidelines, p.Ile138Thr can be classified as a likely pathogenic variant. Variant p.Ile138Thr was found only in Asian people in the Genome Aggregation Database and Exome Aggregation Consortium database at an averaged frequency of 0.026%. An estimated 1.15 million Asian people carry the variant and might be at risk of AF. The KI rats may provide an inflammation-based, genetic animal model for AF valuable for testing anti-inflammation or other therapies for AF.-Cheng, C., Liu, H., Tan, C., Tong, D., Zhao, Y., Liu, X., Si, W., Wang, L., Liang, L., Li, J., Wang, C., Chen, Q., Du, Y., Wang, Q. K., Ren, X. Mutation in NPPA causes atrial fibrillation by activating inflammation and cardiac fibrosis in a knock-in rat model.

Autism risk genes are evolutionarily ancient and maintain a unique feature landscape that echoes their function.

Previous research on autism risk (ASD), developmental regulatory (DevReg), and central nervous system (CNS) genes suggests they tend to be large in size, enriched in nested repeats, and mutation intolerant. The relevance of these genomic features is intriguing yet poorly understood. In this study, we investigated the feature landscape of these gene groups to discover structural themes useful in interpreting their function, developmental patterns, and evolutionary history. ASD, DevReg, CNS, housekeeping, and whole genome control (WGC) groups were compiled using various resources. Multiple gene features of interest were extracted from NCBI/UCSC Bioinformatics. Residual variation intolerance scores, Exome Aggregation Consortium pLI scores, and copy number variation data from Decipher were used to estimate variation intolerance. Gene age and protein-protein interactions (PPI) were estimated using Ensembl and EBI Intact databases, respectively. Compared to WGC: ASD, DevReg, and CNS genes are longer, produce larger proteins, maintain greater numbers/density of conserved noncoding elements and transposable elements, produce more transcript variants, and are comparatively variation intolerant. After controlling for gene size, mutation tolerance, and clinical association, ASD genes still retain many of these same features. In addition, we also found that ASD genes that are extremely mutation intolerant have larger PPI networks. These data support many of the recent findings within the field of autism genetics but also expand our understanding of the evolution of these broad gene groups, their potential regulatory complexity, and the extent to which they interact with the cellular network. Autism Res 2019, 12: 860-869. © 2019 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: Autism risk genes are more ancient compared to other genes in the genome. As such, they exhibit physical features related to their age, including long gene and protein size and regulatory sequences that help to control gene expression. They share many of these same features with other genes that are expressed in the brain and/or are associated with prenatal development.

A 17-Beta-Hydroxysteroid Dehydrogenase 13 Variant Protects From Hepatocellular Carcinoma Development in Alcoholic Liver Disease.

Recently, a loss of function variant (rs72613567) in 17-beta-hydroxysteroid dehydrogenase 13 (HSD17B13) has been identified as protective of nonalcoholic (NAFLD) and alcoholic liver disease (ALD). However, the role of this single-nucleotide polymorphism (SNP) in the development of hepatocellular carcinoma (HCC) is currently unknown. A total of 3,315 European patients with HCC (n=1,109) or without HCC, but with chronic liver disease (CLD; n=2,206), from four centers were analyzed either by whole-exome sequencing (WES; exploratory cohort, 285 HCC) or genotyped for HSD17B13 rs72613567 (validation cohort, 824 HCC and all CLD). We included a control group of 33,337 healthy European individuals from the Exome Aggregation Consortium. We compared distribution of genotype using the chi-square test and logistic regression. In the exploratory cohort analyzed by WES, frequency of the TA allele of HSD17B13 rs72613567 was significantly decreased in HCC patients compared to healthy controls (P=1.52×10-06 ). In the validation cohort, frequency of TA allele carriers was also decreased in patients with CLD and without HCC (39%) compared to healthy individuals (47%; P<0.0001). The protective effect of the TA allele of HSD17B13 rs72613567 was identified in patients with ALD (odds ratio [OR]=0.73; 95% confidence interval [CI], 0.65-0.82; P<0.0001), NAFLD (OR=0.64; 95% CI, 0.49-0.83; P=0.0007), and hepatitis C (OR=0.71; 95% CI, 0.60-0.85; P=0.0002). In patients with ALD, the proportion of TA allele carriers with HCC was significantly lower (32%) than in CLD patients without HCC (40%), even after adjustment for age, sex, and fibrosis (OR=0.64; 95% CI, 0.46-0.87; P=0.005). Conclusion: The HSD17B13 rs72613567 loss of function variant is protective of HCC development in patients with ALD.

MON-335 Phenocopy of Multiple Endocrine Neoplasia Type 1 (MEN1) Due to a Germline Cell Division Cycle 73 (CDC73) Variant

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant disorder characterized by the combined occurrence of parathyroid tumors, and neuroendocrine tumors (NETs) of the pituitary and pancreas. MEN1 is caused by germline mutations of the tumor suppressor gene MEN1, which are found in >75% of MEN1 patients. The remaining patients may have mutations involving: as yet unidentified genes; or other genes with known involvement in endocrine tumors, e.g. cyclin dependent kinase inhibitor 1B (CDKN1B) which is associated with MEN4. Here, we report a patient who had primary hyperparathyroidism, acromegaly, and a pancreatic NET that showed immuno-staining for glucagon and chromogranin A, which are consistent with MEN1. However, mutational analysis did not identify a MEN1 mutation, and analysis of other genes was undertaken, after informed consent was obtained from the patient. These genes include: CDKN1B; CDKN1A; CDKN2B; CDKN2C;cell cycle division 73 (CDC73), causative for hyperparathyroidism-jaw tumour (HPT-JT) syndrome, an autosomal dominant disorder characterised by occurrence of parathyroid tumors, ossifying fibromas of the jaw, renal tumours and uterine tumors; calcium sensing receptor (CASR), causative for familial hypocalciuric hypercalcemia type 1 (FHH1); and aryl-hydrocarbon receptor-interacting protein (AIP), causative for familial pituitary tumors. This revealed a heterozygous c.1138C>T (p.Leu380Phe) missense mutation of CDC73. This is likely a disease-causing mutation as: 1) the p.Leu380Phe substitution is not present in >120,000 alleles of The Exome Aggregation Consortium (ExAc) database; 2) involves an evolutionary conserved Leu380 residue; and 3) is situated within the parafibromin domain predicted to interact with the polymerase associated factor 1 (PAF1) and RNA polymerase II, (POLR2A) complex. In addition, RNA-Scope analysis, to detect specific CDC73 mRNA transcripts in paraffin embedded sections, of the patient’s pancreatic NET revealed that CDC73 mRNA expression was significantly decreased by >13% (p<0.005) in the tumor cells, compared to normal adjacent islets. Thus, our results suggest the spectrum of tumours associated with CDC73 mutations also includes pancreatic NETs, and demonstrate that CDC73 mutations may result in a phenocopy of MEN1.

Rare BANF1 Alleles and Relatively Frequent EMD Alleles Including 'Healthy Lipid' Emerin p.D149H in the ExAC Cohort.

Emerin (EMD) and barrier to autointegration factor 1 (BANF1) each bind A-type lamins (LMNA) as fundamental components of nuclear lamina structure. Mutations in LMNA, EMD and BANF1 are genetically linked to many tissue-specific disorders including Emery-Dreifuss muscular dystrophy and cardiomyopathy (LMNA, EMD), lipodystrophy, insulin resistance and type 2 diabetes (LMNA) and progeria (LMNA, BANF1). To explore human genetic variation in these genes, we analyzed EMD and BANF1 alleles in the Exome Aggregation Consortium (ExAC) cohort of 60,706 unrelated individuals. We identified 13 rare heterozygous BANF1 missense variants (p.T2S, p.H7Y, p.D9N, p.S22R, p.G25E, p.D55N, p.D57Y, p.L63P, p.N70T, p.K72R, p.R75W, p.R75Q, p.G79R), and one homozygous variant (p.D9H). Several variants are known (p.G25E) or predicted (e.g., p.D9H, p.D9N, p.L63P) to perturb BANF1 and warrant further study. Analysis of EMD revealed two previously identified variants associated with adult-onset cardiomyopathy (p.K37del, p.E35K) and one deemed ‘benign’ in an Emery-Dreifuss patient (p.D149H). Interestingly p.D149H was the most frequent emerin variant in ExAC, identified in 58 individuals (overall allele frequency 0.06645%), of whom 55 were East Asian (allele frequency 0.8297%). Furthermore, p.D149H associated with four ‘healthy’ traits: reduced triglycerides (-0.336; p = 0.0368), reduced waist circumference (-0.321; p = 0.0486), reduced cholesterol (-0.572; p = 0.000346) and reduced LDL cholesterol (-0.599; p = 0.000272). These traits are distinct from LMNA-associated metabolic disorders and provide the first insight that emerin influences metabolism. We also identified one novel in-frame deletion (p.F39del) and 62 novel emerin missense variants, many of which were relatively frequent and potentially disruptive including p.N91S and p.S143F (∼0.041% and ∼0.034% of non-Finnish Europeans, respectively), p.G156S (∼0.39% of Africans), p.R204G (∼0.18% of Latinx), p.R207P (∼0.08% of South Asians) and p.R221L (∼0.15% of Latinx). Many novel BANF1 variants are predicted to disrupt dimerization or binding to DNA, histones, emerin or A-type lamins. Many novel emerin variants are predicted to disrupt emerin filament dynamics or binding to BANF1, HDAC3, A-type lamins or other partners. These new human variants provide a foundational resource for future studies to test the molecular mechanisms of BANF1 and emerin function, and to understand the link between emerin variant p.D149H and a ‘healthy’ lipid profile.

High prevalence of cancer‐associated TP53 variants in the gnomAD database: A word of caution concerning the use of variant filtering

The 1,000 genome project, the Exome Aggregation Consortium (ExAC) or the Genome Aggregation database (gnomAD) datasets, were developed to provide large-scale reference data of genetic variations for various populations to filter out common benign variants and identify rare variants of clinical importance based on their frequency in the human population. Using a TP53 repository of 80,000 cancer variants, as well as TP53 variants from multiple cancer genome projects, we have defined a set of certified oncogenic TP53 variants. This specific set has been independently validated by functional and in silico predictive analysis. Here we show that a significant number of these variants are included in gnomAD and ExAC. Most of them correspond to TP53 hotspot variants occurring as somatic and germline events in human cancer. Similarly, disease-associated variants for five other tumor suppressor genes, including BRCA1, BRCA2, APC, PTEN, and MLH1, have also been identified. This study demonstrates that germline TP53 variants in the human population are more frequent than previously thought. Furthermore, population databases such as gnomAD or ExAC must be used with caution and need to be annotated for the presence of oncogenic variants to improve their clinical utility.

A novel homozygous FBXO43 mutation associated with male infertility and teratozoospermia in a consanguineous Chinese family

To identify the genetic causes of male infertility characterized by teratozoospermia. Genetic studies. Medical university. Two infertile brothers with teratozoospermia in a consanguineous Chinese family, another 124 sporadic infertile male patients presenting with teratozoospermia, and 200 male controls with normal fertility. None. Whole exome sequencing and genotype analysis to identify the potential pathogenic mutation, Sanger sequencing to validate the mutation in family members, in silico structural modeling to predict the functional consequences of mutation, and targeted next-generation sequencing to validate the mutation in sporadic cases. A novel homozygous nonsynonymous mutation (C1991T, p.G664D) in FBXO43 (F-box only protein 43) was observed in two brothers from a consanguineous Chinese family. The mutation was segregated with the disease phenotype and was predicted to be a disease causing protein by SIFT, PolyPhen-2, and Mutation Taster. An in silico mutant FBXO43 model predicts that the mutation p.G664D causes shortening of two β-sheets, an additional α-helix, and change in loops, which may result in loss of function of the protein. The homozygous mutation of FBXO43 was absent in the 1000 Genomes Project (1000G) and the Exome Aggregation Consortium (ExAC) databases. Subsequent mutation screening of FBXO43 in a cohort of 124 cases identified four additional cases with heterozygous FBXO43 mutations. No mutations were found in FBXO43 in 200 fertile controls. The mutation in FBXO43 is a causative factor of male infertility and teratozoospermia.

Opening a window on lysosomal acid lipase deficiency: Biochemical, molecular, and epidemiological insights.

Lysosomal acid lipase deficiency (LAL-D) is a multi-organ autosomal recessive disease caused by mutations in LIPA. We reviewed data from 681 samples (white blood cells [WBC] n = 625, fibroblasts = 30, liver = 4, amniocytes = 13, chorionic villus = 9) received for analysis of lysosomal acid lipase (LAL) activity over a 15-year period. LIPA sequencing was performed in 49 patients with reduced (n = 26) or deficient (n = 23) LAL activity. The Exome Aggregation Consortium and Genome Aggregation Database dataset were used for LAL-D prevalence calculations. LAL WBC activity was reduced in 67 patients (10.72%) and deficient in 37 (5.92%). The average of LAL activity ± margin of error (CI 95%) was 19.32 ± 0.86 pmol/min/mg for reduced activity patients and 5.90 ± 1.42 pmol/min/mg for deficient patients. The average age at diagnosis for LAL-D was 23.6 years with several patients older than age 30. The correlation between the age at diagnosis and LAL activity showed a significant moderate direct correlation (Pearson's r = 0.46, P < 0.005). Homozygous or compound heterozygous mutations were identified in 9 out of 23 patients with deficient results (detection rate 39.1%). The average LAL activity in molecularly confirmed patients was 4.02 ± 2.02 pmol/min/mg protein, while in molecularly negative patients was 13.886 ± 1.49 pmol/min/mg (P < 0.0001). Twenty-two different mutations were identified including two novel variants (c.309C>A and c.856G>C). A carrier frequency of approximately 1 in 350 was inferred. LAL activity in WBC is a validated tool for LAL-D diagnosis. Higher residual enzymatic activity might result in a milder phenotype leading to diagnosis delay. A cut-off below 12 pmol/min/mg protein might be useful to discriminate patients with LIPA mutations.

3127 The effect of common genetic variants in the oxytocin receptor gene on oxytocin response.

OBJECTIVES/SPECIFIC AIMS: Previous studies suggest that genetic variants in the oxytocin receptor (OXTR) may alter oxytocin dose requirement for labor induction and may increase risk for preterm labor and neurodevelopmental disorders. However, the mechanisms of actions of these variants remain unknown. The goal of this study was to functionally characterize common missense and noncoding variants in OXTR. First, we aimed to determine the effects of missense variants on two major aspects of receptor function: calcium signaling and β-arrestin recruitment. Second, we used allelic expression imbalance assays in an effort to identify regulatory single nucleotide polymorphisms (SNPs) in noncoding regions of OXTR that alter OXTR mRNA expression. METHODS/STUDY POPULATION: We used the Exome Aggregation Consortium database to identify the 12 most prevalent missense single nucleotide variants in OXTR. To determine the functional effects of these variants, we transfected human embryonic kidney cells (a common model system used to study receptor function) with wild type OXTR, variant OXTR, or empty vector control. We used the calcium-sensitive dye Fluo4 to quantify intracellular calcium flux in response to oxytocin treatment, and used bioluminescence resonance energy transfer assays to measure recruitment of the signaling partner β-arrestin to the receptor. To investigate potential effects of noncoding SNPs on OXTR mRNA expression, we quantified allele-specific expression of OXTR in human uterine tissue obtained from participants at the time of Cesarean section. We used next-generation sequencing (Illumina MiSeq) to count alleles of a reporter SNP in OXTR exon 3. RESULTS/ANTICIPATED RESULTS: Of the 12 most prevalent missense single nucleotide variants, four were predicted to be deleterious by PolyPhen variant annotation software. We anticipate that these variants will alter receptor signaling through calcium or β-arrestin pathways. We further observed that a reporter SNP in OXTR exon 3 exhibits significant allelic expression imbalance in a subset of our myometrial tissue samples, indicating that OXTR expression may be regulated by a functional SNP. Our current work focuses on discovering the functional SNPs in OXTR responsible for the pattern of allelic expression imbalance seen in mRNA. In the future, we will seek to explore the effects of these variants on uterine function by using genome editing of uterine smooth muscle cells. DISCUSSION/SIGNIFICANCE OF IMPACT: Our results suggest that both missense and noncoding variants may affect OXTR expression and function. Future studies may suggest that OXTR sequencing, genotyping, or expression analysis would be useful to identify individuals likely to respond or fail to respond to safe doses of oxytocin for labor induction. Personalizing approaches for labor induction in this way would increase the safety of oxytocin and potentially reduce maternal morbidity and mortality.

How good are pathogenicity predictors in detecting benign variants?

Computational tools are widely used for interpreting variants detected in sequencing projects. The choice of these tools is critical for reliable variant impact interpretation for precision medicine and should be based on systematic performance assessment. The performance of the methods varies widely in different performance assessments, for example due to the contents and sizes of test datasets. To address this issue, we obtained 63,160 common amino acid substitutions (allele frequency ≥1% and <25%) from the Exome Aggregation Consortium (ExAC) database, which contains variants from 60,706 genomes or exomes. We evaluated the specificity, the capability to detect benign variants, for 10 variant interpretation tools. In addition to overall specificity of the tools, we tested their performance for variants in six geographical populations. PON-P2 had the best performance (95.5%) followed by FATHMM (86.4%) and VEST (83.5%). While these tools had excellent performance, the poorest method predicted more than one third of the benign variants to be disease-causing. The results allow choosing reliable methods for benign variant interpretation, for both research and clinical purposes, as well as provide a benchmark for method developers.

Variation in the repulsive guidance molecule family in human populations.

Repulsive guidance molecules, RGMA, RGMB, and RGMC, are related proteins discovered independently through different experimental paradigms. They are encoded by single copy genes in mammalian and other vertebrate genomes, and are ~50% identical in amino acid sequence. The importance of RGM actions in human physiology has not been realized, as most research has focused on non‐human models, although mutations in RGMC are the cause of the severe iron storage disorder, juvenile hemochromatosis. Here I show that repositories of human genomic and population genetic data can be used as starting points for discovery and for developing new testable hypotheses about each of these paralogs in human biology and disease susceptibility. Information was extracted, aggregated, and analyzed from the Ensembl and UCSC Genome Browsers, the Exome Aggregation Consortium, the Genotype‐Tissue Expression project portal, the cBio portal for Cancer Genomics, and the National Cancer Institute Genomic Data Commons data site. Results identify extensive variation in gene expression patterns, substantial alternative RNA splicing, and possible missense alterations and other modifications in the coding regions of each of the three genes, with many putative mutations being detected in individuals with different types of cancers. Moreover, selected amino acid substitutions are highly prevalent in the world population, with minor allele frequencies of up to 37% for RGMA and up to 8% for RGMB. These results indicate that protein sequence variation is common in the human RGM family, and raises the possibility that individual variants will have a significant population impact on human physiology and/or disease predisposition.

Prevalence of Germline Mutations Associated With Cancer Risk in Patients With Intraductal Papillary Mucinous Neoplasms

Many patients with pancreatic adenocarcinoma carry germline mutations associated with increased risk of cancer. It is not clear whether patients with intraductal papillary mucinous neoplasms (IPMNs), which are precursors to some pancreatic cancers, also carry these mutations. We assessed the prevalence of germline mutations associated with cancer risk in patients with histologically confirmed IPMN. We obtained nontumor tissue samples from 315 patients with surgically resected IPMNs from 1997 through 2017, and we sequenced 94 genes with variants associated with cancer risk. Mutations associated with increased risk of cancer were identified and compared with individuals from the Exome Aggregation Consortium. We identified 23 patients with a germline mutation associated with cancer risk (7.3%; 95% confidence interval, 4.9-10.8). Nine patients had a germline mutation associated with pancreatic cancer susceptibility (2.9%; 95% confidence interval, 1.4-5.4). More patients with IPMNs carried germline mutations in ATM (P < .0001), PTCH1 (P < .0001), and SUFU (P < .0001) compared with controls. Patients with IPMNs and germline mutations associated with pancreatic cancer were more like to have concurrent invasive pancreatic carcinoma compared with patients with IPMNs without these mutations (P < .0320). In sequence analyses of 315 patients with surgically resected IPMNs, we found that almost 3% to carry mutations associated with pancreatic cancer risk. More patients with IPMNs and germline mutations associated with pancreatic cancer had concurrent invasive pancreatic carcinoma compared with patients with IPMNs without these mutations. Genetic analysis of patients with IPMNs might identify those at greatest risk for cancer.

Association of Inherited Pathogenic Variants in Checkpoint Kinase 2 (CHEK2) With Susceptibility to Testicular Germ Cell Tumors

ImportanceApproximately 50% of the risk for the development of testicular germ cell tumors (TGCTs) is estimated to be heritable, but no mendelian TGCT predisposition genes have yet been identified. It is hypothesized that inherited pathogenic DNA repair gene (DRG) alterations may drive susceptibility to TGCTs.ObjectiveTo systematically evaluate the enrichment of germline pathogenic variants in the mendelian cancer predisposition DRGs in patients with TGCTs vs healthy controls.Design, Setting, and ParticipantsA case-control enrichment analysis was performed from January 2016 to May 2018 to screen for 48 DRGs in 205 unselected men with TGCT and 27 173 ancestry-matched cancer-free individuals from the Exome Aggregation Consortium cohort in the discovery stage. Significant findings were selectively replicated in independent cohorts of 448 unselected men with TGCTs and 442 population-matched controls, as well as 231 high-risk men with TGCTs and 3090 ancestry-matched controls. Statistical analysis took place from January to May 2018.Main Outcomes and MeasuresGene-level enrichment analysis of germline pathogenic variants in individuals with TGCTs relative to cancer-free controls.ResultsAmong 205 unselected men with TGCTs (mean [SD] age, 33.04 [9.67] years), 22 pathogenic germline DRG variants, one-third of which were in CHEK2 (OMIM 604373), were identified in 20 men (9.8%; 95% CI, 6.1%-14.7%). Unselected men with TGCTs were approximately 4 times more likely to carry germline loss-of-function CHEK2 variants compared with cancer-free individuals from the Exome Aggregation Consortium cohort (odds ratio [OR], 3.87; 95% CI, 1.65-8.86; nominal P = .006; q = 0.018). Similar enrichment was also seen in an independent cohort of 448 unselected Croatian men with TGCTs (mean [SD] age, 31.98 [8.11] years) vs 442 unselected Croatian men without TGCTs (at least 50 years of age at time of sample collection) (OR, >1.4; P = .03) and 231 high-risk men with TGCTs (mean [SD] age, 31.54 [9.24] years) vs 3090 men (all older than 50 years) from the Penn Medicine Biobank (OR, 6.30; 95% CI, 2.34-17.31; P = .001). The low-penetrance CHEK2 variant (p.Ile157Thr) was found to be a Croatian founder TGCT risk variant (OR, 3.93; 95% CI, 1.53-9.95; P = .002). Individuals with the pathogenic CHEK2 loss-of-function variants developed TGCTs 6 years earlier than individuals with CHEK2 wild-type alleles (5.95 years; 95% CI, 1.48-10.42; P = .009).Conclusions and RelevanceThis multicenter case-control analysis of men with or without TGCTs provides evidence for CHEK2 as a novel moderate-penetrance TGCT susceptibility gene, with potential clinical utility. In addition to highlighting DNA-repair deficiency as a potential mechanism driving TGCT susceptibility, this analysis also provides new avenues to explore management strategies and biological investigations for high-risk individuals.

Rare mutations of ADAM17 from TOFs induce hypertrophy in human embryonic stem cell-derived cardiomyocytes via HB-EGF signaling.

Tetralogy of Fallot (TOF) is the most common cyanotic form of congenital heart defects (CHDs). The right ventricular hypertrophy is associated with the survival rate of patients with repaired TOF. However, very little is known concerning its genetic etiology. Based on mouse model studies, a disintergrin and metalloprotease 10/17 (ADAM10 and ADAM17) are the key enzymes for the NOTCH and ErbB pathways, which are critical pathways for heart development. Mutations in these two genes have not been previously reported in human TOF patients. In this study, we sequenced ADAM10 and ADAM17 in a Han Chinese CHD cohort comprised of 80 TOF patients, 286 other CHD patients, and 480 matched healthy controls. Three missense variants of ADAM17 were only identified in 80 TOF patients, two of which (Y42D and L659P) are novel and not found in the Exome Aggregation Consortium (ExAC) database. Point mutation knock-in (KI) and ADAM17 knock-out (KO) human embryonic stem cells (hESCs) were generated by CRISPR/Cas9 and programmed to differentiate into cardiomyocytes (CMs). Y42D or L659P KI cells or complete KO cells all developed hypertrophy with disorganized sarcomeres. RNA-seq results showed that phosphatidylinositide 3-kinases/protein kinase B (PI3K/Akt), which is downstream of epidermal growth factor receptor (EGFR) signaling, was affected in both ADAM17 KO and KI hESC-CMs. In vitro experiments showed that these two mutations are loss-of-function mutations in shedding heparin-binding EGF-like growth factor (HB-EGF) but not NOTCH signaling. Our results revealed that CM hypertrophy in TOF could be the result of mutations in ADAM17 which affects HB-EGF/ErbB signaling.

A comprehensive overview of common polymorphic variants that cause missense mutations in human CYPs and UGTs

Genetic polymorphisms in cytochrome P450 (CYP) or UDP-glucuronosyltransferase (UGT) genes can lead to changes in endocrine regulation or drug metabolism. Based on the recently published allele frequency data of the Exome Aggregation Consortium (ExAC), we extracted all common SNP variants that lead to missense mutations in CYPs and UGTs or in their helping proteins CPR, AdR, Adx, and UGDH. With a total of 737 alleles from 83 genes we provide the most comprehensive overview of missense variation distributions in drug metabolizing enzymes published to date. In all variants the most common allele was always considered to be the wild-type (WT), even if it was not identical to the *1-allele and/or the reference standard sequence of the RefSeq project. Surprisingly, in 15 cases (AdR, CYP2A7, CYP2C19, CYP2D6, CYP4A22, CYP4F11, CYP4F12, CYP4V2, CYP8B1, CYP20A1, UGT1A3, UGT1A7, UGT2A3, UGT2B7, and UGT2B15), the WT protein sequences were found to differ from reference standard sequences in up to four amino acids. We expect that these findings will have an impact on the definition of reference sequence standards for these genes, on the corresponding naming of alleles, and on the definition of reference standard activities.

Identification of putative pathogenic single nucleotide variants (SNVs) in genes associated with heart disease in 290 cases of stillbirth.

The incidence of stillbirth in Sweden has essentially remained constant since the 1980’s, and despite thorough investigation, many cases remain unexplained. It has been suggested that a proportion of stillbirth cases is caused by heart disease, mainly channelopathies. The aim of this study was to analyze DNA from 290 stillbirth cases without chromosomal abnormalities for pathogenic single nucleotide variants (SNVs) in 70 genes associated with cardiac channelopathies and cardiomyopathies. The HaloPlex Target Enrichment System (Agilent Technologies) was utilized to prepare sequencing libraries which were sequenced on the Illumina NextSeq platform. We found that 12.1% of the 290 investigated stillbirth cases had one (n = 31) or two (n = 4) variants with evidence supporting pathogenicity, i.e. loss-of-function variants (nonsense, frameshift, splice site substitutions), evidence from functional studies, or previous identification of the variants in affected individuals. Regarding identified putative pathogenic variants in genes associated with channelopathies, the prevalence was significantly higher in the stillbirth cohort (n = 23, 7.93%) than the corresponding prevalence of the same variants in the non-Finnish European population of the Exome Aggregation Consortium (2.70%, p<0.001) and SweGen, (2.30%, p<0.001). Our results give further support to the hypothesis that cardiac channelopathies might contribute to stillbirth. Screening for pathogenic SNVs in genes associated with heart disease might be a valuable complement for stillbirth cases where today’s conventional investigation does not reveal the underlying cause of fetal demise.

Familial Sinus Node Disease Caused by a Gain of GIRK (G-Protein Activated Inwardly Rectifying K+ Channel) Channel Function.

Inherited forms of sinus node dysfunction (SND) clinically include bradycardia, sinus arrest, and chronotropic incompetence and may serve as disease models to understand sinus node physiology and impulse generation. Recently, a gain-of-function mutation in the G-protein gene GNB2 led to enhanced activation of the GIRK (G-protein activated inwardly rectifying K+ channel). Thus, human cardiac GIRK channels are important for heart rate regulation and subsequently, genes encoding their subunits Kir3.1 and Kir3.4 ( KCNJ3 and KCNJ5) are potential candidates for inherited SND in human. We performed a combined approach of targeted sequencing of KCNJ3 and KCNJ5 in 52 patients with idiopathic SND and subsequent whole exome sequencing of additional family members in a genetically affected patient. A putative novel disease-associated gene variant was functionally analyzed by voltage-clamp experiments using various heterologous cell expression systems (Xenopus oocytes, CHO cells, and rat atrial cardiomyocytes). In a 3-generation family with SND we identified a novel variant in KCNJ5 which leads to an amino acid substitution (p.Trp101Cys) in the first transmembrane domain of the Kir3.4 subunit of the cardiac GIRK channel. The identified variant cosegregated with the disease in the family and was absent in the Exome Variant Server and Exome Aggregation Consortium databases. Expression of mutant Kir3.4 (±native Kir3.1) in different heterologous cell expression systems resulted in increased GIRK currents ( IK,ACh) and a reduced inward rectification which was not compensated by intracellular spermidine. Moreover, in silico modeling of heterotetrameric mutant GIRK channels indicates a structurally altered binding site for spermine. For the first time, an inherited gain-of-function mutation in the human GIRK3.4 causes familial human SND. The increased activity of GIRK channels is likely to lead to a sustained hyperpolarization of pacemaker cells and thereby reduces heart rate. Modulation of human GIRK channels may pave a way for further treatment of cardiac pacemaking.