Dissecting the polygenic basis of atherosclerosis via disease-associated cell state signatures.

Coronary artery disease (CAD) is a pandemic disease where up to half of the risk is explained by genetic factors. Advanced insights into the genetic basis of CAD require deeper understanding of the contributions of different cell types, molecular pathways, and genes to disease heritability. Here, we investigate the biological diversity of atherosclerosis-associated cell states and interrogate their contribution to the genetic risk of CAD by using single-cell and bulk RNA sequencing (RNA-seq) of mouse and human lesions. We identified 12 disease-associated cell states that we characterized further by gene set functional profiling, ligand-receptor prediction, and transcription factor inference. Importantly, Vcam1+ smooth muscle cell state genes contributed most to SNP-based heritability of CAD. In line with this, genetic variants near smooth muscle cell state genes and regulatory elements explained the largest fraction of CAD-risk variance between individuals. Using this information for variant prioritization, we derived a hybrid polygenic risk score (PRS) that demonstrated improved performance over a classical PRS. Our results provide insights into the biological mechanisms associated with CAD risk, which could make a promising contribution to precision medicine and tailored therapeutic interventions in the future.

Late-onset and classic phenotypes of Fabry disease in males with the GLA-Thr410Ala mutation.

OBJECTIVE: To present phenotypic characteristics and biomarkers of a family with the rare mutation Thr410Ala of the α-galactosidase A gene (T410A/GLA) causing Fabry disease (FD). METHODS AND RESULTS: In a woman in her 60s with hypertrophic cardiomyopathy, T410A/GLA was found in screening for variants in 59 cardiomyopathy-related genes. Her son in his 40s, two granddaughters and two great grandsons carried T410A/GLA. The son had a history of hypertension and paroxysmal AF but no microalbuminuria or classic symptoms or signs of FD. Baseline α-galactosidase A enzyme (α-Gal A) activity varied from 0% to 26.5%. Cardiac MRI showed mild Fabry cardiomyopathy (FC). During 11 years of enzyme replacement therapy (ERT), FC progressed and he suffered sudden cardiac death in his 50s. The great grandsons with T410A/GLA had no active α-Gal A, high lyso-Gb3 levels and normal cardiac imaging. They suffered from neuropathic pain and gastrointestinal symptoms and were started with ERT at the age under 10. Granddaughters with T410A/GLA had α-Gal A activities of 8-18 and 10% of normal. The older granddaughter in her 30s was diagnosed with incipient FC. Plasma lyso-Gb3 analogues were elevated, markedly in the elder male with FC and moderately in the elder granddaughter. In young males with classic phenotype, plasma lyso-Gb3 analogues were only slightly elevated. CONCLUSIONS: The T410A/GLA mutation caused late-onset FD with progressive cardiomyopathy in elder male, and classic FD in young males of the same family. Varying levels of α-Gal A and lyso-Gb3 analogues reflected variable phenotype of FD in the family.

Role of 19 SNPs in 10 genes with type 2 diabetes in the Pakistani population.

Meta-GWAS report numerous variants associated with type 2 diabetes (T2D), however, for diabetic retinopathy (DR) no loci achieved genome-wide significance. There are limited candidate gene analyses for T2D and/or DR reported from the Pakistani population. Therefore, the current study was designed to evaluate the genetic association of 10 loci with T2D, non-proliferative DR (NPDR), and proliferative DR (PDR). In total 375 T2D cases and 205 controls were collected. The T2D cases included diabetic no retinopathy (n = 196), NPDR (n = 95), and PDR (n = 84). Genomic DNA was isolated, and 19 SNPs were genotyped. To determine association of SNPs with T2D, logistic regression analyses were performed adjusting for age and sex. Moreover, for association of SNPs with NPDR and PDR logistic regression analyses adjusting for diabetes duration and age of T2D onset were performed. In multivariate analysis, the minor alleles of rs1043618 [G > C, odds ratio (OR) 95% confidence interval (CI) = 1.45 (1.13-1.87), p = 4.00E-3], rs3807987 [G > A, 1.87 (1.22-2.94), p = 0.01], rs12672038 [G > A, 1.53 (1.04-2.30), p = 0.03] and rs2055858 [G > C, 1.70 (1.20-2.43), p = 3.00E-3] were associated with higher risk while rs1801133 (C > T, 0.59 (0.42-0.83), p = 2.28E-3) was associated with a lower risk of T2D. Moreover, minor alleles of rs2055858 [G > C, 1.77 (1.17-2.68), p = 0.02], and rs3759890 [C > G, 2.17 (1.39-3.39), p = 4.00E-3] showed an association with PDR when compared with DNR. However, only the association of rs1801133 survived multiple test correction. Hence, we report that rs1801133 is associated with T2D in the Pakistani population. In addition, out of studied 10 genes 8 proteins had higher interactions among themselves that are predicted to be partially biologically connected, as a group.

Metabolite Signature in the Carriers of Pathogenic Genetic Variants for Cardiomyopathy: A Population-Based METSIM Study.

Hypertrophic (HCM) and dilated (DCM) cardiomyopathies are among the leading causes of sudden cardiac death. We identified 38 pathogenic or likely pathogenic variant carriers for HCM in three sarcomere genes (MYH7, MYBPC3, TPMI) among 9.928 participants of the METSIM Study having whole exome sequencing data available. Eight of them had a clinical diagnosis of HCM. We also identified 20 pathogenic or likely pathogenic variant carriers for DCM in the TTN gene, and six of them had a clinical diagnosis of DCM. The aim of our study was to investigate the metabolite signature in the carriers of the pathogenic or likely pathogenic genetic variants for HCM and DCM, compared to age- and body-mass-index-matched controls. Our novel findings were that the carriers of pathogenic or likely pathogenic variants for HCM had significantly increased concentrations of bradykinin (des-arg 9), vanillactate, and dimethylglycine and decreased concentrations of polysaturated fatty acids (PUFAs) and lysophosphatidylcholines compared with the controls without HCM. Additionally, our novel findings were that the carriers of pathogenic or likely pathogenic variants for DCM had significantly decreased concentrations of 1,5-anhydrogluticol, histidine betaine, N-acetyltryptophan, and methylsuccinate and increased concentrations of trans-4-hydroxyproline compared to the controls without DCM. Our population-based study shows that the metabolite signature of the genetic variants for HCM and DCM includes several novel metabolic pathways not previously described.

The Power of Single-Cell RNA Sequencing in eQTL Discovery.

Genome-wide association studies have successfully mapped thousands of loci associated with complex traits. During the last decade, functional genomics approaches combining genotype information with bulk RNA-sequencing data have identified genes regulated by GWAS loci through expression quantitative trait locus (eQTL) analysis. Single-cell RNA-Sequencing (scRNA-Seq) technologies have created new exciting opportunities for spatiotemporal assessment of changes in gene expression at the single-cell level in complex and inherited conditions. A growing number of studies have demonstrated the power of scRNA-Seq in eQTL mapping across different cell types, developmental stages and stimuli that could be obscured when using bulk RNA-Seq methods. In this review, we outline the methodological principles, advantages, limitations and the future experimental and analytical considerations of single-cell eQTL studies. We look forward to the explosion of single-cell eQTL studies applied to large-scale population genetics to take us one step closer to understanding the molecular mechanisms of disease.

Clinical and Genetic Characterization of 153 Patients with Persistent or Transient Congenital Hyperinsulinism.

CONTEXT: Major advances have been made in the genetics and classification of congenital hyperinsulinism (CHI). OBJECTIVE: To examine the genetics and clinical characteristics of patients with persistent and transient CHI. DESIGN: A cross-sectional study with the register data and targeted sequencing of 104 genes affecting glucose metabolism. PATIENTS: Genetic and phenotypic data were collected from 153 patients with persistent (n = 95) and transient (n = 58) CHI diagnosed between 1972 and 2015. Of these, 86 patients with persistent and 58 with transient CHI participated in the analysis of the selected 104 genes affecting glucose metabolism, including 10 CHI-associated genes, and 9 patients with persistent CHI were included because of their previously confirmed genetic diagnosis. MAIN OUTCOME MEASURES: Targeted next-generation sequencing results and genotype-phenotype associations. RESULTS: Five novel and 21 previously reported pathogenic or likely pathogenic variants in ABCC8, KCNJ11, GLUD1, GCK, HNF4A, and SLC16A1 genes were found in 68% (n = 65) and 0% of the patients with persistent and transient CHI, respectively. KATP channel mutations explained 82% of the mutation positive cases. CONCLUSIONS: The genetic variants found in this nationwide CHI cohort are in agreement with previous studies, mutations in the KATP channel genes being the major causes of the disease. Pathogenic CHI-associated variants were not identified in patients who were both diazoxide responsive and able to discontinue medication within the first 4 months. Therefore, our results support the notion that genetic testing should be focused on patients with inadequate response or prolonged need for medication.

Cardiomyopathy associated with the Ala143Thr variant of the α-galactosidase A gene.

OBJECTIVE: To investigate whether the Ala143Thr variant of the α-galactosidase A gene (A143T/GLA), with conflicting interpretations of pathogenicity, is associated with Fabry cardiomyopathy. METHODS: The index patient, a woman in her 60s with cardiomyopathy, was screened for variants in 59 cardiomyopathy-related genes. A143T/GLA, the only rare variant found, was screened in 10 relatives. GLA activity and lyso-Gb3 levels were measured and echocardiography was performed in 8 of 9 subjects carrying A143T/GLA. Cardiac magnetic resonance (CMR) imaging and 18F-fluorodeoxyglucose (FDG) positron emission tomography/CT (PET/CT) were performed in four adult A143T/GLA carriers. Endomyocardial biopsy was obtained from two adult A143T/GLA carrying sons of the index patient. RESULTS: The index patient and her elder son had a pacemaker implantation because of sick sinus syndrome and atrioventricular block. GLA activities were decreased to 25%-40% of normal in both sons and one granddaughter. Lyso-Gb3 levels were elevated in both sons. In CMR, the index patient and her two sons had left ventricular (LV) hypertrophy and/or dilatation. The elder son had late gadolinium enhancement, high CMR-derived T1 time and positive FDG signal in PET/CT in the basal inferolateral LV wall. The younger son had low T1 time and the mother had positive FDG signal in PET/CT in the basal inferolateral LV wall. Endomyocardial biopsy of both sons showed myocardial accumulation compatible with glycolipids in light and electron microscopy, staining with anti-Gb3 antibody available for the younger son. Five female relatives with A143T/GLA had no cardiomyopathy in cardiac imaging. CONCLUSIONS: A143T/GLA is likely a late-onset Fabry cardiomyopathy causing variant with incomplete penetrance.

Genetic and clinical characterization of Pakistani families with Bardet-Biedl syndrome extends the genetic and phenotypic spectrum.

Bardet-Biedl syndrome (BBS) is an autosomal recessive disorder that is both genetically and clinically heterogeneous. To date 19 genes have been associated with BBS, which encode proteins active at the primary cilium, an antenna-like organelle that acts as the cell's signaling hub. In the current study, a combination of mutation screening, targeted sequencing of ciliopathy genes associated with BBS, and whole-exome sequencing was used for the genetic characterization of five families including four with classic BBS symptoms and one BBS-like syndrome. This resulted in the identification of novel mutations in BBS genes ARL6 and BBS5, and recurrent mutations in BBS9 and CEP164. In the case of CEP164, this is the first report of two siblings with a BBS-like syndrome with mutations in this gene. Mutations in this gene were previously associated with nephronophthisis 15, thus the current results expand the CEP164-associated phenotypic spectrum. The clinical and genetic spectrum of BBS and BBS-like phenotypes is not fully defined in Pakistan. Therefore, genetic studies are needed to gain insights into genotype-phenotype correlations, which will in turn improve the clinician's ability to make an early and accurate diagnosis, and facilitate genetic counseling, leading to directly benefiting families with affected individuals.

Homozygosity mapping and targeted sanger sequencing reveal genetic defects underlying inherited retinal disease in families from pakistan.

BACKGROUND: Homozygosity mapping has facilitated the identification of the genetic causes underlying inherited diseases, particularly in consanguineous families with multiple affected individuals. This knowledge has also resulted in a mutation dataset that can be used in a cost and time effective manner to screen frequent population-specific genetic variations associated with diseases such as inherited retinal disease (IRD). METHODS: We genetically screened 13 families from a cohort of 81 Pakistani IRD families diagnosed with Leber congenital amaurosis (LCA), retinitis pigmentosa (RP), congenital stationary night blindness (CSNB), or cone dystrophy (CD). We employed genome-wide single nucleotide polymorphism (SNP) array analysis to identify homozygous regions shared by affected individuals and performed Sanger sequencing of IRD-associated genes located in the sizeable homozygous regions. In addition, based on population specific mutation data we performed targeted Sanger sequencing (TSS) of frequent variants in AIPL1, CEP290, CRB1, GUCY2D, LCA5, RPGRIP1 and TULP1, in probands from 28 LCA families. RESULTS: Homozygosity mapping and Sanger sequencing of IRD-associated genes revealed the underlying mutations in 10 families. TSS revealed causative variants in three families. In these 13 families four novel mutations were identified in CNGA1, CNGB1, GUCY2D, and RPGRIP1. CONCLUSIONS: Homozygosity mapping and TSS revealed the underlying genetic cause in 13 IRD families, which is useful for genetic counseling as well as therapeutic interventions that are likely to become available in the near future.

Screening of a large cohort of leber congenital amaurosis and retinitis pigmentosa patients identifies novel LCA5 mutations and new genotype-phenotype correlations.

This study was undertaken to investigate the prevalence of sequence variants in LCA5 in patients with Leber congenital amaurosis (LCA), early-onset retinal dystrophy (EORD), and autosomal recessive retinitis pigmentosa (arRP); to delineate the ocular phenotypes; and to provide an overview of all published LCA5 variants in an online database. Patients underwent standard ophthalmic evaluations after providing informed consent. In selected patients, optical coherence tomography (OCT) and fundus autofluorescence imaging were possible. DNA samples from 797 unrelated patients with LCA and 211 with the various types of retinitis pigmentosa (RP) were screened by Sanger sequence analysis of all LCA5 exons and intron/exon junctions. Some LCA patients were prescreened by APEX technology or selected based on homozygosity mapping. In silico analyses were performed to assess the pathogenicity of the variants. Segregation analysis was performed where possible. Published and novel LCA5 variants were collected, amended for their correct nomenclature, and listed in a Leiden Open Variation Database (LOVD). Sequence analysis identified 18 new probands with 19 different LCA5 variants. Seventeen of the 19 LCA5 variants were novel. Except for two missense variants and one splice site variant, all variants were protein-truncating mutations. Most patients expressed a severe phenotype, typical of LCA. However, some LCA subjects had better vision and intact inner segment/outer segment (IS/OS) junctions on OCT imaging. In two families with LCA5 variants, the phenotype was more compatible with EORD with affected individuals displaying preserved islands of retinal pigment epithelium. One of the families with a milder phenotype harbored a homozygous splice site mutation; a second family was found to have a combination of a stop mutation and a missense mutation. This is the largest LCA5 study to date. We sequenced 1,008 patients (797 with LCA, 211 with arRP) and identified 18 probands with LCA5 mutations. Mutations in LCA5 are a rare cause of childhood retinal dystrophy accounting for ∼2% of disease in this cohort, and the majority of LCA5 mutations are likely null. The LCA5 protein truncating mutations are predominantly associated with LCA. However, in two families with the milder EORD, the LCA5 gene analysis revealed a homozygous splice site mutation in one and a stop mutation in combination with a missense mutation in a second family, suggesting that this milder phenotype is due to residual function of lebercilin and expanding the currently known phenotypic spectrum to include the milder early onset RP. Some patients have remaining foveal cone structures (intact IS/OS junctions on OCT imaging) and remaining visual acuities, which may bode well for upcoming treatment trials.

A nonsense mutation in S-antigen (p.Glu306*) causes Oguchi disease.

PURPOSE: Genetic studies were performed to identify the causative mutation in a 15-year-old girl diagnosed with congenital stationary night blindness (CSNB) presenting Mizuo-Nakamura phenomenon, a typical Oguchi disease symptom. The patient also had dural sinus thrombosis (DST), thrombocytopenia, and systemic lupus erythematosus (SLE). METHODS: Mutation analysis was done by sequencing two candidate genes, S-antigen (SAG; arrestin 1), associated with Oguchi type 1, and rhodopsin kinase (GRK1), associated with Oguchi type 2. In addition, the C677T variation in the methylenetetrahydrofolate reductase (MTHFR) gene was also screened in the family, to determine its probable association with hyperhomocysteinemia in the patient. RESULTS: Sequencing of the SAG and GRK1 resulted in identifying a novel homozygous nonsense mutation (c.916G>T; p.Glu306*) in SAG, which in unaffected siblings either was present in a heterozygous state or absent. The C677T heterozygous allele in the MTHFR gene was found to be associated with hyperhomocysteinemia in the patient and other family members. CONCLUSIONS: This is the first report of Oguchi type 1 in a Pakistani patient due to a nonsense mutation (c.916G>T; p.Glu306*) in SAG. The neurologic and hematological abnormalities likely are not associated with the SAG variant.