Trans-ethnic and Ancestry-Specific Blood-Cell Genetics in 746,667 Individuals from 5 Global Populations.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

The Polygenic and Monogenic Basis of Blood Traits and Diseases.

Blood cells play essential roles in human health, underpinning physiological processes such as immunity, oxygen transport, and clotting, which when perturbed cause a significant global health burden. Here we integrate data from UK Biobank and a large-scale international collaborative effort, including data for 563,085 European ancestry participants, and discover 5,106 new genetic variants independently associated with 29 blood cell phenotypes covering a range of variation impacting hematopoiesis. We holistically characterize the genetic architecture of hematopoiesis, assess the relevance of the omnigenic model to blood cell phenotypes, delineate relevant hematopoietic cell states influenced by regulatory genetic variants and gene networks, identify novel splice-altering variants mediating the associations, and assess the polygenic prediction potential for blood traits and clinical disorders at the interface of complex and Mendelian genetics. These results show the power of large-scale blood cell trait GWAS to interrogate clinically meaningful variants across a wide allelic spectrum of human variation.

High-throughput splicing assays identify missense and silent splice-disruptive POU1F1 variants underlying pituitary hormone deficiency.

Pituitary hormone deficiency occurs in ∼1:4,000 live births. Approximately 3% of the cases are due to mutations in the alpha isoform of POU1F1, a pituitary-specific transcriptional activator. We found four separate heterozygous missense variants in unrelated individuals with hypopituitarism that were predicted to affect a minor isoform, POU1F1 beta, which can act as a transcriptional repressor. These variants retain repressor activity, but they shift splicing to favor the expression of the beta isoform, resulting in dominant-negative loss of function. Using a high-throughput splicing reporter assay, we tested 1,070 single-nucleotide variants in POU1F1. We identified 96 splice-disruptive variants, including 14 synonymous variants. In separate cohorts, we found two additional synonymous variants nominated by this screen that co-segregate with hypopituitarism. This study underlines the importance of evaluating the impact of variants on splicing and provides a catalog for interpretation of variants of unknown significance in POU1F1.

Ischemic stroke associated with amyloid-related imaging abnormalities in a patient treated with lecanemab.

INTRODUCTION: Anti-amyloid antibody therapies such as lecanemab are increasingly being used to treat Alzheimer's disease (AD). These therapies are associated with a high rate of amyloid-related imaging abnormalities (ARIA). METHODS: We review the case history of a patient who developed ARIA associated with lecanemab treatment. RESULTS: In addition to microhemorrhages and cerebral edema that are recognized features of ARIA, the patient developed several ischemic strokes. The patient also experienced frequent electrographic seizures without overt clinical seizures. The patient demonstrated clinical and radiographic improvement after steroid treatment. DISCUSSION: Our case suggests that ischemic strokes may be a feature of ARIA and highlights the importance of having a high clinical suspicion for seizures in ARIA. As anti-amyloid therapies are likely going to be increasingly used to treat AD, it is important to appreciate the spectrum of clinical and radiographic findings that can result as side effects from this class of therapies. HIGHLIGHTS: We report a patient who developed severe amyloid-related imaging abnormalities (ARIA) after treatment with lecanemab. Our report suggests that ischemic strokes may be a novel imaging feature of ARIA. Our report highlights the need for high clinical suspicion for seizures in ARIA.

Long-term effect of childhood pandemic experience on medical major choice: Evidence from the 2003 severe acute respiratory syndrome outbreak in China.

This study examines the long-term effect of a pandemic on a crucial human capital decision, namely college major choice. Using China's 2008-2016 major-level National College Entrance Examination (Gaokao) entry grades, we find that the 2003 severe acute respiratory syndrome (SARS) had a substantial deterrent effect on the choice of majoring in medicine among high school graduates who experienced the pandemic in their childhood. In provinces with larger intensities of SARS impact, medical majors become less popular as the average Gaokao grades of enrolled students decline. Further evidence from a nationally representative survey shows that the intensity of the SARS impact significantly decreases children's aspirations to pursue medical occupations, but does not affect their parents' expectations for their children to enter the medical profession. Our discussion on the effect mechanism suggests that the adverse influence of SARS on the popularity of medical majors likely originates from students' childhood experiences.

Burden Testing of Rare Variants Identified through Exome Sequencing via Publicly Available Control Data.

The genetic causes of many Mendelian disorders remain undefined. Factors such as lack of large multiplex families, locus heterogeneity, and incomplete penetrance hamper these efforts for many disorders. Previous work suggests that gene-based burden testing-where the aggregate burden of rare, protein-altering variants in each gene is compared between case and control subjects-might overcome some of these limitations. The increasing availability of large-scale public sequencing databases such as Genome Aggregation Database (gnomAD) can enable burden testing using these databases as controls, obviating the need for additional control sequencing for each study. However, there exist various challenges with using public databases as controls, including lack of individual-level data, differences in ancestry, and differences in sequencing platforms and data processing. To illustrate the approach of using public data as controls, we analyzed whole-exome sequencing data from 393 individuals with idiopathic hypogonadotropic hypogonadism (IHH), a rare disorder with significant locus heterogeneity and incomplete penetrance against control subjects from gnomAD (n = 123,136). We leveraged presumably benign synonymous variants to calibrate our approach. Through iterative analyses, we systematically addressed and overcame various sources of artifact that can arise when using public control data. In particular, we introduce an approach for highly adaptable variant quality filtering that leads to well-calibrated results. Our approach "re-discovered" genes previously implicated in IHH (FGFR1, TACR3, GNRHR). Furthermore, we identified a significant burden in TYRO3, a gene implicated in hypogonadotropic hypogonadism in mice. Finally, we developed a user-friendly software package TRAPD (Test Rare vAriants with Public Data) for performing gene-based burden testing against public databases.

The Genetic Landscape of Diamond-Blackfan Anemia.

Diamond-Blackfan anemia (DBA) is a rare bone marrow failure disorder that affects 7 out of 1,000,000 live births and has been associated with mutations in components of the ribosome. In order to characterize the genetic landscape of this heterogeneous disorder, we recruited a cohort of 472 individuals with a clinical diagnosis of DBA and performed whole-exome sequencing (WES). We identified relevant rare and predicted damaging mutations for 78% of individuals. The majority of mutations were singletons, absent from population databases, predicted to cause loss of function, and located in 1 of 19 previously reported ribosomal protein (RP)-encoding genes. Using exon coverage estimates, we identified and validated 31 deletions in RP genes. We also observed an enrichment for extended splice site mutations and validated their diverse effects using RNA sequencing in cell lines obtained from individuals with DBA. Leveraging the size of our cohort, we observed robust genotype-phenotype associations with congenital abnormalities and treatment outcomes. We further identified rare mutations in seven previously unreported RP genes that may cause DBA, as well as several distinct disorders that appear to phenocopy DBA, including nine individuals with biallelic CECR1 mutations that result in deficiency of ADA2. However, no new genes were identified at exome-wide significance, suggesting that there are no unidentified genes containing mutations readily identified by WES that explain >5% of DBA-affected case subjects. Overall, this report should inform not only clinical practice for DBA-affected individuals, but also the design and analysis of rare variant studies for heterogeneous Mendelian disorders.

Temporal trends and yield of clinical diagnostic genetic testing in adult neurology.

While genetics evaluation is increasingly utilized in adult neurology patients, its usage and efficacy are not well characterized. Here, we report our experience with 1461 consecutive patients evaluated in an adult neurogenetics clinic at a large academic medical center between January 2015 and March 2020. Of the 1461 patients evaluated, 1215 patients were referred for the purposes of identifying a genetic diagnosis for an undiagnosed condition, 90.5% of whom underwent genetic testing. The modalities of genetic testing utilized varied across referral diagnostic categories, including a range of utilization of whole exome sequencing (WES) as an initial test in 13.9% of neuromuscular patients to 52.9% in white matter disorder patients. The usage of WES increased over time, from 7.7% of initial testing in 2015 to a peak of 27.3% in 2019. Overall, genetic testing yielded a causal genetic diagnosis in 30.7% of patients. This yield was higher in certain referring diagnosis categories, such as neuromuscular (39.0%) and epilepsy (29.8%). Our study demonstrates that evaluation at an adult neurogenetics referral center can yield diagnoses in a substantial fraction of patients. Additional research will be needed to determine optimal genetic testing strategies and cost effectiveness of adult neurogenetics evaluation.

Risk of sexual dysfunction with progestin-based contraceptives in women of child-bearing age.

PURPOSE: To examine whether there is a positive association between sexual dysfunction (SD) and different types of progestin-based contraceptives. METHODS: Nested case-control study in women of child-bearing age (15-45 years) from the IQVIA® Ambulatory electronic medical record database from 2008 to 2018. Cases defined by diagnosis of sexual dysfunction identified by international classification for disease clinical modification code 9th and 10th. Each case was matched to four controls and rates of prescriptions of the following were compared: levonorgestrel intra-uterine device (IUD), progestin, and ethinyl estradiol (EE) combined oral contraceptive (COC) formulations including levonorgestrel, norgestimate, drospirenone, desogestrel, norethindrone, and norgestrel; etonogestrel vaginal ring; and medroxyprogesterone injection. RESULTS: Overall, 6689 cases of patients with SD were matched to 26,756 matched controls. Compared with matched controls, more subjects with SD used levonorgestrel IUD (OR 1.24, 95% CI 1.08-1.44), EE-levonorgestrel COC (OR 1.18, 95% CI 1.00-1.41), EE-drospirenone (OR 1.28, 95% CI 1.00-1.67), and medroxyprogesterone (OR 1.38, 95% CI 1.12-1.70). The use of norgestrel exhibited a protective effect (OR 0.83, 95% CI 0.73-0.95). When using the EE-levonorgestrel COC as a comparator, norgestrel users exhibited a protective effect (OR 0.70, 95% CI 0.57-0.87) while no other contraceptives showed a statistically significant difference in association with SD. CONCLUSION: Our study found an increase in the use of levonorgestrel (COC and IUD), drospirenone, and medroxyprogesterone in subjects with SD. The risk of contraceptives did not differ when compared with oral levonorgestrel. The small association size and lack of difference between drug formulations suggest a minimal impact of progestin-based contraceptives on sexual dysfunction.

Perioperative cardiac investigations for chest pain after parathyroidectomy rarely yield a cardiac diagnosis.

BACKGROUND: The incidence of adverse perioperative cardiac complications after parathyroidectomy has not been well described. This study aimed to evaluate the incidence of perioperative chest pain and cardiac complications after parathyroidectomy and to evaluate risk factors that may identify patients who are more likely to benefit from a cardiac workup. METHODS: We performed a retrospective study of all patients undergoing parathyroidectomy for primary hyperparathyroidism by a single endocrine surgeon at a tertiary endocrine centre between 2011 and 2018. Patient demographics, clinicopathologic variables, operative and postoperative details (reported chest pain, performance of a cardiac workup and new postoperative cardiac diagnosis) were reviewed. Patients with chest pain were compared to those without chest pain using the Fisher exact test and Student t test. RESULTS: Fourteen of 295 patients (4.7%) reported chest pain in the immediate postoperative period. Most patients were investigated with a 12-lead electrocardiogram and troponin (n = 12/14), yet none were diagnosed with a cardiac event. When comparing patients with and without chest pain, there was no significant difference in age, gender, body mass index, presence of cardiovascular risk factors, American Society of Anesthesiologists score or length of surgery. CONCLUSION: Postoperative chest pain after parathyroidectomy is not an uncommon event and leads to a cardiac workup in most cases; however, the risk of significant postoperative cardiac events is minimal. In the "choosing wisely" era, one should evaluate each patient's pretest probability of such events and avoid extensive workup in low-risk patients to avoid unnecessary costs to the health care system.

Development and Prospective Validation of a Risk Calculator That Predicts a Low Risk Cohort for Atypical Ductal Hyperplasia Upstaging to Malignancy: Evidence for a Watch and Wait Strategy of a High-Risk Lesion.

BACKGROUND: Guidelines recommend surgical excision of atypical ductal hyperplasia (ADH) due to the concern of undersampling a potential malignancy on core needle biopsy (CNB). The purpose of this study was to determine clinical, radiological and pathological variables associated with ADH upstaging to cancer and to develop a predictive risk calculator capable of identifying women who have a low oncological risk of upstaging. METHODS: A prospectively collected database from a tertiary breast referral center was analyzed for women diagnosed with ADH on CNB between January 2013 to December 2017 who underwent surgical excision. CNB and surgical pathology reports were examined to determine rate of upstaging. The association between clinical, radiological and pathological variables were evaluated using regression analysis to determine predictors of ADH upstaging to cancer. Significant variables (p ≤ 0.05) identified on univariate analysis were assigned a score of "1" and were included in the ADH upstaging risk calculator. RESULTS: A total of 1986 patients underwent surgery for a high-risk lesion. We identified 318 (16.0%) patients who had ADH identified on their CNB who underwent surgery-of which 290 were included in our study. The upstage rate was 24.8%. Five variables were associated with upstaging and included in our calculator: (1) lesion > 5 mm on ultrasound; (2) lesion > 5 mm on mammogram; (3) one or more "high-risk" lesion(s) on CNB; (4) pathological suspicion for cancer and; (5) incomplete removal of calcifications on CNB. Patients with a score of 0 had a 2% risk of being upstaged to cancer and were deemed low risk with 17.2% of patients falling within this category. CONCLUSIONS: Patients with ADH on CNB can be stratified into a low oncological cohort who have a 2% risk of being upstaged to carcinoma. In the future, these select patients may be counselled and potentially offered observation as an alternative to surgery.

TMEM106B Effect on cognition in Parkinson disease and frontotemporal dementia.

OBJECTIVE: Common variants near TMEM106B associate with risk of developing frontotemporal dementia (FTD). Emerging evidence suggests a role for TMEM106B in neurodegenerative processes beyond FTD. We evaluate the effect of TMEM106B genotype on cognitive decline across multiple neurogenerative diseases. METHODS: We longitudinally followed 870 subjects with diagnoses of Parkinson disease (PD; n = 179), FTD (n = 179), Alzheimer disease (AD; n = 300), memory-predominant mild cognitive impairment (MCI; n = 75), or neurologically normal control subjects (NC; n = 137) at the University of Pennsylvania (UPenn). All participants had annual Mini-Mental State Examination (MMSE; median follow-up duration = 3.0 years) and were genotyped at TMEM106B index single nucleotide polymorphism rs1990622. Genotype effects on cognition were confirmed by extending analyses to additional cognitive instruments (Mattis Dementia Rating Scale-2 [DRS-2] and Montreal Cognitive Assessment [MoCA]) and to an international validation cohort (Parkinson's Progression Markers Initiative [PPMI], N = 371). RESULTS: The TMEM106B rs1990622T allele, linked to increased risk of FTD, associated with greater MMSE decline over time in PD subjects but not in AD or MCI subjects. For FTD subjects, rs1990622T associated with more rapid decrease in MMSE only under the minor-allele, rs1990622C , dominant model. Among PD patients, rs1990622T carriers from the UPenn cohort demonstrated more rapid longitudinal decline in DRS-2 scores. Finally, in the PPMI cohort, TMEM106B risk allele carriers demonstrated more rapid longitudinal decline in MoCA scores. INTERPRETATION: Irrespective of cognitive instrument or cohort assessed, TMEM106B acts as a genetic modifier for cognitive trajectory in PD. Our results implicate lysosomal dysfunction in the pathogenesis of cognitive decline in 2 different proteinopathies. ANN NEUROL 2019;85:801-811.

Genetics of rheumatoid arthritis contributes to biology and drug discovery.

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological data sets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA). Here we performed a genome-wide association study meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ∼10 million single-nucleotide polymorphisms. We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 101 (refs 2 - 4). We devised an in silico pipeline using established bioinformatics methods based on functional annotation, cis-acting expression quantitative trait loci and pathway analyses--as well as novel methods based on genetic overlap with human primary immunodeficiency, haematological cancer somatic mutations and knockout mouse phenotypes--to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery.

Comprehensive population-based genome sequencing provides insight into hematopoietic regulatory mechanisms.

Genetic variants affecting hematopoiesis can influence commonly measured blood cell traits. To identify factors that affect hematopoiesis, we performed association studies for blood cell traits in the population-based Estonian Biobank using high-coverage whole-genome sequencing (WGS) in 2,284 samples and SNP genotyping in an additional 14,904 samples. Using up to 7,134 samples with available phenotype data, our analyses identified 17 associations across 14 blood cell traits. Integration of WGS-based fine-mapping and complementary epigenomic datasets provided evidence for causal mechanisms at several loci, including at a previously undiscovered basophil count-associated locus near the master hematopoietic transcription factor CEBPA The fine-mapped variant at this basophil count association near CEBPA overlapped an enhancer active in common myeloid progenitors and influenced its activity. In situ perturbation of this enhancer by CRISPR/Cas9 mutagenesis in hematopoietic stem and progenitor cells demonstrated that it is necessary for and specifically regulates CEBPA expression during basophil differentiation. We additionally identified basophil count-associated variation at another more pleiotropic myeloid enhancer near GATA2, highlighting regulatory mechanisms for ordered expression of master hematopoietic regulators during lineage specification. Our study illustrates how population-based genetic studies can provide key insights into poorly understood cell differentiation processes of considerable physiologic relevance.

Determinants of Power in Gene-Based Burden Testing for Monogenic Disorders.

Whole-exome sequencing has enabled new approaches for discovering genes associated with monogenic disorders. One such approach is gene-based burden testing, in which the aggregate frequency of "qualifying variants" is compared between case and control subjects for each gene. Despite substantial successes of this approach, the genetic causes for many monogenic disorders remain unknown or only partially known. It is possible that particular genetic architectures lower rates of discovery, but the influence of these factors on power has not been rigorously evaluated. Here, we leverage large-scale exome-sequencing data to create an empirically based simulation framework to evaluate the impact of key parameters (background variation rates, locus heterogeneity, mode of inheritance, penetrance) on power in gene-based burden tests in the context of monogenic disorders. Our results demonstrate that across genes, there is a wide range in sample sizes needed to achieve power due to differences in the background rate of rare variants in each gene. Increasing locus heterogeneity results in rapid increases in sample sizes needed to achieve adequate power, particularly when individual genes contribute to less than 5% of cases under a dominant model. Interestingly, incomplete penetrance as low as 10% had little effect on power due to the low prevalence of monogenic disorders. Our results suggest that moderate incomplete penetrance is not an obstacle in this gene-based burden testing approach but that dominant disorders with high locus heterogeneity will require large sample sizes. Our simulations also provide guidance on sample size needs and inform study design under various genetic architectures.

Estimating yields of prenatal carrier screening and implications for design of expanded carrier screening panels.

PURPOSE: Prenatal genetic carrier screening can identify parents at risk of having a child affected by a recessive condition. However, the conditions/genes most appropriate for screening remain a matter of debate. Estimates of carrier rates across genes are needed to guide construction of carrier screening panels. METHOD: We leveraged an exome sequencing database (n = 123,136) to estimate carrier rates across six major ancestries for 415 genes associated with severe recessive conditions. RESULTS: We found that 32.6% (East Asian) to 62.9% (Ashkenazi Jewish) of individuals are variant carriers in at least one of the 415 genes. For couples, screening all 415 genes would identify 0.17-2.52% of couples as being at risk for having a child affected by one of these conditions. Screening just the 40 genes with carrier rate >1.0% would identify more than 76% of these at-risk couples. An ancestry-specific panel designed to capture genes with carrier rates >1.0% would include 5 to 28 genes, while a comparable panethnic panel would include 40 genes. CONCLUSION: Our work guides the design of carrier screening panels and provides data to assist in counseling prospective parents. Our results highlight a high cumulative carrier rate across genes, underscoring the need for careful selection of genes for screening.

A Novel Olfactometer for Efficient and Flexible Odorant Delivery.

Understanding how sensory space maps to neural activity in the olfactory system requires efficiently and flexibly delivering numerous odorants within single experimental preparations. Such delivery is difficult with current olfactometer designs, which typically include limited numbers of stimulus channels and are subject to intertrial and interchannel contamination of odorants. Here, we present a novel olfactometer design that is easily constructed, modular, and capable of delivering an unlimited number of odorants in air with temporal precision and no detectable intertrial or interchannel contamination. The olfactometer further allows for the flexible generation of odorant mixtures and flexible timing of odorant sequences. Odorant delivery from the olfactometer is turbulent but reliable from trial to trial, supporting operant conditioning of mice in an odorant discrimination task and permitting odorants and concentrations to be mapped to neural activity with a level of precision equivalent to that obtained with a flow dilution olfactometer. This novel design thus provides several unique advantages for interrogating olfactory perception and for mapping sensory space to neural activity in the olfactory system.

Genome-wide Analysis of Body Proportion Classifies Height-Associated Variants by Mechanism of Action and Implicates Genes Important for Skeletal Development.

Human height is a composite measurement, reflecting the sum of leg, spine, and head lengths. Many common variants influence total height, but the effects of these or other variants on the components of height (body proportion) remain largely unknown. We studied sitting height ratio (SHR), the ratio of sitting height to total height, to identify such effects in 3,545 African Americans and 21,590 individuals of European ancestry. We found that SHR is heritable: 26% and 39% of the total variance of SHR can be explained by common variants in European and African Americans, respectively, and global European admixture is negatively correlated with SHR in African Americans (r(2) ≈ 0.03). Six regions reached genome-wide significance (p < 5 × 10(-8)) for association with SHR and overlapped biological candidate genes, including TBX2 and IGFBP3. We found that 130 of 670 height-associated variants are nominally associated (p < 0.05) with SHR, more than expected by chance (p = 5 × 10(-40)). At these 130 loci, the height-increasing alleles are associated with either a decrease (71 loci) or increase (59 loci) in SHR, suggesting that different height loci disproportionally affect either leg length or spine/head length. Pathway analyses via DEPICT revealed that height loci affecting SHR, and especially those affecting leg length, show enrichment of different biological pathways (e.g., bone/cartilage/growth plate pathways) than do loci with no effect on SHR (e.g., embryonic development). These results highlight the value of using a pair of related but orthogonal phenotypes, in this case SHR with height, as a prism to dissect the biology underlying genetic associations in polygenic traits and diseases.

Association of Antidepressant Use With Drug-Related Extrapyramidal Symptoms: A Pharmacoepidemiological Study.

BACKGROUND: Antidepressants are one of the most prescribed classes of medications. A number of case reports have linked these drugs to extrapyramidal symptoms (EPSs), but no large epidemiologic study to date has examined this association. We sought to quantify the association of EPSs with different antidepressants by undertaking a large pharmacoepidemiologic study. METHODS: A nested case-control study was conducted using a large health claims database in the United States from June 2006 to December 2015. Subjects with a diagnosis of primary Parkinson disease and those who received prescriptions of levodopa, ropinirole, pramipexole, domperidone, metoclopramide, entacapone, benztropine, selegiline, rasagiline, diphenhydramine, trihexyphenidyl, typical and atypical antipsychotics, and tricyclic antidepressants were excluded. Cases were followed to the first billing code for an extrapyramidal event or last date of enrollment in the cohort. For each case, 10 control subjects were matched by follow-up time, calendar time, and age through density-based sampling. Rate ratios were computed using conditional logistic regression adjusting for other covariates. RESULTS: We identified 3,838 subjects with EPSs compared with 38,380 age-matched control subjects. Rate ratios with respect to EPSs were as follows: duloxetine, 5.68 (95% confidence interval [CI], 4.29-7.53); mirtazapine, 3.78 (95% CI, 1.71-8.32); citalopram, 3.47 (95% CI, 2.68-4.50); escitalopram, 3.23 (95% CI, 2.44-4.26); paroxetine, 3.07 (95% CI, 2.15-4.40); sertraline, 2.57 (95% CI, 2.02-3.28); venlafaxine, 2.37 (95% CI, 1.71-3.29); bupropion, 2.31 (95% CI, 1.67-3.21); and fluoxetine, 2.03 (95% CI, 1.48-2.78). CONCLUSIONS: This observational study demonstrates a harmful association between the incidence of Parkinson disease or associated EPSs and use of the antidepressants duloxetine, mirtazapine, citalopram, escitalopram, paroxetine, sertraline, venlafaxine, bupropion, and fluoxetine.