Next Article in Journal
Corin Overexpression Reduces Myocardial Infarct Size and Modulates Cardiomyocyte Apoptotic Cell Death
Next Article in Special Issue
CRISPR/Cas9-Mediated Gene Correction to Understand ALS
Previous Article in Journal
Rubrofusarin Attenuates Chronic Restraint Stress-Induced Depressive Symptoms
Previous Article in Special Issue
Modeling Cardiovascular Diseases with hiPSC-Derived Cardiomyocytes in 2D and 3D Cultures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Iroquois Homeobox Protein 2 Identified as a Potential Biomarker for Parkinson’s Disease

1
Stem Cell Convergence Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon 34141, Korea
2
Department of Functional Genomics, KRIBB School of Bioscience, University of Science and Technology, Daejeon 34113, Korea
3
Group for Biometrology, Korea Research Institute of Standards and Science (KRISS), Daejeon 34113, Korea
4
Department of Molecular & Life Sciences, College of Science & Technology, Hanyang University, Ansan 15588, Korea
5
Mucosal Immunology Laboratory, Department of Convergence Medicine, University of Ulsan College of Medicine/Asan Medical Center, Seoul 05505, Korea
6
Department of Pharmacology, College of Medicine, Dankook University, Cheonan 31116, Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2020, 21(10), 3455; https://doi.org/10.3390/ijms21103455
Submission received: 27 March 2020 / Revised: 5 May 2020 / Accepted: 12 May 2020 / Published: 14 May 2020
(This article belongs to the Special Issue Disease Modeling Using Human Induced Pluripotent Stem Cells 2.0)

Abstract

:
The diagnosis of Parkinson’s disease (PD) is initiated after the occurrence of motor symptoms, such as resting tremors, rigidity, and bradykinesia. According to previous reports, non-motor symptoms, notably gastrointestinal dysfunction, could potentially be early biomarkers in PD patients as such symptoms occur earlier than motor symptoms. However, connecting PD to the intestine is methodologically challenging. Thus, we generated in vitro human intestinal organoids from PD patients and ex vivo mouse small intestinal organoids from aged transgenic mice. Both intestinal organoids (IOs) contained the human LRRK2 G2019S mutation, which is the most frequent genetic cause of familial and sporadic PD. By conducting comprehensive genomic comparisons with these two types of IOs, we determined that a particular gene, namely, Iroquois homeobox protein 2 (IRX2), showed PD-related expression patterns not only in human pluripotent stem cell (PSC)-derived neuroectodermal spheres but also in human PSC-derived neuronal cells containing dopaminergic neurons. We expected that our approach of using various cell types presented a novel technical method for studying the effects of multi-organs in PD pathophysiology as well as for the development of diagnostic markers for PD.

1. Introduction

Parkinson’s disease (PD) is a common and complex neurodegenerative disorder that results from the progressive loss of midbrain dopaminergic neurons in the substantia nigra (SN). The main cause of PD may be a complicated interplay of genetic and environmental factors [1]. However, it is now recognized that the disease is a multisystem disorder involving both motor and non-motor symptoms associated with the central and peripheral nervous systems [2]. Throughout the past few decades, several genes associated with PD have been discovered through genome-wide association studies (GWAS), such as the alpha-synuclein (SNCA), glucocerebrosidase (GBA), and leucine-rich repeat kinase 2 (LRRK2) genes [3,4]. For the diagnosis of PD, clinical criteria based on the cardinal signs of PD and brain imaging of structural, functional, and molecular changes have been used. However, typical motor symptoms occur when up to 60% of the dopaminergic neurons of the SN have already degenerated [5]. It is, therefore, essential to develop an early diagnosis marker for the successful treatment for PD. According to previous studies, the onset of motor symptoms can be preceded by a premotor or prodromal phase characterized by specific non-motor symptoms [6]. These include emotional problems, cognitive dysfunction, sleep disturbances, sensory manipulation, as well as autonomic dysregulations [7]. Gastrointestinal (GI) dysfunction is a common non-motor symptom of PD and is observed at all stages through GI symptoms, including dysphagia, gastroparesis, and constipation [8]. The results of several studies suggest that infrequent bowel movements or constipation severity are associated with an elevated risk of PD in the future [9,10]. Moreover, the pathogenic accumulation of phosphorylated alpha-synuclein protein, the neuropathological hallmark of PD, in the bowel can occur during the early stages of the disease, prior to the onset of motor symptoms [11]. The risk of PD was lower in patients who underwent truncal vagotomy, which indicates the linkage between the brain and the GI tract is influential in the development of PD [12,13].
The pathological mutations of LRRK2 affect protein functions and molecular pathways involved in PD [14]. In particular, Gly2019Ser substitution (G2019S) in LRRK2, a common genetic variation of familial and sporadic PD, increases kinase activity and induces toxicity [15]. A recent study has identified the mutated LRRK2 gene through exome sequencing of patients with Crohn’s disease (CD) and ulcerative colitis (UC), which indicates that LRRK2 mutations are associated with intestinal inflammatory disorders [16]. Nevertheless, the molecular association between PD-related mutations and intestinal dysfunction remains unclear. Thus, the study of gene expression profiles in PD patient-derived intestinal samples could provide clues to the pathogenesis of PD. In our previous study, in vitro intestinal organoids (IOs) and neuroectodermal spheres (NESs) derived from PD patient’s pluripotent stem cells (PSCs) showed a highly distinct gene expression patterns [17].
Here, we generated PD-relevant three-dimensional (3D) human IOs from PD patient-specific PSCs and mouse IOs from aged mouse intestines harboring the LRRK2 G2019S mutation. Furthermore, we performed transcriptome profiles of human and mouse IOs to determine the molecular association and found a candidate gene: IRX2. This newly identified gene exhibited PD-related expression patterns in PD patient-derived neural cells, including NESs, and neuronal cells containing dopaminergic neurons.

2. Results

2.1. Generation of Mouse Intestinal Organoids from Human LRRK2 G2019S Transgenic Mice and Littermate Controls

In our previous study, we observed the comprehensive molecular changes caused by the LRRK2 G2019S mutation in neural and intestinal 3D culture model systems based on the same PD-specific PSCs [17]. One of the most interesting aspects of the results was that the gene expression difference was more distinct in human intestinal organoids (hIOs) than in human neuroectodermal spheres (hNESs), even though PD is a well-known neurodegenerative disorder. To explore the potential of PD patient-derived hIOs and investigate molecular targets from LRRK2 mutations in PD, we compared differentiated PD patient-derived IOs with ex vivo IOs from PD model mice through microarray analysis. Differentially expressed genes (DEGs) were collected to identify common factors to study the signature of PD concomitants with the intestine.
To reflect age-dependent accumulative PD characteristics, we generated mouse small intestinal organoids (mIOs) from human LRRK2 G2019S transgenic (TG) mice and normal control littermates at 15–18 months of age. Ex vivo cultured mouse’s small intestinal epithelial cells exhibited a bud-like structure after 4 to 7 days of culture on Matrigel (Figure 1A). There were no differences in the forming efficiency, size, and budding structures between the human LRRK2 G2019S mutant mouse small intestinal organoids (GS mIOs) and non-transgenic littermate mouse small intestinal organoids (WT mIOs). GS mIOs were morphologically indistinguishable from WT mIOs after several passages. Moreover, the qRT-PCR analysis showed that GS mIOs and WT mIOs contained similar expression levels of intestinal cell markers, including Lgr5 (an intestinal stem cell marker), Vil1 (villin 1 for enterocytes), Lyz1 (lysozyme for Paneth cells), ChgA (chromogranin A for enteroendocrine cells), and Muc2 (mucin 2 for goblet cells) (Figure 1B). We assessed the expression and localization of the intestinal cell markers via immunostaining (Figure 1C). Despite the fact that the behavioral symptoms of PD induced by the expressions of mutant LRRK2 are found in mice of 12 months and older [18], ex vivo cultured GS mIOs showed similar characteristics to WT mIOs.

2.2. Comparative Transcriptome Analysis of Ex Vivo Cultured mIOs and In Vitro Differentiated hIOs

To understand the alterations in PD using two types of PD-intestine models, we decided to perform microarray analysis to determine the genes with expression levels that were affected by the LRRK2 G2019S mutation. First of all, we selected DEGs that were significantly altered by the LRRK2 G2019S mutation and analyzed transcriptome alterations (Figure 2). Compared to WT mIOs, 1225 genes were significantly up- or down-regulated in GS mIOs. According to a fold-change threshold of 2.0 between WT mIOs and GS mIOs, 148 genes were up-regulated, and 127 genes were down-regulated (Figure 2A, left). In addition, we re-analyzed our previous microarray data associated with pluripotent stem cell (PSC)-derived differentiated hIOs of LRRK2 G2019S PD patients [17] to select the common factors between the two types of models. There was a substantially greater number of significantly altered genes (16,067 genes) in hIOs with the LRRK2 G2019S mutation compared to in healthy control PSC-derived hIOs. Accordingly, greater numbers of up-regulated (3464) and down-regulated (4418) genes were identified in hIOs with the LRRK2 G2019S mutation when we applied a fold-change threshold of 2.0 (Figure 2A, right). The principal component analysis showed that the LRRK2 G2019S mutant and the control were separated in both mIOs and hIOs (Figure 2B). Heatmap analysis also showed partial segregation of the expression patterns (Figure 2C). These results suggested that the gene expression profiles were altered by the LRRK2 G2019S mutation in ex vivo cultured and PSC-derived differentiated IOs.

2.3. Identification and Validation of Common Factors between Two Types of Intestinal Organoids of PD

Through transcriptome analysis on the two types of PD-intestine models, we attempted to select commonly detectable up-regulated or down-regulated genes with a fold-change threshold of 5.0 between the LRRK2 G2019S mutant and the control (Figure S1). We identified four common factors, including IRX2 (Iroquois homeobox protein 2), LMO3 (LIM domain only protein 3), ANXA10 (AnnexinA10), and TFF2 (Trefoil factor 2) (Figure 3A). Microarray data indicated that IRX2 and LMO3 were up-regulated in the LRRK2 G2019S mutant, whereas ANXA10 and TFF2 were down-regulated (Figure 3B). The qRT-PCR analysis confirmed that the expression levels of selected genes were altered in mIOs and hIOs with similar patterns that are represented by the microarray analysis (Figure 3C).

2.4. Confirmation of Commonly Regulated Factors in NESs and Neuronal Cells Containing Dopaminergic Neurons

Previous studies suggest that the gut-brain connection is associated with the initiation and progression of PD [19]. Therefore, we investigated the common factors of intestinal organoids on neural cells. First, we checked the expression patterns of IRX2, LMO3, ANXA10, and TFF2 in human PSC-derived NESs (Figure 4A) [17]. Only IRX2 showed a similar pattern, which indicated that the IRX2 expression was significantly up-regulated in human LRRK2 G2019S PD patient-derived NESs compared to in wild-type control NESs. We also validated the expression levels through qRT-PCR analysis (Figure 4B).
According to previous reports, familial PD patients harboring the G2019S mutation in the LRRK2 gene were clinically indistinguishable from idiopathic PD, showing similarities in terms of motor deficits [20], clinical features [21,22], and neurochemical phenotypes [23]. Thus, we differentiated human PSCs from familial PD patients harboring the LRRK2 G2019S mutation and from sporadic PD patients into neuronal cells containing dopaminergic neurons, which are known to be affected by PD progression. Thirty-five days differentiated neuronal cells exhibited the expression of the dopaminergic neuron-specific marker tyrosine hydroxylase (TH) and the mature neuron marker microtubule-associated protein 2 (MAP2) (Figure 5A). We confirmed cell survival and analyzed TH and MAP2 double-positive neuronal cells from familial/sporadic PD and healthy control (WT)-derived induced pluripotent stem cells (Figure S2). We found that there were no differences in neuronal cell death between WT and PD through the expression levels of cleaved caspase-3 (Figure S2A,B). In addition, there were no significant differences between healthy control (WT) and the PD-group in terms of the expressions of TH and MAP2 double-positive cells (Figure S2C). We detected elevated expressions of IRX2 in familial LRRK2 G2019S and sporadic PD patient-derived neuronal cells in comparison with wild-type control cells (Figure 5B), as shown in human LRRK2 G2019S PD patient-derived IOs and NESs. Finally, we also checked the expression levels of IRX2 in familial LRRK2 G2019S and sporadic PD patient-derived fibroblasts to determine the feasibility of IRX2 as a potential biomarker (Figure 5C). Interestingly, IRX2 showed similar expression patterns in both familial and sporadic PD patient-derived fibroblasts, which are known as appropriate samples for diagnostic purposes. In summary, these data suggested that significantly increased expression levels of IRX2 could be used as a PD-specific biomarker.

3. Discussion

Human development and disease are challenging topics of study due to the low experimental availability of in vivo systems and the complex biological networks that are involved [24]. Although experimental animal models have been used to study the pathogenesis of the disease, the use of model organisms often does not translate into the replication of human biology as such models possess different structures and physiology [25,26]. For these reasons, 3D organoid culture systems are more preferred in recent years and serve as an excellent platform for modeling human tissue or organs.
Organoids have been widely adopted in disease pathology, regenerative medicine, and drug development, including in pharmacokinetics analysis, efficacy evaluation, and toxicity testing [27]. In particular, 3D organotypic cultures from disease-specific PSCs represent a greater number of relevant characteristics and symptoms of complex human diseases [28,29]. Accordingly, patient-specific 3D organoid models represent a powerful tool for investigating disease pathology.
In our previous study, we observed the comprehensive molecular changes caused by the LRRK2 G2019S mutation in 3D culture model systems based on PD PSCs. The expression profiles of PD patient-derived cells were distinct from those of wild-type controls in hIOs as well as hNESs containing neural progenitor cells. This patient-derived 3D culture model system resulted in extensive alterations in global gene expression, such as synaptic transmission, specifically synaptic vesicle trafficking, and defects of which are known to be related to PD [17]. As a result, IOs harboring the LRRK2 G2019S mutation could serve as a versatile model to explain the association between PD and intestinal alterations.
A number of studies have shown that LRRK2 affects pathophysiological modifications in the GI pathogenesis of PD [6,7,10,17]. The mechanism of LRRK2-associated GI pathogenesis has been examined in LRRK2 TG mice and knockout mice [30,31,32]. Thus, in this study, we developed mIOs as an ex vivo system and compared them with differentiated hIOs through transcriptome analysis to identify molecular targets that are regulated by LRRK2 mutation. The results showed that the two types of IOs exhibited distinctive gene expressions induced by the LRRK2 G2019S mutation. We subsequently selected the genes that were significantly up-regulated or down-regulated in the two types of intestinal organoids harboring LRRK2 G2019S to identify the common factors. As a result, four genes were selected: IRX2, LMO3, ANXA10, and TFF2. Two genes (IRX2, LMO3) were up-regulated, and the remaining two (ANXA10, TFF2) were down-regulated by the expression of mutant LRRK2. Taken together, our results suggested that these genes might be potential factors that could possibly emerge in PD pathology.
According to a recent report, the gut-to-brain propagation of pathogenic proteins induces PD through dopaminergic neurodegeneration [33]. Thus, it is worthwhile to verify the expression levels of the selected genes from IOs models of PD in neural cells. Notably, IRX2, one of the selected genes from the LRRK2 G2019S IOs models, exhibited a similar expression pattern in hNESs from LRRK2 G2019S PD patient-derived PSCs. IRX2 is a member of the Iroquois homeobox genes that are expressed in the cerebellum and applies multiple functions during animal development [34,35]. Vertebrate Irx genes participate in the regionalization of the brain and have been reported to exhibit functions in neural development [35]. On the other hand, Irx2-deficient mice have shown normal developmental phenotypes through the dynamic expression of Irx2 in the developing heart, nervous system, and other organs [36]. Most vertebrates include six genes and two clusters of three genes, with the IrxA cluster comprised of Irx1, Irx2, and Irx4, and the IrxB cluster comprised of Irx3, Irx5, and Irx6 [37,38]. Recently, GWAS of PD has discovered the potential functionality of risk correlated with single nucleotide polymorphisms (SNPs) based on non-coding DNA, including IRX2, 3, and 6 [39]. An additional transcriptome mapper (TRAM) profiling has been conducted in post–mortem whole SN tissue and single dopamine neurons from familial and sporadic PD patients to identify chromosomal segments and gene loci diverse expression. TRAM has determined 10 significant genes, including IRX2, that have been up-regulated or down-regulated in SN and midbrain dopaminergic neurons of PD patients [40]. In addition, mass spectrometry-based phosphoproteomics of genetically engineered LRRK2 mice has led to the discovery of several key substrates, including the Irx1 gene [41]. Accordingly, we were able to confirm the elevated expressions of IRX2, which was also discovered in LRRK2 G2019S PD patient PSC-derived neuronal cells containing dopaminergic neurons. Even though previous studies have shown conflicting results [34,35,36], the evolutionarily conserved IRX2 gene, which showed elevated expression patterns not only in neural cells but also in IOs, would conceivably be associated with PD progression. In addition, we also investigated the expression patterns of the Iroquois family of homeobox genes in hNESs and IOs (Figure S3). In hIOs, six genes were up-regulated in the LRRK2 G2019S PD patient-derived cells when compared to the wild-type control cells. However, the IrxA cluster (IRX1, IRX2, and IRX4) showed a significantly greater increase in the LRRK2 G2019S PD patient-derived cells over the wild-type control cells compared to the IrxB cluster (IRX3, IRX5, and IRX6) in hNESs. We further examined the expression of IRX2 in fibroblasts of familial and sporadic PD patients. IRX2 showed an expression pattern similar to those shown in the intestinal and neural models of PD. Through these results, we suggested that increased levels of IRX2 could be used to detect PD in various cell types, including fibroblasts as well as intestinal and neuronal cells containing dopaminergic neurons. Although we were not able to determine the role of the increase in IRX2 in neural and intestinal cells during the progression of PD and underline the mechanisms in this study, the newly discovered IRX2 is expected to be an interesting factor for upcoming research. It may be worthwhile to check the expression levels of IRX2 with a larger group of PD patients.
All in all, IRX2 could be a promising factor in studying molecular alterations in PD and could be used as a diagnostic marker for PD. Furthermore, intestinal organoids could be appropriate PD models even though they are composed of the only epithelium. Our multi-organ culture system could lead to a better understanding of PD.

4. Materials and Methods

4.1. Animals

LRRK2 transgenic mice expressing the G2019S mutation [FVB/N-Tg (LRRK2*G2019S) 1Cjli/J] were purchased from Jackson Laboratory (Bar Harbor, ME, USA). The mice were maintained under a 12 h light/dark cycle and given free access to food and water. All experimental procedures followed guidelines approved by the Institutional Animal Care and Use Committee of Hanyang University (HY-IACUC-12-018).

4.2. Generation of Mouse Small Intestinal Organoids

Mouse small intestinal organoids (mIOs) were isolated and established according to the manufacturer’s instructions using the IntestiCultTM Organoid Growth Medium (STEMCELL Technologies, Vancouver, BC, Canada). Mouse small intestines were harvested about 10–20 cm in length. Fine forceps were used to remove the membrane, blood vessels, and fat from the exterior of the tissue. The intestinal segments were subsequently flushed with cold Dulbecco’s Phosphate-Buffered Saline (DPBS, Welgene, Gyeongsan-si, Korea) and were opened using small scissors. The intestinal sheets were thoroughly rinsed with clean DPBS, and the washed intestine sheets were chopped into 2 mm pieces. Tissue pieces were pipetted up and down gently in DPBS using a pre-wetted 10 mL serological pipette. For the dissociation of tissue pieces and isolation of crypts, the tissue pieces were treated with Gentle Cell Dissociation Reagent (STEMCELL Technologies, Vancouver, BC, Canada) for 15–20 min at room temperature. The pellets were resuspended in cold PBS containing 0.1% BSA (Sigma-Aldrich, St. Louis, MO, USA) and were subsequently filtered through a 70 μm strainer (Corning Inc., Corning, NY, USA) to obtain crypts of appropriate sizes. The selected intestinal crypts were mixed with cold DMEM/F-12 (Thermo Fisher Scientific, Waltham, MA, USA). Through counting, 3000 crypts were collected. To produce the solid dome, the crypts were resuspended in the IntestiCultTM Organoid Growth Medium and mixed with Matrigel Matrix (BD Biosciences, San Jose, CA, USA) with a 1:1 ratio. The 50 µL of the crypts and Matrigel matrix mixture was incubated at 37 °C for 10 min. IntestiCultTM Organoid Growth Medium was added to the sidewall of the well. The culture medium was half-exchanged every three days. Organoids were passaged using the Gentle Cell Dissociation Reagent after 10 to 14 days of culture.

4.3. Quantitative Real-Time PCR (qRT-PCR)

Total RNA was extracted from the organoids using the RNeasy Plus Mini Kit (QIAGEN, Hilden, Germany) and reverse-transcribed using the iScriptTM cDNA Synthesis Kit (Bio-Rad Laboratories, Inc., Hercules, CA, USA). The qRT-PCR was performed using the 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The primers used in this study are listed in Table S1. All experiments were performed in triplicate.

4.4. Western Blot Analysis

Cells were treated with cold sample lysis buffer (1% Triton X-100, Xpert Protease Inhibitor Cocktail Solution (GenDEPOT, Katy, TX, USA), 5 mM Ethylenediaminetetraacetic acid (EDTA, Thermo Fisher Scientific, Waltham, MA, USA), and 1 mM Phenylmethanesulfonyl fluoride (PMSF, Thermo Fisher Scientific, Waltham, MA, USA)). Protein Assay Dye Reagent Concentrate (Bio-Rad Laboratories Inc., Hercules, CA, USA) was used for determining the protein concentration. An equal amount of protein was loaded for SDS-PAGE and transferred onto PVDF membrane using Wet/Tank Blotting Systems (Bio-Rad Laboratories Inc., Hercules, CA, USA). The membranes were incubated with 5% DifcoTM Skim Milk (BD Biosciences, San Jose, CA, USA) for 30 min at room temperature. The membranes were incubated with cleaved caspase-3 (Asp175) antibody (Cat. 9661S, Cell Signaling Technology, Danvers, MA, USA) and GAPDH antibody (Cat. sc-47724, Santacruz Biotechnology Inc., Dallas, TX, USA) overnight at 4 °C, followed by incubation with horseradish peroxidase (HRP)-conjugated secondary antibodies (Cell Signaling Technology, Danvers, MA, USA). We added ECLTM Select Western Blotting Detection Reagent (GE Healthcare, Chicago, IL, USA) on the membrane for the detection of the signals from HRP. The images of protein bands were obtained by a LAS-3000 imaging system (Fujifilm, Minato, Tokyo, Japan).

4.5. Immunofluorescence of Mouse Small Intestinal Organoids and Human Neuronal Cells Containing Dopaminergic Neurons

Whole mIOs were fixed overnight in fresh 4% paraformaldehyde. Fixed organoids were unmasked in PBS through shaking at room temperature and were blocked and permeabilized with 5% BSA and 0.1% Triton X-100 in PBS. The primary antibodies that were used were Villin (Cat. sc-58897, Santa Cruz Biotechnology Inc., Dallas, TX, USA), Mucin 2 (Cat. sc-15334, Santa Cruz Biotechnology Inc., Dallas, TX, USA), Chromogranin A (Cat. ab15160, Abcam, Cambridge, UK), and Lysozyme (Cat. ab108508, Abcam, Cambridge, UK). Human neuronal cells containing dopaminergic neurons were fixed with 80% Methanol for 10 min. After blocking with 5% BSA, the primary antibodies that were used were MAP2 (Cat. ab5392, Abcam, Cambridge, UK) and TH (Cat. Ab152, Millipore Sigma, Burlington, MA, USA). The primary antibodies were bound overnight at 4 °C. The secondary antibodies that were used were donkey anti-mouse and anti-rabbit Alexa Fluor 488 and 594 (Invitrogen, Carlsbad, CA, USA). The secondary antibodies were bound for 1 h at room temperature. DNA was stained with Hoechst 33342 (Cat. H3570, Thermo Fisher Scientific, Waltham, MA, USA). The images of the crypt organoids and human neuronal cells containing dopaminergic neurons were acquired using an LSM 800 confocal microscope with Zeiss software.

4.6. Microarray and Transcriptome Analysis

Microarray analysis was performed using the Low Input Quick Amp Labeling Kit (Agilent Technologies, Santa Clara, CA, USA), RNase Mini-Column (QIAGEN, Hilden, Germany), and Mouse V2 Gene expression 4 × 44K (Agilent Technologies, Santa Clara, CA, USA), according to the manufacturers’ protocols. Data quantification was achieved by using the Feature Extraction software 10.7, and global normalization was conducted using GeneSpringGX 7.3.1 (Agilent Technologies, Santa Clara, CA, USA). Heatmap clustering was visualized using the MeV 4.9.0 software. Principal component analysis (PCA) of the differentially expressed genes was reflected using scikit-learn, Python package. Functional annotation clustering and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were analyzed using David Bioinformatics Resources 6.7 (https://david-d.ncifcrf.gov). We used relative values to compare the expression levels of the selected genes, meaning that the global normalized value of WT1 was used as the standard point of the mouse results, whereas WT#1 was used as the standard point of the human results.

4.7. Differentiation of Human-Induced Pluripotent Cells into Human Intestinal Organoids, Human Neuroectodermal Sphere, and Culture

Human wild-type (WT) and Parkinson’s disease patient (PD) somatic cells (AG02261, GM01680, GM01706, GM02036, GM02623, GM23967, ND14317, ND29492, ND33879, ND38262, AG20439, AG20442, AG20443, AG20445, and AG20446) were purchased from the Coriell Institute (Camden, NJ, USA). The human embryonic stem cell line (H9) was purchased from the WiCell Research Institute (Madison, WI, USA). Healthy fibroblasts (IMR90, MRC5, and CRL-2097) were purchased from the American Type Culture Collection (Manassas, VA, USA). Based on these cell lines, human intestinal organoids (hIOs) and human neuroectodermal spheres (hNESs) were generated, as previously described [17]. To differentiate hindgut spheroids, human pluripotent stem cells (PSCs) were treated with 100 ng/mL Activin A (Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), 500 ng/mL Fibroblast growth factor 4 (FGF4, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), and 500 ng/mL Wnt Family Member 3A (WNT3A, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA). The spheroids were cultured in an hIO media containing 1× B27 (Invitrogen, Carlsbad, CA, USA), 100 ng/mL Epidermal Growth Factor (EGF, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), 500 ng/mL R-Spondin 1 (Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), and 100 ng/mL Noggin (Research and Diagnostic Systems, Inc., Minneapolis, MN, USA). To generate hNESs, Embryonic bodies (EBs) were cultured in a neural progenitor culture medium that includes DMEM/F12, 1X N2/B27 (Invitrogen, Carlsbad, CA, USA), 20 ng/mL basic Fibroblast Growth Factor (bFGF, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), 20 ng/mL EGF (PeproTech, Inc., Rocky Hill, NJ, USA), and 10 ng/mL human Leukemia inhibitory factor (LIF, PeproTech, Inc., Rocky Hill, NJ, USA). Our reference samples were purchased from the Coriell Institute (Camden, NJ, USA), WiCell Research Institute (Madison, WI, USA), and American Type Culture Collection (Manassas, VA, USA) and have been determined to be exempt from IRB review.

4.8. Generation of Human Neuronal Cells Containing Dopaminergic Neurons

Human neuronal cells containing dopaminergic neurons were differentiated from human PSCs, as previously described [42,43], with some modifications. At first, human PSCs were differentiated into human primitive neural stem cells (hpNSCs) with 10 ng/mL human LIF (PeproTech, Inc., Rocky Hill, NJ, USA), 4 μM CHIR99021 (Tocris Bioscience, Bristol, UK), 3 μM SB431542 (Tocris Bioscience, Bristol, UK), and 0.1 μM Compound E (Millipore Sigma, Burlington, MA, USA) for 7 days. Treatment with 2 μM Dorsomorphin (Sigma-Aldrich, St. Louis, MO, USA) was performed during the first 2 days of differentiation. The cells were subsequently passaged and maintained in 10 ng/mL human LIF, 3 μM CHIR99021, and 2 μM SB431542. The hpNSCs were treated with 100 ng/mL Sonic hedgehog (SHH, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), 100 ng/mL Fibroblast Growth Factor 8 (FGF8, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), and 2 μM Purmorphamine (Tocris Bioscience, Bristol, UK) for 7 days. Treatment with 20 ng/mL Brain-derived Neurotrophic Factor (BDNF, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), Glial cell-derived Neurotrophic Factor (GDNF, Research and Diagnostic Systems, Inc., Minneapolis, MN, USA), and 1 ng/mL Transforming Growth Factor beta-3 (TGF-b3, PeproTech, Inc., Rocky Hill, NJ, USA) was performed for 20 days to differentiate into neuronal cells containing dopaminergic neurons.

4.9. Statistical Analysis

Statistical analysis was performed on GraphPad Prism V5 (San Diego, CA, USA). All data are presented as means ± standard error of means (SEM). The statistical results were analyzed through the unpaired, two-tailed Student’s t-test. p-values less than 0.05 were considered as statistically significant. All assays were performed with three-five technical replicates. The number of independent experiments and standard error for each experiment are represented in the relevant manuscript and figure legends.

5. Conclusions

We proposed several genes as potential therapeutic targets from 3D organoid culture systems based on PD-specific mIOs and hIOs containing the LRRK2 G2019S mutation. We observed the molecular changes caused by the LRRK2 G2019S mutation and identified the IRX2 gene. IRX2 showed similar expression patterns in PSC-derived hNES and PSC-derived human neuronal cells containing dopaminergic neurons. Moreover, we could also detect the elevated expression of IRX2 in sporadic PD patient-derived neural cells. Based on our results, we expected that the newly identified genes from our comparative analysis could be used in studying PD pathology and developing diagnostic markers for familial and sporadic PD.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/10/3455/s1, Figure S1. Hierarchically clustered heatmap of the differentially expressed genes between the LRRK2 G2019S mutant and the control with a fold-change threshold of 5.0 in (A) mIOs and (B) hIOs. Figure S2. Analysis of WT and PD-derived neuronal cells containing dopaminergic neurons. Figure S3. Expression patterns of the Iroquois family of homeobox genes in hIOs and hNESs. Table S1. Primer sequences used in this study.

Author Contributions

H.S. (Hyuna Sim), J.-E.L., M.-Y.S., and J.K. (Jisun Kim) conceived and designed the experiments, analyzed the data, and wrote the manuscript. H.M.Y., J.K. (Jisun Kim), S.C., H.L., A.B., H.S. (Hyemyung Seo), M.-N.K., H.L., Y.-J.J., and H.G.K. analyzed the data. M.-Y.S. and J.K. (Janghwan Kim) revised the manuscript and gave the final approval for submission. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the grant from the National Research Foundation of Korea (2015M3A9C7030128), the National Research Foundation of Korea (NRF) funded by the Korean government (MSIT) (NRF-2018M3A9H3023077), and the KRIBB research initiative program funded by the Ministry of Science and ICT.

Acknowledgments

Our reference samples were purchased from the Coriell Institute and ATCC and have been determined to be exempt from IRB review.

Conflicts of Interest

The authors declare that they have no competing interests.

Abbreviations

3Dthree-dimensional
ANXA10AnnexinA10
CDCrohn’s disease
DEGsDifferentially expressed genes
G2019SGly2019Ser substitution
GBAGlucocerebrosidase
GIGastrointestinal
GS mIOsHuman LRRK2 G2019S mutant mouse small intestinal organoids
GWASGenome-wide association studies
hIOsHuman intestinal organoids
hNESsHuman neuroectodermal spheres
IOsIntestinal organoids
iPSCsInduced pluripotent stem cells
IRX2Iroquois homeobox protein 2
LMO3LIM domain only protein 3
LRRK2Leucine-rich repeat kinase 2
MAP2Microtubule-associated protein 2
mIOsMouse small intestinal organoids
NESsNeuroectodermal spheres
PDParkinson’s disease
PSCPluripotent stem cells
SNSubstantia nigra
SNCAAlpha-synuclein
SNPsSingle nucleotide polymorphisms
TFF2Trefoil factor 2
TGTransgenic
THTyrosine hydroxylase
TRAMTranscriptome Mapper
UCUlcerative colitis
WT mIOsNon-transgenic littermate mouse small intestinal organoids

References

  1. Kalia, L.V.; Lang, A.E. Parkinson’s disease. Lancet 2015, 386, 896–912. [Google Scholar] [CrossRef]
  2. Poewe, W.; Seppi, K.; Tanner, C.M.; Halliday, G.M.; Brundin, P.; Volkmann, J.; Schrag, A.E.; Lang, A.E. Parkinson disease. Nat. Rev. Dis. Prim. 2017, 3, 1–21. [Google Scholar] [CrossRef] [PubMed]
  3. Nalls, M.A.; Pankratz, N.; Lill, C.M.; Do, C.B.; Hernandez, D.G.; Saad, M.; Destefano, A.L.; Kara, E.; Bras, J.; Sharma, M.; et al. Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat. Genet. 2014, 46, 989–993. [Google Scholar] [CrossRef]
  4. Chang, D.; Nalls, M.A.; Hallgrímsdóttir, I.B.; Hunkapiller, J.; van der Brug, M.; Cai, F.; Kerchner, G.A.; Ayalon, G.; Bingol, B.; Sheng, M.; et al. A meta-analysis of genome-wide association studies identifies 17 new Parkinson’s disease risk loci. Nat. Genet. 2017, 49, 1511–1516. [Google Scholar] [CrossRef] [PubMed]
  5. Berg, D. In vivo detection of iron and neuromelanin by transcranial sonography - A new approach for early detection of substantia nigra damage. J. Neural Transm. 2006, 113, 775–780. [Google Scholar] [CrossRef] [PubMed]
  6. Goldman, J.G.; Postuma, R. Premotor and nonmotor features of Parkinson’s disease. Curr. Opin. Neurol. 2014, 27, 434–441. [Google Scholar] [CrossRef] [Green Version]
  7. Poewe, W. Non-motor symptoms in Parkinson’s disease. Eur. J. Neurol. 2008, 15, 14–20. [Google Scholar] [CrossRef]
  8. Chen, H.; Zhao, E.J.; Zhang, W.; Lu, Y.; Liu, R.; Huang, X.; Ciesielski-Jones, A.J.; Justice, M.A.; Cousins, D.S.; Peddada, S. Meta-analyses on prevalence of selected Parkinson’s nonmotor symptoms before and after diagnosis. Transl. Neurodegener. 2015, 4, 1. [Google Scholar] [CrossRef] [Green Version]
  9. Lin, C.H.; Lin, J.W.; Liu, Y.C.; Chang, C.H.; Wu, R.M. Risk of Parkinson’s disease following severe constipation: A nationwide population-based cohort study. Park. Relat. Disord. 2014, 20, 1371–1375. [Google Scholar] [CrossRef]
  10. Abbott, R.D. Frequency of bowel movements and future risk of Parkinson’s disease. Neurology 2001, 57, 456–462. [Google Scholar] [CrossRef]
  11. Hilton, D.; Stephens, M.; Kirk, L.; Edwards, P.; Potter, R.; Zajicek, J.; Broughton, E.; Hagan, H.; Carroll, C. Accumulation of α-synuclein in the bowel of patients in the pre-clinical phase of Parkinson’s disease. Acta Neuropathol. 2014, 127, 235–241. [Google Scholar] [CrossRef] [PubMed]
  12. Svensson, E.; Horváth-Puhó, E.; Thomsen, R.W.; Djurhuus, J.C.; Pedersen, L.; Borghammer, P.; Sørensen, H.T. Vagotomy and subsequent risk of Parkinson’s disease. Ann. Neurol. 2015, 78, 522–529. [Google Scholar] [CrossRef] [PubMed]
  13. Liu, B.; Fang, F.; Pedersen, N.L.; Tillander, A.; Ludvigsson, J.F.; Ekbom, A.; Svenningsson, P.; Chen, H.; Wirdefeldt, K. Vagotomy and Parkinson disease. Neurology 2017, 88, 1996–2002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cookson, M.R. The role of leucine-rich repeat kinase 2 (LRRK2) in Parkinson’s disease. Nat. Rev. Neurosci. 2010, 11, 791–797. [Google Scholar] [CrossRef]
  15. Di Maio, R.; Hoffman, E.K.; Rocha, E.M.; Keeney, M.T.; Sanders, L.H.; De Miranda, B.R.; Zharikov, A.; Van Laar, A.; Stepan, A.F.; Lanz, T.A.; et al. LRRK2 activation in idiopathic Parkinson’s disease. Sci. Transl. Med. 2018, 10, 1–13. [Google Scholar] [CrossRef] [Green Version]
  16. Hui, K.Y.; Fernandez-Hernandez, H.; Hu, J.; Schaffner, A.; Pankratz, N.; Hsu, N.Y.; Chuang, L.S.; Carmi, S.; Villaverde, N.; Li, X.; et al. Functional variants in the LRRK2 gene confer shared effects on risk for Crohn’s disease and Parkinson’s disease. Sci. Transl. Med. 2018, 10. [Google Scholar] [CrossRef] [Green Version]
  17. Son, M.Y.; Sim, H.; Son, Y.S.; Jung, K.B.; Lee, M.O.; Oh, J.H.; Chung, S.K.; Jung, C.R.; Kim, J. Distinctive genomic signature of neural and intestinal organoids from familial Parkinson’s disease patient-derived induced pluripotent stem cells. Neuropathol. Appl. Neurobiol. 2017, 43, 584–603. [Google Scholar] [CrossRef]
  18. Han, K.A.; Shin, W.H.; Jung, S.; Seol, W.; Seo, H.; Ko, C.M.; Chung, K.C. Leucine-rich repeat kinase 2 exacerbates neuronal cytotoxicity through phosphorylation of histone deacetylase 3 and histone deacetylation. Hum. Mol. Genet. 2017, 26, 1–18. [Google Scholar] [CrossRef] [Green Version]
  19. Holmqvist, S.; Chutna, O.; Bousset, L.; Aldrin-Kirk, P.; Li, W.; Björklund, T.; Wang, Z.Y.; Roybon, L.; Melki, R.; Li, J.Y. Direct evidence of Parkinson pathology spread from the gastrointestinal tract to the brain in rats. Acta Neuropathol. 2014, 128, 805–820. [Google Scholar] [CrossRef] [Green Version]
  20. Aasly, J.O.; Toft, M.; Fernandez-Mata, I.; Kachergus, J.; Hulihan, M.; White, L.R.; Farrer, M. Clinical features of LRRK2-associated Parkinson’s disease in Central Norway. Ann. Neurol. 2005, 57, 762–765. [Google Scholar] [CrossRef]
  21. Kay, D.M.; Zabetian, C.P.; Factor, S.A.; Nut, J.G.; Samii, A.; Griffith, A.; Brid, T.D.; Kramer, P.; Higgins, D.S.; Payami, H. Parkinson’s disease and LRRK2: Frequency of a common mutation in U.S. movement disorder clinics. Mov. Disord. 2006, 21, 519–523. [Google Scholar] [CrossRef] [PubMed]
  22. Goldwurm, S.; Zini, M.; Di Fonzo, A.; De Gaspari, D.; Siri, C.; Simons, E.J.; van Doeselaar, M.; Tesei, S.; Antonini, A.; Canesi, M.; et al. LRRK2 G2019S mutation and Parkinson’s disease: A clinical, neuropsychological and neuropsychiatric study in a large Italian sample. Park. Relat. Disord. 2006, 12, 410–419. [Google Scholar] [CrossRef] [PubMed]
  23. Adams, J.R.; Van Netten, H.; Schulzer, M.; Mak, E.; Mckenzie, J.; Strongosky, A.; Sossi, V.; Ruth, T.J.; Lee, C.S.; Farrer, M.; et al. PET in LRRK2 mutations: Comparison to sporadic Parkinson’s disease and evidence for presymptomatic compensation. Brain 2005, 128, 2777–2785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Jackson, E.L.; Lu, H. Three-dimensional models for studying development and disease: Moving on from organisms to organs-on-a-chip and organoids. Integr. Biol. 2016, 8, 672–683. [Google Scholar] [CrossRef] [Green Version]
  25. Seok, J.; Warren, H.S.; Cuenca, A.G.; Mindrinos, M.N.; Baker, H.V.; Xu, W.; Richards, D.R.; Mcdonald-Smith, G.P.; Gao, H.; Hennessy, L.; et al. Genomic responses in mouse models poorly mimic human inflammatory diseases, and the Inflammation and Host Response to Injury, Large Scale Collaborative Research Program 4. Proc. Natl. Acad. Sci. USA 2013, 110, 3507–3512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Potashkin, J.A.; Blume, S.R.; Runkle, N.K. Limitations of animal models of Parkinson’s disease. Parkinsons. Dis. 2011, 2011. [Google Scholar] [CrossRef] [Green Version]
  27. Xu, H.; Jiao, Y.; Qin, S.; Zhao, W.; Chu, Q.; Wu, K. Organoid technology in disease modelling, drug development, personalized treatment and regeneration medicine. Exp. Hematol. Oncol. 2018, 7, 1–12. [Google Scholar] [CrossRef] [Green Version]
  28. Horvath, P.; Aulner, N.; Bickle, M.; Davies, A.M.; Del Nery, E.; Ebner, D.; Montoya, M.C.; Östling, P.; Pietiäinen, V.; Price, L.S.; et al. Screening out irrelevant cell-based models of disease. Nat. Rev. Drug Discov. 2016, 15, 751–769. [Google Scholar] [CrossRef] [Green Version]
  29. Fatehullah, A.; Tan, S.H.; Barker, N. Organoids as an in vitro model of human development and disease. Nat. Cell Biol. 2016, 18, 246–254. [Google Scholar] [CrossRef]
  30. Waterston, R.H.; Lindblad-Toh, K.; Birney, E.; Rogers, J.; Abril, J.F.; Agarwal, P.; Agarwala, R.; Ainscough, R.; Alexandersson, M.; An, P.; et al. Initial sequencing and comparative analysis of the mouse genome. Nature 2002, 420, 520–562. [Google Scholar]
  31. Bichler, Z.; Lim, H.C.; Zeng, L.; Tan, E.K. Non-Motor and Motor Features in LRRK2 Transgenic Mice. PLoS ONE 2013, 8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Maekawa, T.; Tsushima, H.; Kawakami, F.; Kawashima, R.; Kodo, M.; Imai, M.; Ichikawa, T. Leucine-rich repeat kinase 2 is associated with activation of the paraventricular nucleus of the hypothalamus and stress-related gastrointestinal dysmotility. Front. Neurosci. 2019, 13, 1–10. [Google Scholar] [CrossRef] [PubMed]
  33. Kim, S.; Kwon, S.H.; Kam, T.I.; Panicker, N.; Karuppagounder, S.S.; Lee, S.; Lee, J.H.; Kim, W.R.; Kook, M.; Foss, C.A.; et al. Transneuronal Propagation of Pathologic α-Synuclein from the Gut to the Brain Models Parkinson’s Disease. Neuron 2019, 103, 627–641.e7. [Google Scholar] [CrossRef] [PubMed]
  34. Matsumoto, K.; Nishihara, S.; Kamimura, M.; Shiraishi, T.; Otoguro, T.; Uehara, M.; Maeda, Y.; Ogura, K.; Lumsden, A.; Ogura, T. The prepattern transcription factor Irx2, a target of the FGF8/MAP kinase cascade, is involved in cerebellum formation. Nat. Neurosci. 2004, 7, 605–612. [Google Scholar] [CrossRef] [PubMed]
  35. Itoh, M.; Kudoh, T.; Dedekian, M.; Kim, C.H.; Chitnis, A.B. A role for iro1 and iro7 in the establishment of an anteposterior compartment of the ectoderm adjacent to the midbrain-hindbrain boundary. Development 2002, 129, 2317–2327. [Google Scholar] [PubMed]
  36. Agarwal, P.; Cheng, C.W.; Kabir, M.G.; Chan, T.Y.; Thanabalasingham, V.; Zhang, X.; Cohen, D.R.; Husain, M.; Cheng, S.H.; Bruneau, B.G.; et al. The Iroquois Homeobox Gene Irx2 Is Not Essential for Normal Development of the Heart and Midbrain-Hindbrain Boundary in Mice Me. Mol. Cell. Biol. 2003, 23, 8216–8225. [Google Scholar]
  37. De La Calle-Mustienes, E.; Feijóo, C.G.; Manzanares, M.; Tena, J.J.; Rodríguez-Seguel, E.; Letizia, A.; Allende, M.L.; Gómez-Skarmeta, J.L. A functional survey of the enhancer activity of conserved non-coding sequences from vertebrate Iroquois cluster gene deserts. Genome Res. 2005, 15, 1061–1072. [Google Scholar] [CrossRef] [Green Version]
  38. Bosse, A.; Zülch, A.; Becker, M.B.; Torres, M.; Gómez-Skarmeta, J.L.; Modolell, J.; Gruss, P. Identification of the vertebrate Iroquois homeobox gene family with overlapping expression during early development of the nervous system. Mech. Dev. 1997, 69, 169–181. [Google Scholar] [CrossRef] [Green Version]
  39. Coetzee, S.G.; Pierce, S.; Brundin, P.; Brundin, L.; Hazelett, D.J.; Coetzee, G.A. Enrichment of risk SNPs in regulatory regions implicate diverse tissues in Parkinson’s disease etiology. Sci. Rep. 2016, 6, 1–11. [Google Scholar] [CrossRef] [Green Version]
  40. Mariani, E.; Frabetti, F.; Tarozzi, A.; Pelleri, M.C.; Pizzetti, F.; Casadei, R. Meta-analysis of Parkinson’s disease transcriptome data using TRAM software: Whole substantia nigra tissue and single dopamine neuron differential gene expression. PLoS ONE 2016, 11, 1–21. [Google Scholar] [CrossRef]
  41. Steger, M.; Tonelli, F.; Ito, G.; Davies, P.; Trost, M.; Vetter, M.; Wachter, S.; Lorentzen, E.; Duddy, G.; Wilson, S.; et al. Phosphoproteomics reveals that Parkinson’s disease kinase LRRK2 regulates a subset of Rab GTPases: (A) Experimental setup of PS1. LRRK2-G2019SGSK mouse embryonic fibroblasts (MEFs, n = 5) were treated with DMSO or each of two structurally distinct LRRK2 inh. Elife 2016, 5, 1–28. [Google Scholar]
  42. Fedele, S.; Collo, G.; Behr, K.; Bischofberger, J.; Müller, S.; Kunath, T.; Christensen, K.; Gündner, A.L.; Graf, M.; Jagasia, R.; et al. Expansion of human midbrain floor plate progenitors from induced pluripotent stem cells increases dopaminergic neuron differentiation potential. Sci. Rep. 2017, 7, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Liu, G.H.; Qu, J.; Suzuki, K.; Nivet, E.; Li, M.; Montserrat, N.; Yi, F.; Xu, X.; Ruiz, S.; Zhang, W.; et al. Progressive degeneration of human neural stem cells caused by pathogenic LRRK2. Nature 2012, 491, 603–607. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Establishment of intestinal organoids as an ex vivo model. (A) Isolated mouse intestinal crypts efficiently formed small intestinal organoids. (scale bar 200 μm) Wild-type (WT), Non-transgenic littermate mouse small intestinal organoids; Leucine-rich repeat kinase 2 (LRRK2) Gly2019Ser substitution (G2019S) (GS), The human LRRK2 G2019S mutant mouse small intestinal organoids. (B) mRNA expression level of Leucine-rich repeat-containing G-protein coupled receptor 5 (Lgr5) (WT, 1.134 ± 0.1343, n = 2; GS, 1.027 ± 0.1325, n = 3, p-value 0.6272), Villin 1 (Vil1) (WT, 1.161 ± 0.1613, n = 2; GS, 1.115 ± 0.8476, n = 3, p-value 0.7929), Lysozyme 1 (Lyz1) (WT, 0.8842 ± 0.1158, n = 2; GS, 1.206 ± 0.2031, n = 3, p-value 0.3275), Chromogranin A (ChgA) (WT, 0.9741 ± 0.02588, n = 2; GS, 1.507 ± 0.3067, n = 3, p-value 0.2712), and Mucin 2 (Muc2) (WT, 1.235 ± 0.2354, n = 2; GS, 1.150 ± 0.2728, n = 3, p-value 0.8427). Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (N.S, not significant). (C) Confocal image for ChgA (enteroendocrine cells), Lysozyme (paneth cells), Mucin2 (goblet cells), and Villin (enterocytes). Counterstain, Hoechst 33342 (scale bar 200 μm).
Figure 1. Establishment of intestinal organoids as an ex vivo model. (A) Isolated mouse intestinal crypts efficiently formed small intestinal organoids. (scale bar 200 μm) Wild-type (WT), Non-transgenic littermate mouse small intestinal organoids; Leucine-rich repeat kinase 2 (LRRK2) Gly2019Ser substitution (G2019S) (GS), The human LRRK2 G2019S mutant mouse small intestinal organoids. (B) mRNA expression level of Leucine-rich repeat-containing G-protein coupled receptor 5 (Lgr5) (WT, 1.134 ± 0.1343, n = 2; GS, 1.027 ± 0.1325, n = 3, p-value 0.6272), Villin 1 (Vil1) (WT, 1.161 ± 0.1613, n = 2; GS, 1.115 ± 0.8476, n = 3, p-value 0.7929), Lysozyme 1 (Lyz1) (WT, 0.8842 ± 0.1158, n = 2; GS, 1.206 ± 0.2031, n = 3, p-value 0.3275), Chromogranin A (ChgA) (WT, 0.9741 ± 0.02588, n = 2; GS, 1.507 ± 0.3067, n = 3, p-value 0.2712), and Mucin 2 (Muc2) (WT, 1.235 ± 0.2354, n = 2; GS, 1.150 ± 0.2728, n = 3, p-value 0.8427). Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (N.S, not significant). (C) Confocal image for ChgA (enteroendocrine cells), Lysozyme (paneth cells), Mucin2 (goblet cells), and Villin (enterocytes). Counterstain, Hoechst 33342 (scale bar 200 μm).
Ijms 21 03455 g001
Figure 2. Comparative analysis of the transcriptome of ex vivo cultured mouse intestinal organoids (mIOs) and in vitro differentiated human pluripotent stem cell (PSC)-derived intestinal organoids (hIOs). (A) Distribution of transcriptome expression ratios in mIOs and hIOs. The Log2 ratio was calculated for each transcript. The frequencies of the Log2 ratio were plotted as a bar chart. (B) Principal component analysis (PCA) of the differentially expressed genes (DEGs) based on the microarray data. (C) Hierarchically clustered heatmap of mIOs and hIOs. Genes with expressions that were not significantly altered were removed. WT, normal control; GS, Parkinson’s disease-relevant cells, carrying LRRK2 G2019S mutation. Both WT1 mIOs and WT#1 hIOs were used as controls.
Figure 2. Comparative analysis of the transcriptome of ex vivo cultured mouse intestinal organoids (mIOs) and in vitro differentiated human pluripotent stem cell (PSC)-derived intestinal organoids (hIOs). (A) Distribution of transcriptome expression ratios in mIOs and hIOs. The Log2 ratio was calculated for each transcript. The frequencies of the Log2 ratio were plotted as a bar chart. (B) Principal component analysis (PCA) of the differentially expressed genes (DEGs) based on the microarray data. (C) Hierarchically clustered heatmap of mIOs and hIOs. Genes with expressions that were not significantly altered were removed. WT, normal control; GS, Parkinson’s disease-relevant cells, carrying LRRK2 G2019S mutation. Both WT1 mIOs and WT#1 hIOs were used as controls.
Ijms 21 03455 g002
Figure 3. Identification and validation of commonly up-regulated or down-regulated factors based on the transcriptome analysis. (A) Venn diagram of the identification of commonly regulated factors between mIOs and hIOs. DEGs were selected on the basis of a >5-fold change (FC). Four factors (Iroquois homeobox 2; IRX2, LIM domain only 3; LMO3, AnnexinA10; ANXA10, and Trefoil factor 2; TFF2) were affiliated to the intersection. (B) Transcriptome analysis data of commonly regulated genes. (C) Quantitative Real-Time PCR (qRT-PCR) verification of commonly regulated genes, including Irx2 (WT, 0.8132 ± 0.1868, n = 2; GS, 3.953 ± 0.4730, n = 3, p-value 0.0153), Lmo3 (WT, 1.414 ± 0.4138, n = 2; GS, 1.526 ± 0.4898, n = 3, p-value 0.8830), Anxa10 (WT, 0.9925 ± 0.007453, n = 2; GS, 0.5974 ± 0.06272, n = 3, p-value 0.0165), and Tff2 (WT, 0.6377 ± 0.3623, n = 2; GS, 0.1783 ± 0.02823, n = 3, p-value 0.1903) in mIOs and IRX2 (WT, 0.5759 ± 0.2178, n = 3; GS, 6.148 ± 1.224, n = 4, p-value 0.0124), LMO3 (WT, 0.5873 ± 0.2296, n = 3; GS, 36.96 ± 2.965, n = 4, p-value 0.0001), ANXA10 (WT, 0.7185 ± 0.1408, n = 3; GS, 0.02325 ± 0.01537, n = 4, p-value 0.0021). and TFF2 (WT, 0.8333 ± 0.1256, n = 3; GS, 0.02061 ± 0.01213, n = 4, p-value 0.0006) in hIOs. FCs in the expression levels were relative to WT1 mIOs and WT#1 hIOs, respectively. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05, ** p < 0.005, *** p < 0.001, N.S, not significant).
Figure 3. Identification and validation of commonly up-regulated or down-regulated factors based on the transcriptome analysis. (A) Venn diagram of the identification of commonly regulated factors between mIOs and hIOs. DEGs were selected on the basis of a >5-fold change (FC). Four factors (Iroquois homeobox 2; IRX2, LIM domain only 3; LMO3, AnnexinA10; ANXA10, and Trefoil factor 2; TFF2) were affiliated to the intersection. (B) Transcriptome analysis data of commonly regulated genes. (C) Quantitative Real-Time PCR (qRT-PCR) verification of commonly regulated genes, including Irx2 (WT, 0.8132 ± 0.1868, n = 2; GS, 3.953 ± 0.4730, n = 3, p-value 0.0153), Lmo3 (WT, 1.414 ± 0.4138, n = 2; GS, 1.526 ± 0.4898, n = 3, p-value 0.8830), Anxa10 (WT, 0.9925 ± 0.007453, n = 2; GS, 0.5974 ± 0.06272, n = 3, p-value 0.0165), and Tff2 (WT, 0.6377 ± 0.3623, n = 2; GS, 0.1783 ± 0.02823, n = 3, p-value 0.1903) in mIOs and IRX2 (WT, 0.5759 ± 0.2178, n = 3; GS, 6.148 ± 1.224, n = 4, p-value 0.0124), LMO3 (WT, 0.5873 ± 0.2296, n = 3; GS, 36.96 ± 2.965, n = 4, p-value 0.0001), ANXA10 (WT, 0.7185 ± 0.1408, n = 3; GS, 0.02325 ± 0.01537, n = 4, p-value 0.0021). and TFF2 (WT, 0.8333 ± 0.1256, n = 3; GS, 0.02061 ± 0.01213, n = 4, p-value 0.0006) in hIOs. FCs in the expression levels were relative to WT1 mIOs and WT#1 hIOs, respectively. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05, ** p < 0.005, *** p < 0.001, N.S, not significant).
Ijms 21 03455 g003
Figure 4. Verification of intestine-based four common factors in human neuroectodermal spheres (NESs). (A) Microarray data of commonly regulated genes in human NESs from a previous study. (B) qRT-PCR confirmation of IRX2 (WT, 0.4329 ± 0.2898, n = 3; GS, 1.628 ± 0.1466, n = 4, p-value 0.0103) expression in human NESs. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05).
Figure 4. Verification of intestine-based four common factors in human neuroectodermal spheres (NESs). (A) Microarray data of commonly regulated genes in human NESs from a previous study. (B) qRT-PCR confirmation of IRX2 (WT, 0.4329 ± 0.2898, n = 3; GS, 1.628 ± 0.1466, n = 4, p-value 0.0103) expression in human NESs. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05).
Ijms 21 03455 g004
Figure 5. Confirmation of the expression of the newly identified gene, IRX2, in neuronal cells containing dopaminergic neurons. (A) Characterization of neuronal cells containing dopaminergic neurons from healthy control (WT), LRRK2 G2019S familial (GS), and sporadic (SPD) PSCs. The neuronal marker (Microtubule-associated protein 2; MAP2) and the dopaminergic neuronal marker (Tyrosine hydroxylase; TH) were co-stained (scale bar 20 μm). (B) Relative expression levels of IRX2 (WT, 0.7100 ± 0.2900, n = 2; PD, 1.598 ± 0.1473, n = 5, p-value 0.0281) in familial/sporadic PD patient-derived neuronal cells containing dopaminergic neurons (PD) in comparison with normal control-derived cells (WT). (C) Determination of IRX2 in human fibroblast cells. Examination of IRX2 (WT, 1.033 ± 0.2582, n = 8; PD, 2.047 ± 0.3931, n = 8, p-value 0.0489) expression in familial/sporadic PD patient-derived fibroblast (PD) and healthy control fibroblast (WT) through qRT-PCR analysis. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05).
Figure 5. Confirmation of the expression of the newly identified gene, IRX2, in neuronal cells containing dopaminergic neurons. (A) Characterization of neuronal cells containing dopaminergic neurons from healthy control (WT), LRRK2 G2019S familial (GS), and sporadic (SPD) PSCs. The neuronal marker (Microtubule-associated protein 2; MAP2) and the dopaminergic neuronal marker (Tyrosine hydroxylase; TH) were co-stained (scale bar 20 μm). (B) Relative expression levels of IRX2 (WT, 0.7100 ± 0.2900, n = 2; PD, 1.598 ± 0.1473, n = 5, p-value 0.0281) in familial/sporadic PD patient-derived neuronal cells containing dopaminergic neurons (PD) in comparison with normal control-derived cells (WT). (C) Determination of IRX2 in human fibroblast cells. Examination of IRX2 (WT, 1.033 ± 0.2582, n = 8; PD, 2.047 ± 0.3931, n = 8, p-value 0.0489) expression in familial/sporadic PD patient-derived fibroblast (PD) and healthy control fibroblast (WT) through qRT-PCR analysis. Data are mean ± standard error of means (SEM) of at least three independent experiments. p-values were analyzed using the unpaired two-tailed Student’s t-test (* p < 0.05).
Ijms 21 03455 g005

Share and Cite

MDPI and ACS Style

Sim, H.; Lee, J.-E.; Yoo, H.M.; Cho, S.; Lee, H.; Baek, A.; Kim, J.; Seo, H.; Kweon, M.-N.; Kim, H.G.; et al. Iroquois Homeobox Protein 2 Identified as a Potential Biomarker for Parkinson’s Disease. Int. J. Mol. Sci. 2020, 21, 3455. https://doi.org/10.3390/ijms21103455

AMA Style

Sim H, Lee J-E, Yoo HM, Cho S, Lee H, Baek A, Kim J, Seo H, Kweon M-N, Kim HG, et al. Iroquois Homeobox Protein 2 Identified as a Potential Biomarker for Parkinson’s Disease. International Journal of Molecular Sciences. 2020; 21(10):3455. https://doi.org/10.3390/ijms21103455

Chicago/Turabian Style

Sim, Hyuna, Joo-Eun Lee, Hee Min Yoo, Sunwha Cho, Hana Lee, Aruem Baek, Jisun Kim, Hyemyung Seo, Mi-Na Kweon, Hyung Gun Kim, and et al. 2020. "Iroquois Homeobox Protein 2 Identified as a Potential Biomarker for Parkinson’s Disease" International Journal of Molecular Sciences 21, no. 10: 3455. https://doi.org/10.3390/ijms21103455

APA Style

Sim, H., Lee, J. -E., Yoo, H. M., Cho, S., Lee, H., Baek, A., Kim, J., Seo, H., Kweon, M. -N., Kim, H. G., Jeon, Y. -J., Son, M. -Y., & Kim, J. (2020). Iroquois Homeobox Protein 2 Identified as a Potential Biomarker for Parkinson’s Disease. International Journal of Molecular Sciences, 21(10), 3455. https://doi.org/10.3390/ijms21103455

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop