Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients.

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

Genome-wide data from medieval German Jews show that the Ashkenazi founder event pre-dated the 14th century.

We report genome-wide data from 33 Ashkenazi Jews (AJ), dated to the 14th century, obtained following a salvage excavation at the medieval Jewish cemetery of Erfurt, Germany. The Erfurt individuals are genetically similar to modern AJ, but they show more variability in Eastern European-related ancestry than modern AJ. A third of the Erfurt individuals carried a mitochondrial lineage common in modern AJ and eight carried pathogenic variants known to affect AJ today. These observations, together with high levels of runs of homozygosity, suggest that the Erfurt community had already experienced the major reduction in size that affected modern AJ. The Erfurt bottleneck was more severe, implying substructure in medieval AJ. Overall, our results suggest that the AJ founder event and the acquisition of the main sources of ancestry pre-dated the 14th century and highlight late medieval genetic heterogeneity no longer present in modern AJ.

Cancer Risk C (CR-C), a functional genomics test is a sensitive and rapid test for germline mismatch repair deficiency.

PURPOSE: Heritable pathogenic variants in the DNA mismatch repair (MMR) pathway cause Lynch syndrome, a condition that significantly increases risk of colorectal and other cancers. At least half of individuals tested using gene panel sequencing have a variant of uncertain significance or no variant identified leading to no diagnosis. To fill this diagnostic gap, we developed Cancer Risk C (CR-C), a flow variant assay test. METHODS: In response to treatment with an alkylating agent, individual assays of the nuclear translocation of MLH1, MSH2, BARD1, PMS2, and BRCA2 proteins and the nuclear phosphorylation of the ATM and ATR proteins distinguished pathogenic/likely pathogenic (P/LP) from benign/likely benign variants in MMR genes. RESULTS: A risk classification score based on MLH1, MSH2, and ATR assays was 100% sensitive and 98% specific. Causality of MMR P/LP variants was shown through gene editing and rescue. In individuals with suspected Lynch syndrome but no P/LP, CR-C identified most (73%) as having germline MMR defects. Direct comparison of CR-C on matched blood samples and lymphoblastoid cell lines yielded comparable results (r2 > 0.9). CONCLUSION: For identifying germline MMR defects, CR-C provides augmentation to traditional panel sequencing through greater accuracy, shorter turnaround time (48 hours), and performance on blood with minimal sample handling.

Pathogenic Variants in MAP3K1 Cause 46,XY Gonadal Dysgenesis: A Review.

Pathogenic variants in the MAP3K1 gene are an important cause of 46,XY non-syndromic partial and complete gonadal dysgenesis, accounting for at least 4% of cases. Inheritance occurs in a sex-limited, autosomal dominant fashion with virtually complete penetrance in 46,XY individuals. 46,XX carriers appear to have normal fertility and no developmental abnormalities. Pathogenic variants occur almost exclusively within known domains of the MAP3K1 protein, facilitating annotation when identified. Where studied, these variants have been modeled to alter the local MAP3K1 folding and surface domains and have been shown to alter interactions with known binding partners. The net effect of these variants is to increase phosphorylation of downstream targets ERK1, ERK2, and p38, resulting in multiple gain-of-function effects interfering with testis determination and enabling ovarian determination.

Prediction of breast cancer risk based on flow variant analysis of circulating peripheral blood mononuclear cells.

Identifying women at high risk for developing breast cancer is potentially lifesaving. Patients with pathogenic genetic variants can embark on a program of surveillance for early detection, chemoprevention, and/or prophylactic surgery. Newly diagnosed cancer patients can also use the results of gene panel sequencing to make decisions about surgery; therefore, rapid turnaround time for results is critical. Cancer Risk B (CR-B), a test that uses flow variant assays to assess the effects of variants in the DNA double-strand break repair, was applied to two groups of subjects who underwent coincidental gene panel testing, thereby allowing an assessment of sensitivity, specificity and accuracy, and utility for annotating variants of uncertain significance (VUS). The test was compared in matched peripheral blood mononuclear cells (PBMCs) and lymphoblastoid cells (LCLs) and tested for rescue in LCLs with gene transfer. The CR-B phenotype demonstrated a bimodal distribution: CR-B+ indicative of DSB repair defects, and CR-B-, indicative of wild-type repair. When comparing matched LCLs and PBMCs and inter-day tests, CR-B yielded highly reproducible results. The CR-B- phenotype was rescued by gene transfer using wild-type cDNA expression plasmids. The CR-B- phenotype predicted VUS as benign or likely benign. CR-B could represent a rapid alternative to panel sequencing for women with cancer and identifying women at high risk for cancer and is a useful adjunct for annotating VUS.

Novel ultra-rare exonic variants identified in a founder population implicate cadherins in schizophrenia.

The identification of rare variants associated with schizophrenia has proven challenging due to genetic heterogeneity, which is reduced in founder populations. In samples from the Ashkenazi Jewish population, we report that schizophrenia cases had a greater frequency of novel missense or loss of function (MisLoF) ultra-rare variants (URVs) compared to controls, and the MisLoF URV burden was inversely correlated with polygenic risk scores in cases. Characterizing 141 "case-only" genes (MisLoF URVs in ≥3 cases with none in controls), the cadherin gene set was associated with schizophrenia. We report a recurrent case mutation in PCDHA3 that results in the formation of cytoplasmic aggregates and failure to engage in homophilic interactions on the plasma membrane in cultured cells. Modeling purifying selection, we demonstrate that deleterious URVs are greatly overrepresented in the Ashkenazi population, yielding enhanced power for association studies. Identification of the cadherin/protocadherin family as risk genes helps specify the synaptic abnormalities central to schizophrenia.

Survey of Radiation Oncologists to Assess Interest and Potential Use of a Genetic Test Predicting Susceptibility for the Development of Toxicities After Prostate Cancer Radiation Therapy.

PURPOSE: A genetic test predicting susceptibility for the development of toxicities after prostate cancer radiation therapy is in development. This test intends to help physicians with treatment decision making. METHODS AND MATERIALS: Radiation oncologists were surveyed using a web-based questionnaire to gauge their interest in using a genetic test predictive of increased risk of radiation therapy toxicities as an aid in determining therapy for men with prostate cancer. Responses were summarized using frequencies, and a χ2 test compared responses among participants. Multivariable ordinal regression identified factors associated with anticipated adoption or nonadoption of such a genetic test by radiation oncologists. RESULTS: Among 204 radiation oncologists (64% from the United States, 36% from other countries), 86.3% would order a genetic test and 80.2% said the test would be useful for treatment discussions. There was wide acceptance (76.7%) to offer a genetic test to all patients considering radiation therapy for prostate cancer. Additionally, 98.1% indicated that patients would be receptive to the test information. There were no significant differences in the likelihood of ordering a genetic test based on practice setting, familiarity with scientific literature, time spent on research, or geographic location (all P > .05). CONCLUSIONS: Radiation oncologists who treat prostate cancer are interested in and willing to order a genetic test predictive of susceptibility to radiation therapy toxicity to aid their treatment decision making.

Radiogenomics Consortium Genome-Wide Association Study Meta-Analysis of Late Toxicity After Prostate Cancer Radiotherapy.

BACKGROUND: A total of 10%-20% of patients develop long-term toxicity following radiotherapy for prostate cancer. Identification of common genetic variants associated with susceptibility to radiotoxicity might improve risk prediction and inform functional mechanistic studies. METHODS: We conducted an individual patient data meta-analysis of six genome-wide association studies (n = 3871) in men of European ancestry who underwent radiotherapy for prostate cancer. Radiotoxicities (increased urinary frequency, decreased urinary stream, hematuria, rectal bleeding) were graded prospectively. We used grouped relative risk models to test associations with approximately 6 million genotyped or imputed variants (time to first grade 2 or higher toxicity event). Variants with two-sided Pmeta less than 5 × 10-8 were considered statistically significant. Bayesian false discovery probability provided an additional measure of confidence. Statistically significant variants were evaluated in three Japanese cohorts (n = 962). All statistical tests were two-sided. RESULTS: Meta-analysis of the European ancestry cohorts identified three genomic signals: single nucleotide polymorphism rs17055178 with rectal bleeding (Pmeta = 6.2 × 10-10), rs10969913 with decreased urinary stream (Pmeta = 2.9 × 10-10), and rs11122573 with hematuria (Pmeta = 1.8 × 10-8). Fine-scale mapping of these three regions was used to identify another independent signal (rs147121532) associated with hematuria (Pconditional = 4.7 × 10-6). Credible causal variants at these four signals lie in gene-regulatory regions, some modulating expression of nearby genes. Previously identified variants showed consistent associations (rs17599026 with increased urinary frequency, rs7720298 with decreased urinary stream, rs1801516 with overall toxicity) in new cohorts. rs10969913 and rs17599026 had similar effects in the photon-treated Japanese cohorts. CONCLUSIONS: This study increases the understanding of the architecture of common genetic variants affecting radiotoxicity, points to novel radio-pathogenic mechanisms, and develops risk models for testing in clinical studies. Further multinational radiogenomics studies in larger cohorts are worthwhile.

Ectopic Otoconin 90 expression in triple negative breast cancer cell lines is associated with metastasis functions.

Triple negative breast cancer (TNBC) is an aggressive tumor with propensity to metastasize and poor treatment options. Improving treatment options would be impactful; thus, finding a tumor-specific cell surface protein with metastasis promoting functions that could be knocked out was the goal of this study. The Otoconin 90 gene (OC90), frequently amplified in tumors on chromosome 8q24.22, was identified as a potential therapeutic candidate. Normally OC90 is expressed in the cochlea with no known function in other normal tissues. In silico analysis of The Cancer Genome Atlas (TCGA) multi-tumor RNAseq cohorts revealed that OC90 is expressed in many tumor types at high prevalence and genomic amplification is associated with the elevated mRNA expression. In vitro assays in TNBC cell lines revealed OC90 expression with control over cell viability, apoptosis and invasion. RNA-seq analysis of OC90-siRNA knockdown and OC90-overexpression in BT20, BT549, HCC38 cell lines identified co-expressed transcripts, HMGA2, POLE2 and TRIB3. Altered expression of HMGA2, POLE2 and TRIB3 was predictive of survival among members of the Metabric breast cancer cohort. Thus, OC90 represents a potential therapeutic target whose knockdown could improve the treatment of TNBC.

Mutations in MAP3K1 that cause 46,XY disorders of sex development disrupt distinct structural domains in the protein.

Missense mutations in the gene, MAP3K1, are a common cause of 46,XY gonadal dysgenesis, accounting for 15-20% of cases [Ostrer, 2014, Disorders of sex development (DSDs): an update. J. Clin. Endocrinol. Metab., 99, 1503-1509]. Functional studies demonstrated that all of these mutations cause a protein gain-of-function that alters co-factor binding and increases phosphorylation of the downstream MAP kinase pathway targets, MAPK11, MAP3K and MAPK1. This dysregulation of the MAP kinase pathway results in increased CTNNB1, increased expression of WNT4 and FOXL2 and decreased expression of SRY and SOX9. Unique and recurrent pathogenic mutations cluster in three semi-contiguous domains outside the kinase region of the protein, a newly identified N-terminal domain that shares homology with the Guanine Exchange Factor (residues Met164 to Glu231), a Plant HomeoDomain (residues Met442 to Trp495) and an ARMadillo repeat domain (residues Met566 to Glu862). Despite the presence of the mutation clusters and clinical data, there exists a dearth of mechanistic insights behind the development imbalance. In this paper, we use structural modeling and functional data of these mutations to understand alterations of the MAP3K1 protein and the effects on protein folding, binding and downstream target phosphorylation. We show that these mutations have differential effects on protein binding depending on the domains in which they occur. These mutations increase the binding of the RHOA, MAP3K4 and FRAT1 proteins and generally decrease the binding of RAC1. Thus, pathologies in MAP3K1 disrupt the balance between the pro-kinase activities of the RHOA and MAP3K4 binding partners and the inhibitory activity of RAC1.

Rapid Next-Generation Sequencing Method for Prediction of Prostate Cancer Risks.

Prostate cancer is the most commonly diagnosed male cancer and the second leading cause of cancer deaths among men in the United States, with approximately 220,000 new diagnoses and approximately 27,000 deaths each year. Men with clinical low-risk disease can receive active surveillance to safely preserve quality of life, provided that the risk of an undetected aggressive cancer can be managed. Thus, prediction of a tumor's metastatic potential, ideally using only a biopsy sample, is critical to choosing appropriate treatment. We previously proposed and verified a metastasis potential score (MPS) based on regions prone to copy number alterations in metastatic prostate cancer; MPS is highly predictive of metastatic potential in primary tumors. We developed a novel, targeted postligation amplification sequencing approach, which we call the next-generation copy number alteration assay, to efficiently interrogate 902 genomic sites that belong to 194 genomic regions used in the MPS calculation. The assay is designed to work with the latest generation of sequencing platforms to produce estimates of copy number alteration events. The assay's technical reproducibility, robustness to low starting genomic material, and accuracy have been verified. The assay performed very well on cell lines, a cohort of prostate cancer surgical research samples, and matched punched biopsy samples, making it a significant step toward incorporating sequencing techniques for prostate cancer evaluation.

Radiation biology and oncology in the genomic era.

Radiobiology research is building the foundation for applying genomics in precision radiation oncology. Advances in high-throughput approaches will underpin increased understanding of radiosensitivity and the development of future predictive assays for clinical application. There is an established contribution of genetics as a risk factor for radiotherapy side effects. An individual's radiosensitivity is an inherited polygenic trait with an architecture that includes rare mutations in a few genes that confer large effects and common variants in many genes with small effects. Current thinking is that some will be tissue specific, and future tests will be tailored to the normal tissues at risk. The relationship between normal and tumor cell radiosensitivity is poorly understood. Data are emerging suggesting interplay between germline genetic variation and epigenetic modification with growing evidence that changes in DNA methylation regulate the radiosensitivity of cancer cells and histone acetyltransferase inhibitors have radiosensitizing effects. Changes in histone methylation can also impair DNA damage response signaling and alter radiosensitivity. An important effort to advance radiobiology in the genomic era was establishment of the Radiogenomics Consortium to enable the creation of the large radiotherapy cohorts required to exploit advances in genomics. To address challenges in harmonizing data from multiple cohorts, the consortium established the REQUITE project to collect standardized data and genotyping for ~5,000 patients. The collection of detailed dosimetric data is important to produce validated multivariable models. Continued efforts will identify new genes that impact on radiosensitivity to generate new knowledge on toxicity pathogenesis and tests to incorporate into the clinical decision-making process.

Robust genomic copy number predictor of pan cancer metastasis.

Copy number alterations(CNAs) are the most common genetic changes observed in many cancers, reflecting the innate chromosomal instability of this disorder. Yet, how these alterations affect gene function to promote metastases across different tumor types has not been established. In this study, we developed a pan-cancer metastasis potential score (panMPS) based on observed CNAs. panMPS predicts metastasis and metastasis-free survival in cohorts of patients with prostate cancer, triple negative breast cancer and lung adenocarcinoma, and overall survival in the Metabric breast cancer cohort and three cohorts from The Cancer Genome Atlas (TCGA), including prostate, breast and lung adenocarcinoma. These CNAs are present in cell lines of metastatic tumors from eight different origins, reflected by an elevated panMPS for all cell lines. Many copy number alterations involve large chromosomal segments that encompass multiple genes ("clumps"). We show that harnessing this structural information to select only one gene per clump captures the contributions of other genes within the clump, resulting in a robust predictor of metastasis outcome. These sets of selected genes are distinct from cancer drivers that undergo mutation, and in fact, metastasis-related functions have been published for over half of them.

High-depth whole genome sequencing of an Ashkenazi Jewish reference panel: enhancing sensitivity, accuracy, and imputation.

While increasingly large reference panels for genome-wide imputation have been recently made available, the degree to which imputation accuracy can be enhanced by population-specific reference panels remains an open question. Here, we sequenced at full-depth (≥ 30×), across two platforms (Illumina X Ten and Complete Genomics, Inc.), a moderately large (n = 738) cohort of samples drawn from the Ashkenazi Jewish population. We developed a series of quality control steps to optimize sensitivity, specificity, and comprehensiveness of variant calls in the reference panel, and then tested the accuracy of imputation against target cohorts drawn from the same population. Quality control (QC) thresholds for the Illumina X Ten platform were identified that permitted highly accurate calling of single nucleotide variants across 94% of the genome. QC procedures also identified numerous regions that are poorly mapped using current reference or alternate assemblies. After stringent QC, the population-specific reference panel produced more accurate and comprehensive imputation results relative to publicly available, large cosmopolitan reference panels, especially in the range of rare variants that may be most critical to further progress in mapping of complex phenotypes. The population-specific reference panel also permitted enhanced filtering of clinically irrelevant variants from personal genomes.

Machine Learning on a Genome-wide Association Study to Predict Late Genitourinary Toxicity After Prostate Radiation Therapy.

PURPOSE: Late genitourinary (GU) toxicity after radiation therapy limits the quality of life of prostate cancer survivors; however, efforts to explain GU toxicity using patient and dose information have remained unsuccessful. We identified patients with a greater congenital GU toxicity risk by identifying and integrating patterns in genome-wide single nucleotide polymorphisms (SNPs). METHODS AND MATERIALS: We applied a preconditioned random forest regression method for predicting risk from the genome-wide data to combine the effects of multiple SNPs and overcome the statistical power limitations of single-SNP analysis. We studied a cohort of 324 prostate cancer patients who were self-assessed for 4 urinary symptoms at 2 years after radiation therapy using the International Prostate Symptom Score. RESULTS: The predictive accuracy of the method varied across the symptoms. Only for the weak stream endpoint did it achieve a significant area under the curve of 0.70 (95% confidence interval 0.54-0.86; P = .01) on hold-out validation data that outperformed competing methods. Gene ontology analysis highlighted key biological processes, such as neurogenesis and ion transport, from the genes known to be important for urinary tract functions. CONCLUSIONS: We applied machine learning methods and bioinformatics tools to genome-wide data to predict and explain GU toxicity. Our approach enabled the design of a more powerful predictive model and the determination of plausible biomarkers and biological processes associated with GU toxicity.

Functional variants in the LRRK2 gene confer shared effects on risk for Crohn's disease and Parkinson's disease.

Crohn's disease (CD), a form of inflammatory bowel disease, has a higher prevalence in Ashkenazi Jewish than in non-Jewish European populations. To define the role of nonsynonymous mutations, we performed exome sequencing of Ashkenazi Jewish patients with CD, followed by array-based genotyping and association analysis in 2066 CD cases and 3633 healthy controls. We detected association signals in the LRRK2 gene that conferred risk for CD (N2081D variant, P = 9.5 × 10-10) or protection from CD (N551K variant, tagging R1398H-associated haplotype, P = 3.3 × 10-8). These variants affected CD age of onset, disease location, LRRK2 activity, and autophagy. Bayesian network analysis of CD patient intestinal tissue further implicated LRRK2 in CD pathogenesis. Analysis of the extended LRRK2 locus in 24,570 CD cases, patients with Parkinson's disease (PD), and healthy controls revealed extensive pleiotropy, with shared genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The LRRK2 N2081D CD risk allele is located in the same kinase domain as G2019S, a mutation that is the major genetic cause of familial and sporadic PD. Like the G2019S mutation, the N2081D variant was associated with increased kinase activity, whereas neither N551K nor R1398H variants on the protective haplotype altered kinase activity. We also confirmed that R1398H, but not N551K, increased guanosine triphosphate binding and hydrolyzing enzyme (GTPase) activity, thereby deactivating LRRK2. The presence of shared LRRK2 alleles in CD and PD provides refined insight into disease mechanisms and may have major implications for the treatment of these two seemingly unrelated diseases.

Assessing risk for Mendelian disorders in a Bronx population.

BACKGROUND: To identify variants likely responsible for Mendelian disorders among the three major ethnic groups in the Bronx that might be useful to include in genetic screening panels or whole exome sequencing filters and to estimate their likely prevalence in these populations. METHODS: Variants from a high-density oligonucleotide screen of 192 members from each of the three ethnic-national populations (African Americans, Puerto Ricans, and Dominicans) were evaluated for overlap with next generation sequencing data. Variants were curated manually for clinical validity and utility using the American College of Medical Genetics (ACMG) scoring system. Additional variants were identified through literature review. RESULTS: A panel of 75 variants displaying autosomal dominant, autosomal recessive, autosomal recessive/digenic recessive, X-linked recessive, and X-linked dominant inheritance patterns representing 39 Mendelian disorders were identified among these populations. CONCLUSION: Screening for a broader range of disorders could offer the benefits of early or presymptomatic diagnosis and reproductive choice.

MAP3K1-related gonadal dysgenesis: Six new cases and review of the literature.

Investigation of disorders of sex development (DSD) has resulted in the discovery of multiple sex-determining genes. MAP3K1 encodes a signal transduction regulator in the sex determination pathway and is emerging as one of the more common genes responsible for 46,XY DSD presenting as complete or partial gonadal dysgenesis. Clinical assessment, endocrine evaluation, and genetic analysis were performed in six individuals from four unrelated families with 46,XY DSD. All six individuals were found to have likely pathogenic MAP3K1 variants. Three of these individuals presented with complete gonadal dysgenesis, characterized by bilateral streak gonads with typical internal and external female genitalia, while the other three presented with partial gonadal dysgenesis, characterized by incomplete testicular development, resulting in clitoral hypertrophy with otherwise typical female external genitalia. Testing for MAP3K1 variants should be considered in patients with 46,XY complete or partial gonadal dysgenesis, particularly in families with multiple members affected with 46,XY DSD. Identification of a MAP3K1 variant should prompt an evaluation for DSD in female siblings of the proband.