Search

Psychosocial trauma has been hypothesized to influence breast cancer risk, but little is known about how co-occurring traumas-particularly during early life-may impact incidence. We examine the relationship between multiple measures of early-life trauma and incident breast cancer.The Sister Study is a prospective cohort study of US women (n = 50,884; enrollment 2003-2009; ages 35-74). Of 45,961 eligible participants, 3,070 developed invasive breast cancer or ductal carcinoma in situ through 2017. We assessed trauma before age 18 using previously studied measures (cumulative score, individual trauma type, and substantive domain) and a six-class latent variable to evaluate co-occurring traumas. We accounted for missing data using multiple imputation and estimated hazard ratios (HRs) and 95% confidence intervals (CIs) using Cox proportional-hazards models.Approximately 49% of participants reported early-life trauma. Using the latent class variable approach, breast cancer hazard was higher among participants who had sexual trauma or household dysfunction (HR = 1.1; CI = 0.93, 1.3) or moderate (HR = 1.2; CI = 0.99, 1.4) but not high trauma (HR = 0.66; CI = 0.44, 0.99) compared to low trauma. Breast cancer HRs associated with sexual early-life trauma or household dysfunction were elevated for pre- and postmenopausal breast cancer and by estrogen receptor status. We found no effect modification by race-ethnicity. Estimated effects were attenuated with report of constant childhood social support.Breast cancer incidence varied by latent patterns of co-occurring early-life trauma. Models capturing childhood social support and trauma patterning, rather than cumulative or discrete indicators, may be more meaningful in breast cancer risk assessment.

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

Background and Purpose: Stroke is the leading cause of death and long-term disability worldwide. Previous genome-wide association studies identified 51 loci associated with stroke (mostly ischemic) and its subtypes among predominantly European populations. Using whole-genome sequencing in ancestrally diverse populations from the Trans-Omics for Precision Medicine (TOPMed) Program, we aimed to identify novel variants, especially low-frequency or ancestry-specific variants, associated with all stroke, ischemic stroke and its subtypes (large artery, cardioembolic, and small vessel), and hemorrhagic stroke and its subtypes (intracerebral and subarachnoid). Methods: Whole-genome sequencing data were available for 6833 stroke cases and 27 116 controls, including 22 315 European, 7877 Black, 2616 Hispanic/Latino, 850 Asian, 54 Native American, and 237 other ancestry participants. In TOPMed, we performed single variant association analysis examining 40 million common variants and aggregated association analysis focusing on rare variants. We also combined TOPMed European populations with over 28 000 additional European participants from the UK BioBank genome-wide array data through meta-analysis. Results: In the single variant association analysis in TOPMed, we identified one novel locus 13q33 for large artery at whole-genome-wide significance ( P −9 ) and 4 novel loci at genome-wide significance ( P −8 ), all of which need confirmation in independent studies. Lead variants in all 5 loci are low-frequency but are more common in non-European populations. An aggregation of synonymous rare variants within the gene C6orf26 demonstrated suggestive evidence of association for hemorrhagic stroke ( P −6 ). By meta-analyzing European ancestry samples in TOPMed and UK BioBank, we replicated several previously reported stroke loci including PITX2 , HDAC9 , ZFHX3 , and LRCH1 . Conclusions: We represent the first association analysis for stroke and its subtypes using whole-genome sequencing data from ancestrally diverse populations. While our findings suggest the potential benefits of combining whole-genome sequencing data with populations of diverse genetic backgrounds to identify possible low-frequency or ancestry-specific variants, they also highlight the need to increase genome coverage and sample sizes.

In this study, the asymptotic distributions of the likelihood ratio test (LRT), the restricted likelihood ratio test (RLRT), the F and the sequence kernel association test (SKAT) statistics for testing an additive effect of the expected familial relatedness (FR) in a linear mixed model are examined based on an eigenvalue approach. First, the covariance structure for modeling the FR effect in a LMM is presented. Then, the multiplicity of eigenvalues for the log-likelihood and restricted log-likelihood is established under a replicate family setting and extended to a more general replicate family setting (GRFS) as well. After that, the asymptotic null distributions of LRT, RLRT, F and SKAT statistics under GRFS are derived. The asymptotic null distribution of SKAT for testing genetic rare variants is also constructed. In addition, a simple formula for sample size calculation is provided based on the restricted maximum likelihood estimate of the effect size for the expected FR. Finally, a power comparison of these test statistics on hypothesis test of the expected FR effect is made via simulation. The four test statistics are also applied to a data set from the UK Biobank.

Many common and rare variants associated with hematologic traits have been discovered through imputation on large-scale reference panels. However, the majority of genome-wide association studies (GWASs) have been conducted in Europeans, and determining causal variants has proved challenging. We performed a GWAS of total leukocyte, neutrophil, lymphocyte, monocyte, eosinophil, and basophil counts generated from 109,563,748 variants in the autosomes and the X chromosome in the Trans-Omics for Precision Medicine (TOPMed) program, which included data from 61,802 individuals of diverse ancestry. We discovered and replicated 7 leukocyte trait associations, including (1) the association between a chromosome X, pseudo-autosomal region (PAR), noncoding variant located between cytokine receptor genes (CSF2RA and CLRF2) and lower eosinophil count; and (2) associations between single variants found predominantly among African Americans at the S1PR3 (9q22.1) and HBB (11p15.4) loci and monocyte and lymphocyte counts, respectively. We further provide evidence indicating that the newly discovered eosinophil-lowering chromosome X PAR variant might be associated with reduced susceptibility to common allergic diseases such as atopic dermatitis and asthma. Additionally, we found a burden of very rare FLT3 (13q12.2) variants associated with monocyte counts. Together, these results emphasize the utility of whole-genome sequencing in diverse samples in identifying associations missed by European-ancestry-driven GWASs.

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

Mosaic mutations in blood are common with increasing age and are prognostic markers for cancer, cardiovascular dysfunction and other diseases. This group of acquired mutations include megabase-scale mosaic chromosomal alterations (mCAs). These large mutations have mainly been surveyed using SNP array data from individuals of European (EA) or Japanese genetic ancestry. To gain a better understanding of mCA rates and associated risk factors in genetically diverse populations, we surveyed whole genome sequencing data from 67,390 individuals, including 20,132 individuals of African ancestry (AA), and 7,608 of Hispanic ancestry (HA) with deep (30X) whole genome sequencing data from the NHLBI Trans Omics for Precision Medicine (TOPMed) program. We adapted an existing mCA calling algorithm for application to WGS data, and observed higher sensitivity with WGS data, compared with array-based data, in uncovering mCAs at low mutant cell fractions. As in previous reports, we observed a strong association with age and a non-uniform distribution of mCAs across the genome. The presence of autosomal (but not chromosome X) mCAs was associated with an increased risk of both lymphoid and myeloid malignancies. After adjusting for age, we found that individuals of European ancestry have the highest rates of autosomal mCAs, mirroring the higher rate of leukemia in this group. Our analysis also uncovered higher rates of chromosome X mCAs in AA and HA compared to EA, again after adjusting for age. Germline variants inATMandMPLshowed strong associations with mCAs incis, including ancestry specific variants. And rare variant gene-burden analysis confirmed the association of putatively protein altering variants inATMandMPLwith mCAs incis. Individual rare variants inDCPS, ADM17, PPP1R16B, andTET2were all associated with autosomal mCAs and rare variants inOR4C16were associated with chromosome X mCAs in females. There was significant enrichment of co-occurrence of CHIP mutations and mCAs both altering cancer associated genesTET2, DNMT3A, JAK2, CUX1, andTP53. Overall, our study demonstrates that rates of mCAs differ across populations and that rare inherited germline variants are strongly associated with mCAs across genetically diverse populations. These results strongly motivate further studies of mCAs in under-represented populations to better understand the causes and consequences of this class of somatic variation.

Polygenic risk scores (PRS) have shown successes in clinics, but most PRS methods have focused only on individuals with one primary continental ancestry, thus poorly accommodating recently-admixed individuals. Here, we develop GAUDI, a novel penalized-regression-based method specifically designed for admixed individuals by explicitly modeling ancestry-specific effects and jointly estimating ancestry-shared effects. We demonstrate marked advantages of GAUDI over other methods through comprehensive simulation and real data analyses.

PURPOSE This study assessed the joint association of pathogenic variants (PVs) in breast cancer (BC) predisposition genes and polygenic risk scores (PRS) with BC in the general population. METHODS A total of 26,798 non-Hispanic white BC cases and 26,127 controls from predominately population-based studies in the Cancer Risk Estimates Related to Susceptibility consortium were evaluated for PVs in BRCA1, BRCA2, ATM, CHEK2, PALB2, BARD1, BRIP1, CDH1, and NF1. PRS based on 105 common variants were created using effect estimates from BC genome-wide association studies; the performance of an overall BC PRS and estrogen receptor–specific PRS were evaluated. The odds of BC based on the PVs and PRS were estimated using penalized logistic regression. The results were combined with age-specific incidence rates to estimate 5-year and lifetime absolute risks of BC across percentiles of PRS by PV status and first-degree family history of BC. RESULTS The estimated lifetime risks of BC among general-population noncarriers, based on 10th and 90th percentiles of PRS, were 9.1%-23.9% and 6.7%-18.2% for women with or without first-degree relatives with BC, respectively. Taking PRS into account, more than 95% of BRCA1, BRCA2, and PALB2 carriers had > 20% lifetime risks of BC, whereas, respectively, 52.5% and 69.7% of ATM and CHEK2 carriers without first-degree relatives with BC, and 78.8% and 89.9% of those with a first-degree relative with BC had > 20% risk. CONCLUSION PRS facilitates personalization of BC risk among carriers of PVs in predisposition genes. Incorporating PRS into BC risk estimation may help identify > 30% of CHEK2 and nearly half of ATM carriers below the 20% lifetime risk threshold, suggesting the addition of PRS may prevent overscreening and enable more personalized risk management approaches.

Background Circulating white blood cell and platelet traits are clinically linked to various disease outcomes and differ across individuals and ancestry groups. Genetic factors play an important role in determining these traits and many loci have been identified. However, most of these findings were identified in populations of European ancestry (EA), with African Americans (AA), Hispanics/Latinos (HL), and other races/ethnicities being severely underrepresented. Results We performed ancestry-combined and ancestry-specific genome-wide association studies (GWAS) for white blood cell and platelet traits in the ancestrally diverse Population Architecture using Genomics and Epidemiology (PAGE) Study, including 16,201 AA, 21,347 HL, and 27,236 EA participants. We identified six novel findings at suggestive significance (P < 5E-8), which need confirmation, and independent signals at six previously established regions at genome-wide significance (P < 2E-9). We confirmed multiple previously reported genome-wide significant variants in the single variant association analysis and multiple genes using PrediXcan. Evaluation of loci reported from a Euro-centric GWAS indicated attenuation of effect estimates in AA and HL compared to EA populations. Conclusions Our results highlighted the potential to identify ancestry-specific and ancestry-agnostic variants in participants with diverse backgrounds and advocate for continued efforts in improving inclusion of racially/ethnically diverse populations in genetic association studies for complex traits.

Background: Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. Methods: We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values-8 as genome-wide significant, and p-values-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants.Results: None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values5. The strongest evidence was found for rs4674019 (p-value=2.27x10-7), which showed genome-wide significant interaction (p-value=3.8x10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. Conclusions: In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association.

Summary Whole-genome sequencing (WGS), a powerful tool for detecting novel coding and non-coding disease-causing variants, has largely been applied to clinical diagnosis of inherited disorders. Here we leveraged WGS data in up to 62,653 ethnically diverse participants from the NHLBI Trans-Omics for Precision Medicine (TOPMed) program and assessed statistical association of variants with seven red blood cell (RBC) quantitative traits. We discovered 14 single variant-RBC trait associations at 12 genomic loci, which have not been reported previously. Several of the RBC trait-variant associations (RPN1, ELL2, MIDN, HBB, HBA1, PIEZO1, and G6PD) were replicated in independent GWAS datasets imputed to the TOPMed reference panel. Most of these discovered variants are rare/low frequency, and several are observed disproportionately among non-European Ancestry (African, Hispanic/Latino, or East Asian) populations. We identified a 3 bp indel p.Lys2169del (g.88717175_88717177TCT[4]) (common only in the Ashkenazi Jewish population) of PIEZO1, a gene responsible for the Mendelian red cell disorder hereditary xerocytosis (MIM: 194380 ), associated with higher mean corpuscular hemoglobin concentration (MCHC). In stepwise conditional analysis and in gene-based rare variant aggregated association analysis, we identified several of the variants in HBB, HBA1, TMPRSS6, and G6PD that represent the carrier state for known coding, promoter, or splice site loss-of-function variants that cause inherited RBC disorders. Finally, we applied base and nuclease editing to demonstrate that the sentinel variant rs112097551 (nearest gene RPN1) acts through a cis-regulatory element that exerts long-range control of the gene RUVBL1 which is essential for hematopoiesis. Together, these results demonstrate the utility of WGS in ethnically diverse population-based samples and gene editing for expanding knowledge of the genetic architecture of quantitative hematologic traits and suggest a continuum between complex trait and Mendelian red cell disorders.

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%., The goals, resources and design of the NHLBI Trans-Omics for Precision Medicine (TOPMed) programme are described, and analyses of rare variants detected in the first 53,831 samples provide insights into mutational processes and recent human evolutionary history.

BACKGROUND: Population-based estimates of the risk of breast cancer associated with germline pathogenic variants in cancer-predisposition genes are critically needed for risk assessment and management in women with inherited pathogenic variants. METHODS: In a population-based case–control study, we performed sequencing using a custom multigene amplicon-based panel to identify germline pathogenic variants in 28 cancer-predisposition genes among 32,247 women with breast cancer (case patients) and 32,544 unaffected women (controls) from population-based studies in the Cancer Risk Estimates Related to Susceptibility (CARRIERS) consortium. Associations between pathogenic variants in each gene and the risk of breast cancer were assessed. RESULTS: Pathogenic variants in 12 established breast cancer–predisposition genes were detected in 5.03% of case patients and in 1.63% of controls. Pathogenic variants in BRCA1 and BRCA2 were associated with a high risk of breast cancer, with odds ratios of 7.62 (95% confidence interval [CI], 5.33 to 11.27) and 5.23 (95% CI, 4.09 to 6.77), respectively. Pathogenic variants in PALB2 were associated with a moderate risk (odds ratio, 3.83; 95% CI, 2.68 to 5.63). Pathogenic variants in BARD1, RAD51C, and RAD51D were associated with increased risks of estrogen receptor–negative breast cancer and triple-negative breast cancer, whereas pathogenic variants in ATM, CDH1, and CHEK2 were associated with an increased risk of estrogen receptor–positive breast cancer. Pathogenic variants in 16 candidate breast cancer–predisposition genes, including the c.657_661del5 founder pathogenic variant in NBN, were not associated with an increased risk of breast cancer. CONCLUSIONS: This study provides estimates of the prevalence and risk of breast cancer associated with pathogenic variants in known breast cancer–predisposition genes in the U.S. population. These estimates can inform cancer testing and screening and improve clinical management strategies for women in the general population with inherited pathogenic variants in these genes. (Funded by the National Institutes of Health and the Breast Cancer Research Foundation.)

Telomere length, a biomarker of cell division and cellular aging, has been associated with multiple chronic disease endpoints. Experienced trauma over the life course may contribute to telomere shortening via mechanisms of stress embodiment. However, it is unclear how patterns of co-occurring trauma during sensitive periods (e.g., early life) throughout the life course may influence telomere shortening. We examine the relationship between co-occurring early life trauma on adult telomere length and the extent to which adulthood trauma, socioeconomic position, and health and lifestyle factors may mediate this relationship.We use data from a sample of participants in the Sister Study (N = 740, analytic sample: n = 602), a prospective cohort of U.S. self-identified females aged 35-74 years at enrollment (2003-2009) for whom leukocyte telomere length was measured in baseline blood samples. Participants reported their experience of 20 different types of trauma, from which we identified patterns of co-occurring early life trauma (before age 18) using latent class analysis. We estimated the direct and indirect effects of early life trauma on leukocyte telomere length using structural equation modeling, allowing for mediating adult pathways.Approximately 47 % of participants reported early life trauma. High early life trauma was associated with shorter telomere length compared to low early life trauma (β = -0.11; 95 % CI: -0.22, -0.004) after adjusting for age and childhood socioeconomic position. The inverse association between early life trauma and adult leukocyte telomere length was largely attributable to the direct effect of early life trauma on telomere length (β = -0.12; 95 %CI: -0.23, -0.01). Mediating indirect pathways via adult trauma, socioeconomic position, and health metrics did not substantively contribute the overall association.These findings highlight the role of patterns of co-occurring early life trauma on shortened telomere length independent of adult pathways.

International audience; Mendelian randomization (MR) is an established approach for assessing the causal effects of heritable exposures on outcomes. Outcomes of interest often include binary clinical endpoints, but may also include censored survival times. We explore the implications of both the Cox proportional hazard model and the additive hazard model in the context of MR, with a specific emphasis on two-stage methods. We show that naive application of standard MR approaches to censored survival times may induce significant bias. Through simulations and analysis of data from the Women's Health Initiative, we provide practical advice on modeling survival outcomes in MRs.

Background Prospective studies have suggested higher factor VIII (FVIII) levels are an independent risk factor for coronary heart disease (CHD) and stroke. However, limited information, including on genetic and epigenetic contributors to FVIII variation, is available specifically among African Americans (AAs), who have higher FVIII levels than Europeans. Objectives We measured FVIII levels in ~3400 AAs from the community-based Jackson Heart Study and assessed genetic, epigenetic, and epidemiological correlates of FVIII, as well as incident cardiovascular disease (CVD) associations. Methods We assessed cross-sectional associations of FVIII with CVD risk factors as well as incident CHD, stroke, heart failure, and mortality associations. We additionally assessed associations with TOPMed whole genome sequencing data and an epigenome-wide methylation array. Results Our results confirmed associations between FVIII and risk of incident CHD events and total mortality in AAs; mortality associations were largely independent of traditional risk factors. We also demonstrate an association of FVIII with incident heart failure, independent of B-type natriuretic peptide. Two genomic regions were strongly associated with FVIII (ABO and VWF). The index variant at VWF is specific to individuals of African descent and is distinct from the previously reported European VWF association signal. Epigenome-wide association analysis showed significant FVIII associations with several CpG sites in the ABO region. However, after adjusting for ABO genetic variants, ABO CpG sites were not significant. Conclusions Larger sample sizes of AAs will be required to discover additional genetic and epigenetic contributors to FVIII phenotypic variation, which may have consequences for CVD health disparities.

The objective of the current study is to use comparative and functional genomic analysis to help to understand the biological mechanism mediating the effect of single nucleotide polymorphisms (SNPs) on blood pressure. We mapped 26 585 SNPs that are in linkage disequilibrium with 1071 human blood pressure–associated sentinel SNPs to 9447 syntenic regions in the mouse genome. Approximately 21.8% of the 1071 linkage disequilibrium regions are located at least 10 kb from any protein-coding gene. Approximately 300 blood pressure–associated SNPs are expression quantitative trait loci for a few dozen known blood pressure physiology genes in tissues including specific kidney regions. Blood pressure–associated sentinel SNPs are significantly enriched for expression quantitative trait loci for blood pressure physiology genes compared with randomly selected SNPs ( P P =4.90×10 −7 , Pearson χ 2 test). One example synonymous SNP rs9337951 was predicted to influence the secondary structure of its host mRNA JCAD (junctional cadherin 5 associated) and was experimentally validated to influence JCAD protein expression. These findings provide an extensive comparative and functional genomic resource for developing experiments to test the functional significance of human blood pressure–associated SNPs in human cells and animal models.

Genome-wide analysis identifies 32 loci associated with breast cancer susceptibility, accounting for estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype(1-3). To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 x 10(-8)), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores.

Most loci identified by GWASs have been found in populations of European ancestry (EUR). In trans-ethnic meta-analyses for 15 hematological traits in 746,667 participants, including 184,535 non-EUR individuals, we identified 5,552 trait-variant associations at p < 5 × 10-9, including 71 novel associations not found in EUR populations. We also identified 28 additional novel variants in ancestry-specific, non-EUR meta-analyses, including an IL7 missense variant in South Asians associated with lymphocyte count in vivo and IL-7 secretion levels in vitro. Fine-mapping prioritized variants annotated as functional and generated 95% credible sets that were 30% smaller when using the trans-ethnic as opposed to the EUR-only results. We explored the clinical significance and predictive value of trans-ethnic variants in multiple populations and compared genetic architecture and the effect of natural selection on these blood phenotypes between populations. Altogether, our results for hematological traits highlight the value of a more global representation of populations in genetic studies.

Current publicly available tools that allow rapid exploration of linkage disequilibrium (LD) between markers (e.g., HaploReg and LDlink) are based on whole-genome sequence (WGS) data from 2,504 individuals in the 1000 Genomes Project. Here, we present TOP-LD, an online tool to explore LD inferred with high-coverage (∼30×) WGS data from 15,578 individuals in the NHLBI Trans-Omics for Precision Medicine (TOPMed) program. TOP-LD provides a significant upgrade compared to current LD tools, as the TOPMed WGS data provide a more comprehensive representation of genetic variation than the 1000 Genomes data, particularly for rare variants and in the specific populations that we analyzed. For example, TOP-LD encompasses LD information for 150.3, 62.2, and 36.7 million variants for European, African, and East Asian ancestral samples, respectively, offering 2.6- to 9.1-fold increase in variant coverage compared to HaploReg 4.0 or LDlink. In addition, TOP-LD includes tens of thousands of structural variants (SVs). We demonstrate the value of TOP-LD in fine-mapping at the GGT1 locus associated with gamma glutamyltransferase in the African ancestry participants in UK Biobank. Beyond fine-mapping, TOP-LD can facilitate a wide range of applications that are based on summary statistics and estimates of LD. TOP-LD is freely available online.

Background: suPAR (Soluble urokinase plasminogen activator receptor) has emerged as an important biomarker of coagulation, inflammation, and cardiovascular disease (CVD) risk. The contribution of suPAR to CVD risk and its genetic influence in Black populations have not been evaluated. Methods: We measured suPAR in 3492 Black adults from the prospective, community-based JHS (Jackson Heart Study). Cross-sectional associations of suPAR with lifestyle and CVD risk factors were assessed, whole-genome sequence data were used to evaluate genetic associations of suPAR, and relationships of suPAR with incident CVD outcomes and overall mortality were estimated over follow-up. Results: In Cox models adjusted for traditional CVD risk factors, estimated glomerular filtration rate, and CRP (C-reactive protein), each 1-SD higher suPAR was associated with a 21% to 31% increased risk of incident coronary heart disease, heart failure, stroke, and mortality. In the genome-wide association study, 2 missense (rs399145 encoding p.Thr86Ala, rs4760 encoding p.Phe272Leu) and 2 noncoding regulatory variants (rs73935023 within an enhancer element and rs4251805 within the promoter) of PLAUR on chromosome 19 were each independently associated with suPAR and together explained 14% of suPAR phenotypic variation. The allele frequencies of each of the four suPAR-associated genetic variants differ considerably across African and European populations. We further show that PLAUR rs73935023 can alter transcriptional activity in vitro. We did not find any association between genetically determined suPAR and CVD in JHS or a larger electronic medical record-based analyses of Blacks or Whites. Conclusions: Our results demonstrate the importance of ancestry-differentiated genetic variation on suPAR levels and indicate suPAR is a CVD biomarker in Black adults.

Soluble urokinase plasminogen activator receptor (suPAR) is the circulating form of urokinase plasminogen activator receptor, a glycosyl-phosphatidylinositol–anchored membrane protein expressed in various cell types including kidney podocytes and endothelial cells ([1][1]). suPAR has been

PURPOSE: Understanding the molecular mediators of breast cancer survival is critical for accurate disease prognosis and improving therapies. Here, we identified Neuronatin (NNAT) as a novel antiproliferative modifier of estrogen receptor-alpha (ER+) breast cancer. EXPERIMENTAL DESIGN: Genomic regions harboring breast cancer modifiers were identified by congenic mapping in a rat model of carcinogen-induced mammary cancer. Tumors from susceptible and resistant congenics were analyzed by RNAseq to identify candidate genes. Candidates were prioritized by correlation with outcome, using a consensus of three breast cancer patient cohorts. NNAT was transgenically expressed in ER+ breast cancer lines (T47D and ZR75), followed by transcriptomic and phenotypic characterization. RESULTS: We identified a region on rat chromosome 3 (142–178 Mb) that modified mammary tumor incidence. RNAseq of the mammary tumors narrowed the candidate list to three differentially expressed genes: NNAT, SLC35C2, and FAM210B. NNAT mRNA and protein also correlated with survival in human breast cancer patients. Quantitative immunohistochemistry of NNAT protein revealed an inverse correlation with survival in a univariate analysis of patients with invasive ER+ breast cancer (training cohort: n = 444, HR = 0.62, p = 0.031; validation cohort: n = 430, HR = 0.48, p = 0.004). NNAT also held up as an independent predictor of survival after multivariable adjustment (HR = 0.64, p = 0.038). NNAT significantly reduced proliferation and migration of ER+ breast cancer cells, which coincided with altered expression of multiple related pathways. CONCLUSIONS: Collectively, these data implicate NNAT as a novel mediator of cell proliferation and migration, which correlates with decreased tumorigenic potential and prolonged patient survival.

Breast cancer risk is 31% heritable, yet the majority of the underlying risk factors remain poorly defined. Here, we used F2-linkage analysis in a rat mammary tumor model to identify a novel 11.2 Mb modifier locus of tumor incidence and burden on rat chromosome 5 (chr5: 15.4 – 26.6 Mb). Genomic and RNA sequencing analysis identified four differentially expressed candidates: TMEM68, IMPAD1, SDCBP, and RBM12B. Analysis of the human syntenic candidate region revealed that SDCBP is in close proximity to a previously reported genetic risk locus for human breast cancer. Moreover, analysis of the candidate genes in The Cancer Genome Atlas (TCGA) revealed that they fall within the commonly amplified 8q12.1 and 8q22.1 regions in human breast cancer patients and are correlated with worse overall survival. Collectively, this study presents novel evidence suggesting that TMEM68, IMPAD1, SDCBP, and RBM12B are potential modifiers of human breast cancer risk and outcome.

Objective: sCD14 (soluble CD14) is a circulating pattern recognition receptor involved in inflammatory signaling. Although it is known that sCD14 levels vary by race, information on the genetic and cardiovascular disease (CVD) risk relationships of sCD14 in Black participants is limited. Approach and Results: We measured sCD14 in plasma at the baseline exam from n=3492 Black participants from the JHS (Jackson Heart Study). We evaluated associations between sCD14 and subclinical CVD measures, incident CVD events, and mortality under 3 levels of covariate adjustment. We used whole-genome sequence data from the Trans-Omics for Personalized Medicine program to identify genetic associations with sCD14. Adjusting for CVD risk factors and C-reactive protein, higher sCD14 was significantly associated with increased risk of mortality (hazard ratio, 1.25 [95% CI, 1.17–1.32]), incident coronary heart disease (hazard ratio, 1.28 [95% CI, 1.11–1.47]), and incident heart failure (hazard ratio, 1.27 [95% CI, 1.15–1.41]), but not stroke (hazard ratio, 0.96 [95% CI, 0.84–1.09]). Some of these relationships differed by age or sex: the association between sCD14 and heart failure was only observed in females; there was an association between sCD14 and stroke only at younger ages (in the lowest tertile of age, CD14 region of chromosome 5q31. We additionally identified 2 novel statistically distinct genetic associations with sCD14 represented by index variants rs770147646 and rs57599368, also in the chromosome 5q31 region. Conclusions: sCD14 independently predicts CVD-related outcomes and mortality in Black participants.

Summary Genome-wide association studies (GWASs) have identified hundreds of thousands of genetic variants associated with complex diseases and traits. However, most variants are noncoding and not clearly linked to genes, making it challenging to interpret these GWAS signals. We present a systematic variant-to-function study, prioritizing the most likely functional elements of the genome for experimental follow-up, for >148,000 variants identified for hematological traits. Specifically, we developed VAMPIRE: Variant Annotation Method Pointing to Interesting Regulatory Effects, an interactive web application implemented in R Shiny. This tool efficiently integrates and displays information from multiple complementary sources, including epigenomic signatures from blood-cell-relevant tissues or cells, functional and conservation summary scores, variant impact on protein and gene expression, chromatin conformation information, as well as publicly available GWAS and phenome-wide association study (PheWAS) results. Leveraging data generated from independently performed functional validation experiments, we demonstrate that our prioritized variants, genes, or variant-gene links are significantly more likely to be experimentally validated. This study not only has important implications for systematic and efficient revelation of functional mechanisms underlying GWAS variants for hematological traits but also provides a prototype that can be adapted to many other complex traits, paving the path for efficient variant-to-function (V2F) analyses., We developed VAMPIRE: Variant Annotation Method Pointing to Interesting Regulatory Effects, an interactive web application implemented in R Shiny. VAMPIRE efficiently integrates and displays information from multiple complementary sources. Using data from various functional experiments, we demonstrated that VAMPIRE-prioritized variants are more likely to play functional roles.