Human rare variant and zebrafish CRISPR/Cas9-mediated mutant analyses reveal novel functions for API5, HSPB7, and LMO2 in heart failure

The clinical heterogeneity of heart failure has challenged our understanding of the underlying genetic mechanisms of this disease. To gain further insights into this complex pathophysiology we combined human rare variant analysis and in vivo CRISPR/Cas9-mediated mutant phenotyping in zebrafish to identify and investigate the role of 3 genes. Whole-exome sequencing of patients identified API5, HSPB7, and LMO2 as causally associated with heart failure and these genes were further investigated, alongside the positive control gata5, using CRISPR/Cas9-mediated multi-locus in vivo mutation in zebrafish. Following effective somatic mutation, we observed multiple impacts on cardiovascular development and function in F0 embryos including reductions in ventricle size, pericardial oedema, and chamber malformation. In the case of lmo2, there was also a significant impact on heart function. Our analysis suggests novel functions for API5, HSPB7, and LMO2 in human cardiovascular disease and identifies them as potential drug targets. Our data also supports in vivo CRISPR/Cas9-mediated multi-locus gene mutation analysis in F0 zebrafish as a rapid and effective primary screen for assessing gene function, as part of an integrated multi-level drug target validation strategy.

STAAR Workflow: A cloud-based workflow for scalable and reproducible rare variant analysis

SummaryWe developed the STAAR WDL workflow to facilitate the analysis of rare variants in whole genome sequencing association studies. The open-access STAAR workflow written in the workflow description language (WDL) allows a user to perform rare variant testing for both gene-centric and genetic region approaches, enabling genome-wide, candidate, and conditional analyses. It incorporates functional annotations into the workflow as introduced in the STAAR method in order to boost the rare variant analysis power. This tool was specifically developed and optimized to be implemented on cloud-based platforms such as BioData Catalyst Powered by Terra. It provides easy-to-use functionality for rare variant analysis that can be incorporated into an exhaustive whole genome sequencing analysis pipeline.Availability and implementationThe workflow is freely available from https://dockstore.org/workflows/github.com/sheilagaynor/STAAR_workflow.

Rare variants in IFFO1, DTNB and NLRC3 associate with Alzheimer's disease CSF profile of neuronal injury and inflammation

Alzheimer's disease (AD) biomarkers represent several neurodegenerative processes, such as synaptic dysfunction, neuronal inflammation and injury, as well as amyloid pathology. We performed an exome-wide rare variant analysis of six AD biomarkers (β-amyloid, total/phosphorylated tau, Nfl, YKL-40, and Neurogranin) to discover genes associated with these markers. Genetic and biomarker information was available for 480 participants from two studies: EMIF-AD and ADNI. We applied a principal component (PC) analysis to derive biomarkers combinations, which represent statistically independent biological processes. We then tested whether rare variants in 9,576 protein-coding genes associate with these PCs using a Meta-SKAT test. We also tested whether the PCs are intermediary to gene effects on AD symptoms with a SMUT test. One PC loaded on Nfl and YKL-40, indicators of neuronal injury and inflammation. Three genes were associated with this PC: IFFO1, DTNB and NLRC3. Mediation tests suggest, that these genes also affect dementia symptoms via inflammation/injury. We also observed an association between a PC loading on Neurogranin, a marker for synaptic functioning, with GABBR2 and CASZ1, but no mediation effects. The results suggest that rare variants in IFFO1, DTNB and NLRC3 heighten susceptibility to neuronal injury and inflammation, potentially by altering cytoskeleton structure and immune activity disinhibition, resulting in an elevated dementia risk. GABBR2 and CASZ1 were associated with synaptic functioning, but mediation analyses suggest that the effect of these two genes on synaptic functioning is not consequential for AD development.

Targeted sequencing of NOTCH signaling pathway genes and association analysis of variants correlated with mandibular prognathism

Introduction The purpose of this study was to systematically identify variants in NOTCH signaling pathway genes that correlate with mandibular prognathism (MP) in the general Chinese population. Methods Targeted sequencing of NOTCH signaling pathway genes was conducted in 199 MP individuals and 197 class I malocclusion control individuals. The associations of common and rare variants with MP, cephalometric parameters, and continuous cephalometric phenotypes were analyzed by principal component (PC) analysis. The associations between rare variants and MP were tested for each gene. Results Six SNPs, including rs415929, rs520688, and rs423023 in an exonic region of NOTCH4; rs1044006 in an exonic region of NOTCH3; rs1051415 in an exonic region of JAG1; and rs75236173 in the 3′-untranslated region (3′-UTR) of NUMB were associated with MP (P < 0.05). One common variant, rs1051415, in an exonic region of JAG1 was significantly related to PC1 (P = 3.608 × 10− 4), which explained 24.3% of the overall phenotypic variation observed and corresponded to the sagittal mandibular position towards the maxilla, ranging from a posterior positioned mandible to an anterior positioned mandible. Additionally, 41 other variants were associated with PC1–5 (P < 0.05). With respect to rare variant analysis, variants within the EP300, NCOR2, and PSEN2 gene showed an association with MP (t < 0 .05). Conclusions An association between NOTCH signaling pathway genes and MP has been identified.

Rare Missense Functional Variants at COL4A1 and COL4A2 in Sporadic Intracerebral Hemorrhage

ObjectiveTo test the genetic contribution of rare missense variants in COL4A1 and COL4A2 in which common variants are genetically associated with sporadic intracerebral hemorrhage (ICH), we performed rare variant analysis in multiple sequencing data for the risk for sporadic ICH.MethodsWe performed sequencing across 559Kbp at 13q34 including COL4A1 and COL4A2 among 2,133 individuals (1,055 ICH cases; 1,078 controls) in US-based and 1,492 individuals (192 ICH cases; 1,189 controls) from Scotland-based cohorts, followed by sequence annotation, functional impact prediction, genetic association testing, and in silico thermodynamic modeling.ResultsWe identified 107 rare nonsynonymous variants in sporadic ICH, of which two missense variants, rs138269346 (COL4A1I110T) and rs201716258 (COL4A2H203L), were predicted to be highly functional and occurred in multiple ICH cases but not in controls from the US-based cohort. The minor allele of rs201716258 was also present in Scottish ICH patients, and rs138269346 was observed in two ICH-free controls with a history of hypertension and myocardial infarction. Rs138269346 was nominally associated with non-lobar ICH risk (P=0.05), but not with lobar ICH (P=0.08), while associations between rs201716258 and ICH subtypes were non-significant (P>0.12). Both variants were considered pathogenic based on minor allele frequency (<0.00035 in EUR), predicted functional impact (deleterious or probably damaging), and in silico modeling studies (substantially altered physical length and thermal stability of collagen).ConclusionsWe identified rare missense variants in COL4A1/A2 in association with sporadic ICH. Our annotation and simulation studies suggest that these variants are highly functional and may represent targets for translational follow-up.

Rare variant analysis of 4241 pulmonary arterial hypertension cases from an international consortium implicates FBLN2, PDGFD, and rare de novo variants in PAH

Background Pulmonary arterial hypertension (PAH) is a lethal vasculopathy characterized by pathogenic remodeling of pulmonary arterioles leading to increased pulmonary pressures, right ventricular hypertrophy, and heart failure. PAH can be associated with other diseases (APAH: connective tissue diseases, congenital heart disease, and others) but often the etiology is idiopathic (IPAH). Mutations in bone morphogenetic protein receptor 2 (BMPR2) are the cause of most heritable cases but the vast majority of other cases are genetically undefined. Methods To identify new risk genes, we utilized an international consortium of 4241 PAH cases with exome or genome sequencing data from the National Biological Sample and Data Repository for PAH, Columbia University Irving Medical Center, and the UK NIHR BioResource – Rare Diseases Study. The strength of this combined cohort is a doubling of the number of IPAH cases compared to either national cohort alone. We identified protein-coding variants and performed rare variant association analyses in unrelated participants of European ancestry, including 1647 IPAH cases and 18,819 controls. We also analyzed de novo variants in 124 pediatric trios enriched for IPAH and APAH-CHD. Results Seven genes with rare deleterious variants were associated with IPAH with false discovery rate smaller than 0.1: three known genes (BMPR2, GDF2, and TBX4), two recently identified candidate genes (SOX17, KDR), and two new candidate genes (fibulin 2, FBLN2; platelet-derived growth factor D, PDGFD). The new genes were identified based solely on rare deleterious missense variants, a variant type that could not be adequately assessed in either cohort alone. The candidate genes exhibit expression patterns in lung and heart similar to that of known PAH risk genes, and most variants occur in conserved protein domains. For pediatric PAH, predicted deleterious de novo variants exhibited a significant burden compared to the background mutation rate (2.45×, p = 2.5e−5). At least eight novel pediatric candidate genes carrying de novo variants have plausible roles in lung/heart development. Conclusions Rare variant analysis of a large international consortium identified two new candidate genes—FBLN2 and PDGFD. The new genes have known functions in vasculogenesis and remodeling. Trio analysis predicted that ~ 15% of pediatric IPAH may be explained by de novo variants.

GWAS and ExWAS of blood Mitochondrial DNA copy number identifies 73 loci and highlights a potential causal role in dementia

Mitochondrial DNA copy number (mtDNA-CN) is an accessible blood-based measurement believed to capture underlying mitochondrial function. The specific biological processes underpinning its regulation, and whether those processes are causative for disease, is an area of active investigation. We developed a novel method for array-based mtDNA-CN estimation suitable for biobank-scale studies, called “AutoMitoC”. We applied AutoMitoC to 395,781 UKBiobank study participants and performed genome and exome-wide association studies, identifying novel common and rare genetic determinants. Overall, genetic analyses identified 73 loci for mtDNA-CN, which implicated several genes involved in rare mtDNA depletion disorders, dNTP metabolism, and the mitochondrial central dogma. Rare variant analysis identified SAMHD1 mutation carriers as having higher mtDNA-CN (beta=0.23 SDs; 95% CI, 0.18-0.29; P=2.6×10−19), a potential therapeutic target for patients with mtDNA depletion disorders, but at increased risk of breast cancer (OR=1.91; 95% CI, 1.52-2.40; P=2.7×10−8). Finally, Mendelian randomization analyses suggest a causal effect of low mtDNA-CN on dementia risk (OR=1.94 per 1 SD decrease in mtDNA-CN; 95% CI, 1.55-2.32; P=7.5×10−4). Altogether, our genetic findings indicate that mtDNA-CN is a complex biomarker reflecting specific mitochondrial processes related to mtDNA regulation, and that these processes are causally related to human diseases.

Discovery of predictors of sudden cardiac arrest in diabetes: rationale and outline of the RESCUED (REcognition of Sudden Cardiac arrest vUlnErability in Diabetes) project

IntroductionEarly recognition of individuals with increased risk of sudden cardiac arrest (SCA) remains challenging. SCA research so far has used data from cardiologist care, but missed most SCA victims, since they were only in general practitioner (GP) care prior to SCA. Studying individuals with type 2 diabetes (T2D) in GP care may help solve this problem, as they have increased risk for SCA, and rich clinical datasets, since they regularly visit their GP for check-up measurements. This information can be further enriched with extensive genetic and metabolic information.AimTo describe the study protocol of the REcognition of Sudden Cardiac arrest vUlnErability in Diabetes (RESCUED) project, which aims at identifying clinical, genetic and metabolic factors contributing to SCA risk in individuals with T2D, and to develop a prognostic model for the risk of SCA.MethodsThe RESCUED project combines data from dedicated SCA and T2D cohorts, and GP data, from the same region in the Netherlands. Clinical data, genetic data (common and rare variant analysis) and metabolic data (metabolomics) will be analysed (using classical analysis techniques and machine learning methods) and combined into a prognostic model for risk of SCA.ConclusionThe RESCUED project is designed to increase our ability at early recognition of elevated SCA risk through an innovative strategy of focusing on GP data and a multidimensional methodology including clinical, genetic and metabolic analyses.