The next-generation sequencing revolution and its impact on genomics.

Genomics is a relatively new scientific discipline, having DNA sequencing as its core technology. As technology has improved the cost and scale of genome characterization over sequencing's 40-year history, the scope of inquiry has commensurately broadened. Massively parallel sequencing has proven revolutionary, shifting the paradigm of genomics to address biological questions at a genome-wide scale. Sequencing now empowers clinical diagnostics and other aspects of medical care, including disease risk, therapeutic identification, and prenatal testing. This Review explores the current state of genomics in the massively parallel sequencing era.

Mono-allelic KCNB2 variants lead to a neurodevelopmental syndrome caused by altered channel inactivation.

Ion channels mediate voltage fluxes or action potentials that are central to the functioning of excitable cells such as neurons. The KCNB family of voltage-gated potassium channels (Kv) consists of two members (KCNB1 and KCNB2) encoded by KCNB1 and KCNB2, respectively. These channels are major contributors to delayed rectifier potassium currents arising from the neuronal soma which modulate overall excitability of neurons. In this study, we identified several mono-allelic pathogenic missense variants in KCNB2, in individuals with a neurodevelopmental syndrome with epilepsy and autism in some individuals. Recurrent dysmorphisms included a broad forehead, synophrys, and digital anomalies. Additionally, we selected three variants where genetic transmission has not been assessed, from two epilepsy studies, for inclusion in our experiments. We characterized channel properties of these variants by expressing them in oocytes of Xenopus laevis and conducting cut-open oocyte voltage clamp electrophysiology. Our datasets indicate no significant change in absolute conductance and conductance-voltage relationships of most disease variants as compared to wild type (WT), when expressed either alone or co-expressed with WT-KCNB2. However, variants c.1141A>G (p.Thr381Ala) and c.641C>T (p.Thr214Met) show complete abrogation of currents when expressed alone with the former exhibiting a left shift in activation midpoint when expressed alone or with WT-KCNB2. The variants we studied, nevertheless, show collective features of increased inactivation shifted to hyperpolarized potentials. We suggest that the effects of the variants on channel inactivation result in hyper-excitability of neurons, which contributes to disease manifestations.

MARK2 variants cause autism spectrum disorder via the downregulation of WNT/ß-catenin signaling pathway.

Microtubule affinity-regulating kinase 2 (MARK2) contributes to establishing neuronal polarity and developing dendritic spines. Although large-scale sequencing studies have associated MARK2 variants with autism spectrum disorder (ASD), the clinical features and variant spectrum in affected individuals with MARK2 variants, early developmental phenotypes in mutant human neurons, and the pathogenic mechanism underlying effects on neuronal development have remained unclear. Here, we report 31 individuals with MARK2 variants and presenting with ASD, other neurodevelopmental disorders, and distinctive facial features. Loss-of-function (LoF) variants predominate (81%) in affected individuals, while computational analysis and in vitro expression assay of missense variants supported the effect of MARK2 loss. Using proband-derived and CRISPR-engineered isogenic induced pluripotent stem cells (iPSCs), we show that MARK2 loss leads to early neuronal developmental and functional deficits, including anomalous polarity and dis-organization in neural rosettes, as well as imbalanced proliferation and differentiation in neural progenitor cells (NPCs). Mark2+/- mice showed abnormal cortical formation and partition and ASD-like behavior. Through the use of RNA sequencing (RNA-seq) and lithium treatment, we link MARK2 loss to downregulation of the WNT/ß-catenin signaling pathway and identify lithium as a potential drug for treating MARK2-associated ASD.

The origin and evolution of mutations in acute myeloid leukemia.

Most mutations in cancer genomes are thought to be acquired after the initiating event, which may cause genomic instability and drive clonal evolution. However, for acute myeloid leukemia (AML), normal karyotypes are common, and genomic instability is unusual. To better understand clonal evolution in AML, we sequenced the genomes of M3-AML samples with a known initiating event (PML-RARA) versus the genomes of normal karyotype M1-AML samples and the exomes of hematopoietic stem/progenitor cells (HSPCs) from healthy people. Collectively, the data suggest that most of the mutations found in AML genomes are actually random events that occurred in HSPCs before they acquired the initiating mutation; the mutational history of that cell is "captured" as the clone expands. In many cases, only one or two additional, cooperating mutations are needed to generate the malignant founding clone. Cells from the founding clone can acquire additional cooperating mutations, yielding subclones that can contribute to disease progression and/or relapse.

Monoallelic de novo variants in DDX17 cause a neurodevelopmental disorder.

DDX17 is an RNA helicase shown to be involved in critical processes during the early phases of neuronal differentiation. Globally, we compiled a case-series of 11 patients with neurodevelopmental phenotypes harbouring de novo monoallelic variants in DDX17. All 11 patients in our case series had a neurodevelopmental phenotype, whereby intellectual disability, delayed speech and language, and motor delay predominated. We performed in utero cortical electroporation in the brain of developing mice, assessing axon complexity and outgrowth of electroporated neurons, comparing wild-type and Ddx17 knockdown. We then undertook ex vivo cortical electroporation on neuronal progenitors to quantitatively assess axonal development at a single cell resolution. Mosaic ddx17 crispants and heterozygous knockouts in Xenopus tropicalis were generated for assessment of morphology, behavioural assays, and neuronal outgrowth measurements. We further undertook transcriptomic analysis of neuroblastoma SH-SY5Y cells, to identify differentially expressed genes in DDX17-KD cells compared to controls. Knockdown of Ddx17 in electroporated mouse neurons in vivo showed delayed neuronal migration as well as decreased cortical axon complexity. Mouse primary cortical neurons revealed reduced axon outgrowth upon knockdown of Ddx17 in vitro. The axon outgrowth phenotype was replicated in crispant ddx17 tadpoles and in heterozygotes. Heterozygous tadpoles had clear neurodevelopmental defects and showed an impaired neurobehavioral phenotype. Transcriptomic analysis identified a statistically significant number of differentially expressed genes involved in neurodevelopmental processes in DDX17-KD cells compared to control cells. We have identified potential neurodevelopment disease-causing variants in a gene not previously associated with genetic disease, DDX17. We provide evidence for the role of the gene in neurodevelopment in both mammalian and non-mammalian species and in controlling the expression of key neurodevelopment genes.

Exome sequencing of Finnish isolates enhances rare-variant association power.

Exome-sequencing studies have generally been underpowered to identify deleterious alleles with a large effect on complex traits as such alleles are mostly rare. Because the population of northern and eastern Finland has expanded considerably and in isolation following a series of bottlenecks, individuals of these populations have numerous deleterious alleles at a relatively high frequency. Here, using exome sequencing of nearly 20,000 individuals from these regions, we investigate the role of rare coding variants in clinically relevant quantitative cardiometabolic traits. Exome-wide association studies for 64 quantitative traits identified 26 newly associated deleterious alleles. Of these 26 alleles, 19 are either unique to or more than 20 times more frequent in Finnish individuals than in other Europeans and show geographical clustering comparable to Mendelian disease mutations that are characteristic of the Finnish population. We estimate that sequencing studies of populations without this unique history would require hundreds of thousands to millions of participants to achieve comparable association power.

Author Correction: Exome sequencing of Finnish isolates enhances rare-variant association power.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

A novel IKZF1 variant in a family with autosomal dominant CVID: A case for expanding exon coverage in inborn errors of immunity.

Common variable immune deficiency (CVID) is a heterogenous group of disorders characterized by varying degrees of hypogammaglobulinemia, recurrent infections, and autoimmunity. Currently, pathogenic variants are identified in approximately 20-30% of CVID cases. Here we report a 3-generation family with autosomal dominant Common Variable Immunodeficiency (CVID) diagnosed in 9 affected individuals. Although primary immune deficiency panels and exome sequencing were non-diagnostic, whole genome sequencing revealed a novel, pathogenic c.499C > T: p.His167Tyr variant in IKZF1, a critical regulator of B cell development. Functional testing done through pericentromeric heterochromatin localization and light shift chemiluminescent electrophoretic mobility shift assay confirmed the variant's deleterious effect via a haploinsufficiency mechanism. Our findings expand the spectrum of known IKZF1 mutations and contribute to a more comprehensive understanding of CVID's genetic heterogeneity. Furthermore, this case underscores the importance of considering whole genome sequencing for comprehensive genetic diagnosis when concern for a monogenic inborn errors of immunity is high.

Cerebral organoids containing an AUTS2 missense variant model microcephaly.

Variants in the AUTS2 gene are associated with a broad spectrum of neurological conditions characterized by intellectual disability, microcephaly, and congenital brain malformations. Here, we use a human cerebral organoid model to investigate the pathophysiology of a heterozygous de novo missense AUTS2 variant identified in a patient with multiple neurological impairments including primary microcephaly and profound intellectual disability. Proband cerebral organoids exhibit reduced growth, deficits in neural progenitor cell (NPC) proliferation and disrupted NPC polarity within ventricular zone-like regions compared to control cerebral organoids. We used CRISPR-Cas9-mediated gene editing to correct this variant and demonstrate rescue of impaired organoid growth and NPC proliferative deficits. Single-cell RNA sequencing revealed a marked reduction of G1/S transition gene expression and alterations in WNT-ß-catenin signalling within proband NPCs, uncovering a novel role for AUTS2 in NPCs during human cortical development. Collectively, these results underscore the value of cerebral organoids to investigate molecular mechanisms underlying AUTS2 syndrome.

Infantile metastatic ependymoma with a novel molecular profile and favorable outcome to intensive chemotherapy without irradiation: Case-based review.

Ependymal tumors are the third most common brain tumor under 14 years old. Even though metastatic disease is a rare event, it affects mostly young children and carries an adverse prognosis. The factors associated with dissemination and the best treatment approach have not yet been established and there is limited published data on how to manage metastatic disease, especially in patients under 3 years of age. We provide a review of the literature on clinical characteristics and radiation-sparing treatments for metastatic ependymoma in children under 3 years of age treated. The majority (73%) of the identified cases were above 12 months old and had the PF as the primary site at diagnosis. Chemotherapy-based approaches, in different regimens, were used with radiation reserved for progression or relapse. The prognosis varied among the studies, with an average of 50%-58% overall survival. This study also describes the case of a 7-month-old boy with metastatic posterior fossa (PF) ependymoma, for whom we identified a novel SPECC1L-RAF1 gene fusion using a patient-centric comprehensive molecular profiling protocol. The patient was successfully treated with intensive induction chemotherapy followed by high-dose chemotherapy and autologous hematopoietic progenitor cell rescue (AuHSCR). Currently, the patient is in continuous remission 5 years after his diagnosis, without radiation therapy. The understanding of the available therapeutic approaches may assist physicians in their management of such patients. This report also opens the perspective of newly identified molecular alterations in metastatic ependymomas that might drive more chemo-sensitive tumors.

Mutations in the Kinesin-2 Motor KIF3B Cause an Autosomal-Dominant Ciliopathy.

Kinesin-2 enables ciliary assembly and maintenance as an anterograde intraflagellar transport (IFT) motor. Molecular motor activity is driven by a heterotrimeric complex comprised of KIF3A and KIF3B or KIF3C plus one non-motor subunit, KIFAP3. Using exome sequencing, we identified heterozygous KIF3B variants in two unrelated families with hallmark ciliopathy phenotypes. In the first family, the proband presents with hepatic fibrosis, retinitis pigmentosa, and postaxial polydactyly; he harbors a de novo c.748G>C (p.Glu250Gln) variant affecting the kinesin motor domain encoded by KIF3B. The second family is a six-generation pedigree affected predominantly by retinitis pigmentosa. Affected individuals carry a heterozygous c.1568T>C (p.Leu523Pro) KIF3B variant segregating in an autosomal-dominant pattern. We observed a significant increase in primary cilia length in vitro in the context of either of the two mutations while variant KIF3B proteins retained stability indistinguishable from wild type. Furthermore, we tested the effects of KIF3B mutant mRNA expression in the developing zebrafish retina. In the presence of either missense variant, rhodopsin was sequestered to the photoreceptor rod inner segment layer with a concomitant increase in photoreceptor cilia length. Notably, impaired rhodopsin trafficking is also characteristic of recessive KIF3B models as exemplified by an early-onset, autosomal-recessive, progressive retinal degeneration in Bengal cats; we identified a c.1000G>A (p.Ala334Thr) KIF3B variant by genome-wide association study and whole-genome sequencing. Together, our genetic, cell-based, and in vivo modeling data delineate an autosomal-dominant syndromic retinal ciliopathy in humans and suggest that multiple KIF3B pathomechanisms can impair kinesin-driven ciliary transport in the photoreceptor.

Dominant-negative variants in CBX1 cause a neurodevelopmental disorder.

PURPOSE: This study aimed to establish variants in CBX1, encoding heterochromatin protein 1ß (HP1ß), as a cause of a novel syndromic neurodevelopmental disorder. METHODS: Patients with CBX1 variants were identified, and clinician researchers were connected using GeneMatcher and physician referrals. Clinical histories were collected from each patient. To investigate the pathogenicity of identified variants, we performed in vitro cellular assays and neurobehavioral and cytological analyses of neuronal cells obtained from newly generated Cbx1 mutant mouse lines. RESULTS: In 3 unrelated individuals with developmental delay, hypotonia, and autistic features, we identified heterozygous de novo variants in CBX1. The identified variants were in the chromodomain, the functional domain of HP1ß, which mediates interactions with chromatin. Cbx1 chromodomain mutant mice displayed increased latency-to-peak response, suggesting the possibility of synaptic delay or myelination deficits. Cytological and chromatin immunoprecipitation experiments confirmed the reduction of mutant HP1ß binding to heterochromatin, whereas HP1ß interactome analysis demonstrated that the majority of HP1ß-interacting proteins remained unchanged between the wild-type and mutant HP1ß. CONCLUSION: These collective findings confirm the role of CBX1 in developmental disabilities through the disruption of HP1ß chromatin binding during neurocognitive development. Because HP1ß forms homodimers and heterodimers, mutant HP1ß likely sequesters wild-type HP1ß and other HP1 proteins, exerting dominant-negative effects.

Missense MED12 variants in 22 males with intellectual disability: From nonspecific symptoms to complete syndromes.

We describe the phenotype of 22 male patients (20 probands) carrying a hemizygous missense variant in MED12. The phenotypic spectrum is very broad ranging from nonspecific intellectual disability (ID) to the three well-known syndromes: Opitz-Kaveggia syndrome, Lujan-Fryns syndrome, or Ohdo syndrome. The identified variants were randomly distributed throughout the gene (p = 0.993, χ2 test), but mostly outside the functional domains (p = 0.004; χ2 test). Statistical analyses did not show a correlation between the MED12-related phenotypes and the locations of the variants (p = 0.295; Pearson correlation), nor the protein domain involved (p = 0.422; Pearson correlation). In conclusion, establishing a genotype-phenotype correlation in MED12-related diseases remains challenging. Therefore, we think that patients with a causative MED12 variant are currently underdiagnosed due to the broad patients' clinical presentations.

Long-read whole genome sequencing reveals HOXD13 alterations in synpolydactyly.

Synpolydactyly 1, also called syndactyly type II (SDTY2), is a genetic limb malformation characterized by polydactyly with syndactyly involving the webbing of the third and fourth fingers, and the fourth and fifth toes. It is caused by heterozygous alterations in HOXD13 with incomplete penetrance and phenotypic variability. In our study, a five-generation family with an SPD phenotype was enrolled in our Rare Disease Genomics Protocol. A comprehensive examination of three generations using Illumina short-read whole-genome sequencing (WGS) did not identify any causative variants. Subsequent WGS using Pacific Biosciences (PacBio) long-read HiFi Circular Consensus Sequencing (CCS) revealed a heterozygous 27-bp duplication in the polyalanine tract of HOXD13. Sanger sequencing of all available family members confirmed that the variant segregates with affected individuals. Reanalysis of an unrelated family with a similar SPD phenotype uncovered a 21-bp (7-alanine) duplication in the same region of HOXD13. Although ExpansionHunter identified these events in most individuals in a retrospective analysis, low sequence coverage due to high GC content in the HOXD13 polyalanine tract makes detection of these events challenging. Our findings highlight the value of long-read WGS in elucidating the molecular etiology of congenital limb malformation disorders.

Expanding the phenotypic and molecular spectrum of NFS1-related disorders that cause functional deficiencies in mitochondrial and cytosolic iron-sulfur cluster containing enzymes.

Iron-sulfur cluster proteins are involved in critical functions for gene expression regulation and mitochondrial bioenergetics including the oxidative phosphorylation system. The c.215G>A p.(Arg72Gln) variant in NFS1 has been previously reported to cause infantile mitochondrial complex II and III deficiency. We describe three additional unrelated patients with the same missense variant. Two infants with the same homozygous variant presented with hypotonia, weakness and lactic acidosis, and one patient with compound heterozygous p.(Arg72Gln) and p.(Arg412His) variants presented as a young adult with gastrointestinal symptoms and fatigue. Skeletal muscle biopsy from patients 1 and 3 showed abnormal mitochondrial morphology, and functional analyses demonstrated decreased activity in respiratory chain complex II and variably in complexes I and III. We found decreased mitochondrial and cytosolic aconitase activities but only mildly affected lipoylation of pyruvate dehydrogenase and 2-oxoglutarate dehydrogenase enzymes. Our studies expand the phenotypic spectrum and provide further evidence for the pathogenicity and functional sequelae of NFS1-related disorders with disturbances in both mitochondrial and cytosolic iron-sulfur cluster containing enzymes.

Case report and review of the literature: immune dysregulation in a large familial cohort due to a novel pathogenic RELA variant.

OBJECTIVE: To explore and define the molecular cause(s) of a multi-generational kindred affected by Bechet's-like mucocutaneous ulcerations and immune dysregulation. METHODS: Whole genome sequencing and confirmatory Sanger sequencing were performed. Components of the NFκB pathway were quantified by immunoblotting, and function was assessed by cytokine production (IL-6, TNF-α, IL-1ß) after lipopolysaccharide (LPS) stimulation. Detailed immunophenotyping of T-cell and B-cell subsets was performed in four patients from this cohort. RESULTS: A novel variant in the RELA gene, p. Tyr349LeufsTer13, was identified. This variant results in premature truncation of the protein before the serine (S) 536 residue, a key phosphorylation site, resulting in enhanced degradation of the p65 protein. Immunoblotting revealed significantly decreased phosphorylated [p]p65 and pIκBα. The decrease in [p]p65 may suggest reduced heterodimer formation between p50/p65 (NFκB1/RelA). Immunophenotyping revealed decreased naïve T cells, increased memory T cells, and expanded senescent T-cell populations in one patient (P1). P1 also had substantially higher IL-6 and TNF-α levels post-stimulation compared with the other three patients. CONCLUSION: Family members with this novel RELA variant have a clinical phenotype similar to other reported RELA cases with predominant chronic mucocutaneous ulceration; however, the clinical phenotype broadens to include Behçet's syndrome and IBD. Here we describe the clinical, immunological and genetic evaluation of a large kindred to further expand identification of patients with autosomal dominant RELA deficiency, facilitating earlier diagnosis and intervention. The functional impairment of the canonical NFκB pathway suggests that this variant is causal for the clinical phenotype in these patients.

Detection of brain somatic variation in epilepsy-associated developmental lesions.

OBJECTIVE: Epilepsy-associated developmental lesions, including malformations of cortical development and low-grade developmental tumors, represent a major cause of drug-resistant seizures requiring surgical intervention in children. Brain-restricted somatic mosaicism has been implicated in the genetic etiology of these lesions; however, many contributory genes remain unidentified. METHODS: We enrolled 50 children who were undergoing epilepsy surgery into a translational research study. Resected tissue was divided for clinical neuropathologic evaluation and genomic analysis. We performed exome and RNA sequencing to identify somatic variation and we confirmed our findings using high-depth targeted DNA sequencing. RESULTS: We uncovered candidate disease-causing somatic variation affecting 28 patients (56%), as well as candidate germline variants affecting 4 patients (8%). In agreement with previous studies, we identified somatic variation affecting solute carrier family 35 member A2 (SLC35A2) and mechanistic target of rapamycin kinase (MTOR) pathway genes in patients with focal cortical dysplasia. Somatic gains of chromosome 1q were detected in 30% (3 of 10) of patients with Type I focal cortical dysplasia (FCD)s. Somatic variation in mitogen-activated protein kinase (MAPK) pathway genes (i.e., fibroblast growth factor receptor 1 [FGFR1], FGFR2, B-raf proto-oncogene, serine/threonine kinase [BRAF], and KRAS proto-oncogene, GTPase [KRAS]) was associated with low-grade epilepsy-associated developmental tumors. RNA sequencing enabled the detection of somatic structural variation that would have otherwise been missed, and which accounted for more than one-half of epilepsy-associated tumor diagnoses. Sampling across multiple anatomic regions revealed that somatic variant allele fractions vary widely within epileptogenic tissue. Finally, we identified putative disease-causing variants in genes not yet associated with focal cortical dysplasia. SIGNIFICANCE: These results further elucidate the genetic basis of structural brain abnormalities leading to focal epilepsy in children and point to new candidate disease genes.

PTEN somatic mutations contribute to spectrum of cerebral overgrowth.

Phosphatase and tensin homologue (PTEN) regulates cell growth and survival through inhibition of the mammalian target of rapamycin (MTOR) signalling pathway. Germline genetic variation of PTEN is associated with autism, macrocephaly and PTEN hamartoma tumour syndromes. The effect of developmental PTEN somatic mutations on nervous system phenotypes is not well understood, although brain somatic mosaicism of MTOR pathway genes is an emerging cause of cortical dysplasia and epilepsy in the paediatric population. Here we report two somatic variants of PTEN affecting a single patient presenting with intractable epilepsy and hemimegalencephaly that varied in clinical severity throughout the left cerebral hemisphere. High-throughput sequencing analysis of affected brain tissue identified two somatic variants in PTEN. The first variant was present in multiple cell lineages throughout the entire hemisphere and associated with mild cerebral overgrowth. The second variant was restricted to posterior brain regions and affected the opposite PTEN allele, resulting in a segmental region of more severe malformation, and the only neurons in which it was found by single-nuclei RNA-sequencing had a unique disease-related expression profile. This study reveals brain mosaicism of PTEN as a disease mechanism of hemimegalencephaly and furthermore demonstrates the varying effects of single- or bi-allelic disruption of PTEN on cortical phenotypes.

Discovery of clinically relevant fusions in pediatric cancer.

BACKGROUND: Pediatric cancers typically have a distinct genomic landscape when compared to adult cancers and frequently carry somatic gene fusion events that alter gene expression and drive tumorigenesis. Sensitive and specific detection of gene fusions through the analysis of next-generation-based RNA sequencing (RNA-Seq) data is computationally challenging and may be confounded by low tumor cellularity or underlying genomic complexity. Furthermore, numerous computational tools are available to identify fusions from supporting RNA-Seq reads, yet each algorithm demonstrates unique variability in sensitivity and precision, and no clearly superior approach currently exists. To overcome these challenges, we have developed an ensemble fusion calling approach to increase the accuracy of identifying fusions. RESULTS: Our Ensemble Fusion (EnFusion) approach utilizes seven fusion calling algorithms: Arriba, CICERO, FusionMap, FusionCatcher, JAFFA, MapSplice, and STAR-Fusion, which are packaged as a fully automated pipeline using Docker and Amazon Web Services (AWS) serverless technology. This method uses paired end RNA-Seq sequence reads as input, and the output from each algorithm is examined to identify fusions detected by a consensus of at least three algorithms. These consensus fusion results are filtered by comparison to an internal database to remove likely artifactual fusions occurring at high frequencies in our internal cohort, while a "known fusion list" prevents failure to report known pathogenic events. We have employed the EnFusion pipeline on RNA-Seq data from 229 patients with pediatric cancer or blood disorders studied under an IRB-approved protocol. The samples consist of 138 central nervous system tumors, 73 solid tumors, and 18 hematologic malignancies or disorders. The combination of an ensemble fusion-calling pipeline and a knowledge-based filtering strategy identified 67 clinically relevant fusions among our cohort (diagnostic yield of 29.3%), including RBPMS-MET, BCAN-NTRK1, and TRIM22-BRAF fusions. Following clinical confirmation and reporting in the patient's medical record, both known and novel fusions provided medically meaningful information. CONCLUSIONS: The EnFusion pipeline offers a streamlined approach to discover fusions in cancer, at higher levels of sensitivity and accuracy than single algorithm methods. Furthermore, this method accurately identifies driver fusions in pediatric cancer, providing clinical impact by contributing evidence to diagnosis and, when appropriate, indicating targeted therapies.

Hypomorphic alleles pose challenges in rare disease genomic variant interpretation.

Exon skipping associated with an ATP7B intronic variant in a patient with Wilson's disease. (A) Sashimi plot visualization of aligned RNA sequencing data from proband liver tissue at ATP7B exons 14-13-12. The red track shows traditional RNA-seq data; the blue track shows RNA-seq enriched with exon capture (cDNA-cap) which achieves higher depth of protein-coding transcripts. The histogram indicates overall sequencing depth while arcs tabulate the number of junction-spanning reads supporting exon pairs. (B) The domain structure (top) and exon structure (bottom) of ATP7B. Loss of exon 13 (dashed box) would remove a transmembrane domain and disrupt the first phosphorylation domain.