Integration of rare expression outlier-associated variants improves polygenic risk prediction.

Polygenic risk scores (PRSs) quantify the contribution of multiple genetic loci to an individual's likelihood of a complex trait or disease. However, existing PRSs estimate this likelihood with common genetic variants, excluding the impact of rare variants. Here, we report on a method to identify rare variants associated with outlier gene expression and integrate their impact into PRS predictions for body mass index (BMI), obesity, and bariatric surgery. Between the top and bottom 10%, we observed a 20.8% increase in risk for obesity (p = 3 × 10-14), 62.3% increase in risk for severe obesity (p = 1 × 10-6), and median 5.29 years earlier onset for bariatric surgery (p = 0.008), as a function of expression outlier-associated rare variant burden when controlling for common variant PRS. We show that these predictions were more significant than integrating the effects of rare protein-truncating variants (PTVs), observing a mean 19% increase in phenotypic variance explained with expression outlier-associated rare variants when compared with PTVs (p = 2 × 10-15). We replicated these findings by using data from the Million Veteran Program and demonstrated that PRSs across multiple traits and diseases can benefit from the inclusion of expression outlier-associated rare variants identified through population-scale transcriptome sequencing.

A genome-wide association study in a large community-based cohort identifies multiple loci associated with susceptibility to bacterial and viral infections.

There is limited data on host-specific genetic determinants of susceptibility to bacterial and viral infections. Genome-wide association studies using large population cohorts can be a first step towards identifying patients prone to infectious diseases and targets for new therapies. Genetic variants associated with clinically relevant entities of bacterial and viral infections (e.g., abdominal infections, respiratory infections, and sepsis) in 337,484 participants of the UK Biobank cohort were explored by genome-wide association analyses. Cases (n = 81,179) were identified based on ICD-10 diagnosis codes of hospital inpatient and death registries. Functional annotation was performed using gene expression (eQTL) data. Fifty-seven unique genome-wide significant loci were found, many of which are novel in the context of infectious diseases. Some of the detected genetic variants were previously reported associated with infectious, inflammatory, autoimmune, and malignant diseases or key components of the immune system (e.g., white blood cells, cytokines). Fine mapping of the HLA region revealed significant associations with HLA-DQA1, HLA-DRB1, and HLA-DRB4 locus alleles. PPP1R14A showed strong colocalization with abdominal infections and gene expression in sigmoid and transverse colon, suggesting causality. Shared significant loci across infections and non-infectious phenotypes in the UK Biobank cohort were found, suggesting associations for example between SNPs identified for abdominal infections and CRP, rheumatoid arthritis, and diabetes mellitus. We report multiple loci associated with bacterial and viral infections. A better understanding of the genetic determinants of bacterial and viral infections can be useful to identify patients at risk and in the development of new drugs.

Population-scale tissue transcriptomics maps long non-coding RNAs to complex disease.

Long non-coding RNA (lncRNA) genes have well-established and important impacts on molecular and cellular functions. However, among the thousands of lncRNA genes, it is still a major challenge to identify the subset with disease or trait relevance. To systematically characterize these lncRNA genes, we used Genotype Tissue Expression (GTEx) project v8 genetic and multi-tissue transcriptomic data to profile the expression, genetic regulation, cellular contexts, and trait associations of 14,100 lncRNA genes across 49 tissues for 101 distinct complex genetic traits. Using these approaches, we identified 1,432 lncRNA gene-trait associations, 800 of which were not explained by stronger effects of neighboring protein-coding genes. This included associations between lncRNA quantitative trait loci and inflammatory bowel disease, type 1 and type 2 diabetes, and coronary artery disease, as well as rare variant associations to body mass index.

Impact of admixture and ancestry on eQTL analysis and GWAS colocalization in GTEx.

BACKGROUND: Population structure among study subjects may confound genetic association studies, and lack of proper correction can lead to spurious findings. The Genotype-Tissue Expression (GTEx) project largely contains individuals of European ancestry, but the v8 release also includes up to 15% of individuals of non-European ancestry. Assessing ancestry-based adjustments in GTEx improves portability of this research across populations and further characterizes the impact of population structure on GWAS colocalization. RESULTS: Here, we identify a subset of 117 individuals in GTEx (v8) with a high degree of population admixture and estimate genome-wide local ancestry. We perform genome-wide cis-eQTL mapping using admixed samples in seven tissues, adjusted by either global or local ancestry. Consistent with previous work, we observe improved power with local ancestry adjustment. At loci where the two adjustments produce different lead variants, we observe 31 loci (0.02%) where a significant colocalization is called only with one eQTL ancestry adjustment method. Notably, both adjustments produce similar numbers of significant colocalizations within each of two different colocalization methods, COLOC and FINEMAP. Finally, we identify a small subset of eQTL-associated variants highly correlated with local ancestry, providing a resource to enhance functional follow-up. CONCLUSIONS: We provide a local ancestry map for admixed individuals in the GTEx v8 release and describe the impact of ancestry and admixture on gene expression, eQTLs, and GWAS colocalization. While the majority of the results are concordant between local and global ancestry-based adjustments, we identify distinct advantages and disadvantages to each approach.

CRISPR-Cas9-mediated knockout of SPRY2 in human hepatocytes leads to increased glucose uptake and lipid droplet accumulation.

BACKGROUND: The prevalence of obesity and its comorbidities, including type 2 diabetes mellitus (T2DM), is dramatically increasing throughout the world; however, the underlying aetiology is incompletely understood. Genome-wide association studies (GWAS) have identified hundreds of genec susceptibility loci for obesity and T2DM, although the causal genes and mechanisms are largely unknown. SPRY2 is a candidate gene identified in GWAS of body fat percentage and T2DM, and has recently been linked to insulin production in pancreatic ß-cells. In the present study, we aimed to further understand SPRY2 via functional characterisation in HepG2 cells, an in vitro model of human hepatocytes widely used to investigate T2DM and insulin resistance. METHODS: CRISPR-Cas9 genome editing was used to target SPRY2 in HepG2 cells, and the functional consequences of SPRY2 knockout (KO) and overexpression subsequently assessed using glucose uptake and lipid droplet assays, measurement of protein kinase phosphorylation and RNA sequencing. RESULTS: The major functional consequence of SPRY2 KO was a significant increase in glucose uptake, along with elevated lipid droplet accumulation. These changes were attenuated, but not reversed, in cells overexpressing SPRY2. Phosphorylation of protein kinases across key signalling pathways (including Akt and mitogen activated protein kinases) was not altered after SPRY2 KO. Transcriptome profiling in SPRY2 KO and mock (control) cells revealed a number of differentially expressed genes related to cholesterol biosynthesis, cell cycle regulation and cellular signalling pathways. Phospholipase A2 group IIA (PLA2G2A) mRNA level was subsequently validated as significantly upregulated following SPRY2 KO, highlighting this as a potential mediator downstream of SPRY2. CONCLUSION: These findings suggest a role for SPRY2 in glucose and lipid metabolism in hepatocytes and contribute to clarifying the function of this gene in the context of metabolic diseases.

Polygenic risk scores in coronary artery disease.

PURPOSE OF REVIEW: Large genome-wide association studies (GWAS) have identified variants accounting for a substantial portion of the heritable risk for coronary artery disease (CAD). These studies have catalyzed drug discovery and generated the possibility of improved risk prediction and stratification. Here, we review the current state-of-the art in polygenic risk scores (PRSs) and look to the future, as these scores move towards clinical application. RECENT FINDINGS: Over the last decade, multilocus PRSs for CAD have expanded to include millions of variants and demonstrated strong association with CAD outcomes, even when adjusted for traditional risk factors. Recently, PRSs have shown better prediction of CAD outcomes than any single traditional risk factor alone. Advances in statistical methods used to generate PRSs have improved their predictive ability and transferability between populations with varied ancestries. Initial clinical studies have also demonstrated the potential of genetic information to impact shared decision-making between patients and providers, leading to improved outcomes. SUMMARY: PRSs can improve risk stratification for CAD especially in white/European populations and have the potential to alter routine clinical care. However, unlocking this potential will require additional research in PRSs in nonwhite populations and substantial investment in clinical implementation studies.

Large-Scale Phenome-Wide Association Study of PCSK9 Variants Demonstrates Protection Against Ischemic Stroke.

BACKGROUND: PCSK9 inhibition is a potent new therapy for hypercholesterolemia and cardiovascular disease. Although short-term clinical trial results have not demonstrated major adverse effects, long-term data will not be available for some time. Genetic studies in large biobanks offer a unique opportunity to predict drug effects and provide context for the evaluation of future clinical trial outcomes. METHODS: We tested the association of the PCSK9 missense variant rs11591147 with predefined phenotypes and phenome-wide, in 337 536 individuals of British ancestry in the UK Biobank, with independent discovery and replication. Using a Bayesian statistical method, we leveraged phenotype correlations to evaluate the phenome-wide impact of PCSK9 inhibition with higher power at a finer resolution. RESULTS: The T allele of rs11591147 showed a protective effect on hyperlipidemia (odds ratio, 0.63±0.04; P=2.32×10-38), coronary heart disease (odds ratio, 0.73±0.09; P=1.05×10-6), and ischemic stroke (odds ratio, 0.61±0.18; P=2.40×10-3) and was associated with increased type 2 diabetes mellitus risk adjusted for lipid-lowering medication status (odds ratio, 1.24±0.10; P=1.98×10-7). We did not observe associations with cataracts, heart failure, atrial fibrillation, and cognitive dysfunction. Leveraging phenotype correlations, we observed evidence of a protective association with cerebral infarction and vascular occlusion. These results explore the effects of direct PCSK9 inhibition; off-target effects cannot be predicted using this approach. CONCLUSIONS: This result represents the first genetic evidence in a large cohort for the protective effect of PCSK9 inhibition on ischemic stroke and corroborates exploratory evidence from clinical trials. PCSK9 inhibition was not associated with variables other than those related to LDL (low-density lipoprotein) cholesterol, atherosclerosis, and type 2 diabetes mellitus, suggesting that other effects are either small or absent.

Low-cost blast wave generator for studies of hearing loss and brain injury: blast wave effects in closed spaces.

BACKGROUND: Military personnel and civilians living in areas of armed conflict have increased risk of exposure to blast overpressures that can cause significant hearing loss and/or brain injury. The equipment used to simulate comparable blast overpressures in animal models within laboratory settings is typically very large and prohibitively expensive. NEW METHOD: To overcome the fiscal and space limitations introduced by previously reported blast wave generators, we developed a compact, low-cost blast wave generator to investigate the effects of blast exposures on the auditory system and brain. RESULTS: The blast wave generator was constructed largely from off the shelf components, and reliably produced blasts with peak sound pressures of up to 198dB SPL (159.3kPa) that were qualitatively similar to those produced from muzzle blasts or explosions. Exposure of adult rats to 3 blasts of 188dB peak SPL (50.4kPa) resulted in significant loss of cochlear hair cells, reduced outer hair cell function and a decrease in neurogenesis in the hippocampus. COMPARISON TO EXISTING METHODS: Existing blast wave generators are typically large, expensive, and are not commercially available. The blast wave generator reported here provides a low-cost method of generating blast waves in a typical laboratory setting. CONCLUSIONS: This compact blast wave generator provides scientists with a low cost device for investigating the biological mechanisms involved in blast wave injury to the rodent cochlea and brain that may model many of the damaging effects sustained by military personnel and civilians exposed to intense blasts.

Presurgical planning using image-based in silico anatomical and functional characterization of Tetralogy of Fallot with associated anomalies.

OBJECTIVES: The ideal vascular reconstruction strategy for anomalies associated with Tetralogy of Fallot (ToF) is often driven by observations made at the operating table. A method to conduct accurate studies to assess the virtues of a certain surgical technique to guide surgical decisions is found wanting. We hypothesize that patient-specific computed tomography (CT)-based morphometry followed by in silico reconstruction of viable surgical options with haemodynamic function assessment using computational fluid dynamics (CFD) can guide surgical decisions and help forecast functional outcomes without invasive measurements. METHODS: A ToF patient associated with additional left pulmonary artery (LPA) stenosis and a patent ductus arteriosus (PDA) who underwent a successful correction using a single pericardial patch (SPP) was selected as a reference for morphological characterization after 3D anatomical reconstruction from CT images. A second patient with morphological similarities established after scaled, co-registration with the reference patient was selected for virtual correction using the same strategy (i.e. SPP repair). CFD was employed for functional analysis of pulmonary artery (PA) pressure gradients in the baseline preoperative and virtually corrected models, using patient-specific cardiac output and Qp/Qs information. RESULTS: SPP repair was modelled in silico following surgical steps of PDA ligation, creation of an incision along the LPA and main PA (MPA) and finally suturing a rectangular SPP, effectively reducing MPA to LPA angle. Analysis of SPP repair revealed significant reduction in right ventricular outflow tract-LPA pressure gradient with improved left-right PA flow distribution in both patients. CONCLUSIONS: In silico surgery followed by CFD evaluation has the potential in augmenting morphology-guided decisions on surgical strategy and holds promise in preoperatively determining optimal intervention strategy. This is a paradigm-shifting concept to evaluate patient-specific anomalies in a manner more objective than mere visual inspection of anatomical traits from radiology images. Present studies are focused on an analysis with a larger patient cohort to establish a library of ToF patients' successful surgical outcomes to inform morphology-based selection of surgical strategy.