- 1Center of Excellence in Environmental Toxicology, Department of Systems Pharmacology and Translational Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
- 2Penn/CHOP Center of Excellence in Friedreich’s Ataxia, Philadelphia, PA, United States
- 3Departments of Pediatrics and Neurology, Children’s Hospital of Philadelphia, Philadelphia, PA, United States
- 4Departments of Pediatrics and Neurology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
Friedreich’s ataxia (FRDA) is an autosomal recessive disease caused by an intronic guanine-adenine-adenine (GAA) triplet expansion in the frataxin (FXN) gene, which leads to reduced expression of full-length frataxin (1–210) also known as isoform 1. Full-length frataxin has a mitochondrial targeting sequence, which facilitates its translocation into mitochondria where it is processed through cleavage at G41-L42 and K80-S81 by mitochondrial processing (MPP) to release mitochondrial mature frataxin (81–210). Alternative splicing of FXN also leads to expression of N-terminally acetylated extra-mitochondrial frataxin (76–210) named isoform E because it was discovered in erythrocytes. Frataxin isoforms are undetectable in serum or plasma, and originally whole blood could not be used as a biomarker in brief therapeutic trials because it is present in erythrocytes, which have a half-life of 115-days and so frataxin levels would remain unaltered. Therefore, an assay was developed for analyzing frataxin in platelets, which have a half-life of only 10-days. However, our discovery that isoform E is only present in erythrocytes, whereas, mature frataxin is present primarily in short-lived peripheral blood mononuclear cells (PBMCs), granulocytes, and platelets, meant that both proteins could be quantified in whole blood samples. We now report a quantitative assay for frataxin proteoforms in whole blood from healthy controls and FRDA patients. The assay is based on stable isotope dilution coupled with immunoprecipitation (IP) and two-dimensional-nano-ultrahigh performance liquid chromatography/parallel reaction monitoring/high resolution mass spectrometry (2D-nano-UHPLC-PRM/HRMS). The lower limit of quantification was 0.5 ng/mL for each proteoform and the assays had 100% sensitivity and specificity for discriminating between healthy controls (n = 11) and FRDA cases (N = 100 in year-1, N = 22 in year-2,3). The mean levels of mature frataxin in whole blood from healthy controls and homozygous FRDA patients were significantly different (p < 0.0001) at 7.5 ± 1.5 ng/mL and 2.1 ± 1.2 ng/mL, respectively. The mean levels of isoform E in whole blood from healthy controls and homozygous FRDA patients were significantly different (p < 0.0001) at 26.8 ± 4.1 ng/mL and 4.7 ± 3.3 ng/mL, respectively. The mean levels of total frataxin in whole blood from healthy controls and homozygous FRDA patients were significantly different (p < 0.0001) at 34.2 ± 4.3 ng/mL and 6.8 ± 4.0 ng/mL, respectively. The assay will make it possible to rigorously monitor the natural history of the disease and explore the potential role of isoform E in etiology of the disease. It will also facilitate the assessment of therapeutic interventions (including gene therapy approaches) that attempt to increase frataxin protein expression as a treatment for this devastating disease.
Introduction
Friedreich’s ataxia (FRDA) is an autosomal recessive neurodegenerative disorder caused by the deficient expression of mitochondrial mature frataxin protein. FRDA is (Marmolino, 2011) considered a rare disease, as it has an estimated prevalence of one in 50,000 in the Caucasian population (Delatycki et al., 2000), which makes it the most common hereditary ataxia in this population. FRDA is characterized by progressive sensory ataxia, areflexia, dysarthria, loss of position and vibratory sense and extensor plantar response (Lynch and Seyer, 2014). As a result of its progressive nature, most patients are wheelchair bound by 15.5 ± 7.4 years (mean age ± SD) after the onset of disease (Harding, 1981; Evans-Galea et al., 2014). However, heart disease is the major cause of death, which typically occurs in the fourth decade of life (Pousset et al., 2015). At present, there is no effective treatment for FRDA, although omaveloxolone was found to improve neurological function, when compared to placebo, and so it provides a potential future therapeutic strategy (Lynch et al., 2021).
The genetic basis of most cases (approximately 97%) is a homozygous guanine-adenine-adenine (GAA) triplet repeat expansion in the first intron of the frataxin (FXN) on both alleles (GAA1 and GAA2) (Chamberlain et al., 1988; Campuzano et al., 1996), which causes epigenetic transcriptional silencing and reduced expression of full-length frataxin protein (Chamberlain et al., 1988; Santos et al., 2010). A minority of FRDA patients (<3%) are compound heterozygotes possessing point or small mutations in the non-expanded allele and a GAA repeat expansion in the other (Cossee et al., 1999; Gellera et al., 2007). This leads to a loss of functional protein by non-transcriptional mechanisms. In typical FRDA patients, the length of the shortest GAA expansion (GAA1) correlates with the severity of the disease and conversely correlates with the age of onset; longer GAA expansions result in earlier onset and a faster progression (Filla et al., 1996; Sacca et al., 2011b), The presence of compound heterozygotes can complicate confirmation of a clinical diagnosis (Anheim et al., 2012) and the phenotype of compound heterozygotes cannot be predicted with certainty (Forrest et al., 1998).
Human full-length frataxin (isoform 1, MW = 23,135 Da) is expressed as a 210 amino acid precursor protein with a N-terminal mitochondrial targeting sequence (Figure 1). A two-step proteolytic cleavage by mitochondrial processing peptidase (MPP) result in the initial formation of an intermediate form (42–210; MW 18,826 Da), which is then converted to mitochondrial mature frataxin (81–210; MW = 14,268 Da) (Gakh et al., 2002; Condo et al., 2007; Schmucker et al., 2008). In contrast, extra-mitochondrial frataxin E protein (Figure 1; 76–210; MW = 14,953) in erythrocytes lacks a mitochondrial targeting sequence and arises through alternative splicing followed by N-terminal acetylation during translation (Guo et al., 2018a). No function has yet been ascribed to frataxin isoform E, and the exact function of mature frataxin has not yet been completely defined. However, several lines of evidence strongly suggest that mature frataxin is a functional component in a series of pathways including iron-sulfur cluster assembly, iron storage, heme biosynthesis, the respiratory chain, and cellular response to oxidative stress (Campuzano et al., 1996; Marmolino, 2011).
Figure 1. Sequence of full-length, intermediate, mature form, and new isoform of frataxin (FXN isoform E). Green highlights are quantitative peptides after Asp-N digestion which represent the total amount of frataxin in the whole blood. Yellow highlights are unique Asp-N digested peptide in mature form (SGTLGHPGSL) and frataxin isoform E (Ac-MNLRKSGTLGHPGSL).
Frataxin levels have been measured by western blot analysis, electro-chemiluminescence (ECL), enzyme-linked immunosorbent assay (ELISA), and western blot from FRDA lymphoblasts, lymphocytes, and peripheral blood mononuclear cells (PBMCs) (Table 1) (Campuzano et al., 1997; Boesch et al., 2007; Boehm et al., 2011). However, these studies were limited by small samples sizes and/or by their time-consuming nature. Therefore, lateral flow immunoassay (LFI), was developed as a sensitive, rapid and accurate method for quantification of frataxin from a variety of cell types (Table 1; Deutsch et al., 2010; Sacca et al., 2011b; Plasterer et al., 2013). There is no need to purify mitochondria allowing fast one-step sample preparation and rapid assay performance (normally less than 45 min). Another immunoassay-based method for frataxin with multiplex capabilities using Luminex xMAP® technology was developed and validated (Deutsch et al., 2014). The nature of the Luminex assay allows multiple purposes for the investigation of multiple markers of disease or markers for multiple diseases from a single sample. Only a small amount of blood is required for the Luminex assay and dried blood spots were typically used for frataxin determination. In most studies, the relative frataxin levels reported show consistently reduced levels in all types of specimens compared to controls (mean = 28.6% of controls, Table 1). Interestingly, a large variation of frataxin levels was observed in heterozygous FRDA patients compared to homozygous FRDA patients (Sacca et al., 2011b; Deutsch et al., 2014; Lazaropoulos et al., 2015; Clark et al., 2017). This was most likely due to cross-reaction of the antibody with the increased levels of isoform E found in some of the heterozytes (Guo et al., 2018a).
In contrast to relative levels, absolute frataxin levels showed considerable differences between laboratories. A mean frataxin levels of 438 pg/μg (n = 5) was reported in control lymphoblasts from one study (Willis et al., 2008) while another study reported a mean level using an ECL assay in the same cells that was almost 30-times lower at 11.9 pg/μg (n = 5) (Steinkellner et al., 2010; Table 1). The concentration of frataxin in normal erythrocytes was reported to be 78.6 ng/mL or 85% of the whole blood frataxin level using LFI (Plasterer et al., 2013). In contrast, a level of 28.5 ng/mL was reported from dried blood spots by Luminex xMAP (Deutsch et al., 2014). Unfortunately, few of the studies have reported details of recombinant frataxin used to construct the standard curves (Steinkellner et al., 2010; Boehm et al., 2011; Oglesbee et al., 2013) and even fewer studies have performed vigorous method validation before clinical sample analysis (Oglesbee et al., 2013; Deutsch et al., 2014). Importantly, different forms of frataxin could not be distinguished by any of these assays so they would quantify both mitochondrial and extra-mitochondrial forms (Nachbauer et al., 2011).
Mass spectrometry (MS) has the potential to discriminate different isoforms of frataxin protein by carefully choosing signature peptides which are released from different isoforms. We previously reported a sensitive assay of quantifying mitochondrial mature frataxin (81–210) from human platelets using immunoprecipitation (IP) coupled with two dimensional-nano-ultrahigh performance liquid chromatography-parallel reaction monitoring (2D-nano-UHPLC-PRM)/MS (Guo et al., 2018b). The characterization of isoform E as the frataxin proteoform that is present only in long-lived erythrocytes (Guo et al., 2018a) has now made it possible to use whole blood (rather than much less convenient blood platelets) for frataxin quantification. We report a validated, sensitive, specific, and robust assay for simultaneous determination of the two isoforms in blood using IP coupled with Asp-N endoproteinase to cleave peptide bonds on the N-terminal side of aspartic acid residues, followed by 2D-nano-UPLC-PRM/MS analysis of selected Asp-N peptides (Figure 2). The method employed stable isotope labeling by amino acids in cell culture (SILAC)-labeled human mature frataxin as the internal standard.
Figure 2. Schematic showing IP of frataxin isoforms from whole blood and analysis by 2D-nano-UHPLC-PRM/MS.
Experimental Section
Chemicals and Materials
All reagents and solvents were LC-MS grade quality unless otherwise noted. [13C615N2]-lysine and [13C615N1]-leucine were from Cambridge Isotope Laboratories (Andover, MA, United States). Anti-frataxin mouse monoclonal antibody [17A11] Ab113691 was from Abcam (Cambridge, MA, United States). Dimethyl pimelimidate dihydrochloride (DMP), Ethylenediaminetetraacetic acid (EDTA), cOmplete™ Mini EDTA-free Easypack protease inhibitor cocktail tablets, endoproteinase Asp-N sequencing grade, DL-dithiothreitol (DTT), bovine serum albumin (BSA), human lysozyme, imidazole, glycerol, phenylmethylsulfonyl fluoride (PMSF), triethanolamine, ethanolamine, and M9, minimal salts, 5X powder, minimal microbial growth medium (M9 media) were purchased from MilliporeSigma (Billerica, MA, United States). Ni-NTA agarose resin was purchased from Qiagen (Germantown, MD, United States). LC grade water and acetonitrile were from Burdick and Jackson (Muskegon, MI, United States). Ammonium bicarbonate and acetic acid were purchased from Fisher Scientific (Pittsburgh, PA, United States). Protein G magnetic beads were obtained from Life Technologies Corporation (Grand Island, NY, United States).
Clinical Samples
Blood samples were obtained from 11-unaffected healthy control subjects and 100-homozygous FRDA patients (Table 2) and 6-heterozygous FRDA patients (Table 3). A blood pool was made from 6 FRDA patients with > 800 GAA1 repeat (Supplementary Table 1). All were enrolled in parallel in an ongoing natural history study at the Children’s Hospital of Philadelphia (Lynch et al., 2006). Written informed consent was obtained from each donor participating in the study. If subjects were under the age of 18, written informed consent was obtained from a parent and/or legal guardian. The study was approved by the Institutional Review Board (IRB) of the Children Hospital of Philadelphia (IRB Protocol # 01–002609). Venous blood was drawn in 8.5 mL purple cap Vacutainer EDTA tubes and gently invert to mix. All samples were immediately aliquoted to Eppendorf tubes and frozen at −80°C until analysis.
Table 2. Demographics for study cohort of healthy controls and study cohort of homozygous FRDA patients.
Expression and Purification of Unlabeled and Stable Isotope Labeling by Amino Acids in Cell Culture-Labeled Mature Frataxin
The expression of unlabeled and SILAC-labeled mature frataxin was performed in Escherichia coli BL21 DE3 as described previously (Guo et al., 2018b). Briefly, the coding sequence of human mature frataxin (81–210) was amplified from FXN cDNA plasmid (pTL1), then cloned into a pET21b plasmid and linked to the 6× histidine (His) sequence. The 6× His-tag fusion of frataxin was expressed in Escherichia coli BL21 DE3 in M9 media containing 1 mM MgSO4, 10 μM CaCl2, and 0.5% glucose with 100 mg/L ampicillin. For expressing unlabeled frataxin, the M9 medium was supplemented with 0.025% leucine and lysine. For expressing SILAC-labeled frataxin, the M9 medium was supplemented with 0.025% [13C615N1]-leucine and [13C615N2]-lysine. The cell pellets were collected and lysed in lysis buffer [50 mM Tris-HCl (pH 8.0), 500 mM NaCl, 10 mM imidazole, 10% glycerol, 2 mM β-mercaptoethanol, 2× protease inhibitor mix, 1 mM PMSF] containing 100 μg/mL human lysozyme. The lysate was centrifuged at 20,000 × g for 30 min at 4°C, and the supernatant was purified with Ni-NTA resin. The purity of the unlabeled mature frataxin and SILAC-labeled mature frataxin was confirmed to be > 95% by SDS-PAGE and Coomassie blue staining.
Whole Blood Sample Preparation Before Immunoprecipitation
All blood samples were thawed at room temperature, and 500 μL of each sample was mixed with 750 μL NP-40 lysis buffer (150 mM NaCl, 50 mM Tris/HCl pH 7.5, 0.5% Triton X-100, 0.5% NP-40, 1 mM DTT, 1 mM EDTA) containing protease inhibitor cocktail. The same amount of SILAC-labeled mature frataxin (20 ng) was spiked in each sample (calibrator, QC and whole blood) as an internal standard. Samples were lysed by probe sonication on ice for 30 pulses at power 5 using a sonic dismembranator (Fisher, Pittsburgh, PA, United States), followed by centrifugation at 17,000 g for 15 min at 4°C. The supernatant was transferred from the pellet and incubated with pre-made DMP-crosslinked anti-frataxin protein G beads for immunoprecipitation.
Immunoprecipitation and Asp-N Digestion
Mouse monoclonal anti-frataxin antibody (4 μg) was cross-linked to protein G beads (0.5 mg) through DMP as described previously (Guo et al., 2018b). Briefly, mouse monoclonal anti-frataxin antibody was firstly incubated with protein G beads overnight at 4°C to form the antibody coupled protein G beads. The antibody coupled beads were incubated with 13 mg/mL DMP solution 1 h at room temperature to form the stable cross-linked anti-frataxin protein G beads. The cross-linked protein G beads can be kept in PBS at 4°C for a week.
The processed whole blood samples (1.25 mL) were added into 0.5 mg anti-frataxin protein G beads to carry out immunoprecipitation at 4°C overnight under rotary agitation. The beads were washed with PBS with 0.02% Tween-20 three times and frataxin and its isoforms were eluted with 100 μL of elution buffer (90% 100 mM acetic acid/10% acetonitrile). Elutes were transferred to deactivated glass inserts (Waters, Milford, MA, United States) and dried in a vacuum concentrator (Jouan RC 10.22, Fisher, Pittsburgh, PA, United States). Samples were dissolved in 50 μL of 25 mM aqueous ammonium bicarbonate containing 100 ng of Asp-N and incubated at 37 overnight before LC-MS analysis.
Method Validation by Stable Isotope Labeling by Amino Acids in Cell Culture-Labeled Mature Frataxin
A 5% BSA and FRDA whole blood pool from six FRDA patients (Supplementary Table 1) were used for preparation of calibration standards and quality controls (QCs). Calibration standards and QCs were prepared in two types of matrix, 500 μL of 5% BSA and 500 μL of FRDA whole blood pool for the evaluation of matrix effect. In 5% BSA, calibration standards (0.5, 2.0, 4.0, 10.0, 20.0, 40.0, and 80.0 ng/mL) and QCs (1.5, 3.0, 30, and 60 ng/mL) were prepared from spiking in additional stock solutions of mature frataxin at 2, 0.2, or 0.02 ng/mL. In FRDA blood pool, calibration standards (baseline + 0, 2.0, 4.0, 10.0, 20.0, 40.0, and 80.0 ng/mL) and QCs (baseline + 0, 3.0, 30, and 60 ng/mL) were prepared from spiking in additional stock solutions of mature frataxin at 2, 0.2, or 0.02 ng/mL. The extra spiked volume was less than 90% of the matrix. The accuracy and precision were determined on four levels of QCs, lower limit of quantification (LLOQ, 1.5 ng/mL), low QC (LQC, 3.0 ng/mL), middle QC (MQC, 30 ng/mL), and high QC (HQC, 80 ng/mL) in 5% BSA. The QCs were prepared in pooled FRDA blood where the frataxin level had been determined. This value was subtracted from the values that were obtained. The intra-batch accuracy and precision were determined on QCs (n = 5) were performed and analyzed on the same day. The inter-batch accuracy and precision were determined on QCs (n = 3) were performed and analyzed on three different days. The acceptance criteria of ± 15% for calibration stands and QCs (± 20% at lowest QC) at protein level were applied.
Stability Assessment
Frataxin stability in whole blood at room temperature was evaluated in healthy control whole blood over a 24-h period. After blood sample collection, 500 μL of blood were aliquoted to eight Eppendorf tubes. The aliquots were kept at room temperature on the benchtop for 1, 2, 4, 6, 8, 12, and 24 h, then frozen at −80°C until analysis. Room temperature stability was evaluated by comparing the incubated samples to freshly prepared sample (T = 0). The freeze-thaw stability was evaluated (−80°C/room temperature) with three replicates at MQC (n = 3, 30 ng/mL).
2D-Nano-UHPLC-PRM/MS
Mass spectrometry was conducted using a Q Exactive HF Orbitrap high resolution mass spectrometer coupled to Dionex Ultimate 3000 RSLCnano with capillary flowmeter chromatographic systems (Thermo Fisher Scientific, San Jose, CA, United States). The 2D system was setup as a pre-concentration mode which was composed of a 10-port valve, one nanopump and a micropump. The LC trapping column was an Acclaim PepMap C18 cartridge (0.3 μm × 5 mm, 100 Å, Thermo Scientific) and the analytical column was a C18 AQ capillary column with a 10 μm pulled tip (75 μm × 25 cm, 3 μm particle size; Columntip, New Haven, CT, United States). The 2D-nano-UHPLC system was controlled by Xcalibur software from the Q-Exactive mass spectrometer.
Loading solvent was water/acetonitrile (99.7:0.3; v/v) containing 0.2% formic acid. Solvent A was water/acetonitrile (99.5:0.5; v/v) containing 0.1% formic acid, and solvent B was acetonitrile/water (98:2, v/v) containing 0.1% formic acid. The valve stayed at loading position (1–2) with loading solvent on the trapping column at 10 μL/min for 4 min. Then the valve changed to injection position (1–10) at which the trapping column was connected with the analytical column, and samples were back-flushed into the analytical column. The valve was changed back to the loading position (1–2) at 50 min at which the trapping column was equilibrated with loading solvent and analytical column was equilibrated with 2% B. Samples were eluted with a linear gradient at a flow rate of 0.4 μL/min. The gradient on the analytical column was as follows: 2% B at the start, 5% B at 10 min, 20% B at 15 min, 60% B at 35 min, 80% B at 40–48 min, and 2% B at 48–60 min. Nanospray was conducted using Nanospray Flex™ ion source (Thermo Scientific). MS operating conditions were as follows: spray voltage 2,500 V, ion transfer capillary temperature 275°C, ion polarity positive, S-lens RF level 55, in-source CID 2.0 eV. Both full scan and PRM were used. The full scan parameters were resolution 60,000, AGC target 2 e5, maximum IT 80 ms, scan range 350–1200 m/z. The PRM parameters were resolution 60,000, AGC target 2 e5, maximum IT 80 ms, loop count 5, MSC count 4, isolation window 1.0 m/z, NCE 25. The PRMs were scheduled for 20.5 to 22.5 min for S81GTLGHPGSL90, 24.0 to 28.0 min for D124VSFGSGVLTVKLGG138, 22.00–24.2 min for D167WTGKNWVYSH177, and 21.5–23.5 min for D199LSSLAYSGK208.
Data Analysis
Data analysis was performed using Skyline (MacCoss Laboratory, University of Washington, Seattle, WA, United States). The peak area ratio of each PRM transition for each unlabeled/light (L) peptide to labeled/heavy (H) peptide was calculated by the Skyline software and used for absolute quantification. The peptide ratios were calculated by the average L/H ratios of the three PRM transitions. The total frataxin levels were calculated from the average of the four selected quantitative peptides, the mature frataxin levels were calculated from SGTLGHPGSL peptide, and the frataxin isoform E levels were calculated by subtracting SGTLGHPGSL peptide from total frataxin. Statistical analyses were conducted using Prism 9 for macOS Version 9.3.1 (GraphPad Software, LLC.). Comparisons between groups were made using an unpaired two-tailed t-test. Correlations between GAA repeats and frataxin levels as well as between age of onset and frataxin levels were conducted using linear regression models.
Results
Selection of Quantitative Peptides
There is an arginine (R79) and a lysine (K80) before serine (S81) in the full-length frataxin sequence. Therefore, trypsin would potentially digest other all frataxin isoforms to release the same N-terminal tryptic S81GTLGHPGSLDETTYER97 peptide. In order to obtain the unique N-terminal peptides for differentiating mature frataxin and frataxin isoform E, Asp-N was chosen as the endoproteinase. Mature frataxin and frataxin isoform E release two different N-terminal peptides with Asp-N digestion, S81GTLGHPGSL90 and Ac-M76NLRKSGTLGHPGSL90 (Figure 1). The expressed and purified mature frataxin (81–210) and SILAC-labeled mature frataxin were digested by Asp-N and the proteolytic peptides were optimized under full scan and PRM/MS modes. The uniqueness of the Asp-N peptides was verified by using the human proteome downloaded from NIST as the background proteome. Initially, seven potentially useful proteolytic peptides with length at 5–12 amino acids were selected, and their precursor and all product ions were assessed. The three most intense Asp-N digested peptides chosen for quantifying total frataxin levels were: D124VSFGSGVLTVKLGG138, D157WTGKNWVYSH177, and D199LSSLAYSGK208 (Supplementary Table 2). The precursor ions exhibited 10–15-fold higher intensities than their corresponding product ions in the neat solution. However, the precursor peaks had interference peaks around or even not detected due to the high interfering background in the whole blood matrix. In contrast, product ions showed clean peaks as fewer interfering signals were present in the channels of product ions (Figure 3). Therefore, the three most intense product ions for corresponding peptides were used for quantifying both endogenous frataxin and SILAC-labeled mature frataxin in whole blood samples. Direct quantification of frataxin isoform E is not accurate since the oxidation of methionine happens during sample preparation step. Due to a variety of technical challenges (Arnesen, 2011; Aksnes et al., 2016), fully acetylated frataxin isoform E standard is not available. Therefore, isoform E levels were calculated by subtracting the amount of mature frataxin (determined from S81GTLGHPGSL90) from total frataxin determined as the mean of amounts determined from the three common peptides (D124VSFGSGVLTVKLGG138, D157WTGKNWVYSH177, and D199LSSLAYSGK208).
Figure 3. Typical chromatograms of three PRM ion transitions of four quantitative Asp-N peptides under PRM mode. Left panel is light peptide, right panel is heavy peptide. Y-axis is response intensity under PRM mode. X-axis is retention time (min). Three PRM transitions shown in blue, purple and red.
Calibration Curves and Lower Limit of Quantitation
A 5% BSA and a whole blood pool from six FRDA patients (Supplementary Table 1) were used as a surrogate matrix and the biological matrix respectively in method validation. Calibration curves were constructed at six different concentrations ranging from 0.5 to 80.0 ng/mL in 5% BSA solution, and from baseline to baseline with spiked-in 80 ng/mL frataxin standard in the FRDA blood pool. Linear standard curves were obtained for each of the four peptides with r2 values range between 0.9986 to 0.9998 in 5% BSA matrix and between 0.9865 to 0.9983 in FRDA blood pool (Supplementary Figure 1). The LLOQ was set at 1.5 ng/mL in 5% BSA, which is below the mean concentration of total frataxin reported in whole blood from healthy controls and FRDA patients determined previously (Table 1).
Accuracy and Precision
Method accuracy and precision were obtained for the analysis of four levels of QCs in whole blood for the LLOQ (1.5 ng/mL), LQC (3.0 ng/mL), MQC (30 ng/mL), and HQC (60 ng/mL). The intra-day precision (N = 5) for LLOQ, LQC, MQC, and HQC (Supplementary Table 3) and inter-day precision for LQC, MQC, and HQC (Supplementary Table 4) precision was better than ± 15% for total frataxin. The accuracy was better than 90–110%. Endogenous frataxin isoforms in the FRDA whole blood that was used to prepare QCs ere determined and subtracted from the determined amounts.
Matrix Effect
A matrix effect could alter the absolute responses of the analyte and internal standard but would not change the L/H ratio of the response. Calibration curves generated from 5% BSA and from FRDA whole blood were compared (Supplementary Figure 1 and Supplementary Table 5). Two sets of calibration curves were parallel to each other. However, the matrix effect did not affect assay accuracy and precision of the assay.
In vitro Sample Stability in Whole Blood
The frataxin levels in the whole blood were stable over a 24-h period at room temperature before sample preparation (Supplementary Table 6). The frataxin levels were the same after two freeze-thaw cycles in three replicates of whole blood (Supplementary Table 7). LQC, MQC, and HQC samples that were re-analyzed after 48 h standing in the autosampler (4°C) gave essentially identical data to that obtained from the original analyses (data not shown).
Incurred Sample Reanalysis (ISR)
Twelve FRDA whole blood samples were re-analyzed. The deviations of the first analysis and repeat analysis for each individual for mature frataxin, isoform E and total frataxin ranged from −26.6 to 18.7% (Supplementary Table 8). The mean value of the first analysis and repeat analysis for mature frataxin, frataxin isoform E and total frataxin was −2.2, −13.3, and −10.8%, respectively (Supplementary Table 8). This demonstrated that data were reproducible in the whole blood samples.
Frataxin Levels in Controls and Friedreich’s Ataxia Cases
The demographic information for the homozygous and heterogenous FRDA patients enrolled in this study are summarized in Tables 2, 3. All frataxin levels were within the linear range of standard curves. The mean concentrations (± SD) of mature frataxin in healthy controls and homozygous FRDA cases in blood samples taken at year-1 were significantly different (p < 0.0001) at 7.5 ± 1.5 ng/mL and 2.1 ± 1.2 ng/mL, respectively (Figure 4A). The mean concentrations (± SD) of isoform E in healthy controls and homozygous FRDA cases in blood samples taken at year-1 were significantly different (p < 0.0001) at 26.8 ± 4.1 ng/mL and 4.7 ± 3.3 ng/mL, respectively (Figure 4B). The mean concentrations (± SD) of total frataxin in healthy controls and homozygous FRDA cases in blood samples taken at year-1 were significantly different (p < 0.001) at 34.2 ± 4.2 ng/mL and 6.8 ± 4.0 ng/mL, respectively (Figure 4C).
Figure 4. Whole blood frataxin levels in healthy controls and homozygous FRDA patients as means ± standard deviation. (A) Mature frataxin. Healthy controls 7.5 ± 1.5 ng/mL blood (n = 11), homozygous cases in year-1 2.1 ± 1.2 ng/mL blood (n = 100, p < 0.0001), homozygous FRDA patients in year-2 2.1 ± 1.4 ng/mL blood (n = 22, p < 0.0001), homozygous FRDA patients in year-3 2.0 ± 1.1 ng/mL blood (n = 22, p < 0.0001). (B) Frataxin isoform E. Healthy controls 26.8 ± 4.1 ng/mL blood (n = 11), homozygous FRDA patients in year-1 4.7 ± 3.3 ng/mL blood (n = 100, p < 0.0001), homozygous FRDA patients in year-2 2.5 ± 2.0 ng/mL blood (n = 22, p < 0.0001), homozygous FRDA patients in year-3 2.5 ± 1.9 ng/mL blood (n = 22, p < 0.0001). (C) Total frataxin. Healthy controls 34.2 ± 4.2 ng/mL blood (n = 11), homozygous FRDA patients in year-1 6.8 ± 4.0 ng/mL blood (n = 100, p < 0.0001), homozygous FRDA patients in year-2 4.5 ± 2.6 ng/mL blood (n = 22, p < 0.0001), homozygous FRDA patients in year-3 4.5 ± 2.8 ng/mL blood (n = 22, p < 0.0001). Homozygous FRDA patients were compared with healthy controls using an unpaired two-tailed t-test.
Discussion
It has been reported that recombinant yeast and human frataxin are able to self-associate in large molecular assemblies as a mechanism for detoxifying redox-active iron (Cavadini et al., 2002; O’Neill et al., 2005). Therefore, to effectively release frataxin proteoforms from whole blood cells, different types of lysis buffers were screened. We found that the non-ionic detergents, NP-40 and Triton X-100 were effective and also less harsh than ionic detergents such as sodium dodecyl sulfate (SDS). Although these two detergents could potentially have affected chromatographic separations and contaminated the mass spectrometer, they were efficiently removed by washing the beads used for IP with PBS. This made it possible to quantify endogenous frataxin levels in whole blood from 11 controls and 100 FRDA patients using the validated stable isotope dilution 2D-nano-UHPLC-PRM/HRMS method. There was no overlap of mature frataxin, isoform E, or total frataxin levels between controls and homozygous FRDA patients using this assay (Figure 4). Therefore, the assay had 100% specificity and 100% sensitivity for distinguishing healthy controls from homozygous FRDA cases. Assays on a sub-set of patients conducted on samples obtained at year-2 and year-3 from the original year sampling were found to give very similar results for mature frataxin, isoform E and total frataxin (Figure 4). This confirmed that the assay is very robust and that the frataxin levels do not change significantly over a 2-year period. The levels of isoform E frataxin in whole blood from healthy controls were similar to those reported previously for frataxin in erythrocyte fractions obtained using a dipstick LFI (Plasterer et al., 2013) most likely because the dipstick assay could not distinguish mature frataxin from isoform E.
A majority of FRDA patients (homozygous, 97%) carry GAA repeat expansions on both alleles, while a subgroup of patients carries a missense mutation on one allele and a GAA repeat expansion on the other (heterozygous, 3%). (Chamberlain et al., 1988; Campuzano et al., 1996). The phenotype of heterozygous patients cannot be predicted with certainty and the frataxin levels vary with specific point mutations. In our cohort, six of the patients were heterozygotes (Table 3). The c.317 T > C mutation (L106S) was reported as a type of missense mutation in the core of the protein, which was predicted to be a requirement for maintenance of a stable and compact frataxin structure (Dhe-Paganon et al., 2000). Missense mutations on this site could lead to folding abnormalities and yield essentially no protein (Lazaropoulos et al., 2015). In agreement with this prediction, a heterozygous FRDA subject with the L106S point mutation had lower levels of mature frataxin, isoform E, and total frataxin than expected based on the GAA1 repeat length of 832 (Figure 5, red circle). In contrast, the FRDA patient with a c.100delG point mutation (A34P) had higher levels of total frataxin, mature frataxin and frataxin isoform E than predicted based on a GAA1 repeat length of 733 (Figure 5, blue square). The c.100delG point mutation was found in exon 1 nt100, was predicted to cause a frameshift and premature termination of frataxin at codon 7 to theoretically give M1S at the amino terminus of the protein (Gellera et al., 2007). A heterozygous FRDA subject who carried a c.2delT mutation (M1S) had highly discrepant levels of frataxin between buccal cells (22% of controls) and whole blood (157% of controls) that were reported in a previous study (Lazaropoulos et al., 2015). The M1S mutation would result in a failure of initiation of translation and together with the effect of GAA1 repeat expansion reduced levels of mature frataxin would result (Zhu et al., 2002). We found a low level of mature frataxin (Figure 5 green square) for a subject with the same M1S point mutation. However, total frataxin was highly elevated because of the high levels of isoform E (Figure 5, green square). This subject was the only FRDA patient with total frataxin and frataxin isoform E levels, that fell into the range typically found for healthy controls (Figure 4).
Figure 5. Whole blood frataxin levels in controls and compound heterozygous cases in year 1 (n = 6). M1S (c.2 del T); A34P (c.100 del G); L106S (c.317 T > C); XX (g.1 1-5 G > C); G130V (p.389 G > T); ● unknown deletion.
The c.100 del G point mutation (A34P) resulted in levels of isoform E and total frataxin that were higher than was observed in most of the homozygous FRDA patients even though this patient had 733 GAA1 repeats (Figure 5, blue square). In contrast, the levels of mature frataxin were typical of those found in homozygous cases with this number of GAA1 repeats (Figure 5, blue square). Blood samples from the other three heterozygous FRDA patients that were analyzed including the patient with a c.100 del G mutation (G130V), had levels typical of homozygous FRDA cases. Furthermore, none of frataxin proteins that were analyzed in blood samples from the heterozygous FRDA cases had the mutations predicted from their mRNA sequences. Therefore, these proteins were either unstable and not present in the blood cells or they were lost during our IP procedure.
In agreement with previous studies (Filla et al., 1996; Deutsch et al., 2010; Plasterer et al., 2013), we observed that in homozygous FRDA patients, mature frataxin, isoform E, and total frataxin levels were inversely correlated with GAA1 repeat lengths (Figure 6). Furthermore, there a direct correlation between mature frataxin, isoform E, and total frataxin levels with age of onset of the disease (Figure 7). Therefore, the availability this validated robust, specific, and sensitive UHPLC-MS assay will make it possible to monitor mature frataxin, isoform E, and total frataxin levels for future studies on the natural history of the disease. It will also facilitate studies to explore the potential role of isoform E in etiology of the disease, particularly in heterozygous FRDA cases. Importantly, it will be possible with this assay to rigorously assess the efficacy of future therapeutic interventions (including gene therapy approaches) that attempt to increase frataxin protein expression as a treatment for this devastating disease.
Figure 6. Correlations of GAA repeat length on the shorter allele (GAA1) with whole blood frataxin levels. (A) Mature frataxin (p < 0.0001). (B) Isoform E (p < 0.0001). (C) Total frataxin (p < 0.0001).
Figure 7. Correlations of age of onset with whole blood frataxin levels. (A) Mature frataxin (p < 0.0001). (B) Isoform E (p < 0.0001). (C) Total frataxin (p < 0.0001).
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Children’s Hospital of Philadelphia. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
QW, LL, LW, NE, and CM conducted the analyses. LH coordinated patient recruitment and sample acquisition. DL and IB directed the project and conceived the bioanalytical approach. QW, CM, DL, and IB drafted the original manuscript. All authors edited the final version of the manuscript.
Funding
We gratefully acknowledge support of the Hamilton and Finneran families, the Penn Medicine/CHOP Friedreich’s Ataxia Center of Excellence, and NIH grants P30ES013508 (IB), 5U01NS114143 (IB), 1R01FD006029 (DL), and R21TR003035 (DL).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnins.2022.874768/full#supplementary-material
References
Aksnes, H., Drazic, A., Marie, M., and Arnesen, T. (2016). First things first: vital protein marks by N-Terminal acetyltransferases. Trends Biochem. Sci. 41, 746–760. doi: 10.1016/j.tibs.2016.07.005
Anheim, M., Mariani, L. L., Calvas, P., Cheuret, E., Zagnoli, F., Odent, S., et al. (2012). Exonic deletions of FXN and early-onset friedreich ataxia. Arch. Neurol. 69, 912–916. doi: 10.1001/archneurol.2011.834
Arnesen, T. (2011). Towards a functional understanding of protein N-terminal acetylation. PLoS Biol. 9:e1001074. doi: 10.1371/journal.pbio.1001074
Boehm, T., Scheiber-Mojdehkar, B., Kluge, B., Goldenberg, H., Laccone, F., and Sturm, B. (2011). Variations of frataxin protein levels in normal individuals. Neurol. Sci. 32, 327–330. doi: 10.1007/s10072-010-0326-1
Boesch, S., Sturm, B., Hering, S., Goldenberg, H., Poewe, W., and Scheiber-Mojdehkar, B. (2007). Friedreich’s ataxia: clinical pilot trial with recombinant human erythropoietin. Ann. Neurol. 62, 521–524. doi: 10.1002/ana.21177
Campuzano, V., Montermini, L., Lutz, Y., Cova, L., Hindelang, C., Jiralerspong, S., et al. (1997). Frataxin is reduced in Friedreich ataxia patients and is associated with mitochondrial membranes. Hum. Mol. Genet. 6, 1771–1780. doi: 10.1093/hmg/6.11.1771
Campuzano, V., Montermini, L., Molto, M. D., Pianese, L., Cossee, M., Cavalcanti, F., et al. (1996). Friedreich’s ataxia: autosomal recessive disease caused by an intronic GAA triplet repeat expansion. Science 271, 1423–1427. doi: 10.1126/science.271.5254.1423
Cavadini, P., O’neill, H. A., Benada, O., and Isaya, G. (2002). Assembly and iron-binding properties of human frataxin, the protein deficient in friedreich ataxia. Hum. Mol. Genet. 11, 217–227. doi: 10.1093/hmg/11.3.217
Chamberlain, S., Shaw, J., Rowland, A., Wallis, J., South, S., Nakamura, Y., et al. (1988). Mapping of mutation causing friedreich’s ataxia to human chromosome 9. Nature 334, 248–250. doi: 10.1038/334248a0
Clark, E., Butler, J. S., Isaacs, C. J., Napierala, M., and Lynch, D. R. (2017). Selected missense mutations impair frataxin processing in friedreich ataxia. Ann. Clin. Transl. Neurol. 4, 575–584. doi: 10.1002/acn3.433
Condo, I., Ventura, N., Malisan, F., Rufini, A., Tomassini, B., and Testi, R. (2007). In vivo maturation of human frataxin. Hum. Mol. Genet. 16, 1534–1540. doi: 10.1093/hmg/ddm102
Cossee, M., Durr, A., Schmitt, M., Dahl, N., Trouillas, P., Allinson, P., et al. (1999). Friedreich’s ataxia: point mutations and clinical presentation of compound heterozygotes. Ann. Neurol. 45, 200–206. doi: 10.1002/1531-8249(199902)45:2<200::aid-ana10>3.0.co;2-u
Delatycki, M. B., Williamson, R., and Forrest, S. M. (2000). Friedreich ataxia: an overview. J. Med. Genet. 37, 1–8. doi: 10.1136/jmg.37.1.1
Deutsch, E. C., Oglesbee, D., Greeley, N. R., and Lynch, D. R. (2014). Usefulness of frataxin immunoassays for the diagnosis of Friedreich ataxia. J. Neurol. Neurosurg. Psychiatry 85, 994–1002. doi: 10.1136/jnnp-2013-306788
Deutsch, E. C., Santani, A. B., Perlman, S. L., Farmer, J. M., Stolle, C. A., Marusich, M. F., et al. (2010). A rapid, noninvasive immunoassay for frataxin: utility in assessment of Friedreich ataxia. Mol. Genet. Metab. 101, 238–245. doi: 10.1016/j.ymgme.2010.07.001
Dhe-Paganon, S., Shigeta, R., Chi, Y. I., Ristow, M., and Shoelson, S. E. (2000). Crystal structure of human frataxin. J. Biol. Chem. 275, 30753–30756. doi: 10.1074/jbc.C000407200
Evans-Galea, M. V., Pebay, A., Dottori, M., Corben, L. A., Ong, S. H., Lockhart, P. J., et al. (2014). Cell and gene therapy for Friedreich ataxia: progress to date. Hum. Gene. Ther. 25, 684–693. doi: 10.1089/hum.2013.180
Filla, A., De Michele, G., Cavalcanti, F., Pianese, L., Monticelli, A., Campanella, G., et al. (1996). The relationship between trinucleotide (GAA) repeat length and clinical features in Friedreich ataxia. Am. J. Hum. Genet. 59, 554–560.
Forrest, S. M., Knight, M., Delatycki, M. B., Paris, D., Williamson, R., King, J., et al. (1998). The correlation of clinical phenotype in Friedreich ataxia with the site of point mutations in the FRDA gene. Neurogenetics 1, 253–257. doi: 10.1007/s100480050037
Gakh, O., Cavadini, P., and Isaya, G. (2002). Mitochondrial processing peptidases. Biochim. Biophys. Acta 1592, 63–77. doi: 10.1016/s0167-4889(02)00265-3
Gellera, C., Castellotti, B., Mariotti, C., Mineri, R., Seveso, V., Didonato, S., et al. (2007). Frataxin gene point mutations in Italian friedreich ataxia patients. Neurogenetics 8, 289–299. doi: 10.1007/s10048-007-0101-5
Guo, L., Wang, Q., Weng, L., Hauser, L. A., Strawser, C. J., Mesaros, C., et al. (2018a). Characterization of a new N-terminally acetylated extra-mitochondrial isoform of frataxin in human erythrocytes. Sci. Rep. 8:17043. doi: 10.1038/s41598-018-35346-y
Guo, L., Wang, Q., Weng, L., Hauser, L. A., Strawser, C. J., Rocha, A. G., et al. (2018b). Liquid chromatography-high resolution mass spectrometry analysis of platelet frataxin as a protein biomarker for the rare disease friedreich’s ataxia. Anal. Chem. 90, 2216–2223. doi: 10.1021/acs.analchem.7b04590
Harding, A. E. (1981). Friedreich’s ataxia: a clinical and genetic study of 90 families with an analysis of early diagnostic criteria and intrafamilial clustering of clinical features. Brain 104, 589–620. doi: 10.1093/brain/104.3.589
Lazaropoulos, M., Dong, Y., Clark, E., Greeley, N. R., Seyer, L. A., Brigatti, K. W., et al. (2015). Frataxin levels in peripheral tissue in friedreich ataxia. Ann. Clin. Transl. Neurol. 2, 831–842. doi: 10.1002/acn3.225
Lynch, D. R., Chin, M. P., Delatycki, M. B., Subramony, S. H., Corti, M., Hoyle, J. C., et al. (2021). Safety and efficacy of omaveloxolone in friedreich ataxia (MOXIe Study). Ann. Neurol. 89, 212–225. doi: 10.1002/ana.25934
Lynch, D. R., Farmer, J. M., Tsou, A. Y., Perlman, S., Subramony, S. H., Gomez, C. M., et al. (2006). Measuring Friedreich ataxia: complementary features of examination and performance measures. Neurology 66, 1711–1716. doi: 10.1212/01.wnl.0000218155.46739.90
Lynch, D. R., and Seyer, L. (2014). Friedreich ataxia: new findings, new challenges. Ann. Neurol. 76, 487–488. doi: 10.1002/ana.24259
Marmolino, D. (2011). Friedreich’s ataxia: past, present and future. Brain Res. Rev. 67, 311–330. doi: 10.1016/j.brainresrev.2011.04.001
Nachbauer, W., Wanschitz, J., Steinkellner, H., Eigentler, A., Sturm, B., Hufler, K., et al. (2011). Correlation of frataxin content in blood and skeletal muscle endorses frataxin as a biomarker in Friedreich ataxia. Mov. Disord. 26, 1935–1938. doi: 10.1002/mds.23789
Oglesbee, D., Kroll, C., Gakh, O., Deutsch, E. C., Lynch, D. R., Gavrilova, R., et al. (2013). High-throughput immunoassay for the biochemical diagnosis of Friedreich ataxia in dried blood spots and whole blood. Clin. Chem. 59, 1461–1469. doi: 10.1373/clinchem.2013.207472
O’Neill, H. A., Gakh, O., and Isaya, G. (2005). Supramolecular assemblies of human frataxin are formed via subunit-subunit interactions mediated by a non-conserved amino-terminal region. J. Mol. Biol. 345, 433–439. doi: 10.1016/j.jmb.2004.10.074
Plasterer, H. L., Deutsch, E. C., Belmonte, M., Egan, E., Lynch, D. R., and Rusche, J. R. (2013). Development of frataxin gene expression measures for the evaluation of experimental treatments in Friedreich’s ataxia. PLoS One 8:e63958. doi: 10.1371/journal.pone.0063958
Pousset, F., Legrand, L., Monin, M. L., Ewenczyk, C., Charles, P., Komajda, M., et al. (2015). A 22-year follow-up study of long-term cardiac outcome and predictors of survival in Friedreich Ataxia. JAMA Neurol. 72, 1334–1341. doi: 10.1001/jamaneurol.2015.1855
Sacca, F., Piro, R., De Michele, G., Acquaviva, F., Antenora, A., Carlomagno, G., et al. (2011a). Epoetin alfa increases frataxin production in Friedreich’s ataxia without affecting hematocrit. Mov. Disord. 26, 739–742. doi: 10.1002/mds.23435
Sacca, F., Puorro, G., Antenora, A., Marsili, A., Denaro, A., Piro, R., et al. (2011b). A combined nucleic acid and protein analysis in Friedreich ataxia: implications for diagnosis, pathogenesis and clinical trial design. PLoS One 6:e17627. doi: 10.1371/journal.pone.0017627
Santos, R., Lefevre, S., Sliwa, D., Seguin, A., Camadro, J. M., and Lesuisse, E. (2010). Friedreich ataxia: molecular mechanisms, redox considerations, and therapeutic opportunities. Antioxid. Redox Signal. 13, 651–690. doi: 10.1089/ars.2009.3015
Schmucker, S., Argentini, M., Carelle-Calmels, N., Martelli, A., and Puccio, H. (2008). The in vivo mitochondrial two-step maturation of human frataxin. Hum. Mol. Genet. 17, 3521–3531. doi: 10.1093/hmg/ddn244
Selak, M. A., Lyver, E., Micklow, E., Deutsch, E. C., Onder, O., Selamoglu, N., et al. (2011). Blood cells from Friedreich ataxia patients harbor frataxin deficiency without a loss of mitochondrial function. Mitochondrion 11, 342–350. doi: 10.1016/j.mito.2010.12.003
Steinkellner, H., Scheiber-Mojdehkar, B., Goldenberg, H., and Sturm, B. (2010). A high throughput electrochemiluminescence assay for the quantification of frataxin protein levels. Anal. Chim. Acta 659, 129–132. doi: 10.1016/j.aca.2009.11.036
Willis, J. H., Isaya, G., Gakh, O., Capaldi, R. A., and Marusich, M. F. (2008). Lateral-flow immunoassay for the frataxin protein in Friedreich’s ataxia patients and carriers. Mol. Genet. Metab. 94, 491–497. doi: 10.1016/j.ymgme.2008.03.019
Yiu, E. M., Tai, G., Peverill, R. E., Lee, K. J., Croft, K. D., Mori, T. A., et al. (2015). An open-label trial in Friedreich ataxia suggests clinical benefit with high-dose resveratrol, without effect on frataxin levels. J. Neurol. 262, 1344–1353. doi: 10.1007/s00415-015-7719-2
Keywords: assay validations, mass spectrometry, stable isotope dilution, neurodegenerative disease, liquid chromatography, homozygous Friedreich’s ataxia, heterozygous Friedreich’s ataxia
Citation: Wang Q, Laboureur L, Weng L, Eskenazi NM, Hauser LA, Mesaros C, Lynch DR and Blair IA (2022) Simultaneous Quantification of Mitochondrial Mature Frataxin and Extra-Mitochondrial Frataxin Isoform E in Friedreich’s Ataxia Blood. Front. Neurosci. 16:874768. doi: 10.3389/fnins.2022.874768
Received: 12 February 2022; Accepted: 22 March 2022;
Published: 28 April 2022.
Edited by:
Xin Qi, Case Western Reserve University, United StatesReviewed by:
Sara Anjomani Virmouni, Brunel University London, United KingdomDario Ronchi, University of Milan, Italy
Copyright © 2022 Wang, Laboureur, Weng, Eskenazi, Hauser, Mesaros, Lynch and Blair. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ian A. Blair, aWFuYmxhaXJAdXBlbm4uZWR1