Clinical model for Hereditary Transthyretin Amyloidosis age of onset prediction

Introduction: Hereditary transthyretin amyloidosis (ATTRv amyloidosis) is a rare neurological hereditary disease clinically characterized as severe, progressive, and life-threatening while the age of onset represents the moment in time when the first symptoms are felt. In this study, we present and discuss our results on the study, development, and evaluation of an approach that allows for time-to-event prediction of the age of onset, while focusing on genealogical feature construction.

Materials and methods: This research was triggered by the need to answer the medical problem of when will an asymptomatic ATTRv patient show symptoms of the disease. To do so, we defined and studied the impact of 77 features (ranging from demographic and genealogical to familial disease history) we studied and compared a pool of prediction algorithms, namely, linear regression (LR), elastic net (EN), lasso (LA), ridge (RI), support vector machines (SV), decision tree (DT), random forest (RF), and XGboost (XG), both in a classification as well as a regression setting; we assembled a baseline (BL) which corresponds to the current medical knowledge of the disease; we studied the problem of predicting the age of onset of ATTRv patients; we assessed the viability of predicting age of onset on short term horizons, with a classification framing, on localized sets of patients (currently symptomatic and asymptomatic carriers, with and without genealogical information); and we compared the results with an out-of-bag evaluation set and assembled in a different time-frame than the original data in order to account for data leakage.

Results: Currently, we observe that our approach outperforms the BL model, which follows a set of clinical heuristics and represents current medical practice. Overall, our results show the supremacy of SV and XG for both the prediction tasks although impacted by data characteristics, namely, the existence of missing values, complex data, and small-sized available inputs.

Discussion: With this study, we defined a predictive model approach capable to be well-understood by medical professionals, compared with the current practice, namely, the baseline approach (BL), and successfully showed the improvement achieved to the current medical knowledge.

1. Introduction

Hereditary amyloidosis related to variant transthyretin (ATTRv amyloidosis) (1) is a multisystemic disease, with predominant involvement of the heart and the peripheral nervous system. Point mutations in transthyretin, a plasma transport protein, mainly produced in the liver lead to instability of the tetrameric structure of the native protein and disassembly of the monomers that undergo a misfolding process resulting in a widespread deposition of amyloid in the extracellular space of many tissues with consequent organ dysfunction (2). Genetic heterogeneity (more than 130 pathogenic mutations have been described) and phenotypic heterogeneity (variable organ involvement and variable age of onset) are characteristics of this condition, but they are not completely understood. Phenotypic heterogeneity, either the variability of the age of onset or the predominance of a given organ involvement, is driven both by mutation and by genetic background (3). Nonetheless, patients with the same mutation originating in the same region, and even in the same family, may present significant clinical diversity. The variability of the age of onset has been extensively studied and its causes are the object of different types of research, including different strategies to identify genetic modifiers (4, 5).

In recent years, researchers developed several approaches to study and predict future symptom occurrence in patients with different, complex, and irregular diseases. A few examples are breast cancer occurrence and renal artery stenosis (6), Framingham risk functions for cardiovascular diseases (7), and, most recently, genetic frontotemporal dementia (8). In these cases, most of the research focuses on applying statistical-based analysis that is dependent on the dimension of the patient cohorts and sometimes suffers from a lack of enough follow-up temporal data.

In the case of ATTRv amyloidosis, where we have variable penetrance of the gene, the probability of a carrier of a given mutation developing the disease has been extensively studied. This variable penetrance translates into a variable age of onset from the late teens until very old ages and also explains many cases of sporadic presentation of the disease in probands, later shown to have relatives with the same mutation remaining silent throughout their lifespan (9). This is a phenomenon called non-penetrance of the gene. It has been noted in this case that penetrance has been found to be influenced by the mutation, the genetic background, the gender of the patient, and the gender of the transmitting parent (10).

The disease is present worldwide (11), but it was initially described as occurring in clusters or foci related to the same mutation, sometimes known to originate from a common founder, many centuries ago. Most of these clusters are related to the TTRVal30Met mutation (also known as p. TTRVal50Met), but other mutations also show the same behavior (12, 13). The study of large clusters of TTRVal30Met mutation such as those present in Portugal, Sweden, and Japan showed a particular relation between the symptomatic presentation and the early or later age of onset, a positive correlation of the age of onset inside the same family but also an unexpected (14) finding of anticipation of the age of onset, related to the gender of the patient, and the gender of the transmitting parent (15). For subjects belonging to an affected family, the study of the pedigree allows for the definition of the risk of being a carrier of the mutation and to proceed to a blood genetic test (preferably under genetic counseling conditions) to identify those who are carriers of the mutation causing the disease in their relatives (16). However, after this definition, it would be important to predict the age of onset of the disease, the moment when the first symptoms will appear, prompting the need to start one of the presently available disease-modifying treatments, and the moment problems commonly disabling will start impacting their personal, familial, professional, and social life. The development of a sound method to predict the age of onset of a carrier could also contribute to the organization of medical care, defining that the moment carriers should be regularly observed to detect the first symptoms or signs of the disease.

Several studies emphasized the importance of early onset of disease-modifying treatments (17) and we can even consider that pre-symptomatic treatment, when the pathogenic process has started but the functional reserve of affected tissues, such as the heart and the peripheral nerves prevent symptomatic manifestations of this silent evolution would have the best results in the long term. The design of a clinical trial targeting this pre-clinical stage would benefit from the possibility of predicting the age of onset of the asymptomatic carriers based on the knowledge of the age of onset of their relatives. Recently, in Conceição et al. (18), authors discussed the prediction of the age of onset in ATTRv patients and proposed a method currently known as PADO which is presented as based in the clinical practice. Further discussion can be found in Section 4.

We define the age of onset as the moment when a set of characteristic symptoms progressively confluent and more frequent and severe start, disregarding isolated, inconstant, and unspecific complaints. The diagnosis must always be confirmed with the identification of the mutation (when not done previously) and the confirmation of amyloid deposition in any tissue, with the least invasive biopsy, such as abdominal fat, salivary gland, or skin punch biopsies). In the rare cases of recurrent biopsies without evidence of amyloid deposition, the diagnosis is accepted only if objective and unequivocal symptoms of the disease are present and other potential causes have been excluded (19).

Predicting the age of onset for ATTRv amyloidosis patients is difficult as the disease is rare, with genetic and phenotypic heterogeneity and has a wide range of possible age of onset values. Family history can help this evaluation if other family members' ages of onset are known, but we must understand how to valorize this information according to what we know about the variability of the age of onset in many families.

In this study, we follow-up on our previous study on genealogical feature construction (20–22). Our objective is to study, develop, and evaluate an approach to allow time to event the age of onset prediction, in a robust manner, while focusing on genealogical feature construction. This methodology lays the foundation to assess the disease evolution over time.

Our main contributions are (i) the definition of a large set of genealogical-based features constructed from heterogeneous clinical data sets, (ii) the study and comparison of a pool of machine learning algorithms optimized for the age of onset prediction, and (iii) the comparison of the results with an out-of-bag evaluation set, assembled in a different time-frame than the original data.

2. Materials and methods

2.1. Data

Our data comes from Unidade Corino de Andrade, a disease specialized unit officially recognized as a Portuguese National Reference Center. It is integrated into one of the largest Portuguese public hospitals, Centro Hospitalar Universitário de Santo António, and its medical professionals treat and maintain accurate patient records since the 1940's. From their records, we started by selecting information from 2,259 patients, diagnosed and followed before 2018. Recently, in order to account and evaluate for possible leakage problems, we added a validation data set and included data from 213 extra patients, diagnosed between 2018 and the end of 2021 (see Table 1). Of note currently, there are 511 patients classified as asymptomatic carriers that are expected to show symptoms in the next few years. Current patients are spread over 932 families, and these numbers tend to increase over time, as new families are diagnosed.

TABLE 1

Table 1. Two thousand four hundred and seventy-two patients diagnosed with ATTRv amyloidosis in Unidade Corino de Andrade as of December 2021.

2.2. Baseline model

Currently, for ATTRv amyloidosis, the incorporation of a predicted age of onset to insure a more accurate asymptomatic patient follow-up is not an established practice. One of the most well-known studies that reflect on this problem can be found in Lemos et al. (15).

By taking Lemos et al. (15), we defined a baseline rule model that assumes a linear trend between the age of onset of a patient and that of the transmitting ancestor. In this case, we have

$\text{if}(fatherandson)\text{then}ŷ=parent(ageonset)-6.06,$

$\text{if}(fatheranddaughter)\text{then}ŷ=parent(ageonset)-1.23,$

$\text{if}(motherandson)\text{then}ŷ=parent(ageonset)-10.43,$

$\text{if}(motheranddaughter)\text{then}ŷ=parent(ageonset)-7.43.$

To produce comparable results, we added to the baseline a rule for cases where we do not have information regarding the gender or the age of onset of the transmitting ancestor. Therefore, for simplicity purposes, we state that the patient will most likely feel symptoms of the disease in the 2 years after the prediction time point. Thus, the last rule states that

$\text{if}(dataunknown)\text{then}ŷ=patient(current\text{\_}age)+2.$

These rules, which have different constant values that represent the average difference in the age of onset of all previously diagnosed known pairs of (parent and child), are based on the existence of an anticipation inheritance mechanism. This genetic anticipation, a common factor in most hereditary diseases, indicates that a patient has a higher risk to show symptoms of the disease prior to that of the parent(s). In this case and for each known different pair of (patient, ancestor) gender, we averaged the overall differences of each pair. These values were obtained from the current pool of patients to allow for a sound comparison between the baseline and our machine learning based approach.

We will use this baseline (BL) as a reference in the validation (results in Section 3) of a more powerful predictive approach.

2.3. Approach

In this study, we focused on answering two clinical questions formulated in a classification and a regression setting.

In the classification-based setting, our purpose is to infer “if a patient previously diagnosed with the genetic error will show symptoms of the disease in the next 4 years.” In this case, we will not be interested in the patient's point age of onset but rather in whether we will be capable of accurately diagnosing them as showing symptoms of the disease. The results will answer how accurate we are capable to diagnose a patient as symptomatic in the period between 2018 and 2021.

In the regression-based setting, our purpose is to infer “when will a patient show symptoms of the disease.” In this study, we are interested in the difference between the ŷ and the y values to accurately model this time event.

2.3.1. Patient representation

To answer these problems we start by representing a patient in specific moments in time (age marks). This is not an easy task since we are interested in a formulation that captures the highly variable risk at which a patient will show symptoms of the disease, while focusing on different sets of characteristics, namely, demographic, genealogical, and based on familial disease history. In this case, we chose to construct features that aggregate information from familial connections between patients and their family members. To do this, we implemented a version of the lowest common ancestor algorithm to have a set of valid relations for each patient (23), which helps to supply the aforementioned features. Therefore, overall, for a specific patient and for each age mark, we calculate three sets of features (see Table 2): patient feature set, first-level feature set, and extended feature set. By taking into account these sets of features we characterize our patients in different age marks and later use them to train different sets of models. These are specialized on an age mark and will generate different risks for each patient. In this case, the age marks for which we group our patients are 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58, 61, and 64 years old.

TABLE 2

Table 2. List of 77 original and generated features classified between patient, first level, and extended level.

The representation of a patient at a specific age mark is only valid if they were still alive both for the classification and regression case and if they were still considered asymptomatic by the medical team in the regression case. As an example, in Figure 1, we show, for a simulated patient, its valid age marks, which range from 22 to 30 years; that the patient is symptomatic, with an age of onset of 30 years; that the asymptomatic ages range from 22 to 28 years old; and that the patient died with 39 years. Therefore, this patient cannot be used in models for ages >39 years (age of death).

FIGURE 1

Figure 1. Conceptual model of virtual age scenarios.

2.3.2. Model construction

After the construction of the set of valid age mark instances for each patient, both for the training/ testing and validation sets, we assembled a set of functions to assist the creation of our model construction and evaluation phases (see Figure 1). To evaluate our approach, we selected a temporal validation block (6, 24). To do so, we sliced the data set into two parts: one part containing early treated patients that were used to develop the model, and another part containing the most recently treated patients. We started with a non-random split (6) in which patients that showed symptoms and were diagnosed prior to 2018 were used in the training/ testing sets to check the robustness of the model, while the patients diagnosed after 2018 and the patients that currently have not shown symptoms, but are expected to, were bagged together and used in the validation phase. The definition of training, validation, and testing is done next.

2.3.3. Model evaluation

To compare our approach with the medical baseline, we separated between the patients with known parent(s) age of onset and those without. For that purpose, we developed three types of models: data set A, containing patients with previously diagnosed ancestor(s); data set B, which contains those for whom we do not have a measurable age of onset for the ancestor; and data set F, composed of all available records.

When building a machine learning model, it is essential to evaluate its performance accurately. To do so, each data set (data set A, data set B, and data set F) was further split into three groups namely, training, testing, and validation sets.

In broad terms, the training data set constitutes the reference samples by which an algorithm learns a specific concept. Once the model has been trained, the testing data set is used to assess the performance of the model and to allow for a good tuning of specific hyperparameters capable to balance between over (when the model performs well on the training set but poorly on unseen data) and under fitting (when the model is too simple or has poor performance on both the training and test data). Finally, the validation data set is used to evaluate the model's final performance on unseen data after it has been optimized using the training and testing sets. This set should only be used once, after the model has been optimized, to provide an unbiased estimate of the model's performance on new data and to provide a better indication of how well the model will perform in the real world.

After dealing with the horizontal and vertical break of the data set, i.e., the separation between different sets of features and instances to train/ test and validate our study, we aligned (see Figure 1) each patient's instance according to the age it represented. This alignment allowed for the construction and training of different models according to the age instance of each patient. With this process, each patient has a set of prediction results, which correspond to a risk factor for each age instance. In this case, the prediction of the age of onset of a patient at each age is an independent task which allows for the bagging of sets of patients from different temporal decades. This way, we end up with an independently tuned model for each time horizon and are able to select the best set of hyperparameters and algorithms for relatively small data sets (25).

3. Results

In this section, we present the experimental set-up and our results both for the classification as well as regression based approaches.

3.1. Experimental setup

As discussed previously, for each asymptomatic patient, we considered the information known at the time at which they were 22, 25, 28, 31, 34, 37, 40, 43, 46, 49, 52, 55, 58, 61, and 64 years old (as long as they were still asymptomatic) and instantiate an individual model for each group of patients. This means that for each age group, and each data set, we have a single experiment with each of the nine different algorithms (i.e., LR, EN, LA, RI, SV, DT, RF, XG, and BL). In total, we ran 405 experiments to compare the results for different patients' ages.

As for our results, these are presented in two settings.

Regarding the question if it is viable to check if a patient will feel symptoms of the disease in the next 4 years (classification setting), we present and analyze the validation set both for asymptomatic and symptomatic patients. Therefore, for each patient, we analyzed each prediction age of the last known instance of each patient before 2018 and verify if the calculated age of onset fell in the period set between 2018 and 2021 or not. This resulted in each patient having a single prediction value, which allows us to group all the patients together, independently of the disease stage they were at (see Table 3 and Figure 2). In this case, our metrics of choice will be recall (R) also referenced as true positive rate, which indicates how well a model is capable to identify all symptomatic patients and precision (P), also referred to as positive predictive value, which indicates how precise or accurate a model's symptomatic predictions are.

TABLE 3

Table 3. Classification-based metrics evaluation.

FIGURE 2

Figure 2. Evolution of precision and recall for V(A), V(B), and V(F) data set (these and other metrics, aggregated in the age dimension, are referenced on Table 3). In the early ages, we are mostly able to identify the patients that will remain asymptomatic. This decreases in the later onset ages (e.g., before 50 years old as opposed to after 50 years old). The best results are given by SV.

On the contrary, regarding the prediction of age of onset (regression setting), we analyze the average position (rank) of the mean absolute error (MAE) and the root mean squared error (RMSE) of each pair of (algorithm and data set) results. We then compare how each model performed over the full set of experiments (see Table 4 and Figure 3). We also analyze and compare the average rank in the training/ testing sets with the one in the validation set. This allows us to infer as to the robustness of our approach and to distinguish our models between top and worst performers.

TABLE 4

Table 4. Average MAE ranks of each tuple (algorithm and data set) for model-based imputation with parameter k = 10.

FIGURE 3

Figure 3. Evolution of MAE for Val(A), Val(B), and Val(F) data sets. For the data set A case, and at the later ages, our approach ranks top when compared with the BL. This also happens in data set B and data set F, notwithstanding the fact that the average error rates of the BL decrease as current age increases. There is some instability in the LR results in the early ages, for data set B resulting of the large number of features. Overall, the SV and XG correspond to the best models for a high majority of age models created.

The mean absolute error (MAE) is an evaluation metric used to measure the average absolute difference between the actual and the predicted values in a set of data. It measures how far, on average, the predictions are from the true values, and it is obtained by taking the absolute value of the difference between each actual value and its predicted value and then taking the average of all those absolute differences. As for root mean squared error (RMSE), it corresponds to a metric used to measure the average difference between the actual and predicted values in a set of data. The main difference between MAE and RMSE is that the later (RMSE) is more sensitive to larger differences (outliers) than the first (MAE).

3.2. Is it viable to check when will a patient feel symptoms?

For the classification-based evaluation, and by considering the interval between 2018 and 2021, we had to align each patient at the beginning of 2018 and transform the predicted age of onset to a year instance. Due to the lack of a valid prediction model prior to 22 years of age, here we had to exclude nine patients, that are currently asymptomatic.

In the study, we observe that in the early ages, we are mostly able to identify the negative examples (patients that will remain asymptomatic). This result decreases in the later onset ages (e.g., before 50 years old as opposed to after 50 years old). We can see the impact by the peek of the onset age dispersing the results following the period of 30–35 years.

Overall, these results are impacted by the rather small prediction period that ranges between 2018 and 2021. The results of the top and worst performers show that although there is some fluctuation, the overall top performer is the support vector classifier (SV) algorithm for most of the chosen metrics. In the case of recall (R), most algorithms obtain good results, and SV is usually well above 0.9 except in one point for data set A. In fact, when having a closer look at the different ages, it is clear the consistency of the SV algorithm.

3.3. When will an asymptomatic patient feel symptoms of the disease?

When comparing the validation results of the MAE error distribution among the different ages (see Figure 3), it is clear that for the data set A case, and at the early ages, our approach ranks the top when compared with the BL although in later years this is not as clear. This fact is likely affected by the decrease in the validation set size throughout the age groups. It is important to reference that this approach suffers from a lack of new asymptomatic patients as the ages increase (presently, the average age of onset of the training/ testing + validation sets is 36.32 years) as well as in the change between the training/testing and validation sets (see Table 1).

As for data set B and data set F, it is clear that our approach shows the best results, with SV being the best-chosen algorithm and XG being the second-best.

Regarding the average ranks (see Tables 4, 5), and when comparing the dispersion and difference between the top two performers, we have two main observations. First, on the one hand, the different MAE ranks are not stable between the training/ testing and validation sets. When taking a closer look at the RMSE average ranks, we see that this does not happen. These rank difference between training/testing (T) and validation (V) sets in MAE results suggest that there is a difference in large error distribution in SV results. This can indicate some future problems since we want to have some assurances as to the future performance of the top selected approaches.

TABLE 5

Table 5. Average RMSE ranks of each tuple (algorithm and data set) for model-based imputation with parameter k = 10.

Second, that, there is a good degree of stability of the ranks between each vertical data set separation (e.g., between data set A, data set B, and data set F).

4. Discussion

This study explores genealogical features in modeling and predicting the age of onset of ATTRv amyloidosis patients. The definition of different sets of data explores different levels of pre-existing information regarding the patients' ancestry.

To compare our approach with current medical practice, we studied works focused on ATTRv amyloidosis guidelines regarding the age of onset of asymptomatic patients, namely, Lemos et al. (15) and Conceição et al. (18).

Our purpose was to provide a concise and clear baseline method, which follows the current medical practice. The baseline needs to be defined in a manner that enables its implementation without major adaptations.

In the first study, authors studied the evolution of the age of onset in a group of families and found the existence of an anticipation mechanism (26), a common factor in hereditary diseases. This fact has been studied for ATTRv amyloidosis in other geographic areas with the studies of Drugge et al. (27), Tashima et al. (28), and Cisneros-Barroso et al. (29). This study, which used a large enough pool set of patient records, presents a statistically based foundation that can foster the development of a wide range of analytical inferences focused on the distribution of the age of onset of ATTRv patients and its genealogical influence.

In a later study (18), a team of doctors with a wide knowledge in treating and following patients diagnosed with ATTRv amyloidosis debated the importance of assessing the predicted age of onset of asymptomatic patients (PADO). As the authors acknowledge, the definition of a prediction value can insure the correct monitoring of known asymptomatic subjects until such time that the onset is medically observed with fewer costs for the public health system. Authors mention that PADO depends on patient mutation, the typical age of onset for that mutation, and the age of onset in family members with ATTRv amyloidosis, specifically its proband (e.g., family first diagnosed case).

When comparing both studies, a few differences emerge. In Lemos et al. (15), researchers study the impact of genealogical inheritance factors and validate their work with data results. In the case of subjects in the asymptomatic status, they only consider information known from their affected relatives and present their results in a clear and concise manner as to allow for the implementation of the resulting rules in a non-biased manner. On the contrary, in Conceição et al. (18), we have a few concerns.

Although we know in different mutations and in many regions, there is a trend in favor of early or late onset with correspondence on the predominant neurological or cardiac phenotype, we also know that there is a major variability in the age of onset and the clinical presentation, even when we consider such subgroups of the disease.

The authors of Conceição et al. (18) indicate that their approach depends on the particular mutation, the typical age of onset for that mutation, and the age of onset in family members with ATTR amyloidosis, without elaborating as to in what manner should the age of onset of family members be used. By this, we mean that authors did not specify if a patient with known family will have a PADO value that gives the average AOO of patients in his/her family tree or will have a PADO that corresponds to the same value, or the result of an anticipation factor when considering the proband and/or immediate family members. These are only a few of the possible considerations that can be made regarding PADO and why it is unclear when it comes to the usage of family data.

We would also like to refer that the authors of Cisneros-Barroso et al. (29), while assessing the hypothesis of an anticipation factor in the AOO in a cohort of Mallorca patients, examine PADO. They state that while in a previous international consensus on ATTRV30M amyloidosis authors recommended to start monitoring asymptomatic carriers 10 years before PADO, their findings suggest that this should be done with caution, specifically when patients come from endemic areas where ATTRv amyloidosis is common.

In summary, when comparing Lemos et al. (15) and Conceição et al. (18), it is clear that while in the first study, authors reflect on the fact that there is an evolution related with the AOO of patients when their family data is taken into consideration, and that despite of the obvious experience of Conceição et al. (18) authors' in treating and following patients, their study lacks a clear mathematical description of the rules that correspond to the V30M variant, which allows for their translation to computational operations. As such, the study of Lemos et al. (15) represents a more suitable baseline (BL) to compare with our study.

Regarding our results for the prediction questions, we have the following issues. Our results in the classification problem show the supremacy of the SV algorithm even though they are impacted by the rather small prediction period. As for the regression problem, our results show the supremacy of the machine learning approaches when compared to the baseline results: in particular, SV and XG are ranked as top performers. In specific, XG results tend to be the most stable, when comparing training/ testing and validation results, namely, with the MAE rank values. Indeed, although SV ranked first in the validation set, its inconsistency between testing and validation sets leads us to rely mainly on XG. In the case of XG, its good position was already expected in part due to the algorithm's strengths and the intrinsic characteristics of the problem of predicting the age of onset of ATTRv amyloidosis patients, with clinical and genealogical data, namely, the existence of missing values, complex data, and small sized available inputs. It is also important to reference that these results might well be influenced by XG built-in capabilities to prevent overfitting while working well with small data, as well as with data with subgroups. The differences in size between the training/ testing and validation set can impact the generalization of each algorithm and can lead to differences between training/ testing and validation results (see Table 1).

Since some of the tried methods showed some instability, namely, in later ages (see Figure 3), from a cautionary perspective, we also suggest the usage of DT, as a more intelligent, but still human perceptible, baseline. In practice, this method represents a simpler approach that, regardless of age group and data set, can function as a computational baseline.

To the best of our knowledge, this is the first work that studies and presents a data-oriented approach that can report the risk of a patient showing symptoms before they appear. We base our results and conclusions on data from a single reference center but with the largest number of patients and families with the same mutation and genetic background known worldwide. The registers we analyzed go back many decades and have been collected and sanitized according to similar criteria. Although this is an important work for guiding the clinical follow-up of Portuguese families, from a cautionary perspective, it would be important to verify its replicability in other disease foci, namely, with different mutations and clinical characteristics such as France, Spain, Sweden, and Japan (30). By this, we mean that even though it may be difficult to apply our conclusions to other centers and groups of patients, we think this study can be a first step in defining the much-needed guidelines for a data-oriented prediction of the age of onset of a carrier of a TTR mutation.

Overall, with this study, we successfully achieved our goals as we defined a predictive model capable to be well-understood by medical professionals since in the top performance, we have tree-based algorithms; compared with the current medical practice, namely, the baseline approach (BL) and successfully showed the improvement with good predictive results. We also defined an adequate experimental approach that enables us to explore different sets of genealogical features and compare their results in a time-to-event prediction setting.

Regarding future work directions, we intend to (i) work on the definition of a modeling strategy to predict the age of onset independently of patients' current age (ensemble guided approach); (ii) work on a modeling strategy that uses genotypic information of the patients, i.e., the genetic mutation of each patient; (iii) study the age limits for which to develop statistical validated models; (iv) attempt to accurately and amply model the costs of asymptomatic patients as well as missed treatments for symptomatic patients; (v) study different approaches to extract valuable information from genealogical data sets; and (vi) explore the definition of another medical baseline that follows the inferences of Conceição et al. (18), although in this case, we would need to integrate information from several different and geographically disperse medical centers. Finally, we add that we believe that there is also a need to focus on studying explainability approaches, such as LIME, in order to understand a few of the irregular trends found in our results and supply medical professionals with a set of workable results for their daily patient assessment.

Data availability statement

Patient data is not publicly available due to information that could compromise research participant privacy/consent. Requests should be directed to MP, bWFyaWEuai5wZWRyb3RvJiN4MDAwNDA7aW5lc2N0ZWMucHQ= or TC, dGNvZWxobyYjeDAwMDQwO25ldGNhYm8ucHQ=.

Ethics statement

The studies involving human participants were reviewed and approved by Centro Hospitalar Universitário de Santo António, Porto, Portugal. The patients/participants provided their written informed consent to participate in this study.

Author contributions

MP, AJ, and JM-M: conceptualization and methodology. MP: software and data curation. MP, TC, AJ, and JM-M: writing—original draft preparation and writing—review and editing. TC, AJ, and JM-M: supervision and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This study was financed by National Funds through the Portuguese Funding Agency, FCT—Fundação para a Ciência e a Tecnologia within project UID/ EEA/ 50014/ 2020.

Conflict of interest

TC received support from Pfizer, Ionis, Alnylam, and Biogen to attend to scientific meetings.

The remaining 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.

References

1. Benson MD, Buxbaum JN, Eisenberg DS, Merlini G, Saraiva MJM, Sekijima Y, et al. Amyloid nomenclature 2020: update and recommendations by the International Society of Amyloidosis (ISA) nomenclature committee. Amyloid. (2020) 27:217–22. doi: 10.1080/13506129.2020.1835263

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Ando Y, Coelho T, Berk JL, Cruz MW, Ericzon BG, Ikeda SI, et al. Guideline of transthyretin-related hereditary amyloidosis for clinicians. Orph J Rare Dis. (2013) 8:31. doi: 10.1186/1750-1172-8-31

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Rapezzi C, Quarta CC, Obici L, Perfetto F, Longhi S, Salvi F, et al. Disease profile and differential diagnosis of hereditary transthyretin-related amyloidosis with exclusively cardiac phenotype: an Italian perspective. Eur Heart J. (2013) 34:ehs123. doi: 10.1093/eurheartj/ehs123

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Alves-Ferreira M, Azevedo A, Coelho T, Santos D, Sequeiros J, Alonso I, et al. Beyond Val30Met transthyretin (TTR): variants associated with age-at-onset in hereditary ATTRv amyloidosis. Amyloid. (2021) 28:100–6. doi: 10.1080/13506129.2020.1857236

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Dardiotis E, Koutsou P, Zamba-Papanicolaou E, Vonta I, Hadjivassiliou M, Hadjigeorgiou G, et al. Complement C1Q polymorphisms modulate onset in familial amyloidotic polyneuropathy TTR Val30Met. J Neurol Sci. (2009) 284:158–62. doi: 10.1016/j.jns.2009.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Steyerberg EW. Clinical Prediction Models: A Practical Approach to Development, Validation, and Updating, 2nd Edn. Cham: Springer (2019). doi: 10.1007/978-3-030-16399-0

CrossRef Full Text | Google Scholar

7. Suvila K, Langén V, Cheng S, Niiranen TJ. Age of hypertension onset: overview of research and how to apply in practice. Tech Rep. (2020) 9:1071. doi: 10.1007/s11906-020-01071-z

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Moore KM, Nicholas J, Grossman M, Boeve BF, Knopman DS, Graff-Radford NR, et al. Age at symptom onset and death and disease duration in genetic frontotemporal dementia: an international retrospective cohort study. Lancet Neurol. (2020) 19:145–56. doi: 10.1016/S1474-4422(19)30394-1

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sousa AMBC. A variabilidade fenotípica da Polineuropatia Amiloidótica Familiar: um estudo de genética quantitativa em Portugal e na Suécia (Phenotypic Variability of Familial Amyloid Polyneuropathy: A Quantitative Genetics Study in Portugal and in Sweden). (Ph. D. Thesis), University of Porto, Porto, Portugal (1995).

Google Scholar

10. Gorram F, Alarcon F, Perdry H, Hébrard B, Damy T, Fanen P, et al. Refine penetrance estimates in the main pathogenic variants of transthyretin hereditary (familial) amyloid polyneuropathy (TTR-FAP) using a new non-parametric approach (NPSE). Amyloid. (2017) 24(Suppl.1):115–6. doi: 10.1080/13506129.2017.1293515

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Waddington-Cruz M, Schmidt H, Botteman MF, Carter JA, Stewart M, Hopps M, et al. Epidemiological and clinical characteristics of symptomatic hereditary transthyretin amyloid polyneuropathy: a global case series. Orph J Rare Dis. (2019) 14:1. doi: 10.1186/s13023-019-1000-1

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Jacobson DR, Alexander AA, Tagoe C, Buxbaum JN. Prevalence of the amyloidogenic transthyretin (TTR) V122I allele in 14–333 African–Americans. Amyloid. (2015) 22:171–4. doi: 10.3109/13506129.2015.1051219

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Reilly MM, Staunton H, Harding AE. Familial amyloid polyneuropathy (TTR ala 60) in north west Ireland: a clinical, genetic, and epidemiological study. J Neurol Neurosurg Psychiatr. (1995) 59:45–9. doi: 10.1136/jnnp.59.1.45

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Koike H, Misu KI, Ikeda S, Ando Y, Nakazato M, Ando E, et al. Type I (transthyretin Met30) familial amyloid polyneuropathy in Japan: early- vs. late-onset form. Archiv Neurol. (2002) 11:1771–6. doi: 10.1001/archneur.59.11.1771

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Lemos C, Coelho T, Alves-Ferreira M, Martins-Da-Silva A, Sequeiros J, Mendonça D, et al. Overcoming artefact: anticipation in 284 Portuguese kindreds with familial amyloid polyneuropathy (FAP) ATTRV30M. J Neurol Neurosurg Psychiatry. (2014) 2013:305383. doi: 10.1136/jnnp-2013-305383

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Obici L, Kuks JB, Buades J, Adams D, Suhr OB, Coelho T, et al. Recommendations for presymptomatic genetic testing and management of individuals at risk for hereditary transthyretin amyloidosis. Curr Opin Neurol. (2016) 29:S27–35. doi: 10.1097/WCO.0000000000000290

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Ando Y, Adams D, Benson MD, Berk JL, Planté-Bordeneuve V, Coelho T, et al. Guidelines and new directions in the therapy and monitoring of ATTRv amyloidosis. Amyloid. (2022) 29:143–55. doi: 10.1080/13506129.2022.2052838

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Conceição I, Damy T, Romero M, Galán L, Attarian S, Luigetti M, et al. Early diagnosis of ATTR amyloidosis through targeted follow-up of identified carriers of TTR gene mutations. Amyloid. (2019) 26:3–9. doi: 10.1080/13506129.2018.1556156

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Planté-Bordeneuve V, Gorram F, Olsson M, Anan I, Mazzeo A, Gentile L, et al. A multicentric study of the disease risks and first manifestations in Hereditary transthyretin amyloidosis (ATTRv): insights for an earlier diagnosis. Amyloid. (2023) 2023:1–8. doi: 10.1080/13506129.2023.2178891

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Pedroto M, Jorge A, Mendes-Moreira J, Coelho T. Predicting age of onset in TTR-FAP patients with genealogical features. In:J Hollmén, C McGregor, P Soda, B Kane, , editors, 31st IEEE International Symposium on Computer-Based Medical Systems,CBMS 2018, Karlstad, Sweden, June 18-21, 2018. Karlstad: IEEE Computer Society (2018). p. 199–204.

Google Scholar

21. Pedroto M, Jorge A, Mendes-Moreira J, Coelho T. Impact of genealogical features in transthyretin familial amyloid polyneuropathy age of onset prediction. In:FR Florentino, M Saberi, A Mohamad, R Miguel, F De Paz Juan, G Pascual, , editors, Practical Applications of Computational Biology and Bioinformatics, 12th International Conference. Cham: Springer International Publishing (2019). p. 35–42.

Google Scholar

22. Petroto M, Jorge A, Mendes-Moreira J, Coelho T. Improving the prediction of age of onset of TTR-FAP patients using graph-embedding features. In:Marreiros G, Martins B, Paiva A, Ribeiro B, Sardinha A, , editors. Progress in Artificial Intelligence. EPIA 2022. Lecture Notes in Computer Science, Vol. 13566. Cham: Springer (2022). doi: 10.1007/978-3-031-16474-3_16

CrossRef Full Text | Google Scholar

23. Bender MA, Farach-Colton M, Pemmasani G, Skiena S, Sumazin P. Lowest common ancestors in trees and directed acyclic graphs 1. J Algorithms. (2005) 57:75–94. doi: 10.1016/j.jalgor.2005.08.001

CrossRef Full Text | Google Scholar

24. Beaulieu D, Berry JD, Paganoni S, Glass JD, Fournier C, Cuerdo J, et al. Development and validation of a machine-learning ALS survival model lacking vital capacity (VC-Free) for use in clinical trials during the COVID-19 pandemic. Amyotr Lat Scler Frontotemp Degener. (2021) 22:22–32. doi: 10.1080/21678421.2021.1924207

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Petersohn C. Temporal Video Segmentation. Vogt (2010).

Google Scholar

26. Friedman JE. Anticipation in hereditary disease: the history of a biomedical concept. Hum Genet. (2011) 130:705–14. doi: 10.1007/s00439-011-1022-9

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Drugge U, Andersson R, Chizari F, Danielsson M, Holmgren G, Sandgren O, et al. Familial amyloidotic polyneuropathy in Sweden: a pedigree analysis. J Med Genet. (1993) 30:388–92.

PubMed Abstract | Google Scholar

28. Tashima K, Ando Y, Tanaka Y, Uchino M, Ando M. Change in the age of onset in patients with familial amyloidotic polyneuropathy type I. Japan Soc Intern Medicine. (1995) 34:748–50.

PubMed Abstract | Google Scholar

29. Cisneros-Barroso E, González-Moreno J, Rodríguez A, Ripoll-Vera T, Álvarez J, Usón M, et al. Anticipation on age at onset in kindreds with hereditary ATTRV30M amyloidosis from the Majorcan cluster. Amyloid. (2020) 2020:254–8. doi: 10.1080/13506129.2020.1789580

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Parman Y, Adams D, Obici L, Galán L, Guergueltcheva V, Suhr OB, et al. Sixty years of transthyretin familial amyloid polyneuropathy (TTR-FAP) in Europe. Curr Opin Neurol. (2016) 2016:288. doi: 10.1097/WCO.0000000000000288

PubMed Abstract | CrossRef Full Text | Google Scholar