Whole genome sequencing based analysis of inflammation biomarkers in the Trans-Omics for Precision Medicine (TOPMed) consortium.

Inflammation biomarkers can provide valuable insight into the role of inflammatory processes in many diseases and conditions. Sequencing based analyses of such biomarkers can also serve as an exemplar of the genetic architecture of quantitative traits. To evaluate the biological insight, which can be provided by a multi-ancestry, whole-genome based association study, we performed a comprehensive analysis of 21 inflammation biomarkers from up to 38 465 individuals with whole-genome sequencing from the Trans-Omics for Precision Medicine (TOPMed) program (with varying sample size by trait, where the minimum sample size was n = 737 for MMP-1). We identified 22 distinct single-variant associations across 6 traits-E-selectin, intercellular adhesion molecule 1, interleukin-6, lipoprotein-associated phospholipase A2 activity and mass, and P-selectin-that remained significant after conditioning on previously identified associations for these inflammatory biomarkers. We further expanded upon known biomarker associations by pairing the single-variant analysis with a rare variant set-based analysis that further identified 19 significant rare variant set-based associations with 5 traits. These signals were distinct from both significant single variant association signals within TOPMed and genetic signals observed in prior studies, demonstrating the complementary value of performing both single and rare variant analyses when analyzing quantitative traits. We also confirm several previously reported signals from semi-quantitative proteomics platforms. Many of these signals demonstrate the extensive allelic heterogeneity and ancestry-differentiated variant-trait associations common for inflammation biomarkers, a characteristic we hypothesize will be increasingly observed with well-powered, large-scale analyses of complex traits.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown1. The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer2-4 and coronary heart disease5-this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP)6. Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Whole-exome sequencing study identifies four novel gene loci associated with diabetic kidney disease.

Diabetic kidney disease (DKD) is recognized as an important public health challenge. However, its genomic mechanisms are poorly understood. To identify rare variants for DKD, we conducted a whole-exome sequencing (WES) study leveraging large cohorts well-phenotyped for chronic kidney disease and diabetes. Our two-stage WES study included 4372 European and African ancestry participants from the Chronic Renal Insufficiency Cohort and Atherosclerosis Risk in Communities studies (stage 1) and 11 487 multi-ancestry Trans-Omics for Precision Medicine participants (stage 2). Generalized linear mixed models, which accounted for genetic relatedness and adjusted for age, sex and ancestry, were used to test associations between single variants and DKD. Gene-based aggregate rare variant analyses were conducted using an optimized sequence kernel association test implemented within our mixed model framework. We identified four novel exome-wide significant DKD-related loci through initiating diabetes. In single-variant analyses, participants carrying a rare, in-frame insertion in the DIS3L2 gene (rs141560952) exhibited a 193-fold increased odds [95% confidence interval (CI): 33.6, 1105] of DKD compared with noncarriers (P = 3.59 × 10-9). Likewise, each copy of a low-frequency KRT6B splice-site variant (rs425827) conferred a 5.31-fold higher odds (95% CI: 3.06, 9.21) of DKD (P = 2.72 × 10-9). Aggregate gene-based analyses further identified ERAP2 (P = 4.03 × 10-8) and NPEPPS (P = 1.51 × 10-7), which are both expressed in the kidney and implicated in renin-angiotensin-aldosterone system modulated immune response. In the largest WES study of DKD, we identified novel rare variant loci attaining exome-wide significance. These findings provide new insights into the molecular mechanisms underlying DKD.

Social and scientific motivations to move beyond groups in allele frequencies: The TOPMed experience.

For the genomics community, allele frequencies within defined groups (or "strata") are useful across multiple research and clinical contexts. Benefits include allowing researchers to identify populations for replication or "look up" studies, enabling researchers to compare population-specific frequencies to validate findings, and facilitating assessment of variant pathogenicity in clinical contexts. However, there are potential concerns with stratified allele frequencies. These include potential re-identification (determining whether or not an individual participated in a given research study based on allele frequencies and individual-level genetic data), harm from associating stigmatizing variants with specific groups, potential reification of race as a biological rather than a socio-political category, and whether presenting stratified frequencies-and the downstream applications that this presentation enables-is consistent with participants' informed consents. The NHLBI Trans-Omics for Precision Medicine (TOPMed) program considered the scientific and social implications of different approaches for adding stratified frequencies to the TOPMed BRAVO (Browse All Variants Online) variant server. We recommend a novel approach of presenting ancestry-specific allele frequencies using a statistical method based upon local genetic ancestry inference. Notably, this approach does not require grouping individuals by either predominant global ancestry or race/ethnicity and, therefore, mitigates re-identification and other concerns as the mixture distribution of ancestral allele frequencies varies across the genome. Here we describe our considerations and approach, which can assist other genomics research programs facing similar issues of how to define and present stratified frequencies in publicly available variant databases.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Polygenic Risk Score Associates With Atherosclerotic Plaque Characteristics at Autopsy.

BACKGROUND: Polygenic risk scores (PRSs) for coronary artery disease (CAD) potentially improve cardiovascular risk prediction. However, their relationship with histopathologic features of CAD has never been examined systematically. METHODS: From 4327 subjects referred to CVPath by the State of Maryland Office Chief Medical Examiner for sudden death between 1994 and 2015, 2455 cases were randomly selected for genotyping. We generated PRS from 291 known CAD risk loci. Detailed histopathologic examination of the coronary arteries was performed in all subjects. The primary study outcome measurements were histopathologic plaque features determining severity of atherosclerosis, including %stenosis, calcification, thin-cap fibroatheromas, and thrombotic CAD. RESULTS: After exclusion of cases with insufficient DNA sample quality or with missing data, 954 cases (mean age, 48.8±14.7 years; 75.7% men) remained in the final study cohort. Subjects in the highest PRS quintile exhibited more severe atherosclerosis compared with subjects in the lowest quintile, with greater %stenosis (80.3%±27.0% versus 50.4%±38.7%; adjusted P<0.001) and a higher frequency of calcification (69.6% versus 35.8%; adjusted P=0.004) and thin-cap fibroatheroma (26.7% versus 9.5%; adjusted P=0.007). Even after adjustment for traditional CAD risk factors, subjects within the highest PRS quintile had higher odds of severe atherosclerosis (ie, ≥75% stenosis; adjusted odds ratio, 3.77 [95% CI, 2.10-6.78]; P<0.001) and plaque rupture (adjusted odds ratio, 4.05 [95% CI, 2.26-7.24]; P<0.001). Moreover, subjects within the highest quintile had higher odds of CAD-associated cause of death, especially among those aged ≤50 years (adjusted odds ratio, 4.08 [95% CI, 2.01-8.30]; P<0.001). No statistically significant associations were observed with plaque erosion after adjusting for covariates. CONCLUSIONS: This is the first autopsy study investigating associations between PRS and atherosclerosis severity at the histopathologic level in subjects with sudden death. Our pathological analysis suggests PRS correlates with plaque burden and features of advanced atherosclerosis and may be useful as a method for CAD risk stratification, especially in younger subjects.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

An in-depth association analysis of genetic variants within nicotine-related loci: Meeting in middle of GWAS and genetic fine-mapping.

In the last two decades of Genome-wide association studies (GWAS), nicotine-dependence-related genetic loci (e.g., nicotinic acetylcholine receptor - nAChR subunit genes) are among the most replicable genetic findings. Although GWAS results have reported tens of thousands of SNPs within these loci, further analysis (e.g., fine-mapping) is required to identify the causal variants. However, it is computationally challenging for existing fine-mapping methods to reliably identify causal variants from thousands of candidate SNPs based on the posterior inclusion probability. To address this challenge, we propose a new method to select SNPs by jointly modeling the SNP-wise inference results and the underlying structured network patterns of the linkage disequilibrium (LD) matrix. We use adaptive dense subgraph extraction method to recognize the latent network patterns of the LD matrix and then apply group LASSO to select causal variant candidates. We applied this new method to the UK biobank data to identify the causal variant candidates for nicotine addiction. Eighty-one nicotine addiction-related SNPs (i.e.,-log(p) > 50) of nAChR were selected, which are highly correlated (average r2>0.8) although they are physically distant (e.g., >200 kilobase away) and from various genes. These findings revealed that distant SNPs from different genes can show higher LD r2 than their neighboring SNPs, and jointly contribute to a complex trait like nicotine addiction.

Missense variants in SORT1 are associated with LDL-C in an Amish population.

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

Monogenic and Polygenic Contributions to QTc Prolongation in the Population.

BACKGROUND: Rare sequence variation in genes underlying cardiac repolarization and common polygenic variation influence QT interval duration. However, current clinical genetic testing of individuals with unexplained QT prolongation is restricted to examination of monogenic rare variants. The recent emergence of large-scale biorepositories with sequence data enables examination of the joint contribution of rare and common variations to the QT interval in the population. METHODS: We performed a genome-wide association study of the QTc in 84 630 UK Biobank participants and created a polygenic risk score (PRS). Among 26 976 participants with whole-genome sequencing and ECG data in the TOPMed (Trans-Omics for Precision Medicine) program, we identified 160 carriers of putative pathogenic rare variants in 10 genes known to be associated with the QT interval. We examined QTc associations with the PRS and with rare variants in TOPMed. RESULTS: Fifty-four independent loci were identified by genome-wide association study in the UK Biobank. Twenty-one loci were novel, of which 12 were replicated in TOPMed. The PRS composed of 1 110 494 common variants was significantly associated with the QTc in TOPMed (ΔQTc/decile of PRS=1.4 ms [95% CI, 1.3 to 1.5]; P=1.1×10-196). Carriers of putative pathogenic rare variants had longer QTc than noncarriers (ΔQTc=10.9 ms [95% CI, 7.4 to 14.4]). Of individuals with QTc>480 ms, 23.7% carried either a monogenic rare variant or had a PRS in the top decile (3.4% monogenic, 21% top decile of PRS). CONCLUSIONS: QTc duration in the population is influenced by both rare variants in genes underlying cardiac repolarization and polygenic risk, with a sizeable contribution from polygenic risk. Comprehensive assessment of the genetic determinants of QTc prolongation includes incorporation of both polygenic and monogenic risk.

Associations of genome-wide and regional autozygosity with 96 complex traits in old order Amish.

BACKGROUND: Autozygosity, the proportion of the genome that is homozygous by descent, has been associated with variation in multiple health-related traits impacting evolutionary fitness. Autozygosity (FROH) is typically measured from runs of homozygosity (ROHs) that arise when identical-by-descent (IBD) haplotypes are inherited from each parent. Population isolates with a small set of common founders have elevated autozygosity relative to outbred populations. METHODS: In this study, we examined whether degree of autozygosity was associated with variation in 96 cardiometabolic traits among 7221 Old Order Amish individuals residing in Lancaster County, PA. We estimated the average length of an ROH segment to be 6350 KB, with each individual having on average 17.2 segments 1.5 KB or larger. Measurements of genome-wide and regional FROH were used as the primary predictors of trait variation in association analysis. RESULTS: In genome-wide FROH analysis, we did not identify any associations that withstood Bonferroni-correction (p = 0.0005). However, on regional FROH analysis, we identified associations exceeding genome-wide thresholds for two traits: serum bilirubin levels, which were significantly associated with a region on chromosome 2 localized to a region surrounding UGT1A10 (p = 1 × 10- 43), and HbA1c levels, which were significantly associated with a region on chromosome 8 localized near CHRNB3 (p = 8 × 10- 10). CONCLUSIONS: These analyses highlight the potential value of autozygosity mapping in founder populations.

Genome-wide association study of circulating interleukin 6 levels identifies novel loci.

Interleukin 6 (IL-6) is a multifunctional cytokine with both pro- and anti-inflammatory properties with a heritability estimate of up to 61%. The circulating levels of IL-6 in blood have been associated with an increased risk of complex disease pathogenesis. We conducted a two-staged, discovery and replication meta genome-wide association study (GWAS) of circulating serum IL-6 levels comprising up to 67 428 (ndiscovery = 52 654 and nreplication = 14 774) individuals of European ancestry. The inverse variance fixed effects based discovery meta-analysis, followed by replication led to the identification of two independent loci, IL1F10/IL1RN rs6734238 on chromosome (Chr) 2q14, (Pcombined = 1.8 × 10-11), HLA-DRB1/DRB5 rs660895 on Chr6p21 (Pcombined = 1.5 × 10-10) in the combined meta-analyses of all samples. We also replicated the IL6R rs4537545 locus on Chr1q21 (Pcombined = 1.2 × 10-122). Our study identifies novel loci for circulating IL-6 levels uncovering new immunological and inflammatory pathways that may influence IL-6 pathobiology.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts.

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

Deciphering the causal relationship between blood pressure and regional white matter integrity: A two-sample Mendelian randomization study.

Elevated arterial blood pressure (BP) is a common risk factor for cerebrovascular and cardiovascular diseases, but no causal relationship has been established between BP and cerebral white matter (WM) integrity. In this study, we performed a two-sample Mendelian randomization (MR) analysis with individual-level data by defining two nonoverlapping sets of European ancestry individuals (genetics-exposure set: N = 203,111; mean age = 56.71 years, genetics-outcome set: N = 16,156; mean age = 54.61 years) from UK Biobank to evaluate the causal effects of BP on regional WM integrity, measured by fractional anisotropy of diffusion tensor imaging. Two BP traits: systolic and diastolic blood pressure were used as exposures. Genetic variant was carefully selected as instrumental variable (IV) under the MR analysis assumptions. We existing large-scale genome-wide association study summary data for validation. The main method used was a generalized version of inverse-variance weight method while other MR methods were also applied for consistent findings. Two additional MR analyses were performed to exclude the possibility of reverse causality. We found significantly negative causal effects (FDR-adjusted p < .05; every 10 mmHg increase in BP leads to a decrease in FA value by .4% ~ 2%) of BP traits on a union set of 17 WM tracts, including brain regions related to cognitive function and memory. Our study extended the previous findings of association to causation for regional WM integrity, providing insights into the pathological processes of elevated BP that might chronically alter the brain microstructure in different regions.

Frontal white matter association with sleep quality and the role of stress.

An important measure of brain health is the integrity of white matter connectivity structures that link brain regions. Studies have found an association between poorer sleep quality and decreased white matter integrity. Stress is among the strongest predictors of sleep quality. This study aimed to evaluate the association between sleep quality and white matter and to test if the relationship persisted after accounting for stress. White matter microstructures were measured by diffusion tensor imaging in a population of Old Order Amish/Mennonite (N = 240). Sleep quality was determined by the Pittsburgh Sleep Quality Index. Current stress levels were measured by the perceived stress scale. Exposure to lifetime stress was measured by the lifetime stressor inventory. Microstructures of four white matter tracts: left and right anterior limbs of internal capsule, left anterior corona radiata, and genu of corpus callosum were significantly correlated with sleep quality (all p ≤ 0.001). The current stress level was a significant predictor of sleep quality (p ≤ 0.001) while lifetime stress was not. PSQI remained significantly associated with white matter integrity in these frontal tracts (all p < 0.01) after accounting for current stress and lifetime stress, while current and lifetime stress were not significant predictors of white matter in any of the four models. Sleep quality did not have any substantial mediation role between stress and white matter integrity. Sleep quality was significantly associated with several frontal white matter tracts that connect brain structures important for sleep regulation regardless of current or past stress levels.

Pharmacogenetics of sodium-glucose co-transporter-2 inhibitors: Validation of a sex-agnostic pharmacodynamic biomarker.

AIM: To validate pharmacodynamic responses to sodium-glucose co-transporter-2 (SGLT2) inhibitors and test for association with genetic variants in SLC5A4, SLC5A9, and SLC2A9. METHODS: Canagliflozin (300 mg), a SGLT2 inhibitor, was administered to 30 healthy volunteers. Several endpoints were measured to assess clinically relevant responses, including drug-induced increases in urinary excretion of glucose, sodium and uric acid. RESULTS: This pilot study confirmed that canagliflozin (300 mg) triggered acute changes in mean levels of several biomarkers: fasting plasma glucose (-4.1 mg/dL; P = 6 × 10-5 ), serum creatinine (+0.05 mg/dL; P = 8 × 10-4 ) and serum uric acid (-0.90 mg/dL; P = 5 × 10-10 ). The effects of sex on glucosuria depended upon how data were normalized. Whereas males' responses were ~60% greater when data were normalized to body surface area, males and females exhibited similar responses when glucosuria was expressed as grams of urinary glucose per gram-creatinine. The magnitude of glucosuria was not significantly correlated with fasting plasma glucose, estimated glomerular filtration rate or age in those healthy individuals without diabetes with an estimated glomerular filtration rate of more than 60 mL/min/1.73m2 . CONCLUSIONS: Normalizing data relative to creatinine excretion will facilitate including data from males and females in a single analysis. Furthermore, because our ongoing pharmacogenomic study (NCT02891954) is conducted in healthy individuals, this will facilitate detection of genetic associations with limited confounding by other factors such as HbA1c and renal function.

Acute pharmacodynamic responses to exenatide: Drug-induced increases in insulin secretion and glucose effectiveness.

AIM: Glucagon-like peptide-1 receptor agonists provide multiple benefits to patients with type 2 diabetes, including improved glycaemic control, weight loss and decreased risk of major adverse cardiovascular events. Because drug responses vary among individuals, we initiated investigations to identify genetic variants associated with the magnitude of drug responses. METHODS: Exenatide (5 µg, subcutaneously) or saline (0.2 ml, subcutaneously) was administered to 62 healthy volunteers. Frequently sampled intravenous glucose tolerance tests were conducted to assess the impact of exenatide on insulin secretion and insulin action. This pilot study was a crossover design in which participants received exenatide and saline in random order. RESULTS: Exenatide increased first phase insulin secretion 1.9-fold (p = 1.9 × 10-9 ) and accelerated the rate of glucose disappearance 2.4-fold (p = 2 × 10-10 ). Minimal model analysis showed that exenatide increased glucose effectiveness (Sg ) by 32% (p = .0008) but did not significantly affect insulin sensitivity (Si ). The exenatide-induced increase in insulin secretion made the largest contribution to interindividual variation in exenatide-induced acceleration of glucose disappearance while interindividual variation in the drug effect on Sg contributed to a lesser extent (ß = 0.58 or 0.27, respectively). CONCLUSIONS: This pilot study provides validation for the value of a frequently sampled intravenous glucose tolerance test (including minimal model analysis) to provide primary data for our ongoing pharmacogenomic study of pharmacodynamic effects of semaglutide (NCT05071898). Three endpoints provide quantitative assessments of the effects of glucagon-like peptide-1 receptor agonists on glucose metabolism: first phase insulin secretion, glucose disappearance rates and glucose effectiveness.

De novo mutations across 1,465 diverse genomes reveal mutational insights and reductions in the Amish founder population.

De novo mutations (DNMs), or mutations that appear in an individual despite not being seen in their parents, are an important source of genetic variation whose impact is relevant to studies of human evolution, genetics, and disease. Utilizing high-coverage whole-genome sequencing data as part of the Trans-Omics for Precision Medicine (TOPMed) Program, we called 93,325 single-nucleotide DNMs across 1,465 trios from an array of diverse human populations, and used them to directly estimate and analyze DNM counts, rates, and spectra. We find a significant positive correlation between local recombination rate and local DNM rate, and that DNM rate explains a substantial portion (8.98 to 34.92%, depending on the model) of the genome-wide variation in population-level genetic variation from 41K unrelated TOPMed samples. Genome-wide heterozygosity does correlate with DNM rate, but only explains <1% of variation. While we are underpowered to see small differences, we do not find significant differences in DNM rate between individuals of European, African, and Latino ancestry, nor across ancestrally distinct segments within admixed individuals. However, we did find significantly fewer DNMs in Amish individuals, even when compared with other Europeans, and even after accounting for parental age and sequencing center. Specifically, we found significant reductions in the number of C→A and T→C mutations in the Amish, which seem to underpin their overall reduction in DNMs. Finally, we calculated near-zero estimates of narrow sense heritability (h2), which suggest that variation in DNM rate is significantly shaped by nonadditive genetic effects and the environment.