Search

Mitochondrial DNA copy number (mtDNA-CN) measured from blood specimens is a minimally invasive marker of mitochondrial function that exhibits both inter-individual and intercellular variation. To identify genes involved in regulating mitochondrial function, we performed a genome-wide association study (GWAS) in 465,809 White individuals from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium and the UK Biobank (UKB). We identified 133 SNPs with statistically significant, independent effects associated with mtDNA-CN across 100 loci. A combination of fine-mapping, variant annotation, and co-localization analyses was used to prioritize genes within each of the 133 independent sites. Putative causal genes were enriched for known mitochondrial DNA depletion syndromes (p = 3.09 × 10-15) and the gene ontology (GO) terms for mtDNA metabolism (p = 1.43 × 10-8) and mtDNA replication (p = 1.2 × 10-7). A clustering approach leveraged pleiotropy between mtDNA-CN associated SNPs and 41 mtDNA-CN associated phenotypes to identify functional domains, revealing three distinct groups, including platelet activation, megakaryocyte proliferation, and mtDNA metabolism. Finally, using mitochondrial SNPs, we establish causal relationships between mitochondrial function and a variety of blood cell-related traits, kidney function, liver function and overall (p = 0.044) and non-cancer mortality (p = 6.56 × 10-4).

Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10 percentage points higher in African populations. Three signals (SERPINA1, ZFP36L2, and TLR10) contain predicted deleterious missense variants. Two loci, SOCS3 and HPN, each harbor two conditionally distinct, non-coding variants. The gene region encoding the fibrinogen protein chain subunits (FGG;FGB;FGA), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common in African ancestry populations but extremely rare in Europeans (MAFAFR=0.180; MAFEUR=0.008). Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation., (Copyright © 2024 American Society of Hematology.)

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium genome-wide association studies meta-analyses of European- (71 771 cases and 1 059 740 controls) and African-ancestry samples (7482 cases and 129 975 controls). We used LDpred2 and PRS-CSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6781 cases and 103 016 controls) and African-ancestry sample (1385 cases and 12 569 controls). Multi-ancestry PRSs with weights tuned in European-ancestry samples slightly outperformed ancestry-specific PRSs in European-ancestry test samples (e.g. the area under the receiver operating curve [AUC] was 0.609 for PRS-CSx_combinedEUR and 0.608 for PRS-CSxEUR [P = 0.00029]). Multi-ancestry PRSs with weights tuned in African-ancestry samples also outperformed ancestry-specific PRSs in African-ancestry test samples (PRS-CSxAFR: AUC = 0.58, PRS-CSx_combined AFR: AUC = 0.59), although this difference was not statistically significant (P = 0.34). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS might be used to improve performance across diverse populations to identify individuals at highest risk for VTE., (© The Author(s) 2024. Published by Oxford University Press.)

African American (AA) kidney transplant recipients (KTRs) have poor outcomes, which may in-part be due to tacrolimus (TAC) sub-optimal immunosuppression. We previously determined the common genetic regulators of TAC pharmacokinetics in AAs which were CYP3A5 *3, *6, and *7. To identify low-frequency variants that impact TAC pharmacokinetics, we used extreme phenotype sampling and compared individuals with extreme high (n = 58) and low (n = 60) TAC troughs (N = 515 AA KTRs). Targeted next generation sequencing was conducted in these two groups. Median TAC troughs in the high group were 7.7 ng/ml compared with 6.3 ng/ml in the low group, despite lower daily doses of 5 versus 12 mg, respectively. Of 34,542 identified variants across 99 genes, 1406 variants were suggestively associated with TAC troughs in univariate models (p-value < 0.05), however none were significant after multiple testing correction. We suggest future studies investigate additional sources of TAC pharmacokinetic variability such as drug-drug-gene interactions and pharmacomicrobiome., (© 2024. The Author(s), under exclusive licence to Springer Nature Limited.)

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., (© The Author(s) 2024. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Genome-wide association studies (GWAS) have become well-powered to detect loci associated with telomere length. However, no prior work has validated genes nominated by GWAS to examine their role in telomere length regulation. We conducted a multi-ancestry meta-analysis of 211,369 individuals and identified five novel association signals. Enrichment analyses of chromatin state and cell-type heritability suggested that blood/immune cells are the most relevant cell type to examine telomere length association signals. We validated specific GWAS associations by overexpressing KBTBD6 or POP5 and demonstrated that both lengthened telomeres. CRISPR/Cas9 deletion of the predicted causal regions in K562 blood cells reduced expression of these genes, demonstrating that these loci are related to transcriptional regulation of KBTBD6 and POP5. Our results demonstrate the utility of telomere length GWAS in the identification of telomere length regulation mechanisms and validate KBTBD6 and POP5 as genes affecting telomere length regulation., (© 2024. The Author(s).)

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 (R 2  = 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., (© 2024. The Author(s).)

Abstract: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Background and Objectives: Nearly all genetic analyses of Parkinson disease (PD) have been in populations of European ancestry. We sought to test the ability of a machine learning method to extract accurate PD diagnoses from an electronic medical record (EMR) system, to see whether genetic variants identified in European populations generalize to individuals of African and Hispanic ancestries, and to compare the rates of PD across ancestries., Methods: A machine learning method using natural language processing was applied to EMRs of US veterans participating in the VA Million Veteran Program (MVP) to identify individuals with PD. These putative cases were vetted via blind chart review by a movement disorder specialist. A polygenic risk score (PRS) of 90 established genetic variants whose genotypes were imputed from a customized Axiom Biobank Array was evaluated in different case groups., Results: The EMR prediction scores had a distinct trimodal distribution, with 97% of the high group and only 30% of the middle group having a credible diagnosis of PD. Using the 3,542 cases from the high group matched 4:1 to controls, the PRS was highly predictive in individuals of European ancestry (n = 3,137 cases; OR = 1.82; p = 8.01E-48), and nearly identical effect sizes were seen in individuals of African (n = 184; OR = 2.07; p = 3.4E-4) and Hispanic ancestries (n = 221; OR = 2.13; p = 3.9E-6). The PRS was much less predictive for the 2,757 European ancestry cases who had an ICD code for PD but for whom the machine learning method had a lower confidence in their diagnosis. No novel ancestry-specific genetic variants were identified. Individuals with African ancestry had one-quarter the rate of PD compared with European or Hispanic ancestries aged 60-70 years and one half the rate in the 70-80 years age range. African American cases had a higher proportion of their DNA originating in Europe compared with African American controls., Discussion: Machine learning can reliably classify PD using data from a large EMR. Larger studies of non-European populations are required to confirm the generalizability of PD risk variants identified in populations of European ancestry and the increased risk coming from a higher proportion of European DNA in African Americans., Competing Interests: The authors report no relevant disclosures. Go to Neurology.org/NG for full disclosures., (Copyright © 2023 The Author(s). Published by Wolters Kluwer Health, Inc. on behalf of the American Academy of Neurology.)

Abdominal aortic aneurysm (AAA) is a common disease with substantial heritability. In this study, we performed a genome-wide association meta-analysis from 14 discovery cohorts and uncovered 141 independent associations, including 97 previously unreported loci. A polygenic risk score derived from meta-analysis explained AAA risk beyond clinical risk factors. Genes at AAA risk loci indicate involvement of lipid metabolism, vascular development and remodeling, extracellular matrix dysregulation and inflammation as key mechanisms in AAA pathogenesis. These genes also indicate overlap between the development of AAA and other monogenic aortopathies, particularly via transforming growth factor β signaling. Motivated by the strong evidence for the role of lipid metabolism in AAA, we used Mendelian randomization to establish the central role of nonhigh-density lipoprotein cholesterol in AAA and identified the opportunity for repurposing of proprotein convertase, subtilisin/kexin-type 9 (PCSK9) inhibitors. This was supported by a study demonstrating that PCSK9 loss of function prevented the development of AAA in a preclinical mouse model., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Megabase-scale mosaic chromosomal alterations (mCAs) in blood are prognostic markers for a host of human diseases. Here, to gain a better understanding of mCA rates in genetically diverse populations, we analyzed whole-genome sequencing data from 67,390 individuals from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine program. We observed higher sensitivity with whole-genome sequencing data, compared with array-based data, in uncovering mCAs at low mutant cell fractions and found that individuals of European ancestry have the highest rates of autosomal mCAs and the lowest rates of chromosome X mCAs, compared with individuals of African or Hispanic ancestry. Although further studies in diverse populations will be needed to replicate our findings, we report three loci associated with loss of chromosome X, associations between autosomal mCAs and rare variants in DCPS, ADM17, PPP1R16B and TET2 and ancestry-specific variants in ATM and MPL with mCAs in cis., (© 2023. The Author(s).)

Background The relationship between mitochondrial DNA copy number (mtDNA CN) and cardiovascular disease remains elusive. Methods and Results We performed cross-sectional and prospective association analyses of blood-derived mtDNA CN and cardiovascular disease outcomes in 27 316 participants in 8 cohorts of multiple racial and ethnic groups with whole-genome sequencing. We also performed Mendelian randomization to explore causal relationships of mtDNA CN with coronary heart disease (CHD) and cardiometabolic risk factors (obesity, diabetes, hypertension, and hyperlipidemia). P <0.01 was used for significance. We validated most of the previously reported associations between mtDNA CN and cardiovascular disease outcomes. For example, 1-SD unit lower level of mtDNA CN was associated with 1.08 (95% CI, 1.04-1.12; P <0.001) times the hazard for developing incident CHD, adjusting for covariates. Mendelian randomization analyses showed no causal effect from a lower level of mtDNA CN to a higher CHD risk (β=0.091; P =0.11) or in the reverse direction (β=-0.012; P =0.076). Additional bidirectional Mendelian randomization analyses revealed that low-density lipoprotein cholesterol had a causal effect on mtDNA CN (β=-0.084; P <0.001), but the reverse direction was not significant ( P =0.059). No causal associations were observed between mtDNA CN and obesity, diabetes, and hypertension, in either direction. Multivariable Mendelian randomization analyses showed no causal effect of CHD on mtDNA CN, controlling for low-density lipoprotein cholesterol level ( P =0.52), whereas there was a strong direct causal effect of higher low-density lipoprotein cholesterol on lower mtDNA CN, adjusting for CHD status (β=-0.092; P <0.001). Conclusions Our findings indicate that high low-density lipoprotein cholesterol may underlie the complex relationships between mtDNA CN and vascular atherosclerosis.

Psychosis is a severe mental disorder characterized by abnormal thoughts and perceptions (e.g., hallucinations) occurring quintessentially in schizophrenia and in several other neuropsychiatric disorders. Schizophrenia is widely considered as a neurodevelopmental disorder that onsets during teenage/early adulthood. A multiplex consanguineous Pakistani family was afflicted with severe psychosis and apparent autosomal recessive transmission. The first-cousin parents and five children were healthy, whereas two teenage daughters were severely affected. Structured interviews confirmed the diagnosis of DSM-V schizophrenia. Probands and father underwent next-generation sequencing. All available relatives were subjected to confirmatory Sanger sequencing. Homozygosity mapping and directed a priori filtering identified only one rare variant [MAF < 5(10) -5 ] at a residue conserved across vertebrates. The variant was a non-catalytic deubiquitinase, USP53 (p.Cys228Arg), predicted in silico as damaging. Genome sequencing did not identify any other potentially pathogenic single nucleotide variant or structural variant. Since the literature on USP53 lacked relevance to mental illness or CNS expression, studies were conducted which revealed USP53 localization in regions of the hippocampus (CA 1-3) and granular dentate. The staining pattern was like that seen with GRIA2/GluA2 and GRIP2 antibodies. All three proteins coimmunoprecipitated. These findings support the glutamate hypothesis of schizophrenia as part of the AMPA-R interactome. If confirmed, USP53 appears to be one of the few Mendelian variants potentially causal to a common-appearing mental disorder that is a rare genetic form of schizophrenia.

Maintenance of telomere length has long been established to play a role in the biology of cancer and several studies suggest that it may be especially important in myeloid malignancies. To overcome potential bias in confounding and reverse causation of observational studies, we use both a polygenic risk score (PRS) and inverse-variance weighted (IVW) Mendelian randomization (MR) analyses to estimate the relationship between genetically predicted leukocyte telomere length (LTL) and acute myeloid leukemia (AML) risk in 498 cases and 2099 controls and myelodysplastic syndrome (MDS) risk in 610 cases and 1759 controls. Genetic instruments derived from four recent studies explaining 1.23-4.57% of telomere variability were considered. We used multivariable logistic regression to estimate odds ratios (OR, 95% confidence intervals [CI]) as the measure of association between individual single-nucleotide polymorphisms and myeloid malignancies. We observed a significant association between a PRS of longer predicted LTL and AML using three genetic instruments (OR = 4.03 per ~1200 base pair [bp] increase in LTL, 95% CI: 1.65, 9.85 using Codd et al. [Codd, V., Nelson, C.P., Albrecht, E., Mangino, M., Deelen, J., Buxton, J.L., Hottenga, J.J., Fischer, K., Esko, T., Surakka, I. et al. (2013) Identification of seven loci affecting mean telomere length and their association with disease. Nat. Genet., 45, 422-427 427e421-422.], OR = 3.48 per one-standard deviation increase in LTL, 95% CI: 1.74, 6.97 using Li et al. [Li, C., Stoma, S., Lotta, L.A., Warner, S., Albrecht, E., Allione, A., Arp, P.P., Broer, L., Buxton, J.L., Alves, A.D.S.C. et al. (2020) Genome-wide association analysis in humans links nucleotide metabolism to leukocyte telomere length. Am. J. Hum. Genet., 106, 389-404.] and OR = 2.59 per 1000 bp increase in LTL, 95% CI: 1.03, 6.52 using Taub et al. [Taub, M.A., Conomos, M.P., Keener, R., Iyer, K.R., Weinstock, J.S., Yanek, L.R., Lane, J., Miller-Fleming, T.W., Brody, J.A., Raffield, L.M. et al. (2022) Genetic determinants of telomere length from 109,122 ancestrally diverse whole-genome sequences in TOPMed. Cell Genom., 2.] genetic instruments). MR analyses further indicated an association between LTL and AML risk (PIVW ≤ 0.049) but not MDS (all PIVW ≥ 0.076). Findings suggest variation in genes relevant to telomere function and maintenance may be important in the etiology of AML but not MDS., (© The Author(s) 2023. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Background: Annular pancreas (AP) is a congenital defect of unknown cause in which the pancreas encircles the duodenum. Theories include abnormal migration and rotation of the ventral bud, persistence of ectopic pancreatic tissue, and inappropriate fusion of the ventral and dorsal buds before rotation. The few reported familial cases suggest a genetic contribution., Methods: We conducted exome sequencing in 115 affected infants from the California birth defects registry., Results: Seven cases had a single heterozygous missense variant in IQGAP1, five of them with CADD scores >20; seven other infants had a single heterozygous missense variant in NRCAM, five of them with CADD scores >20. We also looked at genes previously associated with AP and found two rare heterozygous missense variants, one each in PDX1 and FOXF1., Conclusion: IQGAP1 and NRCAM are crucial in cell polarization and migration. Mutations result in decreased motility which could possibly cause the ventral bud to not migrate normally. To our knowledge, this is the first study reporting a possible association for IQGAP1 and NRCAM with AP. Our findings of rare genetic variants involved in cell migration in 15% of our population raise the possibility that AP may be related to abnormal cell migration., (© 2023 The Authors. Molecular Genetics & Genomic Medicine published by Wiley Periodicals LLC.)

Mitochondria carry their own circular genome and disruption of the mitochondrial genome is associated with various aging-related diseases. Unlike the nuclear genome, mitochondrial DNA (mtDNA) can be present at 1000 s to 10,000 s copies in somatic cells and variants may exist in a state of heteroplasmy, where only a fraction of the DNA molecules harbors a particular variant. We quantify mtDNA heteroplasmy in 194,871 participants in the UK Biobank and find that heteroplasmy is associated with a 1.5-fold increased risk of all-cause mortality. Additionally, we functionally characterize mtDNA single nucleotide variants (SNVs) using a constraint-based score, mitochondrial local constraint score sum (MSS) and find it associated with all-cause mortality, and with the prevalence and incidence of cancer and cancer-related mortality, particularly leukemia. These results indicate that mitochondria may have a functional role in certain cancers, and mitochondrial heteroplasmic SNVs may serve as a prognostic marker for cancer, especially for leukemia., (© 2023. Springer Nature Limited.)

Background: The MHC class I chain-related protein A (MICA) and protein B (MICB) participate in tumor immunosurveillance and may be important in colorectal cancer, but have not been examined in colorectal cancer development., Methods: sMICA and sMICB blood levels were measured by SomaScan in Visit 2 (1990-92, baseline) and Visit 3 (1993-95) samples in cancer-free participants in the Atherosclerosis Risk in Communities Study. We selected rs1051792, rs1063635, rs2516448, rs3763288, rs1131896, rs2596542, and rs2395029 that were located in or in the vicinity of MICA or MICB and were associated with cancer or autoimmune diseases in published studies. SNPs were genotyped by the Affymetrix Genome-Wide Human SNP Array. We applied linear and Cox proportional hazards regressions to examine the associations of preselected SNPs with sMICA and sMICB levels and colorectal cancer risk (236 colorectal cancers, 8,609 participants) and of sMICA and sMICB levels with colorectal cancer risk (312 colorectal cancers, 10,834 participants). In genetic analyses, estimates adjusted for ancestry markers were meta-analyzed., Results: Rs1051792-A, rs1063635-A, rs2516448-C, rs3763288-A, rs2596542-T, and rs2395029-G were significantly associated with decreased sMICA levels. Rs2395029-G, in the vicinity of MICA and MICB, was also associated with increased sMICB levels. Rs2596542-T was significantly associated with decreased colorectal cancer risk. Lower sMICA levels were associated with lower colorectal cancer risk in males (HR = 0.68; 95% confidence interval, 0.49-0.96) but not in females (Pinteraction = 0.08)., Conclusions: Rs2596542-T associated with lower sMICA levels was associated with decreased colorectal cancer risk. Lower sMICA levels were associated with lower colorectal cancer risk in males., Impact: These findings support an importance of immunosurveillance in colorectal cancer., (©2023 American Association for Cancer Research.)

Successful treatment of endogenous Cushing disease (CD) is often followed by a period of adrenal insufficiency (AI). We performed an exploratory study on genetic factors potentially involved in the hypothalamic-pituitary-adrenal (HPA) axis recovery in patients with CD after remission. We identified 90 patients who achieved remission after surgery and had a minimum of 3 months follow-up. Variants in a selected panel of genes that were rare in the general population and predicted as damaging in silico were retrieved from whole exome sequencing analysis. We did not identify any variant with significant correlation with recovery time after adjusting for multiple comparisons. On gene-specific analysis the BAG1 gene showed a correlation with shorter duration of postsurgical AI, but both patients with BAG1 variants later experienced a recurrence. After excluding patients with recurrence, no statistical association was recorded. To conclude, we did not identify a strong genetic modifier of HPA recovery in this exploratory study., (Published by Oxford University Press on behalf of the Endocrine Society 2023.)

Background: Risk for venous thromboembolism has a strong genetic component. Whole genome sequencing from the TOPMed program (Trans-Omics for Precision Medicine) allowed us to look for new associations, particularly rare variants missed by standard genome-wide association studies., Methods: The 3793 cases and 7834 controls (11.6% of cases were individuals of African, Hispanic/Latino, or Asian ancestry) were analyzed using a single variant approach and an aggregate gene-based approach using our primary filter (included only loss-of-function and missense variants predicted to be deleterious) and our secondary filter (included all missense variants)., Results: Single variant analyses identified associations at 5 known loci. Aggregate gene-based analyses identified only PROC (odds ratio, 6.2 for carriers of rare variants; P =7.4×10 -14 ) when using our primary filter. Employing our secondary variant filter led to a smaller effect size at PROC (odds ratio, 3.8; P =1.6×10 -14 ), while excluding variants found only in rare isoforms led to a larger one (odds ratio, 7.5). Different filtering strategies improved the signal for 2 other known genes: PROS1 became significant (minimum P =1.8×10 -6 with the secondary filter), while SERPINC1 did not (minimum P =4.4×10 -5 with minor allele frequency <0.0005). Results were largely the same when restricting the analyses to include only unprovoked cases; however, one novel gene, MS4A1 , became significant ( P =4.4×10 -7 using all missense variants with minor allele frequency <0.0005)., Conclusions: Here, we have demonstrated the importance of using multiple variant filtering strategies, as we detected additional genes when filtering variants based on their predicted deleteriousness, frequency, and presence on the most expressed isoforms. Our primary analyses did not identify new candidate loci; thus larger follow-up studies are needed to replicate the novel MS4A1 locus and to identify additional rare variation associated with venous thromboembolism.

Introduction: Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects., Methods: Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs., Results: CNVs in schizophrenia cases were twice as large on average than CNVs in controls ( p  = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions., Discussion: Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Abumadini, Al Ghamdi, Alqahtani, Almedallah, Callans, Jarad, Cyrus, Koeleman, Keating, Pankratz and Al-Ali.)

Individuals with cystic fibrosis (CF), caused by biallelic germline mutations in the cystic fibrosis transmembrane conductance regulator ( CFTR) , have higher risk and earlier onset of colorectal cancer (CRC). A subset of CRC patients in the non-CF population expresses low levels of tumoral CFTR mRNA which may also cause decreased CFTR activity. To determine the consequences of reduced CFTR expression in this population, we investigated association of tumoral CFTR expression with overall and disease-specific mortality in CRC patients. CFTR mRNA expression, clinical factors and survival data from 1177 CRC patients reported in the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus studies GSE39582 and GSE17538 were included. Log-transformed and z-normalized [mean = 0, standard deviation (SD) = 1] CFTR expression values were modeled as quartiles or dichotomized at the median. Univariate and multivariable Cox proportional hazards regression models were used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for overall and disease-specific mortality in individual studies and meta-analyses. Analyses of each of the three individual datasets showed a robust association of decreased CFTR expression with increased mortality. In meta-analyses adjusted for stage at diagnosis, age and sex, CFTR expression was inversely associated with risk of overall death [pooled HR (95% CI): 0.70 (0.57-0.86)] and disease-specific death [pooled HR (95% CI): 0.68 (0.47-0.99)]. Associations did not differ by stage at diagnosis, age, or sex. Meta-analysis of overall death stratified by microsatellite instable (MSI) versus microsatellite stable (MSS) status indicated potential interaction between MSI/MSS status and CFTR expression, ( p -interaction: 0.06). The findings from these three datasets support the hypothesis that low CFTR expression is associated with increased CRC mortality.

Background: Genome-wide association studies have reported 23 gene loci related to abdominal aortic aneurysm (AAA)-a potentially lethal condition characterized by a weakened dilated vessel wall. This study aimed to identify proteomic signatures and pathways related to these risk loci to better characterize AAA genetic susceptibility., Methods: Plasma concentrations of 4870 proteins were determined using a DNA aptamer-based array. Linear regression analysis estimated the associations between the 23 risk alleles and plasma protein levels with adjustments for potential confounders in a race-stratified analysis of 1671 Black and 7241 White participants. Significant proteins were then evaluated for their prediction of clinical AAA (454 AAA events in 11 064 individuals), and those significantly associated with AAA were further interrogated using Mendelian randomization analysis., Results: Risk variants proximal to PSRC1-CELSR2-SORT1 , PCIF1-ZNF335-MMP9, RP11-136O12.2/TRIB1 , ZNF259/APOA5, IL6R , PCSK9 , LPA , and APOE were associated with 118 plasma proteins in Whites and 59 were replicated in Black participants. Novel associations with clinical AAA incidence were observed for kit ligand (HR, 0.59 [95% CI, 0.42-0.82] for top versus first quintiles) and neogenin (HR, 0.64 [95% CI, 0.46-0.88]) over a median 21.2-year follow-up; neogenin was also associated with ultrasound-detected asymptomatic AAA (N=4295; 57 asymptomatic AAA cases). Mendelian randomization inverse variance weighted estimates suggested that AAA risk is promoted by lower levels of kit ligand (OR per SD=0.67; P =1.4×10 -5 ) and neogenin (OR per SD=0.50; P =0.03)., Conclusions: Low levels of neogenin and kit ligand may be novel risk factors for AAA development in potentially causal pathways. These findings provide insights and potential targets to reduce AAA susceptibility.

Introduction: Most strategies for prevention of venous thromboembolism focus on preventing recurrent events. Yet, primary prevention might be possible through approaches targeting the whole population or high-risk patients. To inform possible prevention strategies, population-based information on the ability of genetic risk scores to identify risk of incident venous thromboembolism is needed., Materials and Methods: We used proportional hazards regression to relate two published genetic risk scores (273-variants versus 5-variants) with venous thromboembolism incidence in the Atherosclerosis Risk in Communities Study (ARIC) cohort (n = 11,292), aged 45-64 at baseline, drawn from 4 US communities., Results: Over a median of 28 years, ARIC identified 788 incident venous thromboembolism events. Incidence rates rose more than two-fold across quartiles of the 273-variant genetic risk score: 1.7, 2.7, 3.4 and 4.0 per 1,000 person-years. For White participants, age, sex, and ancestry-adjusted hazard ratios (95% confidence intervals) across quartiles were strong [1 (reference), 1.30 (0.99,1.70), 1.85 (1.43,2.40), and 2.58 (2.04,3.28)] but weaker for Black participants [1, 1.05 (0.63,1.75), 1.37 (0.84,2.22), and 1.32 (0.80,2.20)]. The 5-variant genetic risk score showed a less steep gradient, with hazard ratios in Whites of 1, 1.17 (0.89,1.54), 1.48 (1.14,1.92), and 2.18 (1.71,2.79). Models including the 273-variant genetic risk score plus lifestyle and clinical factors had a c-statistic of 0.67., Conclusions: In the general population, middle-aged adults in the highest quartile of either genetic risk score studied have approximately two-fold higher risk of an incident venous thromboembolism compared with the lowest quartile. The genetic risk scores show a weaker association with venous thromboembolism for Black people., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Folsom et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Aims/hypothesis: Genetic predisposition to type 2 diabetes is well-established, and genetic risk scores (GRS) have been developed that capture heritable liabilities for type 2 diabetes phenotypes. However, the proteins through which these genetic variants influence risk have not been thoroughly investigated. This study aimed to identify proteins and pathways through which type 2 diabetes risk variants may influence pathophysiology., Methods: Using a proteomics data-driven approach in a discovery sample of 7241 White participants in the Atherosclerosis Risk in Communities Study (ARIC) cohort and a replication sample of 1674 Black ARIC participants, we interrogated plasma levels of 4870 proteins and four GRS of specific type 2 diabetes phenotypes related to beta cell function, insulin resistance, lipodystrophy, BMI/blood lipid abnormalities and a composite score of all variants combined., Results: Twenty-two plasma proteins were identified in White participants after Bonferroni correction. Of the 22 protein-GRS associations that were statistically significant, 10 were replicated in Black participants and all but one were directionally consistent. In a secondary analysis, 18 of the 22 proteins were found to be associated with prevalent type 2 diabetes and ten proteins were associated with incident type 2 diabetes. Two-sample Mendelian randomisation indicated that complement C2 may be causally related to greater type 2 diabetes risk (inverse variance weighted estimate: OR 1.65 per SD; p=7.0 × 10 -3 ), while neuropilin-2 was inversely associated (OR 0.44 per SD; p=8.0 × 10 -3 )., Conclusions/interpretation: Identified proteins may represent viable intervention or pharmacological targets to prevent, reverse or slow type 2 diabetes progression, and further research is needed to pursue these targets., (© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Purpose: Cushing's disease (CD) is often explained by a single somatic sequence change. Germline defects, however, often go unrecognized. We aimed to determine the frequency and associated phenotypes of genetic drivers of CD in a large cohort., Methods: We studied 245 unrelated patients with CD (139 female, 56.7%), including 230 (93.9%) pediatric and 15 (6.1%) adult patients. Germline exome sequencing was performed in 184 patients; tumor exome sequencing was also done in 27 of them. A total of 43 germline samples and 92 tumor samples underwent Sanger sequencing of specific genes. Rare variants of uncertain significance, likely pathogenic (LP), or pathogenic variants in CD-associated genes, were identified., Results: Germline variants (13 variants of uncertain significance, 8 LP, and 11 pathogenic) were found in 8 of 19 patients (42.1%) with positive family history and in 23 of 226 sporadic patients (10.2%). Somatic variants (1 LP and 7 pathogenic) were found in 20 of 119 tested individuals (16.8%); one of them had a coexistent germline defect. Altogether, variants of interest were identified at the germline level in 12.2% of patients, at the somatic level in 7.8%, and coexisting germline and somatic variants in 0.4%, accounting for one-fifth of the cohort., Conclusion: We report an estimate of the contribution of multiple germline and somatic genetic defects underlying CD in a single cohort., Competing Interests: Conflict of Interest The authors declare no conflict of interest., (Copyright © 2022 American College of Medical Genetics and Genomics. All rights reserved.)

Background: Although micronutrients modulate immunity and inflammation, it remains elusive whether they are implicated in the development and progression of chronic pancreatitis (CP). This study aimed to investigate differences in the circulating levels of selected carotenoids and vitamins between CP and controls and trends in the levels of these micronutrients across controls, early CP, and definite CP., Methods: Demographic and lifestyle data were extracted from medical records for 53 patients with CP (13 early and 38 definite) and obtained using a questionnaire for 52 controls. Plasma β-carotene, lycopene, cryptoxanthin, zeaxanthin, and α-tocopherol and serum 25(OH)D, folate, IL-6, TNF-α, and MCP-1 were measured with state-of-the-art methods., Results: The levels of all micronutrients (except folate) were significantly lower in CP than in controls. There was a progressive decrease in the levels of these micronutrients across controls, early CP, and definite CP (all p values for trend: ≤0.0012); e.g., plasma lycopene was 36.6, 21.5, and 14.5 μg/dL for controls, early CP, and definite CP, respectively. After adjustment for confounders, there were strong, inverse associations between the levels of all micronutrients (except folate) and CP (e.g., OR (95% CI) for ≥ median vs.

Background: Venous thromboembolism (VTE) is a life-threatening vascular event with environmental and genetic determinants. Recent VTE genome-wide association studies (GWAS) meta-analyses involved nearly 30 000 VTE cases and identified up to 40 genetic loci associated with VTE risk, including loci not previously suspected to play a role in hemostasis. The aim of our research was to expand discovery of new genetic loci associated with VTE by using cross-ancestry genomic resources., Methods: We present new cross-ancestry meta-analyzed GWAS results involving up to 81 669 VTE cases from 30 studies, with replication of novel loci in independent populations and loci characterization through in silico genomic interrogations., Results: In our genetic discovery effort that included 55 330 participants with VTE (47 822 European, 6320 African, and 1188 Hispanic ancestry), we identified 48 novel associations, of which 34 were replicated after correction for multiple testing. In our combined discovery-replication analysis (81 669 VTE participants) and ancestry-stratified meta-analyses (European, African, and Hispanic), we identified another 44 novel associations, which are new candidate VTE-associated loci requiring replication. In total, across all GWAS meta-analyses, we identified 135 independent genomic loci significantly associated with VTE risk. A genetic risk score of the significantly associated loci in Europeans identified a 6-fold increase in risk for those in the top 1% of scores compared with those with average scores. We also identified 31 novel transcript associations in transcriptome-wide association studies and 8 novel candidate genes with protein quantitative-trait locus Mendelian randomization analyses. In silico interrogations of hemostasis and hematology traits and a large phenome-wide association analysis of the 135 GWAS loci provided insights to biological pathways contributing to VTE, with some loci contributing to VTE through well-characterized coagulation pathways and others providing new data on the role of hematology traits, particularly platelet function. Many of the replicated loci are outside of known or currently hypothesized pathways to thrombosis., Conclusions: Our cross-ancestry GWAS meta-analyses identified new loci associated with VTE. These findings highlight new pathways to thrombosis and provide novel molecules that may be useful in the development of improved antithrombosis treatments.

Inter-individual variation in the number of copies of the mitochondrial genome, called mitochondrial DNA copy number (mtDNA-CN), reflects mitochondrial function and has been associated with various aging-related diseases. We examined 415,422 exomes of self-reported White ancestry individuals from the UK Biobank and tested the impact of rare variants, at the level of single variants and through aggregate variant-set tests, on mtDNA-CN. A survey across nine variant sets tested enrichment of putatively causal variants and identified 14 genes at experiment-wide significance and three genes at marginal significance. These included associations at known mtDNA depletion syndrome genes (mtDNA helicase TWNK , p = 1.1 × 10 -30 ; mitochondrial transcription factor TFAM , p = 4.3 × 10 -15 ; mtDNA maintenance exonuclease MGME1 , p = 2.0 × 10 -6 ) and the V617F dominant gain-of-function mutation in the tyrosine kinase JAK2 (p = 2.7 × 10 -17 ), associated with myeloproliferative disease. Novel genes included the ATP-dependent protease CLPX (p = 8.4 × 10 -9 ), involved in mitochondrial proteome quality, and the mitochondrial adenylate kinase AK2 (p = 4.7 × 10 -8 ), involved in hematopoiesis. The most significant association was a missense variant in SAMHD1 (p = 4.2 × 10 -28 ), found on a rare, 1.2-Mb shared ancestral haplotype on chromosome 20. SAMHD1 encodes a cytoplasmic host restriction factor involved in viral defense response and the mitochondrial nucleotide salvage pathway, and is associated with Aicardi-Goutières syndrome 5, a childhood encephalopathy and chronic inflammatory response disorder. Rare variants were enriched in Mendelian mtDNA depletion syndrome loci, and these variants implicated core processes in mtDNA replication, nucleoid structure formation, and maintenance. These data indicate that strong-effect mutations from the nuclear genome contribute to the genetic architecture of mtDNA-CN., Competing Interests: The authors declare no competing interests., (© 2022 The Author(s).)

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors., (© The Author(s) 2022. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Context: Ectopic posterior pituitary (EPP), a condition in which the posterior pituitary gland is displaced due to defective neuronal migration, is frequently associated with hypopituitarism. Genetic variants play a role, but many cases remain unexplained., Objective: A large EPP cohort was studied to explore the importance of genetic variants and how they correlate with clinical findings., Methods: Whole exome sequencing was performed on a discovery sample of 27 cases to identify rare variants. The variants that met the criteria for rarity and biological relevance, or that were previously associated with EPP ( ROBO1 and HESX1 ), were then resequenced in the 27 cases plus a replication sample of 51 cases., Results: We identified 16 different variants in 12 genes in 15 of the 78 cases (19.2%). Complete anterior pituitary deficiency was twice as common in cases with variants of interest compared to cases without variants (9/15 [60%] vs 19/63 [30.1%], respectively; Z test, P  = 0.06). Breech presentation was more frequent in the variant positive group (5/15 vs 1/63; Z test, P  = 0.003). Four cases had variants in ROBO1 and 1 in HESX1 , genes previously associated with EPP. The ROBO1 p.S18* variant has not been reported previously; ROBO1 p.Q1227H has not been associated with EPP previously., Conclusion: EPP cases with variants of interest identified in this study were more likely to present with severe clinical disease. Several variants were identified in genes not previously associated with EPP. Our findings confirm that EPP is a multigenic disorder. Future studies are needed to identify additional genes., (Published by Oxford University Press on behalf of the Endocrine Society 2022.)

Background: Genetically predicted leukocyte telomere length (LTL) has been evaluated in several studies of childhood and adult cancer. We test whether genetically predicted longer LTL is associated with germ cell tumours (GCT) in children and adults., Methods: Paediatric GCT samples were obtained from a Children's Oncology Group study and state biobank programs in California and Michigan (N = 1413 cases, 1220 biological parents and 1022 unrelated controls). Replication analysis included 396 adult testicular GCTs (TGCT) and 1589 matched controls from the UK Biobank. Mendelian randomisation was used to look at the association between genetically predicted LTL and GCTs and TERT variants were evaluated within GCT subgroups., Results: We identified significant associations between TERT variants reported in previous adult TGCT GWAS in paediatric GCT: TERT/rs2736100-C (OR = 0.82; P = 0.0003), TERT/rs2853677-G (OR = 0.80; P = 0.001), and TERT/rs7705526-A (OR = 0.81; P = 0.003). We also extended these findings to females and tumours outside the testes. In contrast, we did not observe strong evidence for an association between genetically predicted LTL by other variants and GCT risk in children or adults., Conclusion: While TERT is a known susceptibility locus for GCT, our results suggest that LTL predicted by other variants is not strongly associated with risk in either children or adults., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Mitochondrial diseases are a heterogeneous group of disorders that can be caused by mutations in the nuclear or mitochondrial genome. Mitochondrial DNA (mtDNA) variants may exist in a state of heteroplasmy, where a percentage of DNA molecules harbor a variant, or homoplasmy, where all DNA molecules have the same variant. The relative quantity of mtDNA in a cell, or copy number (mtDNA-CN), is associated with mitochondrial function, human disease, and mortality. To facilitate accurate identification of heteroplasmy and quantify mtDNA-CN, we built a bioinformatics pipeline that takes whole genome sequencing data and outputs mitochondrial variants, and mtDNA-CN. We incorporate variant annotations to facilitate determination of variant significance. Our pipeline yields uniform coverage by remapping to a circularized chrM and by recovering reads falsely mapped to nuclear-encoded mitochondrial sequences. Notably, we construct a consensus chrM sequence for each sample and recall heteroplasmy against the sample's unique mitochondrial genome. We observe an approximately 3-fold increased association with age for heteroplasmic variants in non-homopolymer regions and, are better able to capture genetic variation in the D-loop of chrM compared to existing software. Our bioinformatics pipeline more accurately captures features of mitochondrial genetics than existing pipelines that are important in understanding how mitochondrial dysfunction contributes to disease., (© The Author(s) 2022. Published by Oxford University Press on behalf of NAR Genomics and Bioinformatics.)

Human genetic studies support an inverse causal relationship between leukocyte telomere length (LTL) and coronary artery disease (CAD), but directionally mixed effects for LTL and diverse malignancies. Clonal hematopoiesis of indeterminate potential (CHIP), characterized by expansion of hematopoietic cells bearing leukemogenic mutations, predisposes both hematologic malignancy and CAD. TERT (which encodes telomerase reverse transcriptase) is the most significantly associated germline locus for CHIP in genome-wide association studies. Here, we investigated the relationship between CHIP, LTL, and CAD in the Trans-Omics for Precision Medicine (TOPMed) program ( n = 63,302) and UK Biobank ( n = 47,080). Bidirectional Mendelian randomization studies were consistent with longer genetically imputed LTL increasing propensity to develop CHIP, but CHIP then, in turn, hastens to shorten measured LTL (mLTL). We also demonstrated evidence of modest mediation between CHIP and CAD by mLTL. Our data promote an understanding of potential causal relationships across CHIP and LTL toward prevention of CAD.

Background: Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA., Methods: Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced., Results: Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children., Conclusions: To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies., (© 2022 Wiley Periodicals LLC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.)

Many individual genetic risk loci have been associated with multiple common human diseases. However, the molecular basis of this pleiotropy often remains unclear. We present an integrative approach to reveal the molecular mechanism underlying the PROCR locus, associated with lower coronary artery disease (CAD) risk but higher venous thromboembolism (VTE) risk. We identify PROCR-p.Ser219Gly as the likely causal variant at the locus and protein C as a causal factor. Using genetic analyses, human recall-by-genotype and in vitro experimentation, we demonstrate that PROCR-219Gly increases plasma levels of (activated) protein C through endothelial protein C receptor (EPCR) ectodomain shedding in endothelial cells, attenuating leukocyte-endothelial cell adhesion and vascular inflammation. We also associate PROCR-219Gly with an increased pro-thrombotic state via coagulation factor VII, a ligand of EPCR. Our study, which links PROCR-219Gly to CAD through anti-inflammatory mechanisms and to VTE through pro-thrombotic mechanisms, provides a framework to reveal the mechanisms underlying similar cross-phenotype associations., (© 2022. The Author(s).)