Phenotype Risk Score Study of Rare Variants
Ontology highlight
ABSTRACT: This study was designed to assess the phenotypic effects of rare variants. Rare variants are difficult to study in a high-throughput manner because most cohorts are underpowered to detect associations. In order to gain power to test rare variants, we use Phenotype Risk Scores (PheRS) based on features of Mendelian diseases. PheRS is calculated using claims data from an EHR that is mapped to clinical features from OMIM's clinical descriptions of Mendelian disease. We calculated PheRS for 1,204 Mendelians diseases in a cohort of 21,701 individuals genotyped on the HumanExome BeadChip. We then tested for association between these diseases and rare variants in causal genes. A phenotype risk score is calculated as the weighted sum of features for a given disease. The features are defined as a set of phecodes (consolidated ICD codes) associated with a particular Mendelian disease as described in OMIM. The weights for each phecode are calculated as the log inverse prevalence of the phecode in our cohort. For an individual, their score equals the sum of the weights for each phecode that is present in their medical record. Our study was based on a cohort of individuals with genotype data linked to de-identified electronic health records (EHR) from Vanderbilt's BioVU resource. We hope that this method will help generate phenotype-genotype correlation data on rare variants which in turn may inform rare variant interpretations use for diagnostics.
SECONDARY ACCESSION(S): PRJNA428408PRJNA428409
REPOSITORIES: dbGaP
Dataset's files
| Action | DRS | |||
|---|---|---|---|---|
| README.txt | Txt | |||
| phs001516.eMERGE_EMRPhenome.analysis-PI.MULTI.tar.gz | Other | |||
| phs001516.pha004655.txt | Txt | |||
| GapExchange_phs001516.v1.p1.xml | Xml | |||
| dbGaPEx2.1.5.xsd | Other |
Items per page: 5 1 - 5 of 19 |
Similar Datasets
Project description:Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well-established that both common and rare sequence variants contribute to the formation of CL/P, however, the contribution of copy number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed, however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our case cohort compared to controls, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR/Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.
2022-12-01 | GSE212166 | GEO
Project description:Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well-established that both common and rare sequence variants contribute to the formation of CL/P, however, the contribution of copy number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed, however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our case cohort compared to controls, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR/Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.
2022-12-01 | GSE212165 | GEO
Project description:<p>The OPPERA project has two sub-studies investigating effects of genetic, physiological, psychological and clinical factors on the risk of developing painful temporomandibular disorder (TMD).</p> <ol> <li>OPPERA I baseline case-control study of chronic TMJD (dbGaP study accession <a href="./study.cgi?study_id=phs000762">phs000762</a>). This sub-study uses a case-control study design to investigate characteristics associated with chronic TMD. The parent OPPERA phenotype file contains an indicator variable (<a href="javascript:getPage(this, 'variable.cgi', 222740);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222740);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222740);return true;">PHASE1CASECONTROL</a>) identifying participants in this sub-study. Data users should use the relevant case-classification variable (<a href="javascript:getPage(this, 'variable.cgi', 222742);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222742);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222742);return true;">CHRNONICTMDCASE</a>) and apply statistical methods appropriate for case-control studies when analyzing TMJD in this sub-study.</li> <li>OPPERA I prospective cohort for first onset of TMJD (dbGaP study accession <a href="./study.cgi?study_id=phs000850">phs000850</a>). This sub-study uses a prospective cohort study design to investigate characteristics that predict risk of developing first-onset TMD. The parent OPPERA phenotype file contains an indicator variable (<a href="javascript:getPage(this, 'variable.cgi', 222743);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222743);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222743);return true;">PROSPECTIVECOHORT</a>) identifying participants in this sub-study. Data users should use the relevant case-classification variable (<a href="javascript:getPage(this, 'variable.cgi', 222746);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222746);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222746);return true;">INCIDENTCASE</a>) and the relevant time-in-study variable (<a href="javascript:getPage(this, 'variable.cgi', 222745);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222745);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222745);return true;">FOLLOWUPYEARS</a>) when analyzing TMJD in this sub-study. It is recommended that data users apply statistical methods appropriate for cohort studies of a dynamic population (for example, Cox regression models, where <a href="javascript:getPage(this, 'variable.cgi', 222746);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222746);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222746);return true;">INCIDENTCASE</a> is the censoring variable and <a href="javascript:getPage(this, 'variable.cgi', 222745);return true;" onclick="javascript:getPage(this, 'variable.cgi', 222745);return true;" oncontextmenu="javascript:getPage(this, 'variable.cgi', 222745);return true;">FOLLOWUPYEARS</a> is time in study). </li> </ol>
| phs000761 | dbGaP
Project description:Intellectual disability (ID) is a heterogeneous clinical entity and includes an excess of males who harbor variants on the X-chromosome (XLID). We report rare FAM50A missense variants in the original Armfield XLID syndrome family localized in Xq28 and four additional unrelated males with overlapping features. Our fam50a knockout (KO) zebrafish model exhibits abnormal neurogenesis and craniofacial patterning, and in vivo complementation assays indicate that the patient-derived variants are hypomorphic. RNA sequencing analysis from fam50a KO zebrafish show dysregulation of the transcriptome, with augmented spliceosome mRNAs and depletion of transcripts involved in neurodevelopment. Zebrafish RNA-seq datasets show a preponderance of 3’ alternative splicing events in fam50a KO, suggesting a role in the spliceosome C complex. These data are supported with transcriptomic signatures from cell lines derived from affected individuals and FAM50A protein-protein interaction data. In sum, Armfield XLID syndrome is a spliceosomopathy associated with aberrant mRNA processing during development.
2020-05-02 | GSE145711 | GEO
Project description:Intellectual disability (ID) is a heterogeneous clinical entity and includes an excess of males who harbor variants on the X-chromosome (XLID). We report rare FAM50A missense variants in the original Armfield XLID syndrome family localized in Xq28 and four additional unrelated males with overlapping features. Our fam50a knockout (KO) zebrafish model exhibits abnormal neurogenesis and craniofacial patterning, and in vivo complementation assays indicate that the patient-derived variants are hypomorphic. RNA sequencing analysis from fam50a KO zebrafish show dysregulation of the transcriptome, with augmented spliceosome mRNAs and depletion of transcripts involved in neurodevelopment. Zebrafish RNA-seq datasets show a preponderance of 3’ alternative splicing events in fam50a KO, suggesting a role in the spliceosome C complex. These data are supported with transcriptomic signatures from cell lines derived from affected individuals and FAM50A protein-protein interaction data. In sum, Armfield XLID syndrome is a spliceosomopathy associated with aberrant mRNA processing during development.
2020-05-02 | GSE145710 | GEO
Project description:Pathogenic variants in KMT5B, a lysine methyltransferase, are associated with global developmental delay, macrocephaly, autism, and congenital anomalies (OMIM# 617788). Given the novelty of this disorder, it has not been fully characterized. Deep phenotyping of the largest (n=43) patient cohort to date identified that hypotonia and congenital heart defects are novel prominent features. Both missense variants and putative loss-of-function variants resulted in slow growth in patient derived cell lines. KMT5B homozygous knockout mice were smaller in size than their wildtype littermates but did not have significantly smaller brains, suggesting relative macrocephaly, also noted as a prominent clinical feature. RNA-sequencing of patient lymphoblasts and Kmt5b-haploinsufficient mouse brains identified differentially expressed pathways associated with nervous system development and function including axon guidance signaling. Overall, we identified novel pathogenic variants and clinical features in KMT5B-related neurodevelopmental disorder and provide insights into the molecular mechanisms of the disorder using multiple model systems.
2023-01-01 | GSE184953 | GEO
Project description:Clinical whole genome sequencing has enabled the discovery of potentially pathogenic noncoding variants in the genomes of rare disease patients with a prior history of negative genetic testing. However, interpreting the functional consequences of noncoding variants and distinguishing those that contribute to disease etiology remains a challenge. Here we address this challenge by experimentally profiling the functional consequences of rare noncoding variants detected in a cohort of undiagnosed rare disease patients at scale using a massively parallel reporter assay. We demonstrate that this approach successfully identifies rare noncoding variants that alter the regulatory capacity of genomic sequences. In addition, we describe an integrative analysis that utilizes genomic features alongside patient clinical data to further prioritize candidate variants with an increased likelihood of pathogenicity. This works represents an important step towards establishing a framework for the clinical interpretation of noncoding variants.
2022-05-11 | GSE185795 | GEO
Project description:A mutation within FASTKD2 causes a rare form of Mendelian mitochondrial encephalomyopathy. To investigate whether and how RNA binding of FASTKD2 contributes to the disease phenotype, we identified the RNA targets of FASTKD2 by iCLIP.
2015-11-09 | E-MTAB-4010 | biostudies-arrayexpress
Project description:Although Mendelian disorders are overwhelmingly attributed to protein-coding pathogenic variants, a majority of unsolved cases do not harbor obvious causal pathogenic variants in the coding sequence, suggesting a potential non-coding etiology. However, classification of pathogenicity in non-coding sequence remains prohibitive due to a vastly increased search space and the lack of a standardized rubric for interpretation. Here, we present an integrated single cell multiomic framework to nominate pathogenic non-coding variants for the congenital cranial dysinnervation disorders (CCDDs). The CCDDs are Mendelian neurodevelopmental disorders that result from aberrant development of cranial motor neurons in the embryonic brainstem. We created a non-coding reference atlas of single cell chromatin accessibility profiles for 86,089 embryonic mouse cranial motor neurons (cMNs). We found that high-quality single cell ATAC-seq (scATAC) profiles alone were a strong predictor of enhancement (64% in vivo validation rate). To further aid in interpretation, we integrated single cell histone modification and gene expression information to distinguish individual enhancers and their cognate genes. Relatively subtle differences in cellular composition of input data often led to substantial differences in predicted enhancer strength, cognate gene, and tissue of activity. Next, we mapped candidate non-coding variants from 899 whole genome sequences from 270 CCDD pedigrees to the murine cMN-specific regulatory elements and trained a machine learning classifier to accurately predict the functional effects of patient variants within these elements. We then performed high coverage scATACseq and site-specific footprinting analysis on an allelic series of CRISPR-humanised mice to validate our machine learning predictions and render important clues to the mode of pathogenicity. Finally, we performed peak- and gene-centric allelic aggregation to nominate non-coding variants, including those regulating MN1 and EBF3, respectively. Altogether this work extends non-coding variant analysis to Mendelian disease and presents a generalizable framework for nominating novel non-coding variants in other rare disorders.
2024-08-06 | GSE254089 | GEO
Project description:Although Mendelian disorders are overwhelmingly attributed to protein-coding pathogenic variants, a majority of unsolved cases do not harbor obvious causal pathogenic variants in the coding sequence, suggesting a potential non-coding etiology. However, classification of pathogenicity in non-coding sequence remains prohibitive due to a vastly increased search space and the lack of a standardized rubric for interpretation. Here, we present an integrated single cell multiomic framework to nominate pathogenic non-coding variants for the congenital cranial dysinnervation disorders (CCDDs). The CCDDs are Mendelian neurodevelopmental disorders that result from aberrant development of cranial motor neurons in the embryonic brainstem. We created a non-coding reference atlas of single cell chromatin accessibility profiles for 86,089 embryonic mouse cranial motor neurons (cMNs). We found that high-quality single cell ATAC-seq (scATAC) profiles alone were a strong predictor of enhancement (64% in vivo validation rate). To further aid in interpretation, we integrated single cell histone modification and gene expression information to distinguish individual enhancers and their cognate genes. Relatively subtle differences in cellular composition of input data often led to substantial differences in predicted enhancer strength, cognate gene, and tissue of activity. Next, we mapped candidate non-coding variants from 899 whole genome sequences from 270 CCDD pedigrees to the murine cMN-specific regulatory elements and trained a machine learning classifier to accurately predict the functional effects of patient variants within these elements. We then performed high coverage scATACseq and site-specific footprinting analysis on an allelic series of CRISPR-humanised mice to validate our machine learning predictions and render important clues to the mode of pathogenicity. Finally, we performed peak- and gene-centric allelic aggregation to nominate non-coding variants, including those regulating MN1 and EBF3, respectively. Altogether this work extends non-coding variant analysis to Mendelian disease and presents a generalizable framework for nominating novel non-coding variants in other rare disorders.
2024-08-06 | GSE254084 | GEO