Combined Annotation Dependent Depletion Score Research Articles

Are rare heterozygous SYNJ1 variants associated with Parkinson’s disease?

Previous studies have established that rare biallelic SYNJ1 mutations cause autosomal recessive parkinsonism and Parkinson’s disease (PD). We analyzed 8165 PD cases, 818 early-onset-PD (EOPD, < 50 years) and 70,363 controls. Burden meta-analysis revealed an association between rare nonsynonymous variants and variants with high Combined Annotation-Dependent Depletion score (> 20) in the Sac1 SYNJ1 domain and PD (Pfdr = 0.040). A meta-analysis of EOPD patients demonstrated an association between all rare heterozygous SYNJ1 variants and PD (Pfdr = 0.029).

Clinical relevance of protein-truncating variants of germline DNA repair genes in prostate cancer

BackgroundInterpreting genetic variants remains a challenge in prostate cancer (PCa). Although many annotation tools are available for prioritizing causal variants, the clinical relevance of these variants is rarely studied.MethodsWe collected a cohort study that included 274 PCa patients from June 2017 to December 2020 and sequenced 19 DNA damage repair (DDR) genes in these patients and explored the clinical consequence of these different approaches. We also examined all-cause and PCa-specific survival in DDR gene mutation carriers compared to non-carriers after androgen receptor (AR)-directed therapy.ResultsWe identified 13 variants from 19 DDR genes in a total of 14 (5.1%) patients who had at least one presumed pathogenic mutation using different annotation methods. Four variants were annotated as pathogenic, 11 variants were predicted as protein-truncating variants (PTVs), four variants received proxy-deleterious (Combined Annotation-Dependent Depletion scores of > 30), and only one variant was identified as a pathogenic variant or as having a functional effect by all three methods. PCa patients with PTVs were significantly associated with early onset, high cancer stage, and a worse response to AR-directed treatment. However, patients carrying a proxy-deleterious variant were only associated with a higher T (tumor) stage and N (node) stage than those without such a variant, but not associated with other clinical characteristics. In patients treated with AR-directed therapy, patients with a PTV showed an increased risk of all-cause death (adjusted hazard ratio (aHR) = 3.51, 95% confidence interval (CI): 1.06 ~ 11.56) and PCa-specific death (aHR = 4.49, 95% CI: 1.87 ~ 10.77) compared to non-PTV carriers after adjustment. We were unable to examine gene-specific risks due to the small number of patients.ConclusionsPTVs may assist in guiding treatment and early screening in PCa, while population-specific data for pathogenic variants are still being amassed.

Pediatric Epilepsy Genetic Testing Results and Long-term Seizure Freedom.

Objective: To determine whether there is a correlation of genetic diagnosis/result with long-term seizure freedom in pediatric epilepsy patients. Methods: This was a prospective and retrospective cohort study of children with epilepsy referred for genetic testing at a single center. The primary outcomes were presence and type of genetic diagnosis (pathogenic, benign, or variant of uncertain significance) and patient epilepsy status (seizure free, treatment failure, uncertain). Epilepsy gene panels were the primary method of genetic testing. Results: The prospective cohort had 22 patients followed for >11 years and for whom genetic testing was then performed; the retrospective cohort had 78 patients with previous genetic testing followed for >8 years. In the prospective cohort, one patient each of the seizure free or treatment failure groups had a pathogenic genetic variant; mean Combined Annotation Dependent Depletion (CADD) scores 22 and 24, respectively (P = .62). In the retrospective cohort, there was no difference in the number of variants (P = .97), the variant interpretations (P = .29 ClinVar, P = .39 lab interpretation) or mean CADD scores (P = .29) between the seizure-free, treatment failure, and uncertain epilepsy patients. Whole exome and genome sequencing identified pathogenic variants in 70% of patients with treatment failure but were not performed in seizure-free patients. Significance: Our findings show no correlation of the presence or type of epilepsy gene panel result with long-term seizure freedom in pediatric patients. The yield and specificity of pathogenic variants may be higher using whole exome and whole genome sequencing in patients with treatment-resistant epilepsy. Whole exome and whole genome sequencing, or more targeted understanding of specific variants, may be needed to improve the utility of pediatric epilepsy genetic testing.

Protein-truncating and rare missense variants in ATM and CHEK2 and associations with cancer in UK Biobank whole-exome sequence data

BackgroundDeleterious germline variants in ATM and CHEK2 have been associated with a moderately increased risk of breast cancer. Risks for other cancers remain unclear.MethodsCancer associations for coding variants in ATM...

Whole genome sequencing identifies associations for nonsyndromic sagittal craniosynostosis with the intergenic region of BMP2 and noncoding RNA gene LINC01428

Craniosynostosis (CS) is a major birth defect resulting from premature fusion of cranial sutures. Nonsyndromic CS occurs more frequently than syndromic CS, with sagittal nonsyndromic craniosynostosis (sNCS) presenting as the most common CS phenotype. Previous genome-wide association and targeted sequencing analyses of sNCS have identified multiple associated loci, with the strongest association on chromosome 20. Herein, we report the first whole-genome sequencing study of sNCS using 63 proband-parent trios. Sequencing data for these trios were analyzed using the transmission disequilibrium test (TDT) and rare variant TDT (rvTDT) to identify high-risk rare gene variants. Sequencing data were also examined for copy number variants (CNVs) and de novo variants. TDT analysis identified a highly significant locus at 20p12.3, localized to the intergenic region between BMP2 and the noncoding RNA gene LINC01428. Three variants (rs6054763, rs6054764, rs932517) were identified as potential causal variants due to their probability of being transcription factor binding sites, deleterious combined annotation dependent depletion scores, and high minor allele enrichment in probands. Morphometric analysis of cranial vault shape in an unaffected cohort validated the effect of these three single nucleotide variants (SNVs) on dolichocephaly. No genome-wide significant rare variants, de novo loci, or CNVs were identified. Future efforts to identify risk variants for sNCS should include sequencing of larger and more diverse population samples and increased omics analyses, such as RNA-seq and ATAC-seq.

Genomic Biomarkers Can Provide a Deeper Understanding of Recurrent Pressure Injuries.

To identify genetic biomarkers predisposing individuals with spinal cord injury (SCI) to recurrent pressure injuries (PIs). Repeated measures of the transcriptome profile of veterans with SCI at three Veterans Spinal Cord Injuries and Disorders Centers. Exclusion criteria included having significant active systemic disease at time of enrollment. Researchers obtained comprehensive profiles of clinical and health factors and demographic information relevant to PI history at enrollment and at each follow-up visit by reviewing patients' medical charts. Whole blood samples were collected at 6- to 12-month intervals for 2 to 4 years. In addition to DNA profiling with whole genome sequencing of the patients, RNA sequencing was performed to assess pathways associated with PI risk. Whole genome sequencing analysis identified 260 genes that showed increased prevalence of single-nucleotide variations in exonic regions with high (>20) combined annotation-dependent depletion scores between persons with high versus low intramuscular adipose tissue levels when cross-referenced with persons who had recurrent PIs. Gene set enrichment analysis using Hallmark and KEGG (Kyoto Encyclopedia of Genes and Genomes) gene sets of these candidate genes revealed enrichment in genes encoding proteins involved in fatty acid metabolism (P < .01). Further, RNA sequencing revealed upregulated activity in biological senescence pathways and downregulated activity in antimicrobial protection pathways. Genomic biomarkers may complement electronic health records to support management of complex interactive health issues such as risk of recurrent PIs in people with SCI. These findings may also be leveraged for homogeneous phenotypic grouping of higher-risk individuals.

IFIH1 deficiency causing neonatal herpes simplex virus encephalitis and recurrent retinitis

Homozygous loss-of-function variants in Interferon induced with helicase C domain 1 (IFIH1) leads to immunodeficiency with early onset recurrent and severe viral respiratory infections (Immunodeficiency 95; MIM 619773). Severe respiratory infections with RNA viruses such as rhinovirus and respiratory syncytial virus (RSV) are most characteristic of IFIH1 deficiency. IFIH1 encodes melanoma differentiation-associated protein 5 (MDA5), an intracellular sensor of mainly long double-stranded RNA (ds-RNA), a viral replication intermediate. When activated, MDA5 triggers antiviral pathways including increased transcription of Type 1-interferon. The antiviral role of MDA5 is best described for positive (+) RNA viruses such as Picornavirus family which produce long ds-RNA as a replication intermediate. We present a patient with neonatal herpes simplex virus (HSV) encephalitis and recurrent retinitis with homozygous variant of unknown significance (VUS) in IFIH1, c.2062A>G (p.Lys688Glu).Patient is a 10-year-old male, born of consanguineous parents from the United Arab Emirates. Past medical history is notable for neonatal HSV encephalitis, grade 3 intraventricular hemorrhage, cerebral palsy, epilepsy, optic atrophy and recurrent HSV retinitis. Patient has no additional history of severe or recurrent infections. No family history of recurrent or severe infections, malignancy, autoimmunity or autoinflammation.Immune evaluation demonstrated anemia with elevated mean corpuscular volume and mild elevation of platelets and monocytes. Patient had normal lymphocyte subsets, normal T and B cell subset distributions, protective vaccine titers and normal immunoglobins apart from mildly elevated immunoglobulin G. An Invitae primary immunodeficiency panel demonstrated homozygous VUS in IFIH1, c.2062A>G (p.Lys688Glu). This variant is not found in gnomAD or ClinVar and has a Combined Annotation Dependent Depletion score of 22.6. Functional testing demonstrated decreased interferon induction after stimulation with intracellular RNA, intracellular DNA and agonists of both Toll-like receptor 7 and 9. These findings are consistent with MDA5 deficiency given decreased interferon response with known MDA5 agonists.In summary, we describe a patient with recurrent and severe infections with a DNA virus and no known infections with respiratory RNA viruses. Both DNA and RNA viruses produce long ds-RNA as replication intermediates, an MDA5 agonist. This patient highlights the anti-viral role of MDA5 for both DNA and RNA viruses.

Rare TMEM230 Variants are Potentially Associated With Early Onset Parkinson’s Disease (P3-11.013)

Rare TMEM230 Variants are Potentially Associated With Early Onset Parkinson’s Disease (P3-11.013)

Genome-Wide Association and Inheritance-Based Analyses Implicate Unconventional Myosin Genes in Hypoplastic Left Heart Syndrome.

Deciphering hypoplastic left heart syndrome (HLHS) pathogenesis is confounded by its genetic heterogeneity and oligogenic underpinnings. Whole genome sequences were analyzed by 3 independent strategies to identify HLHS gene candidates, ranked by variant, gene, and disease-level metrics. First, a genome-wide association study of 174 cases and 853 controls revealed suggestive association with a MYO18B intron 33 variant (rs2269628-G; frequency=0.55 versus 0.39; OR, 1.97 [95% CI, 1.54-2.52]; P=6.70×10-8). Second, transmission disequilibrium testing of 161 HLHS proband-parent trios revealed overrepresentation of a MYO18B intron 42 variant (rs73154186-A; frequency=0.05; OR, 24 [95% CI, 3.2-177.4]; P=4.23×10-6). Third, rare, predicted-damaging variants were filtered in 2 multiplex families. In 141H, 2 fifth-degree relatives with HLHS shared a paternally-inherited MYO5A missense variant (p.Arg801Trp; frequency=0.00003; combined annotation-dependent depletion score=29), each with a maternally-inherited or de novo candidate modifier variant in a MYO5A-interacting conventional myosin. In 442H, a HLHS proband was compound heterozygous for MYO15A variants-a maternally-inherited pathogenic stop-gain variant co-segregating with tetralogy of Fallot and bicuspid aortic valve in maternal relatives (p.Tyr2819Ter; frequency=0.00003) and a paternally-inherited intronic variant altering a canonical transcription factor binding site (rs1277068603; frequency=0.00001; position weight matrix score=0.98). Collectively, these findings suggest that common and rare alleles within unconventional myosin genes are associated with HLHS susceptibility. The identified candidate MYO18B regulates cardiac sarcomerogenesis, supporting the hypothesis of intrinsic myogenic perturbation in arrested left heart development.

Targeted gene panel analysis of Japanese patients with maturity-onset diabetes of the young-like diabetes mellitus: Roles of inactivating variants in the ABCC8 and insulin resistance genes.

To investigate the genetic background of Japanese patients with suspected maturity-onset diabetes of the young (MODY). On 340 proband patients referred from across Japan, genomic variants were analyzed using a targeted multigene panel analysis combined with the multiplex ligation probe amplification (MLPA) analysis, mitochondrial m.3243A > G analysis and methylation-specific polymerase chain reaction of the imprinted 6q24 locus. Pathogenic/likely pathogenic variants were listed according to the 2015 American College of Medical Genetics and Genomics and the Association for Molecular Pathology criteria. Additionally, variants with a population frequency <0.001 and Combined Annotation Dependent Depletion score >20 (CS >20) were listed as rare variants of uncertain significance-CS >20. A total of 157 pathogenic/likely pathogenic variants and 44 rare variants of uncertain significance-CS >20 were identified. In the pathogenic/likely pathogenic variants, alterations in the GCK gene were the most common (82, 52.2%) followed by HNF1A (29,18.5%), HNF4A (13, 8.3%) and HNF1B (13, 8.3%). One patient was a 29.5% mosaic with a truncating INSR variant. In the rare variants of uncertain significance-CS >20, 20 (45.5%) were in the genes coding for the adenosine triphosphate-sensitive potassium channel, KCNJ11 or ABCC8, and four were in the genes of the insulin-signaling pathway, INSR and PIK3R1. Four variants in ABCC8 were previously reported in patients with congenital hyperinsulinism, suggesting the inactivating nature of these variants, and at least two of our patients had a history of congenital hyperinsulinism evolving into diabetes. In two patients with INSR or PIK3R1 variants, insulin resistance was evident at diagnosis. Causative genomic variants could be identified in at least 46.2% of clinically suspected MODY patients. ABCC8-MODY with inactivating variants could represent a distinct category of MODY. Genes of insulin resistance should be included in the sequencing panel for MODY.

Linking single nucleotide polymorphisms to signaling blueprints in abdominal aortic aneurysms

Abdominal aortic aneurysms (AAA) is a multifactorial complex disease with life-threatening consequences. While Genome-wide association studies (GWAS) have revealed several single nucleotide polymorphisms (SNPs) located in the genome of individuals with AAA, the link between SNPs with the associated pathological signals, the influence of risk factors on their distribution and their combined analysis is not fully understood. We integrated 86 AAA SNPs from GWAS and clinical cohorts from the literature to determine their phenotypical vulnerabilities and association with AAA risk factors. The SNPs were annotated using snpXplorer AnnotateMe tool to identify their chromosomal position, minor allele frequency, CADD (Combined Annotation Dependent Depletion), annotation-based pathogenicity score, variant consequence, and their associated gene. Gene enrichment analysis was performed using Gene Ontology and clustered using REVIGO. The plug-in GeneMANIA in Cytoscape was applied to identify network integration with associated genes and functions. 15 SNPs affecting 20 genes with a CADD score above ten were identified. AAA SNPs were predominantly located on chromosome 3 and 9. Stop-gained rs5516 SNP obtained high frequency in AAA and associated with proinflammatory and vascular remodeling phenotypes. SNPs presence positively correlated with hypertension, dyslipidemia and smoking history. GO showed that AAA SNPs and their associated genes could regulate lipid metabolism, extracellular matrix organization, smooth muscle cell proliferation, and oxidative stress, suggesting that part of these AAA traits could stem from genetic abnormalities. We show a library of inborn SNPs and associated genes that manifest in AAA. We uncover their pathological signaling trajectories that likely fuel AAA development.

Clinical features and pathogenicity assessment in patients with HTRA1-autosomal dominant disease.

Heterozygous mutations in HTRA1 were recently found to cause autosomal dominant cerebral small vessel disease (CSVD), and it was named HTRA1-autosomal dominant disease (AD-HTRA1) in the consensus recommendations of the European Academy of Neurology. This study aimed to investigate the clinical features of a mutation in HTRA1 and the effect of HTRA1 mutation on white matter hyperintensity (WMH). A proband's brain magnetic resonance imaging (MRI) showed multiple lacunar infarctions and multiple WMH in the lateral ventricle, external capsule, frontal lobe and corpus callosum. The proband and family members were tested for CSVD-related genes by next-generation sequencing and the clinical data of the patients were collected. The published literature on AD-HTRA1 was collected, and the clinical characteristics and pathogenicity of the patients were summarized. Combined Annotation Dependent Depletion (CADD) is a tool for scoring the deleteriousness of single-nucleotide variants and insertion/deletion variants in the human genome. The relationship between the degree of WMH and the pathogenicity of the mutation was further analyzed. It was found that the proband and her family members had a heterozygous missense mutation of c.854C > T (p.P285L) in the 4 exon of HTRA1 gene. A retrospective analysis of 5 families with c.854C > T mutation found that the patients had an early age of onset, cognitive impairment was more common, and alopecia and spondylosis could be combined at the same time. By univariate analysis, the severity of WMH was found to be significantly associated with the mutated CADD score (p < 0.05, Spearman's rho = 0.266). The clinical manifestations of AD-HTRA1 with mutation site c.854C > T (p.P285L) are similar to CARASIL, and brain MRI are mainly moderate or severe WMH and lacunar infarction (LI). WMH are affected by mutation sites. Therefore, our pathogenicity score for mutations can predict the severity of WMH.

Molecular Genetic Screening in Patients With ACE Inhibitor/Angiotensin Receptor Blocker-Induced Angioedema to Explore the Role of Hereditary Angioedema Genes.

Angioedema is a relatively rare but potentially life-threatening adverse reaction to angiotensin-converting enzyme inhibitors (ACEi) and angiotensin receptor blockers (ARBs). As with hereditary forms of angioedema (HAE), this adverse reaction is mediated by bradykinin. Research suggests that ACEi/ARB-induced angioedema has a multifactorial etiology. In addition, recent case reports suggest that some ACEi/ARB-induced angioedema patients may carry pathogenic HAE variants. The aim of the present study was to investigate the possible association between ACEi/ARB-induced angioedema and HAE genes via systematic molecular genetic screening in a large cohort of ACEi/ARB-induced angioedema cases. Targeted re-sequencing of five HAE-associated genes (SERPING1, F12, PLG, ANGPT1, and KNG1) was performed in 212 ACEi/ARB-induced angioedema patients recruited in Germany/Austria, Sweden, and Denmark, and in 352 controls from a German cohort. Among patients, none of the identified variants represented a known pathogenic variant for HAE. Moreover, no significant association with ACEi/ARB-induced angioedema was found for any of the identified common [minor allele frequency (MAF) >5%] or rare (MAF < 5%) variants. However, several non-significant trends suggestive of possible protective effects were observed. The lowest p-value for an individual variant was found in PLG (rs4252129, p.R523W, p = 0.057, p.adjust > 0.999, Fisher’s exact test). Variant p.R523W was found exclusively in controls and has previously been associated with decreased levels of plasminogen, a precursor of plasmin which is part of a pathway directly involved in bradykinin production. In addition, rare, potentially functional variants (MAF < 5%, Phred-scaled combined annotation dependent depletion score >10) showed a nominally significant enrichment in controls both: 1) across all five genes; and 2) in the F12 gene alone. However, these results did not withstand correction for multiple testing. In conclusion, our results suggest that HAE-associated mutations are, at best, a rare cause of ACEi/ARB-induced angioedema. Furthermore, we were unable to identify a significant association between ACEi/ARB-induced angioedema and other variants in the investigated genes. Further studies with larger sample sizes are warranted to draw more definite conclusions concerning variants with limited effect sizes, including protective variants.

The Phenotypic Continuum of ATP1A3-Related Disorders.

Background and ObjectivesATP1A3 is associated with a broad spectrum of predominantly neurologic disorders, which continues to expand beyond the initially defined phenotypes of alternating hemiplegia of childhood, rapid-onset dystonia parkinsonism, and cerebellar ataxia, areflexia, pes cavus, optic atrophy, sensorineural hearing loss syndrome. This phenotypic variability makes it challenging to assess the pathogenicity of an ATP1A3 variant found in an undiagnosed patient. We describe the phenotypic features of individuals carrying a pathogenic/likely pathogenic ATP1A3 variant and perform a literature review of all ATP1A3 variants published thus far in association with human neurologic disease. Our aim is to demonstrate the heterogeneous clinical spectrum of the gene and look for phenotypic overlap between patients that will streamline the diagnostic process.MethodsUndiagnosed individuals with ATP1A3 variants were identified within the cohort of the Deciphering Developmental Disorders study with additional cases contributed by collaborators internationally. Detailed clinical data were collected with consent through a questionnaire completed by the referring clinicians. PubMed was searched for publications containing the term “ATP1A3” from 2004 to 2021.ResultsTwenty-four individuals with a previously undiagnosed neurologic phenotype were found to carry 21 ATP1A3 variants. Eight variants have been previously published. Patients experienced on average 2–3 different types of paroxysmal events. Permanent neurologic features were common including microcephaly (7; 29%), ataxia (13; 54%), dystonia (10; 42%), and hypotonia (7; 29%). All patients had cognitive impairment. Neuropsychiatric diagnoses were reported in 16 (66.6%) individuals. Phenotypes were extremely varied, and most individuals did not fit clinical criteria for previously published phenotypes. On review of the literature, 1,108 individuals have been reported carrying 168 different ATP1A3 variants. The most common variants are associated with well-defined phenotypes, while more rare variants often result in very rare symptom correlations, such as are seen in our study. Combined Annotation-Dependent Depletion (CADD) scores of pathogenic and likely pathogenic variants were significantly higher and variants clustered within 6 regions of constraint.DiscussionOur study shows that looking for a combination of paroxysmal events, hyperkinesia, neuropsychiatric symptoms, and cognitive impairment and evaluating the CADD score and variant location can help identify an ATP1A3-related condition, rather than applying diagnostic criteria alone.

The Diverse Phenotype of Intestinal Dysmotility Secondary to ACTG2-related Disorders.

The initial description of a heterozygous dominant ACTG2 variant in familial visceral myopathy was followed by the identification of additional variants in other forms of intestinal dysmotility disorders. we aimed to describe the diverse phenotype of this newly reported and rare disease. Report of 4 new patients, and a systematic review of ACTG2-related disorders. we analyzed the population frequency and used in silico gene damaging predictions. Genotype-phenotype correlations were explored. One hundred three patients (52% girls), from 14 publications, were included. Twenty-eight unique variants were analyzed, all exceedingly rare, and 27 predicted to be highly damaging. The median Combined Annotation Dependent Depletion (CADD) score was 29.2 (Interquartile range 26.3-29.4). Most patients underwent abdominal surgery (66%), about half required intermittent bladder catheterization (48.5%), and more than half were parenteral nutrition (PN)-dependent (53%). One-quarter of the patients died (25.7%), and 6 required transplant (5.8%). Girls had a higher rate of microcolon (P = 0.009), PN dependency (P = 0.003), and death/transplant (P = 0.029) compared with boys, and early disease onset (<2 years of age) was associated with megacystis-microcolon-intestinal hypoperistalsis syndrome (MMIHS) features. There was no statistical association between disease characteristics and CADD scores. Damaging ACTG2 variants are rare, often associated with MMIHS phenotype, and overall have a wide phenotypic variation. Symptoms usually present in the perinatal period but can also appear at a later age. The course of the disease is marked by frequent need for surgical interventions, PN support, and mortality. Poor outcomes are more common among girls with ACTG2 variants.

Whole exome sequencing and analysis of hypohidrotic ectodermal dysplasia patients

Objective: To detect gene mutation in patients with hypohidrotic ectodermal dysplasia (HED) by using whole exome sequencing, to analyze the pathogenicity of the mutations, and to provide reference for the genetic diagnosis of HED patients. Methods: Peripheral blood genomic DNA was extracted from each of the HED patients and their family members collected in Peking University School and Hospital of Stomatology from August 2016 to August 2021. Whole exome sequencing and sanger sequencing were performed to detect gene mutations. Functions of the rare variants after the database filtering were analyzed by bioinformatics tools. Results: Three reported mutations of ectodysplasin A (EDA) gene (c.2T>C, c.161A>G, c.467G>A) and a mutation of ectodysplasin A receptor (EDAR) gene (c.871G>A) were detected by whole genome sequencing in four HED patients, and were verified by Sanger sequencing in four HED families. The EDAR gene mutation founded in this research was reported in HED patients for the first time. Bioinformatics tools predicted that the mutations of EDA gene detected in this study were highly species conserved and disease-causing. The combined annotation dependent depletion (CADD) scores of EDA gene mutations c.2T>C, c.161A>G and c.467G>A were 22.5, 26.3 and 25.5 respectively, and the genomic evolutionary rate profiling (GERP) scores were 2.16, 2.26 and 2.18 respectively. The EDAR gene mutation c.871G>A detected in this study was species conserved and possibly disease-causing. The CADD and GERP scores of EDAR gene mutation c.871G>A were 22.0 and 1.93 respectively. Conclusions: Three reported mutations of EDA gene and a previously unreported mutation of EDAR gene were detected in four HED families. Different mutations of EDA gene and EDAR gene could make different influence on the protein function and lead to the occurrence of HED.

A Machine Learning Approach to Identifying Causal Monogenic Variants in Inflammatory Bowel Disease.

Diagnosis of monogenic disease is increasingly important for patient care and personalizing therapy. However, the current process is nonstandardized, expensive, and time consuming. There is currently no accepted strategy to help identify disease-causing variants in monogenic inflammatory bowel disease (IBD). The aim of the study is to develop a prioritization strategy for monogenic IBD variant discovery through detailed analysis of a whole-exome sequencing (WES) data set. All consenting pediatric patients with IBD presenting to our tertiary care hospital during the study period were enrolled and underwent WES (n= 1005). Available family members also underwent WES. Variants were analyzed en masse using the GEMINI framework and were further annotated using data from dbNSFP, Combined Annotation Dependent Depletion, and gnomAD. Known disease-causing variants (n= 36) were used as positive controls. Machine learning algorithms were optimized and then compared to assist with identifying monogenic IBD case characteristics. Initial gene-level analysis identified 11 genes not previously linked to IBD that could potentially harbor IBD-causing variants. Machine learning algorithms identified 4 primary variant characteristics (Combined Annotation Dependent Depletion score, dbNSFP score, relationship with a known immunodeficiency gene, and alternate allele frequency), and optimal threshold values for each were determined to assist with identifying monogenic IBD variants. Based on these characteristics, an automated variant prioritization pipeline was then created that filters and prioritizes variants from >100,000 variants per patient down to a mean of 15. This pipeline is available online for all to use. Leveraging a large WES data set, we demonstrate a statistically rigorous strategy for prioritization of variants for monogenic IBD diagnosis.

A Floppy Infant with Facial Dysmorphism.

A Floppy Infant with Facial Dysmorphism.

In-vitro characterization of coding variants with predicted functional implications in the efflux transporter multidrug resistance protein 4 (MRP4, ABCC4).

MRP4 (gene ABCC4) is a polymorphic efflux transporter that has been implicated in drug-induced toxicity. We selected ten commonly observed MRP4 coding variants among Europeans for experimental characterization including nine variants predicted to be deleterious or functional (combined annotation-dependent depletion score >15). We assessed protein localization and activity by quantifying intracellular accumulation of two prototypic substrates, taurocholic acid (TCA) and estradiol 17-β-glucuronide (E217βG), in HEK293T over-expressing MRP4 wildtype or variant where cellular substrate loading was optimized through co-transfection with an uptake transporter. V458M, a novel variant not previously studied, and T1142M, showed reduced activity compared to MRP4 wildtype for E217βG and TCA (P < 0.01), while L18I, G187W, K293E, and R531Q moderately increased activity in a substrate-dependent manner. Protein expression analysis indicated reduced cell surface expression for V458M (P < 0.01) but not T1142M compared to wildtype. Reduced activity may result from altered surface expression (V458M) or intrinsic activity as both variants map within the nucleotide-binding domains of MRP4. G187W showed a trend for reduced surface expression (P = 0.054) despite transport comparable or increased to wildtype suggesting enhanced intrinsic activity. Our findings suggest moderately altered MRP4 activity in six out of nine predicted functional variants with likely different mechanisms and substrate-specific effects. Cell-based studies using multiple known substrates are warranted to more accurately predict functional variants in this clinically important transporter.

Integration of Scandinavian genetic data with UK biobank data implicates the RBM20 gene with atrial fibrillation pathogenesis

Abstract Purpose Atrial fibrillation (AF) is the most common sustained arrhythmia. It carries a large healthcare burden and is associated with serious complications. The arrhythmia has a substantial genetic component and is associated with several structural genes, including the gene TTN. A recent large genome-wide association study on AF found an association to RBM20. The RBM20 gene is a splicing factor targeting TTN, RYR2 and CAMK2D among other cardiac genes. Using Next-Generation Sequencing and data derived from the UK Biobank, we aimed to reveal the role of RBM20 in AF. Methods and results We examined the burden of rare (Minor allele frequency (MAF)&amp;lt;0.01%) RBM20 loss-of-function (LOF) variants in whole-exome sequencing data from the UK Biobank (n=175,280). AF was defined by ICD9/10, while individuals without AF were used as controls. Association tests aggregating rare variants in RBM20 using the Efficient Variant-Set Mixed Model Association Test (SMMAT) were performed to assess the effect of LOF RBM20 variants, adjusted for age, sex and principal components. We identified 33 LOF variants in RBM20, which were significantly enriched in AF (P=0.0087). To examine the effect of rare missense RBM20 variants in the splicing of TTN, we screened an in-house cohort of 531 Scandinavian early-onset AF patients using targeted sequencing. We filtered for rare (MAF&amp;lt;0.1%) and deleterious (defined as combined annotation dependent depletion score &amp;gt;20) variants and identified nine missense variants and three novel LOF variants in RBM20. To evaluate the effect of these RBM20 variants, we constructed a series of human RBM20 single nucleotide base exchange mutants. The splicing activity of the variants was measured with RT-qPCR on HEK293 cells transfected with a TTN241–3 splicing reporter. Four of these variants resulted in a significantly altered splicing activity in TTN, with the largest effect observed for LOF variants. In order to examine the biological effect of RBM20 variants on structural changes in atrial tissue, we used a Norwegian Brown rat animal model with loss of RBM20. In this model, Transmission Electron Microscopy revealed altered sarcomere and mitochondrial structure in its atrial cardiomyocytes. Furthermore, nanopore RNA sequencing of atrial tissue from the aforementioned animal model indicated altered expression in several key cardiac genes, including TTN and PITX2. Conclusion Rare RBM20 LOF variants are significantly enriched in AF cases, seen in a large population of 175,000 individuals. We demonstrated that the effect of LOF RBM20 on alternative TTN splicing can be detected on an individual level in patients with AF. Studies using an animal model indicates that LOF in RBM20 may affect atrial function through altered expression of several genes in the atria, and may cause structural changes in the atrial cardiomyocytes. This suggests that RBM20 may be involved in AF pathogenesis mediated through an atrial cardiomyopathy. Funding Acknowledgement Type of funding sources: Foundation. Main funding source(s): Novo Nordisk Foundation Pre-Graduate Scholarships (NNF18OC0053094)The Hallas Møller Emerging Investigator grant (Novo Nordisk Foundation (NNF17OC0031204))