Analysis of rare genetic variation underlying cardiometabolic diseases and traits among 200,000 individuals in the UK Biobank

Prioritizing disease-related rare variants by integrating gene expression data

Rare variants, comprising a vast majority of human genetic variations, are likely to have more deleterious impact on human diseases compared to common variants. Here we present carrier statistic, a statistical framework to prioritize disease-related rare variants by integrating gene expression data. By quantifying the impact of rare variants on gene expression, carrier statistic can prioritize those rare variants that have large functional consequence in the diseased patients. Through simulation studies and analyzing real multi-omics dataset, we demonstrated that carrier statistic is applicable in studies with limited sample size (a few hundreds) and achieves substantially higher sensitivity than existing rare variants association methods. Application to Alzheimer's disease reveals 16 rare variants within 15 genes with extreme carrier statistics. We also found strong excess of rare variants among the top prioritized genes in diseased patients compared to that in healthy individuals. The carrier statistic method can be applied to various rare variant types and is adaptable to other omics data modalities, offering a powerful tool for investigating the molecular mechanisms underlying complex diseases.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Phenotypic overlap between rare disease patients and variant carriers in a large population cohort informs biological mechanisms

The biological mechanisms giving rise to the extreme symptoms exhibited by rare disease patients are complex, heterogenous, and difficult to discern. Understanding these mechanisms is critical for developing treatments that address the underlying causes of diseases rather than merely the presenting symptoms. Moreover, the same dysfunctional biological mechanisms implicated in rare recessive diseases may also lead to milder and potentially preventable symptoms in carriers in the general population. Seizures are a common, extreme phenotype that can result from diverse and often elusive biological pathways in patients with ultrarare or undiagnosed disorders. In this pilot study, we present an approach to understand the biological pathways leading to seizures in patients from the Undiagnosed Diseases Network (UDN) by analyzing aggregated genotype and phenotype data from the UK Biobank (UKB). Specifically, we look for enriched phenotypes across UKB participants who harbor rare variants in the same gene known or suspected to be causally implicated in a UDN patient's recessively manifesting disorder. Analyzing these milder but related associated phenotypes in UKB participants can provide insight into the disease-causing molecular mechanisms at play in the rare disease UDN patient. We present six vignettes of undiagnosed patients experiencing seizures as part of their recessive genetic condition, and we discuss the potential mechanisms underlying the spectrum of symptoms associated with UKB participants to the severe presentations exhibited by UDN patients. We find that in our set of rare disease patients, seizures may result from diverse, multi-step pathways that involve multiple body systems. Analyses of large-scale population cohorts such as the UKB can be a critical tool to further our understanding of rare diseases in general.

Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing

We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.

Atrial fibrillation: age at diagnosis, incident cardiovascular events, and mortality

BACKGROUND AND AIMS

Patients with atrial fibrillation (AF) are at increased risks of cardiovascular diseases and mortality, but risks according to age at diagnosis have not been reported. This study investigated age-specific risks of outcomes among patients with AF and the background population.

METHODS

This nationwide population-based cohort study included patients with AF and controls without outcomes by the application of exposure density matching on the basis of sex, year of birth, and index date. The absolute risks and hazard rates were stratified by age groups and assessed using competing risk survival analyses and Cox regression models, respectively. The expected differences in residual life years among participants were estimated.

RESULTS

The study included 216 579 AF patients from year 2000 to 2020 and 866 316 controls. The mean follow-up time was 7.9 years. Comparing AF patients with matched controls, the hazard ratios among individuals ≤50 years was 8.90 [95% confidence interval (CI), 7.17-11.0] for cardiomyopathy, 8.64 (95% CI, 7.74-9.64) for heart failure, 2.18 (95% CI, 1.89-2.52) for ischaemic stroke, and 2.74 (95% CI, 2.53-2.96) for mortality. The expected average loss of life years among individuals ≤50 years was 9.2 years (95% CI, 9.0-9.3) years. The estimates decreased with older age.

CONCLUSIONS

The findings show that earlier diagnosis of AF is associated with a higher hazard ratio of subsequent myocardial disease and shorter life expectancy. Further studies are needed to determine causality and whether AF could be used as a risk marker among particularly younger patients.

Large-scale whole-exome sequencing of neuropsychiatric diseases and traits in 350,770 adults

While numerous genomic loci have been identified for neuropsychiatric conditions, the contribution of protein-coding variants has yet to be determined. Here we conducted a large-scale whole-exome-sequencing study to interrogate the impact of protein-coding variants on 46 neuropsychiatric diseases and 23 traits in 350,770 adults from the UK Biobank. Twenty new genes were associated with neuropsychiatric diseases through coding variants, among which 16 genes had impacts on the longitudinal risks of diseases. Thirty new genes were associated with neuropsychiatric traits, with SYNGAP1 showing pleiotropic effects across cognitive function domains. Pairwise estimation of genetic correlations at the coding-variant level highlighted shared genetic associations among pairs of neurodegenerative diseases and mental disorders. Lastly, a comprehensive multi-omics analysis suggested that alterations in brain structures, blood proteins and inflammation potentially contribute to the gene-phenotype linkages. Overall, our findings characterized a compendium of protein-coding variants for future research on the biology and therapeutics of neuropsychiatric phenotypes.

Multi-ancestry Whole-exome Sequencing Study of Alcohol Use Disorder in Two Cohorts

Alcohol use disorder (AUD) is a leading cause of death and disability worldwide. There has been substantial progress in identifying genetic variants underlying AUD. However, there are few whole-exome sequencing (WES) studies of AUD. We analyzed WES of 4,530 samples from the Yale-Penn cohort and 469,835 samples from the UK Biobank (UKB). After quality control, 1,420 AUD cases and 619 controls of European ancestry (EUR) and 1,142 cases and 608 controls of African ancestry (AFR) from Yale-Penn were retained for subsequent analyses. WES data from 415,617 EUR samples (12,861 cases), 6,142 AFR samples (130 cases) and 4,607 South Asian (SAS) samples (130 cases) from UKB were also analyzed. Single-variant association analysis identified the well-known functional variant rs1229984 in ADH1B ( P =4.88×10 -31 ) and several other common variants in ADH1C . Gene-based tests identified ADH1B ( P =1.00×10 -31 ), ADH1C ( P =5.23×10 -7 ), CNST ( P =1.19×10 -6 ), and IFIT5 (3.74×10 -6 ). This study extends our understanding of the genetic basis of AUD.

Angiopoietin-like protein 8: a multifaceted protein instrumental in regulating triglyceride metabolism

PURPOSE OF REVIEW

The angiopoietin-like (ANGPTL) proteins ANGPTL3 and ANGPTL4 are critical lipoprotein lipase (LPL) inhibitors. This review discusses the unique ability of the insulin-responsive protein ANGPTL8 to regulate triglyceride (TG) metabolism by forming ANGPTL3/8 and ANGPTL4/8 complexes that control tissue-specific LPL activities.

RECENT FINDINGS

After feeding, ANGPTL4/8 acts locally in adipose tissue, has decreased LPL-inhibitory activity compared to ANGPTL4, and binds tissue plasminogen activator (tPA) and plasminogen to generate plasmin, which cleaves ANGPTL4/8 and other LPL inhibitors. This enables LPL to be fully active postprandially to promote efficient fatty acid (FA) uptake and minimize ectopic fat deposition. In contrast, liver-derived ANGPTL3/8 acts in an endocrine manner, has markedly increased LPL-inhibitory activity compared to ANGPTL3, and potently inhibits LPL in oxidative tissues to direct TG toward adipose tissue for storage. Circulating ANGPTL3/8 levels are strongly correlated with serum TG, and the ANGPTL3/8 LPL-inhibitory epitope is blocked by the TG-lowering protein apolipoprotein A5 (ApoA5).

SUMMARY

ANGPTL8 plays a crucial role in TG metabolism by forming ANGPTL3/8 and ANGPTL4/8 complexes that differentially modulate LPL activities in oxidative and adipose tissues respectively. Selective ANGPTL8 inhibition in the context of the ANGPTL3/8 complex has the potential to be a promising strategy for treating dyslipidemia.

Genetic associations of protein-coding variants in venous thromboembolism

Previous genetic studies of venous thromboembolism (VTE) have been largely limited to common variants, leaving the genetic determinants relatively incomplete. We performed an exome-wide association study of VTE among 14,723 cases and 334,315 controls. Fourteen known and four novel genes (SRSF6, PHPT1, CGN, and MAP3K2) were identified through protein-coding variants, with broad replication in the FinnGen cohort. Most genes we discovered exhibited the potential to predict future VTE events in longitudinal analysis. Notably, we provide evidence for the additive contribution of rare coding variants to known genome-wide polygenic risk in shaping VTE risk. The identified genes were enriched in pathways affecting coagulation and platelet activation, along with liver-specific expression. The pleiotropic effects of these genes indicated the potential involvement of coagulation factors, blood cell traits, liver function, and immunometabolic processes in VTE pathogenesis. In conclusion, our study unveils the valuable contribution of protein-coding variants in VTE etiology and sheds new light on its risk stratification.

Prioritizing disease-related rare variants by integrating gene expression data

Rare variants, comprising a vast majority of human genetic variations, are likely to have more deleterious impact on human diseases compared to common variants. Here we present carrier statistic, a statistical framework to prioritize disease-related rare variants by integrating gene expression data. By quantifying the impact of rare variants on gene expression, carrier statistic can prioritize those rare variants that have large functional consequence in the diseased patients. Through simulation studies and analyzing real multi-omics dataset, we demonstrated that carrier statistic is applicable in studies with limited sample size (a few hundreds) and achieves substantially higher sensitivity than existing rare variants association methods. Application to Alzheimer's disease reveals 16 rare variants within 15 genes with extreme carrier statistics. The carrier statistic method can be applied to various rare variant types and is adaptable to other omics data modalities, offering a powerful tool for investigating the molecular mechanisms underlying complex diseases.

Clonal haematopoiesis of indeterminate potential predicts incident cardiac arrhythmias

BACKGROUND AND AIMS

Clonal haematopoiesis of indeterminate potential (CHIP), the age-related expansion of blood cells with preleukemic mutations, is associated with atherosclerotic cardiovascular disease and heart failure. This study aimed to test the association of CHIP with new-onset arrhythmias.

METHODS

UK Biobank participants without prevalent arrhythmias were included. Co-primary study outcomes were supraventricular arrhythmias, bradyarrhythmias, and ventricular arrhythmias. Secondary outcomes were cardiac arrest, atrial fibrillation, and any arrhythmia. Associations of any CHIP [variant allele fraction (VAF) ≥ 2%], large CHIP (VAF ≥10%), and gene-specific CHIP subtypes with incident arrhythmias were evaluated using multivariable-adjusted Cox regression. Associations of CHIP with myocardial interstitial fibrosis [T1 measured using cardiac magnetic resonance (CMR)] were also tested.

RESULTS

This study included 410 702 participants [CHIP: n = 13 892 (3.4%); large CHIP: n = 9191 (2.2%)]. Any and large CHIP were associated with multi-variable-adjusted hazard ratios of 1.11 [95% confidence interval (CI) 1.04-1.18; P = .001] and 1.13 (95% CI 1.05-1.22; P = .001) for supraventricular arrhythmias, 1.09 (95% CI 1.01-1.19; P = .031) and 1.13 (95% CI 1.03-1.25; P = .011) for bradyarrhythmias, and 1.16 (95% CI, 1.00-1.34; P = .049) and 1.22 (95% CI 1.03-1.45; P = .021) for ventricular arrhythmias, respectively. Associations were independent of coronary artery disease and heart failure. Associations were also heterogeneous across arrhythmia subtypes and strongest for cardiac arrest. Gene-specific analyses revealed an increased risk of arrhythmias across driver genes other than DNMT3A. Large CHIP was associated with 1.31-fold odds (95% CI 1.07-1.59; P = .009) of being in the top quintile of myocardial fibrosis by CMR.

CONCLUSIONS

CHIP may represent a novel risk factor for incident arrhythmias, indicating a potential target for modulation towards arrhythmia prevention and treatment.

Exome sequencing identifies genes associated with sleep-related traits

Sleep is vital for human health and has a moderate heritability. Previous genome-wide association studies have limitations in capturing the role of rare genetic variants in sleep-related traits. Here we conducted a large-scale exome-wide association study of eight sleep-related traits (sleep duration, insomnia symptoms, chronotype, daytime sleepiness, daytime napping, ease of getting up in the morning, snoring and sleep apnoea) among 450,000 participants from UK Biobank. We identified 22 new genes associated with chronotype (ADGRL4, COL6A3, CLK4 and KRTAP3-3), daytime sleepiness (ST3GAL1 and ANKRD12), daytime napping (PLEKHM1, ANKRD12 and ZBTB21), snoring (WDR59) and sleep apnoea (13 genes). Notably, 20 of these genes were confirmed to be significantly associated with sleep disorders in the FinnGen cohort. Enrichment analysis revealed that these discovered genes were enriched in circadian rhythm and central nervous system neurons. Phenotypic association analysis showed that ANKRD12 was associated with cognition and inflammatory traits. Our results demonstrate the value of large-scale whole-exome analysis in understanding the genetic architecture of sleep-related traits and potential biological mechanisms.

Discovery and SAR Study of Boronic Acid-Based Selective PDE3B Inhibitors from a Novel DNA-Encoded Library

Human genetic evidence shows that PDE3B is associated with metabolic and dyslipidemia phenotypes. A number of PDE3 family selective inhibitors have been approved by the FDA for various indications; however, given the undesirable proarrhythmic effects in the heart, selectivity for PDE3B inhibition over closely related family members (such as PDE3A; 48% identity) is a critical consideration for development of PDE3B therapeutics. Selectivity for PDE3B over PDE3A may be achieved in a variety of ways, including properties intrinsic to the compound or tissue-selective targeting. The high (>95%) active site homology between PDE3A and B represents a massive obstacle for obtaining selectivity at the active site; however, utilization of libraries with high molecular diversity in high throughput screens may uncover selective chemical matter. Herein, we employed a DNA-encoded library screen to identify PDE3B-selective inhibitors and identified potent and selective boronic acid compounds bound at the active site.

Inborn Errors of Immunity Contribute to the Burden of Skin Disease and Create Opportunities for Improving the Practice of Dermatology

Opportunities to improve the clinical management of skin disease are being created by advances in genomic medicine. Large-scale sequencing increasingly challenges notions about single-gene disorders. It is now apparent that monogenic etiologies make appreciable contributions to the population burden of disease and that they are underrecognized in clinical practice. A genetic diagnosis informs on molecular pathology and may direct targeted treatments and tailored prevention strategies for patients and family members. It also generates knowledge about disease pathogenesis and management that is relevant to patients without rare pathogenic variants. Inborn errors of immunity are a large class of monogenic etiologies that have been well-studied and contribute to the population burden of inflammatory diseases. To further delineate the contributions of inborn errors of immunity to the pathogenesis of skin disease, we performed a set of analyses that identified 316 inborn errors of immunity associated with skin pathologies, including common skin diseases. These data suggest that clinical sequencing is underutilized in dermatology. We next use these data to derive a network that illuminates the molecular relationships of these disorders and suggests an underlying etiological organization to immune-mediated skin disease. Our results motivate the further development of a molecularly derived and data-driven reorganization of clinical diagnoses of skin disease.

The genetic architecture of the human hypothalamus and its involvement in neuropsychiatric behaviours and disorders

Despite its crucial role in the regulation of vital metabolic and neurological functions, the genetic architecture of the hypothalamus remains unknown. Here we conducted multivariate genome-wide association studies (GWAS) using hypothalamic imaging data from 32,956 individuals to uncover the genetic underpinnings of the hypothalamus and its involvement in neuropsychiatric traits. There were 23 significant loci associated with the whole hypothalamus and its subunits, with functional enrichment for genes involved in intracellular trafficking systems and metabolic processes of steroid-related compounds. The hypothalamus exhibited substantial genetic associations with limbic system structures and neuropsychiatric traits including chronotype, risky behaviour, cognition, satiety and sympathetic-parasympathetic activity. The strongest signal in the primary GWAS, the ADAMTS8 locus, was replicated in three independent datasets (N = 1,685-4,321) and was strengthened after meta-analysis. Exome-wide association analyses added evidence to the association for ADAMTS8, and Mendelian randomization showed lower ADAMTS8 expression with larger hypothalamic volumes. The current study advances our understanding of complex structure-function relationships of the hypothalamus and provides insights into the molecular mechanisms that underlie hypothalamic formation.

Cardiometabolic health across menopausal years is linked to white matter hyperintensities up to a decade later

Introduction

The menopause transition is associated with several cardiometabolic risk factors. Poor cardiometabolic health is further linked to microvascular brain lesions, which can be detected as white matter hyperintensities (WMHs) using T2-FLAIR magnetic resonance imaging (MRI) scans. Females show higher risk for WMHs post-menopause, but it remains unclear whether changes in cardiometabolic risk factors underlie menopause-related increase in brain pathology.

Methods

In this study, we assessed whether cross-sectional measures of cardiometabolic health, including body mass index (BMI) and waist-to-hip ratio (WHR), blood lipids, blood pressure, and long-term blood glucose (HbA1c), as well as longitudinal changes in BMI and WHR, differed according to menopausal status at baseline in 9,882 UK Biobank females (age range 40-70 years, n premenopausal = 3,529, n postmenopausal = 6,353). Furthermore, we examined whether these cardiometabolic factors were associated with WMH outcomes at the follow-up assessment, on average 8.78 years after baseline.

Results

Postmenopausal females showed higher levels of baseline blood lipids (HDL β  = 0.14, p < 0.001, LDL β  = 0.20, p < 0.001, triglycerides β  = 0.12, p < 0.001) and HbA1c (β  = 0.24, p < 0.001) compared to premenopausal women, beyond the effects of age. Over time, BMI increased more in the premenopausal compared to the postmenopausal group (β  = -0.08, p < 0.001), while WHR increased to a similar extent in both groups (β  = -0.03, p = 0.102). The change in WHR was however driven by increased waist circumference only in the premenopausal group. While the group level changes in BMI and WHR were in general small, these findings point to distinct anthropometric changes in pre- and postmenopausal females over time. Higher baseline measures of BMI, WHR, triglycerides, blood pressure, and HbA1c, as well as longitudinal increases in BMI and WHR, were associated with larger WMH volumes (β range = 0.03-0.13, p ≤ 0.002). HDL showed a significant inverse relationship with WMH volume (β  = -0.27, p < 0.001).

Discussion

Our findings emphasise the importance of monitoring cardiometabolic risk factors in females from midlife through the menopause transition and into the postmenopausal phase, to ensure improved cerebrovascular outcomes in later years.

Genetic perspectives on childhood monogenic diabetes: Diagnosis, management, and future directions

Monogenic diabetes is caused by one or even more genetic variations, which may be uncommon yet have a significant influence and cause diabetes at an early age. Monogenic diabetes affects 1 to 5% of children, and early detection and gene-tically focused treatment of neonatal diabetes and maturity-onset diabetes of the young can significantly improve long-term health and well-being. The etiology of monogenic diabetes in childhood is primarily attributed to genetic variations affecting the regulatory genes responsible for beta-cell activity. In rare instances, mutations leading to severe insulin resistance can also result in the development of diabetes. Individuals diagnosed with specific types of monogenic diabetes, which are commonly found, can transition from insulin therapy to sulfonylureas, provided they maintain consistent regulation of their blood glucose levels. Scientists have successfully devised materials and methodologies to distinguish individuals with type 1 or 2 diabetes from those more prone to monogenic diabetes. Genetic screening with appropriate findings and interpretations is essential to establish a prognosis and to guide the choice of therapies and management of these interrelated ailments. This review aims to design a comprehensive literature summarizing genetic insights into monogenetic diabetes in children and adolescents as well as summarizing their diagnosis and mana-gement.

Selenoprotein deficiency disorder predisposes to aortic aneurysm formation

Aortic aneurysms, which may dissect or rupture acutely and be lethal, can be a part of multisystem disorders that have a heritable basis. We report four patients with deficiency of selenocysteine-containing proteins due to selenocysteine Insertion Sequence Binding Protein 2 (SECISBP2) mutations who show early-onset, progressive, aneurysmal dilatation of the ascending aorta due to cystic medial necrosis. Zebrafish and male mice with global or vascular smooth muscle cell (VSMC)-targeted disruption of Secisbp2 respectively show similar aortopathy. Aortas from patients and animal models exhibit raised cellular reactive oxygen species, oxidative DNA damage and VSMC apoptosis. Antioxidant exposure or chelation of iron prevents oxidative damage in patient's cells and aortopathy in the zebrafish model. Our observations suggest a key role for oxidative stress and cell death, including via ferroptosis, in mediating aortic degeneration.

Lack of Evidence for the Role of the p.(Ser96Ala) Polymorphism in Histidine-Rich Calcium Binding Protein as a Secondary Hit in Cardiomyopathies

Inherited forms of arrhythmogenic and dilated cardiomyopathy (ACM and DCM) are characterized by variable disease expression and age-related penetrance. Calcium (Ca2+) is crucially important for proper cardiac function, and dysregulation of Ca2+ homeostasis seems to underly cardiomyopathy etiology. A polymorphism, c.286T>G p.(Ser96Ala), in the gene encoding the histidine-rich Ca2+ binding (HRC) protein, relevant for sarcoplasmic reticulum Ca2+ cycling, has previously been associated with a marked increased risk of life-threatening arrhythmias among idiopathic DCM patients. Following this finding, we investigated whether p.(Ser96Ala) affects major cardiac disease manifestations in carriers of the phospholamban (PLN) c.40_42delAGA; p.(Arg14del) pathogenic variant (cohort 1); patients diagnosed with, or predisposed to, ACM (cohort 2); and DCM patients (cohort 3). We found that the allele frequency of the p.(Ser96Ala) polymorphism was similar across the general European-American population (control cohort, 40.3-42.2%) and the different cardiomyopathy cohorts (cohorts 1-3, 40.9-43.9%). Furthermore, the p.(Ser96Ala) polymorphism was not associated with life-threatening arrhythmias or heart failure-related events across various patient cohorts. We therefore conclude that there is a lack of evidence supporting the important role of the HRC p.(Ser96Ala) polymorphism as a modifier in cardiomyopathy, refuting previous findings. Further research is required to identify bona fide genomic predictors for the stratification of cardiomyopathy patients and their risk for life-threatening outcomes.

Decoding Genetics, Ancestry, and Geospatial Context for Precision Health

Mass General Brigham, an integrated healthcare system based in the Greater Boston area of Massachusetts, annually serves 1.5 million patients. We established the Mass General Brigham Biobank (MGBB), encompassing 142,238 participants, to unravel the intricate relationships among genomic profiles, environmental context, and disease manifestations within clinical practice. In this study, we highlight the impact of ancestral diversity in the MGBB by employing population genetics, geospatial assessment, and association analyses of rare and common genetic variants. The population structures captured by the genetics mirror the sequential immigration to the Greater Boston area throughout American history, highlighting communities tied to shared genetic and environmental factors. Our investigation underscores the potency of unbiased, large-scale analyses in a healthcare-affiliated biobank, elucidating the dynamic interplay across genetics, immigration, structural geospatial factors, and health outcomes in one of the earliest American sites of European colonization.

Exploring the complex spectrum of dominance and recessiveness in genetic cardiomyopathies

Discrete categorization of Mendelian disease genes into dominant and recessive models often oversimplifies their underlying genetic architecture. Cardiomyopathies (CMs) are genetic diseases with complex etiologies for which an increasing number of recessive associations have recently been proposed. Here, we comprehensively analyze all published evidence pertaining to biallelic variation associated with CM phenotypes to identify high-confidence recessive genes and explore the spectrum of monoallelic and biallelic variant effects in established recessive and dominant disease genes. We classify 18 genes with robust recessive association with CMs, largely characterized by dilated phenotypes, early disease onset and severe outcomes. Several of these genes have monoallelic association with disease outcomes and cardiac traits in the UK Biobank, including LMOD2 and ALPK3 with dilated and hypertrophic CM, respectively. Our data provide insights into the complex spectrum of dominance and recessiveness in genetic heart disease and demonstrate how such approaches enable the discovery of unexplored genetic associations.

Novel MYBPC3 Mutations in Indian Population with Cardiomyopathies

Background

Mutations in Myosin Binding Protein C (MYBPC3) are one of the most frequent causes of cardiomyopathies in the world, but not much data are available in India.

Methods

We carried out targeted direct sequencing of MYBPC3 in 115 hypertrophic (HCM) and 127 dilated (DCM) cardiomyopathies against 197 ethnically matched healthy controls from India.

Results

We detected 34 single nucleotide variations in MYBPC3, of which 19 were novel. We found a splice site mutation [(IVS6+2T) T>G] and 16 missense mutations in Indian cardiomyopathies [5 in HCM; E258K, T262S, H287L, R408M, V483A: 4 in DCM; T146N, V321L, A392T, E393K and 7 in both HCM and DCM; L104M, V158M, S236G, R272C, T290A, G522E, A626V], but those were absent in 197 normal healthy controls. Interestingly, we found 7 out of 16 missense mutations (V158M, E258K, R272C, A392T, V483A, G522E, and A626V) in MYBPC3 were altering the evolutionarily conserved native amino acids, accounted for 8.7% and 6.3% in HCM and DCM, respectively. The bioinformatic tools predicted that those 7 missense mutations were pathogenic. Moreover, the co-segregation of those 7 mutations in families further confirmed their pathogenicity. Remarkably, we also identified compound mutations within the MYBPC3 gene of 6 cardiomyopathy patients (5%) with more severe disease phenotype; of which, 3 were HCM (2.6%) [(1. K244K + E258K + (IVS6+2T) T>G); (2. L104M + G522E + A626V); (3. P186P + G522E + A626V]; and 3 were DCM (2.4%) [(1. 5'UTR + A392T; 2. V158M+G522E; and 3.V158M + T262T + A626V].

Conclusion

The present comprehensive study on MYBPC3 has revealed both single and compound mutations in MYBPC3 and their association with disease in Indian Population with Cardiomyopathies. Our findings may perhaps help in initiating diagnostic strategies and eventually recognizing the targets for therapeutic interventions.

Genetic Burden of TNNI3K in Diagnostic Testing of Patients With Dilated Cardiomyopathy and Supraventricular Arrhythmias

BACKGROUND

Genetic variants in TNNI3K (troponin-I interacting kinase) have previously been associated with dilated cardiomyopathy (DCM), cardiac conduction disease, and supraventricular tachycardias. However, the link between TNNI3K variants and these cardiac phenotypes shows a lack of consensus concerning phenotype and protein function.

METHODS

We describe a systematic retrospective study of a cohort of patients undergoing genetic testing for cardiac arrhythmias and cardiomyopathy including TNNI3K. We further performed burden testing of TNNI3K in the UK Biobank. For 2 novel TNNI3K variants, we tested cosegregation. TNNI3K kinase function was estimated by TNNI3K autophosphorylation assays.

RESULTS

We demonstrate enrichment of rare coding TNNI3K variants in DCM patients in the Amsterdam cohort. In the UK Biobank, we observed an association between TNNI3K missense (but not loss-of-function) variants and DCM and atrial fibrillation. Furthermore, we demonstrate genetic segregation for 2 rare variants, TNNI3K-p.Ile512Thr and TNNI3K-p.His592Tyr, with phenotypes consisting of DCM, cardiac conduction disease, and supraventricular tachycardia, together with increased autophosphorylation. In contrast, TNNI3K-p.Arg556_Asn590del, a likely benign variant, demonstrated depleted autophosphorylation.

CONCLUSIONS

Our findings demonstrate an increased burden of rare coding TNNI3K variants in cardiac patients with DCM. Furthermore, we present 2 novel likely pathogenic TNNI3K variants with increased autophosphorylation, suggesting that enhanced autophosphorylation is likely to drive pathogenicity.

Common and rare variants associated with cardiometabolic traits across 98,622 whole-genome sequences in the All of Us research program

All of Us is a biorepository aiming to advance biomedical research by providing various types of data in diverse human populations. Here we present a demonstration project validating the program's genomic data in 98,622 participants. We sought to replicate known genetic associations for three diseases (atrial fibrillation [AF], coronary artery disease, type 2 diabetes [T2D]) and two quantitative traits (height and low-density lipoprotein [LDL]) by conducting common and rare variant analyses. We identified one known risk locus for AF, five loci for T2D, 143 loci for height, and nine loci for LDL. In gene-based burden tests for rare loss-of-function variants, we replicated associations between TTN and AF, GIGYF1 and T2D, ADAMTS17, ACAN, NPR2 and height, APOB, LDLR, PCSK9 and LDL. Our results are consistent with previous literature, indicating that the All of Us program is a reliable resource for advancing the understanding of complex diseases in diverse human populations.

Genetic and epigenetic background of diabetic kidney disease

Diabetic kidney disease (DKD) is a severe diabetic complication that affects up to half of the individuals with diabetes. Elevated blood glucose levels are a key underlying cause of DKD, but DKD is a complex multifactorial disease, which takes years to develop. Family studies have shown that inherited factors also contribute to the risk of the disease. During the last decade, genome-wide association studies (GWASs) have emerged as a powerful tool to identify genetic risk factors for DKD. In recent years, the GWASs have acquired larger number of participants, leading to increased statistical power to detect more genetic risk factors. In addition, whole-exome and whole-genome sequencing studies are emerging, aiming to identify rare genetic risk factors for DKD, as well as epigenome-wide association studies, investigating DNA methylation in relation to DKD. This article aims to review the identified genetic and epigenetic risk factors for DKD.

The impact of rare protein coding genetic variation on adult cognitive function

Compelling evidence suggests that human cognitive function is strongly influenced by genetics. Here, we conduct a large-scale exome study to examine whether rare protein-coding variants impact cognitive function in the adult population (n = 485,930). We identify eight genes (ADGRB2, KDM5B, GIGYF1, ANKRD12, SLC8A1, RC3H2, CACNA1A and BCAS3) that are associated with adult cognitive function through rare coding variants with large effects. Rare genetic architecture for cognitive function partially overlaps with that of neurodevelopmental disorders. In the case of KDM5B we show how the genetic dosage of one of these genes may determine the variability of cognitive, behavioral and molecular traits in mice and humans. We further provide evidence that rare and common variants overlap in association signals and contribute additively to cognitive function. Our study introduces the relevance of rare coding variants for cognitive function and unveils high-impact monogenic contributions to how cognitive function is distributed in the normal adult population.

Systematic elucidation of genetic mechanisms underlying cholesterol uptake

Genetic variation contributes greatly to LDL cholesterol (LDL-C) levels and coronary artery disease risk. By combining analysis of rare coding variants from the UK Biobank and genome-scale CRISPR-Cas9 knockout and activation screening, we substantially improve the identification of genes whose disruption alters serum LDL-C levels. We identify 21 genes in which rare coding variants significantly alter LDL-C levels at least partially through altered LDL-C uptake. We use co-essentiality-based gene module analysis to show that dysfunction of the RAB10 vesicle transport pathway leads to hypercholesterolemia in humans and mice by impairing surface LDL receptor levels. Further, we demonstrate that loss of function of OTX2 leads to robust reduction in serum LDL-C levels in mice and humans by increasing cellular LDL-C uptake. Altogether, we present an integrated approach that improves our understanding of the genetic regulators of LDL-C levels and provides a roadmap for further efforts to dissect complex human disease genetics.

Investigating Gene-Diet Interactions Impacting the Association Between Macronutrient Intake and Glycemic Traits

Few studies have demonstrated reproducible gene-diet interactions (GDIs) impacting metabolic disease risk factors, likely due in part to measurement error in dietary intake estimation and insufficient capture of rare genetic variation. We aimed to identify GDIs across the genetic frequency spectrum impacting the macronutrient-glycemia relationship in genetically and culturally diverse cohorts. We analyzed 33,187 participants free of diabetes from 10 National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program cohorts with whole-genome sequencing, self-reported diet, and glycemic trait data. We fit cohort-specific, multivariable-adjusted linear mixed models for the effect of diet, modeled as an isocaloric substitution of carbohydrate for fat, and its interactions with common and rare variants genome-wide. In main effect meta-analyses, participants consuming more carbohydrate had modestly lower glycemic trait values (e.g., for glycated hemoglobin [HbA1c], -0.013% HbA1c/250 kcal substitution). In GDI meta-analyses, a common African ancestry-enriched variant (rs79762542) reached study-wide significance and replicated in the UK Biobank cohort, indicating a negative carbohydrate-HbA1c association among major allele homozygotes only. Simulations revealed that >150,000 samples may be necessary to identify similar macronutrient GDIs under realistic assumptions about effect size and measurement error. These results generate hypotheses for further exploration of modifiable metabolic disease risk in additional cohorts with African ancestry.

ARTICLE HIGHLIGHTS

We aimed to identify genetic modifiers of the dietary macronutrient-glycemia relationship using whole-genome sequence data from 10 Trans-Omics for Precision Medicine program cohorts. Substitution models indicated a modest reduction in glycemia associated with an increase in dietary carbohydrate at the expense of fat. Genome-wide interaction analysis identified one African ancestry-enriched variant near the FRAS1 gene that may interact with macronutrient intake to influence hemoglobin A1c. Simulation-based power calculations accounting for measurement error suggested that substantially larger sample sizes may be necessary to discover further gene-macronutrient interactions.

Clinical case study meets population cohort: identification of a BRCA1 pathogenic founder variant in Orcadians

We multiply ascertained the BRCA1 pathogenic missense variant c.5207T > C; p.Val1736Ala (V1736A) in clinical investigation of breast and ovarian cancer families from Orkney in the Northern Isles of Scotland, UK. We sought to investigate the frequency and clinical relevance of this variant in those of Orcadian ancestry as an exemplar of the value of population cohorts in clinical care, especially in isolated populations. Oral history and birth, marriage and death registrations indicated genealogical linkage of the clinical cases to ancestors from the Isle of Westray, Orkney. Further clinical cases were identified through targeted testing for V1736A in women of Orcadian ancestry attending National Health Service (NHS) genetic clinics for breast and ovarian cancer family risk assessments. The variant segregates with female breast and ovarian cancer in clinically ascertained cases. Separately, exome sequence data from 2088 volunteer participants with three or more Orcadian grandparents, in the ORCADES research cohort, was interrogated to estimate the population prevalence of V1736A in Orcadians. The effects of the variant were assessed using Electronic Health Record (EHR) linkage. Twenty out of 2088 ORCADES research volunteers (~1%) carry V1736A, with a common haplotype around the variant. This allele frequency is ~480-fold higher than in UK Biobank participants. Cost-effectiveness of population screening for BRCA1 founder pathogenic variants has been demonstrated at a carrier frequency below the ~1% observed here. Thus we suggest that Orcadian women should be offered testing for the BRCA1 V1736A founder pathogenic variant, starting with those with known Westray ancestry.

Genomic analysis of lean individuals with NAFLD identifies monogenic disorders in a prospective cohort study

Background & Aims

Lean patients with non-alcoholic fatty liver disease (NAFLD) represent 10-20% of the affected population and may have heterogeneous drivers of disease. We have recently proposed the evaluation of patients with lean NAFLD without visceral adiposity for rare monogenic drivers of disease. Here, we aimed to validate this framework in a well-characterised cohort of patients with biopsy-proven NAFLD by performing whole exome sequencing.

Methods

This prospective study included 124 patients with biopsy-proven NAFLD and paired liver biopsies who underwent standardised research visits including advanced magnetic resonance imaging (MRI) assessment of liver fat and stiffness.

Results

Six patients with lean NAFLD were identified and underwent whole exome sequencing. Two lean patients (33%) were identified to have monogenic disorders. The lean patients with monogenic disorders had similar age, and anthropometric and MRI characteristics to lean patients without a monogenic disorder. Patient 1 harbours a rare homozygous pathogenic mutation in ALDOB (aldolase B) and was diagnosed with hereditary fructose intolerance. Patient 2 harbours a rare heterozygous mutation in apolipoprotein B (APOB). The pathogenicity of this APOB variant (p.Val1856CysfsTer2) was further validated in the UK Biobank and associated with lower circulating APOB levels (beta = -0.51 g/L, 95% CI -0.65 to -0.36 g/L, p = 1.4 × 10-11) and higher liver fat on MRI (beta = +10.4%, 95% CI 4.3-16.5%, p = 8.8 × 10-4). Hence, patient 2 was diagnosed with heterozygous familial hypobetalipoproteinaemia.

Conclusions

In this cohort of well-characterised patients with lean NAFLD without visceral adiposity, 33% (2/6) had rare monogenic drivers of disease, highlighting the importance of genomic analysis in this NAFLD subtype.

Impact and Implications

Although most people with non-alcoholic fatty liver disease (NAFLD) are overweight or obese, a subset are lean and may have unique genetic mutations that cause their fatty liver disease. We show that 33% of study participants with NAFLD who were lean harboured unique mutations that cause their fatty liver, and that these mutations had effects beyond the liver. This study demonstrates the value of genetic assessment of NAFLD in lean individuals to identify distinct subtypes of disease.

Adjusting for common variant polygenic scores improves yield in rare variant association analyses

With the emergence of large-scale sequencing data, methods for improving power in rare variant association tests are needed. Here we show that adjusting for common variant polygenic scores improves yield in gene-based rare variant association tests across 65 quantitative traits in the UK Biobank (up to 20% increase at α = 2.6 × 10-6), without marked increases in false-positive rates or genomic inflation. Benefits were seen for various models, with the largest improvements seen for efficient sparse mixed-effects models. Our results illustrate how polygenic score adjustment can efficiently improve power in rare variant association discovery.

Clinical and genetic associations of deep learning-derived cardiac magnetic resonance-based left ventricular mass

Left ventricular mass is a risk marker for cardiovascular events, and may indicate an underlying cardiomyopathy. Cardiac magnetic resonance is the gold-standard for left ventricular mass estimation, but is challenging to obtain at scale. Here, we use deep learning to enable genome-wide association study of cardiac magnetic resonance-derived left ventricular mass indexed to body surface area within 43,230 UK Biobank participants. We identify 12 genome-wide associations (1 known at TTN and 11 novel for left ventricular mass), implicating genes previously associated with cardiac contractility and cardiomyopathy. Cardiac magnetic resonance-derived indexed left ventricular mass is associated with incident dilated and hypertrophic cardiomyopathies, and implantable cardioverter-defibrillator implant. An indexed left ventricular mass polygenic risk score ≥90th percentile is also associated with incident implantable cardioverter-defibrillator implant in separate UK Biobank (hazard ratio 1.22, 95% CI 1.05-1.44) and Mass General Brigham (hazard ratio 1.75, 95% CI 1.12-2.74) samples. Here, we perform a genome-wide association study of cardiac magnetic resonance-derived indexed left ventricular mass to identify 11 novel variants and demonstrate that cardiac magnetic resonance-derived and genetically predicted indexed left ventricular mass are associated with incident cardiomyopathy.

Effects of protein-coding variants on blood metabolite measurements and clinical biomarkers in the UK Biobank

Genome-wide association studies (GWASs) have established the contribution of common and low-frequency variants to metabolic blood measurements in the UK Biobank (UKB). To complement existing GWAS findings, we assessed the contribution of rare protein-coding variants in relation to 355 metabolic blood measurements-including 325 predominantly lipid-related nuclear magnetic resonance (NMR)-derived blood metabolite measurements (Nightingale Health Plc) and 30 clinical blood biomarkers-using 412,393 exome sequences from four genetically diverse ancestries in the UKB. Gene-level collapsing analyses were conducted to evaluate a diverse range of rare-variant architectures for the metabolic blood measurements. Altogether, we identified significant associations (p < 1 × 10^-8) for 205 distinct genes that involved 1,968 significant relationships for the Nightingale blood metabolite measurements and 331 for the clinical blood biomarkers. These include associations for rare non-synonymous variants in PLIN1 and CREB3L3 with lipid metabolite measurements and SYT7 with creatinine, among others, which may not only provide insights into novel biology but also deepen our understanding of established disease mechanisms. Of the study-wide significant clinical biomarker associations, 40% were not previously detected on analyzing coding variants in a GWAS in the same cohort, reinforcing the importance of studying rare variation to fully understand the genetic architecture of metabolic blood measurements.

Association of Rare Protein-Truncating DNA Variants in APOB or PCSK9 With Low-density Lipoprotein Cholesterol Level and Risk of Coronary Heart Disease

Importance

Protein-truncating variants (PTVs) in apolipoprotein B (APOB) and proprotein convertase subtilisin/kexin type 9 (PCSK9) are associated with significantly lower low-density lipoprotein (LDL) cholesterol concentrations. The association of these PTVs with coronary heart disease (CHD) warrants further characterization in large, multiracial prospective cohort studies.

Objective

To evaluate the association of PTVs in APOB and PCSK9 with LDL cholesterol concentrations and CHD risk.

Design, Setting, and Participants

This studied included participants from 5 National Heart, Lung, and Blood Institute (NHLBI) studies and the UK Biobank. NHLBI study participants aged 5 to 84 years were recruited between 1971 and 2002 across the US and underwent whole-genome sequencing. UK Biobank participants aged 40 to 69 years were recruited between 2006 and 2010 in the UK and underwent whole-exome sequencing. Data were analyzed from June 2021 to October 2022.

Exposures

PTVs in APOB and PCSK9.

Main Outcomes and Measures

Estimated untreated LDL cholesterol levels and CHD.

Results

Among 19 073 NHLBI participants (10 598 [55.6%] female; mean [SD] age, 52 [17] years), 139 (0.7%) carried an APOB or PCSK9 PTV, which was associated with 49 mg/dL (95% CI, 43-56) lower estimated untreated LDL cholesterol level. Over a median (IQR) follow-up of 21.5 (13.9-29.4) years, incident CHD was observed in 12 of 139 carriers (8.6%) vs 3029 of 18 934 noncarriers (16.0%), corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.28-0.89; P = .02). Among 190 464 UK Biobank participants (104 831 [55.0%] female; mean [SD] age, 57 [8] years), 662 (0.4%) carried a PTV, which was associated with 45 mg/dL (95% CI, 42-47) lower estimated untreated LDL cholesterol level. Estimated CHD risk by age 75 years was 3.7% (95% CI, 2.0-5.3) in carriers vs 7.0% (95% CI, 6.9-7.2) in noncarriers, corresponding to an adjusted hazard ratio of 0.51 (95% CI, 0.32-0.81; P = .004).

Conclusions and Relevance

Among 209 537 individuals in this study, 0.4% carried an APOB or PCSK9 PTV that was associated with less exposure to LDL cholesterol and a 49% lower risk of CHD.

BIGKnock: fine-mapping gene-based associations via knockoff analysis of biobank-scale data

We propose BIGKnock (BIobank-scale Gene-based association test via Knockoffs), a computationally efficient gene-based testing approach for biobank-scale data, that leverages long-range chromatin interaction data, and performs conditional genome-wide testing via knockoffs. BIGKnock can prioritize causal genes over proxy associations at a locus. We apply BIGKnock to the UK Biobank data with 405,296 participants for multiple binary and quantitative traits, and show that relative to conventional gene-based tests, BIGKnock produces smaller sets of significant genes that contain the causal gene(s) with high probability. We further illustrate its ability to pinpoint potential causal genes at [Formula: see text] of the associated loci.

Systematic elucidation of genetic mechanisms underlying cholesterol uptake

Genetic variation contributes greatly to LDL cholesterol (LDL-C) levels and coronary artery disease risk. By combining analysis of rare coding variants from the UK Biobank and genome-scale CRISPR-Cas9 knockout and activation screening, we have substantially improved the identification of genes whose disruption alters serum LDL-C levels. We identify 21 genes in which rare coding variants significantly alter LDL-C levels at least partially through altered LDL-C uptake. We use co-essentiality-based gene module analysis to show that dysfunction of the RAB10 vesicle transport pathway leads to hypercholesterolemia in humans and mice by impairing surface LDL receptor levels. Further, we demonstrate that loss of function of OTX2 leads to robust reduction in serum LDL-C levels in mice and humans by increasing cellular LDL-C uptake. Altogether, we present an integrated approach that improves our understanding of genetic regulators of LDL-C levels and provides a roadmap for further efforts to dissect complex human disease genetics.

Cardiomyopathy prevalence exceeds 30% in individuals with TTN variants and early atrial fibrillation

PURPOSE

TTN truncating variants (TTNtvs) represent the largest known genetic cause of dilated cardiomyopathies (DCMs), however their penetrance for DCM in general populations is low. More broadly, patients with cardiomyopathies (CMs) often exhibit other cardiac conditions, such as atrial fibrillation (Afib), which has also been linked to TTNtvs. This retrospective analysis aims to characterize the relationship between different cardiac conditions in those with TTNtvs and identify individuals with the highest risk of DCM.

METHODS

In this work we leverage longitudinal electronic health record and exome sequencing data from approximately 450,000 individuals in 2 health systems to statistically confirm and pinpoint the genetic footprint of TTNtv-related diagnoses aside from CM, such as Afib, and determine whether vetting additional significantly associated phenotypes better stratifies CM risk across those with TTNtvs. We focused on TTNtvs in exons with a percentage spliced in >90% (hiPSI TTNtvs), a representation of constitutive cardiac expression.

RESULTS

When controlling for CM and Afib, other cardiac conditions retained only nominal association with TTNtvs. A sliding window analysis of TTNtvs across the locus confirms that the association is specific to hiPSI exons for both CM and Afib, with no meaningful associations in percent spliced in ≤90% exons (loPSI TTNtvs). The combination of hiPSI TTNtv status and early Afib diagnosis (before age 60) found a subset of TTNtv individuals at high risk for CM. The prevalence of CM in this subset was 33%, a rate that was 3.5 fold higher than that in individuals with hiPSI TTNtvs (9% prevalence), 5-fold higher than that in individuals without TTNtvs with early Afib (6% prevalence), and 80-fold higher than that in the general population.

CONCLUSION

Our retrospective analyses revealed that those with hiPSI TTNtvs and early Afib (∼1/2900) have a high prevalence of CM (33%), far exceeding that in other individuals with TTNtvs and in those without TTNtvs with an early Afib diagnosis. These results show that combining phenotypic information along with genomic population screening can identify patients at higher risk for progressing to symptomatic heart failure.

Cardiac Repolarization in Health and Disease

Abnormal cardiac repolarization is at the basis of life-threatening arrhythmias in various congenital and acquired cardiac diseases. Dysfunction of ion channels involved in repolarization at the cellular level are often the underlying cause of the repolarization abnormality. The expression pattern of the gene encoding the affected ion channel dictates its impact on the shape of the T-wave and duration of the QT interval, thereby setting the stage for both the occurrence of the trigger and the substrate for maintenance of the arrhythmia. Here we discuss how research into the genetic and electrophysiological basis of repolarization has provided us with insights into cardiac repolarization in health and disease and how this in turn may provide the basis for future improved patient-specific management.

Precision and genomic medicine for dilated and hypertrophic cardiomyopathy

Cardiomyopathy develops through an interaction of genetic and environmental factors. The clinical manifestations of both dilated cardiomyopathy and hypertrophic cardiomyopathy are diverse, but genetic testing defines the causative genes in about half of cases and can predict clinical prognosis. It has become clear that cardiomyopathy is caused not only by single rare variants but also by combinations of multiple common variants, and genome-wide genetic research is important for accurate disease risk assessment. Single-cell analysis research aimed at understanding the pathophysiology of cardiomyopathy is progressing rapidly, and it is expected that genomic analysis and single-cell molecular profiling will be combined to contribute to more detailed stratification of cardiomyopathy.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Damaging missense variants in IGF1R implicate a role for IGF-1 resistance in the etiology of type 2 diabetes

Type 2 diabetes (T2D) is a heritable metabolic disorder. While population studies have identified hundreds of common genetic variants associated with T2D, the role of rare (frequency < 0.1%) protein-coding variation is less clear. We performed exome sequence analysis in 418,436 (n = 32,374 T2D cases) individuals in the UK Biobank. We identified previously reported genes (GCK, GIGYF1, HNF1A) in addition to missense variants in ZEB2 (n = 31 carriers; odds ratio [OR] = 5.5 [95% confidence interval = 2.5-12.0]; p = 6.4 × 10^-7), MLXIPL (n = 245; OR = 2.3 [1.6-3.2]; p = 3.2 × 10^-7), and IGF1R (n = 394; OR = 2.4 [1.8-3.2]; p = 1.3 × 10^-10). Carriers of damaging missense variants within IGF1R were also shorter (-2.2 cm [-1.8 to -2.7]; p = 1.2 × 10^-19) and had higher circulating insulin-like growth factor-1 (IGF-1) protein levels (2.3 nmol/L [1.7-2.9]; p = 2.8 × 10^-14), indicating relative IGF-1 resistance. A likely causal role of IGF-1 resistance was supported by Mendelian randomization analyses using common variants. These results increase understanding of the genetic architecture of T2D and highlight the growth hormone/IGF-1 axis as a potential therapeutic target.

The Regulation of Triacylglycerol Metabolism and Lipoprotein Lipase Activity

Triacylglycerol (TG) metabolism is tightly regulated to maintain a pool of TG within circulating lipoproteins that can be hydrolyzed in a tissue-specific manner by lipoprotein lipase (LPL) to enable the delivery of fatty acids to adipose or oxidative tissues as needed. Elevated serum TG concentrations, which result from a deficiency of LPL activity or, more commonly, an imbalance in the regulation of tissue-specific LPL activities, have been associated with an increased risk of atherosclerotic cardiovascular disease through multiple studies. Among the most critical LPL regulators are the angiopoietin-like (ANGPTL) proteins ANGPTL3, ANGPTL4, and ANGPTL8, and a number of different apolipoproteins including apolipoprotein A5 (ApoA5), apolipoprotein C2 (ApoC2), and apolipoprotein C3 (ApoC3). These ANGPTLs and apolipoproteins work together to orchestrate LPL activity and therefore play pivotal roles in TG partitioning, hydrolysis, and utilization. This review summarizes the mechanisms of action, epidemiological findings, and genetic data most relevant to these ANGPTLs and apolipoproteins. The interplay between these important regulators of TG metabolism in both fasted and fed states is highlighted with a holistic view toward understanding key concepts and interactions. Strategies for developing safe and effective therapeutics to reduce circulating TG by selectively targeting these ANGPTLs and apolipoproteins are also discussed.

Systematic single-variant and gene-based association testing of thousands of phenotypes in 394,841 UK Biobank exomes

Genome-wide association studies have successfully discovered thousands of common variants associated with human diseases and traits, but the landscape of rare variations in human disease has not been explored at scale. Exome-sequencing studies of population biobanks provide an opportunity to systematically evaluate the impact of rare coding variations across a wide range of phenotypes to discover genes and allelic series relevant to human health and disease. Here, we present results from systematic association analyses of 4,529 phenotypes using single-variant and gene tests of 394,841 individuals in the UK Biobank with exome-sequence data. We find that the discovery of genetic associations is tightly linked to frequency and is correlated with metrics of deleteriousness and natural selection. We highlight biological findings elucidated by these data and release the dataset as a public resource alongside the Genebass browser for rapidly exploring rare-variant association results.

Identifying interpretable gene-biomarker associations with functionally informed kernel-based tests in 190,000 exomes

Here we present an exome-wide rare genetic variant association study for 30 blood biomarkers in 191,971 individuals in the UK Biobank. We compare gene-based association tests for separate functional variant categories to increase interpretability and identify 193 significant gene-biomarker associations. Genes associated with biomarkers were ~ 4.5-fold enriched for conferring Mendelian disorders. In addition to performing weighted gene-based variant collapsing tests, we design and apply variant-category-specific kernel-based tests that integrate quantitative functional variant effect predictions for missense variants, splicing and the binding of RNA-binding proteins. For these tests, we present a computationally efficient combination of the likelihood-ratio and score tests that found 36% more associations than the score test alone while also controlling the type-1 error. Kernel-based tests identified 13% more associations than their gene-based collapsing counterparts and had advantages in the presence of gain of function missense variants. We introduce local collapsing by amino acid position for missense variants and use it to interpret associations and identify potential novel gain of function variants in PIEZO1. Our results show the benefits of investigating different functional mechanisms when performing rare-variant association tests, and demonstrate pervasive rare-variant contribution to biomarker variability.

Endophenotype effect sizes support variant pathogenicity in monogenic disease susceptibility genes

Accurate and efficient classification of variant pathogenicity is critical for research and clinical care. Using data from three large studies, we demonstrate that population-based associations between rare variants and quantitative endophenotypes for three monogenic diseases (low-density-lipoprotein cholesterol for familial hypercholesterolemia, electrocardiographic QTc interval for long QT syndrome, and glycosylated hemoglobin for maturity-onset diabetes of the young) provide evidence for variant pathogenicity. Effect sizes are associated with pathogenic ClinVar assertions (P < 0.001 for each trait) and discriminate pathogenic from non-pathogenic variants (area under the curve 0.82-0.84 across endophenotypes). An effect size threshold of ≥ 0.5 times the endophenotype standard deviation nominates up to 35% of rare variants of uncertain significance or not in ClinVar in disease susceptibility genes with pathogenic potential. We propose that variant associations with quantitative endophenotypes for monogenic diseases can provide evidence supporting pathogenicity.