Long-Read Sequencing Identifies the First Retrotransposon Insertion and Resolves Structural Variants Causing Antithrombin Deficiency

Next-generation sequencing strategies in venous thromboembolism: in whom and for what purpose?

This invited review follows the oral presentation "To Sequence or Not to Sequence, That Is Not the Question; But 'When, Who, Which and What For?' Is" given during the State of the Art session "Translational Genomics in Thrombosis: From OMICs to Clinics" of the International Society on Thrombosis and Haemostasis 2023 Congress. Emphasizing the power of next-generation sequencing technologies and the diverse strategies associated with DNA variant analysis, this review highlights the unresolved questions and challenges in their implementation both for the clinical diagnosis of venous thromboembolism and in translational research.

Long-range and real-time PCR identification of a large SERPINC1 deletion in a patient with antithrombin deficiency

Congenital antithrombin (AT) or serpin C1 deficiency, caused by a SERPINC1 abnormality, is a high-risk factor for venous thrombosis. SERPINC1 is prone to genetic rearrangement, because it contains numerous Alu elements. In this study, a Japanese patient who developed deep vein thrombosis during pregnancy and exhibited low AT activity underwent SERPINC1 gene analysis using routine methods: long-range polymerase chain reaction (PCR) and real-time PCR. Sequencing using long-range PCR products revealed no pathological variants in SERPINC1 exons or exon-intron junctions, and all the identified variants were homozygous, suggesting a deletion in one SERPINC1 allele. Copy number quantification for each SERPINC1 exon using real-time PCR revealed half the number of exon 1 and 2 copies compared with controls. Moreover, a deletion region was deduced by quantifying the 5'-upstream region copy number of SERPINC1 for each constant region. Direct long-range PCR sequencing with primers for the 5'-end of each presumed deletion region revealed a large Alu-mediated deletion (∼13 kb) involving SERPINC1 exons 1 and 2. Thus, a large deletion was identified in SERPINC1 using conventional PCR methods.

Targeted long-read sequencing identifies and characterizes structural variants in cases of inherited platelet disorders

BACKGROUND

Genetic diagnosis of inherited platelet disorders (IPDs) is mainly performed by high-throughput sequencing (HTS). These short-read-based sequencing methods sometimes fail to characterize the genetics of the disease.

OBJECTIVES

To evaluate nanopore long-read DNA sequencing for characterization of structural variants (SVs) in patients with IPDs.

METHODS

Four patients with a clinical and laboratory diagnosis of Glanzmann thrombasthenia (GT) (P1 and P2) and Hermansky-Pudlak syndrome (HPS) (P3 and P4) in whom HTS missed the underlying molecular cause were included. DNA was analyzed by both standard HTS and nanopore sequencing on a MinION device (Oxford Nanopore Technologies) after enrichment of DNA spanning regions covering GT and HPS genes.

RESULTS

In patients with GT, HTS identified only 1 heterozygous ITGB3 splice variant c.2301+1G>C in P2. In patients with HPS, a homozygous deletion in HPS5 was suspected in P3, and 2 heterozygous HPS3 variants, c.2464C>T (p.Arg822∗) and a deletion affecting 2 exons, were reported in P4. Nanopore sequencing revealed a complex SV affecting exons 2 to 6 in ITGB3 (deletion-inversion-duplication) in homozygosity in P1 and compound heterozygosity with the splice variant in P2. In the 2 patients with HPS, nanopore defined the length of the SVs, which were characterized at nucleotide resolution. This allowed the identification of repetitive Alu elements at the breakpoints and the design of specific polymerase chain reactions for family screening.

CONCLUSION

The nanopore technology overcomes the limitations of standard short-read sequencing techniques in SV characterization. Using nanopore, we characterized novel defects in ITGB3, HPS5, and HPS3, highlighting the utility of long-read sequencing as an additional diagnostic tool in IPDs.

The landscape of human SVA retrotransposons

SINE-VNTR-Alu (SVA) retrotransposons are evolutionarily young and still-active transposable elements (TEs) in the human genome. Several pathogenic SVA insertions have been identified that directly mutate host genes to cause neurodegenerative and other types of diseases. However, due to their sequence heterogeneity and complex structures as well as limitations in sequencing techniques and analysis, SVA insertions have been less well studied compared to other mobile element insertions. Here, we identified polymorphic SVA insertions from 3646 whole-genome sequencing (WGS) samples of >150 diverse populations and constructed a polymorphic SVA insertion reference catalog. Using 20 long-read samples, we also assembled reference and polymorphic SVA sequences and characterized the internal hexamer/variable-number-tandem-repeat (VNTR) expansions as well as differing SVA activity for SVA subfamilies and human populations. In addition, we developed a module to annotate both reference and polymorphic SVA copies. By characterizing the landscape of both reference and polymorphic SVA retrotransposons, our study enables more accurate genotyping of these elements and facilitate the discovery of pathogenic SVA insertions.

Detection and annotation of transposable element insertions and deletions on the human genome using nanopore sequencing

Repetitive sequences represent about 45% of the human genome. Some are transposable elements (TEs) with the ability to change their position in the genome, creating genetic variability both as insertions or deletions, with potential pathogenic consequences. We used long-read nanopore sequencing to identify TE variants in the genomes of 24 patients with antithrombin deficiency. We identified 7 344 TE insertions and 3 056 TE deletions, 2 926 were not previously described in publicly available databases. The insertions affected 3 955 genes, with 6 insertions located in exons, 3 929 in introns, and 147 in promoters. Potential functional impact was evaluated with gene annotation and enrichment analysis, which suggested a strong relationship with neuron-related functions and autism. We conclude that this study encourages the generation of a complete map of TEs in the human genome, which will be useful for identifying new TEs involved in genetic disorders.

Genome sequencing and comprehensive rare-variant analysis of 465 families with neurodevelopmental disorders

Despite significant progress in unraveling the genetic causes of neurodevelopmental disorders (NDDs), a substantial proportion of individuals with NDDs remain without a genetic diagnosis after microarray and/or exome sequencing. Here, we aimed to assess the power of short-read genome sequencing (GS), complemented with long-read GS, to identify causal variants in participants with NDD from the National Institute for Health and Care Research (NIHR) BioResource project. Short-read GS was conducted on 692 individuals (489 affected and 203 unaffected relatives) from 465 families. Additionally, long-read GS was performed on five affected individuals who had structural variants (SVs) in technically challenging regions, had complex SVs, or required distal variant phasing. Causal variants were identified in 36% of affected individuals (177/489), and a further 23% (112/489) had a variant of uncertain significance after multiple rounds of re-analysis. Among all reported variants, 88% (333/380) were coding nuclear SNVs or insertions and deletions (indels), and the remainder were SVs, non-coding variants, and mitochondrial variants. Furthermore, long-read GS facilitated the resolution of challenging SVs and invalidated variants of difficult interpretation from short-read GS. This study demonstrates the value of short-read GS, complemented with long-read GS, in investigating the genetic causes of NDDs. GS provides a comprehensive and unbiased method of identifying all types of variants throughout the nuclear and mitochondrial genomes in individuals with NDD.

A Genome-Wide Screen for the Exonisation of Reference SINE-VNTR-Alus and Their Expression in CNS Tissues of Individuals with Amyotrophic Lateral Sclerosis

The hominid-specific retrotransposon SINE-VNTR-Alu (SVA) is a composite element that has contributed to the genetic variation between individuals and influenced genomic structure and function. SVAs are involved in modulating gene expression and splicing patterns, altering mRNA levels and sequences, and have been associated with the development of disease. We evaluated the genome-wide effects of SVAs present in the reference genome on transcript sequence and expression in the CNS of individuals with and without the neurodegenerative disorder Amyotrophic Lateral Sclerosis (ALS). This study identified SVAs in the exons of 179 known transcripts, several of which were expressed in a tissue-specific manner, as well as 92 novel exonisation events occurring in the motor cortex. An analysis of 65 reference genome SVAs polymorphic for their presence/absence in the ALS consortium cohort did not identify any elements that were significantly associated with disease status, age at onset, and survival. However, there were transcripts, such as transferrin and HLA-A, that were differentially expressed between those with or without disease, and expression levels were associated with the genotype of proximal SVAs. This study demonstrates the functional consequences of several SVA elements altering mRNA splicing patterns and expression levels in tissues of the CNS.

Impact of genetic structural variants in factor XI deficiency: identification, accurate characterization, and inferred mechanism by long-read sequencing

BACKGROUND

Congenital factor XI (FXI) deficiency is a probably underestimated coagulopathy that confers antithrombotic protection. Characterization of genetic defects in F11 is mainly focused on the identification of single-nucleotide variants and small insertion/deletions because they represent up to 99% of the alterations accounting for factor deficiency, with only 3 gross gene defects of structural variants (SVs) having been described.

OBJECTIVES

To identify and characterize the SVs affecting F11.

METHODS

The study was performed in 93 unrelated subjects with FXI deficiency recruited in Spanish hospitals over a period of 25 years (1997-2022). F11 was analyzed by next-generation sequencing, multiplex ligand probe amplification, and long-read sequencing.

RESULTS

Our study identified 30 different genetic variants. Interestingly, we found 3 SVs, all heterozygous: a complex duplication affecting exons 8 and 9, a tandem duplication of exon 14, and a large deletion affecting the whole gene. Nucleotide resolution obtained by long-read sequencing revealed Alu repetitive elements involved in all breakpoints. The large deletion was probably generated de novo in the paternal allele during gametogenesis, and despite affecting 30 additional genes, no syndromic features were described.

CONCLUSION

SVs may account for a high proportion of F11 genetic defects implicated in the molecular pathology of congenital FXI deficiency. These SVs, likely caused by a nonallelic homologous recombination involving repetitive elements, are heterogeneous in both type and length and may be de novo. These data support the inclusion of methods to detect SVs in this disorder, with long-read-based methods being the most appropriate because they detect all SVs and achieve adequate nucleotide resolution.

Usefulness and Limitations of Multiple Ligation-Dependent Probe Amplification in Antithrombin Deficiency

Multiplex ligation-dependent probe amplification (MLPA) identifies genetic structural variants in SERPINC1 in 5% of cases with antithrombin deficiency (ATD), the most severe congenital thrombophilia. Our aim was to unravel the utility and limitations of MLPA in a large cohort of unrelated patients with ATD (N = 341). MLPA identified 22 structural variants (SVs) causing ATD (6.5%). MLPA did not detect SVs affecting introns (four cases), and the diagnosis was inaccurate in two cases according to long-range PCR or nanopore sequencing. MLPA was used to detect possible hidden SVs in 61 cases with type I deficiency with single nucleotide variations (SNVs) or small insertion/deletion (INDEL). One case had a false deletion of exon 7, as the 29-bp deletion affected an MLPA probe. We evaluated 32 variants affecting MLPA probes: 27 SNVs and 5 small INDELs. In three cases, MLPA gave false-positive results, all diagnosed as deletions of the affected exon: a small INDEL complex, and two SNVs affecting MLPA probes. Our study confirms the utility of MLPA to detect SVs in ATD, but also shows some limitations in detecting intronic SVs. MLPA renders imprecise and false-positive results for genetic defects which affect MLPA probes. Our results encourage the validation of MLPA results.

Novel Alu-mediated deletions of the SMN1 gene were identified by ultra-long read sequencing technology in patients with spinal muscular atrophy

Spinal muscular atrophy (SMA) is a neuromuscular disease caused by biallelic variants of the survival motor neuron 1 (SMN1) gene. In this study, our aim was to make a molecular diagnosis in two patients with SMA carrying only one SMN1 copy number. Using ultra-long read sequencing (Ultra-LRS), 1415 bp deletion and 3348 bp deletion of the SMN1 gene were identified in patient 1 and the father of patient 2, respectively. Ultra-LRS revealed two novel deletions, starting from the SMN1 promoter to intron 1. It also accurately provided the location of the deletion breakpoints in the SMN1 gene: chr5 g.70,924,798-70,926,212 for a 1415 bp deletion; chr5 g.70,922,695-70,926,042 for a 3348 bp deletion. By analyzing the breakpoint junctions, we identified that these genomic sequences were composed of Alu sequences, including AluJb, AluYm1, AluSq, and AluYm1, indicating that Alu-mediated rearrangements are a mechanism of SMN1 deletion events. In addition, full-length SMN1 transcripts and SMN protein in patient 1 were significantly decreased (p < 0.01), suggesting that a 1415 bp deletion that included the transcription and translation initiation sites of the SMN1 gene had severe consequences for SMN expression. Ultra-LRS can easily distinguish highly homozygous genes compared to other detection technologies, which is useful for detecting SMN1 intragenic mutations, to quickly discover structural rearrangements and to precisely present the breakpoint positions.

Improvement of large copy number variant detection by whole genome nanopore sequencing

INTRODUCTION

Whole-genome sequencing using nanopore technologies can uncover structural variants, which are DNA rearrangements larger than 50 base pairs. Nanopore technologies can also characterize their boundaries with single-base accuracy, owing to the kilobase-long reads that encompass either full variants or their junctions. Other methods, such as next-generation short read sequencing or PCR assays, are limited in their capabilities to detect or characterize structural variants. However, the existing software for nanopore sequencing data analysis still reports incomplete variant sets, which also contain erroneous calls, a considerable obstacle for the molecular diagnosis or accurate genotyping of populations.

METHODS

We compared multiple factors affecting variant calling, such as reference genome version, aligner (minimap2, NGMLR, and lra) choice, and variant caller combinations (Sniffles, CuteSV, SVIM, and NanoVar), to find the optimal group of tools for calling large (>50 kb) deletions and duplications, using data from seven patients exhibiting gross gene defects on SERPINC1 and from a reference variant set as the control. The goal was to obtain the most complete, yet reasonably specific group of large variants using a single cell of PromethION sequencing, which yielded lower depth coverage than short-read sequencing. We also used a custom method for the statistical analysis of the coverage value to refine the resulting datasets.

RESULTS

We found that for large deletions and duplications (>50 kb), the existing software performed worse than for smaller ones, in terms of both sensitivity and specificity, and newer tools had not improved this. Our novel software, disCoverage, could polish variant callers' results, improving specificity by up to 62% and sensitivity by 15%, the latter requiring other data or samples.

CONCLUSION

We analyzed the current situation of >50-kb copy number variants with nanopore sequencing, which could be improved. The methods presented in this work could help to identify the known deletions and duplications in a set of patients, while also helping to filter out erroneous calls for these variants, which might aid the efforts to characterize a not-yet well-known fraction of genetic variability in the human genome.