Molecular Dissection of Structural Variations Involved in Antithrombin Deficiency

Long-read genome sequencing reveals a novel intronic retroelement insertion in NR5A1 associated with 46,XY differences of sexual development

Despite significant advancements in rare genetic disease diagnostics, many patients with rare genetic disease remain without a molecular diagnosis. Novel tools and methods are needed to improve the detection of disease-associated variants and understand the genetic basis of many rare diseases. Long-read genome sequencing provides improved sequencing in highly repetitive, homologous, and low-complexity regions, and improved assessment of structural variation and complex genomic rearrangements compared to short-read genome sequencing. As such, it is a promising method to explore overlooked genetic variants in rare diseases with a high suspicion of a genetic basis. We therefore applied PacBio HiFi sequencing in a large multi-generational family presenting with autosomal dominant 46,XY differences of sexual development (DSD), for whom extensive molecular testing over multiple decades had failed to identify a molecular diagnosis. This revealed a rare SINE-VNTR-Alu retroelement insertion in intron 4 of NR5A1, a gene in which loss-of-function variants are an established cause of 46,XY DSD. The insertion segregated among affected family members and was associated with loss-of-expression of alleles in cis, demonstrating a functional impact on NR5A1. This case highlights the power of long-read genome sequencing to detect genomic variants that have previously been intractable to detection by standard short-read genomic testing.

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.

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.

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.

Identification and characterization of two SERPINC1 mutations causing congenital antithrombin deficiency

BACKGROUND

Antithrombin (AT) is the main physiological anticoagulant involved in hemostasis. Hereditary AT deficiency is a rare autosomal dominant thrombotic disease mainly caused by mutations in SERPINC1, which was usually manifested as venous thrombosis and pulmonary embolism. In this study, we analyzed the clinical characteristics and screened for mutant genes in two pedigrees with hereditary AT deficiency, and the functional effects of the pathogenic mutations were evaluated.

METHODS

Candidate gene variants were analyzed by next-generation sequencing to screen pathogenic mutations in probands, followed by segregation analysis in families by Sanger sequencing. Mutant and wild-type plasmids were constructed and transfected into HEK293T cells to observe protein expression and cellular localization of SERPINC1. The structure and function of the mutations were analyzed by bioinformatic analyses.

RESULTS

The proband of pedigree A with AT deficiency carried a heterozygous frameshift mutation c.1377delC (p.Asn460Thrfs*20) in SERPINC1 (NM000488.3), a 1377C base deletion in exon 7 resulting in a backward shift of the open reading frame, with termination after translation of 20 residues, and a different residue sequence translated after the frameshift. Bioinformatics analysis suggests that the missing amino acid sequence caused by the frameshift mutation might disrupt the disulfide bond between Cys279 and Cys462 and affect the structural function of the protein. This newly discovered variant is not currently included in the ClinVar and HGMD databases. p.Arg229* resulted in a premature stop codon in exon 4, and bioinformatics analysis suggests that the truncated protein structure lost its domain of interaction with factor IX (Ala414 site) after the deletion of nonsense mutations. However, considering the AT truncation protein resulting from the p.Arg229* variant loss a great proportion of the molecule, we speculate the variant may affect two functional domains HBS and RCL and lack of the corresponding function. The thrombophilia and decreased-AT-activity phenotypes of the two pedigrees were separated from their genetic variants. After lentiviral plasmid transfection into HEK293T cells, the expression level of AT protein decreased in the constructed c.1377delC mutant cells compared to that in the wild-type, which was not only reduced in c.685C > T mutant cells but also showed a significant band at 35 kDa, suggesting a truncated protein. Immunofluorescence localization showed no significant differences in protein localization before and after the mutation.

CONCLUSIONS

The p.Asn460Thrfs*20 and p.Arg229* variants of SERPINC1 were responsible for the two hereditary AT deficiency pedigrees, which led to AT deficiency by different mechanisms. The p.Asn460Thrfs*20 variant is reported for the first time.

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.