Phenotype-Specific Enrichment of Mendelian Disorder Genes near GWAS Regions across 62 Complex Traits

Non-linear transcriptional responses to gradual modulation of transcription factor dosage

Genomic loci associated with common traits and diseases are typically non-coding and likely impact gene expression, sometimes coinciding with rare loss-of-function variants in the target gene. However, our understanding of how gradual changes in gene dosage affect molecular, cellular, and organismal traits is currently limited. To address this gap, we induced gradual changes in gene expression of four genes using CRISPR activation and inactivation. Downstream transcriptional consequences of dosage modulation of three master trans-regulators associated with blood cell traits (GFI1B, NFE2, and MYB) were examined using targeted single-cell multimodal sequencing. We showed that guide tiling around the TSS is the most effective way to modulate cis gene expression across a wide range of fold-changes, with further effects from chromatin accessibility and histone marks that differ between the inhibition and activation systems. Our single-cell data allowed us to precisely detect subtle to large gene expression changes in dozens of trans genes, revealing that many responses to dosage changes of these three TFs are non-linear, including non-monotonic behaviours, even when constraining the fold-changes of the master regulators to a copy number gain or loss. We found that the dosage properties are linked to gene constraint and that some of these non-linear responses are enriched for disease and GWAS genes. Overall, our study provides a straightforward and scalable method to precisely modulate gene expression and gain insights into its downstream consequences at high resolution.

Analysis of Evolutionary Conservation, Expression Level, and Genetic Association at a Genome-wide Scale Reveals Heterogeneity Across Polygenic Phenotypes

Understanding the expression level and evolutionary rate of associated genes with human polygenic diseases provides crucial insights into their disease-contributing roles. In this work, we leveraged genome-wide association studies (GWASs) to investigate the relationship between the genetic association and both the evolutionary rate (dN/dS) and expression level of human genes associated with the two polygenic diseases of schizophrenia and coronary artery disease. Our findings highlight a distinct variation in these relationships between the two diseases. Genes associated with both diseases exhibit a significantly greater variance in evolutionary rate compared to those implicated in monogenic diseases. Expanding our analyses to 4,756 complex traits in the GWAS atlas database, we unraveled distinct trait categories with a unique interplay among the evolutionary rate, expression level, and genetic association of human genes. In most polygenic traits, highly expressed genes were more associated with the polygenic phenotypes compared to lowly expressed genes. About 69% of polygenic traits displayed a negative correlation between genetic association and evolutionary rate, while approximately 30% of these traits showed a positive correlation between genetic association and evolutionary rate. Our results demonstrate the presence of a spectrum among complex traits, shaped by natural selection. Notably, at opposite ends of this spectrum, we find metabolic traits being more likely influenced by purifying selection, and immunological traits that are more likely shaped by positive selection. We further established the polygenic evolution portal (evopolygen.de) as a resource for investigating relationships and generating hypotheses in the field of human polygenic trait evolution.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

Disentangling multiple sclerosis phenotypes through Mendelian disorders: A network approach

BACKGROUND

The increasing knowledge about multiple sclerosis (MS) pathophysiology has reinforced the need for an improved description of disease phenotypes, connected to disease biology. Growing evidence indicates that complex diseases constitute phenotypical and genetic continuums with "simple," monogenic disorders, suggesting shared pathomechanisms.

OBJECTIVES

The objective of this study was to depict a novel MS phenotypical framework leveraging shared physiopathology with Mendelian diseases and to identify phenotype-specific candidate drugs.

METHODS

We performed an enrichment testing of MS-associated variants with Mendelian disorders genes. We defined a "MS-Mendelian network," further analyzed to define enriched phenotypic subnetworks and biological processes. Finally, a network-based drug screening was implemented.

RESULTS

Starting from 617 MS-associated loci, we showed a significant enrichment of monogenic diseases (p < 0.001). We defined an MS-Mendelian molecular network based on 331 genes and 486 related disorders, enriched in four phenotypic classes: neurologic, immunologic, metabolic, and visual. We prioritized a total of 503 drugs, of which 27 molecules active in 3/4 phenotypical subnetworks and 140 in subnetwork pairs.

CONCLUSION

The genetic architecture of MS contains the seeds of pathobiological multiplicities shared with immune, neurologic, metabolic and visual monogenic disorders. This result may inform future classifications of MS endophenotypes and support the development of new therapies in both MS and rare diseases.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.

Speos: an ensemble graph representation learning framework to predict core gene candidates for complex diseases

Understanding phenotype-to-genotype relationships is a grand challenge of 21st century biology with translational implications. The recently proposed "omnigenic" model postulates that effects of genetic variation on traits are mediated by core-genes and -proteins whose activities mechanistically influence the phenotype, whereas peripheral genes encode a regulatory network that indirectly affects phenotypes via core gene products. Here, we develop a positive-unlabeled graph representation-learning ensemble-approach based on a nested cross-validation to predict core-like genes for diverse diseases using Mendelian disorder genes for training. Employing mouse knockout phenotypes for external validations, we demonstrate that core-like genes display several key properties of core genes: Mouse knockouts of genes corresponding to our most confident predictions give rise to relevant mouse phenotypes at rates on par with the Mendelian disorder genes, and all candidates exhibit core gene properties like transcriptional deregulation in disease and loss-of-function intolerance. Moreover, as predicted for core genes, our candidates are enriched for drug targets and druggable proteins. In contrast to Mendelian disorder genes the new core-like genes are enriched for druggable yet untargeted gene products, which are therefore attractive targets for drug development. Interpretation of the underlying deep learning model suggests plausible explanations for our core gene predictions in form of molecular mechanisms and physical interactions. Our results demonstrate the potential of graph representation learning for the interpretation of biological complexity and pave the way for studying core gene properties and future drug development.

Inborn Errors of Immunity in Hidradenitis Suppurativa Pathogenesis and Disease Burden

Hidradenitis suppurativa (HS), also known as Verneuil's disease and acne inversa, is a prevalent, debilitating, and understudied inflammatory skin disease. It is marked by repeated bouts of pathological inflammation causing pain, hyperplasia, aberrant healing, and fibrosis. HS is difficult to manage and has many unmet medical needs. There is clinical and pharmacological evidence for extensive etiological heterogeneity with HS, suggesting that this clinical diagnosis is capturing a spectrum of disease entities. Human genetic studies provide robust insight into disease pathogenesis. They also can be used to resolve etiological heterogeneity and to identify drug targets. However, HS has not been extensively investigated with well-powered genetic studies. Here, we review what is known about its genetic architecture. We identify overlap in molecular, cellular, and clinical features between HS and inborn errors of immunity (IEI). This evidence indicates that HS may be an underrecognized component of IEI and suggests that undiagnosed IEI are present in HS cohorts. Inborn errors of immunity represent a salient opportunity for rapidly resolving the immunological landscape of HS pathogenesis, for prioritizing drug repurposing studies, and for improving the clinical management of HS.

Shared genetic risk loci between Alzheimer's disease and related dementias, Parkinson's disease, and amyotrophic lateral sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have indicated moderate genetic overlap between Alzheimer's disease (AD) and related dementias (ADRD), Parkinson's disease (PD) and amyotrophic lateral sclerosis (ALS), neurodegenerative disorders traditionally considered etiologically distinct. However, the specific genetic variants and loci underlying this overlap remain almost entirely unknown.

METHODS

We leveraged state-of-the-art GWAS for ADRD, PD, and ALS. For each pair of disorders, we examined each of the GWAS hits for one disorder and tested whether they were also significant for the other disorder, applying Bonferroni correction for the number of variants tested. This approach rigorously controls the family-wise error rate for both disorders, analogously to genome-wide significance.

RESULTS

Eleven loci with GWAS hits for one disorder were also associated with one or both of the other disorders: one with all three disorders (the MAPT/KANSL1 locus), five with ADRD and PD (near LCORL, CLU, SETD1A/KAT8, WWOX, and GRN), three with ADRD and ALS (near GPX3, HS3ST5/HDAC2/MARCKS, and TSPOAP1), and two with PD and ALS (near GAK/TMEM175 and NEK1). Two of these loci (LCORL and NEK1) were associated with an increased risk of one disorder but decreased risk of another. Colocalization analysis supported a shared causal variant between ADRD and PD at the CLU, WWOX, and LCORL loci, between ADRD and ALS at the TSPOAP1 locus, and between PD and ALS at the NEK1 and GAK/TMEM175 loci. To address the concern that ADRD is an imperfect proxy for AD and that the ADRD and PD GWAS have overlapping participants (nearly all of which are from the UK Biobank), we confirmed that all our ADRD associations had nearly identical odds ratios in an AD GWAS that excluded the UK Biobank, and all but one remained nominally significant (p < 0.05) for AD.

CONCLUSIONS

In one of the most comprehensive investigations to date of pleiotropy between neurodegenerative disorders, we identify eleven genetic risk loci shared among ADRD, PD, and ALS. These loci support lysosomal/autophagic dysfunction (GAK/TMEM175, GRN, KANSL1), neuroinflammation/immunity (TSPOAP1), oxidative stress (GPX3, KANSL1), and the DNA damage response (NEK1) as transdiagnostic processes underlying multiple neurodegenerative disorders.

A human FLII gene variant alters sarcomeric actin thin filament length and predisposes to cardiomyopathy

To better understand the genetic basis of heart disease, we identified a variant in the Flightless-I homolog (FLII) gene that generates a R1243H missense change and predisposes to cardiac remodeling across multiple previous human genome-wide association studies (GWAS). Since this gene is of unknown function in the mammalian heart we generated gain- and loss-of-function genetically altered mice, as well as knock-in mice with the syntenic R1245H amino acid substitution, which showed that Flii protein binds the sarcomeric actin thin filament and influences its length. Deletion of Flii from the heart, or mice with the R1245H amino acid substitution, show cardiomyopathy due to shortening of the actin thin filaments. Mechanistically, Flii is a known actin binding protein that we show associates with tropomodulin-1 (TMOD1) to regulate sarcomere thin filament length. Indeed, overexpression of leiomodin-2 in the heart, which lengthens the actin-containing thin filaments, partially rescued disease due to heart-specific deletion of Flii. Collectively, the identified FLII human variant likely increases cardiomyopathy risk through an alteration in sarcomere structure and associated contractile dynamics, like other sarcomere gene-based familial cardiomyopathies.

The emergence of genotypic divergence and future precision medicine applications

Genotypic divergence is a term adapted from population genetics and intimately linked to evolution. We use divergence here to emphasize the differences that set individuals apart in any cohort. The history of genetics is filled with descriptions of genotypic differences, but causal inference of interindividual biological variation has been scarce. We suggest that the practice of precision medicine requires a divergent approach, an approach dependent on the causal interpretation of previous convergent (and preliminary) knowledge in the field. This knowledge has relied on convergent descriptive syndromology (lumping), which has overemphasized a reductionistic gene determinism on the quest of seeking associations without causal understanding. Regulatory variants with small effect and somatic mutations are some of the modifying factors that lead to incomplete penetrance and intrafamilial variable expressivity often observed in apparently monogenic clinical disorders. A truly divergent approach to precision medicine requires splitting, that is, the consideration of different layers of genetic phenomena that interact causally in a nonlinear fashion. This chapter reviews convergences and divergences in genetics and genomics, aiming to discuss what can be causally understood to approximate the as-yet utopian lands of Precision Medicine for patients with neurodegenerative disorders.

Common and rare variant associations with latent traits underlying depression, bipolar disorder, and schizophrenia

Genetic studies in psychiatry have primarily focused on the effects of common genetic variants, but few have investigated the role of rare genetic variants, particularly for major depression. In order to explore the role of rare variants in the gap between estimates of single nucleotide polymorphism (SNP) heritability and twin study heritability, we examined the contribution of common and rare genetic variants to latent traits underlying psychiatric disorders using high-quality imputed genotype data from the UK Biobank. Using a pre-registered analysis, we used items from the UK Biobank Mental Health Questionnaire relevant to three psychiatric disorders: major depression (N = 134,463), bipolar disorder (N = 117,376) and schizophrenia (N = 130,013) and identified a general hierarchical factor for each that described participants' responses. We calculated participants' scores on these latent traits and conducted single-variant genetic association testing (MAF > 0.05%), gene-based burden testing and pathway association testing associations with these latent traits. We tested for enrichment of rare variants (MAF 0.05-1%) in genes that had been previously identified by common variant genome-wide association studies, and genes previously associated with Mendelian disorders having relevant symptoms. We found moderate genetic correlations between the latent traits in our study and case-control phenotypes in previous genome-wide association studies, and identified one common genetic variant (rs72657988, minor allele frequency = 8.23%, p = 1.01 × 10^-9) associated with the general factor of schizophrenia, but no other single variants, genes or pathways passed significance thresholds in this analysis, and we did not find enrichment in previously identified genes.

The missing link between genetic association and regulatory function

The genetic basis of most traits is highly polygenic and dominated by non-coding alleles. It is widely assumed that such alleles exert small regulatory effects on the expression of cis-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene-trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue this 'missing regulation.'

Scalable Functional Assays for the Interpretation of Human Genetic Variation

Scalable sequence-function studies have enabled the systematic analysis and cataloging of hundreds of thousands of coding and noncoding genetic variants in the human genome. This has improved clinical variant interpretation and provided insights into the molecular, biophysical, and cellular effects of genetic variants at an astonishing scale and resolution across the spectrum of allele frequencies. In this review, we explore current applications and prospects for the field and outline the principles underlying scalable functional assay design, with a focus on the study of single-nucleotide coding and noncoding variants.

Deconstructing a Syndrome: Genomic Insights Into PCOS Causal Mechanisms and Classification

Polycystic ovary syndrome (PCOS) is among the most common disorders in women of reproductive age, affecting up to 15% worldwide, depending on the diagnostic criteria. PCOS is characterized by a constellation of interrelated reproductive abnormalities, including disordered gonadotropin secretion, increased androgen production, chronic anovulation, and polycystic ovarian morphology. It is frequently associated with insulin resistance and obesity. These reproductive and metabolic derangements cause major morbidities across the lifespan, including anovulatory infertility and type 2 diabetes (T2D). Despite decades of investigative effort, the etiology of PCOS remains unknown. Familial clustering of PCOS cases has indicated a genetic contribution to PCOS. There are rare Mendelian forms of PCOS associated with extreme phenotypes, but PCOS typically follows a non-Mendelian pattern of inheritance consistent with a complex genetic architecture, analogous to T2D and obesity, that reflects the interaction of susceptibility genes and environmental factors. Genomic studies of PCOS have provided important insights into disease pathways and have indicated that current diagnostic criteria do not capture underlying differences in biology associated with different forms of PCOS. We provide a state-of-the-science review of genetic analyses of PCOS, including an overview of genomic methodologies aimed at a general audience of non-geneticists and clinicians. Applications in PCOS will be discussed, including strengths and limitations of each study. The contributions of environmental factors, including developmental origins, will be reviewed. Insights into the pathogenesis and genetic architecture of PCOS will be summarized. Future directions for PCOS genetic studies will be outlined.

Familial Clustering and Genetic Analysis of Severe Thumb Carpometacarpal Joint Osteoarthritis in a Large Statewide Cohort

PURPOSE

Our goals were to identify individuals who required surgery for thumb carpometacarpal (CMC) joint osteoarthritis (OA), determine if CMC joint OA clusters in families, define the magnitude of familial risk of CMC joint OA, identify risk factors associated with CMC joint OA, and identify rare genetic variants that segregate with familial CMC joint OA.

METHODS

We searched the Utah Population Database to identify a cohort of CMC joint OA patients who required surgery. Affected individuals were mapped to pedigrees to identify high-risk families with excess clustering of CMC joint OA. Cox regression models were used to calculate familial risk of CMC joint OA in related individuals. Risk factors were evaluated using logistic regression models. Whole exome sequencing was used to identify rare coding variants associated with familial CMC joint OA.

RESULTS

We identified 550 pedigrees with excess clustering of severe CMC joint OA. The relative risk of CMC joint OA requiring surgical treatment was elevated significantly in first- and third-degree relatives of affected individuals, and significant associations with advanced age, female sex, obesity, and tobacco use were observed. We discovered candidate genes that dominantly segregate with severe CMC joint OA in 4 independent families, including a rare variant in Chondroitin Sulfate Synthase 3 (CHSY3).

CONCLUSIONS

Familial clustering of severe CMC joint OA was observed in a statewide population. Our data indicate that genetic and environmental factors contribute to the disease process, further highlighting the multifactorial nature of the disease. Genomic analyses suggest distinct biological processes are involved in CMC joint OA pathogenesis.

CLINICAL RELEVANCE

Awareness of associated comorbidities may guide the diagnosis of CMC joint OA in at-risk populations and help identify individuals who may not do well with nonoperative treatment. Further pursuit of the genes associated with severe CMC joint OA may lead to assays for detection of early stages of disease and have therapeutic potential.

NOD/RIPK2 signalling pathway contributes to osteoarthritis susceptibility

OBJECTIVES

How inflammatory signalling contributes to osteoarthritis (OA) susceptibility is undetermined. An allele encoding a hyperactive form of the Receptor Interacting Protein Kinase 2 (RIPK2) proinflammatory signalling intermediate has been associated with familial OA. To test whether altered nucleotide-binding oligomerisation domain (NOD)/RIPK2 pathway activity causes heightened OA susceptibility, we investigated whether variants affecting additional pathway components are associated with familial OA. To determine whether the Ripk2104Asp disease allele is sufficient to account for the familial phenotype, we determined the effect of the allele on mice.

METHODS

Genomic analysis of 150 independent families with dominant inheritance of OA affecting diverse joints was used to identify coding variants that segregated strictly with occurrence of OA. Genome editing was used to introduce the OA-associated RIPK2 (p.Asn104Asp) allele into the genome of inbred mice. The consequences of the Ripk2104Asp disease allele on physiology and OA susceptibility in mice were measured by histology, immunohistochemistry, serum cytokine levels and gene expression.

RESULTS

We identified six novel variants affecting components of the NOD/RIPK2 inflammatory signalling pathway that are associated with familial OA affecting the hand, shoulder or foot. The Ripk2104Asp allele acts dominantly to alter basal physiology and response to trauma in the mouse knee. Whereas the knees of uninjured Ripk2Asp104 mice appear normal histologically, the joints exhibit a set of marked gene expression changes reminiscent of overt OA. Although the Ripk2104Asp mice lack evidence of chronically elevated systemic inflammation, they do exhibit significantly increased susceptibility to post-traumatic OA (PTOA).

CONCLUSIONS

Two types of data support the hypothesis that altered NOD/RIPK2 signalling confers susceptibility to OA.

Polygenic enrichment distinguishes disease associations of individual cells in single-cell RNA-seq data

Single-cell RNA sequencing (scRNA-seq) provides unique insights into the pathology and cellular origin of disease. We introduce single-cell disease relevance score (scDRS), an approach that links scRNA-seq with polygenic disease risk at single-cell resolution, independent of annotated cell types. scDRS identifies cells exhibiting excess expression across disease-associated genes implicated by genome-wide association studies (GWASs). We applied scDRS to 74 diseases/traits and 1.3 million single-cell gene-expression profiles across 31 tissues/organs. Cell-type-level results broadly recapitulated known cell-type-disease associations. Individual-cell-level results identified subpopulations of disease-associated cells not captured by existing cell-type labels, including T cell subpopulations associated with inflammatory bowel disease, partially characterized by their effector-like states; neuron subpopulations associated with schizophrenia, partially characterized by their spatial locations; and hepatocyte subpopulations associated with triglyceride levels, partially characterized by their higher ploidy levels. Genes whose expression was correlated with the scDRS score across cells (reflecting coexpression with GWAS disease-associated genes) were strongly enriched for gold-standard drug target and Mendelian disease genes.

Mendel paved the path toward understanding genetic diseases

Background

July 20th, 2022, marks the 200th anniversary of the “Father of Genetics,” Gregor Mendel’s birth. His experiments with pea plants established the fundamental principles of genetic inheritance.

Main text

In this article, a succinct review of literature is hereby done to answer two key questions: (1) How Mendel’s principles of genetic inheritance helped us understand Mendelian disorders? and (2) How the study of Mendelian disorders can help us understand complex diseases?

Conclusion

This literature review concludes that continued effort to understand the genetic basis of Mendelian disorders will improve our understanding and treatment strategies for complex diseases.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

Polygenic risk scores of endo-phenotypes identify the effect of genetic background in congenital heart disease

Congenital heart disease (CHD) is a rare structural defect that occurs in ∼1% of live births. Studies on CHD genetic architecture have identified pathogenic single-gene mutations in less than 30% of cases. Single-gene mutations often show incomplete penetrance and variable expressivity. Therefore, we hypothesize that genetic background may play a role in modulating disease expression. Polygenic risk scores (PRSs) aggregate effects of common genetic variants to investigate whether, cumulatively, these variants are associated with disease penetrance or severity. However, the major limitations in this field have been in generating sufficient sample sizes for these studies. Here we used CHD-phenotype matched genome-wide association study (GWAS) summary statistics from the UK Biobank (UKBB) as our base study and whole-genome sequencing data from the CHD cohort (n1 = 711 trios, n2 = 362 European trios) of the Gabriella Miller Kids First dataset as our target study to develop PRSs for CHD. PRSs estimated using a GWAS for heart valve problems and heart murmur explain 2.5% of the variance in case-control status of CHD (all SNVs, p = 7.90 × 10-3; fetal cardiac SNVs, p = 8.00 × 10-3) and 1.8% of the variance in severity of CHD (fetal cardiac SNVs, p = 6.20 × 10-3; all SNVs, p = 0.015). These results show that common variants captured in CHD phenotype-matched GWASs have a modest but significant contribution to phenotypic expression of CHD. Further exploration of the cumulative effect of common variants is necessary for understanding the complex genetic etiology of CHD and other rare diseases.

SNP-to-gene linking strategies reveal contributions of enhancer-related and candidate master-regulator genes to autoimmune disease

We assess contributions to autoimmune disease of genes whose regulation is driven by enhancer regions (enhancer-related) and genes that regulate other genes in trans (candidate master-regulator). We link these genes to SNPs using several SNP-to-gene (S2G) strategies and apply heritability analyses to draw three conclusions about 11 autoimmune/blood-related diseases/traits. First, several characterizations of enhancer-related genes using functional genomics data are informative for autoimmune disease heritability after conditioning on a broad set of regulatory annotations. Second, candidate master-regulator genes defined using trans-eQTL in blood are also conditionally informative for autoimmune disease heritability. Third, integrating enhancer-related and master-regulator gene sets with protein-protein interaction (PPI) network information magnified their disease signal. The resulting PPI-enhancer gene score produced >2-fold stronger heritability signal and >2-fold stronger enrichment for drug targets, compared with the recently proposed enhancer domain score. In each case, functionally informed S2G strategies produced 4.1- to 13-fold stronger disease signals than conventional window-based strategies.

Translational impact of omics studies in alopecia areata: recent advances and future perspectives

INTRODUCTION

Alopecia areata (AA) is a non-scarring, hair loss disorder and a common autoimmune-mediated disease with an estimated lifetime risk of about 2%. To date, the treatment of AA is mainly based on suppression or stimulation of the immune response. Genomics and transcriptomics studies generated important insights into the underlying pathophysiology, enabled discovery of molecular disease signatures, which were used in some of the recent clinical trials to monitor drug response and substantiated the consideration of new therapeutic modalities for the treatment of AA such as abatacept, dupilumab, ustekinumab and Janus Kinase (JAK) inhibitors.

AREAS COVERED

In this review, genomics and transcriptomics studies in AA are discussed in detail with particular emphasis on their past and prospective translational impacts. Microbiome studies are also briefly introduced.

EXPERT OPINION

The generation of large datasets using the new high-throughput technologies has revolutionized medical research and AA has also benefited from the wave of omics studies. However, the limitations associated with JAK inhibitors and clinical heterogeneity in AA patients underscore the necessity for continuing omics research in AA for discovery of novel therapeutic modalities and development of clinical tools for precision medicine.

Multivariate GWAS of psychiatric disorders and their cardinal symptoms reveal two dimensions of cross-cutting genetic liabilities

Understanding which biological pathways are specific versus general across diagnostic categories and levels of symptom severity is critical to improving nosology and treatment of psychopathology. Here, we combine transdiagnostic and dimensional approaches to genetic discovery for the first time, conducting a novel multivariate genome-wide association study of eight psychiatric symptoms and disorders broadly related to mood disturbance and psychosis. We identify two transdiagnostic genetic liabilities that distinguish between common forms of psychopathology versus rarer forms of serious mental illness. Biological annotation revealed divergent genetic architectures that differentially implicated prenatal neurodevelopment and neuronal function and regulation. These findings inform psychiatric nosology and biological models of psychopathology, as they suggest that the severity of mood and psychotic symptoms present in serious mental illness may reflect a difference in kind rather than merely in degree.

Combining SNP-to-gene linking strategies to identify disease genes and assess disease omnigenicity

Disease-associated single-nucleotide polymorphisms (SNPs) generally do not implicate target genes, as most disease SNPs are regulatory. Many SNP-to-gene (S2G) linking strategies have been developed to link regulatory SNPs to the genes that they regulate in cis. Here, we developed a heritability-based framework for evaluating and combining different S2G strategies to optimize their informativeness for common disease risk. Our optimal combined S2G strategy (cS2G) included seven constituent S2G strategies and achieved a precision of 0.75 and a recall of 0.33, more than doubling the recall of any individual strategy. We applied cS2G to fine-mapping results for 49 UK Biobank diseases/traits to predict 5,095 causal SNP-gene-disease triplets (with S2G-derived functional interpretation) with high confidence. We further applied cS2G to provide an empirical assessment of disease omnigenicity; we determined that the top 1% of genes explained roughly half of the SNP heritability linked to all genes and that gene-level architectures vary with variant allele frequency.

Partitioning gene-level contributions to complex-trait heritability by allele frequency identifies disease-relevant genes

Recent works have shown that SNP heritability-which is dominated by low-effect common variants-may not be the most relevant quantity for localizing high-effect/critical disease genes. Here, we introduce methods to estimate the proportion of phenotypic variance explained by a given assignment of SNPs to a single gene ("gene-level heritability"). We partition gene-level heritability by minor allele frequency (MAF) to find genes whose gene-level heritability is explained exclusively by "low-frequency/rare" variants (0.5% ≤ MAF < 1%). Applying our method to ∼16K protein-coding genes and 25 quantitative traits in the UK Biobank (N = 290K "White British"), we find that, on average across traits, ∼2.5% of nonzero-heritability genes have a rare-variant component and only ∼0.8% (327 gene-trait pairs) have heritability exclusively from rare variants. Of these 327 gene-trait pairs, 114 (35%) were not detected by existing gene-level association testing methods. The additional genes we identify are significantly enriched for known disease genes, and we find several examples of genes that have been previously implicated in phenotypically related Mendelian disorders. Notably, the rare-variant component of gene-level heritability exhibits trends different from those of common-variant gene-level heritability. For example, while total gene-level heritability increases with gene length, the rare-variant component is significantly larger among shorter genes; the cumulative distributions of gene-level heritability also vary across traits and reveal differences in the relative contributions of rare/common variants to overall gene-level polygenicity. While nonzero gene-level heritability does not imply causality, if interpreted in the correct context, gene-level heritability can reveal useful insights into complex-trait genetic architecture.

Partitioning gene-mediated disease heritability without eQTLs

Unknown SNP-to-gene regulatory architecture complicates efforts to link noncoding GWAS associations with genes implicated by sequencing or functional studies. eQTLs are often used to link SNPs to genes, but expression in bulk tissue explains a small fraction of disease heritability. A simple but successful approach has been to link SNPs with nearby genes via base pair windows, but genes may often be regulated by SNPs outside their window. We propose the abstract mediation model (AMM) to estimate (1) the fraction of heritability mediated by the closest or kth-closest gene to each SNP and (2) the mediated heritability enrichment of a gene set (e.g., genes with rare-variant associations). AMM jointly estimates these quantities by matching the decay in SNP enrichment with distance from genes in the gene set. Across 47 complex traits and diseases, we estimate that the closest gene to each SNP mediates 27% (SE: 6%) of heritability and that a substantial fraction is mediated by genes outside the ten closest. Mendelian disease genes are strongly enriched for common-variant heritability; for example, just 21 dyslipidemia genes mediate 25% of LDL heritability (211× enrichment, p = 0.01). Among brain-related traits, genes involved in neurodevelopmental disorders are only about 4× enriched, but gene expression patterns are highly informative, as they have detectable differences in per-gene heritability even among weakly brain-expressed genes.

Priority index: database of genetic targets in immune-mediated disease

We describe a comprehensive and unique database 'Priority index' (Pi; http://pi.well.ox.ac.uk) of prioritized genes encoding potential therapeutic targets that encompasses all major immune-mediated diseases. We provide targets at the gene level, each receiving a 5-star rating supported by: genomic evidence arising from disease genome-wide associations and functional immunogenomics, annotation evidence using ontologies restricted to genes with genomic evidence, and network evidence from protein interactions. Target genes often act together in related molecular pathways. The underlying Pi approach is unique in identifying a network of highly rated genes that mediate pathway crosstalk. In the Pi website, disease-centric pages are specially designed to enable the users to browse a complete list of prioritized genes and also a manageable list of nodal genes at the pathway crosstalk level; both switchable by clicks. Moreover, target genes are cross-referenced and supported using additional information, particularly regarding tractability, including druggable pockets viewed in 3D within protein structures. Target genes highly rated across diseases suggest drug repurposing opportunity, while genes in a particular disease reveal disease-specific targeting potential. To facilitate the ease of such utility, cross-disease comparisons involving multiple diseases are also supported. This facility, together with the faceted search, enhances integrative mining of the Pi resource to accelerate early-stage therapeutic target identification and validation leveraging human genetics.

Bench Research Informed by GWAS Results

Scientifically interesting as well as practically important phenotypes often belong to the realm of complex traits. To the extent that these traits are hereditary, they are usually 'highly polygenic'. The study of such traits presents a challenge for researchers, as the complex genetic architecture of such traits makes it nearly impossible to utilise many of the usual methods of reverse genetics, which often focus on specific genes. In recent years, thousands of genome-wide association studies (GWAS) were undertaken to explore the relationships between complex traits and a large number of genetic factors, most of which are characterised by tiny effects. In this review, we aim to familiarise 'wet biologists' with approaches for the interpretation of GWAS results, to clarify some issues that may seem counterintuitive and to assess the possibility of using GWAS results in experiments on various complex traits.

Integration of genetic, transcriptomic, and clinical data provides insight into 16p11.2 and 22q11.2 CNV genes

BACKGROUND

Deletions and duplications of the multigenic 16p11.2 and 22q11.2 copy number variant (CNV) regions are associated with brain-related disorders including schizophrenia, intellectual disability, obesity, bipolar disorder, and autism spectrum disorder (ASD). The contribution of individual CNV genes to each of these identified phenotypes is unknown, as well as the contribution of these CNV genes to other potentially subtler health implications for carriers. Hypothesizing that DNA copy number exerts most effects via impacts on RNA expression, we attempted a novel in silico fine-mapping approach in non-CNV carriers using both GWAS and biobank data.

METHODS

We first asked whether gene expression level in any individual gene in the CNV region alters risk for a known CNV-associated behavioral phenotype(s). Using transcriptomic imputation, we performed association testing for CNV genes within large genotyped cohorts for schizophrenia, IQ, BMI, bipolar disorder, and ASD. Second, we used a biobank containing electronic health data to compare the medical phenome of CNV carriers to controls within 700,000 individuals in order to investigate the full spectrum of health effects of the CNVs. Third, we used genotypes for over 48,000 individuals within the biobank to perform phenome-wide association studies between imputed expressions of individual 16p11.2 and 22q11.2 genes and over 1500 health traits.

RESULTS

Using large genotyped cohorts, we found individual genes within 16p11.2 associated with schizophrenia (TMEM219, INO80E, YPEL3), BMI (TMEM219, SPN, TAOK2, INO80E), and IQ (SPN), using conditional analysis to identify upregulation of INO80E as the driver of schizophrenia, and downregulation of SPN and INO80E as increasing BMI. We identified both novel and previously observed over-represented traits within the electronic health records of 16p11.2 and 22q11.2 CNV carriers. In the phenome-wide association study, we found seventeen significant gene-trait pairs, including psychosis (NPIPB11, SLX1B) and mood disorders (SCARF2), and overall enrichment of mental traits.

CONCLUSIONS

Our results demonstrate how integration of genetic and clinical data aids in understanding CNV gene function and implicates pleiotropy and multigenicity in CNV biology.

Estimation of regional polygenicity from GWAS provides insights into the genetic architecture of complex traits

The number of variants that have a non-zero effect on a trait (i.e. polygenicity) is a fundamental parameter in the study of the genetic architecture of a complex trait. Although many previous studies have investigated polygenicity at a genome-wide scale, a detailed understanding of how polygenicity varies across genomic regions is currently lacking. In this work, we propose an accurate and scalable statistical framework to estimate regional polygenicity for a complex trait. We show that our approach yields approximately unbiased estimates of regional polygenicity in simulations across a wide-range of various genetic architectures. We then partition the polygenicity of anthropometric and blood pressure traits across 6-Mb genomic regions (N = 290K, UK Biobank) and observe that all analyzed traits are highly polygenic: over one-third of regions harbor at least one causal variant for each of the traits analyzed. Additionally, we observe wide variation in regional polygenicity: on average across all traits, 48.9% of regions contain at least 5 causal SNPs, 5.44% of regions contain at least 50 causal SNPs. Finally, we find that heritability is proportional to polygenicity at the regional level, which is consistent with the hypothesis that heritability enrichments are largely driven by the variation in the number of causal SNPs.

The proportion of SNPs with nonzero effects on a trait, or polygenicity, is a key quantity used to describe the genetic architecture of a complex trait. Furthermore, identifying the specific genomic regions that contribute to trait variation requires an understanding of how the number of causal SNPs varies across regions of the genome (regional polygenicity). In this work, we propose a statistical framework to estimate regional polygenicity for a complex trait using marginal effect sizes from GWAS and LD information. We demonstrate in simulation and empirical data that our approach accurately and efficiently estimates regional polygenicity. We find that SNP-heritability is proportional to polygenicity both on the genome-wide and regional scale, suggesting that the observed differences in heritability across traits stem from differences in the underlying number of causal SNPs.

Finding commonalities in rare diseases through the undiagnosed diseases network

Objective

When studying any specific rare disease, heterogeneity and scarcity of affected individuals has historically hindered investigators from discerning on what to focus to understand and diagnose a disease. New nongenomic methodologies must be developed that identify similarities in seemingly dissimilar conditions.

Materials and Methods

This observational study analyzes 1042 patients from the Undiagnosed Diseases Network (2015-2019), a multicenter, nationwide research study using phenotypic data annotated by specialized staff using Human Phenotype Ontology terms. We used Louvain community detection to cluster patients linked by Jaccard pairwise similarity and 2 support vector classifier to assign new cases. We further validated the clusters’ most representative comorbidities using a national claims database (67 million patients).

Results

Patients were divided into 2 groups: those with symptom onset before 18 years of age (n = 810) and at 18 years of age or older (n = 232) (average symptom onset age: 10 [interquartile range, 0-14] years). For 810 pediatric patients, we identified 4 statistically significant clusters. Two clusters were characterized by growth disorders, and developmental delay enriched for hypotonia presented a higher likelihood of diagnosis. Support vector classifier showed 0.89 balanced accuracy (0.83 for Human Phenotype Ontology terms only) on test data.

Discussions

To set the framework for future discovery, we chose as our endpoint the successful grouping of patients by phenotypic similarity and provide a classification tool to assign new patients to those clusters.

Conclusion

This study shows that despite the scarcity and heterogeneity of patients, we can still find commonalities that can potentially be harnessed to uncover new insights and targets for therapy.

Genetic Influences on Disease Subtypes

Disease classification, or nosology, was historically driven by careful examination of clinical features of patients. As technologies to measure and understand human phenotypes advanced, so too did classifications of disease, and the advent of genetic data has led to a surge in genetic subtyping in the past decades. Although the fundamental process of refining disease definitions and subtypes is shared across diverse fields, each field is driven by its own goals and technological expertise, leading to inconsistent and conflicting definitions of disease subtypes. Here, we review several classical and recent subtypes and subtyping approaches and provide concrete definitions to delineate subtypes. In particular, we focus on subtypes with distinct causal disease biology, which are of primary interest to scientists, and subtypes with pragmatic medical benefits, which are of primary interest to physicians. We propose genetic heterogeneity as a gold standard for establishing biologically distinct subtypes of complex polygenic disease. We focus especially on methods to find and validate genetic subtypes, emphasizing common pitfalls and how to avoid them.

Modeling regulatory network topology improves genome-wide analyses of complex human traits

Genome-wide association studies (GWAS) have cataloged many significant associations between genetic variants and complex traits. However, most of these findings have unclear biological significance, because they often have small effects and occur in non-coding regions. Integration of GWAS with gene regulatory networks addresses both issues by aggregating weak genetic signals within regulatory programs. Here we develop a Bayesian framework that integrates GWAS summary statistics with regulatory networks to infer genetic enrichments and associations simultaneously. Our method improves upon existing approaches by explicitly modeling network topology to assess enrichments, and by automatically leveraging enrichments to identify associations. Applying this method to 18 human traits and 38 regulatory networks shows that genetic signals of complex traits are often enriched in interconnections specific to trait-relevant cell types or tissues. Prioritizing variants within enriched networks identifies known and previously undescribed trait-associated genes revealing biological and therapeutic insights.

Functional dynamic genetic effects on gene regulation are specific to particular cell types and environmental conditions

Genetic effects on gene expression and splicing can be modulated by cellular and environmental factors; yet interactions between genotypes, cell type, and treatment have not been comprehensively studied together. We used an induced pluripotent stem cell system to study multiple cell types derived from the same individuals and exposed them to a large panel of treatments. Cellular responses involved different genes and pathways for gene expression and splicing and were highly variable across contexts. For thousands of genes, we identified variable allelic expression across contexts and characterized different types of gene-environment interactions, many of which are associated with complex traits. Promoter functional and evolutionary features distinguished genes with elevated allelic imbalance mean and variance. On average, half of the genes with dynamic regulatory interactions were missed by large eQTL mapping studies, indicating the importance of exploring multiple treatments to reveal previously unrecognized regulatory loci that may be important for disease.

MendelVar: gene prioritization at GWAS loci using phenotypic enrichment of Mendelian disease genes

Motivation

Gene prioritization at human GWAS loci is challenging due to linkage-disequilibrium and long-range gene regulatory mechanisms. However, identifying the causal gene is crucial to enable identification of potential drug targets and better understanding of molecular mechanisms. Mapping GWAS traits to known phenotypically relevant Mendelian disease genes near a locus is a promising approach to gene prioritization.

Results

We present MendelVar, a comprehensive tool that integrates knowledge from four databases on Mendelian disease genes with enrichment testing for a range of associated functional annotations such as Human Phenotype Ontology, Disease Ontology and variants from ClinVar. This open web-based platform enables users to strengthen the case for causal importance of phenotypically matched candidate genes at GWAS loci. We demonstrate the use of MendelVar in post-GWAS gene annotation for type 1 diabetes, type 2 diabetes, blood lipids and atopic dermatitis.

Availability and implementation

MendelVar is freely available at https://mendelvar.mrcieu.ac.uk

Supplementary information

Supplementary data are available at Bioinformatics online.

The Impact of Environmental Factors on Monogenic Mendelian Diseases

Environmental factors and gene-environment interactions modify the variable expressivity, progression, severity, and onset of some classic (monogenic) Mendelian-inherited genetic diseases. Cystic fibrosis, Huntington disease, Parkinson's disease, and sickle cell disease are examples of well-known Mendelian disorders that are influenced by exogenous exposures. Environmental factors may act by direct or indirect mechanisms to modify disease severity, timing, and presentation, including through epigenomic influences, protein misfolding, miRNA alterations, transporter activity, and mitochondrial effects. Because pathological features of early-onset Mendelian diseases can mimic later onset complex diseases, we propose that studies of environmental exposure vulnerabilities using monogenic model systems of rare Mendelian diseases have high potential to provide insight into complex disease phenotypes arising from multi-genetic/multi-toxicant interactions. Mendelian disorders can be modeled by homologous mutations in animal model systems with strong recapitulation of human disease etiology and natural history, providing an important advantage for study of these diseases. Monogenic high penetrant mutations are ideal for toxicant challenge studies with a wide variety of environmental stressors, because background genetic variability may be less able to alter the relatively strong phenotype driving disease-causing mutations. These models promote mechanistic understandings of gene-environment interactions and biological pathways relevant to both Mendelian and related sporadic complex disease outcomes by creating a sensitized background for relevant environmental risk factors. Additionally, rare disease communities are motivated research participants, creating the potential of strong research allies among rare Mendelian disease advocacy groups and disease registries and providing a variety of translational opportunities that are under-utilized in genetic or environmental health science.

Genomics of hypertension: the road to precision medicine

The known genetic architecture of blood pressure now comprises >30 genes, with rare variants resulting in monogenic forms of hypertension or hypotension and >1,477 common single-nucleotide polymorphisms (SNPs) being associated with the blood pressure phenotype. Monogenic blood pressure syndromes predominantly involve the renin-angiotensin-aldosterone system and the adrenal glucocorticoid pathway, with a smaller fraction caused by neuroendocrine tumours of the sympathetic and parasympathetic nervous systems. The SNPs identified in genome-wide association studies (GWAS) as being associated with the blood pressure phenotype explain only approximately 27% of the 30-50% estimated heritability of blood pressure, and the effect of each SNP on the blood pressure phenotype is small. A paucity of SNPs from GWAS are mapped to known genes causing monogenic blood pressure syndromes. For example, a GWAS signal mapped to the gene encoding uromodulin has been shown to affect blood pressure by influencing sodium homeostasis, and the effects of another GWAS signal were mediated by endothelin. However, the majority of blood pressure-associated SNPs show pleiotropic associations. Unravelling these associations can potentially help us to understand the underlying biological pathways. In this Review, we appraise the current knowledge of blood pressure genomics, explore the causal pathways for hypertension identified in Mendelian randomization studies and highlight the opportunities for drug repurposing and pharmacogenomics for the treatment of hypertension.

Novel Linkage Peaks Discovered for Diabetic Nephropathy in Individuals With Type 1 Diabetes

Genome-wide association studies (GWAS) and linkage studies have had limited success in identifying genome-wide significantly linked regions or risk loci for diabetic nephropathy (DN) in individuals with type 1 diabetes (T1D). As GWAS cohorts have grown, they have also included more documented and undocumented familial relationships. Here we computationally inferred and manually curated pedigrees in a study cohort of >6,000 individuals with T1D and their relatives without diabetes. We performed a linkage study for 177 pedigrees consisting of 452 individuals with T1D and their relatives using a genome-wide genotyping array with >300,000 single nucleotide polymorphisms and PSEUDOMARKER software. Analysis resulted in genome-wide significant linkage peaks on eight chromosomal regions from five chromosomes (logarithm of odds score >3.3). The highest peak was localized at the HLA region on chromosome 6p, but whether the peak originated from T1D or DN remained ambiguous. Of other significant peaks, the chromosome 4p22 region was localized on top of ARHGAP24, a gene associated with focal segmental glomerulosclerosis, suggesting this gene may play a role in DN as well. Furthermore, rare variants have been associated with DN and chronic kidney disease near the 4q25 peak, localized on top of CCSER1.

Transcriptomic Insight Into the Polygenic Mechanisms Underlying Psychiatric Disorders

Over the past decade, large-scale genetic studies have successfully identified hundreds of genetic variants robustly associated with risk for psychiatric disorders. However, mechanistic insight and clinical translation continue to lag the pace of risk variant identification, hindered by the sheer number of targets and their predominant noncoding localization, as well as pervasive pleiotropy and incomplete penetrance. Successful next steps require identification of "causal" genetic variants and their proximal biological consequences; placing variants within biologically defined functional contexts, reflecting specific molecular pathways, cell types, circuits, and developmental windows; and characterizing the downstream, convergent neurobiological impact of polygenicity within an individual. Here, we discuss opportunities and challenges of high-throughput transcriptomic profiling in the human brain, and how transcriptomic approaches can help pinpoint mechanisms underlying genetic risk for psychiatric disorders at a scale necessary to tackle daunting levels of polygenicity. These include transcriptome-wide association studies for risk gene prioritization through integration of genome-wide association studies with expression quantitative trait loci. We outline transcriptomic results that inform our understanding of the brain-level molecular pathology of psychiatric disorders, including autism spectrum disorder, bipolar disorder, major depressive disorder, and schizophrenia. Finally, we discuss systems-level approaches for integration of distinct genetic, genomic, and phenotypic levels, including combining spatially resolved gene expression and human neuroimaging maps. Results highlight the importance of understanding gene expression (dys)regulation across human brain development as a major contributor to psychiatric disease pathogenesis, from common variants acting as expression quantitative trait loci to rare variants enriched for gene expression regulatory pathways.

Improving the informativeness of Mendelian disease-derived pathogenicity scores for common disease

Despite considerable progress on pathogenicity scores prioritizing variants for Mendelian disease, little is known about the utility of these scores for common disease. Here, we assess the informativeness of Mendelian disease-derived pathogenicity scores for common disease and improve upon existing scores. We first apply stratified linkage disequilibrium (LD) score regression to evaluate published pathogenicity scores across 41 common diseases and complex traits (average N = 320K). Several of the resulting annotations are informative for common disease, even after conditioning on a broad set of functional annotations. We then improve upon published pathogenicity scores by developing AnnotBoost, a machine learning framework to impute and denoise pathogenicity scores using a broad set of functional annotations. AnnotBoost substantially increases the informativeness for common disease of both previously uninformative and previously informative pathogenicity scores, implying that Mendelian and common disease variants share similar properties. The boosted scores also produce improvements in heritability model fit and in classifying disease-associated, fine-mapped SNPs. Our boosted scores may improve fine-mapping and candidate gene discovery for common disease.

A phenome-wide association study of 26 mendelian genes reveals phenotypic expressivity of common and rare variants within the general population

The clinical evaluation of a genetic syndrome relies upon recognition of a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. However, individuals displaying a single phenotype of a complex syndrome may not meet criteria for clinical diagnosis or genetic testing. Here, we present a phenome-wide association study (PheWAS) approach to systematically explore the phenotypic expressivity of common and rare alleles in genes associated with four well-described syndromic diseases (Alagille (AS), Marfan (MS), DiGeorge (DS), and Noonan (NS) syndromes) in the general population.

Using human phenotype ontology (HPO) terms, we systematically mapped 60 phenotypes related to AS, MS, DS and NS in 337,198 unrelated white British from the UK Biobank (UKBB) based on their hospital admission records, self-administrated questionnaires, and physiological measurements. We performed logistic regression adjusting for age, sex, and the first 5 genetic principal components, for each phenotype and each variant in the target genes (JAG1, NOTCH2 FBN1, PTPN1 and RAS-opathy genes, and genes in the 22q11.2 locus) and performed a gene burden test.

Overall, we observed multiple phenotype-genotype correlations, such as the association between variation in JAG1, FBN1, PTPN11 and SOS2 with diastolic and systolic blood pressure; and pleiotropy among multiple variants in syndromic genes. For example, rs11066309 in PTPN11 was significantly associated with a lower body mass index, an increased risk of hypothyroidism and a smaller size for gestational age, all in concordance with NS-related phenotypes. Similarly, rs589668 in FBN1 was associated with an increase in body height and blood pressure, and a reduced body fat percentage as observed in Marfan syndrome.

Our findings suggest that the spectrum of associations of common and rare variants in genes involved in syndromic diseases can be extended to individual phenotypes within the general population.

Standard medical evaluation of genetic syndromes relies upon recognizing a characteristic pattern of signs or symptoms to guide targeted genetic testing for confirmation of the diagnosis. This may lead to missing diagnoses in patients with silent or a low expressed form of the syndrome. Here we take advantage of a rich electronic health record, various phenotypic measurements, and genetic information in 337,198 unrelated white British from the UKBB, to study the relation between single syndromic disease phenotypes and genes related to syndromic disease. We show multiple phenotype-genotype associations in concordance with phenotypes variations found in syndromic diseases. For example, we show that a commonly found variant in FBN1 was associated with high standing/sitting height ratio and reduced body fat percentage as observed in individuals with Marfan syndrome. Our findings suggest that common and rare alleles in syndromic disease genes are causative of individual component phenotypes present in a general population; further research is needed to characterize the pleiotropic effect of alleles in syndromic genes in persons without the syndromic disease.

An Imperative Need for Further Genetic Studies of Alopecia Areata

Human genetic studies of diseases that are multifactorial and prevalent have generated a wealth of knowledge about the genetic architecture of chronic diseases. Generalizable attributes are shaping the development of models to explain how the human genome influences our health and can be leveraged to improve it. Importantly, both rare and common genetic variants contribute to disease risk and provide complementary information. Although initial genetic studies of alopecia areata have yielded insight with high clinical impact, there remains a number of important unanswered questions pertaining to disease biology and patient care that could be addressed by further genetic investigations.

Molecular pathways in patients with systemic lupus erythematosus revealed by gene-centred DNA sequencing

Objectives

Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE.

Methods

We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS).

Results

We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes.

Conclusions

Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.

A Single Nucleotide ADA Genetic Variant Is Associated to Central Inflammation and Clinical Presentation in MS: Implications for Cladribine Treatment

In multiple sclerosis (MS), activated T and B lymphocytes and microglial cells release various proinflammatory cytokines, promoting neuroinflammation and negatively affecting the course of the disease. The immune response homeostasis is crucially regulated by the activity of the enzyme adenosine deaminase (ADA), as evidenced in patients with genetic ADA deficiency and in those treated with cladribine tablets. We investigated in a group of patients with MS the associations of a single nucleotide polymorphism (SNP) of ADA gene with disease characteristics and cerebrospinal fluid (CSF) inflammation. The SNP rs244072 of the ADA gene was determined in 561 patients with MS. Disease characteristics were assessed at the time of diagnosis; furthermore, in 258 patients, proinflammatory and anti-inflammatory molecules were measured in the CSF. We found a significant association between rs244072 and both clinical characteristics and central inflammation. In C-carriers, significantly enhanced disability and increased CSF levels of TNF, IL-5 and RANTES was observed. In addition, lower CSF levels of the anti-inflammatory cytokine IL-10 were found. Finally, the presence of the C allele was associated with a tendency of increased lymphocyte count. In MS patients, ADA SNP rs244072 is associated with CSF inflammation and disability. The selective targeting of the ADA pathway through cladribine tablet therapy could be effective in MS by acting on a pathogenically relevant biological mechanism.