Consensus guidance for monitoring individuals with islet autoantibody-positive pre-stage 3 type 1 diabetes.

Given the proven benefits of screening to reduce diabetic ketoacidosis (DKA) likelihood at the time of stage 3 type 1 diabetes diagnosis, and emerging availability of therapy to delay disease progression, type 1 diabetes screening programmes are being increasingly emphasised. Once broadly implemented, screening initiatives will identify significant numbers of islet autoantibody-positive (IAb+) children and adults who are at risk of (confirmed single IAb+) or living with (multiple IAb+) early-stage (stage 1 and stage 2) type 1 diabetes. These individuals will need monitoring for disease progression; much of this care will happen in non-specialised settings. To inform this monitoring, JDRF in conjunction with international experts and societies developed consensus guidance. Broad advice from this guidance includes the following: (1) partnerships should be fostered between endocrinologists and primary-care providers to care for people who are IAb+; (2) when people who are IAb+ are initially identified there is a need for confirmation using a second sample; (3) single IAb+ individuals are at lower risk of progression than multiple IAb+ individuals; (4) individuals with early-stage type 1 diabetes should have periodic medical monitoring, including regular assessments of glucose levels, regular education about symptoms of diabetes and DKA, and psychosocial support; (5) interested people with stage 2 type 1 diabetes should be offered trial participation or approved therapies; and (6) all health professionals involved in monitoring and care of individuals with type 1 diabetes have a responsibility to provide education. The guidance also emphasises significant unmet needs for further research on early-stage type 1 diabetes to increase the rigour of future recommendations and inform clinical care.

Consensus Guidance for Monitoring Individuals With Islet Autoantibody-Positive Pre-Stage 3 Type 1 Diabetes.

Given the proven benefits of screening to reduce diabetic ketoacidosis (DKA) likelihood at the time of stage 3 type 1 diabetes diagnosis, and emerging availability of therapy to delay disease progression, type 1 diabetes screening programs are being increasingly emphasized. Once broadly implemented, screening initiatives will identify significant numbers of islet autoantibody-positive (IAb+) children and adults who are at risk for (confirmed single IAb+) or living with (multiple IAb+) early-stage (stage 1 and stage 2) type 1 diabetes. These individuals will need monitoring for disease progression; much of this care will happen in nonspecialized settings. To inform this monitoring, JDRF, in conjunction with international experts and societies, developed consensus guidance. Broad advice from this guidance includes the following: 1) partnerships should be fostered between endocrinologists and primary care providers to care for people who are IAb+; 2) when people who are IAb+ are initially identified, there is a need for confirmation using a second sample; 3) single IAb+ individuals are at lower risk of progression than multiple IAb+ individuals; 4) individuals with early-stage type 1 diabetes should have periodic medical monitoring, including regular assessments of glucose levels, regular education about symptoms of diabetes and DKA, and psychosocial support; 5) interested people with stage 2 type 1 diabetes should be offered trial participation or approved therapies; and 6) all health professionals involved in monitoring and care of individuals with type 1 diabetes have a responsibility to provide education. The guidance also emphasizes significant unmet needs for further research on early-stage type 1 diabetes to increase the rigor of future recommendations and inform clinical care.

Dysanapsis Genetic Risk Predicts Lung Function Across the Lifespan.

Rationale Dysanapsis refers to a mismatch between airway tree caliber and lung size arising early in life. Dysanapsis assessed by computed tomography (CT) is evident by early adulthood and associated with chronic obstructive pulmonary disease (COPD) risk later in life. Objective By examining the genetic factors associated with CT-assessed dysanapsis, we aimed to elucidate its molecular underpinnings and physiological significance across the lifespan. Methods We performed a genome-wide association study (GWAS) of CT-assessed dysanapsis in 11,951 adults, including individuals from two population-based and two COPD-enriched studies. We applied colocalization analysis to integrate GWAS and gene expression data from whole blood and lung. Genetic variants associated with dysanapsis were combined into a genetic risk score that was applied to examine association with lung function in children from a population-based birth cohort (n=1,278) and adults from the UK Biobank (n=369,157). Measurements and Main Results CT-assessed dysanapsis was associated with genetic variants from 21 independent signals in 19 gene regions, implicating HHIP, DSP, and NPNT as potential molecular targets based on colocalization of their expression. Higher dysanapsis genetic risk score was associated with obstructive spirometry among 5 year old children and among adults in the 5th, 6th and 7th decades of life. Conclusions CT-assessed dysanapsis is associated with variation in genes previously implicated in lung development and dysanapsis genetic risk is associated with obstructive lung function from early life through older adulthood. Dysanapsis may represent an endo-phenotype link between the genetic variations associated with lung function and COPD.

Epigenetic and proteomic signatures associate with clonal hematopoiesis expansion rate.

Clonal hematopoiesis of indeterminate potential (CHIP), whereby somatic mutations in hematopoietic stem cells confer a selective advantage and drive clonal expansion, not only correlates with age but also confers increased risk of morbidity and mortality. Here, we leverage genetically predicted traits to identify factors that determine CHIP clonal expansion rate. We used the passenger-approximated clonal expansion rate method to quantify the clonal expansion rate for 4,370 individuals in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) cohort and calculated polygenic risk scores for DNA methylation aging, inflammation-related measures and circulating protein levels. Clonal expansion rate was significantly associated with both genetically predicted and measured epigenetic clocks. No associations were identified with inflammation-related lab values or diseases and CHIP expansion rate overall. A proteome-wide search identified predicted circulating levels of myeloid zinc finger 1 and anti-Müllerian hormone as associated with an increased CHIP clonal expansion rate and tissue inhibitor of metalloproteinase 1 and glycine N-methyltransferase as associated with decreased CHIP clonal expansion rate. Together, our findings identify epigenetic and proteomic patterns associated with the rate of hematopoietic clonal expansion.

Plasma Proteomics of Exercise Blood Pressure and Incident Hypertension.

Importance: Blood pressure response during acute exercise (exercise blood pressure [EBP]) is associated with the future risk of hypertension and cardiovascular disease (CVD). Biochemical characterization of EBP could inform disease biology and identify novel biomarkers of future hypertension. Objective: To identify protein markers associated with EBP and test their association with incident hypertension. Design, Setting, and Participants: This study assayed 4977 plasma proteins in 681 healthy participants (from 763 assessed) of the Health, Risk Factors, Exercise Training and Genetics (HERITAGE; data collection from January 1993 to December 1997 and plasma proteomics from January 2019 to January 2020) Family Study at rest who underwent 2 cardiopulmonary exercise tests. Individuals were free of CVD at the time of recruitment. Individuals with resting SBP ≥160 mm Hg or DBP ≥100 mm Hg or taking antihypertensive drug therapy were excluded from the study. The association between resting plasma protein levels to both resting BP and EBP was evaluated. Proteins associated with EBP were analyzed for their association with incident hypertension in the Framingham Heart Study (FHS; n = 1177) and validated in the Jackson Heart Study (JHS; n = 772) and Multi-Ethnic Study of Atherosclerosis (MESA; n = 1367). Proteins associated with incident hypertension were tested for putative causal links in approximately 700 000 individuals using cis-protein quantitative loci mendelian randomization (cis-MR). Data were analyzed from January 2023 to January 2024. Exposures: Plasma proteins. Main Outcomes and Measures: EBP was defined as the BP response during a fixed workload (50 W) on a cycle ergometer. Hypertension was defined as BP ≥140/90 mm Hg or taking antihypertensive medication. Results: Among the 681 participants in the HERITAGE Family Study, the mean (SD) age was 34 (13) years; 366 participants (54%) were female; 238 (35%) were self-reported Black and 443 (65%) were self-reported White. Proteomic profiling of EBP revealed 34 proteins that would not have otherwise been identified through profiling of resting BP alone. Transforming growth factor ß receptor 3 (TGFBR3) and prostaglandin D2 synthase (PTGDS) had the strongest association with exercise systolic BP (SBP) and diastolic BP (DBP), respectively (TGFBR3: exercise SBP, ß estimate, -3.39; 95% CI, -4.79 to -2.00; P = 2.33 × 10-6; PTGDS: exercise DBP ß estimate, -2.50; 95% CI, -3.29 to -1.70; P = 1.18 × 10-9). In fully adjusted models, TGFBR3 was inversely associated with incident hypertension in FHS, JHS, and MESA (hazard ratio [HR]: FHS, 0.86; 95% CI, 0.75-0.97; P = .01; JHS, 0.87; 95% CI, 0.77-0.97; P = .02; MESA, 0.84; 95% CI, 0.71-0.98; P = .03; pooled cohort, 0.86; 95% CI, 0.79-0.92; P = 6 × 10-5). Using cis-MR, genetically predicted levels of TGFBR3 were associated with SBP, hypertension, and CVD events (SBP: ß, -0.38; 95% CI, -0.64 to -0.11; P = .006; hypertension: odds ratio [OR], 0.99; 95% CI, 0.98-0.99; P < .001; heart failure with hypertension: OR, 0.86; 95% CI, 0.77-0.97; P = .01; CVD: OR, 0.84; 95% CI, 0.77-0.92; P = 8 × 10-5; cerebrovascular events: OR, 0.77; 95% CI, 0.70-0.85; P = 5 × 10-7). Conclusions and Relevance: Plasma proteomic profiling of EBP identified a novel protein, TGFBR3, which may protect against elevated BP and long-term CVD outcomes.

Genetic mapping across autoimmune diseases reveals shared associations and mechanisms.

Autoimmune and inflammatory diseases are polygenic disorders of the immune system. Many genomic loci harbor risk alleles for several diseases, but the limited resolution of genetic mapping prevents determining whether the same allele is responsible, indicating a shared underlying mechanism. Here, using a collection of 129,058 cases and controls across 6 diseases, we show that ~40% of overlapping associations are due to the same allele. We improve fine-mapping resolution for shared alleles twofold by combining cases and controls across diseases, allowing us to identify more expression quantitative trait loci driven by the shared alleles. The patterns indicate widespread sharing of pathogenic mechanisms but not a single global autoimmune mechanism. Our approach can be applied to any set of traits and is particularly valuable as sample collections become depleted.

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.

Genetic Associations with C-peptide Levels before Type 1 Diabetes Diagnosis in At-Risk Relatives.

OBJECTIVE: We sought to determine whether the type 1 diabetes genetic risk score-2 (T1D-GRS2) and single nucleotide polymorphisms (SNPs) are associated with C-peptide preservation before type 1 diabetes diagnosis. METHODS: We conducted a retrospective analysis of 713 autoantibody-positive participants who developed type 1 diabetes in the TrialNet Pathway to Prevention Study who had T1DExomeChip data. We evaluated the relationships of 16 known SNPs and T1D-GRS2 with area under the curve (AUC) C-peptide levels during oral glucose tolerance tests conducted in the 9 months before diagnosis. RESULTS: Higher T1D-GRS2 was associated with lower C-peptide AUC in the 9 months before diagnosis in univariate (ß=-0.06, P<0.0001) and multivariate (ß=-0.03, P=0.005) analyses. Participants with the JAZF1 rs864745 T allele had lower C-peptide AUC in both univariate (ß=-0.11, P=0.002) and multivariate (ß=-0.06, P=0.018) analyses. CONCLUSIONS: The type 2 diabetes-associated JAZF1 rs864745 T allele and higher T1D-GRS2 are associated with lower C-peptide AUC prior to diagnosis of type 1 diabetes, with implications for the design of prevention trials.

LXR signaling pathways link cholesterol metabolism with risk for prediabetes and diabetes.

BACKGROUNDPreclinical studies suggest that cholesterol accumulation leads to insulin resistance. We previously reported that alterations in a monocyte cholesterol metabolism transcriptional network (CMTN) - suggestive of cellular cholesterol accumulation - were cross-sectionally associated with obesity and type 2 diabetes (T2D). Here, we sought to determine whether the CMTN alterations independently predict incident prediabetes/T2D risk, and correlate with cellular cholesterol accumulation.METHODSMonocyte mRNA expression of 11 CMTN genes was quantified among 934 Multi-Ethnic Study of Atherosclerosis (MESA) participants free of prediabetes/T2D; cellular cholesterol was measured in a subset of 24 monocyte samples.RESULTSDuring a median 6-year follow-up, lower expression of 3 highly correlated LXR target genes - ABCG1 and ABCA1 (cholesterol efflux) and MYLIP (cholesterol uptake suppression) - and not other CMTN genes, was significantly associated with higher risk of incident prediabetes/T2D. Lower expression of the LXR target genes correlated with higher cellular cholesterol levels (e.g., 47% of variance in cellular total cholesterol explained by ABCG1 expression). Further, adding the LXR target genes to overweight/obesity and other known predictors significantly improved prediction of incident prediabetes/T2D.CONCLUSIONThese data suggest that the aberrant LXR/ABCG1-ABCA1-MYLIP pathway (LAAMP) is a major T2D risk factor and support a potential role for aberrant LAAMP and cellular cholesterol accumulation in diabetogenesis.FUNDINGThe MESA Epigenomics and Transcriptomics Studies were funded by NIH grants 1R01HL101250, 1RF1AG054474, R01HL126477, R01DK101921, and R01HL135009. This work was supported by funding from NIDDK R01DK103531 and NHLBI R01HL119962.

Machine learning models for predicting blood pressure phenotypes by combining multiple polygenic risk scores.

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations.

Integrated analysis of blood DNA methylation, genetic variants, circulating proteins, microRNAs, and kidney failure in type 1 diabetes.

Variation in DNA methylation (DNAmet) in white blood cells and other cells/tissues has been implicated in the etiology of progressive diabetic kidney disease (DKD). However, the specific mechanisms linking DNAmet variation in blood cells with risk of kidney failure (KF) and utility of measuring blood cell DNAmet in personalized medicine are not clear. We measured blood cell DNAmet in 277 individuals with type 1 diabetes and DKD using Illumina EPIC arrays; 51% of the cohort developed KF during 7 to 20 years of follow-up. Our epigenome-wide analysis identified DNAmet at 17 CpGs (5'-cytosine-phosphate-guanine-3' loci) associated with risk of KF independent of major clinical risk factors. DNAmet at these KF-associated CpGs remained stable over a median period of 4.7 years. Furthermore, DNAmet variations at seven KF-associated CpGs were strongly associated with multiple genetic variants at seven genomic regions, suggesting a strong genetic influence on DNAmet. The effects of DNAmet variations at the KF-associated CpGs on risk of KF were partially mediated by multiple KF-associated circulating proteins and KF-associated circulating miRNAs. A prediction model for risk of KF was developed by adding blood cell DNAmet at eight selected KF-associated CpGs to the clinical model. This updated model significantly improved prediction performance (c-statistic = 0.93) versus the clinical model (c-statistic = 0.85) at P = 6.62 × 10-14. In conclusion, our multiomics study provides insights into mechanisms through which variation of DNAmet may affect KF development and shows that blood cell DNAmet at certain CpGs can improve risk prediction for KF in T1D.

Determinants of mosaic chromosomal alteration fitness.

Clonal hematopoiesis (CH) is characterized by the acquisition of a somatic mutation in a hematopoietic stem cell that results in a clonal expansion. These driver mutations can be single nucleotide variants in cancer driver genes or larger structural rearrangements called mosaic chromosomal alterations (mCAs). The factors that influence the variations in mCA fitness and ultimately result in different clonal expansion rates are not well understood. We used the Passenger-Approximated Clonal Expansion Rate (PACER) method to estimate clonal expansion rate as PACER scores for 6,381 individuals in the NHLBI TOPMed cohort with gain, loss, and copy-neutral loss of heterozygosity mCAs. Our mCA fitness estimates, derived by aggregating per-individual PACER scores, were correlated (R2 = 0.49) with an alternative approach that estimated fitness of mCAs in the UK Biobank using population-level distributions of clonal fraction. Among individuals with JAK2 V617F clonal hematopoiesis of indeterminate potential or mCAs affecting the JAK2 gene on chromosome 9, PACER score was strongly correlated with erythrocyte count. In a cross-sectional analysis, genome-wide association study of estimates of mCA expansion rate identified a TCL1A locus variant associated with mCA clonal expansion rate, with suggestive variants in NRIP1 and TERT.

Cardiovascular Burden of the V142I Transthyretin Variant.

Importance: Individual cohort studies concur that the amyloidogenic V142I variant of the transthyretin (TTR) gene, present in 3% to 4% of US Black individuals, increases heart failure (HF) and mortality risk. Precisely defining carrier risk across relevant clinical outcomes and estimating population burden of disease are important given established and emerging targeted treatments. Objectives: To better define the natural history of disease in carriers across mid to late life, assess variant modifiers, and estimate cardiovascular burden to the US population. Design, Setting, and Participants: A total of 23 338 self-reported Black participants initially free from HF were included in 4 large observational studies across the US (mean [SD], 15.5 [8.2] years of follow-up). Data analysis was performed between May 2023 and February 2024. Exposure: V142I carrier status (n = 754, 3.2%). Main Outcomes and Measures: Hospitalizations for HF (including subtypes of reduced and preserved ejection fraction) and all-cause mortality. Outcomes were analyzed by generating 10-year hazard ratios for each age between 50 and 90 years. Using actuarial methods, mean survival by carrier status was estimated and applied to the 2022 US population using US Census data. Results: Among the 23 338 participants, the mean (SD) age at baseline was 62 (9) years and 76.7% were women. Ten-year carrier risk increased for HF hospitalization by age 63 years, predominantly driven by HF with reduced ejection fraction, and 10-year all-cause mortality risk increased by age 72 years. Only age (but not sex or other select variables) modified risk with the variant, with estimated reductions in longevity ranging from 1.9 years (95% CI, 0.6-3.1) at age 50 to 2.8 years (95% CI, 2.0-3.6) at age 81. Based on these data, 435 851 estimated US Black carriers between ages 50 and 95 years are projected to cumulatively lose 957 505 years of life (95% CI, 534 475-1 380 535) due to the variant. Conclusions and Relevance: Among self-reported Black individuals, male and female V142I carriers faced similar and substantial risk for HF hospitalization, predominantly with reduced ejection fraction, and death, with steep age-dependent penetrance. Delineating the individual contributions of, and complex interplay among, the V142I variant, ancestry, the social construct of race, and biological or social determinants of health to cardiovascular disease merits further investigation.

Islet autoantibodies as precision diagnostic tools to characterize heterogeneity in type 1 diabetes: a systematic review.

BACKGROUND: Islet autoantibodies form the foundation for type 1 diabetes (T1D) diagnosis and staging, but heterogeneity exists in T1D development and presentation. We hypothesized that autoantibodies can identify heterogeneity before, at, and after T1D diagnosis, and in response to disease-modifying therapies. METHODS: We systematically reviewed PubMed and EMBASE databases (6/14/2022) assessing 10 years of original research examining relationships between autoantibodies and heterogeneity before, at, after diagnosis, and in response to disease-modifying therapies in individuals at-risk or within 1 year of T1D diagnosis. A critical appraisal checklist tool for cohort studies was modified and used for risk of bias assessment. RESULTS: Here we show that 152 studies that met extraction criteria most commonly characterized heterogeneity before diagnosis (91/152). Autoantibody type/target was most frequently examined, followed by autoantibody number. Recurring themes included correlations of autoantibody number, type, and titers with progression, differing phenotypes based on order of autoantibody seroconversion, and interactions with age and genetics. Only 44% specifically described autoantibody assay standardization program participation. CONCLUSIONS: Current evidence most strongly supports the application of autoantibody features to more precisely define T1D before diagnosis. Our findings support continued use of pre-clinical staging paradigms based on autoantibody number and suggest that additional autoantibody features, particularly in relation to age and genetic risk, could offer more precise stratification. To improve reproducibility and applicability of autoantibody-based precision medicine in T1D, we propose a methods checklist for islet autoantibody-based manuscripts which includes use of precision medicine MeSH terms and participation in autoantibody standardization workshops.

Islet autoantibodies are markers found in the blood when insulin-producing cells in the pancreas become damaged and can be used to predict future development of type 1 diabetes. We evaluated published literature to determine whether characteristics of islet antibodies (type, levels, numbers) could improve prediction and help understand differences in how individuals with type 1 diabetes respond to treatments. We found existing evidence shows that islet autoantibody type and number are most useful to predict disease progression before diagnosis. In addition, the age when islet autoantibodies first appear strongly influences rate of progression. These findings provide important information for patients and care providers on how islet autoantibodies can be used to understand future type 1 diabetes development and to identify individuals who have the potential to benefit from intervention or prevention therapy.

Longitudinal change in mitochondrial heteroplasmy exhibits positive selection for deleterious variants.

A common feature of human aging is the acquisition of somatic mutations, and mitochondria are particularly prone to mutation due to their inefficient DNA repair and close proximity to reactive oxygen species, leading to a state of mitochondrial DNA heteroplasmy1,2. Cross-sectional studies have demonstrated that detection of heteroplasmy increases with participant age3, a phenomenon that has been attributed to genetic drift4-7. In this first large-scale longitudinal study, we measured heteroplasmy in two prospective cohorts (combined n=1405) at two timepoints (mean time between visits, 8.6 years), demonstrating that deleterious heteroplasmies were more likely to increase in variant allele fraction (VAF). We further demonstrated that increase in VAF was associated with increased risk of overall mortality. These results challenge the claim that somatic mtDNA mutations arise mainly due to genetic drift, instead demonstrating positive selection for predicted deleterious mutations at the cellular level, despite an negative impact on overall mortality.

Time-to-Event Genome-Wide Association Study for Incident Cardiovascular Disease in People With Type 2 Diabetes.

OBJECTIVE: To identify genetic risk factors for incident cardiovascular disease (CVD) among people with type 2 diabetes (T2D). RESEARCH DESIGN AND METHODS: We conducted a multiancestry time-to-event genome-wide association study for incident CVD among people with T2D. We also tested 204 known coronary artery disease (CAD) variants for association with incident CVD. RESULTS: Among 49,230 participants with T2D, 8,956 had incident CVD events (event rate 18.2%). We identified three novel genetic loci for incident CVD: rs147138607 (near CACNA1E/ZNF648, hazard ratio [HR] 1.23, P = 3.6 × 10-9), rs77142250 (near HS3ST1, HR 1.89, P = 9.9 × 10-9), and rs335407 (near TFB1M/NOX3, HR 1.25, P = 1.5 × 10-8). Among 204 known CAD loci, 5 were associated with incident CVD in T2D (multiple comparison-adjusted P < 0.00024, 0.05/204). A standardized polygenic score of these 204 variants was associated with incident CVD with HR 1.14 (P = 1.0 × 10-16). CONCLUSIONS: The data point to novel and known genomic regions associated with incident CVD among individuals with T2D.

MagicalRsq-X: A cross-cohort transferable genotype imputation quality metric.

Since genotype imputation was introduced, researchers have been relying on the estimated imputation quality from imputation software to perform post-imputation quality control (QC). However, this quality estimate (denoted as Rsq) performs less well for lower-frequency variants. We recently published MagicalRsq, a machine-learning-based imputation quality calibration, which leverages additional typed markers from the same cohort and outperforms Rsq as a QC metric. In this work, we extended the original MagicalRsq to allow cross-cohort model training and named the new model MagicalRsq-X. We removed the cohort-specific estimated minor allele frequency and included linkage disequilibrium scores and recombination rates as additional features. Leveraging whole-genome sequencing data from TOPMed, specifically participants in the BioMe, JHS, WHI, and MESA studies, we performed comprehensive cross-cohort evaluations for predominantly European and African ancestral individuals based on their inferred global ancestry with the 1000 Genomes and Human Genome Diversity Project data as reference. Our results suggest MagicalRsq-X outperforms Rsq in almost every setting, with 7.3%-14.4% improvement in squared Pearson correlation with true R2, corresponding to 85-218 K variant gains. We further developed a metric to quantify the genetic distances of a target cohort relative to a reference cohort and showed that such metric largely explained the performance of MagicalRsq-X models. Finally, we found MagicalRsq-X saved up to 53 known genome-wide significant variants in one of the largest blood cell trait GWASs that would be missed using the original Rsq for QC. In conclusion, MagicalRsq-X shows superiority for post-imputation QC and benefits genetic studies by distinguishing well and poorly imputed lower-frequency variants.

Genome-wide analysis in over 1 million individuals of European ancestry yields improved polygenic risk scores for blood pressure traits.

Hypertension affects more than one billion people worldwide. Here we identify 113 novel loci, reporting a total of 2,103 independent genetic signals (P < 5 × 10-8) from the largest single-stage blood pressure (BP) genome-wide association study to date (n = 1,028,980 European individuals). These associations explain more than 60% of single nucleotide polymorphism-based BP heritability. Comparing top versus bottom deciles of polygenic risk scores (PRSs) reveals clinically meaningful differences in BP (16.9 mmHg systolic BP, 95% CI, 15.5-18.2 mmHg, P = 2.22 × 10-126) and more than a sevenfold higher odds of hypertension risk (odds ratio, 7.33; 95% CI, 5.54-9.70; P = 4.13 × 10-44) in an independent dataset. Adding PRS into hypertension-prediction models increased the area under the receiver operating characteristic curve (AUROC) from 0.791 (95% CI, 0.781-0.801) to 0.826 (95% CI, 0.817-0.836, ∆AUROC, 0.035, P = 1.98 × 10-34). We compare the 2,103 loci results in non-European ancestries and show significant PRS associations in a large African-American sample. Secondary analyses implicate 500 genes previously unreported for BP. Our study highlights the role of increasingly large genomic studies for precision health research.

Methylation patterns associated with C-reactive protein in racially and ethnically diverse populations.

Systemic low-grade inflammation is a feature of chronic disease. C-reactive protein (CRP) is a common biomarker of inflammation and used as an indicator of disease risk; however, the role of inflammation in disease is not completely understood. Methylation is an epigenetic modification in the DNA which plays a pivotal role in gene expression. In this study we evaluated differential DNA methylation patterns associated with blood CRP level to elucidate biological pathways and genetic regulatory mechanisms to improve the understanding of chronic inflammation. The racially and ethnically diverse participants in this study were included as 50% White, 41% Black or African American, 7% Hispanic or Latino/a, and 2% Native Hawaiian, Asian American, American Indian, or Alaska Native (total n = 13,433) individuals. We replicated 113 CpG sites from 87 unique loci, of which five were novel (CADM3, NALCN, NLRC5, ZNF792, and cg03282312), across a discovery set of 1,150 CpG sites associated with CRP level (p < 1.2E-7). The downstream pathways affected by DNA methylation included the identification of IFI16 and IRF7 CpG-gene transcript pairs which contributed to the innate immune response gene enrichment pathway along with NLRC5, NOD2, and AIM2. Gene enrichment analysis also identified the nuclear factor-kappaB transcription pathway. Using two-sample Mendelian randomization (MR) we inferred methylation at three CpG sites as causal for CRP levels using both White and Black or African American MR instrument variables. Overall, we identified novel CpG sites and gene transcripts that could be valuable in understanding the specific cellular processes and pathogenic mechanisms involved in inflammation.

A multi-ancestry genome-wide association study in type 1 diabetes.

Type 1 diabetes (T1D) is an autoimmune disease caused by destruction of the pancreatic ß-cells. Genome-wide association (GWAS) and fine mapping studies have been conducted mainly in European ancestry (EUR) populations. We performed a multi-ancestry GWAS to identify SNPs and HLA alleles associated with T1D risk and age at onset. EUR families (N = 3223), and unrelated individuals of African (AFR, N = 891) and admixed (Hispanic/Latino) ancestry (AMR, N = 308) were genotyped using the Illumina HumanCoreExome BeadArray, with imputation to the TOPMed reference panel. The Multi-Ethnic HLA reference panel was utilized to impute HLA alleles and amino acid residues. Logistic mixed models (T1D risk) and frailty models (age at onset) were used for analysis. In GWAS meta-analysis, seven loci were associated with T1D risk at genome-wide significance: PTPN22, HLA-DQA1, IL2RA, RNLS, INS, IKZF4-RPS26-ERBB3, and SH2B3, with four associated with T1D age at onset (PTPN22, HLA-DQB1, INS, and ERBB3). AFR and AMR meta-analysis revealed NRP1 as associated with T1D risk and age at onset, although NRP1 variants were not associated in EUR ancestry. In contrast, the PTPN22 variant was significantly associated with risk only in EUR ancestry. HLA alleles and haplotypes most significantly associated with T1D risk in AFR and AMR ancestry differed from that seen in EUR ancestry; in addition, the HLA-DRB1*08:02-DQA1*04:01-DQB1*04:02 haplotype was 'protective' in AMR while HLA-DRB1*08:01-DQA1*04:01-DQB1*04:02 haplotype was 'risk' in EUR ancestry, differing only at HLA-DRB1*08. These results suggest that much larger sample sizes in non-EUR populations are required to capture novel loci associated with T1D risk.