Human mutations in high-confidence Tourette disorder genes affect sensorimotor behavior, reward learning, and striatal dopamine in mice.

Tourette disorder (TD) is poorly understood, despite affecting 1/160 children. A lack of animal models possessing construct, face, and predictive validity hinders progress in the field. We used CRISPR/Cas9 genome editing to generate mice with mutations orthologous to human de novo variants in two high-confidence Tourette genes, CELSR3 and WWC1. Mice with human mutations in Celsr3 and Wwc1 exhibit cognitive and/or sensorimotor behavioral phenotypes consistent with TD. Sensorimotor gating deficits, as measured by acoustic prepulse inhibition, occur in both male and female Celsr3 TD models. Wwc1 mice show reduced prepulse inhibition only in females. Repetitive motor behaviors, common to Celsr3 mice and more pronounced in females, include vertical rearing and grooming. Sensorimotor gating deficits and rearing are attenuated by aripiprazole, a partial agonist at dopamine type II receptors. Unsupervised machine learning reveals numerous changes to spontaneous motor behavior and less predictable patterns of movement. Continuous fixed-ratio reinforcement shows that Celsr3 TD mice have enhanced motor responding and reward learning. Electrically evoked striatal dopamine release, tested in one model, is greater. Brain development is otherwise grossly normal without signs of striatal interneuron loss. Altogether, mice expressing human mutations in high-confidence TD genes exhibit face and predictive validity. Reduced prepulse inhibition and repetitive motor behaviors are core behavioral phenotypes and are responsive to aripiprazole. Enhanced reward learning and motor responding occur alongside greater evoked dopamine release. Phenotypes can also vary by sex and show stronger affection in females, an unexpected finding considering males are more frequently affected in TD.

Human mutations in high-confidence Tourette disorder genes affect sensorimotor behavior, reward learning, and striatal dopamine in mice.

Tourette disorder (TD) is poorly understood, despite affecting 1/160 children. A lack of animal models possessing construct, face, and predictive validity hinders progress in the field. We used CRISPR/Cas9 genome editing to generate mice with mutations orthologous to human de novo variants in two high-confidence Tourette genes, CELSR3 and WWC1 . Mice with human mutations in Celsr3 and Wwc1 exhibit cognitive and/or sensorimotor behavioral phenotypes consistent with TD. Sensorimotor gating deficits, as measured by acoustic prepulse inhibition, occur in both male and female Celsr3 TD models. Wwc1 mice show reduced prepulse inhibition only in females. Repetitive motor behaviors, common to Celsr3 mice and more pronounced in females, include vertical rearing and grooming. Sensorimotor gating deficits and rearing are attenuated by aripiprazole, a partial agonist at dopamine type II receptors. Unsupervised machine learning reveals numerous changes to spontaneous motor behavior and less predictable patterns of movement. Continuous fixed-ratio reinforcement shows Celsr3 TD mice have enhanced motor responding and reward learning. Electrically evoked striatal dopamine release, tested in one model, is greater. Brain development is otherwise grossly normal without signs of striatal interneuron loss. Altogether, mice expressing human mutations in high-confidence TD genes exhibit face and predictive validity. Reduced prepulse inhibition and repetitive motor behaviors are core behavioral phenotypes and are responsive to aripiprazole. Enhanced reward learning and motor responding occurs alongside greater evoked dopamine release. Phenotypes can also vary by sex and show stronger affection in females, an unexpected finding considering males are more frequently affected in TD. Significance Statement: We generated mouse models that express mutations in high-confidence genes linked to Tourette disorder (TD). These models show sensorimotor and cognitive behavioral phenotypes resembling TD-like behaviors. Sensorimotor gating deficits and repetitive motor behaviors are attenuated by drugs that act on dopamine. Reward learning and striatal dopamine is enhanced. Brain development is grossly normal, including cortical layering and patterning of major axon tracts. Further, no signs of striatal interneuron loss are detected. Interestingly, behavioral phenotypes in affected females can be more pronounced than in males, despite male sex bias in the diagnosis of TD. These novel mouse models with construct, face, and predictive validity provide a new resource to study neural substrates that cause tics and related behavioral phenotypes in TD.

LD block disorder-specific pleiotropic roles of novel CRHR1 in type 2 diabetes and depression disorder comorbidity.

Major depressive disorder (MDD) and type 2 diabetes (T2D) are complex disorders whose comorbidity can be due to hypercortisolism and may be explained by dysfunction of the corticotropin-releasing hormone receptor 1 (CRHR1) and cortisol feedback within the hypothalamic-pituitary-adrenal axis (HPA axis). To investigate the role of the CRHR1 gene in familial T2D, MDD, and MDD-T2D comorbidity, we tested 152 CRHR1 single-nucleotide-polymorphisms (SNPs), via 2-point parametric linkage and linkage disequilibrium (LD; i.e., association) analyses using 4 models, in 212 peninsular families with T2D and MDD. We detected linkage/LD/association to/with MDD and T2D with 122 (116 novel) SNPs. MDD and T2D had 4 and 3 disorder-specific novel risk LD blocks, respectively, whose risk variants reciprocally confirm one another. Comorbidity was conferred by 3 novel independent SNPs. In silico analyses reported novel functional changes, including the binding site of glucocorticoid receptor-alpha [GR-α] on CRHR1 for transcription regulation. This is the first report of CRHR1 pleiotropic linkage/LD/association with peninsular familial MDD and T2D. CRHR1 contribution to MDD is stronger than to T2D and may antecede T2D onset. Our findings suggest a new molecular-based clinical entity of MDD-T2D and should be replicated in other ethnic groups.

Design, analysis, and interpretation of treatment response heterogeneity in personalized nutrition and obesity treatment research.

It is increasingly assumed that there is no one-size-fits-all approach to dietary recommendations for the management and treatment of chronic diseases such as obesity. This phenomenon that not all individuals respond uniformly to a given treatment has become an area of research interest given the rise of personalized and precision medicine. To conduct, interpret, and disseminate this research rigorously and with scientific accuracy, however, requires an understanding of treatment response heterogeneity. Here, we define treatment response heterogeneity as it relates to clinical trials, provide statistical guidance for measuring treatment response heterogeneity, and highlight study designs that can quantify treatment response heterogeneity in nutrition and obesity research. Our goal is to educate nutrition and obesity researchers in how to correctly identify and consider treatment response heterogeneity when analyzing data and interpreting results, leading to rigorous and accurate advancements in the field of personalized medicine.

Comorbidity of Novel CRHR2 Gene Variants in Type 2 Diabetes and Depression.

The corticotropin-releasing hormone receptor 2 (CRHR2) gene encodes CRHR2, contributing to the hypothalamic-pituitary-adrenal stress response and to hyperglycemia and insulin resistance. CRHR2-/- mice are hypersensitive to stress, and the CRHR2 locus has been linked to type 2 diabetes and depression. While CRHR2 variants confer risk for mood disorders, MDD, and type 2 diabetes, they have not been investigated in familial T2D and MDD. In 212 Italian families with type 2 diabetes and depression, we tested 17 CRHR2 single nucleotide polymorphisms (SNPs), using two-point parametric-linkage and linkage-disequilibrium (i.e., association) analysis (models: dominant-complete-penetrance-D1, dominant-incomplete-penetrance-D2, recessive-complete-penetrance-R1, recessive-incomplete-penetrance-R2). We detected novel linkage/linkage-disequilibrium/association to/with depression (3 SNPs/D1, 2 SNPs/D2, 3 SNPs/R1, 3 SNPs/R2) and type 2 diabetes (3 SNPs/D1, 2 SNPs/D2, 2 SNPs/R1, 1 SNP/R2). All detected risk variants are novel. Two depression-risk variants within one linkage-disequilibrium block replicate each other. Two independent novel SNPs were comorbid while the most significant conferred either depression- or type 2 diabetes-risk. Although the families were primarily ascertained for type 2 diabetes, depression-risk variants showed higher significance than type 2 diabetes-risk variants, implying CRHR2 has a stronger role in depression-risk than type 2 diabetes-risk. In silico analysis predicted variants' dysfunction. CRHR2 is for the first time linked to/in linkage-disequilibrium/association with depression-type 2 diabetes comorbidity and may underlie the shared genetic pathogenesis via pleiotropy.

Childhood maltreatment, serotonin transporter gene, and risk for callous and unemotional traits: A prospective investigation.

Previous studies have reported associations between the serotonin transporter 5-HTTLPR genotype and antisocial and aggressive traits and between child maltreatment and antisocial traits. However, few studies have examined whether 5-HTTLPR moderates the influence of childhood maltreatment on callous and unemotional traits, a hallmark of psychopathy. Using a prospective cohort design, children with documented cases of maltreatment and matched controls were followed up and interviewed in adulthood. DNA was extracted from blood and saliva (N = 414) and callous-unemotional (CU) traits were assessed. Childhood maltreatment predicted higher CU scores in adulthood, whereas the effect of 5-HTTLPR was not significant. The effect of child maltreatment on CU traits did not differ by genetic risk (high or low activity 5-HTTLPR), whereas controls with the LL genotype had higher CU scores than controls with the SS genotype. Similar results were found for females and White, non-Hispanics, but not for males and Blacks. Variations in 5-HTTLPR did not affect the impact of child maltreatment on CU traits in adulthood. Genetic risk had a stronger effect on adults with lower environmental risk (controls). Having a history of child maltreatment or the LL genotype placed participants at risk for higher levels of callous and unemotional trait scores.

A reliable method to determine which candidate chemotherapeutic drugs effectively inhibit tumor growth in patient-derived xenografts (PDX) in single mouse trials.

PURPOSE: We report on a statistical method for grouping anti-cancer drugs (GRAD) in single mouse trials (SMT). The method assigns candidate drugs into groups that inhibit or do not inhibit tumor growth in patient-derived xenografts (PDX). It determines the statistical significance of the group assignments without replicate trials of each drug. METHODS: The GRAD method applies a longitudinal finite mixture model, implemented in the statistical package PROC TRAJ, to analyze a mixture of tumor growth curves for portions of the same tumor in different mice, each single mouse exposed to a different drug. Each drug is classified into an inhibitory or non-inhibitory group. There are several advantages to the GRAD method for SMT. It determines that probability that the grouping is correct, uses the entire longitudinal tumor growth curve data for each drug treatment, can fit different shape growth curves, accounts for missing growth curve data, and accommodates growth curves of different time periods. RESULTS: We analyzed data for 22 drugs for 18 human colorectal tumors provided by researchers in a previous publication. The GRAD method identified 18 drugs that were inhibitory against at least one tumor, and 10 tumors for which there was at least one inhibitory drug. Analysis of simulated data indicated that the GRAD method has a sensitivity of 84% and a specificity of 98%. CONCLUSION: A statistical method, GRAD, can group anti-cancer drugs into those that are inhibitory and those that are non-inhibitory in single mouse trials and provide probabilities that the grouping is correct.

Bone marrow-derived epithelial cells and hair follicle stem cells contribute to development of chronic cutaneous neoplasms.

We used allogeneic bone marrow transplantation (BMT) and a mouse multistage cutaneous carcinogenesis model to probe recruitment of bone marrow-derived epithelial cells (BMDECs) in skin tumors initiated with the carcinogen, dimethylbenz[a]anthracene (DMBA), and promoted with 12-O-tetradecanolyphorbol-13-acetate (TPA). BMDECs clustered in the lesional epithelium, expressed cytokeratins, proliferated, and stratified. We detected cytokeratin induction in plastic-adherent bone marrow cells (BMCs) cultured in the presence of filter-separated keratinocytes (KCs) and bone morphogenetic protein 5 (BMP5). Lineage-depleted BMCs migrated towards High Mobility Group Box 1 (HMGB1) protein and epidermal KCs in ex vivo invasion assays. Naive female mice receiving BMTs from DMBA-treated donors developed benign and malignant lesions after TPA promotion alone. We conclude that BMDECs contribute to the development of papillomas and dysplasia, demonstrating a systemic contribution to these lesions. Furthermore, carcinogen-exposed BMCs can initiate benign and malignant lesions upon tumor promotion. Ultimately, these findings may suggest targets for treatment of non-melanoma skin cancers.

Investigation of previously implicated genetic variants in chronic tic disorders: a transmission disequilibrium test approach.

Genetic studies in Tourette syndrome (TS) are characterized by scattered and poorly replicated findings. We aimed to replicate findings from candidate gene and genome-wide association studies (GWAS). Our cohort included 465 probands with chronic tic disorder (93% TS) and both parents from 412 families (some probands were siblings). We assessed 75 single nucleotide polymorphisms (SNPs) in 465 parent-child trios; 117 additional SNPs in 211 trios; and 4 additional SNPs in 254 trios. We performed SNP and gene-based transmission disequilibrium tests and compared nominally significant SNP results with those from a large independent case-control cohort. After quality control 71 SNPs were available in 371 trios; 112 SNPs in 179 trios; and 3 SNPs in 192 trios. 17 were candidate SNPs implicated in TS and 2 were implicated in obsessive-compulsive disorder (OCD) or autism spectrum disorder (ASD); 142 were tagging SNPs from eight monoamine neurotransmitter-related genes (including dopamine and serotonin); 10 were top SNPs from TS GWAS; and 13 top SNPs from attention-deficit/hyperactivity disorder, OCD, or ASD GWAS. None of the SNPs or genes reached significance after adjustment for multiple testing. We observed nominal significance for the candidate SNPs rs3744161 (TBCD) and rs4565946 (TPH2) and for five tagging SNPs; none of these showed significance in the independent cohort. Also, SLC1A1 in our gene-based analysis and two TS GWAS SNPs showed nominal significance, rs11603305 (intergenic) and rs621942 (PICALM). We found no convincing support for previously implicated genetic polymorphisms. Targeted re-sequencing should fully appreciate the relevance of candidate genes.

An Analytic Solution to the Computation of Power and Sample Size for Genetic Association Studies under a Pleiotropic Mode of Inheritance.

Our motivation here is to calculate the power of 3 statistical tests used when there are genetic traits that operate under a pleiotropic mode of inheritance and when qualitative phenotypes are defined by use of thresholds for the multiple quantitative phenotypes. Specifically, we formulate a multivariate function that provides the probability that an individual has a vector of specific quantitative trait values conditional on having a risk locus genotype, and we apply thresholds to define qualitative phenotypes (affected, unaffected) and compute penetrances and conditional genotype frequencies based on the multivariate function. We extend the analytic power and minimum-sample-size-necessary (MSSN) formulas for 2 categorical data-based tests (genotype, linear trend test [LTT]) of genetic association to the pleiotropic model. We further compare the MSSN of the genotype test and the LTT with that of a multivariate ANOVA (Pillai). We approximate the MSSN for statistics by linear models using a factorial design and ANOVA. With ANOVA decomposition, we determine which factors most significantly change the power/MSSN for all statistics. Finally, we determine which test statistics have the smallest MSSN. In this work, MSSN calculations are for 2 traits (bivariate distributions) only (for illustrative purposes). We note that the calculations may be extended to address any number of traits. Our key findings are that the genotype test usually has lower MSSN requirements than the LTT. More inclusive thresholds (top/bottom 25% vs. top/bottom 10%) have higher sample size requirements. The Pillai test has a much larger MSSN than both the genotype test and the LTT, as a result of sample selection. With these formulas, researchers can specify how many subjects they must collect to localize genes for pleiotropic phenotypes.

Expression and characterization of three Aurora kinase C splice variants found in human oocytes.

Chromosome segregation is an extensively choreographed process yet errors still occur frequently in female meiosis, leading to implantation failure, miscarriage or offspring with developmental disorders. Aurora kinase C (AURKC) is a component of the chromosome passenger complex and is highly expressed in gametes. Studies in mouse oocytes indicate that AURKC is required to regulate chromosome segregation during meiosis I; however, little is known about the functional significance of AURKC in human oocytes. Three splice variants of AURKC exist in testis tissue. To determine which splice variants human oocytes express, we performed quantitative real-time PCR using single oocytes and found expression of all three variants. To evaluate the functional differences between the variants, we created green fluorescent protein-tagged constructs of each variant to express in oocytes from Aurkc(-/-) mice. By quantifying metaphase chromosome alignment, cell cycle progression, phosphorylation of INCENP and microtubule attachments to kinetochores, we found that AURKC_v1 was the most capable of the variants at supporting metaphase I chromosome segregation. AURKC_v3 localized to chromosomes properly and supported cell cycle progression to metaphase II, but its inability to correct erroneous microtubule attachments to kinetochores meant that chromosome segregation was not as accurate compared with the other two variants. Finally, when we expressed the three variants simultaneously, error correction was more robust than when they were expressed on their own. Therefore, oocytes express three variants of AURKC that are not functionally equivalent in supporting meiosis, but fully complement meiosis when expressed simultaneously.

A PAX1 enhancer locus is associated with susceptibility to idiopathic scoliosis in females.

Idiopathic scoliosis (IS) is a common paediatric musculoskeletal disease that displays a strong female bias. By performing a genome-wide association study (GWAS) of 3,102 individuals, we identify significant associations with 20p11.22 SNPs for females (P=6.89 × 10(-9)) but not males (P=0.71). This association with IS is also found in independent female cohorts from the United States of America and Japan (overall P=2.15 × 10(-10), OR=1.30 (rs6137473)). Unexpectedly, the 20p11.22 IS risk alleles were previously associated with protection from early-onset alopecia, another sexually dimorphic condition. The 174-kb associated locus is distal to PAX1, which encodes paired box 1, a transcription factor involved in spine development. We identify a sequence in the associated locus with enhancer activity in zebrafish somitic muscle and spinal cord, an activity that is abolished by IS-associated SNPs. We thus identify a sexually dimorphic IS susceptibility locus, and propose the first functionally defined candidate mutations in an enhancer that may regulate expression in specific spinal cells.

Deficits in Trabecular Bone Microarchitecture in Young Women With Type 1 Diabetes Mellitus.

The pathophysiological mechanism of increased fractures in young adults with type 1 diabetes mellitus (T1DM) is unclear. We conducted a case-control study of trabecular bone microarchitecture and vertebral marrow adiposity in young women with T1DM. Thirty women with T1DM with a median age (range) age of 22.0 years (16.9, 36.1) attending one outpatient clinic with a median age at diagnosis of 9.7 years (0.46, 14.8) were compared with 28 age-matched healthy women who acted as controls. Measurements included MRI-based assessment of proximal tibial bone volume/total volume (appBV/TV), trabecular separation (appTb.Sp), vertebral bone marrow adiposity (BMA), and abdominal adipose tissue and biochemical markers of GH/IGF-1 axis (IGF-1, IGFBP3, ALS) and bone turnover. Median appBV/TV in cases and controls was 0.3 (0.22, 0.37) and 0.33 (0.26, 0.4), respectively (p = 0.018) and median appTb.Sp in T1DM was 2.59 (2.24, 3.38) and 2.32 (2.03, 2.97), respectively (p = 0.012). The median appBV/TV was 0.28 (0.22, 0.33) in those cases with retinopathy (n = 15) compared with 0.33 (0.25, 0.37) in those without retinopathy (p = 0.02). Although median visceral adipose tissue in cases was higher than in controls at 5733 mm(3) (2030, 11,144) and 3460 mm(3) (1808, 6832), respectively (p = 0.012), there was no difference in median BMA, which was 31.1% (9.9, 59.9) and 26.3% (8.5, 49.8) in cases and controls, respectively (p = 0.2). Serum IGF-1 and ALS were also lower in cases, and the latter showed an inverse association to appTbSp (r = -0.30, p = 0.04). Detailed MRI studies in young women with childhood-onset T1DM have shown clear deficits in trabecular microarchitecture of the tibia. Underlying pathophysiological mechanisms may include a microvasculopathy.

Novel cystine ester mimics for the treatment of cystinuria-induced urolithiasis in a knockout mouse model.

OBJECTIVE: To assess the effectiveness of l-cystine dimethyl ester (CDME), an inhibitor of cystine crystal growth, for the treatment of cystine urolithiasis in an Slc3a1 knockout mouse model of cystinuria. MATERIALS AND METHODS: CDME (200 µg per mouse) or water was delivered by gavage daily for 4 weeks. Higher doses by gavage or in the water supply were administered to assess organ toxicity. Urinary amino acids and cystine stones were analyzed to assess drug efficacy using several analytical methods. RESULTS: Treatment with CDME led to a significant decrease in stone size compared with that of the water group (P = .0002), but the number of stones was greater (P = .005). The change in stone size distribution between the 2 groups was evident by micro computed tomography. Overall, cystine excretion in urine was the same between the 2 groups (P = .23), indicating that CDME did not interfere with cystine metabolism. Scanning electron microscopy analysis of cystine stones from the CDME group demonstrated a change in crystal habit, with numerous small crystals. l-cysteine methyl ester was detected by ultra-performance liquid chromatography-mass spectrometer in stones from the CDME group only, indicating that a CDME metabolite was incorporated into the crystal structure. No pathologic changes were observed at the doses tested. CONCLUSION: These data demonstrate that CDME promotes formation of small stones but does not prevent stone formation, consistent with the hypothesis that CDME inhibits cystine crystal growth. Combined with the lack of observed adverse effects, our findings support the use of CDME as a viable treatment for cystine urolithiasis.

Mapping genes with longitudinal phenotypes via Bayesian posterior probabilities.

Most association studies focus on disease risk, with less attention paid to disease progression or severity. These phenotypes require longitudinal data. This paper presents a new method for analyzing longitudinal data to map genes in both population-based and family-based studies. Using simulated systolic blood pressure measurements obtained from Genetic Analysis Workshop 18, we cluster the phenotype data into trajectory subgroups. We then use the Bayesian posterior probability of being in the high subgroup as a quantitative trait in an association analysis with genotype data. This method maintains high power (>80%) in locating genes known to affect the simulated phenotype for most specified significance levels (α). We believe that this method can be useful to aid in the discovery of genes that affect severity or progression of disease.

A meta-analysis identifies adolescent idiopathic scoliosis association with LBX1 locus in multiple ethnic groups.

BACKGROUND: Adolescent idiopathic scoliosis (AIS) is a common rotational deformity of the spine that presents in children worldwide, yet its etiology is poorly understood. Recent genome-wide association studies (GWAS) have identified a few candidate risk loci. One locus near the chromosome 10q24.31 LBX1 gene (OMIM #604255) was originally identified by a GWAS of Japanese subjects and replicated in additional Asian populations. To extend this result, and to create larger AIS cohorts for the purpose of large-scale meta-analyses in multiple ethnicities, we formed a collaborative group called the International Consortium for Scoliosis Genetics (ICSG). METHODS: Here, we report the first ICSG study, a meta-analysis of the LBX1 locus in six Asian and three non-Asian cohorts. RESULTS: We find significant evidence for association of this locus with AIS susceptibility in all nine cohorts. Results for seven cohorts containing both genders yielded P=1.22×10-43 for rs11190870, and P=2.94×10-48 for females in all nine cohorts. Comparing the regional haplotype structures for three populations, we refined the boundaries of association to a ∼25 kb block encompassing the LBX1 gene. The LBX1 protein, a homeobox transcription factor that is orthologous to the Drosophila ladybird late gene, is involved in proper migration of muscle precursor cells, specification of cardiac neural crest cells, and neuronal determination in developing neural tubes. CONCLUSIONS: Our results firmly establish the LBX1 region as the first major susceptibility locus for AIS in Asian and non-Hispanic white groups, and provide a platform for larger studies in additional ancestral groups.

Genetic variants in GPR126 are associated with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is the most common pediatric skeletal disease. We previously reported a locus on chromosome 10q24.31 associated with AIS susceptibility in Japanese using a genome-wide association study (GWAS) consisting of 1,033 cases and 1,473 controls. To identify additional AIS-associated loci, we expanded the study by adding X-chromosome SNPs in the GWAS and increasing the size of the replication cohorts. Through a stepwise association study including 1,819 cases and 25,939 controls, we identified a new susceptibility locus on chromosome 6q24.1 in Japanese (P = 2.25 × 10(-10); odds ratio (OR) = 1.28). The most significantly associated SNP, rs6570507, was in GPR126 (encoding G protein-coupled receptor 126). Its association was replicated in Han Chinese and European-ancestry populations (combined P = 1.27 × 10(-14); OR = 1.27). GPR126 was highly expressed in cartilage, and the knockdown of gpr126 in zebrafish caused delayed ossification of the developing spine. Our results should provide insights into the etiology and pathogenesis of AIS.

A novel method for analyzing genetic association with longitudinal phenotypes.

Knowledge of genes influencing longitudinal patterns may offer information about predicting disease progression. We developed a systematic procedure for testing association between SNP genotypes and longitudinal phenotypes. We evaluated false positive rates and statistical power to localize genes for disease progression. We used genome-wide SNP data from the Framingham Heart Study. With longitudinal data from two real studies unrelated to Framingham, we estimated three trajectory curves from each study. We performed simulations by randomly selecting 500 individuals. In each simulation replicate, we assigned each individual to one of the three trajectory groups based on the underlying hypothesis (null or alternative), and generated corresponding longitudinal data. Individual Bayesian posterior probabilities (BPPs) for belonging to a specific trajectory curve were estimated. These BPPs were treated as a quantitative trait and tested (using the Wald test) for genome-wide association. Empirical false positive rates and power were calculated. Our method maintained the expected false positive rate for all simulation models. Also, our method achieved high empirical power for most simulations. Our work presents a method for disease progression gene mapping. This method is potentially clinically significant as it may allow doctors to predict disease progression based on genotype and determine treatment accordingly.