DNA mismatch and damage patterns revealed by single-molecule sequencing.

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

A common flanking variant is associated with enhanced stability of the FGF14-SCA27B repeat locus.

The factors driving or preventing pathological expansion of tandem repeats remain largely unknown. Here, we assessed the FGF14 (GAA)·(TTC) repeat locus in 2,530 individuals by long-read and Sanger sequencing and identified a common 5'-flanking variant in 70.34% of alleles analyzed (3,463/4,923) that represents the phylogenetically ancestral allele and is present on all major haplotypes. This common sequence variation is present nearly exclusively on nonpathogenic alleles with fewer than 30 GAA-pure triplets and is associated with enhanced stability of the repeat locus upon intergenerational transmission and increased Fiber-seq chromatin accessibility.

Structural polymorphism and diversity of human segmental duplications.

Segmental duplications (SDs) contribute significantly to human disease, evolution, and diversity yet have been difficult to resolve at the sequence level. We present a population genetics survey of SDs by analyzing 170 human genome assemblies where the majority of SDs are fully resolved using long-read sequence assembly. Excluding the acrocentric short arms, we identify 173.2 Mbp of duplicated sequence (47.4 Mbp not present in the telomere-to-telomere reference) distinguishing fixed from structurally polymorphic events. We find that intrachromosomal SDs are among the most variable with rare events mapping near their progenitor sequences. African genomes harbor significantly more intrachromosomal SDs and are more likely to have recently duplicated gene families with higher copy number when compared to non-African samples. A comparison to a resource of 563 million full-length Iso-Seq reads identifies 201 novel, potentially protein-coding genes corresponding to these copy number polymorphic SDs.

Addressing dispersion in mis-measured multivariate binomial outcomes: A novel statistical approach for detecting differentially methylated regions in bisulfite sequencing data.

Motivated by a DNA methylation application, this article addresses the problem of fitting and inferring a multivariate binomial regression model for outcomes that are contaminated by errors and exhibit extra-parametric variations, also known as dispersion. While dispersion in univariate binomial regression has been extensively studied, addressing dispersion in the context of multivariate outcomes remains a complex and relatively unexplored task. The complexity arises from a noteworthy data characteristic observed in our motivating dataset: non-constant yet correlated dispersion across outcomes. To address this challenge and account for possible measurement error, we propose a novel hierarchical quasi-binomial varying coefficient mixed model, which enables flexible dispersion patterns through a combination of additive and multiplicative dispersion components. To maximize the Laplace-approximated quasi-likelihood of our model, we further develop a specialized two-stage expectation-maximization (EM) algorithm, where a plug-in estimate for the multiplicative scale parameter enhances the speed and stability of the EM iterations. Simulations demonstrated that our approach yields accurate inference for smooth covariate effects and exhibits excellent power in detecting non-zero effects. Additionally, we applied our proposed method to investigate the association between DNA methylation, measured across the genome through targeted custom capture sequencing of whole blood, and levels of anti-citrullinated protein antibodies (ACPA), a preclinical marker for rheumatoid arthritis (RA) risk. Our analysis revealed 23 significant genes that potentially contribute to ACPA-related differential methylation, highlighting the relevance of cell signaling and collagen metabolism in RA. We implemented our method in the R Bioconductor package called "SOMNiBUS."

Evaluation of Hi-C sequencing for the detection of gene fusions in hematologic and solid pediatric cancer samples.

HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome and has shown promise in detecting genomic rearrangements in translational research studies. To evaluate HiC as a potential clinical diagnostic platform, analytical concordance with routine laboratory testing was assessed using primary pediatric leukemia and sarcoma specimens previously positive for clinically significant genomic rearrangements. Archived specimen types tested included viable and nonviable frozen leukemic cells, as well as formalin-fixed paraffin-embedded (FFPE) tumor tissues. Initially, pediatric acute myeloid leukemia (AML) and alveolar rhabdomyosarcoma (A-RMS) specimens with known genomic rearrangements were subjected to HiC analysis to assess analytical concordance. Subsequently, a discovery cohort consisting of AML and acute lymphoblastic leukemia (ALL) cases with no known genomic rearrangements based on prior clinical diagnostic testing were evaluated to determine whether HiC could detect rearrangements. Using a standard sequencing depth of 50 million raw read-pairs per sample, or approximately 5X raw genomic coverage, 100% concordance was observed between HiC and previous clinical cytogenetic and molecular testing. In the discovery cohort, a clinically relevant gene fusion was detected in 45% of leukemia cases (5/11). This study demonstrates the value of HiC sequencing to medical diagnostic testing as it identified several clinically significant rearrangements, including those that might have been missed by current clinical testing workflows. Key points: HiC sequencing is a DNA-based next-generation sequencing method that preserves the 3D conformation of the genome, facilitating detection of genomic rearrangements.HiC was 100% concordant with clinical diagnostic testing workflows for detecting clinically significant genomic rearrangements in pediatric leukemia and rhabdomyosarcoma specimens.HiC detected clinically significant genomic rearrangements not previously detected by prior clinical cytogenetic and molecular testing.HiC performed well with archived non-viable and viable frozen leukemic cell samples, as well as archived formalin-fixed paraffin-embedded tumor tissue specimens.

Large-scale mutational analysis identifies UNC93B1 variants that drive TLR-mediated autoimmunity in mice and humans.

Nucleic acid-sensing Toll-like receptors (TLR) 3, 7/8, and 9 are key innate immune sensors whose activities must be tightly regulated to prevent systemic autoimmune or autoinflammatory disease or virus-associated immunopathology. Here, we report a systematic scanning-alanine mutagenesis screen of all cytosolic and luminal residues of the TLR chaperone protein UNC93B1, which identified both negative and positive regulatory regions affecting TLR3, TLR7, and TLR9 responses. We subsequently identified two families harboring heterozygous coding mutations in UNC93B1, UNC93B1+/T93I and UNC93B1+/R336C, both in key negative regulatory regions identified in our screen. These patients presented with cutaneous tumid lupus and juvenile idiopathic arthritis plus neuroinflammatory disease, respectively. Disruption of UNC93B1-mediated regulation by these mutations led to enhanced TLR7/8 responses, and both variants resulted in systemic autoimmune or inflammatory disease when introduced into mice via genome editing. Altogether, our results implicate the UNC93B1-TLR7/8 axis in human monogenic autoimmune diseases and provide a functional resource to assess the impact of yet-to-be-reported UNC93B1 mutations.

Independent expansion, selection and hypervariability of the TBC1D3 gene family in humans.

TBC1D3 is a primate-specific gene family that has expanded in the human lineage and has been implicated in neuronal progenitor proliferation and expansion of the frontal cortex. The gene family and its expression have been challenging to investigate because it is embedded in high-identity and highly variable segmental duplications. We sequenced and assembled the gene family using long-read sequencing data from 34 humans and 11 nonhuman primate species. Our analysis shows that this particular gene family has independently duplicated in at least five primate lineages, and the duplicated loci are enriched at sites of large-scale chromosomal rearrangements on chromosome 17. We find that most humans vary along two TBC1D3 clusters where human haplotypes are highly variable in copy number, differing by as many as 20 copies, and structure (structural heterozygosity 90%). We also show evidence of positive selection, as well as a significant change in the predicted human TBC1D3 protein sequence. Lastly, we find that, despite multiple duplications, human TBC1D3 expression is limited to a subset of copies and, most notably, from a single paralog group: TBC1D3-CDKL. These observations may help explain why a gene potentially important in cortical development can be so variable in the human population.

Genomic Answers for Kids: Toward more equitable access to genomic testing for rare diseases in rural populations.

Next-generation sequencing has revolutionized the speed of rare disease (RD) diagnoses. While clinical exome and genome sequencing represent an effective tool for many RD diagnoses, there is room to further improve the diagnostic odyssey of many RD patients. One recognizable intervention lies in increasing equitable access to genomic testing. Rural communities represent a significant portion of underserved and underrepresented individuals facing additional barriers to diagnosis and treatment. Primary care providers (PCPs) at local clinics, though sometimes suspicious of a potential benefit of genetic testing for their patients, have significant constraints in pursuing it themselves and rely on referrals to specialists. Yet, these referrals are typically followed by long waitlists and significant delays in clinical assessment, insurance clearance, testing, and initiation of diagnosis-informed care management. Not only is this process time intensive, but it also often requires multiple visits to urban medical centers for which distance may be a significant barrier to rural families. Therefore, providing early, "direct-to-provider" (DTP) local access to unrestrictive genomic testing is likely to help speed up diagnostic times and access to care for RD patients in rural communities. In a pilot study with a PCP clinic in rural Kansas, we observed a minimum 5.5 months shortening of time to diagnosis through the DTP exome sequencing program as compared to rural patients receiving genetic testing through the "traditional" PCP-referral-to-specialist scheme. We share our experience to encourage future partnerships beyond our center. Our efforts represent just one step in fostering greater diversity and equity in genomic studies.

Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores.

Recent studies have revealed the pervasive landscape of rare structural variants (rSVs) present in human genomes. rSVs can have extreme effects on the expression of proximal genes and, in a rare disease context, have been implicated in patient cases where no diagnostic single nucleotide variant (SNV) was found. Approaches for integrating rSVs to date have focused on targeted approaches in known Mendelian rare disease genes. This approach is intractable for rare diseases with many causal loci or patients with complex, multi-phenotype syndromes. We hypothesized that integrating trait-relevant polygenic scores (PGS) would provide a substantial reduction in the number of candidate disease genes in which to assess rSV effects. We further implemented a method for ranking PGS genes to define a set of core/key genes where a rSV has the potential to exert relatively larger effects on disease risk. Among a subset of patients enrolled in the Genomic Answers for Kids (GA4K) rare disease program (N=497), we used PacBio HiFi long-read whole genome sequencing (lrWGS) to identify rSVs intersecting genes in trait-relevant PGSs. Illustrating our approach in Autism (N=54 cases), we identified 22, 019 deletions, 2,041 duplications, 87,826 insertions, and 214 inversions overlapping putative core/key PGS genes. Additionally, by integrating genomic constraint annotations from gnomAD, we observed that rare duplications overlapping putative core/key PGS genes were frequently in higher constraint regions compared to controls (P = 1×10-03). This difference was not observed in the lowest-ranked gene set (P = 0.15). Overall, our study provides a framework for the annotation of long-read rSVs from lrWGS data and prioritization of disease-linked genomic regions for downstream functional validation of rSV impacts. To enable reuse by other researchers, we have made SV allele frequencies and gene associations freely available.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads.

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly.

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

EpiVar Browser: advanced exploration of epigenomics data under controlled access.

MOTIVATION: Human epigenomic data has been generated by large consortia for thousands of cell types to be used as a reference map of normal and disease chromatin states. Since epigenetic data contains potentially identifiable information, similarly to genetic data, most raw files generated by these consortia are stored in controlled-access databases. It is important to protect identifiable information, but this should not hinder secure sharing of these valuable datasets. RESULTS: Guided by the Framework for responsible sharing of genomic and health-related data from the Global Alliance for Genomics and Health (GA4GH), we have developed an approach and a tool to facilitate the exploration of epigenomics datasets' aggregate results, while filtering out identifiable information. Specifically, the EpiVar Browser allows a user to navigate an epigenetic dataset from a cohort of individuals and enables direct exploration of genotype-chromatin phenotype relationships. Because individual genotypes and epigenetic signal tracks are not directly accessible, and rather aggregated in the portal output, no identifiable data is released, yet the interface allows for dynamic genotype-epigenome interrogation. This approach has the potential to accelerate analyses that would otherwise require a lengthy multi-step approval process and provides a generalizable strategy to facilitate responsible access to sensitive epigenomics data. AVAILABILITY AND IMPLEMENTATION: Online portal: https://computationalgenomics.ca/tools/epivar; EpiVar Browser source code: https://github.com/c3g/epivar-browser; bw-merge-window tool source code: https://github.com/c3g/bw-merge-window.

Epigenetic variation impacts individual differences in the transcriptional response to influenza infection.

Humans display remarkable interindividual variation in their immune response to identical challenges. Yet, our understanding of the genetic and epigenetic factors contributing to such variation remains limited. Here we performed in-depth genetic, epigenetic and transcriptional profiling on primary macrophages derived from individuals of European and African ancestry before and after infection with influenza A virus. We show that baseline epigenetic profiles are strongly predictive of the transcriptional response to influenza A virus across individuals. Quantitative trait locus (QTL) mapping revealed highly coordinated genetic effects on gene regulation, with many cis-acting genetic variants impacting concomitantly gene expression and multiple epigenetic marks. These data reveal that ancestry-associated differences in the epigenetic landscape can be genetically controlled, even more than gene expression. Lastly, among QTL variants that colocalized with immune-disease loci, only 7% were gene expression QTL, while the remaining genetic variants impact epigenetic marks, stressing the importance of considering molecular phenotypes beyond gene expression in disease-focused studies.

Pangenome graphs improve the analysis of structural variants in rare genetic diseases.

Rare DNA alterations that cause heritable diseases are only partially resolvable by clinical next-generation sequencing due to the difficulty of detecting structural variation (SV) in all genomic contexts. Long-read, high fidelity genome sequencing (HiFi-GS) detects SVs with increased sensitivity and enables assembling personal and graph genomes. We leverage standard reference genomes, public assemblies (n = 94) and a large collection of HiFi-GS data from a rare disease program (Genomic Answers for Kids, GA4K, n = 574 assemblies) to build a graph genome representing a unified SV callset in GA4K, identify common variation and prioritize SVs that are more likely to cause genetic disease (MAF < 0.01). Using graphs, we obtain a higher level of reproducibility than the standard reference approach. We observe over 200,000 SV alleles unique to GA4K, including nearly 1000 rare variants that impact coding sequence. With improved specificity for rare SVs, we isolate 30 candidate SVs in phenotypically prioritized genes, including known disease SVs. We isolate a novel diagnostic SV in KMT2E, demonstrating use of personal assemblies coupled with pangenome graphs for rare disease genomics. The community may interrogate our pangenome with additional assemblies to discover new SVs within the allele frequency spectrum relevant to genetic diseases.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.

A syndromic neurodevelopmental disorder caused by rare variants in PPFIA3.

PPFIA3 encodes the protein-tyrosine phosphatase, receptor-type, F-polypeptide-interacting-protein-alpha-3 (PPFIA3), which is a member of the LAR-protein-tyrosine phosphatase-interacting-protein (liprin) family involved in synapse formation and function, synaptic vesicle transport, and presynaptic active zone assembly. The protein structure and function are evolutionarily well conserved, but human diseases related to PPFIA3 dysfunction are not yet reported in OMIM. Here, we report 20 individuals with rare PPFIA3 variants (19 heterozygous and 1 compound heterozygous) presenting with developmental delay, intellectual disability, hypotonia, dysmorphisms, microcephaly or macrocephaly, autistic features, and epilepsy with reduced penetrance. Seventeen unique PPFIA3 variants were detected in 18 families. To determine the pathogenicity of PPFIA3 variants in vivo, we generated transgenic fruit flies producing either human wild-type (WT) PPFIA3 or five missense variants using GAL4-UAS targeted gene expression systems. In the fly overexpression assays, we found that the PPFIA3 variants in the region encoding the N-terminal coiled-coil domain exhibited stronger phenotypes compared to those affecting the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin-α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 function is partially conserved in the fly. However, two of the tested variants failed to rescue the lethality at the larval stage and one variant failed to rescue lethality at the adult stage. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant-negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

Characterization and visualization of tandem repeats at genome scale.

Tandem repeat (TR) variation is associated with gene expression changes and numerous rare monogenic diseases. Although long-read sequencing provides accurate full-length sequences and methylation of TRs, there is still a need for computational methods to profile TRs across the genome. Here we introduce the Tandem Repeat Genotyping Tool (TRGT) and an accompanying TR database. TRGT determines the consensus sequences and methylation levels of specified TRs from PacBio HiFi sequencing data. It also reports reads that support each repeat allele. These reads can be subsequently visualized with a companion TR visualization tool. Assessing 937,122 TRs, TRGT showed a Mendelian concordance of 98.38%, allowing a single repeat unit difference. In six samples with known repeat expansions, TRGT detected all expansions while also identifying methylation signals and mosaicism and providing finer repeat length resolution than existing methods. Additionally, we released a database with allele sequences and methylation levels for 937,122 TRs across 100 genomes.

Comparing Attitudes About Genomic Privacy and Data Sharing in Adolescents and Parents of Children Enrolled in a Genomic Research Repository.

BACKGROUND: Sharing of genomic data aims to make efficient use of limited resources, which may be particularly valuable in rare disease research. Adult research participants and parents of pediatric research participants have shown support for data sharing with protections, but little is known about adolescent attitudes on genomic privacy and data sharing. METHODS: In-depth interviews were conducted with 10 adolescents and 18 parents of children enrolled in a pediatric genomic research repository. Interview transcripts were analyzed for themes on attitudes toward genomic privacy, restricted-access data sharing, and open-access data sharing. Findings in adolescent and parent participants were compared and contrasted. RESULTS: No adolescents endorsed privacy concerns for restricted-access data sharing. Both adolescents and parents saw value in data sharing for reaching the goals of research and discussed trust in institutions and researchers to protect their data and use it as intended. Adolescents were more likely than parents to accept open-access data sharing, including after risks were discussed. CONCLUSIONS: In this exploratory study, adolescents and parents enrolled in a genomic research repository shared many attitudes about genomic data sharing, but adolescents were less concerned about privacy and more agreeable toward open-access data sharing. Future research is needed to investigate this hypothesis in expanded populations and settings, and to clarify whether adolescent attitudes change with age and experiences.