- 1Division of Basic Neurosciences, Biomedical Research Foundation, Academy of Athens, Athens, Greece
- 2Univ. Lille, Inserm, CHU Lille, U1172 - Lille Neuroscience & Cognition, Lille, France
- 3Inserm, UMR-S 1172, Team “Brain Biology and Chemistry”, Lille, France
- 4Merck & Co., Kenilworth, NJ, United States
- 5Biogen, Cambridge, MA, United States
- 6Denali Therapeutics Inc., South San Francisco, CA, United States
Evidence is mounting that LRRK2 function, particularly its kinase activity, is elevated in multiple forms of Parkinson’s disease, both idiopathic as well as familial forms linked to mutations in the LRRK2 gene. However, sensitive quantitative markers of LRRK2 activation in clinical samples remain at the early stages of development. There are several measures of LRRK2 activity that could potentially be used in longitudinal studies of disease progression, as inclusion/exclusion criteria for clinical trials, to predict response to therapy, or as markers of target engagement. Among these are levels of LRRK2, phosphorylation of LRRK2 itself, either by other kinases or via auto-phosphorylation, its in vitro kinase activity, or phosphorylation of downstream substrates. This is advantageous on many levels, in that multiple indices of elevated kinase activity clearly strengthen the rationale for targeting this kinase with novel therapeutic candidates, and provide alternate markers of activation in certain tissues or biofluids for which specific measures are not detectable. However, this can also complicate interpretation of findings from different studies using disparate measures. In this review we discuss the current state of LRRK2-focused biomarkers, the advantages and disadvantages of the current pallet of outcome measures, the gaps that need to be addressed, and the priorities that the field has defined.
Introduction
Parkinson’s disease (PD) is a debilitating neurodegenerative disorder, affecting millions of people worldwide. The current therapeutic options address symptoms only and there is no approved therapy that slows progression or modifies disease course. PD is a complex disorder influenced by both genetic and environmental factors. The first unequivocal genetic data supporting susceptibility to PD were mutations found in SNCA (encoding α-synuclein) and the subsequent identification of SNCA gene duplications (Polymeropoulos et al., 1997; Singleton et al., 2003). A few years later, mutations in the leucine-rich repeat kinase 2 (LRRK2) gene were found to exhibit significant impact across familial and sporadic PD (Paisan-Ruiz et al., 2004; Zimprich et al., 2004). Hundreds of nonsense or missense genetic variations have been described in the LRRK2 locus (Ross et al., 2011). However, only a few are considered pathogenic: p.Asn1437His, p.Arg1441Gly, p.Arg1441Cys, p.Arg1441His, p.Arg1441Ser, p.Tyr1699Cys, p.Gly2019Ser, and p.Ile2020Thr; with several other risk factors (e.g., p.Gly2385Arg) or variants of unclear pathogenicity (p.Arg1628Pro and p.Ser1761Arg). Their frequency varies markedly depending on the population founder effects of the G2019S-LRRK2 mutation, reaching 30–42% of PD patients in North African Arabic populations as well as 6–30% in Ashkenazi Jewish populations, probably resulting from a mutation arising at least 5,000 years ago (for review of the genetics of LRRK2, please see Monfrini and Di Fonzo, 2017). The collective data strongly suggesting that each of the different point mutations increase kinase activity (Liu et al., 2018). Genome wide association studies (GWAS) analyses also demonstrated that variants at the LRRK2 locus, such as single nucleotide polymorphisms, are among the most important genetic risk factors for PD (Monfrini and Di Fonzo, 2017). Emerging data also suggests that intergenic LRRK2 variants may be associated with increases in LRRK2 gene expression and accelerated PD motor symptom development (Võsa et al., 2018; Iwaki et al., 2019). In Figure 1, we show a schematic of the LRRK2 domain architecture, highlighting both pathogenic as well as other risk factor or functional variants.
Figure 1. Domain architecture of human LRRK2 protein. A schematic of the known functional domains within the LRRK2 protein. Also indicated are the currently identified pathogenic mutations linked to PD (bold magenta), risk factor variants (italicized blue), protective variants (green); and below the schematic are shown several key phosphorylated residues that are auto-phosphorylation sites (pink) or phosphorylated by other kinases (blue).
LRRK2 plays an important role in vesicular trafficking. It impacts endosomal, lysosomal and autophagosomal pathways (Roosen and Cookson, 2016; Alessi and Sammler, 2018), which are also affected by other well-defined PD genes, such as SNCA and GBA1 (Blandini et al., 2019), strongly implicating these fundamental cellular processes in PD pathophysiology. Recent data from post-mortem PD brain and multiple in vivo models, suggest a role for LRRK2 in idiopathic disease as well (Di Maio et al., 2018). Preclinical studies have shown that genetic knock-out of LRRK2, inhibition of LRRK2 with small molecules, or ASO-mediated knockdown reduce pathology and protect from α-synuclein induced dopaminergic neuronal loss in rodents (Daher et al., 2014, 2015; Zhao et al., 2017), also supporting the hypothesis that, even in the absence of familial mutations, LRRK2 can be pathogenic under certain conditions.
Collectively, human genetic studies and preclinical data have led to biopharma initiating drug discovery efforts that have resulted in 2 potential therapeutics progressing into clinical trials (clinicaltrials.gov/ct2/show/study/NCT03710707; clinicaltrials.gov/ct2/show/NCT03976349). There are three potential strategies for clinical development of these LRRK2 therapeutics. Firstly, trials may selectively include genetically defined LRRK2 mutation carriers that have been diagnosed with PD. This would be dependent on patients knowing their own genetic status or focused screening efforts1. However, limitations in enrolling appropriate numbers of suitable LRRK2 mutation carriers will likely provide a significant hurdle in Phase 2 and Phase 3 trials; as the prevalence of LRRK2 mutations, which are estimated at approximately 5% of all PD cases, vary significantly depending on the geographic location, and the relative frequency of specific mutations, which also varies greatly (Monfrini and Di Fonzo, 2017). If there were to be additional stratification, for example, only including G2019S or R1441C/G, this would further reduce this limited patient pool. A second potential clinical design is a prodromal approach: identifying subjects with LRRK2 mutations and determining if disease onset could be prevented by pre-treatment of the potential therapeutic. The major limitations of this approach are the limited genetic penetrance of LRRK2 (Goldwurm et al., 2011; Lee et al., 2018) and, in the cases where the mutation carriers do progress to disease, the unpredictable age of onset, as well as the absence of safety data in subjects undergoing long-term chronic LRRK2 kinase inhibition; providing significant cost/length challenges and an uncharted regulatory path.
Finally, strengthening the link between LRRK2 and idiopathic PD (iPD) could identify cohorts of patients where LRRK2, in the absence of known pathogenic mutations, is driving disease pathophysiology. In this case, the need for LRRK2 biomarkers, i.e., biological measures related to LRRK2 that can identify PD processes or therapeutic response, is absolutely critical given the heterogeneous nature of PD. As introduced above, there is evidence highlighting a link between α-synuclein and LRRK2 (Daher et al., 2014, 2015; Zhao et al., 2017). Similarly, a link has been postulated between LRRK2 and GCase activity (Alcalay et al., 2015; Nguyen and Krainc, 2018); although some controversy still exists concerning the nature of this link. Given these links as well as the prevalence of LRRK2 risk variants in the sporadic PD population, there is significant evidence supporting the therapeutic potential of LRRK2 inhibition in sporadic PD as well as additional familial PD cohorts beyond LRRK2 mutation carriers.
The clinical development of LRRK2 therapeutics will be strongly dependent on biomarkers, as target engagement and pharmacodynamic endpoints are critical for successful progression of clinical candidates (Morgan et al., 2012). This is particularly vital in PD where efficacy trials are long (average length of current trials is ∼2 years), will require significant numbers of subjects (100+ per arm), and will be costly (in the hundreds of millions USD).
There are several measures of LRRK2 function that could potentially be used in longitudinal studies of disease progression, as inclusion/exclusion criteria for clinical trials, as markers of target engagement, and as markers to predict response to therapy. Among these are levels of LRRK2, phosphorylation of LRRK2, either by other kinases or via auto-phosphorylation, in vitro LRRK2 kinase activity, and phosphorylation of downstream substrates or functional endpoints related to elevated (or therapeutic suppression of) LRRK2 function, which will be covered in this review.
LRRK2 Outcome Measures
In probing the function of LRRK2, with the goal of quantifying changes that coalesce around specific disease-stratifying variables (e.g., disease state, LRRK2 mutation status, etc.), there are a number of biochemical outcome measures that are available. These include the quantification of: total LRRK2 levels; phosphorylated LRRK2 (at multiple residues; including S935, S1292, see Figure 1 and below for more details); in vitro kinase activity using model peptide substrates; phosphorylation of endogenous LRRK2 substrates (e.g., Rab10); and others. The specific methodologies employed for each of these measures are discussed in more detail below (see section “Assays Being Employed”). However, to date, most of the early reports (with a few exceptions, see above) assessing these targets have relied largely on Western immunoblotting, which in comparison to ELISA-based approaches for example, is limited in terms of the quantitative range that is possible, and sensitivity. Each of the measures described reveal a distinct, yet equally important, feature of the activation “state” of LRRK2; and importantly, this pattern may also manifest differently depending on the source of the biospecimen examined. Note that a summary overview of LRRK2 related measures and potential applications is given in Table 1.
Total LRRK2 Levels
Total expression levels of LRRK2, depending on the tissue/cell type, can vary in PD, and thus can potentially be a useful tool to assess activation during the different stages of the disease. For example, in the CNS, LRRK2 protein levels are elevated in prefrontal cortex of PD patients (Cho et al., 2013), while CSF levels were only elevated in G2019S PD, but not in iPD or non-manifesting G2019S carriers (Mabrouk et al., 2020). Outside the CNS, immune cells are an ideal source of LRRK2 since they are obtained non-invasively and previous reports have shown elevated levels in iPD compared to healthy controls (Cook et al., 2017). In that study, levels were determined by a novel flow cytometric approach using a LRRK2 knockout validated antibody (rabbit monoclonal; clone c41-2). Specifically, LRRK2 expression was increased in CD16+ monocytes, as well as B and T cells; and this expression was correlated with both intracellular as well as secreted levels of certain cytokines (Cook et al., 2017). The regulation of LRRK2 expression in specific immune cell sub-types is unclear, however, it is known that specific pro-inflammatory mediators, such as IFN-γ, can induce expression of LRRK2 (Gardet et al., 2010); thus, the increased levels of LRRK2 in specific immune cells may be linked to elevated peripheral inflammation, which may or may not be associated with PD (e.g., Dzamko, 2017). In the earlier study of Dzamko (2017), assessing pS935-LRRK2 levels by Western immunoblotting in isolated PBMCs, no difference in total LRRK2 expression was detected between iPD and healthy control subjects in this mixed cell population. Thus, given that the bulk of LRRK2 expression in blood cells is concentrated in a few cellular sub-types (e.g., see Fan et al., 2018), including neutrophils (which were not specifically assessed in either study), it is possible that changes in LRRK2 levels, like phosphorylation of LRRK2 as discussed above, in specific types of peripheral blood cells are heterogeneous.
Heterologous LRRK2 Phosphorylation
Phosphorylation of LRRK2 at a cluster of serine residues located within the N-terminal region of LRRK2 (e.g., S910, S935, S955, and S973), immediately upstream of the namesake leucine-rich repeat domain, represents an additional biochemical readout of LRRK2. The apparent relative abundance of these post-translational modifications (PTMs) in comparison to other sites, such as S1292, has rendered phosphorylation at these sites more easily detected, however, functional interpretation of these findings is complicated by the fact that these are not auto-phosphorylation modifications, as is the case for pS1292. Multiple kinases have been implicated in the phosphorylation of these residues, including: CK1-α1 (Chia et al., 2014), PKA (Muda et al., 2014), TBK-1 (Dzamko et al., 2012), and others. However; analogous to what is observed for auto-phosphorylation sites, phosphorylation at S935 is sensitive to pharmacological LRRK2 kinase inhibition, such that there is a rapid de-phosphorylation at this site (and the other N-terminal serine residues) following treatment of cells, or in vivo, with specific LRRK2 kinase inhibitors (e.g., Dzamko et al., 2010; Vancraenenbroeck et al., 2014). Interestingly, over-expressed kinase inactive mutant LRRK2 [e.g., D1994A and K1906M/R; (Ito et al., 2014)], does not display S935 dephosphorylation relative to WT, indicating that acute (pharmacological) inhibition of LRRK2 alters this regulatory cycle, while chronic genetic ablation of LRRK2 kinase activity does not. This is explained by the fact that the S935-LRRK2 phosphorylation levels do not correlate to kinase activity but rather to the sensitivity of LRRK2 to phosphatases. Indeed, LRRK2 is sensitized to dephosphorylation by LRRK2 kinase inhibitors and for certain LRRK2 mutants with reduced basal S935-LRRK2 phosphorylation. LRRK2 de-phosphorylation at the S935 cluster is mediated by the catalytic subunit of protein phosphatase 1 that is recruited to the LRRK2 complex in conditions of pharmacological inhibition of the LRRK2 kinase (Lobbestael et al., 2013). Conversely, over-expression of pathogenic mutant forms of LRRK2, such as G2019S or R1441C/G, which are known to enhance the kinase activity of LRRK2, does not show enhanced levels of pS935-LRRK2, and in fact have been reported to have decreased levels of phosphorylation at this site (Nichols et al., 2010; Li et al., 2011), including at endogenous levels in immortalized lymphoblasts from G2019S-LRRK2 mutation carriers (Dzamko et al., 2010).
The earliest report of an assay designed to quantify pS935-LRRK2 at endogenous levels came from the group of Delbroek et al. (2013). Using well-validated antibodies (i.e., in knock-out tissue), this group established a quantitative detection method for S935-LRRK2 phosphorylation that demonstrated loss of signal in kinase inhibitor-treated cells and animals. Additionally, as a proof of concept, phosphorylation at this site was detected in human PBMCs from healthy volunteers, that also was sensitive to LRRK2 kinase inhibition (Delbroek et al., 2013). Apart from the ELISA-based detection of this PTM of LRRK2, several studies employing Western immunoblotting have also been reported. The same year as the report from Delbroek et al. (2013), the group of Dzamko et al. (2013) examined pS910-LRRK2 and pS935-LRRK2 levels in PBMCs from healthy controls or iPD patients. While a significant correlation between both phosphorylated sites and total LRRK2 was found in these cells, there was no significant change in pS935 or pS910 levels, when normalized to total LRRK2, in the iPD group (Dzamko et al., 2013). The authors in this study correctly pointed out that in a mixed cellular population of PBMCs, where LRRK2 expression is concentrated in few cell sub-types (Fan et al., 2018), changes in phosphorylation of LRRK2 in distinct cellular types may not be uniform. Additional clinical studies assessing these changes in LRRK2 within specific purified cell types (e.g., monocytes, neutrophils, etc.) are necessary to determine if pS935-LRRK2 levels are detectable in selective cellular populations. In a study of a small cohort of iPD and G2019S mutation carriers, levels of pS935-LRRK2, normalized to total LRRK2 by Western immunoblotting, showed a non-significant decrease in comparison to healthy controls in isolated neutrophils (Fan et al., 2018).
LRRK2 Autophosphorylation
The phosphorylation status of LRRK2 is reflective of its activation in a number of distinct ways. First, and most directly, auto-phosphorylation of LRRK2, for example at the S1292 site, is indicative of its kinase activity in the cell of origin at the time of collection. A number of factors can come into play to determine the level of phosphorylation at this, or other auto-phosphorylation site(s), not just the level of kinase activity of LRRK2 alone. For example, the presence and activity of relevant phosphatases, the sub-cellular localization of LRRK2, the activity of upstream regulators of LRRK2, the status of the ROC GTPase domain, and even the cell type, can all influence the final degree of pS1292 observed. Finally, phosphorylation at S1292 has also been detected in EVs present in CSF (Wang S. et al., 2017), at significantly higher levels, with the signal saturated in many samples, in comparison to pS1292-LRRK2 present in urinary EVs. This saturation effect in the Western immunoblot detection of pS1292-LRRK2 from CSF EVs prevented the stratification of LRRK2 G2019S carriers from non-carriers. This limitation would likely be overcome using an ELISA-based approach (with suitable antibodies for pS1292-LRRK2) in which the usable linear range of detection is typically much broader.
Intrinsic LRRK2 Kinase Activity
Finally, in addition to the cellular indices of LRRK2 kinase activity (LRRK2 phosphorylation, phosphorylation of endogenous substrates such as Rab GTPases), the intrinsic kinase activity of isolated LRRK2 can also be informative. In this approach, LRRK2 is purified from a specific biosample (e.g., PBMCs), and an in vitro kinase reaction is performed using as a substrate model peptides such as LRRKtide or the related NICtide. There are several key differences between assessing kinase activity in this way (i.e., the in vitro activity of the purified enzyme), vs. assessing kinase activation by determining auto-phosphorylation (e.g., pS1292-LRRK2) or phosphorylation of endogenous cellular substrates (e.g., pT73-Rab10). First, performing an in vitro kinase reaction will allow the determination of any changes in the intrinsic activity of purified LRRK2. For example, it is possible that certain PTMs that are known to affect LRRK2 (e.g., phosphorylation and ubiquitination) can alter the intrinsic activity of the kinase domain. If such PTMs are more prevalent in the diseased state, compared to healthy control subjects, the functional consequence of these may be altered kinase activity. This alteration can be detected in an in vitro assay, in an un-encumbered way, without the potential influence of interacting proteins (depending on the stringency of the purification conditions). Secondly, “cellular” assays (measuring phosphorylation of LRRK2 or its substrates) provide a “snap-shot” of kinase activation that is the result of the coordinated action of myriad upstream and downstream regulatory factors, interacting proteins, and the general activation state of the cell. We have employed such an assay, initially in over-expression models (e.g., Leandrou et al., 2019), but more recently in a clinical study assessing LRRK2 in peripheral blood cells (Melachroinou et al., 2020). In this approach, LRRK2 is purified in an ELISA plate, capturing the protein with anti-LRRK2 (or epitope tag) antibodies, followed by an in-well kinase reaction in which the reaction mixture containing the peptide substrate is added directly to the well containing immobilized LRRK2. Possible evolutions of this in vitro kinase activity approach could be to include measures of autophosphorylation (such as the pS1292-LRRK2 measure, as described in Melachroinou et al., 2016) or of Rab substrate phosphorylation (by spiking in recombinant Rab substrate proteins rather than peptide substrates).
Substrates of LRRK2 Kinase
Phosphorylation of LRRK2 substrates represents another potentially informative outcome measure of LRRK2 kinase activation; and like several of the other markers discussed, can also be dependent upon cell or tissue source. In 2016, in a landmark study from the groups of Alessi and Mann, several members of the Rab GTPase family were identified as endogenous kinase substrates of LRRK2 (Steger et al., 2016). A conserved residue within the switch II domain of these GTPases was found to be robustly phosphorylated both in cellular systems as well as in vivo. Several phospho-specific antibodies to certain Rab proteins have since been developed and characterized (Lis et al., 2018), and are now being deployed in studies of LRRK2 activation in clinical samples and as potential markers of target engagement. In another study, using in-house developed phospho-specific antibodies to pT73-Rab10 and pS106-Rab12, Thirstrup et al. (2017) demonstrated that a novel inhibitor of LRRK2 kinase activity could reduce phospho-Rab levels in stimulated PBMCs, but only after 24 h of treatment (Thirstrup et al., 2017). Likewise, similar to other reports, kinase inhibition, at 24 h, also reduced LRRK2 levels in comparison to un-treated cells. The goal in this study, as samples from PD cohorts were not examined, was to assess the utility of assessing pRab10 and pRab12 rates (i.e., phosphorylated Rab as a proportion of total Rab expression) as a marker of target engagement, and as a proof of concept study, this was indeed demonstrated. The principal caveat associated with this study is that LRRK2 levels were artificially induced in isolated PBMCs, following culture for 3 days in the presence of PMA and IFN-γ (Thirstrup et al., 2017). Later evidence demonstrated the translatability of both pS935-LRRK2 and pT73-Rab10 as pharmacodynamic readouts in the clinical setting in unstimulated PBMCs. In human subjects treated with the LRRK2 inhibitor DNL201 for 1–10 days, both readouts showed a robust exposure-dependent reduction in PBMCs (Denali Therapeutics Inc., MJFF PD Therapeutics Conference 2018).
Two studies, thus far, have examined Rab10 phosphorylation (pT73) in peripheral blood cells of PD patients, both with and without the G2019S LRRK2 mutation. A first study by Fan et al. (2018) showed the feasibility of using Rab10 phosphorylation measures by demonstrating good detection levels of Rab10 and pT73-Rab10 in neutrophils and beginning to show increases in Rab10 phosphorylation in small samples of idiopathic PD or PD with LRRK2-G2019S compared to healthy controls. In a larger study, comprised of almost 50 subjects from control or iPD groups, Rab10 phosphorylation in isolated neutrophils or PBMCs was assessed. Consistent with the earlier report, LRRK2 inhibitor treatment significantly reduced pS935 levels as well as pRab10-T73 in both neutrophils and mixed PBMC cellular populations (Atashrazm et al., 2018), with no difference in the degree of response between control and iPD subjects. Interestingly, similar to the report of increased LRRK2 expression in B or T cells, or CD16+ monocytes (Cook et al., 2017), levels of LRRK2 in purified neutrophils (but not PBMCs) are also elevated (Atashrazm et al., 2018). Additionally, neither cell type revealed differences in phosphorylation of Rab between iPD and control subjects. Taken together, while phosphorylation of Rab10, by Western immunoblotting, appears to be a suitable marker for target engagement in clinical studies of LRRK2 inhibition, it remains unclear whether this readout can reliably stratify subjects according to patient group; as thus far, differences between control and iPD or LRRK2 mutation carriers have not been observed. It should be noted, however, that for the LRRK2 mutation carrier study, the sample size and statistical power was low (intentionally, for a proof-of-concept study) precluding the possibility to reach significant conclusions. Further analyses in larger cohorts, ideally with more quantitative approaches, are clearly warranted.
Current Assays Being Employed
An important aspect of the evaluation of LRRK2 and related targets as potential biomarkers of PD is to have a good understanding of the assays used. We present here the assay methods that have been used in recent literature (e.g., ELISA), or are at earlier developmental stages (e.g., PET tracer ligands) to measure LRRK2 status.
Western Immunoblots to Measure LRRK2 Function
Western blots targeting pS935-LRRK2, pS1292-LRRK2, and pT73-Rab10 have been successfully used pre-clinically to detect and measure total levels of LRRK2, activation of LRRK2 kinase, and LRRK2 function. As a pharmacodynamic endpoint reflecting LRRK2 inhibition, the phospho-specific LRRK2 and Rab10 targets are well established pre-clinically. Measurement of pS935-LRRK2 showed a rapid reduction in S935 phosphorylation following LRRK2 inhibition in cellular models (Dzamko et al., 2010) and in vivo pharmacokinetics/dynamics studies (Delbroek et al., 2013; Fell et al., 2015; Fuji et al., 2015), enabling quantification of LRRK2 inhibitor potency in cells and tissues where LRRK2 is endogenously expressed. Similarly, both pS1292-LRRK2 (in HEK cells overexpressing a mutant form of LRRK2) and pT73-Rab10 (in mouse tissues and HEK cells overexpressing Rab10 and LRRK2) are dose-dependently reduced following LRRK2 inhibition as measured by Western blot (Sheng et al., 2012; Atashrazm et al., 2018; Fan et al., 2018; Lis et al., 2018). pS935-LRRK2 and pT73-Rab10 are also measurable by Western blot and reduced following LRRK2 inhibition ex vivo in PBMCs and neutrophils, demonstrating the potential translatability of these markers for clinical use (Perera et al., 2016; Atashrazm et al., 2018; Fan et al., 2018; Lis et al., 2018; and see below in section “Tissue/Biofluid Origin”).
In addition to use as pharmacodynamic readouts, several studies have measured pS1292-LRRK2 and phosphorylated Rab proteins in PD patient samples to test the hypothesis that LRRK2 kinase activity is elevated in all or a subset of PD (Dzamko et al., 2013; Fraser et al., 2016a; Atashrazm et al., 2018; Fan et al., 2018; Lis et al., 2018). pS1292-LRRK2 has not been reproducibly detectable in accessible blood matrices, while the results with pT73-Rab10 have not conclusively demonstrated elevated LRRK2 kinase activity in PD patient samples. Therefore, at this point, pT73-Rab10 (as well as pS935-LRRK2) are more likely to be useful as pharmacodynamic markers than patient selection markers, although additional studies using more sensitive and high throughput assays in additional matrices are ongoing that may change the landscape on this point.
Despite the successes of assessing LRRK2 function via Western blot in many preclinical studies; Western blot has strong disadvantages as a potential biomarker endpoint in the context of a clinical trial. In order to enable clear interpretation and quantitative analysis with rapid turnaround time, clinical assays for pharmacodynamic readouts or patient selection must be highly quantitative, ideally allowing for absolute measurement of the analyte of interest, robust to implement in different locations or over extended periods of time, and relatively high throughput. Western blots are semi-quantitative, differ greatly from user to user, and generally allow for analysis of <100 samples at a time. Therefore, new methods of LRRK2 and Rab measurement must be developed to maximize the utility of these biomarkers in the clinical setting.
ELISA to Measure LRRK2 Function
ELISAs offer a more sensitive and high-throughput method to interrogate LRRK2 kinase activity and pharmacodynamics of LRRK2 kinase inhibitors. Thus far, three assays have been published (Delbroek et al., 2013; Henderson et al., 2015; Scott et al., 2017), each utilizing a sandwich-ELISA approach by capture with a total LRRK2 antibody, followed by detection with a specific pS935-LRRK2 antibody. The latter two studies have facilitated accurate IC50 measurements for LRRK2 kinase inhibitors from treated mouse tissues (brain and kidney) (Henderson et al., 2015), LRRK2 G2019S SH-SY5Y cell lysates (Scott et al., 2017), and human PBMC lysates (Padmanabhan et al., 2020). Notably, Meso Scale Discovery (MSD) has made a pS935-LRRK2 sandwich-ELISA-based assay commercially available, alongside a comparable assay that measures total LRRK2 protein levels. Use of both assays facilitates normalization of pS935-LRRK2 to total LRRK2 levels to account for compound effects on LRRK2 expression or half-life, in addition to standard normalization to tissue weight or total protein levels. These ELISAs offer enhanced sensitivity compared to Western blots, as low as 400 picomolar, as well as the option for high-throughput 384-well assay design.
An emerging alternative is the SIMOA platform offered by Quanterix, which applies digital ELISA technology. The SIMOA platform utilizes a bead-based approach to enable single molecule labeling detected by fluorescence. In addition, the SIMOA assay is able to use sample volumes often <5 μL and run up to 400 samples per shift. A recent MJFF-led study (Padmanabhan et al., 2020) utilizing the SIMOA platform, reported levels as low as 19 pg/mL for total LRRK2 and 4.2 pg/mL for pS935-LRRK2 using full-length recombinant human LRRK2; and subsequently applied this approach to human PBMC lysates. Altogether, ELISA-based assays offer a more sensitive, high-throughput alternative to Western blotting with multiplex potential for measurement of pS935-LRRK2 biomarker levels. It is crucial to note, however, that it will be vital to compare each approach, across platforms and in different centers, with parallel samples to determine if similar estimations of LRRK2 concentration and phosphorylation are obtained by the various assays.
Quantification of reduced pS935-LRRK2 by conventional ELISA and SIMOA assays can accurately reflect pharmacodynamic response following administration of LRRK2 kinase inhibitors and therefore is currently used as a surrogate biomarker, even though pS935-LRRK2 is not a direct measurement of kinase activity (see above). In fact, it has been reported that the ratio of pS935-LRRK2 to total LRRK2 is significantly reduced in human PBMC lysates from PD manifesting LRRK2 G2019S carriers compared to iPD samples and healthy controls [with and without G2019S mutations (Padmanabhan et al., 2020)], although an alternative in-house developed ELISA detected a slight but significant increase in pS935-LRRK2 in PBMCs of iPD, compared to healthy controls (Melachroinou et al., 2020). Measurement of the auto-phosphorylation site pS1292-LRRK2 would be a more ideal marker of LRRK2 kinase activity, but reliance on this biomarker has been hindered by low physiological stoichiometry (Sheng et al., 2012), and limited phospho-specific antibodies. As newer clones of antibodies targeting this site, and pT73-Rab10 as well (see above), are validated for use in more quantitative and sensitive methods such as ELISA, these challenges will likely be overcome. Nonetheless, Di Maio et al. (2018) recently reported a method using proximity ligation to amplify pS1292-LRRK2 immunostaining in the substantia nigra of human iPD tissue, which was increased compared to healthy controls. These exciting data suggest LRRK2 kinase inhibitors may have broader therapeutic potential for the larger PD patient population, beyond those carrying mutations in the LRRK2 gene. There is potential for proximity ligation technology to be converted to more high-throughput qPCR-based platforms, though this has not yet been reported for pS1292-LRRK2. Additional improvements on the quality of reagents available for pS1292-LRRK2 detection will likely enable better utilization of this site as a biomarker. Similar quantification strategies for LRRK2-mediated phosphorylation of Rab substrates, particularly of pT73-Rab10, may also offer additional alternatives for more direct markers of LRRK2 kinase activity in the future.
Mass Spectrometry to Measure LRRK2 Levels and Function
Liquid chromatography-mass spectrometry (LC-MS) has wide ranging applications from exploratory to regulated clinical use, yet quantitative measurements of very low abundance proteins remains challenging due to sensitivity limits and artifacts such as matrix-induced ion suppression. By and large, protein measurements by LC-MS utilize “bottom up” proteomics techniques whereby proteins are digested by proteases into smaller peptides, which are then analyzed for their signature parent and fragment ion mass-to-charge (m/z) ratios. Peptides between 10 and 20 amino acids are in the ideal range for specificity (i.e., not likely to exist in different protein types) and sensitivity (i.e., are more likely to perform better in electrospray ionization-MS). Trypsin, which cleaves proteins at the c-terminus of arginine (R) and lysine (K), is the most commonly used protease in this context. In some instances, trypsin does not yield an appropriate peptide when a specific amino acid sequence is desired. For example, a recent article by Wang S. et al. (2017) showed detection of total LRRK2 and pS1292-LRRK2 by LC-MS using the Glu-C protease since trypsin would not generate a viable peptide containing S1292 - the S1292 site is flanked by K residues (KLSK), thus trypsin would generate a 3 amino acid peptide (LSK). A peptide this short would not necessarily only come from LRRK2 and so assay specificity would be lost. The group instead chose to use the less common protease Glu-C, which cleaves at the C-term of glutamic and aspartic acid residues and this process generated a 14 amino acid peptide between E1287 and E1301 (MGKLSKIWDLPLDE) containing S1292. The mass spectrometer can then distinguish and quantify the un-phosphorylated and phosphorylated peptide species. The group then showed that phosphorylated rLRRK2 (including pS1292) was reliably detected, however, the article stops short of quantifying pS1292 LRRK2 in biological samples. It is likely that an antibody enrichment step would still be required in biological samples to be successful since no cleanup step was applied.
Although there are a number of sample cleanup steps that can reduce sample complexity (including sample fractionation), these steps can be laborious and can introduce variability. Gaining momentum in the field of protein biomarker quantification is the so-called “hybrid ligand binding assay (LBA)-LC-MS” methodology, whereby proteins are isolated from samples using antibodies (similar to ELISA), followed by protease digestion and LC-MS analysis. This methodology has the advantage of greatly reducing sample complexity and improving MS analysis. When a high-resolution mass spectrometer such as an orbitrap or FT-ICR system is used, specificity of signal is encoded by unique peptides that only exist in the targeted protein. This is an advantage over traditional ELISAs where detection specificity must be demonstrated experimentally by analyzing samples in various matrices and testing KO tissues, for example. To our knowledge there are no reported hybrid-LC-MS assays in the literature being used for routine LRRK2-pLRRK2 quantitation in the context of a fit-for-purpose biomarker assay. As this approach becomes more common, LRRK2-pLRRK2 would be well positioned for this type of assay development because of the availability of several high quality LRRK2 antibodies.
Another variation of the hybrid approach is called SISCAPA (stable isotope standard and capture by anti-peptide antibodies). This technique goes even further in reducing sample complexity. In this approach, samples containing proteins of interest are digested using a protease, and then peptides (not proteins) are isolated using anti-peptide antibodies (Anderson et al., 2004). In principle, following elution from an antibody, samples are purified for a single peptide species. In comparison, anti-protein immunocapture eluent will contain peptides from the entire protein as well as peptides from proteases used. A recent initiative by MJFF sought to develop SISCAPA based assays against regions of LRRK2 that would serve as both total LRRK2 and kinase activity endpoints. Specifically, the MJFF-SISCAPA collaboration developed mouse monoclonal antibodies against linear epitopes containing S935 (HSNSLGPIFDHEDLLK) and S1292 (MGKLSKIWDLPLD) capable of detecting both the native and phosphorylated forms of the peptides. Unfortunately, those results showed only nanogram level sensitivity, which was attributed to the performance of the target peptides on the particular LC/MS platforms used as well as the need for a higher affinity rabbit monoclonal antibody (data not published). As such, the existing assays would have limited sensitivity in the context of human CSF.
Elsewhere in this issue (Mabrouk et al., 2020), we will describe a novel SISCAPA assay using commercially available antibodies that function as anti-peptide antibodies to measure total LRRK2 with sensitivity sufficient for CSF detection.
Development of LRRK2 PET Ligands
Positron emission tomography (PET) is a non-invasive and highly sensitive molecular imaging technique that has multiple applications across the CNS drug discovery field. For example, PET imaging with a radiolabeled molecule can be used to assess that molecule’s biodistribution properties thus allowing for the assessment of brain penetration which otherwise cannot be definitively determined in the clinical setting. PET imaging can also be used to quantify CNS target occupancy by a drug molecule and to confirm CNS target engagement. This is an incredibly powerful tool as it can determine if the hypothesis in question has been sufficiently tested in the clinic [i.e., a proof of concept (PoC) trial outcome was negative, but the CNS target was engaged sufficiently such that it rules out a role for that target in the disease/disease stage]. Finally, PET imaging has the potential to serve as a disease state biomarker if the radiolabeled molecule is specific to a target that is associated with disease or a particular stage of disease.
Given the applications of PET imaging to CNS drug discovery, the identification of a LRRK2 PET ligand could significantly enable the clinical development of LRRK2 kinase inhibitors and has been the subject of intense focus from both industry and academic groups alike. Despite the identification of numerous potent and selective LRRK2 kinase inhibitors, from a variety of structural classes, there are limited reports detailing the successful development of radiolabeled LRRK2 kinase inhibitors. In 2013, Roche/Genentech published a patent in which they described the synthesis of 11C- or 18F-labeled LRRK2 inhibitors, which were related to GNE-1023. Similarly, Wang M. et al. (2017) described the radiolabeling of [11C]-HG-10-102-01 but as with the Roche/Genentech probes, no in vitro or in vivo PET characterization of this molecule was described. Malik et al. (2017) reported that they had successfully radiolabeled [3H]-LRRK2-IN-1, however, its use as a CNS PET tracer is limited by poor off-target selectivity and limited brain penetration of the base molecule. Most recently, Chen et al., 2019 reported on the development of [11C]-GNE-1023 and reported excellent in vitro specific binding of [11C]-GNE-1023 to LRRK2 in rat and NHP brain sections (Chen et al., 2019). However, whole-body ex vivo biodistribution studies exhibited limited brain uptake of [11C]-GNE-1023 in mice despite not being a substrate of the brain efflux transporter Pgp. The authors reported that studies in higher species such as NHP and the development of tracers with improved brain penetration were ongoing.
Additionally, GNE-1023 has been labeled with [18F] rather than [11C], however, minimal specific binding in caudate putamen homogenates from rat, rhesus monkey, and human was reported (Zeng et al., 2018). This group also reported on studies with another radiolabeled LRRK2 kinase inhibitor (compound-B) that is derived from the indazole class and is structurally similar to the highly potent and selective LRRK2 kinase inhibitor MLi-2. While [3H]-compound B showed high binding affinity to LRRK2 WT full-length enzyme (Kd = 57 pM), only modest displaceable and saturable binding of [3H]-compound B was observed in rhesus monkey brain CPu homogenates (Kd = 0.09 nM) (Zeng et al., 2018). Importantly, using either [3H]-compound B or [18F]-GNE-1023, they determined that the Bmax for LRRK2 in the NHP and human brain was very low (∼0.4 nM) and that the resulting tracer binding potentials (Bmax/Kd ratio) were far below the desired Bmax/Kd ratio > 10 which is typically required for the successful development of CNS PET tracers (Patel and Gibson, 2008). In summary, a validated PET ligand for monitoring changes in LRRK2 is not currently available and the probability of success for developing a LRRK2 PET tracer is low, based on observed low Bmax (<1 nM) in the CNS regions of interest.
Tissue/Biofluid Origin
LRRK2 and related measures have already been assessed as biomarker in a wide variety of tissues and biofluids (see Figure 2). Here, we provide an overview of key findings for: blood, urinary exosomes, CSF exosomes, and gut/saliva.
Figure 2. Expression of LRRK2 in multiple tissues/cell types. LRRK2 is widely expressed throughout the body in a variety of cell types and tissues, including high levels of expression in the kidney, lung, and cells of the peripheral immune system; but also in multiple brain regions, the intestine, as well as extracellularly via exosomal release.
Measurement of LRRK2 in Blood and Blood Derivatives
Despite LRRK2 being connected most closely with a disorder of the central nervous system, LRRK2 expression levels are highest in the periphery, in particular white blood cells (Fuji et al., 2015). This enables measurement of LRRK2 markers in blood or cells derived from blood such as PBMCs, a practical and accessible matrix in the context of clinical applications where frequent sampling for pharmacokinetics/dynamics analysis will be required. Indeed, many groups have successfully measured LRRK2 inhibition in PBMCs ex vivo from human samples, and in some cases in vivo in cynomolgus monkeys treated with LRRK2 inhibitors, by quantifying LRRK2pS935 reduction (Delbroek et al., 2013; Fuji et al., 2015; Perera et al., 2016). After the discovery that LRRK2 phosphorylates several Rab GTPases, phospho-specific antibodies targeting the LRRK2-dependent Rabs were developed and used to measure LRRK2 inhibition in PBMCs, in particular pT73-Rab10 (Steger et al., 2016; Fan et al., 2018; Lis et al., 2018). PBMCs are commonly isolated in labs and clinical sites for many applications, and are therefore clearly translatable for the purposes of measuring LRRK2 inhibition in human subjects.
LRRK2 expression varies among the different cell types of cells within PBMCs. It is most highly expressed in neutrophils and monocytes, with lower expression in T cells, B cells, dendritic cells, and natural killer cells (Fuji et al., 2015; Fan et al., 2018). This heterogeneity in LRRK2 expression combined with heterogeneity of cell populations from person to person may add to variability of LRRK2 marker quantification in PBMCs. It has therefore been proposed that isolation of neutrophils and/or monocytes for the purposes of measuring LRRK2 markers may reduce inter- and intra-subject variability, and at least in the case of neutrophils, isolation from many donors and measurement of LRRK2 inhibition either by pS935-LRRK2 or pRab10 measurement has been successfully performed (Fan et al., 2018).
In the clinical setting, it is likely that most centers will have more experience isolating PBMCs, compared to specific sub-types such as neutrophils or monocytes, so the practicalities of cell isolation must be balanced with theoretical gains of isolating a pure and homogeneous cell population. With respect to practicality, the most ideal solution for clinical measurement of LRRK2 inhibition in the periphery would be to measure it in whole blood rather than a population of cells isolated from whole blood. This would make the assay more broadly applicable and practical for clinical sites, however, the ability to track changes in LRRK2 activation within specific cell types will be sacrificed. For this reason, we developed an ELISA-based assay of pS935-LRRK2 and total LRRK2 with sufficient sensitivity for detection of these analytes in whole blood (Denali Therapeutics Inc., MJFF PD Therapeutics Conference 2018). This has indeed resulted in more practical and streamlined sample collection processes for clinical sites, compared to PBMC isolation, that are applicable for multi-center, international studies. Alternatively, at sites with such capabilities, immortalization of lymphocytes might be a useful strategy to identify new biomarkers from one type of cell. For instance, we were able to detect centrosomal cohesion deficits in PBMC derived lymphoblastoid cell lines from LRRK2 G2019S Parkinson’s disease patients, as well as in a subset of sporadic PD patients (Fernandez et al., 2019). This approach, however, is better suited for patient stratification purposes in clinical research studies, as compared to rapid and sensitive markers of target engagement required in a clinical trial.
Thus far, we have only been considering measurement of LRRK2 in blood for the purposes of target engagement, but there has been considerable effort put into measurement of LRRK2 levels and LRRK2 function in blood for the purpose of patient stratification or testing the hypothesis that PD patients without LRRK2 mutations have elevated LRRK2 function that contributes to PD pathogenesis. Total LRRK2, pS935-LRRK2 and pT73-Rab10 have all been measured in PBMCs and in neutrophils in sporadic PD patients, non-PD controls, and LRRK2 carriers with and without PD (Dzamko et al., 2013; Atashrazm et al., 2018; Fan et al., 2018; Melachroinou et al., 2020; Padmanabhan et al., 2020). LRRK2 S935 phosphorylation rates decrease in LRRK2 carriers with PD, while all other groups show no significant differences in levels of the tested analytes (Padmanabhan et al., 2020). In another study, however, pS935-LRRK2 levels were reported to be slightly elevated in iPD patients (Melachroinou et al., 2020). To be fair, for this purpose, one must consider the most relevant cell type in which to measure LRRK2 markers. In particular, in at least one report, specific monocyte sub-types have elevated LRRK2 in PD patients and release inflammatory cytokines to a greater extent in PD patients than in healthy controls following stimulation (Bliederhaeuser et al., 2016; Cook et al., 2017). Given this connection with disease, it is possible that in studies focused on patient stratification, purified monocytes may be the most relevant cell population to examine when developing blood-based markers of increased LRRK2 pathway activity in PD. Thus far, a broad characterization of LRRK2 markers or expression in monocytes in well-powered groups of PD, non-PD, and LRRK2 mutation carriers with or without PD has not been undertaken.
Urine-Derived Exosomes
LRRK2 is present in exosomes, i.e., cell-derived extracellular vesicles (EVs) of 30–100 nm in diameter, in several biofluids including urine ((Fraser et al., 2013) and our own results, Mutez et al., 2016). Proteomics screens of exosomes isolated from urine first indicated the presence of LRRK2 in urinary exosomes (Gonzales et al., 2009). Subsequently, the development of sensitive and specific anti-LRRK2 antibodies allowed the confirmation of the presence of phosphorylated LRRK2 in urinary exosomes.
Semi-quantitative western blot analyses of urinary exosomes have determined that LRRK2 is present in the high pg/ml to low ng/ml range (close to 1,000 pg/ml). Double immunofluorescence labeling of extracellular vesicles with anti-LRRK2 and the exosomal marker TSG101, confirmed the identity of the vesicles containing LRRK2 (Fraser et al., 2013). In light of the gain of toxic kinase function hypothesis in Parkinson’s disease, measures of LRRK2 kinase function are of particular interest, for instance the measure of LRRK2 autophosphorylation sites, including the S1292 site that has robustly been confirmed on endogenous LRRK2 in model systems as well as in human samples. Testing of LRRK2-S1292 phosphorylation in urine has revealed significantly elevated pS1292 levels in subjects harboring the G2019S mutation (Fraser et al., 2016a). This study also reported that for subjects with the G2019S mutation, S1292 phosphorylation is elevated in groups with PD symptoms compared to those without. In a separate study, the same group showed that S1292 phosphorylation is significantly increased in idiopathic PD compared with matched healthy controls (Fraser et al., 2016b). Interestingly, this study also revealed that the severity of cognitive impairment correlates with increased S1292 phosphorylation. Furthermore, a third study by the same lab examined LRRK2 in urinary exosomes compared to CSF exosomes of the same individuals and found that S1292-LRRK2 phosphorylation increases observed in urinary exosomes in subjects harboring the G2019S mutation is reflected by a similar increase in CSF (Wang S. et al., 2017). Interestingly, the study also observes that S1292-LRRK2 phosphorylation is significantly higher in CSF compared to urine in all subjects, suggesting a higher activation level of LRRK2 in brain compared to urine, and highlighting the need for more quantitative measures of LRRK2 functions.
These observations are consistent with the hypothesis that LRRK2 in urinary exosomes is modulated in disease, warranting further study of LRRK2 as a biomarker in this biofluid. For instance, the published results show a partial overlap in the distribution of S1292 phosphorylation levels in urinary exosomes in control and mutant/disease groups, suggesting that it is not an absolute predictor of disease. Also, it remains to be elucidated whether some of the observed differences are specific to certain ethnic groups or are (co-) dependent on additional factors such as dietary habits or sleep patterns, or additional lifestyle factors such as occupation. Weaknesses of this approach include the fact that it is impossible to know the cell type(s) or tissues of origin for the recovered EVs present in urine; however, given the high level of LRRK2 expression in the kidney, it is likely that much of the LRRK2 detected in these samples arises from these cells. Additionally, it is possible that more subtle changes in pS1292-LRRK2 levels may be overlooked due to the reduced sensitivity and quantitative limitations inherent to Western immunoblotting.
Besides measuring LRRK2 and its phosphorylation, other proteins in the LRRK2 complex or LRRK2 pathway offer additional possibilities as LRRK2 related biomarkers of disease or pharmacodynamic response. Therefore LRRK2’s substrates, such as the ezrin–radixin–moesin family of proteins (Jaleel et al., 2007), microtubule affinity-regulating kinase 1 (MARK1) (Krumova et al., 2015), endophilin A (Matta et al., 2012), or Rab proteins (Steger et al., 2016), are also potential biomarkers.
LRRK2 or LRRK2 pathway proteins in urinary exosomes also offer the possibility of monitoring pharmacodynamics response to potential LRRK2 targeting therapeutics. According to this hypothesis, pS1292-LRRK2, pS935-LRRK2 or phospho-Rabs would be reduced in urinary exosomes following LRRK2 inhibitor treatment. This hypothesis remains to be confirmed in biofluids. A caveat to the potential use of this biospecimen source in target engagement measures is that it was initially shown that LRRK2 release in exosomes was sensitive to pharmacological kinase inhibition, specifically via its interaction with 14-3-3 (Fraser et al., 2013). Thus, in samples from subjects undergoing LRRK2 kinase inhibitor treatment, the detection of exosomal LRRK2 will likely be impaired.
It should be noted that these studies also revealed sex differences in LRRK2 levels in urinary exosomes. Most notably, total LRRK2 levels are found to be higher in male compared to female subjects (Fraser et al., 2016b). In addition, pS1292-LRRK2 median levels were higher in men compared to women, while the relative elevation in pS1292-LRRK2 levels for PD versus healthy subjects is greater in women than in men. Interestingly, in a different sample set from a Norwegian patient cohort, sex differences displayed a different trend with males harboring the G2019S mutation showing higher pS1292-LRRK2 levels while the opposite holds true for females (Wang S. et al., 2017).
CSF Exosomes
LRRK2 is not thought to exist as a soluble protein in CSF, which presents a challenge when interrogating its function in the CNS. Despite this obstacle a number of studies have demonstrated LRRK2 detection in CSF after isolating small extracellular vesicles through techniques such as differential ultracentrifugation (e.g., Fraser et al., 2013; Wang S. et al., 2017). For instance, Fraser et al. (2013) showed that in neat CSF, LRRK2 is not detectable, nor in the supernatant of ultracentrifuged CSF, but only in the pellet which contains small EVs (exosomes). Following exosome enrichment, this group has successfully applied western blotting techniques to detect total LRRK2 and pS1292 LRRK2 signals and they continue to study the biological mechanism whereby LRRK2 is introduced into these vesicles. An interesting point is that CSF pLRRK2 does not appear to correlate with urinary pLRRK2 levels and CSF levels did not correlate with disease severity while urinary levels did (Wang S. et al., 2017). It should be noted that the CSF pS1292-LRRK2 levels became saturated (within the semi-quantitative linear range of the Western immunoblot approach) compared to urinary exosomes, complicating the analyses of potential correlations with clinical features.
In terms of having a reliable biomarker endpoint that can be used in a clinical trial, exosome enrichment poses several challenges. Differential ultracentrifugation may be difficult to perform in a reproducible manner across different labs and volume requirements are quite high (∼1 ml). In addition, Western blotting analysis techniques are not considered amenable to the throughput and robustness requirements of a clinical trial. Therefore, additional techniques which can isolate LRRK2 in CSF without exosome enrichment/isolation (see Mabrouk et al., 2020) would be beneficial going forward.
Future Directions and Approaches
Nucleic Acid-Based Approaches
Genome wide association studies analyses revealed that LRRK2 polymorphisms are not only associated with PD, but also other disorders including Crohn’s disease and Leprosy pointing out the importance of the immune functions of LRRK2. Thus, one may expect that LRRK2 genotype stratification might help to better classify patients with higher risk to develop prominent immune phenotypes to orientate clinical trials and pharmacogenomics studies.
Genome Wide Methylation
Assuming that environmental factors may have a larger impact on sporadic PD development compared to familial PD, it is surprising that no difference is found between sporadic PD and LRRK2 patients heterozygous for a LRRK2 mutation either in the methylation status of islands of the LRRK2 promoter in patient derived leucocytes (Fernandez-Santiago et al., 2015) or when investigating whole genome methylation of dopaminergic neurons generated from patient derived induced pluripotent stem cells (iPSCs). Several interpretations might be formulated to explain this result. The role of genetics and environmental factors is proposed to explain the reduced penetrance of LRRK2. It is thus possible that patients with or without LRRK2 mutations share a similar influence of environmental factors or that these unknown environmental factors influenced numerous low risk alleles in genes converging to LRRK2 pathways. Moreover, this same study also revealed an important PD associated hypermethylation occurring only upon differentiation into dopaminergic neurons of PD patients, but not somatic cells (Fernandez-Santiago et al., 2015). This shows that the epigenetic control of the differentiation into dopaminergic neurons plays a crucial role for the development of PD phenotypes. They also highlight the need for exploring the transcriptome expression profiles of sporadic PD and LRRK2 patients to identify biomarkers and other pathways of interest to help better understand the pathogenesis of PD.
LRRK2 RNA Expression and Splicing
The LRRK2 gene on chromosome 12p12 is composed of 51 exons. Usually, large genes are more likely to give rise to several transcripts due to alternative splicing events. The Ensembl database showing only one transcript encoding the full-length protein of 2527 AA is supported by strong biological evidence. Other transcripts encoding proteins of 1271 AA, 454 AA, 521 AA, 206 AA, or 78 AA, as well as 3 transcripts not encoding proteins have been proposed based on computational mapping, based on gene identifiers from Ensembl, Ensembl Genomes and model organism databases. With the development of new sequencing technologies, such as RNAseq, several groups have investigated the existence of LRRK2 RNA expression and/or splicing variants in the brain and other tissues. Of interest, association of quantitative trait locus (QTL) involving exons 32–33 have been found in the brain and is associated with the presence of a polymorphism rs3761863 (p.M2397T, involved in Crohn’s disease) together with two additional QTLs in liver and monocytes. Nevertheless, a 2019 study by Vlachakis et al. (2018) recently confirmed the existence of several spliced transcripts in brain occurring at different ratios according the studied brain regions. The development of large transcript sequencing technologies such as PacBio will enable a more robust mapping and reconstitution of each LRRK2 transcript structure. Such analyses have the potential to identify a specific transcript whose expression may be used as an early biomarker of PD and that might then be easily detectable in PBMCs.
Transcriptome Analyses of LRRK2 Patients
Because of the nature of such studies (assessing changes at the transcriptional level rather than the protein level), transcriptomic analyses, in the context of PD biomarker development, are restricted to patient stratification in studies of disease severity and/or progression. These kinds of studies are not applicable to clinical trials of investigational compounds in which markers of target engagement are required. The transcriptome of blood or neurons heterozygous for LRRK2 variants has revealed numerous pathways, similar to idiopathic PD, that differ from controls. In PBMCs and dopaminergic neurons, we found a prevalent common dysregulation of translation, immune system signaling, and vesicular trafficking and endocytosis (Mutez et al., 2014). These results are sustained by other observations showing that LRRK2 controls several steps of these key mechanisms, such as the phosphorylation of several proteins of the translation machinery, the eukaryotic initiation factors 4EBP and ribosomal protein S15 (within drosophila models), and thus deregulating translation (Martin et al., 2014). However, the exact mechanism leading to the deregulation of translation remains poorly understood.
Transcriptome analyses have also highlighted deregulation of intracellular vesicle trafficking and function within the endocytic pathways. Their biological relevance has been confirmed (see above), since we know that LRRK2 phosphorylates at least 10 Rab GTPases regulating such processes as vesicular trafficking and endocytosis. The recent observations of Connor-Robson et al. (2019) study, using both transcriptome and proteome analyses, demonstrated that 25 of the 70 Rabs are deregulated, confirming a major role of LRRK2 in endocytosis.
RNAseq and microarray analyses of both PBMC and iPSC derived dopaminergic neurons have also demonstrated the strong deregulation of the “axon guidance pathway.” Ensemble of Gene Set Enrichment Analyses (EGSEA) of the integrated dataset revealed endocytosis and axon guidance as the two most significantly perturbed pathways, both of which were predicted to be inhibited in the presence of the G2019S mutation. The LRRK2-G2019S mutation has previously been demonstrated to disrupt axon guidance in iPSC-derived dopaminergic neurons (Sanchez-Danes et al., 2012; Reinhardt et al., 2013; Su and Qi, 2013; Borgs et al., 2016). Numerous reports using animal models confirmed deregulation of axonal guidance proteins [for review (Civiero et al., 2018)]. In addition, another analysis using single-cell transcriptional profiles of LRRK2 multipotent neural stem cells revealed neuronal lineages with signature similar to PD. Of note, among these genes, two regulate neurite extension upon down-regulation (NRSN1) or overexpression (SRRM4) (Ohnishi et al., 2017). The authors suggest that it could explain the discrepancies in the results obtained on neurite outgrowth assaying the LRRK2 role in neurite extension (Garcia-Miralles et al., 2015).
Interestingly, deregulation of transcripts linked to mitochondrial dysfunction is also observed by the above studies (Ohnishi et al., 2017). For instance, a significant up-regulation of nine mitochondrial genes was noted, emphasizing the critical role of mitochondria in the disease process. Additionally, the role of LRRK2 mutations in mitochondrial dysfunction is also reported in other PD patient-specific human neuroepithelial stem cells. Aberrations in mitochondrial morphology and functionality were evident in neurons bearing the LRRK2-G2019S mutation compared with isogenic controls (Walter et al., 2019).
Since deregulation of these pathways was also observed in blood cells, further study needs to be performed to establish whether some of these changes might be useful as PD biomarkers, giving clues to the development of novel neuroprotective therapeutics. In this context, Infante et al. (2016) compared the transcriptome of carriers of the LRRK2 G2019S mutation (symptomatic and asymptomatic) as well as PD patients without the G2019S mutation and controls. These comparisons highlighted six deregulated genes that were previously associated to PD risk in Genome-wide association studies (Do et al., 2011; Rhodes et al., 2011; Pankratz et al., 2012; Nalls et al., 2014). Among the 58 genes deregulated in both idiopathic PD and LRRK2 patients, those involved in oxygen transport function or iron metabolism were significantly enriched as we previously noted (Mutez et al., 2011; Mutez et al., 2014). Cell adhesion molecule perturbations were also noted in these latter studies. The deregulation of the extracellular matrix (ECM) was also noted in transcriptome profiles of iPSC derived midbrain-patterned astrocytes from PD patients harboring the LRRK2 G2019S missense mutation (Connor-Robson et al., 2019). These data put forward the involvement of Transforming growth factor beta 1 (TGFB1), an inhibitor of microglial inflammatory processes in murine models of PD (Chen et al., 2017), and matrix metallopeptidase 2 (MMP2), known to degrade α-synuclein aggregates (Oh et al., 2017).
Lipidomics
Another potential alternative LRRK2 related biomarker is BMP [bis(monoacylglycero)phosphate], also known as lysobisphosphatidic acid (LBPA), which is an anionic phospholipid found exclusively on the intra-lumenal vesicles of late endosomes and lysosomes (Bissig and Gruenberg, 2013). BMP can also be secreted into biofluids, where it may be enriched on exosomes (Miranda et al., 2018) or apolipoprotein particles, like HDL (Grabner et al., 2019). BMP promotes electrostatic interactions between intralumenal vesicles and lysosomal lipases and their regulators (e.g., saposins) in order to facilitate glycosphingolipid degradation (Gallala and Sandhoff, 2011). BMP di22:6 levels increase dramatically in urine from patients with the lysosomal storage disorder Niemann-Pick type C, highlighting the translatability of this biomarker as an indicator of changes in lysosome function in vivo (Liu et al., 2014).
Several reports have now firmly demonstrated that LRRK2 activity modulates BMP levels in urine, providing key foundational evidence linking LRRK2 to lysosome function. Cynomolgus monkeys treated with LRRK2 inhibitors GNE-7915 and GNE-0877 for 7 and 29 days showed a dose dependent decrease in urine BMP after 29 days of dosing. This effect was recapitulated in LRRK2 KO mouse urine, demonstrating that the effect is on-target (Fuji et al., 2015). Similarly, the recent study by the group of Alcalay et al. (2019) showed that LRRK2 carriers had elevated urinary di-BMP levels, suggesting a link between LRRK2 and lysosomal function. While BMP reduction in urine represents on-target pharmacology of LRRK2 inhibitors, much work remains to fully understand the dynamics and biological significance of this biomarker. Studies of BMP reductions in urine following LRRK2 inhibitor treatment have focused on time points of maximal inhibition or on long-term recovery time points, so we do not have a good understanding of the timecourse or dose-dependence of BMP reduction relative to other measures of LRRK2 inhibition such as pS935 or pRab10 (Fuji et al., 2015). Additionally, the mechanism by which LRRK2 LOF leads to changes in species of BMP on a cellular level is currently unknown, confounding our understanding of the biological effects of changes in BMP in biofluids. Studies such as this give investigators new directions in understanding LRRK2 biology but also serve as potential biomarkers in clinical trials. Future lipidomic studies examining the relationship between LRRK2, GBA and lysosomal function will help define common mechanisms of genetic PD. In addition, other LRRK2 interactors have been discovered which may have value as biomarkers of LRRK2 function such as 14-3-3 (Nichols et al., 2010).
Alternate Sample Types: Gut and Saliva
Besides improvements in detection or exploitation of additional markers in the LRRK2 pathway, additional avenues can be opened by studying alternate sample types. Besides urine, PBMCs or CSF, other types of human samples may be of interest to monitor LRRK2 or LRRK2 pathway proteins as disease or pharmacodynamic biomarkers. Interestingly, the presence of LRRK2 has been confirmed in both the enteric nervous system (ENS) as well as in the epithelial gut cells. In the enteric nervous system, Maekawa et al. (2017) report LRRK2 expression in the myenteric plexus of the small intestine. These may be of interest in relation to the gut-brain hypothesis of PD pathology whereby the GI tract is considered a trigger site of PD pathological processes (e.g., Santos et al., 2019). In relation to this hypothesis, alpha-synuclein positive structures can be found in neurons of the submucosal plexus of sporadic PD patients and these structures are similar in LRRK2-G2019S PD subjects (Rouaud et al., 2017). Further work will be required to establish whether LRRK2 expression in the ENS is limited to the myenteric plexus or whether LRRK2 is also expressed in the submucosal plexus or in enteric glial cells (Derkinderen, 2017). It also remains to be determined whether LRRK2 may contribute to α-synuclein pathology in the ENS and/or to the transmission of pathological α-synuclein species from the ENS to the CNS.
An additional link of LRRK2 with the gut is the expression of LRRK2 in epithelial gut cells, including Paneth cells (Zhang et al., 2015). This pattern of expression may be put into relation with the finding that genetic association studies have found LRRK2 to be a risk factor for inflammatory bowel disease (Crohn’s disease, CD) (e.g., Derkinderen and Neunlist, 2018; Hui et al., 2018; Ridler, 2018). Studies of LRRK2 KO mice have shown that LRRK2 in Paneth cells is involved in the lysosome sorting process to protect from enteric infection, pointing to a potential pathological mechanism for Crohn’s disease involving LRRK2 in Paneth cells (Rocha et al., 2015). It is also possible that LRRK2 gut expression may affect digestive tract symptoms that are very common in PD such as constipation. From the few studies focusing on non-motor symptoms, the frequency of such GI complications is similar between LRRK2-PD and iPD (e.g., Gaig et al., 2014). It remains to be elucidated whether levels of LRRK2, phospho-LRRK2 or the LRRK2 pathway proteins are affected in CD or PD at the level of the GI tract. A practical consideration here is the invasiveness of collecting gut samples for diagnostic purposes. The procedure performed via endoscopy is considered moderately invasive and is used on a routine basis to diagnose digestive disorders such as colorectal cancer, inflammatory bowel disease and peptic ulcer, therefore its application for Parkinson’s disease is feasible (Corbille et al., 2016).
Saliva is also considered a valuable biofluid for biomarker analysis and has specifically been highlighted for its potential for PD biomarkers. Indeed saliva is an attractive biofluid for diagnostics, especially in the elderly as it is much less invasive than other sample types. Principally, the primary use of saliva is as a source of DNA for genetic testing. Despite the growing interest of saliva as a biomarker fluid, little has been done to analyze LRRK2 protein or LRRK2 pathway proteins in saliva. Recently, proteomics analyses have uncovered that LRRK2 is detected in saliva as one of more than 2,000 confidently identified proteins (Pappa et al., 2018). Further research should now be performed to develop robust and quantifiable detection methods for LRRK2 in saliva and assess LRRK2 and phospho-LRRK2 levels in patient groups compared to controls.
Conclusion
As we have outlined in the sections above, there are a great many options already available for the interrogation of LRRK2 and LRRK2-related pathways as tools in the clinical setting. We have summarized the current state of biomarker development and use in Table 1, and key outstanding issues are highlighted in Box 1. For example, as LC-MS instrumentation manufacturers continue to make gains in terms of sensitivity, ease of use, robustness, more discoveries will be made leading to novel biomarkers to advance clinical stage programs. Mass spectrometry will continue to play an important role both in LRRK2 biomarker discovery and LRRK2 clinical development. These techniques have also been used to identify novel phosphorylation sites on LRRK2 protein (Greggio et al., 2009). From an exploratory perspective, the evolution and adoption of LC-MS techniques has proven to be extremely powerful with the discovery of the Rab proteins as bona fide substrates of LRRK2 kinase activity (Steger et al., 2016), and in just a few short years, Rab10/pRab10 measurements have been introduced as a clinical endpoint in a LRRK2 therapeutic trial.
Box 1. Outstanding issues.
1. In general, studies reporting differences in specific biomarkers in PD patient groups compared to healthy control are still few in number. It remains therefore important to verify whether initial findings can be broadly replicated and extended to longitudinal studies.
2. Biomarker readouts have often been assessed individually, however, it is unclear whether a single biomarker will have sufficient predictive power. One potential path to resolve this issue is to develop a scoring system that would allow researchers to combine several biomarker readouts and thereby enhance predictive power.
3. Assays used to assess biomarker potential of LRRK2 and LRRK2 related measures have often been low-throughput assays in research laboratories (e.g., Western immunoblotting). For the most promising biomarkers, there remains a need for higher throughput robust assays that can be deployed broadly in clinical laboratories.
4. Our understanding of LRRK2 pathways has increased considerably in the last half decade. Besides kinase substrates that have begun to be considered, several other partners in these pathways remain to be assessed as potential PD biomarkers.
5. Similarly, LRRK2 phosphorylation has been intensely studied for a limited number of phosphosites (particularly S935 and S1292), however, it remains to be assessed what added value other less studied sites may have as PD biomarkers.
6. Besides potential LRRK2 related biomarkers that have emerged from proteomics and phosphoproteomic studies, other omics studies including lipidomics, transcriptomics have begun to point to potential additional potential biomarkers that require further assessment.
Finally, as it should be clear from the literature reviewed here, while the field has made great advances in the use of LRRK2-targeted biomarkers as measures of target engagement (i.e., for small molecule inhibitors of LRRK2 kinase), much work remains in optimizing the interpretation of these outcome measures for use in staging the disease, tracking progression, predicting pheno-conversion (in carriers of specific mutations), or as a tool to confirm the diagnosis of PD. For this aspect to be developed, larger multi-cohort longitudinal studies will be required, assessing multiple readouts for the presence of correlations with specific clinical features, at various disease stages.
Author Contributions
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Conflict of Interest
MF and CL are employees of Merck & Co. WH and OM are employees of Biogen. SH-R is an employee of Denali Therapeutics Inc.
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.
The handling editor declared a past co-authorship with one of the authors HR.
Footnotes
References
Alcalay, R. N., Hsieh, F., Tengstrand, E., Padmanabhan, S., Baptista, M., Kehoe, C., et al. (2019). Higher Urine bis(Monoacylglycerol)phosphate levels in LRRK2 G2019S mutation carriers: implications for therapeutic development. Mov. Disord. 35, 134–141. doi: 10.1002/mds.27818
Alcalay, R. N., Levy, O. A., Waters, C. C., Fahn, S., Ford, B., Kuo, S. H., et al. (2015). Glucocerebrosidase activity in Parkinson’s disease with and without GBA mutations. Brain 138, 2648–2658.
Anderson, N. L., Anderson, N. G., Haines, L. R., Hardie, D. B., Olafson, R. W., and Pearson, T. W. (2004). Mass spectrometric quantitation of peptides and proteins using stable isotope standards and capture by Anti-peptide antibodies (SISCAPA). J. Proteome Res. 3, 235–244. doi: 10.1021/pr034086h
Atashrazm, F., Hammond, D., Perera, G., Bolliger, M. F., Matar, E., and Halliday, G. M. (2018). LRRK2-mediated Rab10 phosphorylation in immune cells from Parkinson’s disease patients. Mov. Disord. 34, 406–415. doi: 10.1002/mds.27601
Bissig, C., and Gruenberg, J. (2013). Lipid sorting and multivesicular endosome biogenesis. Cold Spring Harb. Perspect. Biol. 5:a016816. doi: 10.1101/cshperspect.a016816
Blandini, F., Cilia, R., Cerri, S., Pezzoli, G., Schapira, A. H. V., Mullin, S., et al. (2019). Glucocerebrosidase mutations and synucleinopathies: toward a model of precision medicine. Mov. Disord. 34, 9–21. doi: 10.1002/mds.27583
Bliederhaeuser, C., Zondler, L., Grozdanov, V., Ruf, W. P., Brenner, D., Melrose, H. L., et al. (2016). LRRK2 contributes to monocyte dysregulation in Parkinson’s disease. Acta Neuropathol. Commun. 4:123.
Borgs, L., Peyre, E., Alix, P., Hanon, K., Grobarczyk, B., Godin, J. D., et al. (2016). Dopaminergic neurons differentiating from LRRK2 G2019S induced pluripotent stem cells show early neuritic branching defects. Sci. Rep. 6:33377.
Chen, G., Liu, J., Jiang, L., Ran, X., He, D., Li, Y., et al. (2017). Galangin reduces the loss of dopaminergic neurons in an LPS-evoked model of Parkinson’s disease in rats. Int. J. Mol. Sci. 19:12. doi: 10.3390/ijms19010012
Chen, Z., Shao, T., Gao, W., Fu, H., Collier, T. L., Rong, J., et al. (2019). Synthesis and preliminary evaluation of [(11) C]GNE-1023 as a potent PET probe for imaging leucine-rich repeat kinase 2 (LRRK2) in Parkinson’s disease. Chemmedchem 14, 1580–1585. doi: 10.1002/cmdc.201900321
Chia, R., Haddock, S., Beilina, A., Rudenko, I. N., Mamais, A., Kaganovich, A., et al. (2014). Phosphorylation of LRRK2 by casein kinase 1alpha regulates trans-Golgi clustering via differential interaction with ARHGEF7. Nat. Commun. 5:5827.
Cho, H. J., Liu, G., Jin, S. M., Parisiadou, L., Xie, C., Yu, J., et al. (2013). MicroRNA-205 regulates the expression of Parkinson’s disease-related leucine-rich repeat kinase 2 protein. Hum. Mol. Genet. 22, 608–620. doi: 10.1093/hmg/dds470
Civiero, L., Cogo, S., Biosa, A., and Greggio, E. (2018). The role of LRRK2 in cytoskeletal dynamics. Biochem. Soc. Trans. 46, 1653–1663. doi: 10.1042/bst20180469
Connor-Robson, N., Booth, H., Martin, J. G., Gao, B., Li, K., Doig, N., et al. (2019). An integrated transcriptomics and proteomics analysis reveals functional endocytic dysregulation caused by mutations in LRRK2. Neurobiol. Dis. 127, 512–526. doi: 10.1016/j.nbd.2019.04.005
Cook, D. A., Kannarkat, G. T., Cintron, A. F., Butkovich, L. M., Fraser, K. B., Chang, J., et al. (2017). LRRK2 levels in immune cells are increased in Parkinson’s disease. NPJ Parkinsons. Dis. 3:11.
Corbille, A. G., Clairembault, T., Coron, E., Leclair-Visonneau, L., Preterre, C., Neunlist, M., et al. (2016). What a gastrointestinal biopsy can tell us about Parkinson’s disease? Neurogastroenterol. Motil. 28, 966–974. doi: 10.1111/nmo.12797
Daher, J. P., Abdelmotilib, H. A., Hu, X., Volpicelli-Daley, L. A., Moehle, M. S., Fraser, K. B., et al. (2015). Leucine-rich repeat Kinase 2 (LRRK2) pharmacological inhibition abates alpha-synuclein gene-induced neurodegeneration. J. Biol. Chem. 290, 19433–19444. doi: 10.1074/jbc.m115.660001
Daher, J. P., Volpicelli-Daley, L. A., Blackburn, J. P., Moehle, M. S., and West, A. B. (2014). Abrogation of alpha-synuclein-mediated dopaminergic neurodegeneration in LRRK2-deficient rats. Proc. Natl. Acad. Sci. U.S.A. 111, 9289–9294. doi: 10.1073/pnas.1403215111
Delbroek, L., Van Kolen, K., Steegmans, L., Da Cunha, R., Mandemakers, W., Daneels, G., et al. (2013). Development of an enzyme-linked immunosorbent assay for detection of cellular and in vivo LRRK2 S935 phosphorylation. J. Pharm. Biomed. Anal. 76, 49–58. doi: 10.1016/j.jpba.2012.12.002
Derkinderen, P. (2017). LRRK2 expression in the enteric nervous system: ENSuring its significance. Dig. Dis. Sci. 62, 826–827. doi: 10.1007/s10620-017-4500-7
Derkinderen, P., and Neunlist, M. (2018). Crohn’s and Parkinson disease: is LRRK2 lurking around the corner? Nat. Rev. Gastroenterol. Hepatol. 15, 330–331. doi: 10.1038/s41575-018-0006-9
Di Maio, R., Hoffman, E. K., Rocha, E. M., Keeney, M. T., Sanders, L. H., De Miranda, B. R., et al. (2018). LRRK2 activation in idiopathic Parkinson’s disease. Sci. Transl. Med. 10:eaar5429.
Do, C. B., Tung, J. Y., Dorfman, E., Kiefer, A. K., Drabant, E. M., Francke, U., et al. (2011). Web-based genome-wide association study identifies two novel loci and a substantial genetic component for Parkinson’s disease. PLoS Genet. 7:e1002141. doi: 10.1371/journal.pgen.1002141
Dzamko, N., Chua, G., Ranola, M., Rowe, D. B., and Halliday, G. M. (2013). Measurement of LRRK2 and Ser910/935 phosphorylated LRRK2 in peripheral blood mononuclear cells from idiopathic Parkinson’s disease patients. J. Parkinsons. Dis. 3, 145–152. doi: 10.3233/jpd-130174
Dzamko, N., Deak, M., Hentati, F., Reith, A. D., Prescott, A. R., Alessi, D. R., et al. (2010). Inhibition of LRRK2 kinase activity leads to dephosphorylation of Ser(910)/Ser(935), disruption of 14-3-3 binding and altered cytoplasmic localization. Biochem. J. 430, 405–413. doi: 10.1042/bj20100784
Dzamko, N., Inesta-Vaquera, F., Zhang, J., Xie, C., Cai, H., Arthur, S., et al. (2012). The IkappaB kinase family phosphorylates the Parkinson’s disease kinase LRRK2 at Ser935 and Ser910 during Toll-like receptor signaling. PLoS One 7:e39132. doi: 10.1371/journal.pgen.39132
Dzamko, N. L. (2017). LRRK2 and the immune system. Adv. Neurobiol. 14, 123–143. doi: 10.1007/978-3-319-49969-7_7
Fan, Y., Howden, A. J. M., Sarhan, A. R., Lis, P., Ito, G., Martinez, T. N., et al. (2018). Interrogating Parkinson’s disease LRRK2 kinase pathway activity by assessing Rab10 phosphorylation in human neutrophils. Biochem. J. 475, 23–44. doi: 10.1042/bcj20170803
Fell, M. J., Mirescu, C., Basu, K., Cheewatrakoolpong, B., Demong, D. E., Ellis, J. M., et al. (2015). MLi-2, a potent, selective, and centrally active compound for exploring the therapeutic potential and safety of LRRK2 kinase inhibition. J. Pharmacol. Exp. Ther. 355, 397–409. doi: 10.1124/jpet.115.227587
Fernandez, B., Lara Ordonez, A. J., Fdez, E., Mutez, E., Comptdaer, T., Leghay, C., et al. (2019). Centrosomal cohesion deficits as cellular biomarker in lymphoblastoid cell lines from LRRK2 Parkinson’s disease patients. Biochem. J. 476, 2797–2813. doi: 10.1042/bcj20190315
Fernandez-Santiago, R., Carballo-Carbajal, I., Castellano, G., Torrent, R., Richaud, Y., Sanchez-Danes, A., et al. (2015). Aberrant epigenome in iPSC-derived dopaminergic neurons from Parkinson’s disease patients. EMBO Mol. Med. 7, 1529–1546. doi: 10.15252/emmm.201505439
Fraser, K. B., Moehle, M. S., Alcalay, R. N., West, A. B., and Consortium, L. C. (2016a). Urinary LRRK2 phosphorylation predicts parkinsonian phenotypes in G2019S LRRK2 carriers. Neurology 86, 994–999. doi: 10.1212/wnl.0000000000002436
Fraser, K. B., Rawlins, A. B., Clark, R. G., Alcalay, R. N., Standaert, D. G., Liu, N., et al. (2016b). Ser(P)-1292 LRRK2 in urinary exosomes is elevated in idiopathic Parkinson’s disease. Mov. Disord. 31, 1543–1550. doi: 10.1002/mds.26686
Fraser, K. B., Moehle, M. S., Daher, J. P., Webber, P. J., Williams, J. Y., Stewart, C. A., et al. (2013). LRRK2 secretion in exosomes is regulated by 14-3-3. Hum. Mol. Genet. 22, 4988–5000. doi: 10.1093/hmg/ddt346
Fuji, R. N., Flagella, M., Baca, M., Baptista, M. A., Brodbeck, J., Chan, B. K., et al. (2015). Effect of selective LRRK2 kinase inhibition on nonhuman primate lung. Sci. Transl. Med. 7:273ra215.
Gaig, C., Vilas, D., Infante, J., Sierra, M., Garcia-Gorostiaga, I., Buongiorno, M., et al. (2014). Nonmotor symptoms in LRRK2 G2019S associated Parkinson’s disease. PLoS One 9:e108982. doi: 10.1371/journal.pgen.108982
Gallala, H. D., and Sandhoff, K. (2011). Biological function of the cellular lipid BMP-BMP as a key activator for cholesterol sorting and membrane digestion. Neurochem. Res. 36, 1594–1600. doi: 10.1007/s11064-010-0337-6
Garcia-Miralles, M., Coomaraswamy, J., Habig, K., Herzig, M. C., Funk, N., Gillardon, F., et al. (2015). No dopamine cell loss or changes in cytoskeleton function in transgenic mice expressing physiological levels of wild type or G2019S mutant LRRK2 and in human fibroblasts. PLoS One 10:e0118947. doi: 10.1371/journal.pgen.0118947
Gardet, A., Benita, Y., Li, C., Sands, B. E., Ballester, I., Stevens, C., et al. (2010). LRRK2 is involved in the IFN-gamma response and host response to pathogens. J. Immunol. 185, 5577–5585. doi: 10.4049/jimmunol.1000548
Goldwurm, S., Tunesi, S., Tesei, S., Zini, M., Sironi, F., Primignani, P., et al. (2011). Kin-cohort analysis of LRRK2-G2019S penetrance in Parkinson’s disease. Mov. Disord. 26, 2144–2145. doi: 10.1002/mds.23807
Gonzales, P. A., Pisitkun, T., Hoffert, J. D., Tchapyjnikov, D., Star, R. A., Kleta, R., et al. (2009). Large-scale proteomics and phosphoproteomics of urinary exosomes. J. Am. Soc. Nephrol. 20, 363–379. doi: 10.1681/asn.2008040406
Grabner, G. F., Fawzy, N., Pribasnig, M. A., Trieb, M., Taschler, U., Holzer, M., et al. (2019). Metabolic disease and ABHD6 alter the circulating bis(monoacylglycerol)phosphate profile in mice and humans. J. Lipid Res. 60, 1020–1031. doi: 10.1194/jlr.m093351
Greggio, E., Taymans, J. M., Zhen, E. Y., Ryder, J., Vancraenenbroeck, R., Beilina, A., et al. (2009). The Parkinson’s disease kinase LRRK2 autophosphorylates its GTPase domain at multiple sites. Biochem. Biophys. Res. Commun. 389, 449–454. doi: 10.1016/j.bbrc.2009.08.163
Henderson, J. L., Kormos, B. L., Hayward, M. M., Coffman, K. J., Jasti, J., Kurumbail, R. G., et al. (2015). Discovery and preclinical profiling of 3-[4-(morpholin-4-yl)-7H-pyrrolo[2,3-d]pyrimidin-5-yl]benzonitrile (PF-06447475), a highly potent, selective, brain penetrant, and in vivo active LRRK2 kinase inhibitor. J. Med. Chem. 58, 419–432. doi: 10.1021/jm5014055
Hui, K. Y., Fernandez-Hernandez, H., Hu, J., Schaffner, A., Pankratz, N., Hsu, N. Y., et al. (2018). Functional variants in the LRRK2 gene confer shared effects on risk for Crohn’s disease and Parkinson’s disease. Sci. Transl. Med. 10:eaai7795.
Infante, J., Prieto, C., Sierra, M., Sanchez-Juan, P., Gonzalez-Aramburu, I., Sanchez-Quintana, C., et al. (2016). Comparative blood transcriptome analysis in idiopathic and LRRK2 G2019S-associated Parkinson’s disease. Neurobiol. Aging 38, 211–214.
Ito, G., Fujimoto, T., Kamikawaji, S., Kuwahara, T., and Iwatsubo, T. (2014). Lack of correlation between the kinase activity of LRRK2 harboring kinase-modifying mutations and its phosphorylation at Ser910, 935, and Ser955. PLoS One 9:e97988. doi: 10.1371/journal.pgen.0097988
Iwaki, H., Blauwendraat, C., Leonard, H. L., Liu, G., Maple-Grodem, J., Corvol, J. C., et al. (2019). Genetic risk of Parkinson disease and progression: an analysis of 13 longitudinal cohorts. Neurol. Genet. 5:e348.
Jaleel, M., Nichols, R. J., Deak, M., Campbell, D. G., Gillardon, F., Knebel, A., et al. (2007). LRRK2 phosphorylates moesin at threonine-558: characterization of how Parkinson’s disease mutants affect kinase activity. Biochem. J. 405, 307–317. doi: 10.1042/bj20070209
Krumova, P., Reyniers, L., Meyer, M., Lobbestael, E., Stauffer, D., Gerrits, B., et al. (2015). Chemical genetic approach identifies microtubule affinity-regulating kinase 1 as a leucine-rich repeat kinase 2 substrate. FASEB J. 29, 2980–2992. doi: 10.1096/fj.14-262329
Leandrou, E., Markidi, E., Memou, A., Melachroinou, K., Greggio, E., and Rideout, H. J. (2019). Kinase activity of mutant LRRK2 manifests differently in hetero-dimeric vs. homo-dimeric complexes. Biochem. J. 476, 559–579. doi: 10.1042/bcj20180589
Lee, C. J., Kang, D., Lee, S., Lee, S., Kang, J., and Kim, S. (2018). In silico experiment system for testing hypothesis on gene functions using three condition specific biological networks. Methods 145, 10–15. doi: 10.1016/j.ymeth.2018.05.003
Li, X., Wang, Q. J., Pan, N., Lee, S., Zhao, Y., Chait, B. T., et al. (2011). Phosphorylation-dependent 14-3-3 binding to LRRK2 is impaired by common mutations of familial Parkinson’s disease. PLoS One 6:e17153. doi: 10.1371/journal.pgen.17153
Lis, P., Burel, S., Steger, M., Mann, M., Brown, F., Diez, F., et al. (2018). Development of phospho-specific Rab protein antibodies to monitor in vivo activity of the LRRK2 Parkinson’s disease kinase. Biochem. J. 475, 1–22. doi: 10.1042/bcj20170802
Liu, N., Tengstrand, E. A., Chourb, L., and Hsieh, F. Y. (2014). Di-22:6-bis(monoacylglycerol)phosphate: a clinical biomarker of drug-induced phospholipidosis for drug development and safety assessment. Toxicol. Appl. Pharmacol. 279, 467–476. doi: 10.1016/j.taap.2014.06.014
Liu, Z., Bryant, N., Kumaran, R., Beilina, A., Abeliovich, A., Cookson, M. R., et al. (2018). LRRK2 phosphorylates membrane-bound Rabs and is activated by GTP-bound Rab7L1 to promote recruitment to the trans-Golgi network. Hum. Mol. Genet. 27, 385–395. doi: 10.1093/hmg/ddx410
Lobbestael, E., Zhao, J., Rudenko, I. N., Beylina, A., Gao, F., Wetter, J., et al. (2013). Identification of protein phosphatase 1 as a regulator of the LRRK2 phosphorylation cycle. Biochem. J. 456, 119–128. doi: 10.1042/BJ20121772
Mabrouk, O. S., Chen, S., Edwards, A. L., Yang, M., Hirst, W. D., and Graham, D. L. (2020). Quantitative Measurements of LRRK2 in human cerebrospinal fluid demonstrates increased levels in G2019S patients. Front. Neurosci. 14:526. doi: 10.3389/fnins.2020.00526
Maekawa, T., Shimayama, H., Tsushima, H., Kawakami, F., Kawashima, R., Kubo, M., et al. (2017). LRRK2: an emerging new molecule in the enteric neuronal system that quantitatively regulates neuronal peptides and IgA in the gut. Dig. Dis. Sci. 62, 903–912. doi: 10.1007/s10620-017-4476-3
Malik, N., Gifford, A. N., Sandell, J., Tuchman, D., and Ding, Y. S. (2017). Synthesis and in vitro and in vivo evaluation of [(3)H]LRRK2-IN-1 as a novel radioligand for LRRK2. Mol. Imaging Biol. 19, 837–845. doi: 10.1007/s11307-017-1070-1
Martin, I., Kim, J. W., Lee, B. D., Kang, H. C., Xu, J. C., Jia, H., et al. (2014). Ribosomal protein s15 phosphorylation mediates LRRK2 neurodegeneration in Parkinson’s disease. Cell 157, 472–485. doi: 10.1016/j.cell.2014.01.064
Matta, S., Van Kolen, K., Da Cunha, R., Van Den Bogaart, G., Mandemakers, W., Miskiewicz, K., et al. (2012). LRRK2 controls an EndoA phosphorylation cycle in synaptic endocytosis. Neuron 75, 1008–1021. doi: 10.1016/j.neuron.2012.08.022
Melachroinou, K., Kang, M. S., Liong, C., Narayan, S., Levers, N., and Joshi, N., et al. (2020). Elevated in vitro kinase activity in PBMCs of LRRK2 G2019S carriers: a novel ELISA-based method. Mov. Dis. doi: 10.1002/mds.28175 [Epub ahead of print].
Melachroinou, K., Leandrou, E., Valkimadi, P. E., Memou, A., Hadjigeorgiou, G., Stefanis, L., et al. (2016). Activation of FADD-dependent neuronal death pathways as a predictor of pathogenicity for LRRK2 mutations. PLoS One 11:e0166053. doi: 10.1371/journal.pgen.0166053
Miranda, A. M., Lasiecka, Z. M., Xu, Y., Neufeld, J., Shahriar, S., Simoes, S., et al. (2018). Neuronal lysosomal dysfunction releases exosomes harboring APP C-terminal fragments and unique lipid signatures. Nat. Commun. 9:291.
Monfrini, E., and Di Fonzo, A. (2017). Leucine-rich repeat kinase (LRRK2) genetics and Parkinson’s Disease. Adv. Neurobiol. 14, 3–30. doi: 10.1007/978-3-319-49969-7_1
Morgan, P., Van Der Graaf, P. H., Arrowsmith, J., Feltner, D. E., Drummond, K. S., Wegner, C. D., et al. (2012). Can the flow of medicines be improved? Fundamental pharmacokinetic and pharmacological principles toward improving Phase II survival. Drug Discov. Today 17, 419–424. doi: 10.1016/j.drudis.2011.12.020
Muda, K., Bertinetti, D., Gesellchen, F., Hermann, J. S., Von Zweydorf, F., Geerlof, A., et al. (2014). Parkinson-related LRRK2 mutation R1441C/G/H impairs PKA phosphorylation of LRRK2 and disrupts its interaction with 14-3-3. Proc. Natl. Acad. Sci. U.S.A. 111, E34–E43.
Mutez, E., Larvor, L., Lepretre, F., Mouroux, V., Hamalek, D., Kerckaert, J. P., et al. (2011). Transcriptional profile of Parkinson blood mononuclear cells with LRRK2 mutation. Neurobiol. Aging 32, 1839–1848. doi: 10.1016/j.neurobiolaging.2009.10.016
Mutez, E., Nkiliza, A., Belarbi, K., De Broucker, A., Vanbesien-Mailliot, C., Bleuse, S., et al. (2014). Involvement of the immune system, endocytosis and EIF2 signaling in both genetically determined and sporadic forms of Parkinson’s disease. Neurobiol. Dis. 63, 165–170. doi: 10.1016/j.nbd.2013.11.007
Mutez, E., Sibran, W., Comptdaer, T., Semaille, P., Drouyer, M., Bleuse, S., et al. (2016). Detection of leucine-rich repeat kinase 2 in human biofluids. Mov. Disord. 31:S272.
Nalls, M. A., Pankratz, N., Lill, C. M., Do, C. B., Hernandez, D. G., Saad, M., et al. (2014). Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat. Genet. 46, 989–993.
Nguyen, M., and Krainc, D. (2018). LRRK2 phosphorylation of auxilin mediates synaptic defects in dopaminergic neurons from patients with Parkinson’s disease. Proc. Natl. Acad. Sci. U.S.A. 115, 5576–5581. doi: 10.1073/pnas.1717590115
Nichols, R. J., Dzamko, N., Morrice, N. A., Campbell, D. G., Deak, M., Ordureau, A., et al. (2010). 14-3-3 binding to LRRK2 is disrupted by multiple Parkinson’s disease-associated mutations and regulates cytoplasmic localization. Biochem. J. 430, 393–404. doi: 10.1042/bj20100483
Oh, S. H., Kim, H. N., Park, H. J., Shin, J. Y., Kim, D. Y., and Lee, P. H. (2017). The cleavage effect of mesenchymal stem cell and its derived matrix metalloproteinase-2 on extracellular alpha-synuclein aggregates in parkinsonian models. Stem Cells Transl. Med. 6, 949–961. doi: 10.5966/sctm.2016-0111
Ohnishi, T., Shirane, M., and Nakayama, K. I. (2017). SRRM4-dependent neuron-specific alternative splicing of protrudin transcripts regulates neurite outgrowth. Sci. Rep. 7:41130.
Padmanabhan, S., Lanz, T. A., Gorman, D., Wolfe, M., Joyce, A., Cabrera, C., et al. (2020). An assessment of LRRK2 serine 935 phosphorylation in human peripheral blood mononuclear cells in idiopathic parkinson’s disease and G2019S LRRK2 cohorts. J. Parkinsons. Dis. 10, 623–629. doi: 10.3233/jpd-191786
Paisan-Ruiz, C., Jain, S., Evans, E. W., Gilks, W. P., Simon, J., Van Der Brug, M., et al. (2004). Cloning of the gene containing mutations that cause PARK8-linked Parkinson’s disease. Neuron 44, 595–600. doi: 10.1016/j.neuron.2004.10.023
Pankratz, N., Beecham, G. W., Destefano, A. L., Dawson, T. M., Doheny, K. F., Factor, S. A., et al. (2012). Meta-analysis of Parkinson’s disease: identification of a novel locus, RIT2. Ann. Neurol. 71, 370–384. doi: 10.1002/ana.22687
Pappa, E., Vastardis, H., Mermelekas, G., Gerasimidi-Vazeou, A., Zoidakis, J., and Vougas, K. (2018). Saliva proteomics analysis offers insights on type 1 diabetes pathology in a pediatric population. Front. Physiol. 9:444. doi: 10.3389/fnins.2020.00444
Patel, S., and Gibson, R. (2008). In vivo site-directed radiotracers: a mini-review. Nucl. Med. Biol. 35, 805–815. doi: 10.1016/j.nucmedbio.2008.10.002
Perera, G., Ranola, M., Rowe, D. B., Halliday, G. M., and Dzamko, N. (2016). Inhibitor treatment of peripheral mononuclear cells from Parkinson’s disease patients further validates LRRK2 dephosphorylation as a pharmacodynamic biomarker. Sci. Rep. 6:31391.
Polymeropoulos, M. H., Lavedan, C., Leroy, E., Ide, S. E., Dehejia, A., Dutra, A., et al. (1997). Mutation in the alpha-synuclein gene identified in families with Parkinson’s disease. Science 276, 2045–2047. doi: 10.1126/science.276.5321.2045
Reinhardt, P., Schmid, B., Burbulla, L. F., Schondorf, D. C., Wagner, L., Glatza, M., et al. (2013). Genetic correction of a LRRK2 mutation in human iPSCs links parkinsonian neurodegeneration to ERK-dependent changes in gene expression. Cell Stem Cell 12, 354–367. doi: 10.1016/j.stem.2013.01.008
Rhodes, S. L., Sinsheimer, J. S., Bordelon, Y., Bronstein, J. M., and Ritz, B. (2011). Replication of GWAS associations for GAK and MAPT in Parkinson’s disease. Ann. Hum. Genet. 75, 195–200.
Ridler, C. (2018). Parkinson disease: LRRK2 variants linked to PD and Crohn’s disease. Nat. Rev. Neurol. 14:126. doi: 10.1038/nrneurol.2018.10
Rocha, J. D., Schlossmacher, M. G., and Philpott, D. J. (2015). LRRK2 and Nod2 promote lysozyme sorting in Paneth cells. Nat. Immunol. 16, 898–900. doi: 10.1038/ni.3255
Roosen, D. A., and Cookson, M. R. (2016). LRRK2 at the interface of autophagosomes, endosomes and lysosomes. Mol. Neurodegener. 11:73.
Ross, O. A., Soto-Ortolaza, A. I., Heckman, M. G., Aasly, J. O., Abahuni, N., Annesi, G., et al. (2011). Association of LRRK2 exonic variants with susceptibility to Parkinson’s disease: a case-control study. Lancet Neurol. 10, 898–908.
Rouaud, T., Clairembault, T., Coron, E., Neunlist, M., Anheim, M., and Derkinderen, P. (2017). Enteric alpha-synuclein pathology in LRRK2-G2019S Parkinson’s disease. Parkinsonism. Relat. Disord. 40, 83–84. doi: 10.1016/j.parkreldis.2017.05.001
Sanchez-Danes, A., Richaud-Patin, Y., Carballo-Carbajal, I., Jimenez-Delgado, S., Caig, C., Mora, S., et al. (2012). Disease-specific phenotypes in dopamine neurons from human iPS-based models of genetic and sporadic Parkinson’s disease. EMBO Mol. Med. 4, 380–395. doi: 10.1002/emmm.201200215
Santos, S. F., De Oliveira, H. L., Yamada, E. S., Neves, B. C., and Pereira, A. Jr. (2019). The Gut and Parkinson’s disease-a bidirectional pathway. Front. Neurol. 10:574. doi: 10.3389/fnins.2020.00574
Scott, J. D., Demong, D. E., Greshock, T. J., Basu, K., Dai, X., Harris, J., et al. (2017). Discovery of a 3-(4-Pyrimidinyl) Indazole (MLi-2), an orally available and selective leucine-rich repeat kinase 2 (LRRK2) inhibitor that reduces brain kinase activity. J. Med. Chem. 60, 2983–2992.
Sheng, Z., Zhang, S., Bustos, D., Kleinheinz, T., Le Pichon, C. E., Dominguez, S. L., et al. (2012). Ser1292 autophosphorylation is an indicator of LRRK2 kinase activity and contributes to the cellular effects of PD mutations. Sci. Transl. Med. 4:164ra161. doi: 10.1126/scitranslmed.3004485
Singleton, A. B., Farrer, M., Johnson, J., Singleton, A., Hague, S., Kachergus, J., et al. (2003). alpha-Synuclein locus triplication causes Parkinson’s disease. Science 302:841. doi: 10.1126/science.1090278
Steger, M., Tonelli, F., Ito, G., Davies, P., Trost, M., Vetter, M., et al. (2016). Phosphoproteomics reveals that Parkinson’s disease kinase LRRK2 regulates a subset of Rab GTPases. eLife 5:e12813.
Su, Y. C., and Qi, X. (2013). Inhibition of excessive mitochondrial fission reduced aberrant autophagy and neuronal damage caused by LRRK2 G2019S mutation. Hum. Mol. Genet. 22, 4545–4561. doi: 10.1093/hmg/ddt301
Thirstrup, K., Dachsel, J. C., Oppermann, F. S., Williamson, D. S., Smith, G. P., Fog, K., et al. (2017). Selective LRRK2 kinase inhibition reduces phosphorylation of endogenous Rab10 and Rab12 in human peripheral mononuclear blood cells. Sci. Rep. 7:10300.
Vancraenenbroeck, R., De Raeymaecker, J., Lobbestael, E., Gao, F., De Maeyer, M., Voet, A., et al. (2014). In silico, in vitro and cellular analysis with a kinome-wide inhibitor panel correlates cellular LRRK2 dephosphorylation to inhibitor activity on LRRK2. Front. Mol. Neurosci. 7:51. doi: 10.3389/fnmol.2014.00051
Vlachakis, D., Labrou, N. E., Iliopoulos, C., Hardy, J., Lewis, P. A., Rideout, H., et al. (2018). Insights into the influence of specific splicing events on the structural organization of LRRK2. Int. J. Mol. Sci. 19:2784. doi: 10.3390/ijms19092784
Võsa, U., Claringbould, A., Westra, H.-J., Bonder, M. J., Deelen, P., Zeng, B., et al. (2018). Unraveling the polygenic architecture of complex traits using blood eQTL metaanalysis. bioRxiv [Preprint], doi: 10.1101/447367
Walter, J., Bolognin, S., Antony, P. M. A., Nickels, S. L., Poovathingal, S. K., Salamanca, L., et al. (2019). Neural stem cells of Parkinson’s disease patients exhibit aberrant mitochondrial morphology and functionality. Stem Cell Rep. 12, 878–889. doi: 10.1016/j.stemcr.2019.03.004
Wang, M., Gao, M., Xu, Z., and Zheng, Q. H. (2017). Synthesis of [(11)C]HG-10-102-01 as a new potential PET agent for imaging of LRRK2 enzyme in Parkinson’s disease. Bioorg. Med. Chem. Lett. 27, 1351–1355. doi: 10.1016/j.bmcl.2017.02.019
Wang, S., Liu, Z., Ye, T., Mabrouk, O. S., Maltbie, T., Aasly, J., et al. (2017). Elevated LRRK2 autophosphorylation in brain-derived and peripheral exosomes in LRRK2 mutation carriers. Acta Neuropathol. Commun. 5:86.
Zeng, Z., Miller, P., Graham, T., Tong, L., Hesk, D., Li, W., et al. (2018). Study on feasibility to develop a central nervous system (CNS) leucine-rich repeat kinase 2 (LRRK2) PET tracer [abstract]. Mov. Disord. 33 (suppl 2).
Zhang, Q., Pan, Y., Yan, R., Zeng, B., Wang, H., Zhang, X., et al. (2015). Commensal bacteria direct selective cargo sorting to promote symbiosis. Nat. Immunol. 16, 918–926. doi: 10.1038/ni.3233
Zhao, H. T., John, N., Delic, V., Ikeda-Lee, K., Kim, A., Weihofen, A., et al. (2017). LRRK2 Antisense oligonucleotides ameliorate alpha-synuclein inclusion formation in a Parkinson’s disease mouse model. Mol. Ther. Nucleic Acids 8, 508–519. doi: 10.1016/j.omtn.2017.08.002
Keywords: LRRK2, Rab GTPase, biomarker, Parkinson’s disease, kinase
Citation: Rideout HJ, Chartier-Harlin M-C, Fell MJ, Hirst WD, Huntwork-Rodriguez S, Leyns CEG, Mabrouk OS and Taymans J-M (2020) The Current State-of-the Art of LRRK2-Based Biomarker Assay Development in Parkinson’s Disease. Front. Neurosci. 14:865. doi: 10.3389/fnins.2020.00865
Received: 28 May 2020; Accepted: 24 July 2020;
Published: 18 August 2020.
Edited by:
Patrick Lewis, Royal Veterinary College (RVC), United KingdomReviewed by:
Nicolas Dzamko, The University of Sydney, AustraliaDan Berwick, St George’s, University of London, United Kingdom
Copyright © 2020 Rideout, Chartier-Harlin, Fell, Hirst, Huntwork-Rodriguez, Leyns, Mabrouk and Taymans. 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: Hardy J. Rideout, hrideout@bioacademy.gr; hrideout08@gmail.com; Warren D. Hirst, warren.hirst@biogen.com