Skip to main content

LRRK2 kinase activity is necessary for development and regeneration in Nematostella vectensis

Abstract

Background

The starlet sea anemone, Nematostella vectensis, is an emerging model organism with a high regenerative capacity, which was recently found to possess an orthologue to the human Leucine Rich Repeat Kinase 2 (LRRK2) gene. Mutations in this gene are the most common cause of inherited Parkinson’s Disease (PD), highlighting the importance of understanding its function. Despite two decades of research, however, the function of LRRK2 is not well established.

Methods

To investigate the function of LRRKs in Nematostella vectensis, we applied small molecule inhibitors targeting the kinase activity of LRRK2 to examine its function in development, homeostasis and regeneration in Nematostella vectensis.

Results

In vivo analyses inhibiting the kinase function of this enzyme demonstrated a role of nvLRRK2 in development and regeneration of N. vectensis. These findings implicate a developmental role of LRRK2 in Nematostella, adding to the expanding knowledge of its physiological function.

Conclusions

Our work introduces a new model organism with which to study LRRK biology. We report that LRRK kinase activity is necessary for the development and regeneration of Nematostella. Given the short generation time, genetic trackability and in vivo imaging capabilities, this work introduces Nematostella vectensis as a new model in which to study genes linked to neurodegenerative diseases such as Parkinson’s.

Background

The starlet sea anemone Nematostella vectensis is a member of the Cnidarian family [1, 2], and is an emerging model for studying development and regeneration due to its ease of maintenance in the lab, sequenced genome and genetic tractability [1, 3,4,5]. Cnidaria are considered the sister group to the Bilateria, putting them in an excellent position to study the evolutionary trajectory of gene families. Animals from the Cnidarian phylum possess a strikingly similar gene content to humans compared to better studied ecdysozoan models (such as Drosophila melanogaster) and possess syntenic genetic loci [6,7,8,9]. Of particular interest, N. vectensis possesses 4 LRRK genes, one of which being an orthologue to the human leucine rich repeat kinase two gene (LRRK2) [10].

Being composed of two cellular layers, the ectoderm and the endoderm, Nematostella have a simple body plan and tissue organization [11,12,13,14]. Additionally, Nematostella possesses a nervous system that takes the form of a nerve net, composed of sensory cells, glandular cells, multipolar ganglion cells, and cnidocytes (stinging cells) [4, 15, 16].

Unlike bilaterians, the nervous system of cnidarians is not centralized and represents an ancestral, simple organization [17, 18]. Nonetheless, several conserved neurogenic pathways and proteins mirroring what is seen in bilaterians highlights their potential usefulness for studying the underlying neurobiology of neurodegenerative diseases like Parkinson’s Disease (PD) [1]. For example, a recent study by Steger et al. [8] identified a key upstream regulator of Nematostella neuroglandular lineages (composed of neurons, cnidocytes and gland cells): SoxC, which has similar roles in bilaterians [19,20,21]. Prior studies have also demonstrated homologous transcription factors (such as Hox genes) that are involved in neurogenesis, as well as numerous neuropeptides [22, 23]. It is thought that there are 32 potential neuronal cell types in the Nematostella, emphasizing the unusual complexity of the nervous system in this organism [4, 24, 25].

A key feature of N. vectensis is their remarkable ability to regenerate after injury via cellular proliferation [26,27,28,29,30]. Nematostella can be cut into multiple pieces, and each piece will give rise to a new animal, the exception being that if the foot alone is amputated it cannot regenerate [3, 27, 29]. There has long been an interest in whether regeneration reuses the same genes that were developmentally deployed. Recent work from the Rottinger lab has compared gene expression profiles during early development and regeneration and developed a user-friendly website that allows scientists to view the expression of their genes of interest during development versus regeneration in Nematostella vectensis [3].

The human LRRK2 gene is highly studied due to its contribution to the genetic risk of developing PD, with rare autosomal dominant coding mutations as well as more common non-coding variation at this locus linked to disease risk [31]. The LRRK2 gene, located on human chromosome 12p11.2-q13.1, encodes a 2527 amino acid protein member of the ROCO family of Ras-GTPases [32,33,34]. Coding mutations in human LRRK2 (hereafter huLRRK2), including N1437H, R1441C, Y1699C, G2019S and I2020T, cluster in the enzymatic domains of the protein and increase its kinase activity, with dysregulation of a subset of Rab GTPases a key consequence of this [35]. How this alteration in enzymatic function leads to neurodegeneration is unclear, with studies implicating LRRK2 in a range of cellular processes including lysosomal function, mitochondrial biology and synaptic signaling. The tight association between LRRK2 and Parkinson’s disease has led to a number of in-human clinical trials for both small molecule kinase inhibitors and antisense oligonucleotide gene therapy [36]. Several inhibitors of the LRRK2 kinase domain have been developed, the most studied being MLi-2. MLi-2 binds in the ATP binding pocket of the kinase domain and forms hydrogen bonds with residues R1895 and E1948, thereby blocking the kinase activity [37,38,39].

Intriguingly, recent work has identified a role for LRRK proteins in stem cell proliferation and tissue remodeling in intestinal enterocytes during whole-body regeneration in the bilaterian freshwater planaria, the Schmidtea mediterranea, [40] highlighting this as a potential phenotype of interest Nematostella vectensis.

In this study we have examined the cross-species conservation of LRRK2 and investigated the functional importance of LRRK kinase activity in Nematostella development and regeneration, demonstrating a key role for LRRK2 in these processes.

Methods

Animal care

Nematostella vectensis were maintained at 17–20 °C in Pyrex glass bowls kept in the dark in 15 parts per thousand Instant Ocean [41]. Animals were fed 48-hour old artemia, five times per week. Animals were cleaned a few hours after feeding. Spawning was induced by exposure to light and an increase in temperature to 23–25 °C, and embryos were collected immediately after spawning.

Sequence alignment analysis

NCBI reference sequence proteins for human LRRK2 (NP_940980.4) were compared to Nematostella LRRK2 (XP_048584778.1). The domain structure and boundaries for both proteins were derived from information in the conserved domain databases of NCBI [42]. To compare the similarity between the individual domains of nvLRRK2 and hsLRRK2, the FASTA sequences encoding each domain (based on amino acid positions) were identified and subsequently input into the Needleman-Wunch algorithm via the BLAST database [43,44,45]. The most common mutated residues in human LRRK2 pertaining to PD [46,47,48] and the corresponding residues in nvLRRK2 were identified by sequence alignment. The same approach was used for identifying conservation of reported phosphorylated residues.

nvLRRK2 kinase inhibition

To obtain embryos, N. vectensis were spawned by exposure to increased temperature and light. LRRK2 kinase inhibitors GNE-0877 and MLi-2 were purchased from Tocris Bioscience. The kinase inhibitors were made in a stock of 10mM in DMSO and then diluted in seawater. Concentrations in the range of 1–20 μm were first tested, and the concentrations giving a phenotype and high viability were then used for further experiments. The concentration of the kinase inhibitors used for the embryo experiments was 5 μm, and for the adult experiments, it was 10 μm. Some experiments were carried out on Nv-LWamide transgenic animals that were expressing mCherry protein their neurons [18]. Embryos were added to a 6-well plate and were incubated in a solution of 5 ml sea water and 250 µl of each inhibitor. The control embryos were incubated in seawater. For regeneration experiment wild type or transgenic adult Nematostella were relaxed in 7.4% MgCl2 (ThermoFisher) for 15 min and were subsequently amputated below the pharynx on a plastic petri dish using a sterile no.10 disposable scalpel (World Precision Instruments). The animals were monitored every 2 days and the solutions were changed. Whole animals were fed artemia to ensure any differences were not due to lack of food.

For homeostatic experiments, 6-month-old animals were incubated in 10 μm MLi-2 LRRK2 kinase inhibitor for 1 week, the solution was changed every 2 days. After 1 week the animals were relaxed and imaged on the confocal. The animals were then kept in normal sea water for one week and then reimaged. For these experiments the neuronal transgenic line LWamide-mcherry was used.

Upon completion of this experiment, the organisms were relaxed in 7.4% MgCl2 (ThermoFisher), and then in 10% MgCl2, 15% MgCl2 and 20% MgCl2. The animals were then fixed in 4% paraformaldehyde (PFA) (ThermoFisher) at 4 °C and were subsequently stored at 4 °C until imaged. The transgenic animals were imaged using a Zeiss LSM780 microscope.

qRT-PCR analysis

Embryos were flash frozen on liquid nitrogen, approximately 150 embryos were used for each sample. RNA isolation was carried out following Invitrogen phenol-chloroform RNA extraction protocol. cDNA was then synthesized using iScript™ cDNA Synthesis Kit (BioRad). qRT-PCR reactions were prepared using SsoAdvanced Universal SYBR® Green Supermix (BioRad) and carried out on Real-Time PCR (CFX Opus 96) (BioRad).

Primers used:

18SF: CGG CTT AAT TTG ACT CAA CAC G

18SR: TTA GCA TGC CAG AGT CTC GTT C

Wnt4F: CGC CTA ACT ACT GCC ACA AA

Wnt4R: CCT CGC CCA CAA CAA AGA TA

SoxB1F: GTT GAC GGC TGA AGA GAA GG

SoxB1R: AGA ATT TGT CAA CCG CCA TC

LRRK2F: CCC ATA CCT CAC AGC TAC TTT AC

LRRK2R: CGA ATT CCG CCC TGT GTA TAA

In situ hybridization

A 300 bp fragment of the NvLRRK2 coding region was cloned by PCR. The resulting PCR product was used to synthesize in situ probe by the addition of DIG-labeled UTP (Roche) plus the appropriate RNA Polymerase T7 or Sp6 (NEB). Probes were purified with RNA Clean Up Kit (Qiagen) and resuspended in 100uL of hybridization buffer.

Nematostella were collected at different time points after fertilization, fixed in 4%PFA overnight at 4 C. Embryos or young larval animals were dehydrated in methanol series and stored at -20 C. Before starting the in-situ protocol animals were rehydrated again through a methanol series back into phosphate buffered saline (PBS). In situs were carried out using published protocols [49]. In summary, embryos were incubated in 1:1 PBST: Hyb for 30 min and pre-hybridized for 30 min. NvLRRK2 probe was diluted into hybridization buffer and slides were allowed to hybridize overnight at 55 °C. The following day the animals were washed 3 times in Wash Buffer (50% Formamide, 5x SSC and 0.1% Tween) for 30 min each, once in 1:1 Wash Buffer: PBST for 30 min and once in PBST at 55 °C. Samples were rinsed in room temperature PBST 3 times for 5 min each before blocking buffer was added (2% goat serum, 2% BST in PBSTx) for 1 h. Anti-DIG FAB (Roche) was diluted 1:1000 in blocking buffer and slides were incubated for at least 1 h. Samples were washed 3 times for 10 min each before addition of fresh AP Buffer for at least 10 min. Finally, samples were incubated in BM Purple (Roche) until colored reaction was observed. The reaction was stopped by several quick rinses in PBS and were fixed in 4% PFA for 10 min. Samples were mounted in 80% glycerol and images were taken with a Zeiss Discovery V8 microscope using Zen software.

Statistical analysis

Statistical tests were performed using the Graphpad Prism 9 software. A D’Agostino and Pearson test was carried out to test for normality. Upon finding a non-normal distribution, a non-parametric Kruskal Wallis test was used to test for statistical significance between the medians of each measurement in the different experimental groups. A Dunn’s post-hoc multiple comparison test was used to compare the medians of the experimental groups against the control group.

Results

Nematostella vectenis NvLRRK2 exhibits key residues conservation with human LRRK2

LRRK2 contains four protein-protein interaction domains, as well as domains conferring two distinct enzymatic activities. The kinase domain is a serine-threonine kinase capable of autophosphorylating residues within LRRK2, as well as heterologous range of substrates, most notably a subset of Rab GTPases [50,51,52,53,54,55,56,57,58].

To investigate functional and disease relevant conservation between human and Nematostella LRRK2s, the domain structure of N. vectensis and human LRRK2 were aligned, using publicly available data sets from NCBI [42]. Consistent previous reports and with the overall conservation of LRRK proteins across evolution, these orthologues display a similar domain organization [59]. A Needleman-Wunsch algorithm was used to align the protein sequences to observe any biological similarities or differences [45], revealing conservation of 29% identity and 48% similarity across the open reading frames of the two proteins. To investigate conservation at a domain level, the amino acid sequences of individual domains in human and Nematostella LRRK2 were compared under the same global sequence alignment tool, the results of which are shown in Fig. 1A. The kinase domains were highly conserved, including the key kinase DYGI motif [60]. In contrast, the scaffolding and ROC/COR domains are less well conserved.

Fig. 1
figure 1

(A) Domain structure of nvLRRK2 (NCBI RefSeq: XP_048584778.1) compared to human LRRK2 (hsLRRK2) (NCBI RefSeq: NP_940980.4). (B) Global alignment via Needleman Wunsch algorithm illustrating conservation of mutated residues. Arrow indicates most common LRRK2 mutations in PD, G2019S, and I2020T, are conserved in the nvLRRK2 orthologue. (C) Global alignment via Needleman Wunsch algorithm illustrating conservation of residues in nvLRRK2 that are commonly phosphorylated in hsLRRK2. Arrow indicates S910 and S935 which are not conserved in nvLRRK2, as well as S973 which is conserved

Mutated residues are conserved in N. Vectensis

To assess functional and disease relevant conservation across hsLRRK2 and nvLRRK2, residues that are mutated in human LRRK2 contributing to PD development were identified. These include R1441C/G, and N1437H (localized to the Roc domain), Y1699C (in the COR domain), G2019S and I2020T (in the kinase domain) [61,62,63,64,65,66]. Consistent with the high conservation of the kinase domains between these two proteins, Fig. 1B demonstrates that the most common PD associated LRRK2 mutation, G2019S, as well as the I2020T mutation residing in the same motif, is conserved in nvLRRK2.

In contrast, residues that are mutated in the ROC domain are not conserved.

Residues involved in post-translational modifications are not well conserved

To assess conservation of post-translational regulation between human and Nematostella LRRKs, residues previously reported to be involved in signal transduction were compared. Phosphorylation sites on LRRK2 are involved in regulating function and downstream signaling events [67]. Characterized phosphorylation sites in hsLRRK2 include S910, S935, S973, S955 located between the Ank and LRR domain, S1292 between LRR and Roc, and T2031, S2032 and T2035 in the kinase domain [67].

Unlike the Parkinson’s associated mutated residues, phosphorylated residues were not well conserved in either the LRR domain and the kinase domain. Figure 1C shows the sequence alignment demonstrating non-conserved residues as well as a conserved site. Only the S973 residue in the LRR domain and T2035 in the kinase domain were conserved in nvLRRK2, suggesting that the mechanisms involved in regulating signal transduction have diverged between these two proteins over evolution.

Scaffolding and structural domains in LRRK2 are less well conserved throughout species due to the lack of dependence for the functioning of the protein. Thus, is it unsurprising that these residues were not well conserved in nvLRRK2, and it suggests that post-translational modifications have perhaps drifted throughout evolution to suit higher or lower species depending on their functional requirements. This is consistent with global phosphoproteomic analyses that suggest phosphorylation events evolve individually between different species [68].

NvLRRK2 kinase activity is necessary for faithful embryonic development

Based upon the high conservation of the kinase domains of hsLRRK2 and nsLRRK2, to investigate a role of the kinase activity of nvLrrk2 in regulating the development of N.vectensis, anemones were treated with structurally distinct validated inhibitors of LRRK2 kinase activity GNE-0877 and MLi-2. Freshly laid embryos were incubated in 5µM of respective inhibitor diluted in Nematostella water, while control embryos were grown in just water. The concentration of inhibitor to use was determined empirically by first testing concentrations in the range of 1–20µM, 5µM was chosen as for both inhibitors most of the embryos survived and showed a phenotype, while in higher concentrations a large percentage of the embryos died within a few days. Embryos incubated at 5µM developed through the first phases of development and reached the motile stage at relatively the same frequency as the control embryos. However, by 1 week post fertilization clear differences were observed, usually at this timepoint the embryos have developed tentacles and will commence feeding. Animals where the kinase activity of LRRK2 has been inhibited were overall smaller in length (Fig. 2) and most showed no development of tentacles or stunted growth of tentacles (Fig. 2) in comparison to controls that have much longer bodies and clear tentacles. Overall, the inhibitor-treated animals display stunted inhibited growth from around 4 days post-fertilization onwards making it very difficult to determine if they developed normal organs like the pharynx and mesenteries. However, as all LRRK kinase inhibitor treated animals displayed stunted growth of tentacles or no tentacle growth we next examined expression of genes involved in specifying the oral region of the animal, Wnt4. Wnt4 is expressed early in development from the planula stages and is important for induction of genes involved in specifying the oral region of the embryo [30, 69,70,71,72]. Here we find that at 48 h post fertilization (hpf) in embryos exposed to the LRRK inhibitor that levels of Wnt4 are significantly decreased in comparison to the control embryos, suggesting that the embryos are not inducing the genes necessary to specify head and to ultimately direct the embryos towards making tentacles (Fig. 3A). As the tentacles are formed early in development and are highly innervated, we also examined induction and formation of neurons. NvSoxB is well-characterized to be expressed early in development in the cells that will form neurons and nematocysts [73, 74]. Quantitative analysis of NvSoxB levels in control embryos versus the LRRK2 kinase inhibitor treated animal discovered a significant lack of expression of NvSoxB in inhibitor treated animals (Fig. 3B). We next examined the presence of development in the Nematostella embryos taking advantage of the NvLWamide-like::mCherry transgenic reporter line [18]. By one week post-fertilization young animals have differentiated neurons with a complex network of axons in the body and tentacles as seen in Fig. 3C. In comparison in animals exposed to the LRRK kinase inhibitor far fewer neurons are present, and they fail to extend axons and form networks (Fig. 3D, E) suggesting that LRRK2 kinase activity is necessary for specification and proper differentiation of neurons during embryonic development. To determine if NvLRRK2 kinase activity is also necessary for maintenance of neurons in adult animals we placed 3-month and 6-month adult transgenic animals in LRRK2 kinase inhibitors and imaged their neurons after 4 days. We observed in all cases that the axons appeared to degenerate, and no axonal processes could be observed in the inhibitor treated animals versus control animals (Fig. 2S), suggesting that the kinase activity of LRRK2 is necessary to maintain healthy connected axons. Interestingly we carried out further experiments where the animals were exposed to the kinase inhibitor for one week, neurons appeared to degenerate, then the drug was washed out and the animals allowed to recover for 1 week (Fig. 3F-H). After a 1-week recovery period, axons were observed to regenerate in the kinase-treated animals, the axons did not regenerate to the same density of axons observed in the controls after 1-week (Fig. 3F-H), but the initial regeneration suggests that there is the latent potential to recover from the effects of an inactive LRKK2 kinase activity. After the initial 1-week inhibitor treatment neuronal cell bodies could still be observed in the animals, this suggests that the neurons regenerate the axons.

Fig. 2
figure 2

Pharmacological inhibition of the NvLRRK2 kinase domain causes defects in embryonic development. A Control embryos 12 days post-fertilization, showing normal body length and development of internal organs like mesenteries and pharynx and external tentacles. B-C embryos grown in the presence of the Lrrk2 kinase inhibitor GNE-0877 or Mli-2 exhibit stunted body growth and failure to development proper tencles. Data shows mean + SD of body length in each groups; controls (12dpf n = 120), GNE-0877 (12dpf n = 121), and MLi-2 (12dpf n = 130). Non-parametric post-hoc Dunn’s multiple comparison test following a significant Kruskal Wallis test.****P ≤ 0.0001, ***P ≤ 0.001, **P ≤ 0.01, *P ≤ 0.05, ns = not significant. Scale bar = 1 mm

Fig. 3
figure 3

Inhibition of LRRK2 kinase activity leads to defects in neurons. A qRt-PCR of expression of the wnt4 gene that is necessary for induction of oral identity, levels are significantly decreased in inhibitor treated animals. B Quantification of the neuronal specification gene SoxB, kinase inhibition during embryonic development causes a decrease in SoxB expression levels. ****P ≤ 0.0001, ***P ≤ 0.001, **P ≤ 0.01. C-B Confocal images of Nv-Lwamide-mcherry labelled neurons in developing embryos. By 7 days post fertilization a radial network of axons is visible in control embryos (C, n = 150) while in Lrrk2 kinase treated embryos no network of axons is observed (D, E, n = 160, n = 140). Scale bar = 50 μm. F-H Adult animals exposed to the LRRK2 kinase inhibitor for 7 days exhibit degeneration of the radial axons (G, n = 180) in comparison to control animals (F, n = 150). If the inhibitor treated animals are returned to normal water and re-imaged after 7 days recovery of some of the radial nerve network is observed (H, n = 180)

LRRK kinase activity is necessary for regeneration in Nematostella

Nematostella are well-characterized to be capable of regeneration throughout life [27,28,29, 75,76,77]. The animal can be cut into multiple pieces and each fragment is capable of forming a whole new animal, the only exceptions to this is that a piece of foot alone or amputated tenacles are unable to facilitate full body regeneration. To assess whether LRRK2 has a role in regeneration in N. vectensis, the kinase activity of LRRK2 was inhibited after amputation below the pharynx. In adult animals’ pharynx and tentacle regeneration is complete within 7 days as seen here in control animals (Fig. 4A). However, when exposed to the LRRK2 kinase inhibitor animals appeared to heal the wound but failed to regenerate significant tentacle, often both the pharynx and tentacles were missing. (Fig. 4). Additionally, there was a difference in body length for those treated with the kinase inhibitor, suggesting that like in the homeostatic conditions lack of a functional LRRK2 kinase leads to degeneration of the axons and overall shrinking of the body axis (Fig. 2S). Together, these data suggest that nvLRRK2 plays an important role in the development, homeostasis, and regeneration of neurons.

Fig. 4
figure 4

Pharmacological LRRK2 kinase inhibition impacts regeneration of N. vec tensis tentacles. A Control animals regenerate their oral portions including tentacles withing 12 days in adult animals, in comparison animals exposed to the Lrrk2 kinase inhibitor GNE-0877 fail to regenerate tentacles. B Data shows mean + SD of tentacle length in each of the experimental groups; controls (4dpa n = 44, 8dpa n = 45, 12dpa n = 44), GNE-0877 (4dpa n = 45, 8dpa n = 44, 12dpa n = 44). Non-parametric Dunn’s multiple comparison test following a significant Kruskal Wallis. ns = not significant, ****P ≤ 0.0001, ***P ≤ 0.001, **P ≤ 0.01 *P ≤ 0.05

Discussion

In this study, we have introduced a novel organism for studying LRRK biology, and demonstrated that its LRRK2 orthologue, nvLRRK2, shares key features with hsLRRK2. It has previously been noted that Nematostella are more similar to vertebrates on a genomic scale, than more well-used organisms such as the Drosophila are [70]. This study expands on this and demonstrates its great potential for providing information regarding the function of LRRK2 in a disease context.

Here we have taken advantage of the availability of large number of Nematostella embryos and of commercially available validated specific and potent LRRK2 kinase inhibitors to investigate the function of nvLRRK2. We have uncovered a key role for the kinase activity of this gene in promoting normal embryonic growth, in specification of the oral region of the animal and in neuronal differentiation. Additionally, we have demonstrated that the kinase activity of LRRK2 is necessary to maintain a functional radial network of nerves in the Nematostella and to promote head regeneration where a complex set of nerves must be regenerated. Interestingly we found that in animals that have been exposed to the kinase inhibitor and where the axons are observed to appear to degenerate, these animals can in fact at least partially regenerate the nerve network when the drug is removed. In the future it will be interesting to determine is the new nerves comes from the original cell body or if the actual neuron is being replaced.

This study indicates conservation of a gene which has been primarily studied in vertebrates and illustrates its high conservation in invertebrates. In previous studies, LRRK2 has been identified as making an important contribution to regeneration, it was identified in a large-scale transcriptional profiling approach during planaria regeneration, where the authors suggest it may play a role in activated of the neoblast populations that are essential for planaria regeneration [40]. It will be important to determine in the future the identity of the effector pathway used by nvLRRK2 to regulate regeneration and development, providing mechanistic insight into this process. In mammals it is well-established that a subset of Rab proteins – including Rab10 - are LRRK2 substrates [78], undergoing phosphorylation by LRRK2 [79]. We were able to confirm that Nematostella has a Rab10 ortholog, but in the absence of specific tool antibodies to investigate this were unable to test for direct phosphorylation by nvLRRK2. Although the inhibitors used in this study are specific for LRRK2 (inhibiting the human paralog LRRK1 only at much higher concentrations), it is not possible to exclude inhibition of the other Nematostella LRRK proteins as underlying the phenotypes we observe. Future experiments targeting each of the LRRK genes in Nematostella will test whether the inhibitor treatments are nvLRRK2 specific. Likewise, the conservation of key residues that are mutated in human disease (G2019S and I2020T) with defined biochemical consequences provides an opportunity to assess the implications of gain of LRRK2 kinase function on regeneration and development in the Nematostella system.

Our results open up a new avenue of investigation into the physiological role and importance of LRRK genes and provides a novel model platform to test the consequences of modulating LRRK2 function – with the potential to provide insights relevant to human disease.

Data availability

Nematostella stocks used in these studies are available from the corresponding author. The original data files generated and analyzed for this research will be made available by the authors upon request.

Abbreviations

Dpf:

Days post fertilization

Dpa:

Days post amputation

LRRK:

Leucine rich repeat kinase

NvLRRK2:

Nematostella vectensis leucine rich repeat kinase 2

PD:

Parkinsons disease

PBS:

Phosphate buffered saline

PFA:

Paraformaldehyde

References

  1. Layden MJ, Rentzsch F, Rottinger E. The rise of the starlet sea anemone Nematostella Vectensis as a model system to investigate development and regeneration. Wiley Interdiscip Rev Dev Biol. 2016;5(4):408–28.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Miller DJ, Ball EE, Technau U. Cnidarians and ancestral genetic complexity in the animal kingdom. Trends Genet. 2005;21(10):536–9.

    Article  CAS  PubMed  Google Scholar 

  3. Warner JF et al. NvERTx: a gene expression database to compare embryogenesis and regeneration in the sea anemone Nematostella Vectensis. Development, 2018. 145(10).

  4. Sebe-Pedros A et al. Cnidarian Cell Type Diversity and Regulation Revealed by Whole-Organism Single-Cell RNA-Seq Cell, 2018;173(6):1520–e153420.

  5. Ikmi A, et al. TALEN and CRISPR/Cas9-mediated genome editing in the early-branching metazoan Nematostella vectensis. Nat Commun. 2014;5:5486.

    Article  CAS  PubMed  Google Scholar 

  6. Desvignes T, Pontarotti P, Bobe J. Nme gene family evolutionary history reveals pre-metazoan origins and high conservation between humans and the sea anemone, Nematostella vectensis. PLoS ONE. 2010;5(11):e15506.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Putnam NH, et al. Sea anemone genome reveals ancestral eumetazoan gene repertoire and genomic organization. Science. 2007;317(5834):86–94.

    Article  CAS  PubMed  Google Scholar 

  8. Steger J, et al. Single-cell transcriptomics identifies conserved regulators of neuroglandular lineages. Cell Rep. 2022;40(12):111370.

    Article  CAS  PubMed  Google Scholar 

  9. Zimmermann B, et al. Sea anemone genomes reveal ancestral metazoan chromosomal macrosynteny. bioRxiv; 2020.

  10. Marin I. Ancient origin of the Parkinson disease gene LRRK2. J Mol Evol. 2008;67(1):41–50.

    Article  CAS  PubMed  Google Scholar 

  11. Rentzsch F, et al. FGF signalling controls formation of the apical sensory organ in the cnidarian Nematostella vectensis. Development. 2008;135(10):1761–9.

    Article  CAS  PubMed  Google Scholar 

  12. Amiel AR, et al. A bipolar role of the transcription factor ERG for cnidarian germ layer formation and apical domain patterning. Dev Biol. 2017;430(2):346–61.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Technau U. Gastrulation and germ layer formation in the sea anemone Nematostella vectensis and other cnidarians. Mech Dev. 2020;163:103628.

    Article  CAS  PubMed  Google Scholar 

  14. Salinas-Saavedra M, Rock AQ, Martindale MQ. Germ layer-specific Regul cell Polarity Adhes Gives Insight into Evol Mesoderm Elife, 2018. 7.

  15. Nakanishi N, Martindale MQ. CRISPR knockouts reveal an endogenous role for ancient neuropeptides in regulating developmental timing in a sea anemone. Elife, 2018. 7.

  16. Nakanishi N, et al. Nervous systems of the sea anemone Nematostella vectensis are generated by ectoderm and endoderm and shaped by distinct mechanisms. Development. 2012;139(2):347–57.

    Article  CAS  PubMed  Google Scholar 

  17. Havrilak JA, et al. Characterization of the dynamics and variability of neuronal subtype responses during growth, degrowth, and regeneration of Nematostella vectensis. BMC Biol. 2021;19(1):104.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Havrilak JA, et al. Characterization of NvLWamide-like neurons reveals stereotypy in Nematostella nerve net development. Dev Biol. 2017;431(2):336–46.

    Article  CAS  PubMed  Google Scholar 

  19. Bylund M, et al. Vertebrate neurogenesis is counteracted by Sox1-3 activity. Nat Neurosci. 2003;6(11):1162–8.

    Article  CAS  PubMed  Google Scholar 

  20. Graham V, et al. SOX2 functions to maintain neural progenitor identity. Neuron. 2003;39(5):749–65.

    Article  CAS  PubMed  Google Scholar 

  21. Sandberg M, Kallstrom M, Muhr J. Sox21 promotes the progression of vertebrate neurogenesis. Nat Neurosci. 2005;8(8):995–1001.

    Article  CAS  PubMed  Google Scholar 

  22. Galliot B, et al. Origins of neurogenesis, a cnidarian view. Dev Biol. 2009;332(1):2–24.

    Article  CAS  PubMed  Google Scholar 

  23. Watanabe H, et al. Sequential actions of beta-catenin and bmp pattern the oral nerve net in Nematostella vectensis. Nat Commun. 2014;5:5536.

    Article  CAS  PubMed  Google Scholar 

  24. Jegla T, et al. Expanded functional diversity of shaker K(+) channels in cnidarians is driven by gene expansion. PLoS ONE. 2012;7(12):e51366.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Babonis LS, Martindale MQ. Old cell, new trick? Cnidocytes as a model for the evolution of novelty. Integr Comp Biol. 2014;54(4):714–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. DuBuc TQ, Traylor-Knowles N, Martindale MQ. Initiating a regenerative response; cellular and molecular features of wound healing in the cnidarian Nematostella vectensis. BMC Biol. 2014;12:24.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Amiel AR, et al. Characterization of Morphological and Cellular events underlying oral regeneration in the Sea Anemone, Nematostella vectensis. Int J Mol Sci. 2015;16(12):28449–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Amiel AR, et al. [The sea anemone Nematostella Vectensis, an emerging model for biomedical research: Mechano-sensitivity, extreme regeneration and longevity]. Med Sci (Paris). 2021;37(2):167–77.

    Article  PubMed  Google Scholar 

  29. Amiel AR, Röttinger E. Experimental tools to study regeneration in the Sea Anemone Nematostella vectensis. Methods Mol Biol. 2021;2219:69–80.

    Article  CAS  PubMed  Google Scholar 

  30. Schaffer AA, et al. A transcriptional time-course analysis of oral vs. aboral whole-body regeneration in the Sea anemone Nematostella Vectensis. BMC Genomics. 2016;17(1):718.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Zhu C, Herbst S, Lewis PA. Leucine-rich repeat kinase 2 at a glance. J Cell Sci, 2023. 136(17).

  32. Schlitter AM, et al. The LRRK2 gene in Parkinson’s disease: mutation screening in patients from Germany. J Neurol Neurosurg Psychiatry. 2006;77(7):891–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Berwick DC, Harvey K. LRRK2: an eminence grise of wnt-mediated neurogenesis? Front Cell Neurosci. 2013;7:82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Myasnikov A, et al. Structural analysis of the full-length human LRRK2. Cell. 2021;184(13):3519–e352710.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Alessi DR, Sammler E. LRRK2 kinase in Parkinson’s disease. Science. 2018;360(6384):36–7.

    Article  CAS  PubMed  Google Scholar 

  36. Lewis PA. A step forward for LRRK2 inhibitors in Parkinson’s disease. Sci Transl Med. 2022;14(648):eabq7374.

    Article  CAS  PubMed  Google Scholar 

  37. Hatcher JM, et al. Small-molecule inhibitors of LRRK2. Adv Neurobiol. 2017;14:241–64.

    Article  PubMed  Google Scholar 

  38. Soliman A, Cankara FN, Kortholt A. Allosteric inhibition of LRRK2, where are we now. Biochem Soc Trans. 2020;48(5):2185–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Tasegian A, et al. Impact of type II LRRK2 inhibitors on signaling and mitophagy. Biochem J. 2021;478(19):3555–73.

    Article  CAS  PubMed  Google Scholar 

  40. Benham-Pyle BW, et al. Identification of rare, transient post-mitotic cell states that are induced by injury and required for whole-body regeneration in Schmidtea mediterranea. Nat Cell Biol. 2021;23(9):939–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Klein S et al. Common Environmental pollutants negatively affect development and regeneration in the Sea Anemone Nematostella Vectensis Holobiont. Front Ecol Evol 2021. 9.

  42. Sayers EW, et al. Database resources of the national center for biotechnology information. Nucleic Acids Res. 2022;50(D1):D20–6.

    Article  CAS  PubMed  Google Scholar 

  43. Altschul SF, et al. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.

    Article  CAS  PubMed  Google Scholar 

  44. Altschul SF, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Needleman SB, Wunsch CD. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J Mol Biol. 1970;48(3):443–53.

    Article  CAS  PubMed  Google Scholar 

  46. Cookson MR. The role of leucine-rich repeat kinase 2 (LRRK2) in Parkinson’s disease. Nat Rev Neurosci. 2010;11(12):791–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Mills RD, et al. Analysis of LRRK2 accessory repeat domains: prediction of repeat length, number and sites of Parkinson’s disease mutations. Biochem Soc Trans. 2012;40(5):1086–9.

    Article  CAS  PubMed  Google Scholar 

  48. Rui Q, et al. The role of LRRK2 in neurodegeneration of Parkinson Disease. Curr Neuropharmacol. 2018;16(9):1348–57.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. Wolenski FS, et al. Characterizing the spatiotemporal expression of RNAs and proteins in the starlet sea anemone, Nematostella vectensis. Nat Protoc. 2013;8(5):900–15.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Baptista MAS et al. LRRK2 inhibitors induce reversible changes in nonhuman primate lungs without measurable pulmonary deficits. Sci Transl Med, 2020. 12(540).

  51. Dodson MW, et al. Roles of the Drosophila LRRK2 homolog in Rab7-dependent lysosomal positioning. Hum Mol Genet. 2012;21(6):1350–63.

    Article  CAS  PubMed  Google Scholar 

  52. Eguchi T, et al. LRRK2 and its substrate rab GTPases are sequentially targeted onto stressed lysosomes and maintain their homeostasis. Proc Natl Acad Sci U S A. 2018;115(39):E9115–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Fdez E, et al. Pathogenic LRRK2 regulates centrosome cohesion via Rab10/RILPL1-mediated CDK5RAP2 displacement. iScience. 2022;25(6):104476.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Huber LA, et al. Rab8, a small GTPase involved in vesicular traffic between the TGN and the basolateral plasma membrane. J Cell Biol. 1993;123(1):35–45.

    Article  CAS  PubMed  Google Scholar 

  55. Lara Ordonez AJ, et al. RAB8, RAB10 and RILPL1 contribute to both LRRK2 kinase-mediated centrosomal cohesion and ciliogenesis deficits. Hum Mol Genet. 2019;28(21):3552–68.

    Article  PubMed  PubMed Central  Google Scholar 

  56. Lis P, et al. Development of phospho-specific Rab protein antibodies to monitor in vivo activity of the LRRK2 Parkinson’s disease kinase. Biochem J. 2018;475(1):1–22.

    Article  CAS  PubMed  Google Scholar 

  57. MacLeod DA, et al. RAB7L1 interacts with LRRK2 to modify intraneuronal protein sorting and Parkinson’s disease risk. Neuron. 2013;77(3):425–39.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Malik AU, et al. Deciphering the LRRK code: LRRK1 and LRRK2 phosphorylate distinct Rab proteins and are regulated by diverse mechanisms. Biochem J. 2021;478(3):553–78.

    Article  CAS  PubMed  Google Scholar 

  59. Marin I, van Egmond WN, van Haastert PJ. The Roco protein family: a functional perspective. FASEB J. 2008;22(9):3103–10.

    Article  CAS  PubMed  Google Scholar 

  60. Schmidt SH, et al. Conformation and dynamics of the kinase domain drive subcellular location and activation of LRRK2. Volume 118. Proc Natl Acad Sci U S A; 2021. 23.

  61. Funayama M, et al. A new locus for Parkinson’s disease (PARK8) maps to chromosome 12p11.2-q13.1. Ann Neurol. 2002;51(3):296–301.

    Article  CAS  PubMed  Google Scholar 

  62. Goldwurm S, et al. The G6055A (G2019S) mutation in LRRK2 is frequent in both early and late onset Parkinson’s disease and originates from a common ancestor. J Med Genet. 2005;42(11):e65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Greene ID, et al. Evidence that the LRRK2 ROC domain Parkinson’s disease-associated mutants A1442P and R1441C exhibit increased intracellular degradation. J Neurosci Res. 2014;92(4):506–16.

    Article  CAS  PubMed  Google Scholar 

  64. Kalia LV, et al. Clinical correlations with Lewy body pathology in LRRK2-related Parkinson disease. JAMA Neurol. 2015;72(1):100–5.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Puschmann A, et al. First neuropathological description of a patient with Parkinson’s disease and LRRK2 p.N1437H mutation. Parkinsonism Relat Disord. 2012;18(4):332–8.

    Article  PubMed  Google Scholar 

  66. Ruiz-Martinez J, et al. Penetrance in Parkinson’s disease related to the LRRK2 R1441G mutation in the basque country (Spain). Mov Disord. 2010;25(14):2340–5.

    Article  PubMed  Google Scholar 

  67. Price A, et al. The LRRK2 signalling system. Cell Tissue Res. 2018;373(1):39–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Pasquier C, Robichon A. Evolutionary divergence of phosphorylation to regulate interactive protein networks in Lower and higher species. Int J Mol Sci, 2022. 23(22).

  69. Wijesena N, Simmons DK, Martindale MQ. Antagonistic BMP-cWNT signaling in the cnidarian Nematostella vectensis reveals insight into the evolution of mesoderm. Proc Natl Acad Sci U S A. 2017;114(28):E5608–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  70. Sullivan JC, et al. Conserved and novel wnt clusters in the basal eumetazoan Nematostella vectensis. Dev Genes Evol. 2007;217(3):235–9.

    Article  CAS  PubMed  Google Scholar 

  71. Lee PN, et al. A WNT of things to come: evolution of wnt signaling and polarity in cnidarians. Semin Cell Dev Biol. 2006;17(2):157–67.

    Article  PubMed  Google Scholar 

  72. DuBuc TQ et al. Hox and Wnt pattern the primary body axis of an anthozoan cnidarian before gastrulation Nat Commun, 2018;9(1):2007.

  73. Magie CR, Pang K, Martindale MQ. Genomic inventory and expression of Sox and Fox genes in the cnidarian Nematostella vectensis. Dev Genes Evol. 2005;215(12):618–30.

    Article  CAS  PubMed  Google Scholar 

  74. Richards GS, Rentzsch F. Transgenic analysis of a SoxB gene reveals neural progenitor cells in the cnidarian Nematostella vectensis. Development. 2014;141(24):4681–9.

    Article  CAS  PubMed  Google Scholar 

  75. Bossert P, Thomsen GH. Inducing complete polyp regeneration from the Aboral Physa of the Starlet Sea Anemone Nematostella vectensis. J Vis Exp, 2017(119).

  76. Bossert PE, Dunn MP, Thomsen GH. A staging system for the regeneration of a polyp from the aboral physa of the anthozoan Cnidarian Nematostella vectensis. Dev Dyn. 2013;242(11):1320–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  77. Schaffer AA, et al. A transcriptional time-course analysis of oral vs. aboral whole-body regeneration in the Sea anemone Nematostella Vectensis. BMC Genomics. 2016;17:718.

    Article  PubMed  PubMed Central  Google Scholar 

  78. Steger M et al. Phosphoproteomics reveals that Parkinson’s disease kinase LRRK2 regulates a subset of Rab GTPases. Elife, 2016. 5.

  79. Ordóñez A et al. LRRK2 causes centrosomal deficits via phosphorylated Rab10 and RILPL1 at centriolar subdistal appendages. BioRxiv, 2021.

Download references

Acknowledgements

The authors thank all members of the Lewis and Echeverri Labs for their support during this project. We also thank Michael Layden for sharing the Nv-LWamide-mcherry transgenic line with us.

Funding

G.H. was supported by a travel grant from the Company of Biologists. The collaboration between P.A.L. and K.E. were supported by Royal Society International Exchange Grant IES\R3\203078. KE is supported by a grant from NICHD R01 HD092451, start-up funds from the MBL and funding from the Owens Family Foundation. P.A.L. is a Royal Society Industry Fellow (fellowship IF\R2\222002).

Author information

Authors and Affiliations

Authors

Contributions

GH and SF carried out the experiments. GH, KE and PL prepared the figures. All authors contributed to writing and reviewing the manuscript.

Corresponding authors

Correspondence to Patrick Alfryn Lewis or Karen Echeverri.

Ethics declarations

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

: Figure S1. (A-C) In situ hybridization of NvLrrk2 during Nematostella embryonic development. (D) Quantification real-time PCR analysis of Nematostella NvLrrk2 gene during the first week of development in Nematostella embryos. (E) Quantification real-time PCR analysis of Nematostella wnt1 and col1A genes in control versus inhibitor treated embryos, the expression of these genes is unaffected by the LRRK2 inhibitors.

Supplementary Material 2

: Figure S2. LRRK2 kinase inhibition results in neuronal degeneration in adult Nematostella. NvLWamide-mCherry neurons imaged in control animals, nerve networks are clearly visible (A), after 4 days exposure to kinase inhibitors GNE-0887 (B) or Mli-2 (C) very few axons are observed in adult Nematostella, axons appear to be degenerating. (control N = 60, GNE-0887 n = 52, Mli-2 N = 45). Scale bar = 75µm.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it.The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Holmes, G., Ferguson, S.R., Lewis, P.A. et al. LRRK2 kinase activity is necessary for development and regeneration in Nematostella vectensis. Neural Dev 19, 16 (2024). https://doi.org/10.1186/s13064-024-00193-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13064-024-00193-3

Keywords