Search

Identification of therapeutic targets from genome-wide association studies (GWAS) requires insights into downstream functional consequences. We harmonized 8,613 RNA-sequencing samples from 14 brain datasets to create the MetaBrain resource and performed cis- and trans-expression quantitative trait locus (eQTL) meta-analyses in multiple brain region- and ancestry-specific datasets (n ≤ 2,759). Many of the 16,169 cortex cis-eQTLs were tissue-dependent when compared with blood cis-eQTLs. We inferred brain cell types for 3,549 cis-eQTLs by interaction analysis. We prioritized 186 cis-eQTLs for 31 brain-related traits using Mendelian randomization and co-localization including 40 cis-eQTLs with an inferred cell type, such as a neuron-specific cis-eQTL (CYP24A1) for multiple sclerosis. We further describe 737 trans-eQTLs for 526 unique variants and 108 unique genes. We used brain-specific gene-co-regulation networks to link GWAS loci and prioritize additional genes for five central nervous system diseases. This study represents a valuable resource for post-GWAS research on central nervous system diseases.

Genetic colocalisation is an important tool to test for shared genetic aetiology and is commonly used to strengthen causal inference in genetic studies of molecular traits and drug targets. However, the single causal variant assumption of the original colocalization method is a considerable limitation in genomic regions with multiple causal effects.We integrated conditional analyses (GCTA-COJO) and colocalisation analyses (coloc), into a novel analysis tool called Pair-Wise Conditional Colocalization (PWCoCo). PWCoCo performs conditional analyses to identify independent signals for the two tested traits in a genomic region and then conducts colocalisation of each pair of conditionally independent signals for the two traits using summary-level data. This allows for the stringent single-variant assumption to hold for each pair of colocalisation analysis.We found that the computational efficiency of PWCoCo is on average better than colocalisation with Sum of Single Effects Regression using Summary Stats (SuSiE-RSS), with greater gains in efficiency for high-throughput analysis. In a case study using GWAS data for multiple sclerosis and brain cortex-derived eQTLs (MetaBrain), we recapitulated all previously identified genes, which showcased the robustness of the method. We further found colocalisation evidence for secondary signals in nine additional loci, which was not identifiable in conventional GWAS and/or colocalisation.PWCoCo offers key improvements over existing methods, including: (1) robust colocalisation when the single variant assumption is violated; (2) independent colocalisation of secondary signals, which enables identification of novel disease-causing variants; (3) an easy-to-use and computationally efficient tool to test for colocalisation of high-dimensional omics data.

Funding: This work was supported by grants from the Emergency Medicine Foundation Career Development Grant; AHA Mentored Clinical & Population Research Award (14CRP18860027); NIH/NINDS-R01-NS085419 (C.C., J.M.L.); NIH/NINDS-R37-NS107230, NIH/NINDS U24-NS107230 (J.M.L.); NIH/NINDS-K23-NS099487 (L.H.); NIH/NIA-K99-AG062723 (L.I.); Barnes-Jewish Hospital Foundation (J.M.L.); Biogen (C.C., J.M.L.); Bright Focus Foundation, US Department of Defense, Helsinki University Central Hospital; Finnish Medical Foundation; Finland government subsidiary funds; Spanish Ministry of Science and Innovation; Instituto de Salud Carlos III (grants ‘Registro BASICMAR’ Funding for Research in Health (PI051737), ‘GWALA project’ from Fondos de Investigación Sanitaria ISC III (PI10/02064, PI12/01238 and PI15/00451), JR18/00004); Fondos FEDER/EDRF Red de Investigación Cardiovascular (RD12/0042/0020); Fundació la Marató TV3; Genestroke Consortium (76/C/2011); Recercaixa’13 (JJ086116). Tomás Sobrino (CPII17/00027), Francisco Campos (CPII19/00020) and Israel Fernandez are supported by Miguel Servet II Program from Instituto de Salud Carlos III and Fondos FEDER. I.F. is also supported by Maestro project (PI18/01338) and Pre-test project (PMP15/00022) from Instituto de Salud Carlos III and Fondos Feder, Agaur; and Epigenesis project from Marató TV3 Foundation. J.C., J.M., A.D., J.M.-F., J.A. and I.F. are supported by Invictus plus Network (RD16/0019) from Instituto de Salud Carlos III and Fondos Feder. Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP-2013/07559-3) (I.L.-C.), Sigrid Juselius Foundation. The MEGASTROKE project received funding from sources specified at http://www.megastroke.org/acknowledgments.html. B.S., B.A. and F.S. are supported by NIH awards NS097000, NS101718, NS075035, NS079153 and NS106950. During the first hours after stroke onset, neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between the National Institutes of Health Stroke Scale (NIHSS) within 6h of stroke onset and NIHSS at 24h. A total of 5876 individuals from seven countries (Spain, Finland, Poland, USA, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of NIHSS at 24h of variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture from that of stroke risk. Eight loci (1p21.1, 1q42.2, 2p25.1, 2q31.2, 2q33.3, 5q33.2, 7p21.2 and 13q31.1) were genome-wide significant and explained 1.8% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each locus. Expression quantitative trait loci mapping and summary data-based Mendelian randomization indicate that ADAM23 (log Bayes factor = 5.41) was driving the association for 2q33.3. Gene-based analyses suggested that GRIA1 (log Bayes factor = 5.19), which is predominantly expressed in the brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated GNPAT (log Bayes factor = 7.64) ABCB5 (log Bayes factor = 5.97) for the 1p21.1 and 7p21.1 loci. Human brain single-nuclei RNA-sequencing indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a presynaptic protein and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provide the first genetic evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischaemic stroke. Postprint

Biliary atresia (BA) is a severe pediatric liver disease resulting in necroinflammatory obliteration of the extrahepatic biliary tree. BA presents within the first few months of life as either an isolated finding or with additional syndromic features. The etiology of isolated BA is unknown, with evidence for infectious, environmental, and genetic risk factors described. However, to date, there are no definitive causal genes identified for isolated BA in humans, and the question of whether single gene defects play a major role remains open. We performed exome-sequencing in 101 North American patients of European descent with isolated BA (including 30 parent-child trios) and considered several experimental designs to identify potentially deleterious protein-altering variants that may be involved in the disease. In a case-only analysis, we did not identify genes with variants shared among more than two probands, and burden tests of rare variants using a case-case control design did not yield significant results. In the trio analysis of 30 simplex families (patient and parent trios), we identified 66 de novo variants in 66 genes including potentially deleterious variants in STIP1 and REV1. STIP1 is a co-chaperone for the heat-shock protein, HSP90, and has been shown to have diverse functions in yeast, flies and mammals, including stress-responses. REV1 is known to be a key player in DNA repair pathway and to interact with HSP90. In conclusion, our results do not support the hypothesis that a simple genetic model is responsible for the majority of cases of isolated BA. Our finding of de novo variants in genes linked to evolutionarily conserved stress responses (STIP1 and REV1) suggests that exploration of how genetic susceptibility and environmental exposure may interact to cause BA is warranted.

Complex traits are characterized by multiple genes and variants acting simultaneously on a phenotype. However, studying the contribution of individual pairs of genes to complex traits has been challenging since human genetics necessitates very large population sizes, while findings from model systems do not always translate to humans. Here, we combine genetics with combinatorial RNAi (coRNAi) to systematically test for pairwise additive effects (AEs) and genetic interactions (GIs) between 30 lipid genome-wide association studies (GWAS) genes. Gene-based burden tests from 240,970 exomes show that in carriers with truncating mutations in both, APOB and either PCSK9 or LPL (“human double knock-outs”) plasma lipid levels change additively. Genetics and coRNAi identify overlapping AEs for 12 additional gene pairs. Overlapping GIs are observed for TOMM40/APOE with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate plasma and cellular lipid levels primarily via AEs and nominates putative drug target pairs for improved lipid-lowering combination therapies., Studying the contribution of pairs of genes to complex traits has been challenging. Here, the authors combine exome and genotype data with RNAi to screen for genetic interactions between 30 genes identified in lipid GWAS to hint at pairs whose joint modulation may improve lipid-lowering therapies.

Nearly two hundred common-variant depression risk loci have been identified by genome-wide association studies (GWAS)1–4. However, the impact of rare coding variants on depression remains poorly understood. Here, we present the largest to date exome analysis of depression based on 320,356 UK Biobank participants. We show that the burden of rare disruptive coding variants in loss-of-function intolerant genes is significantly associated with depression risk. Among 30 genes with false discovery rate (FDR) SLC2A1, a blood-brain barrier glucose transporter underlying GLUT1 deficiency syndrome5–7, reached exome-wide significance (P=2.96e-7). Gene-set enrichment supports neuron projection development and muscle activities2, 3as implicated in depression. Integrating exomes with polygenic risk revealed additive contributions from common and rare variants to depression risk. The burden of rare disruptive coding variants for depression overlapped with that of developmental disorder, autism and schizophrenia. Our study provides novel insight into the contribution of rare coding variants on depression and genetic relationships across developmental and psychiatric disorders.

Gaining insight into the downstream consequences of non-coding variants is an essential step towards the identification of therapeutic targets from genome-wide association study (GWAS) findings. Here we have harmonized and integrated 8,727 RNA-seq samples with accompanying genotype data from multiple brain-regions from 14 datasets. This sample size enabled us to perform both cis- and trans-expression quantitative locus (eQTL) mapping. Upon comparing the brain cortex cis-eQTLs (for 12,307 unique genes at FDRcis-eQTL analysis (n=31,684 samples), we observed that brain eQTLs are more tissue specific than previously assumed.We inferred the brain cell type for 1,515 cis-eQTLs by using cell type proportion information. We conducted Mendelian Randomization on 31 brain-related traits using cis-eQTLs as instruments and found 159 significant findings that also passed colocalization. Furthermore, two multiple sclerosis (MS) findings had cell type specific signals, a neuron-specific cis-eQTL for CYP24A1 and a macrophage specific cis-eQTL for CLECL1.To further interpret GWAS hits, we performed trans-eQTL analysis. We identified 2,589 trans-eQTLs (at FDRWe also generated a brain-specific gene-coregulation network that we used to predict which genes have brain-specific functions, and to perform a novel network analysis of Alzheimer’s disease (AD), amyotrophic lateral sclerosis (ALS), multiple sclerosis (MS) and Parkinson’s disease (PD) GWAS data. This resulted in the identification of distinct sets of genes that show significantly enriched co-regulation with genes inside the associated GWAS loci, and which might reflect drivers of these diseases.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease with an estimated heritability of around 50%. DNA methylation patterns can serve as biomarkers of (past) exposures and disease progression, as well as providing a potential mechanism that mediates genetic or environmental risk. Here, we present a blood-based epigenome-wide association study (EWAS) meta-analysis in 10,462 samples (7,344 ALS patients and 3,118 controls), representing the largest case-control study of DNA methylation for any disease to date. We identified a total of 45 differentially methylated positions (DMPs) annotated to 42 genes, which are enriched for pathways and traits related to metabolism, cholesterol biosynthesis, and immunity. We show that DNA-methylation-based proxies for HDL-cholesterol, BMI, white blood cell (WBC) proportions and alcohol intake were independently associated with ALS. Integration of these results with our latest GWAS showed that cholesterol biosynthesis was causally related to ALS. Finally, we found that DNA methylation levels at several DMPs and blood cell proportion estimates derived from DNA methylation data, are associated with survival rate in patients, and could represent indicators of underlying disease processes.

Significance Parkinson’s disease (PD) is the second most prevalent neurodegenerative disease of aging, affecting approximately 10 million patients worldwide with no approved therapies to modify progression of disease. Further understanding of the cellular mechanisms contributing to the development of PD is necessary to discover therapies. Here, we characterize the role of a recently identified GWAS hit for sporadic PD, ITPKB, in the aggregation of α-synuclein, the primary pathological feature of disease. These results identify inhibition of inositol-1,4,5,-triphosphate (IP3)-mediated ER-to-mitochondria calcium release as a potential therapeutic approach for reducing neuropathology in PD., Inositol-1,4,5-triphosphate (IP3) kinase B (ITPKB) is a ubiquitously expressed lipid kinase that inactivates IP3, a secondary messenger that stimulates calcium release from the endoplasmic reticulum (ER). Genome-wide association studies have identified common variants in the ITPKB gene locus associated with reduced risk of sporadic Parkinson’s disease (PD). Here, we investigate whether ITPKB activity or expression level impacts PD phenotypes in cellular and animal models. In primary neurons, knockdown or pharmacological inhibition of ITPKB increased levels of phosphorylated, insoluble α-synuclein pathology following treatment with α-synuclein preformed fibrils (PFFs). Conversely, ITPKB overexpression reduced PFF-induced α-synuclein aggregation. We also demonstrate that ITPKB inhibition or knockdown increases intracellular calcium levels in neurons, leading to an accumulation of calcium in mitochondria that increases respiration and inhibits the initiation of autophagy, suggesting that ITPKB regulates α-synuclein pathology by inhibiting ER-to-mitochondria calcium transport. Furthermore, the effects of ITPKB on mitochondrial calcium and respiration were prevented by pretreatment with pharmacological inhibitors of the mitochondrial calcium uniporter complex, which was also sufficient to reduce α-synuclein pathology in PFF-treated neurons. Taken together, these results identify ITPKB as a negative regulator of α-synuclein aggregation and highlight modulation of ER-to-mitochondria calcium flux as a therapeutic strategy for the treatment of sporadic PD.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a private-public partnership between the UK Biobank (UKB) and eight biopharmaceutical companies that will complete the sequencing of exomes for all ~500,000 UKB participants. Here, we describe the early results from ~200,000 UKB participants and the features of this project that enabled its success. The biopharmaceutical industry has increasingly used human genetics to improve success in drug discovery. Recognizing the need for large-scale human genetics data, as well as the unique value of the data access and contribution terms of the UKB, the UKB-ESC was formed. As a result, exome data from 200,643 UKB enrollees are now available. These data include ~10 million exonic variants-a rich resource of rare coding variation that is particularly valuable for drug discovery. The UKB-ESC precompetitive collaboration has further strengthened academic and industry ties and has provided teams with an opportunity to interact with and learn from the wider research community.

The UK Biobank Exome Sequencing Consortium (UKB-ESC) is a unique private/public partnership between the UK Biobank and eight biopharma companies that will sequence the exomes of all ∼500,000 UK Biobank participants. Here we describe early results from the exome sequence data generated by this consortium for the first ∼200,000 UKB subjects and the key features of this project that enabled the UKB-ESC to come together and generate this data.Exome sequencing data from the first 200,643 UKB enrollees are now accessible to the research community. Approximately 10M variants were observed within the targeted regions, including: 8,086,176 SNPs, 370,958 indels and 1,596,984 multi-allelic variants. Of the ∼8M variants observed, 84.5% are coding variants and include 2,139,318 (25.3%) synonymous, 4,549,694 (53.8%) missense, 453,733 (5.4%) predicted loss-of-function (LOF) variants (initiation codon loss, premature stop codons, stop codon loss, splicing and frameshift variants) affecting at least one coding transcript. This open access data provides a rich resource of coding variants for rare variant genetic studies, and is particularly valuable for drug discovery efforts that utilize rare, functionally consequential variants.Over the past decade, the biopharma industry has increasingly leveraged human genetics as part of their drug discovery and development strategies. This shift was motivated by technical advances that enabled cost-effective human genetics research at scale, the emergence of electronic health records and biobanks, and a maturing understanding of how human genetics can increase the probability of successful drug development. Recognizing the need for large-scale human genetics data to drive drug discovery, and the unique value of the open data access policies and contribution terms of the UK Biobank, the UKB-ESC was formed. This precompetitive collaboration has further strengthened the ties between academia and industry and provided teams an unprecedented opportunity to interact with and learn from the wider research community.

During the first hours after stroke onset neurological deficits can be highly unstable: some patients rapidly improve, while others deteriorate. This early neurological instability has a major impact on long-term outcome. Here, we aimed to determine the genetic architecture of early neurological instability measured by the difference between NIH stroke scale (NIHSS) within six hours of stroke onset and NIHSS at 24h (ΔNIHSS). A total of 5,876 individuals from seven countries (Spain, Finland, Poland, United States, Costa Rica, Mexico and Korea) were studied using a multi-ancestry meta-analyses. We found that 8.7% of ΔNIHSS variance was explained by common genetic variations, and also that early neurological instability has a different genetic architecture than that of stroke risk. Seven loci (2p25.1, 2q31.2, 2q33.3, 4q34.3, 5q33.2, 6q26 and 7p21.1) were genome-wide significant and explained 2.1% of the variability suggesting that additional variants influence early change in neurological deficits. We used functional genomics and bioinformatic annotation to identify the genes driving the association from each loci. eQTL mapping and SMR indicate that ADAM23 (log Bayes Factor (LBF)=6.34) was driving the association for 2q33.3. Gene based analyses suggested that GRIA1 (LBF=5.26), which is predominantly expressed in brain, is the gene driving the association for the 5q33.2 locus. These analyses also nominated PARK2 (LBF=5.30) and ABCB5 (LBF=5.70) for the 6q26 and 7p21.1 loci. Human brain single nuclei RNA-seq indicates that the gene expression of ADAM23 and GRIA1 is enriched in neurons. ADAM23, a pre-synaptic protein, and GRIA1, a protein subunit of the AMPA receptor, are part of a synaptic protein complex that modulates neuronal excitability. These data provides the first evidence in humans that excitotoxicity may contribute to early neurological instability after acute ischemic stroke.RESEARCH INTO CONTEXTEvidence before this studyNo previous genome-wide association studies have investigated the genetic architecture of early outcomes after ischemic stroke.Added Value of this studyThis is the first study that investigated genetic influences on early outcomes after ischemic stroke using a genome-wide approach, revealing seven genome-wide significant loci. A unique aspect of this genetic study is the inclusion of all of the major ethnicities by recruiting from participants throughout the world. Most genetic studies to date have been limited to populations of European ancestry.Implications of all available evidenceThe findings provide the first evidence that genes implicating excitotoxicity contribute to human acute ischemic stroke, and demonstrates proof of principle that GWAS of acute ischemic stroke patients can reveal mechanisms involved in ischemic brain injury.

SUMMARYGenetic interactions (GIs), the joint impact of different genes or variants on a phenotype, are foundational to the genetic architecture of complex traits. However, identifying GIs through human genetics is challenging since it necessitates very large population sizes, while findings from model systems not always translate to humans. Here, we combined exome-sequencing and genotyping in the UK Biobank with combinatorial RNA-interference (coRNAi) screening to systematically test for pairwise GIs between 30 lipid GWAS genes. Gene-based protein-truncating variant (PTV) burden analyses from 240,970 exomes revealed additive GIs for APOB with PCSK9 and LPL, respectively. Both, genetics and coRNAi identified additive GIs for 12 additional gene pairs. Overlapping non-additive GIs were detected only for TOMM40 at the APOE locus with SORT1 and NCAN. Our study identifies distinct gene pairs that modulate both, plasma and cellular lipid levels via additive and non-additive effects and nominates drug target pairs for improved lipid-lowering combination therapies.

Genetic predisposition has been shown to contribute substantially to the age at which we die. Genome-wide association studies (GWASs) have linked more than 20 loci to phenotypes related to human lifespan1. However, little is known about how lifespan is impacted by gene loss of function. Through whole-exome sequencing of 352,338 UK Biobank participants of European ancestry, we assessed the relevance of protein-truncating variant (PTV) gene burden on individual and parental survival. We identified four exome-wide significant (P −7) human lifespan genes, BRCA1, BRCA2, ATM and TET2. Gene and gene-set, PTV-burden, phenome-wide association studies support known roles of these genes in cancer to impact lifespan at the population level. The TET2 PTV burden was associated with a lifespan through somatic mutation events presumably due to clonal hematopoiesis. The overlap between PTV burden and common variant-based lifespan GWASs was modest, underscoring the value of exome sequencing in well-powered biobank cohorts to complement GWASs for identifying genes underlying complex traits.

Leukodystrophies are a heterogeneous group of heritable disorders characterized by abnormal brain white matter signal on magnetic resonance imaging (MRI) and primary involvement of the cellular components of myelin. Previous estimates suggest the incidence of leukodystrophies as a whole to be 1 in 7,000 individuals, however the frequency of specific diagnoses relative to others has not been described. Next generation sequencing approaches offer the opportunity to redefine our understanding of the relative frequency of different leukodystrophies. We assessed the relative frequency of all 30 leukodystrophies (associated with 55 genes) in more than 49,000 exomes. We identified a relatively high frequency of disorders previously thought of as very rare, including Aicardi Goutieres Syndrome, TUBB4A-related leukodystrophy, Peroxisomal biogenesis disorders, POLR3-related Leukodystrophy, Vanishing White Matter, and Pelizaeus-Merzbacher Disease. Despite the relative frequency of these conditions, carrier-screening laboratories regularly test only 20 of the 55 leukodystrophy-related genes, and do not test at all, or test only one or a few, genes for some of the higher frequency disorders. Relative frequency of leukodystrophies previously considered very rare suggests these disorders may benefit from expanded carrier screening.

PURPOSE: As massively parallel sequencing is increasingly being used for clinical decision-making, it has become critical to understand parameters that affect sequencing quality and to establish methods for measuring and reporting clinical sequencing standards. In this report, we propose a definition for reduced coverage regions and have established a set of standards for variant calling in clinical sequencing applications. METHODS: To enable sequencing centers to assess the regions of poor sequencing quality in their own data, we optimized and used a tool (ExCid) to identify reduced coverage loci within genes or regions of particular interest. We used this framework to examine sequencing data from 500 patients generated in ten projects from sequencing centers in the NHGRI/NCI Clinical Sequencing Exploratory Research (CSER) Consortium. RESULTS: This approach identified reduced coverage regions in clinically relevant genes, including known clinically relevant loci that were uniquely missed at individual centers, in multiple centers, and in all centers. CONCLUSIONS: This report provides a process roadmap for clinical sequencing centers looking to perform similar analyses on their data.

Summary Background There are more than 300 known red blood cell (RBC) antigens and 33 platelet antigens that differ between individuals. Sensitisation to antigens is a serious complication that can occur in prenatal medicine and after blood transfusion, particularly for patients who require multiple transfusions. Although pre-transfusion compatibility testing largely relies on serological methods, reagents are not available for many antigens. Methods based on single-nucleotide polymorphism (SNP) arrays have been used, but typing for ABO and Rh—the most important blood groups—cannot be done with SNP typing alone. We aimed to develop a novel method based on whole-genome sequencing to identify RBC and platelet antigens. Methods This whole-genome sequencing study is a subanalysis of data from patients in the whole-genome sequencing arm of the MedSeq Project randomised controlled trial (NCT01736566) with no measured patient outcomes. We created a database of molecular changes in RBC and platelet antigens and developed an automated antigen-typing algorithm based on whole-genome sequencing (bloodTyper). This algorithm was iteratively improved to address cis–trans haplotype ambiguities and homologous gene alignments. Whole-genome sequencing data from 110 MedSeq participants (30 × depth) were used to initially validate bloodTyper through comparison with conventional serology and SNP methods for typing of 38 RBC antigens in 12 blood-group systems and 22 human platelet antigens. bloodTyper was further validated with whole-genome sequencing data from 200 INTERVAL trial participants (15 × depth) with serological comparisons. Findings We iteratively improved bloodTyper by comparing its typing results with conventional serological and SNP typing in three rounds of testing. The initial whole-genome sequencing typing algorithm was 99·5% concordant across the first 20 MedSeq genomes. Addressing discordances led to development of an improved algorithm that was 99·8% concordant for the remaining 90 MedSeq genomes. Additional modifications led to the final algorithm, which was 99·2% concordant across 200 INTERVAL genomes (or 99·9% after adjustment for the lower depth of coverage). Interpretation By enabling more precise antigen-matching of patients with blood donors, antigen typing based on whole-genome sequencing provides a novel approach to improve transfusion outcomes with the potential to transform the practice of transfusion medicine. Funding National Human Genome Research Institute, Doris Duke Charitable Foundation, National Health Service Blood and Transplant, National Institute for Health Research, and Wellcome Trust.

Biliary atresia (BA) is a pediatric cholangiopathy with unknown etiology occurring in isolated and syndromic forms. Laterality defects affecting the cardiovascular and gastrointestinal systems are the most common features present in syndromic BA. Most cases are sporadic, although reports of familial cases have led to the hypothesis of genetic susceptibility in some patients. We identified a child with BA, malrotation, and interrupted inferior vena cava whose father presented with situs inversus, polysplenia, panhypopituitarism, and mildly dysmorphic facial features. Chromosomal microarray analysis demonstrated a 277 kb heterozygous deletion on chromosome 20, which included a single gene, FOXA2, in the proband and her father. This deletion was confirmed to be de novo in the father. The proband and her father share a common diagnosis of heterotaxy, but they also each presented with a variety of other issues. Further genetic screening revealed that the proband carried an additional protein-altering polymorphism (rs1904589; p.His165Arg) in the NODAL gene that is not present in the father, and this variant has been shown to decrease expression of the gene. As FOXA2 can be a regulator of NODAL expression, we propose that haploinsufficiency for FOXA2 combined with a decreased expression of NODAL is the likely cause for syndromic BA in this proband.

This report describes an algorithm developed to predict the pathogenicity of copy number variants (CNVs) in large sample cohorts. CNVs (genomic deletions and duplications) are found in healthy individuals and in individuals with genetic diagnoses, and differentiation of these two classes of CNVs can be challenging and usually requires extensive manual curation. We have developed PECONPI, an algorithm to assess the pathogenicity of CNVs based on gene content and CNV frequency. This software was applied to a large cohort of patients with genetically heterogeneous non-syndromic hearing loss to score and rank each CNV based on its relative pathogenicity. Of 636 individuals tested, we identified the likely underlying etiology of the hearing loss in 14 (2%) of the patients (1 with a homozygous deletion, 7 with a deletion of a known hearing loss gene and a point mutation on the trans allele and 6 with a deletion larger than 1 Mb). We also identified two probands with smaller deletions encompassing genes that may be functionally related to their hearing loss. The ability of PECONPI to determine the pathogenicity of CNVs was tested on a second genetically heterogeneous cohort with congenital heart defects (CHDs). It successfully identified a likely etiology in 6 of 355 individuals (2%). We believe this tool is useful for researchers with large genetically heterogeneous cohorts to help identify known pathogenic causes and novel disease genes.

Hearing impairment affects 1 in 650 newborns, making it the most common congenital sensory impairment. Autosomal recessive nonsyndromic sensorineural hearing impairment (ARNSHI) comprises 80% of familial hearing impairment cases. Mutations in GJB2 account for a significant number of ARNSHI (and up to 50% of documented recessive (e.g., more than 1 affected sibling) hearing impairment in some populations). Mutations in the GJB2 gene are amongst the most common causes of hearing impairment in populations of various ethnic backgrounds. Two mutations of this gene, 35delG and 167delT, account for the majority of reported mutations in Caucasian populations, especially those of Mediterranean and Ashkenazi Jewish background. The 235delC mutation is most prevalent in East Asian populations. Some mutations are of less well-characterized significance. The V37I missense mutation, common in Asian populations, was initially described as a polymorphism and later as a potentially pathogenic mutation. We report here on 15 unrelated individuals with ARNSHI and homozygosity for the V37I GJB2 missense mutation. Nine individuals are of Chinese ancestry, two are of unspecified Asian descent, one is of Japanese descent, one individual is of Vietnamese ancestry, one of Philippine background and one of Italian and Cuban/Caucasian background. Homozygosity for the V37I GJB2 mutation may be a more common pathogenic missense mutation in Asian populations, resulting in mild to moderate sensorineural hearing impairment. We report a presumed haplotype block specific to East Asian individuals with the V37I mutation encompassing the GJB2 gene that may account for the high prevalence in East Asian populations.

Summary: We report our new DRAW+SneakPeek software for DNA-seq analysis. DNA resequencing analysis workflow (DRAW) automates the workflow of processing raw sequence reads including quality control, read alignment and variant calling on high-performance computing facilities such as Amazon elastic compute cloud. SneakPeek provides an effective interface for reviewing dozens of quality metrics reported by DRAW, so users can assess the quality of data and diagnose problems in their sequencing procedures. Both DRAW and SneakPeek are freely available under the MIT license, and are available as Amazon machine images to be used directly on Amazon cloud with minimal installation. Availability: DRAW+SneakPeek is released under the MIT license and is available for academic and nonprofit use for free. The information about source code, Amazon machine images and instructions on how to install and run DRAW+SneakPeek locally and on Amazon elastic compute cloud is available at the National Institute on Aging Genetics of Alzheimer’s Disease Data Storage Site (http://www.niagads.org/) and Wang lab Web site (http://wanglab.pcbi.upenn.edu/). Contact: gerardsc@mail.med.upenn.edu or lswang@mail.med.upenn.edu

We established a bioinformatics procedure at our laboratory that processes sequencing data, identifies single nucleotide variants (SNVs) and small insertions/deletions (indels), annotates the variants based on the curated variant databases and algorithms, and filters variants according to custom defined filtration schemas that are selected by geneticists specifically for each case. Since we are operating in the clinical settings, rigorous validation procedures are applied at each step to ensure the accuracy and quality of the data. Here we describe how we implemented our current clinical exome pipeline, and our approaches for validation.

In the United States, biliary atresia (BA) is the most frequent indication for liver transplantation in pediatric patients. BA is a complex disease, with suspected environmental and genetic risk factors. A genome-wide association study in Chinese patients identified association to the 10q24.2 (hg18) genomic region. This signal was upstream of two genes, XPNPEP1 and ADD3, both expressed in intrahepatic bile ducts. We tested association to this region in 171 BA patients and 1,630 controls of European descent and found the strongest signal to be at rs7099604 (p = 2.5 × 10(-3)) in intron 1 of the ADD3 gene. Moreover, expression data suggest that ADD3, but not XPNPEP1, is differentially expressed in BA patients. The role of ADD3 in biliary development is unclear, but our findings suggest that this gene may be functionally relevant for the development of BA.

Background & Aims Biliary atresia (BA) is a progressive fibroinflammatory disorder of infants involving the extrahepatic and intrahepatic biliary tree. Its etiology is unclear but is believed to involve exposure of a genetically susceptible individual to certain environmental factors. BA occurs exclusively in the neonatal liver, so variants of genes expressed during hepatobiliary development could affect susceptibility. Genome-wide association studies previously identified a potential region of interest at 2q37. We continued these studies to narrow the region and identify BA susceptibility genes. Methods We searched for copy number variants that were increased among patients with BA (n = 61) compared with healthy individuals (controls; n=5088). After identifying a candidate gene, we investigated expression patterns of orthologues in zebrafish liver and the effects of reducing expression, with morpholino antisense oligonucleotides, on biliary development, gene expression, and signal transduction. Results We observed a statistically significant increase in deletions at 2q37.3 in patients with BA that resulted in deletion of one copy of GPC1 , which encodes glypican 1, a heparan sulfate proteoglycan that regulates Hedgehog signaling and inflammation. Knockdown of gpc1 in zebrafish led to developmental biliary defects. Exposure of the gpc1 morphants to cyclopamine, a Hedgehog antagonist, partially rescued the gpc1 -knockdown phenotype. Injection of zebrafish with recombinant Sonic Hedgehog led to biliary defects similar to those of the gpc1 morphants. Liver samples from patients with BA had reduced levels of apical GPC1 in cholangiocytes compared with samples from controls. Conclusions Based on genetic analysis of patients with BA and zebrafish, GPC1 appears to be a BA susceptibility gene. These findings also support a role for Hedgehog signaling in the pathogenesis of BA.

Effective implementation of precision medicine will be enhanced by a thorough understanding of each patient’s genetic composition to better treat his or her presenting symptoms or mitigate the onset of disease. This ideally includes the sequence information of a complete genome for each individual. At Partners HealthCare Personalized Medicine, we have developed a clinical process for whole genome sequencing (WGS) with application in both healthy individuals and those with disease. In this manuscript, we will describe our bioinformatics strategy to efficiently process and deliver genomic data to geneticists for clinical interpretation. We describe the handling of data from FASTQ to the final variant list for clinical review for the final report. We will also discuss our methodology for validating this workflow and the cost implications of running WGS.

Background & Aims: Alagille syndrome is an autosomal-dominant, multisystem disorder caused primarily by mutations in JAG1, resulting in bile duct paucity, cholestasis, cardiac disease, and other features. Liver disease severity in Alagille syndrome is highly variable, however, factors influencing the hepatic phenotype are unknown. We hypothesized that genetic modifiers may contribute to the variable expressivity of this disorder. Methods: We performed a genome-wide association study in a cohort of Caucasian subjects with known pathogenic JAG1 mutations, comparing patients with mild vs severe liver disease, followed by functional characterization of a candidate locus. Results: We identified a locus that reached suggestive genome-level significance upstream of the thrombospondin 2 (THBS2) gene. THBS2 codes for a secreted matricellular protein that regulates cell proliferation, apoptosis, and angiogenesis, and has been shown to affect Notch signaling. By using a reporter mouse line, we detected thrombospondin 2 expression in bile ducts and periportal regions of the mouse liver. Examination of Thbs2-null mouse livers showed increased microvessels in the portal regions of adult mice. We also showed that thrombospondin 2 interacts with NOTCH1 and NOTCH2 and can inhibit JAG1âNOTCH2 interactions. Conclusions: Based on the genome-wide association study results, thrombospondin 2 localization within bile ducts, and demonstration of interactions of thrombospondin 2 with JAG1 and NOTCH2, we propose that changes in thrombospondin 2 expression may further perturb JAG1âNOTCH2 signaling in patients harboring a JAG1 mutation and lead to a more severe liver phenotype. These results implicate THBS2 as a plausible candidate genetic modifier of liver disease severity in Alagille syndrome. Keywords: JAG1, NOTCH2, Gene Modifier, Cholestasis, Genome-Wide Association Study