Long-read sequencing improves the genetic diagnosis of retinitis pigmentosa by identifying an Alu retrotransposon insertion in the EYS gene

BACKGROUND

Biallelic variants in EYS are the major cause of autosomal recessive retinitis pigmentosa (arRP) in certain populations, a clinically and genetically heterogeneous disease that may lead to legal blindness. EYS is one of the largest genes (~ 2 Mb) expressed in the retina, in which structural variants (SVs) represent a common cause of disease. However, their identification using short-read sequencing (SRS) is not always feasible. Here, we conducted targeted long-read sequencing (T-LRS) using adaptive sampling of EYS on the MinION sequencing platform (Oxford Nanopore Technologies) to definitively diagnose an arRP family, whose affected individuals (n = 3) carried the heterozygous pathogenic deletion of exons 32-33 in the EYS gene. As this was a recurrent variant identified in three additional families in our cohort, we also aimed to characterize the known deletion at the nucleotide level to assess a possible founder effect.

RESULTS

T-LRS in family A unveiled a heterozygous AluYa5 insertion in the coding exon 43 of EYS (chr6(GRCh37):g.64430524_64430525ins352), which segregated with the disease in compound heterozygosity with the previously identified deletion. Visual inspection of previous SRS alignments using IGV revealed several reads containing soft-clipped bases, accompanied by a slight drop in coverage at the Alu insertion site. This prompted us to develop a simplified program using grep command to investigate the recurrence of this variant in our cohort from SRS data. Moreover, LRS also allowed the characterization of the CNV as a ~ 56.4kb deletion spanning exons 32-33 of EYS (chr6(GRCh37):g.64764235_64820592del). The results of further characterization by Sanger sequencing and linkage analysis in the four families were consistent with a founder variant.

CONCLUSIONS

To our knowledge, this is the first report of a mobile element insertion into the coding sequence of EYS, as a likely cause of arRP in a family. Our study highlights the value of LRS technology in characterizing and identifying hidden pathogenic SVs, such as retrotransposon insertions, whose contribution to the etiopathogenesis of rare diseases may be underestimated.

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.

Full-length antithrombin frameshift variant with aberrant C-terminus causes endoplasmic reticulum retention with a dominant-negative effect

Antithrombin, a major endogenous anticoagulant, is a serine protease inhibitor (serpin). We characterized the biological and clinical impact of variants involving C-terminal antithrombin. We performed comprehensive molecular, cellular, and clinical characterization of patients with C-terminal antithrombin variants from a cohort of 444 unrelated individuals with confirmed antithrombin deficiency. We identified 17 patients carrying 12 C-terminal variants, 5 of whom had the p.Arg445Serfs*17 deletion. Five missense variants caused qualitative deficiency, and 7, including 4 insertion-deletion variants, induced severe quantitative deficiency, particularly p.Arg445Serfs*17 (antithrombin <40%). This +1 frameshift variant had a molecular size similar to that of WT antithrombin but possessed a different C-terminus. Morphologic and cotransfection experiments showed that recombinant p.Arg445Serfs*17 was retained at the endoplasmic reticulum and had a dominant-negative effect on WT antithrombin. Characterization of different 1+ frameshift, aberrant C-terminal variants revealed that protein secretion was determined by frameshift site. The introduction of Pro441 in the aberrant C-terminus, shared by 5 efficiently secreted variants, partially rescued p.Arg445Serfs*17 secretion. C-terminal antithrombin mutants have notable heterogeneity, related to variant type and localization. Aberrant C-terminal variants caused by 1+ frameshift, with similar size as WT antithrombin, may be secreted or not, depending on frameshift site. The severe clinical phenotypes of these genetic changes are consistent with their dominant-negative effects.

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.

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

The identification of inherited antithrombin deficiency (ATD) is critical to prevent potentially life-threatening thrombotic events. Causal variants in SERPINC1 are identified for up to 70% of cases, the majority being single-nucleotide variants and indels. The detection and characterization of structural variants (SVs) in ATD remain challenging due to the high number of repetitive elements in SERPINC1. Here, we performed long-read whole-genome sequencing on 10 familial and 9 singleton cases with type I ATD proven by functional and antigen assays, who were selected from a cohort of 340 patients with this rare disorder because genetic analyses were either negative, ambiguous, or not fully characterized. We developed an analysis workflow to identify disease-associated SVs. This approach resolved, independently of its size or type, all eight SVs detected by multiple ligation-dependent probe amplification, and identified for the first time a complex rearrangement previously misclassified as a deletion. Remarkably, we identified the mechanism explaining ATD in 2 out of 11 cases with previous unknown defect: the insertion of a novel 2.4 kb SINE-VNTR-Alu retroelement, which was characterized by de novo assembly and verified by specific polymerase chain reaction amplification and sequencing in the probands and affected relatives. The nucleotide-level resolution achieved for all SVs allowed breakpoint analysis, which revealed repetitive elements and microhomologies supporting a common replication-based mechanism for all the SVs. Our study underscores the utility of long-read sequencing technology as a complementary method to identify, characterize, and unveil the molecular mechanism of disease-causing SVs involved in ATD, and enlarges the catalogue of genetic disorders caused by retrotransposon insertions.

Molecular Dissection of Structural Variations Involved in Antithrombin Deficiency

Inherited antithrombin deficiency, the most severe form of thrombophilia, is predominantly caused by variants in SERPINC1. Few causal structural variants have been described, usually detected by multiplex ligation-dependent probe amplification or cytogenetic arrays, which only define the gain or loss and the approximate size and location. This study has done a complete dissection of the structural variants affecting SERPINC1 of 39 unrelated patients with antithrombin deficiency using multiplex ligation-dependent probe amplification, comparative genome hybridization array, long-range PCR, and whole genome nanopore sequencing. Structural variants, in all cases only affecting one allele, were deleterious and caused a severe type I deficiency. Most defects were deletions affecting exons of SERPINC1 (82.1%), but the whole cohort was heterogeneous, as tandem duplications, deletion of introns, or retrotransposon insertions were also detected. Their size was also variable, ranging from 193 bp to 8 Mb, and in 54% of the cases involved neighboring genes. All but two structural variants had repetitive elements and/or microhomologies in their breakpoints, suggesting a common mechanism of formation. This study also suggested regions recurrently involved in structural variants causing antithrombin deficiency and found three structural variants with a founder effect: the insertion of a retrotransposon, duplication of exon 6, and a 20-gene deletion. Finally, nanopore sequencing was determined to be the most appropriate method to identify and characterize all structural variants at nucleotide level, independently of their size or type.

High penetrance of inferior vena cava system atresia in severe thrombophilia caused by homozygous antithrombin Budapest 3 variant: Description of a new syndrome

Atresia of inferior vena cava (IVC) is a rare congenital malformation associated with high risk of venous thrombosis that still has unknown etiology, although intrauterine IVC thrombosis has been suggested to be involved. The identification of IVC atresia in a case with early idiopathic venous thrombosis and antithrombin deficiency caused by the homozygous SERPINC1 c.391C > T variant (p.Leu131Phe; antithrombin Budapest 3) encouraged us to evaluate the role of this severe thrombophilia in this vascular abnormality. We have done a cross-sectional study in previously identified cohorts of patients homozygous for the Budapest 3 variant (N = 61) selected from 1118 patients with congenital antithrombin deficiency identified in two different populations: Spain (N = 692) and Hungary (N = 426). Image analysis included computed tomography and phlebography. Atresia of the IVC system was observed in 17/24 cases (70.8%, 95% confidence interval [CI]: 48.9%-87.3%) homozygous for antithrombin Budapest 3 with available computed tomography (5/8 and 12/16 in the Spanish and Hungarian cohorts, respectively), 16 had an absence of infrarenal IVC and one had atresia of the left common iliac vein. All cases with vascular defects had compensatory mechanisms, azygos-hemiazygos continuation or double IVC, and seven also had other congenital anomalies. Short tandem repeat analysis supported the specific association of the IVC system atresia with SERPINC1. We show the first evidence of the association of a severe thrombophilia with IVC system atresia, supporting the possibility that a thrombosis in the developing fetal vessels is the reason for this anomaly. Our hypothesis-generating results encourage further studies to investigate severe thrombophilic states in patients with atresia of IVC.

Prognostic value of thrombin generation parameters in hospitalized COVID-19 patients

SARS-CoV-2 infection increases the risk of thrombosis by different mechanisms not fully characterized. Although still debated, an increase in D-dimer has been proposed as a first-line hemostasis test associated with thromboembolic risk and unfavorable prognosis. We aim to systematically and comprehensively evaluate the association between thrombin generation parameters and the inflammatory and hypercoagulable state, as well as their prognostic value in COVID-19 patients. A total of 127 hospitalized patients with confirmed COVID-19, 24 hospitalized patients with SARS-CoV-2-negative pneumonia and 12 healthy subjects were included. Clinical characteristics, thrombin generation triggered by tissue factor with and without soluble thrombomodulin, and also by silica, as well as other biochemical parameters were assessed. Despite the frequent use of heparin, COVID-19 patients had similar thrombin generation to healthy controls. In COVID-19 patients, the thrombin generation lag-time positively correlated with markers of cell lysis (LDH), inflammation (CRP, IL-6) and coagulation (D-dimer), while the endogenous thrombin potential (ETP) inversely correlated with D-dimer and LDH, and positively correlated with fibrinogen levels. Patients with more prolonged lag-time and decreased ETP had higher peak ISTH-DIC scores, and had more severe disease (vascular events and death). The ROC curve and Kaplan Meier estimate indicated that the D-dimer/ETP ratio was associated with in-hospital mortality (HR 2.5; p = 0.006), and with the occurrence of major adverse events (composite end-point of vascular events and death) (HR 2.38; p = 0.004). The thrombin generation ETP and lag-time variables correlate with thromboinflammatory markers, and the D-dimer/ETP ratio can predict major adverse events in COVID-19.

Antithrombin p.Thr147Ala: The First Founder Mutation in People of African Origin Responsible for Inherited Antithrombin Deficiency

BACKGROUND

 Hereditary antithrombin deficiency is a rare autosomal-dominant disorder predisposing to recurrent venous thromboembolism (VTE). To date, only two founder mutations have been described.

OBJECTIVES

 We investigated the antithrombin p.Thr147Ala variant, found in 12 patients of African origin. This variant is known as rs2227606 with minor allele frequency of 0.5% in Africans and absent in Europeans. A possible founder effect was investigated.

METHODS

 Phenotypical characterization was established through immunological and functional methods, both under basal and stress conditions. Recombinant antithrombin molecules were constructed by site-directed mutagenesis and expressed in HEK-293T cells. Secreted antithrombin was purified and functionally characterized. Structural modeling was performed to predict the impact of the mutation on protein structure. A novel nanopore sequencing approach was used for haplotype investigation.

RESULTS

 Ten patients experienced VTE, stroke, or obstetric complications. Antithrombin antigen levels and anti-IIa activity were normal or slightly reduced while anti-Xa activity was reduced with only one commercial assay. On crossed immunoelectrophoresis, an increase of antithrombin fractions with reduced heparin affinity was observed under high ionic strength conditions but not under physiological conditions. The recombinant p.Thr147Ala protein displayed a reduced anti-Xa activity. Structural modeling revealed that residue Thr147 forms three hydrogen bonds that are abolished when mutated to alanine. The investigated patients shared a common haplotype involving 13 SERPINC1 intragenic single nucleotide polymorphisms.

CONCLUSION

 Antithrombin p.Thr147Ala, responsible for antithrombin type II heparin binding site deficiency, is the first founder mutation reported in people of African ancestry. This study further emphasizes the limitations of commercial methods to diagnose this specific subtype.

ALG12-CDG: An unusual patient without intellectual disability and facial dysmorphism, and with a novel variant

Background

Congenital disorder of glycosylation (CDG) type I is a group of rare disorders caused by recessive mutations in up to 25 genes that impair the N‐glycan precursor formation and its transfer to proteins resulting in hypoglycosylation of multiple proteins. Congenital disorder of glycosylation causes multisystem defects usually with psychomotor delay that is diagnosed in the infancy. We aim to supply further evidences supporting that CDG may be underestimated.

Methods

Antithrombin and factor XI were studied by chromogenic and coagulometric methods. Hypoglycosylation of plasma proteins was evaluated by western blot, HPLC, Q‐TOF, and RP‐LC‐MRM‐MS. Genetic analysis included whole exome, Sanger sequencing, and PCR‐allele specific assay.

Results

We here present an intriguing patient with an exceptional phenotype: 25‐year‐old women with a ventricular septal defect and severe idiopathic scoliosis but no facial dysmorphism, who dances as a professional, and has a University degree. Congenital disorder of glycosylation diagnosis started through the identification of antithrombin deficiency without SERPINC1 defect and the detection of hypoglycosylated forms. Increased levels of hypoglycosylated forms of F XI (also with significant deficiency) and transferrin were also detected. Whole exome analysis showed a novel homozygous ALG12 variant c.77T>A, p.(Val26Asp) supporting an ALG12‐CDG diagnosis. It also showed three new variants in KMT2D, and a mild, known ALG6 variant.

Conclusions

This novel ALG12‐CDG patient (the 13th reported) underlines the heterogeneity of this CDG and broadens its phenotypical spectrum, supports that these disorders are underestimated, and suggests that combination of global hypoglycosylation with specific gene defects might determine the clinical manifestations of CDG patients.

Incidence and features of thrombosis in children with inherited antithrombin deficiency

Pediatric thromboembolism (≤18 years) is very rare (0.07-0.14/10,000/year) but may be more prevalent in children with severe thrombophilia (protein C, protein S or antithrombin deficiency). The aim of this study was to define the prevalence and clinical characteristics of pediatric thrombosis in subjects with inherited antithrombin deficiency. Our observational retrospective multicentric study from two countries recruited 968 patients of any age from 441 unrelated families with genetically, biochemically and functionally characterized antithrombin deficiency. Seventy-three subjects (7.5%) developed thrombosis before 19 years of age. Two high-risk periods for thrombosis were identified: adolescence (12-18 years, n=49) with thrombus localization (lower limb deep venous thrombosis or pulmonary embolism) and triggering factors common to adults (oral contraceptives, surgery or pregnancy); and the neonatal period (<30 days, n=15) with idiopathic thrombosis at unusual sites. The clinical evaluation of pediatric thrombosis in subjects with antithrombin deficiency revealed: i) a high prevalence of cerebral sinovenous thrombosis (n=13, 17.8%), mainly at young age (8 neonates and 4 children <6 years); ii) severe outcome with fatality in six cases (3 neonates, two of them homozygous for p.Leu131Phe). The majority of subjects (76.7%) carried quantitative type I deficiency. This retrospective analysis includes the largest cohort of subjects with inherited antithrombin deficiency so far and provides strong evidence for an increased risk of pediatric thrombosis associated with this thrombophilia (300-fold compared with the general population: 0.41%/year vs. 0.0014%/year, respectively). Our results support testing for antithrombin deficiency in children of affected families, particularly in case of type I deficiency.

High levels of latent antithrombin in plasma from patients with antithrombin deficiency

Antithrombin is an anticoagulant serpin that efficiently inhibits multiple procoagulant proteases. The cost for the structural flexibility required for this function is the vulnerability to mutations that impact its folding pathway. Most conformational mutations identified in serpins cause polymerisation. Only three mutations in SERPINC1 affecting two residues have been found to favour transformation to the latent conformation of antithrombin, another hyperstable non-anticoagulant form with strong antiangiogenic activity that constitutes 3 % of plasma antithrombin in healthy subjects. The analysis of latent antithrombin in 141 unrelated patients with antithrombin deficiency carrying 89 different SERPINC1 mutations identified four cases with higher levels than that of controls: p.Pro439Thr, p.Pro461Ser, p.Met283Val, and p.His401Tyr, the last also with circulating polymers. Heating of plasma at 42ºC exacerbated the transformation to the latent conformation in p.Pro439Thr and p.Pro461Ser. The conformational effect of p.Met283Val, the mutation associated with the highest levels of latent antithrombin detected in four members of a family, was verified in a recombinant model. Antithrombin deficiency in these cases should be classified as pleiotropic based on the impaired reactivity and low heparin affinity of the variant. Despite high levels of latent antithrombin (up to 80 µg/ml in p.Met283Val carriers), no vascular defects were described in carriers of these mutations. In conclusion, our study identifies new residues involved in the structural stability of antithrombin (and potentially of all serpins). High levels of endogenous latent antithrombin seem to play a minor antiangiogenic effect. Finally, pleiotropic deficiencies may be caused by mutations inducing transformation to the latent conformation.