Ellen A. Tsai

Evidence From Human and Zebrafish That GPC1 Is a Biliary Atresia Susceptibility Gene

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.

Heterozygous Deletion of FOXA2 Segregates with Disease in a Family with Heterotaxy, Panhypopituitarism, and Biliary Atresia

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.

Homozygosity for the V37I GJB2 mutation in fifteen probands with mild to moderate sensorineural hearing impairment: Further confirmation of pathogenicity and haplotype analysis in Asian populations

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.

Bioinformatics Workflow for Clinical Whole Genome Sequencing at Partners HealthCare Personalized Medicine

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.

THBS2 Is a Candidate Modifier of Liver Disease Severity in Alagille Syndrome

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.

PECONPI: A Novel Software for Uncovering Pathogenic Copy Number Variations in Non-Syndromic Sensorineural Hearing Loss and Other Genetically Heterogeneous Disorders

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 heterogenous 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.

Automated typing of red blood cell and platelet antigens: a whole-genome sequencing study

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.

Characterizing reduced coverage regions through comparison of exome and genome sequencing data across 10 centers

PurposeAs 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 describe a set of standards for variant calling in clinical sequencing applications.MethodsTo 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 10 projects at sequencing centers in the National Human Genome Research Institute/National Cancer Institute Clinical Sequencing Exploratory Research Consortium.ResultsThis 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.ConclusionThis report provides a process road map for clinical sequencing centers looking to perform similar analyses on their data.

DRAW+SneakPeek: Analysis workflow and quality metric management for DNA-seq experiments

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