The current human reference sequence (GRCh38) is a foundation for large-scale sequencing projects. However, recent studies have suggested that GRCh38 may be incomplete and give a suboptimal representation of specific population groups. Here, we performed a de novo assembly of two Swedish genomes that revealed over 10 Mb of sequences absent from the human GRCh38 reference in each individual. Around 6 Mb of these novel sequences (NS) are shared with a Chinese personal genome. The NS are highly repetitive, have an elevated GC-content, and are primarily located in centromeric or telomeric regions. Up to 1 Mb of NS can be assigned to chromosome Y, and large segments are also missing from GRCh38 at chromosomes 14, 17, and 21. Inclusion of NS into the GRCh38 reference radically improves the alignment and variant calling from short-read whole-genome sequencing data at several genomic loci. A re-analysis of a Swedish population-scale sequencing project yields > 75,000 putative novel single nucleotide variants (SNVs) and removes > 10,000 false positive SNV calls per individual, some of which are located in protein coding regions. Our results highlight that the GRCh38 reference is not yet complete and demonstrate that personal genome assemblies from local populations can improve the analysis of short-read whole-genome sequencing data.

Over the past decade, the Database of Genomic Variants (DGV; http://dgv.tcag.ca/) has provided a publicly accessible, comprehensive curated catalogue of structural variation (SV) found in the genomes of control individuals from worldwide populations. Here, we describe updates and new features, which have expanded the utility of DGV for both the basic research and clinical diagnostic communities. The current version of DGV consists of 55 published studies, comprising >2.5 million entries identified in >22,300 genomes. Studies included in DGV are selected from the accessioned data sets in the archival SV databases dbVar (NCBI) and DGVa (EBI), and then further curated for accuracy and validity. The core visualization tool (gbrowse) has been upgraded with additional functions to facilitate data analysis and comparison, and a new query tool has been developed to provide flexible and interactive access to the data. The content from DGV is regularly incorporated into other large-scale genome reference databases and represents a standard data resource for new product and database development, in particular for copy number variation testing in clinical labs. The accurate cataloguing of variants in DGV will continue to enable medical genetics and genome sequencing research.

The starting material for RNA sequencing (RNA-seq) studies is usually total RNA or polyA+ RNA. Both forms of RNA represent heterogeneous pools of RNA molecules at different levels of maturation and processing. Such heterogeneity, in addition to the biases associated with polyA+ purification steps, may influence the analysis, sensitivity and the interpretation of RNA-seq data. We hypothesize that subcellular fractions of RNA may provide a more accurate picture of gene expression.

Intellectual Disability (ID) is a clinically heterogeneous condition that affects 2-3% of population worldwide. In recent years, exome sequencing has been a successful strategy for studies of genetic causes of ID, providing a growing list of both candidate and validated ID genes. In this study, exome sequencing was performed on 28 ID patients in 27 patient-parent trios with the aim to identify de novo variants (DNVs) in known and novel ID associated genes. We report the identification of 25 DNVs out of which five were classified as pathogenic or likely pathogenic. Among these, a two base pair deletion was identified in the PUF60 gene, which is one of three genes in the critical region of the 8q24.3 microdeletion syndrome (Verheij syndrome). Our result adds to the growing evidence that PUF60 is responsible for the majority of the symptoms reported for carriers of a microdeletion across this region. We also report variants in several genes previously not associated with ID, including a de novo missense variant in NAA15. We highlight NAA15 as a novel candidate ID gene based on the vital role of NAA15 in the generation and differentiation of neurons in neonatal brain, the fact that the gene is highly intolerant to loss of function and coding variation, and previously reported DNVs in neurodevelopmental disorders.

Glycosylphosphatidylinositol (GPI) is a glycolipid that tethers more than 150 different proteins to the cell surface. Aberrations in biosynthesis of GPI anchors cause congenital disorders of glycosylation with clinical features including intellectual disability (ID), seizures, and facial dysmorphism. Here, we present two siblings with ID, cerebellar hypoplasia, cerebellar ataxia, early-onset seizures, and minor facial dysmorphology. Using exome sequencing, we identified a homozygous nonsense variant (NM_001127178.1:c.1640G>A, p.Trp547*) in the gene Phosphatidylinositol Glycan Anchor Biosynthesis, Class G (PIGG) in both the patients. Variants in several other GPI anchor synthesis genes lead to a reduced expression of GPI-anchored proteins (GPI-APs) that can be measured by flow cytometry. No significant differences in GPI-APs could be detected in patient granulocytes, consistent with recent findings. However, fibroblasts showed a reduced global level of GPI anchors and of specific GPI-linked markers. These findings suggest that fibroblasts might be more sensitive to pathogenic variants in GPI synthesis pathway and are well suited to screen for GPI-anchor deficiencies. Based on genetic and functional evidence, we confirm that pathogenic variants in PIGG cause an ID syndrome, and we find that loss of function of PIGG is associated with GPI deficiency.

Tourette Syndrome (TS) is a neuropsychiatric disorder in children characterized by motor and verbal tics. Although several genes have been suggested in the etiology of TS, the genetic mechanisms remain poorly understood.

With a draft genome-sequence assembly for the chimpanzee available, it is now possible to perform genome-wide analyses to identify, at a submicroscopic level, structural rearrangements that have occurred between chimpanzees and humans. The goal of this study was to investigate chromosomal regions that are inverted between the chimpanzee and human genomes. Using the net alignments for the builds of the human and chimpanzee genome assemblies, we identified a total of 1,576 putative regions of inverted orientation, covering more than 154 mega-bases of DNA. The DNA segments are distributed throughout the genome and range from 23 base pairs to 62 mega-bases in length. For the 66 inversions more than 25 kilobases (kb) in length, 75% were flanked on one or both sides by (often unrelated) segmental duplications. Using PCR and fluorescence in situ hybridization we experimentally validated 23 of 27 (85%) semi-randomly chosen regions; the largest novel inversion confirmed was 4.3 mega-bases at human Chromosome 7p14. Gorilla was used as an out-group to assign ancestral status to the variants. All experimentally validated inversion regions were then assayed against a panel of human samples and three of the 23 (13%) regions were found to be polymorphic in the human genome. These polymorphic inversions include 730 kb (at 7p22), 13 kb (at 7q11), and 1 kb (at 16q24) fragments with a 5%, 30%, and 48% minor allele frequency, respectively. Our results suggest that inversions are an important source of variation in primate genome evolution. The finding of at least three novel inversion polymorphisms in humans indicates this type of structural variation may be a more common feature of our genome than previously realized.

Significant advances have been made over the past 5 years in mapping and characterizing structural variation in the human genome. Despite this progress, our understanding of inversion variants is still very restricted. While unbalanced variants such as copy number variations can be mapped using array-based approaches, strategies for characterization of inversion variants have been limited and underdeveloped. Traditional cytogenetic approaches have long been able to identify microscopic inversion events, but discovery of submicroscopic events has remained elusive and largely ignored. With the advent of paired-end sequencing approaches, it is now possible to map inversions across the human genome. Based on the paired-end sequencing studies published to date, it is now feasible to make a first map of inversions across the human genome and to use this map to explore the characteristics and distribution of this form of variation. The current map of inversions indicates that many remain to be identified, especially in the smaller size ranges. This review provides an overview of the current knowledge about human inversions and their contribution to human phenotypes. Further characterization of inversions should be considered as an important step towards a deeper understanding of human variation and genome dynamics.

Several genomes have now been sequenced, with millions of genetic variants annotated. While significant progress has been made in mapping single nucleotide polymorphisms (SNPs) and small (<10 bp) insertion/deletions (indels), the annotation of larger structural variants has been less comprehensive. It is still unclear to what extent a typical genome differs from the reference assembly, and the analysis of the genomes sequenced to date have shown varying results for copy number variation (CNV) and inversions.

Insulin degrading enzyme, encoded by IDE, plays a primary role in the degradation of amyloid beta-protein (A beta), the deposition of which in senile plaques is one of the defining hallmarks of Alzheimer's disease (AD). We recently identified haplotypes in a broad linkage disequilibrium (LD) block encompassing IDE that associate with several AD-related quantitative traits. Here, by examining 32 polymorphic markers extending across IDE and testing quantitative measures of plaque density and cognitive function in three independent Swedish AD samples, we have refined the probable position of pathogenic sequences to a 3' region of IDE, with local maximum effects in the proximity of marker rs1887922. To replicate these findings, a subset of variants were examined against measures of brain A beta load in an independent English AD sample, whereby maximum effects were again observed for rs1887922. For both Swedish and English autopsy materials, variation at rs1887922 explained approximately 10% of the total variance in the respective histopathology traits. However, across all clinical materials studied to date, this variant site does not appear to associate directly with disease, suggesting that IDE may affect AD severity rather than risk. Results indicate that alleles of IDE contribute to variability in A beta deposition in the AD brain and suggest that this relationship may have relevance for the degree of cognitive dysfunction in AD patients.

Circular RNAs (circRNAs) were recently discovered as a class of widely expressed noncoding RNA and have been implicated in regulation of gene expression. However, the function of the majority of circRNAs remains unknown. Studies of circRNAs have been hampered by a lack of essential approaches for detection, quantification and visualization. We therefore developed a target-enrichment sequencing method suitable for screening of circRNAs and their linear counterparts in large number of samples. We also applied padlock probes and in situ sequencing to visualize and determine circRNA localization in human brain tissue at subcellular levels. We measured circRNA abundance across different human samples and tissues. Our results highlight the potential of this RNA class to act as a specific diagnostic marker in blood and serum, by detection of circRNAs from genes exclusively expressed in the brain. The powerful and scalable tools we present will enable studies of circRNA function and facilitate screening of circRNA as diagnostic biomarkers.

In recent years, studies have shown higher prevalence of autoantibodies in patients with schizophrenia compared to healthy individuals. This study applies an untargeted and a targeted affinity proteomics approach to explore and characterize the autoantibody repertoire in brain tissues from 73 subjects diagnosed with schizophrenia and 52 control subjects with no psychiatric or neurological disorders. Selected brain tissue lysates were first explored for IgG reactivity on planar microarrays composed of 11,520 protein fragments representing 10,820 unique proteins. Based on these results of ours and other previous studies of autoantibodies related to psychosis, we selected 226 fragments with an average length of 80 amino acids, representing 127 unique proteins. Tissue-based analysis of IgG reactivities using antigen suspension bead arrays was performed in a multiplex and parallel fashion for all 125 subjects. Among the detected autoantigens, higher IgG reactivity in subjects with schizophrenia, as compared to psychiatrically healthy subjects, was found against the glutamate ionotropic receptor NMDA type subunit 2D (anti-GluN2D). In a separate cohort with serum samples from 395 young adults with a wider spectrum of psychiatric disorders, higher levels of serum autoantibodies targeting GluN2D were found when compared to 102 control individuals. By further validating GluN2D and additional potential autoantigens, we will seek insights into how these are associated with severe mental illnesses.

Copy number variation (CNV) plays an important role in human genetic diversity and has been associated with multiple complex disorders. Here we investigate a CNV on chromosome 10q11.22 that spans NPY4R, the gene for the appetite-regulating pancreatic polypeptide receptor Y4. This genomic region has been challenging to map due to multiple repeated elements and its precise organization has not yet been resolved. Previous studies using microarrays were interpreted to show that the most common copy number was 2 per genome.

Chromatin-remodeling factors are required for a wide range of cellular and biological processes including development and cognition, mainly by regulating gene expression. As these functions would predict, deregulation of chromatin-remodeling factors causes various disease syndromes, including neurodevelopmental disorders. Recent reports have linked mutations in several genes coding for chromatin-remodeling factors to intellectual disability (ID). Here, we used exome sequencing and identified a nonsynonymous de novo mutation in BAZ1A (NM_182648.2:c.4043T > G, p.Phe1348Cys), encoding the ATP-utilizing chromatin assembly and remodeling factor 1 (ACF1), in a patient with unexplained ID. ACF1 has been previously reported to bind to the promoter of the vitamin D receptor (VDR)-regulated genes and suppress their expression. Our results show that the patient displays decreased binding of ACF1 to the promoter of the VDR-regulated gene CYP24A1. Using RNA sequencing, we find that the mutation affects the expression of genes involved in several pathways including vitamin D metabolism, Wnt signaling and synaptic formation. RNA sequencing of BAZ1A knockdown cells and Baz1a knockout mice revealed that BAZ1A carry out distinctive functions in different tissues. We also demonstrate that BAZ1A depletion influence the expression of genes important for nervous system development and function. Our data point to an important role for BAZ1A in neurodevelopment, and highlight a possible link for BAZ1A to ID.

Renewed attention has been directed to the functions of the Y chromosome in the central nervous system during early human male development, due to the recent proposed involvement in neurodevelopmental diseases. PCDH11Y and NLGN4Y are of special interest because they belong to gene families involved in cell fate determination and formation of dendrites and axon.

RNA sequencing has become an important method to perform hypothesis-free characterization of global gene expression. One of the limitations of RNA sequencing is that most sequence reads represent highly expressed transcripts, whereas low level transcripts are challenging to detect. To combine the benefits of traditional expression arrays with the advantages of RNA sequencing, we have used whole exome enrichment prior to sequencing of total RNA. We show that whole exome capture can be successfully applied to cDNA to study the transcriptional landscape in human tissues. By introducing the exome enrichment step, we are able to identify transcripts present at very low levels, which are below the level of detection in conventional RNA sequencing. Although the enrichment increases the ability to detect presence of transcripts, it also lowers the accuracy of quantification of expression levels. Our results yield a large number of novel exons and splice isoforms, suggesting that conventional RNA sequencing methods only detect a small fraction of the full transcript diversity. We propose that whole exome enrichment of RNA is a suitable strategy for genome-wide discovery of novel transcripts, alternative splice variants and fusion genes.

Structural variation (copy number variation [CNV] including deletion and duplication, translocation, inversion) of chromosomes has been identified in some individuals with autism spectrum disorder (ASD), but the full etiologic role is unknown. We performed genome-wide assessment for structural abnormalities in 427 unrelated ASD cases via single-nucleotide polymorphism microarrays and karyotyping. With microarrays, we discovered 277 unbalanced CNVs in 44% of ASD families not present in 500 controls (and re-examined in another 1152 controls). Karyotyping detected additional balanced changes. Although most variants were inherited, we found a total of 27 cases with de novo alterations, and in three (11%) of these individuals, two or more new variants were observed. De novo CNVs were found in approximately 7% and approximately 2% of idiopathic families having one child, or two or more ASD siblings, respectively. We also detected 13 loci with recurrent/overlapping CNV in unrelated cases, and at these sites, deletions and duplications affecting the same gene(s) in different individuals and sometimes in asymptomatic carriers were also found. Notwithstanding complexities, our results further implicate the SHANK3-NLGN4-NRXN1 postsynaptic density genes and also identify novel loci at DPP6-DPP10-PCDH9 (synapse complex), ANKRD11, DPYD, PTCHD1, 15q24, among others, for a role in ASD susceptibility. Our most compelling result discovered CNV at 16p11.2 (p = 0.002) (with characteristics of a genomic disorder) at approximately 1% frequency. Some of the ASD regions were also common to mental retardation loci. Structural variants were found in sufficiently high frequency influencing ASD to suggest that cytogenetic and microarray analyses be considered in routine clinical workup.

Chromoanagenesis is a genomic event responsible for the formation of complex structural chromosomal rearrangements (CCRs). Germline chromoanagenesis is rare and the majority of reported cases are associated with an affected phenotype. Here, we report a healthy female carrying two de novo CCRs involving chromosomes 4, 19, 21 and X and chromosomes 7 and 11, respectively, with a total of 137 breakpoint junctions (BPJs). We characterized the CCRs using a hybrid-sequencing approach, combining short-read sequencing, nanopore sequencing, and optical mapping. The results were validated using multiple cytogenetic methods, including fluorescence in situ hybridization, spectral karyotyping, and Sanger sequencing. We identified 137 BPJs, which to our knowledge is the highest number of reported breakpoint junctions in germline chromoanagenesis. We also performed a statistical assessment of the positioning of the breakpoints, revealing a significant enrichment of BPJ-affecting genes (96 intragenic BPJs, 26 genes, p < 0.0001), indicating that the CCRs formed during active transcription of these genes. In addition, we find that the DNA fragments are unevenly and non-randomly distributed across the derivative chromosomes indicating a multistep process of scattering and re-joining of DNA fragments. In summary, we report a new maximum number of BPJs (137) in germline chromoanagenesis. We also show that a hybrid sequencing approach is necessary for the correct characterization of complex CCRs. Through in-depth statistical assessment, it was found that the CCRs most likely was formed through an event resembling chromoplexy-a catastrophic event caused by erroneous transcription factor binding.

Amplification of DNA is required as a mandatory step during library preparation in most targeted sequencing protocols. This can be a critical limitation when targeting regions that are highly repetitive or with extreme guanine-cytosine (GC) content, including repeat expansions associated with human disease. Here, we used an amplification-free protocol for targeted enrichment utilizing the CRISPR/Cas9 system (No-Amp Targeted sequencing) in combination with single molecule, real-time (SMRT) sequencing for studying repeat elements in the huntingtin (HTT) gene, where an expanded CAG repeat is causative for Huntington disease. We also developed a robust data analysis pipeline for repeat element analysis that is independent of alignment of reads to a reference genome. The method was applied to 11 diagnostic blood samples, and for all 22 alleles the resulting CAG repeat count agreed with previous results based on fragment analysis. The amplification-free protocol also allowed for studying somatic variability of repeat elements in our samples, without the interference of PCR stutter. In summary, with No-Amp Targeted sequencing in combination with our analysis pipeline, we could accurately study repeat elements that are difficult to investigate using PCR-based methods.

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.