Sequencing by avidity enables high accuracy with low reagent consumption.

We present avidity sequencing, a sequencing chemistry that separately optimizes the processes of stepping along a DNA template and that of identifying each nucleotide within the template. Nucleotide identification uses multivalent nucleotide ligands on dye-labeled cores to form polymerase-polymer-nucleotide complexes bound to clonal copies of DNA targets. These polymer-nucleotide substrates, termed avidites, decrease the required concentration of reporting nucleotides from micromolar to nanomolar and yield negligible dissociation rates. Avidity sequencing achieves high accuracy, with 96.2% and 85.4% of base calls having an average of one error per 1,000 and 10,000 base pairs, respectively. We show that the average error rate of avidity sequencing remained stable following a long homopolymer.

An automated 13.5 hour system for scalable diagnosis and acute management guidance for genetic diseases.

While many genetic diseases have effective treatments, they frequently progress rapidly to severe morbidity or mortality if those treatments are not implemented immediately. Since front-line physicians frequently lack familiarity with these diseases, timely molecular diagnosis may not improve outcomes. Herein we describe Genome-to-Treatment, an automated, virtual system for genetic disease diagnosis and acute management guidance. Diagnosis is achieved in 13.5 h by expedited whole genome sequencing, with superior analytic performance for structural and copy number variants. An expert panel adjudicated the indications, contraindications, efficacy, and evidence-of-efficacy of 9911 drug, device, dietary, and surgical interventions for 563 severe, childhood, genetic diseases. The 421 (75%) diseases and 1527 (15%) effective interventions retained are integrated with 13 genetic disease information resources and appended to diagnostic reports ( https://gtrx.radygenomiclab.com ). This system provided correct diagnoses in four retrospectively and two prospectively tested infants. The Genome-to-Treatment system facilitates optimal outcomes in children with rapidly progressive genetic diseases.

ExpansionHunter Denovo: a computational method for locating known and novel repeat expansions in short-read sequencing data.

Repeat expansions are responsible for over 40 monogenic disorders, and undoubtedly more pathogenic repeat expansions remain to be discovered. Existing methods for detecting repeat expansions in short-read sequencing data require predefined repeat catalogs. Recent discoveries emphasize the need for methods that do not require pre-specified candidate repeats. To address this need, we introduce ExpansionHunter Denovo, an efficient catalog-free method for genome-wide repeat expansion detection. Analysis of real and simulated data shows that our method can identify large expansions of 41 out of 44 pathogenic repeats, including nine recently reported non-reference repeat expansions not discoverable via existing methods.

Genome sequencing in persistently unsolved white matter disorders.

Genetic white matter disorders have heterogeneous etiologies and overlapping clinical presentations. We performed a study of the diagnostic efficacy of genome sequencing in 41 unsolved cases with prior exome sequencing, resolving an additional 14 from an historical cohort (n = 191). Reanalysis in the context of novel disease-associated genes and improved variant curation and annotation resolved 64% of cases. The remaining diagnoses were directly attributable to genome sequencing, including cases with small and large copy number variants (CNVs) and variants in deep intronic and technically difficult regions. Genome sequencing, in combination with other methodologies, achieved a diagnostic yield of 85% in this retrospective cohort.

The genome-wide multi-layered architecture of chromosome pairing in early Drosophila embryos.

Genome organization involves cis and trans chromosomal interactions, both implicated in gene regulation, development, and disease. Here, we focus on trans interactions in Drosophila, where homologous chromosomes are paired in somatic cells from embryogenesis through adulthood. We first address long-standing questions regarding the structure of embryonic homolog pairing and, to this end, develop a haplotype-resolved Hi-C approach to minimize homolog misassignment and thus robustly distinguish trans-homolog from cis contacts. This computational approach, which we call Ohm, reveals pairing to be surprisingly structured genome-wide, with trans-homolog domains, compartments, and interaction peaks, many coinciding with analogous cis features. We also find a significant genome-wide correlation between pairing, transcription during zygotic genome activation, and binding of the pioneer factor Zelda. Our findings reveal a complex, highly structured organization underlying homolog pairing, first discovered a century ago in Drosophila. Finally, we demonstrate the versatility of our haplotype-resolved approach by applying it to mammalian embryos.

Highly structured homolog pairing reflects functional organization of the Drosophila genome.

Trans-homolog interactions have been studied extensively in Drosophila, where homologs are paired in somatic cells and transvection is prevalent. Nevertheless, the detailed structure of pairing and its functional impact have not been thoroughly investigated. Accordingly, we generated a diploid cell line from divergent parents and applied haplotype-resolved Hi-C, showing that homologs pair with varying precision genome-wide, in addition to establishing trans-homolog domains and compartments. We also elucidate the structure of pairing with unprecedented detail, observing significant variation across the genome and revealing at least two forms of pairing: tight pairing, spanning contiguous small domains, and loose pairing, consisting of single larger domains. Strikingly, active genomic regions (A-type compartments, active chromatin, expressed genes) correlated with tight pairing, suggesting that pairing has a functional implication genome-wide. Finally, using RNAi and haplotype-resolved Hi-C, we show that disruption of pairing-promoting factors results in global changes in pairing, including the disruption of some interaction peaks.

Publisher Correction: Heterochromatin drives compartmentalization of inverted and conventional nuclei.

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

Heterochromatin drives compartmentalization of inverted and conventional nuclei.

The nucleus of mammalian cells displays a distinct spatial segregation of active euchromatic and inactive heterochromatic regions of the genome1,2. In conventional nuclei, microscopy shows that euchromatin is localized in the nuclear interior and heterochromatin at the nuclear periphery1,2. Genome-wide chromosome conformation capture (Hi-C) analyses show this segregation as a plaid pattern of contact enrichment within euchromatin and heterochromatin compartments3, and depletion between them. Many mechanisms for the formation of compartments have been proposed, such as attraction of heterochromatin to the nuclear lamina2,4, preferential attraction of similar chromatin to each other1,4-12, higher levels of chromatin mobility in active chromatin13-15 and transcription-related clustering of euchromatin16,17. However, these hypotheses have remained inconclusive, owing to the difficulty of disentangling intra-chromatin and chromatin-lamina interactions in conventional nuclei18. The marked reorganization of interphase chromosomes in the inverted nuclei of rods in nocturnal mammals19,20 provides an opportunity to elucidate the mechanisms that underlie spatial compartmentalization. Here we combine Hi-C analysis of inverted rod nuclei with microscopy and polymer simulations. We find that attractions between heterochromatic regions are crucial for establishing both compartmentalization and the concentric shells of pericentromeric heterochromatin, facultative heterochromatin and euchromatin in the inverted nucleus. When interactions between heterochromatin and the lamina are added, the same model recreates the conventional nuclear organization. In addition, our models allow us to rule out mechanisms of compartmentalization that involve strong euchromatin interactions. Together, our experiments and modelling suggest that attractions between heterochromatic regions are essential for the phase separation of the active and inactive genome in inverted and conventional nuclei, whereas interactions of the chromatin with the lamina are necessary to build the conventional architecture from these segregated phases.

ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions.

SUMMARY: We describe a novel computational method for genotyping repeats using sequence graphs. This method addresses the long-standing need to accurately genotype medically important loci containing repeats adjacent to other variants or imperfect DNA repeats such as polyalanine repeats. Here we introduce a new version of our repeat genotyping software, ExpansionHunter, that uses this method to perform targeted genotyping of a broad class of such loci. AVAILABILITY AND IMPLEMENTATION: ExpansionHunter is implemented in C++ and is available under the Apache License Version 2.0. The source code, documentation, and Linux/macOS binaries are available at https://github.com/Illumina/ExpansionHunter/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Measuring the reproducibility and quality of Hi-C data.

BACKGROUND: Hi-C is currently the most widely used assay to investigate the 3D organization of the genome and to study its role in gene regulation, DNA replication, and disease. However, Hi-C experiments are costly to perform and involve multiple complex experimental steps; thus, accurate methods for measuring the quality and reproducibility of Hi-C data are essential to determine whether the output should be used further in a study. RESULTS: Using real and simulated data, we profile the performance of several recently proposed methods for assessing reproducibility of population Hi-C data, including HiCRep, GenomeDISCO, HiC-Spector, and QuASAR-Rep. By explicitly controlling noise and sparsity through simulations, we demonstrate the deficiencies of performing simple correlation analysis on pairs of matrices, and we show that methods developed specifically for Hi-C data produce better measures of reproducibility. We also show how to use established measures, such as the ratio of intra- to interchromosomal interactions, and novel ones, such as QuASAR-QC, to identify low-quality experiments. CONCLUSIONS: In this work, we assess reproducibility and quality measures by varying sequencing depth, resolution and noise levels in Hi-C data from 13 cell lines, with two biological replicates each, as well as 176 simulated matrices. Through this extensive validation and benchmarking of Hi-C data, we describe best practices for reproducibility and quality assessment of Hi-C experiments. We make all software publicly available at http://github.com/kundajelab/3DChromatin_ReplicateQC to facilitate adoption in the community.

Higher-Order Organization Principles of Pre-translational mRNPs.

Compared to noncoding RNAs (ncRNAs), such as rRNAs and ribozymes, for which high-resolution structures abound, little is known about the tertiary structures of mRNAs. In eukaryotic cells, newly made mRNAs are packaged with proteins in highly compacted mRNA particles (mRNPs), but the manner of this mRNA compaction is unknown. Here, we developed and implemented RIPPLiT (RNA immunoprecipitation and proximity ligation in tandem), a transcriptome-wide method for probing the 3D conformations of RNAs stably associated with defined proteins, in this case, exon junction complex (EJC) core factors. EJCs multimerize with other mRNP components to form megadalton-sized complexes that protect large swaths of newly synthesized mRNAs from endonuclease digestion. Unlike ncRNPs, wherein strong locus-specific structures predominate, mRNPs behave more like flexible polymers. Polymer analysis of proximity ligation data for hundreds of mRNA species demonstrates that nascent and pre-translational mammalian mRNAs are compacted by their associated proteins into linear rod-like structures.

Integrative detection and analysis of structural variation in cancer genomes.

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Detection of long repeat expansions from PCR-free whole-genome sequence data.

Identifying large expansions of short tandem repeats (STRs), such as those that cause amyotrophic lateral sclerosis (ALS) and fragile X syndrome, is challenging for short-read whole-genome sequencing (WGS) data. A solution to this problem is an important step toward integrating WGS into precision medicine. We developed a software tool called ExpansionHunter that, using PCR-free WGS short-read data, can genotype repeats at the locus of interest, even if the expanded repeat is larger than the read length. We applied our algorithm to WGS data from 3001 ALS patients who have been tested for the presence of the C9orf72 repeat expansion with repeat-primed PCR (RP-PCR). Compared against this truth data, ExpansionHunter correctly classified all (212/212, 95% CI [0.98, 1.00]) of the expanded samples as either expansions (208) or potential expansions (4). Additionally, 99.9% (2786/2789, 95% CI [0.997, 1.00]) of the wild-type samples were correctly classified as wild type by this method with the remaining three samples identified as possible expansions. We further applied our algorithm to a set of 152 samples in which every sample had one of eight different pathogenic repeat expansions, including those associated with fragile X syndrome, Friedreich's ataxia, and Huntington's disease, and correctly flagged all but one of the known repeat expansions. Thus, ExpansionHunter can be used to accurately detect known pathogenic repeat expansions and provides researchers with a tool that can be used to identify new pathogenic repeat expansions.

X-linked hypomyelination with spondylometaphyseal dysplasia (H-SMD) associated with mutations in AIFM1.

An X-linked condition characterized by the combination of hypomyelinating leukodystrophy and spondylometaphyseal dysplasia (H-SMD) has been observed in only four families, with linkage to Xq25-27, and recent genetic characterization in two families with a common AIFM1 mutation. In our study, 12 patients (6 families) with H-SMD were identified and underwent comprehensive assessment accompanied by whole-exome sequencing (WES). Pedigree analysis in all families was consistent with X-linked recessive inheritance. Presentation typically occurred between 12 and 36 months. In addition to the two disease-defining features of spondylometaphyseal dysplasia and hypomyelination on MRI, common clinical signs and symptoms included motor deterioration, spasticity, tremor, ataxia, dysarthria, cognitive defects, pulmonary hypertension, nystagmus, and vision loss due to retinopathy. The course of the disease was slowly progressive. All patients had maternally inherited or de novo mutations in or near exon 7 of AIFM1, within a region of 70 bp, including synonymous and intronic changes. AIFM1 mutations have previously been associated with neurologic presentations as varied as intellectual disability, hearing loss, neuropathy, and striatal necrosis, while AIFM1 mutations in this small region present with a distinct phenotype implicating bone. Analysis of cell lines derived from four patients identified significant reductions in AIFM1 mRNA and protein levels in osteoblasts. We hypothesize that AIFM1 functions in bone metabolism and myelination and is responsible for the unique phenotype in this condition.

RUNX1 contributes to higher-order chromatin organization and gene regulation in breast cancer cells.

RUNX1 is a transcription factor functioning both as an oncogene and a tumor suppressor in breast cancer. RUNX1 alters chromatin structure in cooperation with chromatin modifier and remodeling enzymes. In this study, we examined the relationship between RUNX1-mediated transcription and genome organization. We characterized genome-wide RUNX1 localization and performed RNA-seq and Hi-C in RUNX1-depleted and control MCF-7 breast cancer cells. RNA-seq analysis showed that RUNX1 depletion led to up-regulation of genes associated with chromatin structure and down-regulation of genes related to extracellular matrix biology, as well as NEAT1 and MALAT1 lncRNAs. Our ChIP-Seq analysis supports a prominent role for RUNX1 in transcriptional activation. About 30% of all RUNX1 binding sites were intergenic, indicating diverse roles in promoter and enhancer regulation and suggesting additional functions for RUNX1. Hi-C analysis of RUNX1-depleted cells demonstrated that overall three-dimensional genome organization is largely intact, but indicated enhanced association of RUNX1 near Topologically Associating Domain (TAD) boundaries and alterations in long-range interactions. These results suggest an architectural role for RUNX1 in fine-tuning local interactions rather than in global organization. Our results provide novel insight into RUNX1-mediated perturbations of higher-order genome organization that are functionally linked with RUNX1-dependent compromised gene expression in breast cancer cells.

Structural organization of the inactive X chromosome in the mouse.

X-chromosome inactivation (XCI) involves major reorganization of the X chromosome as it becomes silent and heterochromatic. During female mammalian development, XCI is triggered by upregulation of the non-coding Xist RNA from one of the two X chromosomes. Xist coats the chromosome in cis and induces silencing of almost all genes via its A-repeat region, although some genes (constitutive escapees) avoid silencing in most cell types, and others (facultative escapees) escape XCI only in specific contexts. A role for Xist in organizing the inactive X (Xi) chromosome has been proposed. Recent chromosome conformation capture approaches have revealed global loss of local structure on the Xi chromosome and formation of large mega-domains, separated by a region containing the DXZ4 macrosatellite. However, the molecular architecture of the Xi chromosome, in both the silent and expressed regions,remains unclear. Here we investigate the structure, chromatin accessibility and expression status of the mouse Xi chromosome in highly polymorphic clonal neural progenitors (NPCs) and embryonic stem cells. We demonstrate a crucial role for Xist and the DXZ4-containing boundary in shaping Xi chromosome structure using allele-specific genome-wide chromosome conformation capture (Hi-C) analysis, an assay for transposase-accessible chromatin with high throughput sequencing (ATAC-seq) and RNA sequencing. Deletion of the boundary disrupts mega-domain formation, and induction of Xist RNA initiates formation of the boundary and the loss of DNA accessibility. We also show that in NPCs, the Xi chromosome lacks active/inactive compartments and topologically associating domains (TADs), except around genes that escape XCI. Escapee gene clusters display TAD-like structures and retain DNA accessibility at promoter-proximal and CTCF-binding sites. Furthermore, altered patterns of facultative escape genes indifferent neural progenitor clones are associated with the presence of different TAD-like structures after XCI. These findings suggest a key role for transcription and CTCF in the formation of TADs in the context of the Xi chromosome in neural progenitors.

SMARCA4 regulates gene expression and higher-order chromatin structure in proliferating mammary epithelial cells.

The packaging of DNA into chromatin plays an important role in transcriptional regulation and nuclear processes. Brahma-related gene-1 SMARCA4 (also known as BRG1), the essential ATPase subunit of the mammalian SWI/SNF chromatin remodeling complex, uses the energy from ATP hydrolysis to disrupt nucleosomes at target regions. Although the transcriptional role of SMARCA4 at gene promoters is well-studied, less is known about its role in higher-order genome organization. SMARCA4 knockdown in human mammary epithelial MCF-10A cells resulted in 176 up-regulated genes, including many related to lipid and calcium metabolism, and 1292 down-regulated genes, some of which encode extracellular matrix (ECM) components that can exert mechanical forces and affect nuclear structure. ChIP-seq analysis of SMARCA4 localization and SMARCA4-bound super-enhancers demonstrated extensive binding at intergenic regions. Furthermore, Hi-C analysis showed extensive SMARCA4-mediated alterations in higher-order genome organization at multiple resolutions. First, SMARCA4 knockdown resulted in clustering of intra- and inter-subtelomeric regions, demonstrating a novel role for SMARCA4 in telomere organization. SMARCA4 binding was enriched at topologically associating domain (TAD) boundaries, and SMARCA4 knockdown resulted in weakening of TAD boundary strength. Taken together, these findings provide a dynamic view of SMARCA4-dependent changes in higher-order chromatin organization and gene expression, identifying SMARCA4 as a novel component of chromatin organization.

Genetics and Genomics of Longitudinal Lung Function Patterns in Individuals with Asthma.

RATIONALE: Patterns of longitudinal lung function growth and decline in childhood asthma have been shown to be important in determining risk for future respiratory ailments including chronic airway obstruction and chronic obstructive pulmonary disease. OBJECTIVES: To determine the genetic underpinnings of lung function patterns in subjects with childhood asthma. METHODS: We performed a genome-wide association study of 581 non-Hispanic white individuals with asthma that were previously classified by patterns of lung function growth and decline (normal growth, normal growth with early decline, reduced growth, and reduced growth with early decline). The strongest association was also measured in two additional cohorts: a small asthma cohort and a large chronic obstructive pulmonary disease metaanalysis cohort. Interaction between the genomic region encompassing the most strongly associated single-nucleotide polymorphism and nearby genes was assessed by two chromosome conformation capture assays. MEASUREMENTS AND MAIN RESULTS: An intergenic single-nucleotide polymorphism (rs4445257) on chromosome 8 was strongly associated with the normal growth with early decline pattern compared with all other pattern groups (P = 6.7 × 10-9; odds ratio, 2.8; 95% confidence interval, 2.0-4.0); replication analysis suggested this variant had opposite effects in normal growth with early decline and reduced growth with early decline pattern groups. Chromosome conformation capture experiments indicated a chromatin interaction between rs4445257 and the promoter of the distal CSMD3 gene. CONCLUSIONS: Early decline in lung function after normal growth is associated with a genetic polymorphism that may also protect against early decline in reduced growth groups. Clinical trial registered with www.clinicaltrials.gov (NCT00000575).

Activation of proto-oncogenes by disruption of chromosome neighborhoods.

Oncogenes are activated through well-known chromosomal alterations such as gene fusion, translocation, and focal amplification. In light of recent evidence that the control of key genes depends on chromosome structures called insulated neighborhoods, we investigated whether proto-oncogenes occur within these structures and whether oncogene activation can occur via disruption of insulated neighborhood boundaries in cancer cells. We mapped insulated neighborhoods in T cell acute lymphoblastic leukemia (T-ALL) and found that tumor cell genomes contain recurrent microdeletions that eliminate the boundary sites of insulated neighborhoods containing prominent T-ALL proto-oncogenes. Perturbation of such boundaries in nonmalignant cells was sufficient to activate proto-oncogenes. Mutations affecting chromosome neighborhood boundaries were found in many types of cancer. Thus, oncogene activation can occur via genetic alterations that disrupt insulated neighborhoods in malignant cells.