Predicting the presence of coronary plaques featuring high-risk characteristics using polygenic risk scores and targeted proteomics in patients with suspected coronary artery disease.

BACKGROUND: The presence of coronary plaques with high-risk characteristics is strongly associated with adverse cardiac events beyond the identification of coronary stenosis. Testing by coronary computed tomography angiography (CCTA) enables the identification of high-risk plaques (HRP). Referral for CCTA is presently based on pre-test probability estimates including clinical risk factors (CRFs); however, proteomics and/or genetic information could potentially improve patient selection for CCTA and, hence, identification of HRP. We aimed to (1) identify proteomic and genetic features associated with HRP presence and (2) investigate the effect of combining CRFs, proteomics, and genetics to predict HRP presence. METHODS: Consecutive chest pain patients (n = 1462) undergoing CCTA to diagnose obstructive coronary artery disease (CAD) were included. Coronary plaques were assessed using a semi-automatic plaque analysis tool. Measurements of 368 circulating proteins were obtained with targeted Olink panels, and DNA genotyping was performed in all patients. Imputed genetic variants were used to compute a multi-trait multi-ancestry genome-wide polygenic score (GPSMult). HRP presence was defined as plaques with two or more high-risk characteristics (low attenuation, spotty calcification, positive remodeling, and napkin ring sign). Prediction of HRP presence was performed using the glmnet algorithm with repeated fivefold cross-validation, using CRFs, proteomics, and GPSMult as input features. RESULTS: HRPs were detected in 165 (11%) patients, and 15 input features were associated with HRP presence. Prediction of HRP presence based on CRFs yielded a mean area under the receiver operating curve (AUC) ± standard error of 73.2 ± 0.1, versus 69.0 ± 0.1 for proteomics and 60.1 ± 0.1 for GPSMult. Combining CRFs with GPSMult increased prediction accuracy (AUC 74.8 ± 0.1 (P = 0.004)), while the inclusion of proteomics provided no significant improvement to either the CRF (AUC 73.2 ± 0.1, P = 1.00) or the CRF + GPSMult (AUC 74.6 ± 0.1, P = 1.00) models, respectively. CONCLUSIONS: In patients with suspected CAD, incorporating genetic data with either clinical or proteomic data improves the prediction of high-risk plaque presence. TRIAL REGISTRATION: https://clinicaltrials.gov/ct2/show/NCT02264717 (September 2014).

Obesity Variants in the GIPR Gene Are not Associated With Risk of Fracture or Bone Mineral Density.

CONTEXT: It is not clear if antagonizing the GIP (glucose-dependent insulinotropic polypeptide) receptor (GIPR) for treatment of obesity is likely to increase the risk of fractures, or to lower bone mineral density (BMD) beyond what is expected with rapid weight loss. OBJECTIVE: The objective of this study was to investigate the risk of fracture and BMD of sequence variants in GIPR that reduce the activity of the GIP receptor and have been associated with reduced body mass index (BMI). METHODS: We analyzed the association of 3 missense variants in GIPR, a common variant, rs1800437 (p.Glu354Gln), and 2 rare variants, rs139215588 (p.Arg190Gln) and rs143430880 (p.Glu288Gly), as well as a burden of predicted loss-of-function (LoF) variants with risk of fracture and with BMD in a large meta-analysis of up to 1.2 million participants. We analyzed associations with fractures at different skeletal sites in the general population: any fractures, hip fractures, vertebral fractures and forearm fractures, and specifically nonvertebral and osteoporotic fractures in postmenopausal women. We also evaluated associations with BMD at the lumbar spine, femoral neck, and total body measured with dual-energy x-ray absorptiometry (DXA), and with BMD estimated from heel ultrasound (eBMD). RESULTS: None of the 3 missense variants in GIPR was significantly associated with increased risk of fractures or with lower BMD. Burden of LoF variants in GIPR was not associated with fractures or with BMD measured with clinically validated DXA, but was associated with eBMD. CONCLUSION: Missense variants in GIPR, or burden of LoF variants in the gene, are not associated with risk of fractures or with lower BMD.

Variant in the synaptonemal complex protein SYCE2 associates with pregnancy loss through effect on recombination.

Two-thirds of all human conceptions are lost, in most cases before clinical detection. The lack of detailed understanding of the causes of pregnancy losses constrains focused counseling for future pregnancies. We have previously shown that a missense variant in synaptonemal complex central element protein 2 (SYCE2), in a key residue for the assembly of the synaptonemal complex backbone, associates with recombination traits. Here we show that it also increases risk of pregnancy loss in a genome-wide association analysis on 114,761 women with reported pregnancy loss. We further show that the variant associates with more random placement of crossovers and lower recombination rate in longer chromosomes but higher in the shorter ones. These results support the hypothesis that some pregnancy losses are due to failures in recombination. They further demonstrate that variants with a substantial effect on the quality of recombination can be maintained in the population.

Loss-of-function variants in ITSN1 confer high risk of Parkinson's disease.

Parkinson's disease (PD) is a debilitating neurodegenerative disorder and its rising global incidence highlights the need for the identification of modifiable risk factors. In a gene-based burden test of rare variants (8647 PD cases and 777,693 controls) we discovered a novel association between loss-of-function variants in ITSN1 and PD. This association was further supported with burden data from the Neurodegenerative Disease Knowledge Portal and the Accelerating Medicines Partnership Parkinson's Disease Knowledge Platform. Our findings show that Rho GTPases and disruptions in synaptic vesicle transport may be involved in the pathogenesis of PD, pointing to the possibility of novel therapeutic approaches.

The correlation between CpG methylation and gene expression is driven by sequence variants.

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Start codon variant in LAG3 is associated with decreased LAG-3 expression and increased risk of autoimmune thyroid disease.

Autoimmune thyroid disease (AITD) is a common autoimmune disease. In a GWAS meta-analysis of 110,945 cases and 1,084,290 controls, 290 sequence variants at 225 loci are associated with AITD. Of these variants, 115 are previously unreported. Multiomics analysis yields 235 candidate genes outside the MHC-region and the findings highlight the importance of genes involved in T-cell regulation. A rare 5'-UTR variant (rs781745126-T, MAF = 0.13% in Iceland) in LAG3 has the largest effect (OR = 3.42, P = 2.2 × 10-16) and generates a novel start codon for an open reading frame upstream of the canonical protein translation initiation site. rs781745126-T reduces mRNA and surface expression of the inhibitory immune checkpoint LAG-3 co-receptor on activated lymphocyte subsets and halves LAG-3 levels in plasma among heterozygotes. All three homozygous carriers of rs781745126-T have AITD, of whom one also has two other T-cell mediated diseases, that is vitiligo and type 1 diabetes. rs781745126-T associates nominally with vitiligo (OR = 5.1, P = 6.5 × 10-3) but not with type 1 diabetes. Thus, the effect of rs781745126-T is akin to drugs that inhibit LAG-3, which unleash immune responses and can have thyroid dysfunction and vitiligo as adverse events. This illustrates how a multiomics approach can reveal potential drug targets and safety concerns.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

POPDC2 encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss of Popdc2 in mice results in sinus pauses and bradycardia and morpholino knockdown of popdc2 in zebrafish results in atrioventricular (AV) block. We identified bi-allelic variants in POPDC2 in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identified POPDC2 variants are predicted to diminish the ability of POPDC2 to bind cAMP. In in vitro electrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1-POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants in POPDC1 and POPDC2 and presence of sinus node disease in POPDC2, but not in POPDC1 related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence for POPDC2 as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functional POPDC2 display sinus and AV node dysfunction.

Integrative common and rare variant analyses provide insights into the genetic architecture of liver cirrhosis.

We report a multi-ancestry genome-wide association study on liver cirrhosis and its associated endophenotypes, alanine aminotransferase (ALT) and γ-glutamyl transferase. Using data from 12 cohorts, including 18,265 cases with cirrhosis, 1,782,047 controls, up to 1 million individuals with liver function tests and a validation cohort of 21,689 cases and 617,729 controls, we identify and validate 14 risk associations for cirrhosis. Many variants are located near genes involved in hepatic lipid metabolism. One of these, PNPLA3 p.Ile148Met, interacts with alcohol intake, obesity and diabetes on the risk of cirrhosis and hepatocellular carcinoma (HCC). We develop a polygenic risk score that associates with the progression from cirrhosis to HCC. By focusing on prioritized genes from common variant analyses, we find that rare coding variants in GPAM associate with lower ALT, supporting GPAM as a potential target for therapeutic inhibition. In conclusion, this study provides insights into the genetic underpinnings of cirrhosis.

Genome-wide analyses identify 21 infertility loci and over 400 reproductive hormone loci across the allele frequency spectrum.

Genome-wide association studies (GWASs) may help inform treatments for infertility, whose causes remain unknown in many cases. Here we present GWAS meta-analyses across six cohorts for male and female infertility in up to 41,200 cases and 687,005 controls. We identified 21 genetic risk loci for infertility (P≤5E-08), of which 12 have not been reported for any reproductive condition. We found positive genetic correlations between endometriosis and all-cause female infertility (rg=0.585, P=8.98E-14), and between polycystic ovary syndrome and anovulatory infertility (rg=0.403, P=2.16E-03). The evolutionary persistence of female infertility-risk alleles in EBAG9 may be explained by recent directional selection. We additionally identified up to 269 genetic loci associated with follicle-stimulating hormone (FSH), luteinising hormone, oestradiol, and testosterone through sex-specific GWAS meta-analyses (N=6,095-246,862). While hormone-associated variants near FSHB and ARL14EP colocalised with signals for anovulatory infertility, we found no rg between female infertility and reproductive hormones (P>0.05). Exome sequencing analyses in the UK Biobank (N=197,340) revealed that women carrying testosterone-lowering rare variants in GPC2 were at higher risk of infertility (OR=2.63, P=1.25E-03). Taken together, our results suggest that while individual genes associated with hormone regulation may be relevant for fertility, there is limited genetic evidence for correlation between reproductive hormones and infertility at the population level. We provide the first comprehensive view of the genetic architecture of infertility across multiple diagnostic criteria in men and women, and characterise its relationship to other health conditions.

Variants at the Interleukin 1 Gene Locus and Pericarditis.

Importance: Recurrent pericarditis is a treatment challenge and often a debilitating condition. Drugs inhibiting interleukin 1 cytokines are a promising new treatment option, but their use is based on scarce biological evidence and clinical trials of modest sizes, and the contributions of innate and adaptive immune processes to the pathophysiology are incompletely understood. Objective: To use human genomics, transcriptomics, and proteomics to shed light on the pathogenesis of pericarditis. Design, Setting, and Participants: This was a meta-analysis of genome-wide association studies of pericarditis from 5 countries. Associations were examined between the pericarditis-associated variants and pericarditis subtypes (including recurrent pericarditis) and secondary phenotypes. To explore mechanisms, associations with messenger RNA expression (cis-eQTL), plasma protein levels (pQTL), and CpG methylation of DNA (ASM-QTL) were assessed. Data from Iceland (deCODE genetics, 1983-2020), Denmark (Copenhagen Hospital Biobank/Danish Blood Donor Study, 1977-2022), the UK (UK Biobank, 1953-2021), the US (Intermountain, 1996-2022), and Finland (FinnGen, 1970-2022) were included. Data were analyzed from September 2022 to August 2023. Exposure: Genotype. Main Outcomes and Measures: Pericarditis. Results: In this genome-wide association study of 4894 individuals with pericarditis (mean [SD] age at diagnosis, 51.4 [17.9] years, 2734 [67.6%] male, excluding the FinnGen cohort), associations were identified with 2 independent common intergenic variants at the interleukin 1 locus on chromosome 2q14. The lead variant was rs12992780 (T) (effect allele frequency [EAF], 31%-40%; odds ratio [OR], 0.83; 95% CI, 0.79-0.87; P = 6.67 × 10-16), downstream of IL1B and the secondary variant rs7575402 (A or T) (EAF, 45%-55%; adjusted OR, 0.89; 95% CI, 0.85-0.93; adjusted P = 9.6 × 10-8). The lead variant rs12992780 had a smaller odds ratio for recurrent pericarditis (0.76) than the acute form (0.86) (P for heterogeneity = .03) and rs7575402 was associated with CpG methylation overlapping binding sites of 4 transcription factors known to regulate interleukin 1 production: PU.1 (encoded by SPI1), STAT1, STAT3, and CCAAT/enhancer-binding protein ß (encoded by CEBPB). Conclusions and Relevance: This study found an association between pericarditis and 2 independent sequence variants at the interleukin 1 gene locus. This finding has the potential to contribute to development of more targeted and personalized therapy of pericarditis with interleukin 1-blocking drugs.