Skip to main content

Enhanced rare disease mapping for phenome-wide genetic association in the UK Biobank

Abstract

Background

Rare diseases collectively affect up to 10% of the population, but often lack effective treatment, and typically little is known about their pathophysiology. Major challenges include suboptimal phenotype mapping and limited statistical power. Population biobanks, such as the UK Biobank, recruit many individuals who can be affected by rare diseases; however, investigation into their utility for rare disease research remains limited. We hypothesized the UK Biobank can be used as a unique population assay for rare diseases in the general population.

Methods

We constructed a consensus mapping between ICD-10 codes and ORPHA codes for rare diseases, then identified individuals with each rare condition in the UK Biobank, and investigated their age at recruitment, sex bias, and comorbidity distributions. Using exome sequencing data from 167,246 individuals of European ancestry, we performed genetic association controlling for case/control imbalance (SAIGE) to identify potential rare pathogenic variants for each disease.

Results

Using our mapping approach, we identified and characterized 420 rare diseases affecting 23,575 individuals in the UK Biobank. Significant genetic associations included JAK2 V617F for immune thrombocytopenic purpura (p = 1.24 × 10−13) and a novel CALR loss of function variant for essential thrombocythemia (p = 1.59 × 10−13). We constructed an interactive resource highlighting demographic information (http://www-personal.umich.edu/~mattpat/rareDiseases.html) and demonstrate transferability by applying our mapping to a medical claims database.

Conclusions

Enhanced disease mapping and increased power from population biobanks can elucidate the demographics and genetic associations for rare diseases.

Background

Rare diseases account for a high proportion of hospital visits [1, 2] and can severely reduce quality of life [3, 4]. A recent study [1] indicates ~10% of hospital discharges were patients who had a rare disease. Many rare diseases have a substantial impact on mortality [2, 5] and are associated with severe physical and mental disability [6]. Legislation has been introduced (e.g., the Rare Diseases Act and the Orphan Drug Act in the USA [7]) to encourage research and drug development for rare diseases, yet < 5% have an FDA approved treatment [5], and available treatments are often highly expensive [8] or inaccessible [9]. Delays in diagnosis are commonplace [3] and patients reported challenges in accessing specialists with experience of their condition [4]. Inference from Orphanet suggests 80% of rare diseases are genetic, yet disease-associated genes have been identified for < 40% of rare diseases [5]. A clearer understanding of rare disease pathogenesis is needed, to develop new treatments and reduce the burden on society.

Addressing challenges in healthcare for rare diseases is challenging due to the vast number and heterogeneity of rare conditions, as well as limited resources available to study each disease. There are believed to be > 10,000 rare diseases [10], with estimates varying due to divergent definitions of rarity [5] and lack of consensus over which diseases constitute distinct entities [11]. It can be difficult to collect sufficient samples for statistical significance, due to low prevalence [12], and clinical trials for rare diseases typically include fewer participants [13]. There are ongoing attempts to establish biobanks for specific rare diseases, such as the EuroBioBank network [14], which provide access to samples for research, and information to guide clinical trials. However, biobanks are not available for every disease, and they cannot be used to understand demographics of rare diseases in the general population.

Population-based biobanks, such as the UK Biobank (UKB) [15] and All of Us research program [16], provide extensive information (both phenotypes and genotypes) on individuals with a wide range of (rare and common) diseases, reflecting the general population. The overall prevalence of having at least one rare disease has been estimated at > 3.5% in Europeans [17], which is higher than many common diseases, such as psoriasis [18] and glaucoma [19]; some sources suggest as many as 1 in 10 Americans are living with a rare disease [4, 10]. One difficulty in harnessing population biobanks is their diagnosis coding systems (e.g., the 10th revision of the International Statistical Classification of Diseases and Related Health Problems, ICD-10) are not specifically designed for rare diseases. ICD-10 codes are typically broader than those from rare disease nomenclatures (such as Orphanet [20]), meaning that while it is possible to map from ORPHA codes (in Orphanet) to ICD-10 for billing purposes, individuals annotated with those ICD-10 codes do not necessarily have a rare disease.

We hypothesize that applying enhanced phenotype mapping to the large sample sizes of population biobanks can provide valuable information regarding the demographics of rare diseases and could help advance the study of previously underserved conditions. We perform a consensus mapping of ICD-10 to ORPHA codes, to identify which ICD-10 codes reliably indicate rare diseases, then use it to catalog and explore the data available on rare diseases through the UKB. To facilitate the dissemination of this information, we have created an online resource (http://www-personal.umich.edu/~mattpat/rareDiseases.html) which will enable researchers to investigate specific rare diseases in the UKB and other systems that rely on ICD-10 codes. Furthermore, our association studies highlight significant genetic-disease associations, shedding light on the importance of refined phenotype mapping.

Methods

Mapping ICD-10 to ORPHA codes

From the ICD-10 mapping provided by Orphanet (which does not require the ICD-10 codes to be as specific as their corresponding ORPHA codes), we conducted refined mapping such that an individual annotated with an ICD-10 code can be expected to have the rare disease indicated by the corresponding ORPHA code (Supplementary Figure 1). For example, Orphanet maps ORPHA:314 (Leiner’s disease) to L21.1 (seborrheic infantile dermatitis); however, seborrheic dermatitis is in general common among newborns [21], so this ICD-10 code is excluded by our approach. We thus required the ORPHA code selected for each ICD-10 code to be as specific as possible (i.e., provide the most narrow and precise description), but no more specific than the ICD-10 code, such that we are able to reliably map ICD-10 codes to ORPHA codes representing specific rare diseases. No parent ICD-10 codes of the mapped ICD-10 codes are used in our mapping.

Clinical and demographic data, from 502,493 individuals in the UKB, was downloaded on 24 August 2020, with diagnoses annotated by ICD-10 codes. We started with 2,044 ICD-10 codes from the UKB mapped to 6762 ORPHA codes by Orphanet and applied a consensus approach to identify the subset of mappings that meet our criteria, selecting a single ORPHA code for each ICD-10 code (Additional file 1: Table S1). It should be noted that multiple ICD-10 codes can map to the same ORPHA code if the ICD-10 codes represent subtypes of the rare disease; for example, multisystemic (C96.0), unisystemic (C96.5), and unifocal (C96.6) Langerhans cell histiocytoses are all mapped to ORPHA:389 (Langerhans cell histiocytosis). Putative ICD-10/ORPHA pairs were assessed independently (by RB and MTP), and then the assessments were compared together. As the third voter in our consensus, we employed the results of a previous study mapping a subset of the Australian modification of ICD-10 (ICD-10-AM) to ORPHA codes [1]. Where assessments differed, further investigation was conducted until a joint decision could be reached. The end result is a set of ICD-10 codes for each ORPHA code that can be reliably mapped, allowing rare diseases to be identified in the UKB and other resources that use ICD-10 codes.

UKB data analysis

We extracted the subset of individuals with ICD-10 codes from primary or secondary diagnoses corresponding to each ORPHA code in our mapping for rare diseases. Disease prevalence was estimated by dividing the number of individuals with ICD-10 codes mapping to those ORPHA codes by the total number of individuals in the UKB (of all ancestries). We compared the estimated prevalence for each rare disease with criteria for considering a disease as rare: < 1 in 2,000 people in Europe (< 0.05%) or < 200,000 people in the USA (< 0.06%) [5]. Statistics describing the proportion of males, age at recruitment, and comorbidities (collected from the entire study period) were compiled for each disease from the available data in the UKB. Diseases were grouped based on the different chapters in the ICD-10 coding system, splitting off “Other Immune” from Chapter III for D80-D89, then combining Chapters V and VI as neurological diseases, VII and VIII as Eye/Ear, chapters XV and XVI as pregnancy/childbirth, and XVIII, XX, XXI, and XXII as Miscellaneous. Access to data from the UK Biobank (https://www.ukbiobank.ac.uk/learn-more-about-uk-biobank/about-us/ethics) was obtained through a material transfer agreement which falls within the UK Biobank’s generic Research Tissue Bank (RTB) approval from the NHS North West REC.

Comorbidities

We investigated the prevalence of different groups of comorbidities among individuals with different groups of rare diseases using the groups we identified from the ICD-10 codes. We also compared the enrichment of these groups of comorbidities against two previous studies (one for common diseases [22] and one for Mendelian diseases [23]) using Fisher exact tests.

Exome sequencing analysis

We performed genetic association analysis on 167,246 individuals of European ancestry from the UKB who self-reported as White British and were genetically confirmed by principal components analysis (UKB field 22,006), using exome sequencing data prepared with the Original Quality Functionally Equivalent (OQFE) protocol. We restricted the genetic analysis to individuals of European ancestry because they constitute the vast majority of the UK Biobank and imbalanced cases/controls for rare diseases could otherwise lead to false positives for variants associated with ancestry. Variant- and gene-level association tests were conducted using SAIGE/SAIGE-GENE [24, 25], a generalized mixed model that provides accurate results in situations where there are far fewer cases than controls and applies a sparse relatedness matrix to account for population stratification. Variants were required to be rare (minor allele frequency, MAF ≤ 1%) in the full population but have a minor allele count of at least three among individuals with the rare disease and more than three overall, annotated as having a high or moderate impact by SnpEff [26] or a score ≥ 0.75 by REVEL [27], resulting in 33,981 variants across 162 rare diseases. Significant variant and gene associations from SAIGE were compared against the genes annotated for each disease by Orphanet and nine different variant effect predictors (CADD, FATHMM, LRTori, MetaSVM, MutationAssessor, MutationTaster Polyphen2 using HDIV or HVAR, PROVEAN and SIFT) annotated in dbNSFP v4.1 [28]. We then evaluated the pathogenicity of variants identified using the 2015 ACMG/AMP classification [29]. We also investigated loss of function variants predicted using SnpEff [26] by filtering the SAIGE associations to those variants. Significant rare disease associations were further filtered to highlight diseases which, in addition to being indicated by Orphanet as rare in Europe (< 1 in 2000 people) are indicated by NIH’s Genetic and Rare Diseases Information Center (GARD) as being rare in the USA (< 200,000 people).

Confirmation using other data sources

To illustrate how our ICD-10/ORPHA mappings can facilitate rare disease research in additional datasets, we applied them to a nationwide database of medical claims from 39 million patients (Optum’s deidentified Clinformatics® Data Mart [30]) and compared the estimated prevalence of each disease against the UKB. Since a large proportion of diagnoses in Clinformatics® are coded using ICD-9, it is necessary to find equivalent ICD-9 codes for the ICD-10/ORPHA mappings. However, the ICD-9 coding system is less detailed than ICD-10, particularly for rare diseases. We therefore used the 2018 General Equivalence Mappings (GEMS) published by the Centers for Medicare and Medicaid Services [31] and excluded any diseases labeled as having approximate mappings. Although this limits the number of diseases that can be mapped, it allows for more accurate disease identification. Access to data from the Clinformatics® Data Mart (https://www.optum.com/content/dam/optum/resources/productSheets/Clinformatics_for_Data_Mart.pdf) was obtained through a data use agreement according to a license agreement at the University of Michigan Institute for Healthcare Policy and Innovation (IHPI).

Statistical analysis

We applied Bonferroni adjustment to identify significant variants and gene level associations from SAIGE and SAIGE-gene, respectively. Pearson correlation and Mann-Whitney tests were used to assess the associations between age/gender and rare disease diagnoses. We also used Pearson correlation and test of proportions to compare the prevalence of rare diseases in Clinformatics® and the UK Biobank. In each case, we present both the effect size and p-value, as appropriate.

Results

Rare diseases identified

Table 1 provides some examples of our disease mapping, along with the number of individuals in the UKB with the rare disease and group (identified from the ICD-10 chapter) the disease belongs to (Additional file 1: Table S2 contains the full list). Using our consensus approach, we mapped 1,176 ICD-10 codes to 720 ORPHA codes. 23,575 individuals in the UKB (~5% prevalence) were found to have at least one of 420 specific rare diseases, with 2,602 individuals having more than one rare disease (~10% prevalence). It is interesting that individuals with a rare disease have higher susceptibility to other rare diseases, as this suggests some of the rare conditions may be related, either causally or by shared molecular mechanisms. This finding confirms a previous study in the USA [4], which revealed that ~13% of individuals with a rare condition have more than one.

Table 1 Sample sizes obtained for some of the rare disorders. For each of the rare diseases we identified in the UK Biobank, we provide the set of ICD-10 codes that map directly to the codes from Orphanet, such that the ORPHA code is no more specific that the ICD-10 codes. We then provide the number of individuals with that rare disease in the UK Biobank (UKB count) along with the group the disease belongs to, based on its ICD-10 chapter

Figure 1a presents the distribution of the number of individuals who are recorded as having each disease, with the highest disease density occurring for 5.4 individuals (1 in 100,000 estimated prevalence). A similarly high density is observed for diseases with up to 30 individuals, and density decreases as the number of individuals increases. The vast majority of diseases (indicated as rare in Orphanet) also meet the criteria for being rare in the UKB, suggesting our mapping approach has successfully identified the rare diseases. However, 24 out of the 420 diseases (6%) have an estimated prevalence greater than 1 in 2000 (indicated by the dashed red line), meaning they are not rare by the European definition. Furthermore, 21 of these diseases do not meet the US criterion (affecting fewer than 200,000 people) when extrapolating to the population of the USA. This is important to consider with respect to the demographics of rare diseases in the UKB.

Fig. 1
figure 1

Rare diseases in the UK Biobank. a Density plot of the 420 rare diseases we identified in the UK Biobank by mapping ICD-10 codes to ORPHA codes. The x-axis shows the log10 number of individuals recorded as having each disease, while the y-axis shows the density of diseases with that number of individuals. The red dashed line indicates the (fewer than 1 in 2000) prevalence criterion for rare diseases in Europe. b Scatter plot comparing the prevalence of rare diseases in the UK Biobank and the Optum dataset. Each point represents a rare disease, and the dotted red line represents the linear regression between the prevalence in Optum and the UK Biobank. c The groups of rare diseases identified in the UK Biobank are shown as a bar plot, with the y-axis indicating the number of diseases in each group. Overlaid is a second bar plot (in red), with the y-axis indicating the number of individuals who have at least one disease in each group. d Hexagon/scatter plot showing the mean age at recruitment and proportion of males for each disease. The x-axis shows the percentage of males recorded as having each disease and the y-axis shows the mean age at recruitment. The hexagons show the density of diseases at each mean age and sex proportion, while the asterisks indicate the actual values for a particular disease. The red dashed line shows the overall mean age and sex proportion in the UK Biobank. e Box/scatter plot of comorbidities for rare diseases. The y-axis shows, for each rare disease, the percentage of individuals who have at least one comorbidity in each group (excluding the rare disease itself)

The total number of rare disease individuals in the UKB differs depending on whether diseases with >1/2,000 prevalence are included (23,575 individuals, 5% of the total) or excluded (10,635 individuals, 2% of the total). Both estimates are within the range of rare disease prevalence previously reported [4, 10], and for the remainder of our analysis, we include all 420 diseases, using Orphanet as the deciding factor over whether a disease should be considered rare. It is also worth noting our approach is conservative, since it only includes rare diseases that can be reliably mapped from ICD-10 to ORPHA. As many rare diseases are not precisely defined within ICD-10, so are excluded by our mapping approach, we should expect the overall prevalence of rare diseases to be higher than our reported number. However, it is also possible due to inaccuracies in diagnosis that some individuals annotated with a particular ICD-10 code might not have the disease it refers to.

Confirmation of applicability to Optum’s deidentified Clinformatics® Data Mart

By converting the ICD-10 codes from our consensus mapping to non-approximate ICD-9 codes, we created ICD-9 to ORPHA mappings for 114 rare diseases (Additional file 2: Table S3). We applied our ICD-9/10 mappings to Clinformatics® and compared the prevalence of individuals with each rare disease to the UKB (Fig. 1b). The prevalence of rare diseases in Clinformatics® and the UKB is highly correlated (r = 0.77, p = 2.2 × 10−23). Interestingly, celiac artery compression syndrome had substantially higher prevalence in Clinformatics® (0.044%) compared to the UKB (0.001%), as did burning mouth syndrome (Clinformatics®: 0.117%, UKB: 0.008%). However, generally rare diseases had higher prevalence in Clinformatics®, with only 11 having higher prevalence in the UKB. We further compared the prevalence estimated from the UKB against those reported in Orphanet and found that of the 186 rare diseases from our mapping that have prevalence information in Orphanet, 165 (89%) are within one order of magnitude difference; 66 (35%) have higher prevalence in the UKB, and 48 (26%) have lower prevalence in the UKB, compared with the range reported by Orphanet.

The relatively high recruitment age in the UKB might be expected to result in more rare diseases arising or being identified. Polymyalgia rheumatica is more common among older adults and does have a slightly higher prevalence in the UKB (0.39%) than Clinformatics® (0.38%). However, trigeminal neuralgia is also more common in older adults and has more than twice the prevalence in Clinformatics® (0.27%) than the UKB (0.12%), fold change (FC = 2.29, p = 9.03 × 10−208, test of proportions). Both Clinformatics® and the UKB may be subject to selection biases. In addition to variation between the UK and USA, there are also data collection differences: Clinformatics® claims are enriched for individuals in contact with a health system, while the UKB is community focused but may be biased towards healthy volunteers [32].

Quantifying groups of rare diseases

Figure 1c presents the number of rare diseases identified in the UKB for each non-overlapping group defined from ICD-10 chapters (congenital, neoplasms and infectious/parasitic rare diseases are the most numerous); although individuals can have diseases from multiple groups, each disease belongs to exactly one group. Overlaid (in red) are the number of individuals with at least one disease from each group (neoplasms, musculoskeletal and digestive system diseases are the most common). Orphanet also provides an overlapping set of categories (Additional file 2: Figure S1); each disease is associated with multiple categories, with genetic, neurological, and hematological categories among the most frequent.

Demographics in the UKB

From the UKB, we extracted the age at recruitment, gender, and comorbidities for individuals who have each disease. Individuals with a rare disease had slightly higher median age (60) than those without (58), Mann-Whitney p = 9.0 × 10−314 FC = 1.03; individuals with more than one rare disease also had a higher median age (61) than those with only one rare disease (60), Mann-Whitney p = 6.5 × 10−6 FC = 1.02 (some individuals with an ICD-10 code might not have the disease due to inaccuracies in diagnosis). Figure 1d and Additional file 2: Figure S2 show a modest (positive) correlation (r = 0.32, p = 4.8 × 10−8) between male proportion and median age at recruitment for different rare diseases, with male dominated diseases having higher age at recruitment. Some male-biased diseases have a low median recruitment age, including cocaine intoxication and Kaposi’s sarcoma (associated with HIV [33]). There are also female-biased diseases with older age, such as malignant tumor of fallopian tubes and vulvar intraepithelial neoplasia.

Although rare diseases affected a similar proportion of males and females overall, approximating the general population, there are considerable differences in sex proportion between rare disease groups (Fig. 2 and Additional file 2: Figure S3). Neoplasms and digestive system rare diseases are more frequent among males, while musculoskeletal and skin/subcutaneous rare diseases are more frequent among females. Figure 1e and Additional file 2: Figure S4 present the comorbidities of rare diseases, determined by ICD-10 codes. The percentage of individuals with comorbidities in each group differs substantially between rare diseases, illustrating diversity in how rare diseases relate to common comorbidities and suggesting the broad impact of rare diseases. Additional file 2: Figure S5 shows the enrichment of comorbidities among individuals with different groups of rare diseases as a heatmap. Most groups were significantly enriched for comorbidities in the same group (i.e., on the diagonal), while individuals with rare blood or genitourinary disorders were more likely to have pregnancy/childbirth comorbidities, and individuals with rare respiratory system disorders were more likely to have infectious/parasitic comorbidities.

Fig. 2
figure 2

Sex of individuals with different groups of rare disease. Each bar plot presents the number of male and female individuals in the UK Biobank who have at least one rare disease from a particular group. Non-overlapping groups of rare diseases were identified from their corresponding ICD-10 chapters. For each group, we conducted a Fisher enrichment test, comparing the number of males and females in the group with the number of males and females in the UK Biobank overall; p-values and odds ratios are provided under each bar plot and the subtitles of groups with significant sex bias are indicated in bold font

We compared our findings with those from a recent work [22] that studied comorbidities of common diseases in the UK Biobank (Additional file 2: Figure S6); having a comorbidity in the same category was generally less common for rare diseases than for common diseases, with the notable exception of rare respiratory system diseases, which have higher same-category and cross-category enrichment. As with our study (Additional file 2: Figure S5), the Dong et al. study reported significant comorbidities between common infectious and respiratory disorders, as well as between common respiratory and metabolic conditions. However, uniquely, we observed individuals with a rare form of neoplasm are more likely to have a blood disorder, whereas individuals who have rare blood disorders are not necessarily more likely to have neoplasms. Overall, the differences in enrichments across the disease categories between our study and the common disease study are minor, further supporting the accuracy of our analysis. We further compared our results with a study that investigated the association between Mendelian and complex genetic diseases [23], focusing on the 15 rare diseases from our mapping in the UK Biobank that it included. Interestingly, we identified 10 significant enrichments (Additional file 2: Figure S7) after Bonferroni correction, all of which were supported by the previous study. Although we did not observe any significantly enriched comorbidities between the Mendelian diseases, we expanded our analysis to all the rare diseases from our mapping and identified 222 significantly co-occurring rare diseases (Additional file 1: Table S4) and 389 significant complex disease comorbidities of rare diseases (Additional file 1: Table S5) after Bonferroni correction.

Interactive web browser for rare diseases in the UKB

To help facilitate the use of our findings for future rare disease research, we created a website (http://www-personal.umich.edu/~mattpat/rareDiseases.html) which allows rare diseases within the UKB to be easily explored. Our website includes an interactive table with details of each rare disease (their name, ORPHA code, ICD-10 codes, number of individuals in the UKB, and disease group). When the user accesses a disease in the table, more information is provided in the tabs above, including the prevalence, age at onset, age at death, inheritance type, genes, and categories from Orphanet (where available), in addition to a table of phenotypes and their expected frequency. On different tabs are figures showing the age at recruitment and sex of individuals, relative to the general population, and the proportion of individuals with different comorbidities, with a link to more information regarding the disease from Orphanet. To ensure the anonymity of participants in the UKB, no age, sex, or comorbidity information is shown for diseases affecting fewer than five individuals.

Exome sequencing analysis

Applying SAIGE/SAIGE-GENE to 33,981 variants across 162 diseases indicated as rare by Orphanet (< 1 in 2000 people) and using Bonferroni adjustment, we identified 19 significant variant-level associations and 20 significant gene-level associations. We further restricted to diseases indicated by NIH’s GARD as rare (< 200,000 people), resulting in 14 variant-level associations (Table 2) and 14 gene-level associations (Additional file 2: Table S6) significant for rare diseases. The full set of gene-level associations and variant level associations with false discovery rate (FDR) ≤ 0.05 are provided in Additional file 1: Table S7 and Additional file 1: Table S8, respectively. Six of the 14 variant-level associations (43%) had previously been reported in ClinVar.

Table 2 Significant variant-level associations. We applied SAIGE’s GLMM test with Bonferroni adjustment to identify significant variant associations for each rare disease and Bonferroni adjustment is applied. Genomic positions are provided in hg38 build and HGVS nomenclature. The ACMG/AMP classification for each association was determined through the use of Varsome and InterVar. Associations are specified as being reported if they have previously been indicated as pathogenic or likely pathogenic in ClinVar for that disease. Abbreviations are as follows: MAC, minor allele count; MAF, minor allele frequency. Although additional significant variant-level associations were identified for systemic lupus erythematosus and von Willebrand disease, these diseases were indicated as not being rare in the USA by NIH’s GARD. Furthermore, significant variant-level associations with interatrial communication, benign epithelial tumor salivary glands, and endophthalmitis were excluded because these diseases were not listed on NIH’s GARD, so it is difficult to confirm their rareness in the USA

The most significant associations involved the JAK2 V617F variant, previously indicated as pathogenic for various diseases in ClinVar. We confirmed associations for polycythemia vera (p = 1.32 × 10−114), chronic myeloproliferative disease (p = 2.89 × 10−67), primary myelofibrosis (p = 5.30 × 10−40), and essential thrombocythemia (p = 1.59 × 10−12). Furthermore, we identified a significant V617F association for immune thrombocytopenic purpura (p = 1.24 × 10−13), which while not indicated in ClinVar or Orphanet, was previously reported in case studies from Poland [34] and Italy [35] and in mouse experiments [36]. Outside of JAK2, significant associations were identified involving three variants indicated as pathogenic or likely pathogenic by ClinVar: HBB with beta-thalassemia (p = 7.34 × 10−12), F11 with congenital factor XI deficiency (p = 3.41 × 10−11), and MYD88 with B-cell chronic lymphocytic leukemia (p = 1.08 × 10−9). The HBB gene encodes for β-globin, an important subunit of hemoglobin [37], while F11 encodes for factor XI, which is needed for coagulation [38], and mutations in MYD88 can induce oncogenesis through their impact on NFκB and JAK signaling regulation [39].

We also identified a significant association between SRSF2 and chronic myelomonocytic leukemia (p = 1.19 × 10−13). Although SRSF2 is a known gene for this disease (potentially causing oncogenesis through its impact on CD4/CD8 T-cells [40]) and was indicated by Orphanet, our variant (rs751713049) is novel. It is indicated as pathogenic (Additional file 1: Table S8) by four predictors (CADD, PROVEAN, MutationTester and SIFT), as well as possibly pathogenic by MutationAssessor and PolyPhen2 (using both HVAR and HDIV). A novel association was identified between MYH1 and osteochondritis dissecans (p = 1.01 × 10−7), involving a variant predicted as pathogenic by seven predictors; MYH1 is involved in skeletal muscle and has been linked to rhabdomyolysis [41].

We performed association analysis on the subset of variants predicted to cause loss of function in SnpEff (Additional file 2: Table S9). Nine significant associations were identified after Bonferroni adjustment, including a frameshift variant in CALR for essential thrombocythemia (p = 1.59 × 10−13) and the previously mentioned stop gained variant in HBB for beta-thalassemia. The HBB variant is recorded in ClinVar as pathogenic for beta-thalassemia, and the gene is annotated in Orphanet for this disease. While CALR is annotated in Orphanet as involved in essential thrombocythemia, our variant is novel.

Discussion

Previous research has typically focused on specific rare conditions. Sun et al. [42] used ICD-9 codes to identify membranous nephropathy patients in the Kaiser Permanente health system, while Dickey et al. [43] used exome sequencing from the UKB to investigate whether erythropoietic protoporphyria may be under-diagnosed. The UKB [44] has also been analyzed more broadly, with rare variants in JAK2 and F11 associated to groups of myeloproliferative disease and congenital coagulation defects, respectively. Another study investigated extreme red blood cell indices in the UKB [45]. More recently, genetic association was performed for 33 rare diseases in data from 23andMe [46], and replication analysis confirmed significant associations for two diseases in the UKB.

Accurate disease identification is essential for expanding population-level biobank research to rare diseases. Although Orphanet indicated ICD-10 codes for > 6000 rare diseases, our consensus mapping approach revealed that only 566 of these were specific. ICD-10 codes are typically broader than their corresponding ORPHA codes and attempts to identify individuals using these codes would not be accurate for 92% of rare diseases. A previous study created a mapping from (the Australian modification of) ICD-10 codes to ORPHA codes [1] to investigate demographics in a public health system; however, it includes fewer rare diseases than our own; of the 14 significant variant-level associations (Table 2), only six involve a disease included in the Western Australia mapping. Other ontologies that have mappings for rare diseases include the UMLS [47], MonDO [10], and OMIM [48], and it is also possible to identify rare diseases by phenotype similarity [49]. We evaluated MonDO in our project by selecting the MonDO codes highest up the hierarchy (most general) for each ICD-10 code, then mapping to the exact match ORPHA code, where available (aiming for the most specific ORPHA code that is no more specific than the ICD-10 code). Of the 4682 ICD-10 codes mapped, 650 (14%) mapped to more than one ORPHA code, and some of the mappings did not meet our criteria. For example, A07.3 (isosporiasis) mapped to both ORPHA:472 (isosporiasis) and ORPHA:210 (cyclosporosis), which are caused by different parasites. A07.8 mapped to ORPHA:54368 (sarcocystosis), but its definition also includes other protozoal intestinal diseases. We therefore determined our consensus mapping approach was needed to identify specific rare diseases.

A single ICD-10 code can describe multiple diseases, so it is not always clear which disease an individual has, nor whether that disease is rare. Twenty-four of the diseases indicated as rare in Orphanet had greater than 1 in 2000 prevalence in the UKB; however, the UKB only recruited individuals aged 40–69, and this could lead to increased representation of rare diseases affecting that age range. To help ensure the diseases we identify are rare in Europe and the USA, we restricted our results to diseases reported by the NIH Genetic and Rare Diseases Information Center (GARD). Nevertheless, GARD leaves out some low prevalence diseases, such as cerebral sinovenous thrombosis (CSVT), while polymyalgia rheumatica is indicated as rare in both Orphanet and GARD but has a relatively high prevalence (~0.4%) in the UKB and Clinformatics®. Due to the challenges of diagnosing rare diseases, it is likely some individuals have a rare disease yet to be diagnosed.

Interestingly, many of the significant genetic associations we identified were for hematological diseases. While the proportion of hematological diseases we identified is large (Additional file 2: Figure S1), it is not the most frequent category; more neurological diseases were identified, yet we only found a significant genetic association for one of them (AA amyloidosis). Since the majority of UKB samples used for DNA extraction were taken from blood, some of the genetic associations we identified may be from de novo rather than hereditary variants. Identification of de novo variants is limited by the number of related individuals with rare diseases in the UKB. Of 24,416 individuals with a rare disease, 852 are estimated to be related to another individual in the UKB through genetics. Restricting to individuals with WES data, there are 22 pairs of related individuals where one has a rare disease. Of these, one pair includes an individual with Immune thrombocytopenic purpura and another includes an individual with B-cell chronic lymphocytic leukemia, however all individuals have the major allele for the associated variants.

Six of the significant (after Bonferroni correction) variant-level associations involve the JAK2 V617F variant, for rare diseases that can be broadly categorized as myeloproliferative disorders. Among the 484 individuals who have at least one of these rare diseases, 72 of them (15%) have more than one of the six diseases, and it is plausible the mechanisms driving these associations may overlap [50]. Individuals with this variant had a slightly higher median age at recruitment (63 vs 58, Mann-Whitney p = 2.32 × 10−12), suggesting that the associated marker is a somatic mutation. Expanding our analysis to include variants with FDR< = 0.05 in genetic association (Additional file 1: Table S8), we found 40 significantly comorbid rare disease pairs that share an associated variant (Additional file 2: Table S10). In addition to the JAK2 V617F variant pairs, immune thrombocytopenia, essential thrombocythemia, and chronic myeloproliferative disease were connected by SELENON G315S, which is indicated as pathogenic by ClinVar and has been associated with congenital myopathies among other diseases [51, 52].

Conclusions

We have shown how consensus mapping between ICD and ORPHA codes can reveal relevant demographics and genetic associations for a wide range of rare diseases in different population-level datasets. By analyzing exome sequencing data from 167,246 individuals of European ancestry in the UKB, we confirmed and identified pathogenic variants for rare diseases, and as sample sizes continue to increase, the power available can lead to important discoveries. We have provided our findings, along with the ICD-10/ORPHA mapping in an interactive website to facilitate investigation of specific rare diseases.

Disease coding systems are constantly evolving, and ICD-11 now includes over 5000 rare diseases [53, 54]. When ICD-11 is widely implemented across clinics, future research will extend and repeat our consensus mapping approach, so that it can continue to enable rare disease research in population-level datasets. Our approach can therefore increase the sample size for rare conditions, especially those currently under-represented. The diverse range of information provided by the UKB (including health records, genetics, drug prescriptions, lifestyle, family/medical history etc.) can improve understanding of the overall burden of rare diseases.

Availability of data and materials

Phenotypic data is provided through our interactive website http://www-personal.umich.edu/~mattpat/rareDiseases.html), and all information related to our mappings and significant genetic associations (adjusting for false discovery rate) are provided in supplementary material. Exome sequencing data of the UK Biobank is available to researchers on successful application and requires an access fee [55] (we downloaded the latest version available on 24 August 2020). The medical claims data used to demonstrate the transferability of our approach is available from the Clinformatics® Datamart available from Optum: optum.org/life-sciences-solutions (we used version 7.2). We obtained this data through the University of Michigan’s Institute for Healthcare Policy & Innovation (IHPI). The Clinformatics® Datamart includes a wide variety of data, including demographics, hospital admissions, procedures, and provider details; however, for the purposes of our study, we focused on the deidentified ICD-9/10 codes. Orphanet’s current mapping of ORPHA codes to ICD-10 codes is available online [56] (we downloaded the latest version available on 2 September 2020). The code from our analyses is available on GitHub [57].

Abbreviations

GARD:

NIH’s Genetic and Rare Diseases Information Center

ICD-10:

10th revision of the International Statistical Classification of Diseases and Related Health Problems

MAF:

Minor allele frequency

UKB:

UK Biobank

References

  1. Walker CE, Mahede T, Davis G, Miller LJ, Girschik J, Brameld K, et al. The collective impact of rare diseases in Western Australia: an estimate using a population-based cohort. Genet Med. 2017;19:546–52. https://doi.org/10.1038/gim.2016.143.

    Article  PubMed  Google Scholar 

  2. Soneda A, Teruya H, Furuya N, Yoshihashi H, Enomoto K, Ishikawa A, et al. Proportion of malformations and genetic disorders among cases encountered at a high-care unit in a children’s hospital. Eur J Pediatr. 2012;171:301–5. https://doi.org/10.1007/s00431-011-1534-2.

    Article  PubMed  Google Scholar 

  3. Schieppati A, Henter JI, Daina E, Aperia A. Why rare diseases are an important medical and social issue. Lancet. 2008;371:2039–41. https://doi.org/10.1016/s0140-6736(08)60872-7.

    Article  PubMed  Google Scholar 

  4. Bogart KR, Irvin VL. Health-related quality of life among adults with diverse rare disorders. Orphanet J Rare Dis. 2017;12:177. https://doi.org/10.1186/s13023-017-0730-1.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Ferreira CR. The burden of rare diseases. Am J Med Genet A. 2019;179:885–92. https://doi.org/10.1002/ajmg.a.61124.

    Article  PubMed  Google Scholar 

  6. Guillem P, Cans C, Robert-Gnansia E, Aymé S, Jouk PS. Rare diseases in disabled children: an epidemiological survey. Arch Dis Child. 2008;93:115–8. https://doi.org/10.1136/adc.2006.104455.

    CAS  Article  PubMed  Google Scholar 

  7. Wästfelt M, Fadeel B, Henter JI. A journey of hope: lessons learned from studies on rare diseases and orphan drugs. J Intern Med. 2006;260:1–10. https://doi.org/10.1111/j.1365-2796.2006.01666.x.

    Article  PubMed  Google Scholar 

  8. Luzzatto L, Hyry HI, Schieppati A, Costa E, Simoens S, Schaefer F, et al. Outrageous prices of orphan drugs: a call for collaboration. Lancet. 2018;392:791–4. https://doi.org/10.1016/s0140-6736(18)31069-9.

    Article  PubMed  Google Scholar 

  9. Detiček A, Locatelli I, Kos M. Patient access to medicines for rare diseases in European countries. Value Health. 2018;21:553–60. https://doi.org/10.1016/j.jval.2018.01.007.

    Article  PubMed  Google Scholar 

  10. Haendel M, Vasilevsky N, Unni D, Bologa C, Harris N, Rehm H, et al. How many rare diseases are there? Nat Rev Drug Discov. 2020;19:77–8. https://doi.org/10.1038/d41573-019-00180-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  11. Herder M. What is the purpose of the orphan drug act? PLoS Med. 2017;14:e1002191. https://doi.org/10.1371/journal.pmed.1002191.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Mitani AA, Haneuse S. Small data challenges of studying rare diseases. JAMA Netw Open. 2020;3:e201965. https://doi.org/10.1001/jamanetworkopen.2020.1965.

    Article  PubMed  Google Scholar 

  13. Bell SA, Tudur SC. A comparison of interventional clinical trials in rare versus non-rare diseases: an analysis of ClinicalTrials.gov. Orphanet J Rare Dis. 2014;9:170. https://doi.org/10.1186/s13023-014-0170-0.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Mora M, Angelini C, Bignami F, Bodin AM, Crimi M, Di Donato JH, et al. The EuroBioBank Network: 10 years of hands-on experience of collaborative, transnational biobanking for rare diseases. Eur J Hum Genet. 2015;23:1116–23. https://doi.org/10.1038/ejhg.2014.272.

    CAS  Article  PubMed  Google Scholar 

  15. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562:203–9. https://doi.org/10.1038/s41586-018-0579-z.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. Sankar PL, Parker LS. The Precision Medicine Initiative’s All of Us Research Program: an agenda for research on its ethical, legal, and social issues. Genet Med. 2017;19:743–50. https://doi.org/10.1038/gim.2016.183.

    Article  PubMed  Google Scholar 

  17. Nguengang Wakap S, Lambert DM, Olry A, Rodwell C, Gueydan C, Lanneau V, et al. Estimating cumulative point prevalence of rare diseases: analysis of the Orphanet database. Eur J Hum Genet. 2020;28:165–73. https://doi.org/10.1038/s41431-019-0508-0.

    Article  PubMed  Google Scholar 

  18. Armstrong AW, Mehta MD, Schupp CW, Gondo GC, Bell SJ, Griffiths CEM. Psoriasis prevalence in adults in the United States. JAMA Dermatol. 2021. https://doi.org/10.1001/jamadermatol.2021.2007.

  19. Gupta P, Zhao D, Guallar E, Ko F, Boland MV, Friedman DS. Prevalence of Glaucoma in the United States: The 2005-2008 National Health and Nutrition Examination Survey. Invest Ophthalmol Vis Sci. 2016;57:2905–13. https://doi.org/10.1167/iovs.15-18469.

    Article  PubMed  Google Scholar 

  20. Rath A, Olry A, Dhombres F, Brandt MM, Urbero B, Ayme S. Representation of rare diseases in health information systems: the Orphanet approach to serve a wide range of end users. Hum Mutat. 2012;33:803–8. https://doi.org/10.1002/humu.22078.

    Article  PubMed  Google Scholar 

  21. Victoire A, Magin P, Coughlan J, van Driel ML. Interventions for infantile seborrhoeic dermatitis (including cradle cap). Cochrane Database Syst Rev. 2019;3:Cd011380. https://doi.org/10.1002/14651858.CD011380.pub2.

    Article  PubMed  Google Scholar 

  22. Dong G, Feng J, Sun F, Chen J, Zhao XM. A global overview of genetically interpretable multimorbidities among common diseases in the UK Biobank. Genome Med. 2021;13:110. https://doi.org/10.1186/s13073-021-00927-6.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Blair DR, Lyttle CS, Mortensen JM, Bearden CF, Jensen AB, Khiabanian H, et al. A nondegenerate code of deleterious variants in Mendelian loci contributes to complex disease risk. Cell. 2013;155:70–80. https://doi.org/10.1016/j.cell.2013.08.030.

    CAS  Article  PubMed  Google Scholar 

  24. Zhou W, Nielsen JB, Fritsche LG, Dey R, Gabrielsen ME, Wolford BN, et al. Efficiently controlling for case-control imbalance and sample relatedness in large-scale genetic association studies. Nat Genet. 2018;50:1335–41. https://doi.org/10.1038/s41588-018-0184-y.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  25. Zhou W, Zhao Z, Nielsen JB, Fritsche LG, LeFaive J, Gagliano Taliun SA, et al. Scalable generalized linear mixed model for region-based association tests in large biobanks and cohorts. Nat Genet. 2020;52:634–9. https://doi.org/10.1038/s41588-020-0621-6.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  26. Cingolani P, Platts A, Wang le L, Coon M, Nguyen T, Wang L, et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly (Austin). 2012;6:80–92. https://doi.org/10.4161/fly.19695.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  27. Ioannidis NM, Rothstein JH, Pejaver V, Middha S, McDonnell SK, Baheti S, et al. REVEL: an ensemble method for predicting the pathogenicity of rare missense variants. Am J Hum Genet. 2016;99:877–85. https://doi.org/10.1016/j.ajhg.2016.08.016.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  28. Liu X, Li C, Mou C, Dong Y, Tu Y. dbNSFP v4: a comprehensive database of transcript-specific functional predictions and annotations for human nonsynonymous and splice-site SNVs. Genome Med. 2020;12:103. https://doi.org/10.1186/s13073-020-00803-9.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  29. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405–24. https://doi.org/10.1038/gim.2015.30.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Gunaseelan V, Kenney B, Lee JS, Hu HM. Databases for surgical health services research: Clinformatics Data Mart. Surgery. 2019;165:669–71. https://doi.org/10.1016/j.surg.2018.02.002.

    Article  PubMed  Google Scholar 

  31. CMS. 2018 ICD-10 CM and GEMs. https://www.cms.gov/Medicare/Coding/ICD10/2018-ICD-10-CM-and-GEMs.

  32. Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, et al. Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population. Am J Epidemiol. 2017;186:1026–34. https://doi.org/10.1093/aje/kwx246.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Krown SE, Lee JY, Dittmer DP. More on HIV-associated Kaposi’s sarcoma. N Engl J Med. 2008;358:535–6; author reply 6. https://doi.org/10.1056/NEJMc072994.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  34. Sobas MA, Wróbel T, Zduniak K, Podolak-Dawidziak M, Rybka J, Biedroń M, et al. Immune thrombocytopenia and JAK2V617F positive essential thrombocythemia: literature review and case report. Case Rep Hematol. 2017;2017:3725089. https://doi.org/10.1155/2017/3725089.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  35. Caocci G, Atzeni S, Usai M, La Nasa G. Essential thrombocytemia following immune thrombocytopenia with JAK2V617F mutation. Leuk Res Rep. 2018;9:14–5. https://doi.org/10.1016/j.lrr.2017.12.002.

    Article  PubMed  Google Scholar 

  36. Grisouard J, Hao-Shen H, Dirnhofer S, Wagner KU, Skoda RC. Selective deletion of Jak2 in adult mouse hematopoietic cells leads to lethal anemia and thrombocytopenia. Haematologica. 2014;99:e52–4. https://doi.org/10.3324/haematol.2013.100016.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. Aldakeel SA, Ghanem NZ, Al-Amodi AM, Osman AK, Al Asoom LI, Ahmed NR, et al. Identification of seven novel variants in the β-globin gene in transfusion-dependent and normal patients. Arch Med Sci. 2020;16:453–9. https://doi.org/10.5114/aoms.2019.84825.

    CAS  Article  PubMed  Google Scholar 

  38. Mohammed BM, Matafonov A, Ivanov I, Sun MF, Cheng Q, Dickeson SK, et al. An update on factor XI structure and function. Thromb Res. 2018;161:94–105. https://doi.org/10.1016/j.thromres.2017.10.008.

    CAS  Article  PubMed  Google Scholar 

  39. Ngo VN, Young RM, Schmitz R, Jhavar S, Xiao W, Lim KH, et al. Oncogenically active MYD88 mutations in human lymphoma. Nature. 2011;470:115–9. https://doi.org/10.1038/nature09671.

    CAS  Article  PubMed  Google Scholar 

  40. Aujla A, Linder K, Iragavarapu C, Karass M, Liu D. SRSF2 mutations in myelodysplasia/myeloproliferative neoplasms. Biomark Res. 2018;6:29. https://doi.org/10.1186/s40364-018-0142-y.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Alsaif HS, Alshehri A, Sulaiman RA, Al-Hindi H, Guzmán-Vega FJ, Arold ST, et al. MYH1 is a candidate gene for recurrent rhabdomyolysis in humans. Am J Med Genet A. 2021. https://doi.org/10.1002/ajmg.a.62188.

  42. Sun AZ, Shu YH, Harrison TN, Hever A, Jacobsen SJ, O’Shaughnessy MM, et al. Identifying patients with rare disease using electronic health record data: the Kaiser Permanente Southern California membranous nephropathy cohort. Perm J. 2020;24. https://doi.org/10.7812/tpp/19.126.

  43. Dickey AK, Quick C, Ducamp S, Zhu Z, Feng YA, Naik H, et al. Evidence in the UK Biobank for the underdiagnosis of erythropoietic protoporphyria. Genet Med. 2021;23:140–8. https://doi.org/10.1038/s41436-020-00951-8.

    CAS  Article  PubMed  Google Scholar 

  44. Zhao Z, Bi W, Zhou W, VandeHaar P, Fritsche LG, Lee S. UK Biobank whole-exome sequence binary phenome analysis with robust region-based rare-variant test. Am J Hum Genet. 2020;106:3–12. https://doi.org/10.1016/j.ajhg.2019.11.012.

    CAS  Article  PubMed  Google Scholar 

  45. Turro E, Astle WJ, Megy K, Gräf S, Greene D, Shamardina O, et al. Whole-genome sequencing of patients with rare diseases in a national health system. Nature. 2020;583:96–102. https://doi.org/10.1038/s41586-020-2434-2.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. Shringarpure SS, Wang W, Jiang Y, Acevedo A, Dhamija D, Cameron B, et al. Large-scale trans-ethnic replication and discovery of genetic associations for rare diseases with self-reported medical data. medRxiv. 2021:2021.06.09.21258643. https://doi.org/10.1101/2021.06.09.21258643.

  47. Bodenreider O. The Unified Medical Language System (UMLS): integrating biomedical terminology. Nucleic Acids Res. 2004;32:D267–70. https://doi.org/10.1093/nar/gkh061.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  48. Amberger JS, Bocchini CA, Scott AF, Hamosh A. OMIM.org: leveraging knowledge across phenotype-gene relationships. Nucleic Acids Res. 2019;47:D1038–d43. https://doi.org/10.1093/nar/gky1151.

    CAS  Article  PubMed  Google Scholar 

  49. Greene D, Richardson S, Turro E. Phenotype similarity regression for identifying the genetic determinants of rare diseases. Am J Hum Genet. 2016;98:490–9. https://doi.org/10.1016/j.ajhg.2016.01.008.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  50. Jones AV, Kreil S, Zoi K, Waghorn K, Curtis C, Zhang L, et al. Widespread occurrence of the JAK2 V617F mutation in chronic myeloproliferative disorders. Blood. 2005;106:2162–8. https://doi.org/10.1182/blood-2005-03-1320.

    CAS  Article  PubMed  Google Scholar 

  51. Bachmann C, Noreen F, Voermans NC, Schär PL, Vissing J, Fock JM, et al. Aberrant regulation of epigenetic modifiers contributes to the pathogenesis in patients with selenoprotein N-related myopathies. Hum Mutat. 2019;40:962–74. https://doi.org/10.1002/humu.23745.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  52. Ferreiro A, Quijano-Roy S, Pichereau C, Moghadaszadeh B, Goemans N, Bönnemann C, et al. Mutations of the selenoprotein N gene, which is implicated in rigid spine muscular dystrophy, cause the classical phenotype of multiminicore disease: reassessing the nosology of early-onset myopathies. Am J Hum Genet. 2002;71:739–49. https://doi.org/10.1086/342719.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Chute CG. The rendering of human phenotype and rare diseases in ICD-11. J Inherit Metab Dis. 2018;41:563–9. https://doi.org/10.1007/s10545-018-0172-5.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Aymé S, Bellet B, Rath A. Rare diseases in ICD11: making rare diseases visible in health information systems through appropriate coding. Orphanet J Rare Dis. 2015;10:35. https://doi.org/10.1186/s13023-015-0251-8.

    Article  PubMed  PubMed Central  Google Scholar 

  55. UKBiobank. Genetic data. Available from: http://www.ukbiobank.ac.uk/scientists-3/genetic-data/.

  56. Orphanet. http://www.orphadata.com/cgi-bin/index.php.

  57. Scripts for rare disease analysis. Available from: https://github.com/CutaneousBioinf/RareDiseaseAnalysis.

Download references

Acknowledgements

We acknowledge all the participants for the cohorts we included in the study.

Funding

This work was supported by the Rare Disease Foundation (M.T.P). L.C.T was supported by awards from the National Institutes of Health (K01AR072129, SBDRC 1P30AR075043), the Taubman Institute Innovation Project, the Dawn and Dudley Holmes Foundation, and the Babcock Memorial Trust. M.T.P was supported by a Precision Health Scholars Award from the University of Michigan and a Career Development Award from the Dermatology Foundation.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization: M.T.P., L.C.T.; data curation: M.T.P., R.B.; formal analysis: M.T.P., W.Z.; funding acquisition: M.T.P., L.C.T.; investigation: M.T.P., R.B. L.C.T.; methodology: M.T.P., W.Z., L.C.T.; software: M.T.P., W.Z.; supervision: L.C.T.; visualization: M.T.P.; writing—original draft: M.T.P., L.C.T.; writing—review and editing: M.T.P., R.B., W.Z., J.T.E., J.E.G., L.C.T. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Matthew T. Patrick or Lam C. Tsoi.

Ethics declarations

Ethics approval and consent to participate

This study conforms to the Declaration of Helsinki. The UK Biobank is advised by its Ethics Advisory Committee (EAC) and consent to take part in UK Biobank has been obtained from the participants. This research has been conducted using the UK Biobank Resource under Application Number 40212. All data used in this study was statistically de-identified by the appropriate covered entity. Access to data from the Clinformatics® Data Mart was obtained through a data use agreement according to a license agreement at the University of Michigan Institute for Healthcare Policy and Innovation (IHPI).

Consent for publication

Not applicable.

Competing interests

JEG has served as a consultant to AbbVie, Eli Lilly, Almirall, Celgene, BMS, Janssen, Prometheus, TimberPharma, Galderma, Novatis, MiRagen, and AnaptysBio and has received research support from AbbVie, SunPharma, Eli Lilly, Kyowa Kirin, Almirall, Celgene, BMS, Janssen, Prometheus, and TimberPharma. LCT has received support from Galderma and Janssen. The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1

: ICD-10 code to Orpha code consensus mapping. Table S2: Sample sizes and demographics obtained for rare disorders. Table S4: Significantly comorbid rare disease pairs (Bonferroni adjustment). Table S5: Complex diseases significantly comorbid with rare diseases (Bonferroni adjustment). Table S7: Significant gene-level associations (FDR<=0.05). Table S8: Significant variant-level associations (FDR<=0.05).

Additional file 2: Figure S1

: Rare disease mapping and frequency of rare disease categories. Figure S2: Hexagon/scatter plot showing the mean age at recruitment and proportion of males for the diseases in each category. Figure S3: Sex of individuals with different groups of rare disease, grouped by age. Figure S4: Box/scatter plot of comorbidities for rare diseases, grouped by age. Figure S5: Heatmap showing the enrichment of comorbidities for individuals with specific groups of rare diseases compared to the full set of individuals with any rare disease. The colors represent odds ratios (OR) from Fisher exact tests, while asterisks indicate enrichments with significant p-values (after Bonferroni correction). Figure S6: Heatmap comparing the enrichment of comorbidities for individuals with specific groups of rare diseases with those from a previous study on comorbidities for individuals with common diseases1. The colors represent differences in odds ratios (OR) from Fisher exact tests. Figure S7: Heatmap showing the enrichment of complex disease comorbidities for individuals with 15 specific rare diseases included in the list of Mendelian diseases from a previous paper2. The colors represent log10 odds ratios (OR) from Fisher exact tests, while asterisks indicate enrichments with significant p-values (after Bonferroni correction). Figure S8: Histogram showing the number of ICD-10 codes mapping to different numbers of ORPHA codes in the original Orphanet mapping (in grey) as well as the number of these codes for which we were able to identify a single ORPHA code, such that individuals with the ICD-10 code should be expected to have the rare disease indicated by the ORPHA code. Some ICD-10 codes originally mapped to a large number of ORPHA codes, but across each of the bins, we were able to identify an appropriate single ORPHA code for a large proportion of ICD-10 codes, through our consensus mapping approach. Table S3: Comparing prevalence of in the UK Biobank and Optum. Table S6: Significant gene-level associations (Bonferroni adjustment). Table S9: Significant associations with loss of function variants (Bonferroni adjustment). Table S10: Shared variants between significantly comorbid rare diseases. Supplementary Note: Improvement in mapping through our consensus approach.

Additional file 3.

Tutorial for using interactive website.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Patrick, M.T., Bardhi, R., Zhou, W. et al. Enhanced rare disease mapping for phenome-wide genetic association in the UK Biobank. Genome Med 14, 85 (2022). https://doi.org/10.1186/s13073-022-01094-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13073-022-01094-y

Keywords

  • Rare disease
  • UK Biobank
  • Genetic associations
  • Phenotyping
  • Demographics