Characterizing the Major Structural Variant Alleles of the Human Genome.

In order to provide a comprehensive resource for human structural variants (SVs), we generated long-read sequence data and analyzed SVs for fifteen human genomes. We sequence resolved 99,604 insertions, deletions, and inversions including 2,238 (1.6 Mbp) that are shared among all discovery genomes with an additional 13,053 (6.9 Mbp) present in the majority, indicating minor alleles or errors in the reference. Genotyping in 440 additional genomes confirms the most common SVs in unique euchromatin are now sequence resolved. We report a ninefold SV bias toward the last 5 Mbp of human chromosomes with nearly 55% of all VNTRs (variable number of tandem repeats) mapping to this portion of the genome. We identify SVs affecting coding and noncoding regulatory loci improving annotation and interpretation of functional variation. These data provide the framework to construct a canonical human reference and a resource for developing advanced representations capable of capturing allelic diversity.

CpG Island Hypermethylation Mediated by DNMT3A Is a Consequence of AML Progression.

DNMT3A mutations occur in ∼25% of acute myeloid leukemia (AML) patients. The most common mutation, DNMT3AR882H, has dominant negative activity that reduces DNA methylation activity by ∼80% in vitro. To understand the contribution of DNMT3A-dependent methylation to leukemogenesis, we performed whole-genome bisulfite sequencing of primary leukemic and non-leukemic cells in patients with or without DNMT3AR882 mutations. Non-leukemic hematopoietic cells with DNMT3AR882H displayed focal methylation loss, suggesting that hypomethylation antedates AML. Although virtually all AMLs with wild-type DNMT3A displayed CpG island hypermethylation, this change was not associated with gene silencing and was essentially absent in AMLs with DNMT3AR882 mutations. Primary hematopoietic stem cells expanded with cytokines were hypermethylated in a DNMT3A-dependent manner, suggesting that hypermethylation may be a response to, rather than a cause of, cellular proliferation. Our findings suggest that hypomethylation is an initiating phenotype in AMLs with DNMT3AR882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes.

We sequenced the MSY (male-specific region of the Y chromosome) of the C57BL/6J strain of the laboratory mouse Mus musculus. In contrast to theories that Y chromosomes are heterochromatic and gene poor, the mouse MSY is 99.9% euchromatic and contains about 700 protein-coding genes. Only 2% of the MSY derives from the ancestral autosomes that gave rise to the mammalian sex chromosomes. Instead, all but 45 of the MSY's genes belong to three acquired, massively amplified gene families that have no homologs on primate MSYs but do have acquired, amplified homologs on the mouse X chromosome. The complete mouse MSY sequence brings to light dramatic forces in sex chromosome evolution: lineage-specific convergent acquisition and amplification of X-Y gene families, possibly fueled by antagonism between acquired X-Y homologs. The mouse MSY sequence presents opportunities for experimental studies of a sex-specific chromosome in its entirety, in a genetically tractable model organism.

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.

Chromothripsis and human disease: piecing together the shattering process.

The unprecedented resolution of high-throughput genomics has enabled the recent discovery of a phenomenon by which specific regions of the genome are shattered and then stitched together via a single devastating event, referred to as chromothripsis. Potential mechanisms governing this process are now emerging, with implications for our understanding of the role of genomic rearrangements in development and disease.

Evolution of human-specific neural SRGAP2 genes by incomplete segmental duplication.

Gene duplication is an important source of phenotypic change and adaptive evolution. We leverage a haploid hydatidiform mole to identify highly identical sequences missing from the reference genome, confirming that the cortical development gene Slit-Robo Rho GTPase-activating protein 2 (SRGAP2) duplicated three times exclusively in humans. We show that the promoter and first nine exons of SRGAP2 duplicated from 1q32.1 (SRGAP2A) to 1q21.1 (SRGAP2B) ∼3.4 million years ago (mya). Two larger duplications later copied SRGAP2B to chromosome 1p12 (SRGAP2C) and to proximal 1q21.1 (SRGAP2D) ∼2.4 and ∼1 mya, respectively. Sequence and expression analyses show that SRGAP2C is the most likely duplicate to encode a functional protein and is among the most fixed human-specific duplicate genes. Our data suggest a mechanism where incomplete duplication created a novel gene function-antagonizing parental SRGAP2 function-immediately "at birth" 2-3 mya, which is a time corresponding to the transition from Australopithecus to Homo and the beginning of neocortex expansion.

Genomic landscape of non-small cell lung cancer in smokers and never-smokers.

We report the results of whole-genome and transcriptome sequencing of tumor and adjacent normal tissue samples from 17 patients with non-small cell lung carcinoma (NSCLC). We identified 3,726 point mutations and more than 90 indels in the coding sequence, with an average mutation frequency more than 10-fold higher in smokers than in never-smokers. Novel alterations in genes involved in chromatin modification and DNA repair pathways were identified, along with DACH1, CFTR, RELN, ABCB5, and HGF. Deep digital sequencing revealed diverse clonality patterns in both never-smokers and smokers. All validated EFGR and KRAS mutations were present in the founder clones, suggesting possible roles in cancer initiation. Analysis revealed 14 fusions, including ROS1 and ALK, as well as novel metabolic enzymes. Cell-cycle and JAK-STAT pathways are significantly altered in lung cancer, along with perturbations in 54 genes that are potentially targetable with currently available drugs.

Author Correction: Elephant shark genome provides unique insights into gnathostome evolution.

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

A human genome structural variation sequencing resource reveals insights into mutational mechanisms.

Understanding the prevailing mutational mechanisms responsible for human genome structural variation requires uniformity in the discovery of allelic variants and precision in terms of breakpoint delineation. We develop a resource based on capillary end sequencing of 13.8 million fosmid clones from 17 human genomes and characterize the complete sequence of 1054 large structural variants corresponding to 589 deletions, 384 insertions, and 81 inversions. We analyze the 2081 breakpoint junctions and infer potential mechanism of origin. Three mechanisms account for the bulk of germline structural variation: microhomology-mediated processes involving short (2-20 bp) stretches of sequence (28%), nonallelic homologous recombination (22%), and L1 retrotransposition (19%). The high quality and long-range continuity of the sequence reveals more complex mutational mechanisms, including repeat-mediated inversions and gene conversion, that are most often missed by other methods, such as comparative genomic hybridization, single nucleotide polymorphism microarrays, and next-generation sequencing.

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.

Full-length isoform concatenation sequencing to resolve cancer transcriptome complexity.

BACKGROUND: Cancers exhibit complex transcriptomes with aberrant splicing that induces isoform-level differential expression compared to non-diseased tissues. Transcriptomic profiling using short-read sequencing has utility in providing a cost-effective approach for evaluating isoform expression, although short-read assembly displays limitations in the accurate inference of full-length transcripts. Long-read RNA sequencing (Iso-Seq), using the Pacific Biosciences (PacBio) platform, can overcome such limitations by providing full-length isoform sequence resolution which requires no read assembly and represents native expressed transcripts. A constraint of the Iso-Seq protocol is due to fewer reads output per instrument run, which, as an example, can consequently affect the detection of lowly expressed transcripts. To address these deficiencies, we developed a concatenation workflow, PacBio Full-Length Isoform Concatemer Sequencing (PB_FLIC-Seq), designed to increase the number of unique, sequenced PacBio long-reads thereby improving overall detection of unique isoforms. In addition, we anticipate that the increase in read depth will help improve the detection of moderate to low-level expressed isoforms. RESULTS: In sequencing a commercial reference (Spike-In RNA Variants; SIRV) with known isoform complexity we demonstrated a 3.4-fold increase in read output per run and improved SIRV recall when using the PB_FLIC-Seq method compared to the same samples processed with the Iso-Seq protocol. We applied this protocol to a translational cancer case, also demonstrating the utility of the PB_FLIC-Seq method for identifying differential full-length isoform expression in a pediatric diffuse midline glioma compared to its adjacent non-malignant tissue. Our data analysis revealed increased expression of extracellular matrix (ECM) genes within the tumor sample, including an isoform of the Secreted Protein Acidic and Cysteine Rich (SPARC) gene that was expressed 11,676-fold higher than in the adjacent non-malignant tissue. Finally, by using the PB_FLIC-Seq method, we detected several cancer-specific novel isoforms. CONCLUSION: This work describes a concatenation-based methodology for increasing the number of sequenced full-length isoform reads on the PacBio platform, yielding improved discovery of expressed isoforms. We applied this workflow to profile the transcriptome of a pediatric diffuse midline glioma and adjacent non-malignant tissue. Our findings of cancer-specific novel isoform expression further highlight the importance of long-read sequencing for characterization of complex tumor transcriptomes.

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.

Sequence analysis in Bos taurus reveals pervasiveness of X-Y arms races in mammalian lineages.

Studies of Y Chromosome evolution have focused primarily on gene decay, a consequence of suppression of crossing-over with the X Chromosome. Here, we provide evidence that suppression of X-Y crossing-over unleashed a second dynamic: selfish X-Y arms races that reshaped the sex chromosomes in mammals as different as cattle, mice, and men. Using super-resolution sequencing, we explore the Y Chromosome of Bos taurus (bull) and find it to be dominated by massive, lineage-specific amplification of testis-expressed gene families, making it the most gene-dense Y Chromosome sequenced to date. As in mice, an X-linked homolog of a bull Y-amplified gene has become testis-specific and amplified. This evolutionary convergence implies that lineage-specific X-Y coevolution through gene amplification, and the selfish forces underlying this phenomenon, were dominatingly powerful among diverse mammalian lineages. Together with Y gene decay, X-Y arms races molded mammalian sex chromosomes and influenced the course of mammalian evolution.

A novel sialic acid-binding adhesin present in multiple species contributes to the pathogenesis of Infective endocarditis.

Bacterial binding to platelets is a key step in the development of infective endocarditis (IE). Sialic acid, a common terminal carbohydrate on host glycans, is the major receptor for streptococci on platelets. So far, all defined interactions between streptococci and sialic acid on platelets are mediated by serine-rich repeat proteins (SRRPs). However, we identified Streptococcus oralis subsp. oralis IE-isolates that bind sialic acid but lack SRRPs. In addition to binding sialic acid, some SRRP- isolates also bind the cryptic receptor ß-1,4-linked galactose through a yet unknown mechanism. Using comparative genomics, we identified a novel sialic acid-binding adhesin, here named AsaA (associated with sialic acid adhesion A), present in IE-isolates lacking SRRPs. We demonstrated that S. oralis subsp. oralis AsaA is required for binding to platelets in a sialic acid-dependent manner. AsaA comprises a non-repeat region (NRR), consisting of a FIVAR/CBM and two Siglec-like and Unique domains, followed by 31 DUF1542 domains. When recombinantly expressed, Siglec-like and Unique domains competitively inhibited binding of S. oralis subsp. oralis and directly interacted with sialic acid on platelets. We further demonstrated that AsaA impacts the pathogenesis of S. oralis subsp. oralis in a rabbit model of IE. Additionally, we found AsaA orthologues in other IE-causing species and demonstrated that the NRR of AsaA from Gemella haemolysans blocked binding of S. oralis subsp. oralis, suggesting that AsaA contributes to the pathogenesis of multiple IE-causing species. Finally, our findings provide evidence that sialic acid is a key factor for bacterial-platelets interactions in a broader range of species than previously appreciated, highlighting its potential as a therapeutic target.

Orthotopic patient-derived xenografts of paediatric solid tumours.

Paediatric solid tumours arise from endodermal, ectodermal, or mesodermal lineages. Although the overall survival of children with solid tumours is 75%, that of children with recurrent disease is below 30%. To capture the complexity and diversity of paediatric solid tumours and establish new models of recurrent disease, here we develop a protocol to produce orthotopic patient-derived xenografts at diagnosis, recurrence, and autopsy. Tumour specimens were received from 168 patients, and 67 orthotopic patient-derived xenografts were established for 12 types of cancer. The origins of the patient-derived xenograft tumours were reflected in their gene-expression profiles and epigenomes. Genomic profiling of the tumours, including detailed clonal analysis, was performed to determine whether the clonal population in the xenograft recapitulated the patient's tumour. We identified several drug vulnerabilities and showed that the combination of a WEE1 inhibitor (AZD1775), irinotecan, and vincristine can lead to complete response in multiple rhabdomyosarcoma orthotopic patient-derived xenografts tumours in vivo.

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.

Long-read sequence and assembly of segmental duplications.

We have developed a computational method based on polyploid phasing of long sequence reads to resolve collapsed regions of segmental duplications within genome assemblies. Segmental Duplication Assembler (SDA; https://github.com/mvollger/SDA ) constructs graphs in which paralogous sequence variants define the nodes and long-read sequences provide attraction and repulsion edges, enabling the partition and assembly of long reads corresponding to distinct paralogs. We apply it to single-molecule, real-time sequence data from three human genomes and recover 33-79 megabase pairs (Mb) of duplications in which approximately half of the loci are diverged (<99.8%) compared to the reference genome. We show that the corresponding sequence is highly accurate (>99.9%) and that the diverged sequence corresponds to copy-number-variable paralogs that are absent from the human reference genome. Our method can be applied to other complex genomes to resolve the last gene-rich gaps, improve duplicate gene annotation, and better understand copy-number-variant genetic diversity at the base-pair level.

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.