Clonal hematopoiesis and risk of prostate cancer in large samples of European ancestry men.

Little is known regarding the potential relationship between clonal hematopoiesis (CH) of indeterminate potential (CHIP), which is the expansion of hematopoietic stem cells with somatic mutations, and risk of prostate cancer, the fifth leading cause of cancer death of men worldwide. We evaluated the association of age-related CHIP with overall and aggressive prostate cancer risk in two large whole-exome sequencing studies of 75 047 European ancestry men, including 7663 prostate cancer cases, 2770 of which had aggressive disease, and 3266 men carrying CHIP variants. We found that CHIP, defined by over 50 CHIP genes individually and in aggregate, was not significantly associated with overall (aggregate HR = 0.93, 95% CI = 0.76-1.13, P = 0.46) or aggressive (aggregate OR = 1.14, 95% CI = 0.92-1.41, P = 0.22) prostate cancer risk. CHIP was weakly associated with genetic risk of overall prostate cancer, measured using a polygenic risk score (OR = 1.05 per unit increase, 95% CI = 1.01-1.10, P = 0.01). CHIP was not significantly associated with carrying pathogenic/likely pathogenic/deleterious variants in DNA repair genes, which have previously been found to be associated with aggressive prostate cancer. While findings from this study suggest that CHIP is likely not a risk factor for prostate cancer, it will be important to investigate other types of CH in association with prostate cancer risk.

Associations of minor histocompatibility antigens with outcomes following allogeneic hematopoietic cell transplantation.

The role of minor histocompatibility antigens (mHAs) in mediating graft versus leukemia and graft versus host disease (GvHD) following allogeneic hematopoietic cell transplantation (alloHCT) is recognized but not well-characterized. By implementing improved methods for mHA prediction in two large patient cohorts, this study aimed to comprehensively explore the role of mHAs in alloHCT by analyzing whether (1) the number of predicted mHAs, or (2) individual mHAs are associated with clinical outcomes. The study population consisted of 2249 donor-recipient pairs treated for acute myeloid leukemia and myelodysplastic syndrome with alloHCT. A Cox proportional hazard model showed that patients with a class I mHA count greater than the population median had an increased hazard of GvHD mortality (hazard ratio [HR] = 1.39, 95% confidence interval [CI] = 1.01, 1.77, p = .046). Competing risk analyses identified the class I mHAs DLRCKYISL (GSTP), WEHGPTSLL (CRISPLD2), and STSPTTNVL (SERPINF2) were associated with increased GVHD mortality (HR = 2.84, 95% CI = 1.52, 5.31, p = .01), decreased leukemia-free survival (LFS) (HR = 1.94, 95% CI = 1.27, 2.95, p = .044), and increased disease-related mortality (DRM) (HR = 2.32, 95% CI = 1.5, 3.6, p = .008), respectively. One class II mHA YQEIAAIPSAGRERQ (TACC2) was associated with increased risk of treatment-related mortality (TRM) (HR = 3.05, 95% CI = 1.75, 5.31, p = .02). WEHGPTSLL and STSPTTNVL were both present within HLA haplotype B*40:01-C*03:04 and showed a positive dose-response relationship with increased all-cause mortality and DRM and decreased LFS, indicating these two mHAs contribute to the risk of mortality in an additive manner. Our study reports the first large-scale investigation of the associations of predicted mHA peptides with clinical outcomes following alloHCT.

Population-specific reference panels are crucial for genetic analyses: an example of the CREBRF locus in Native Hawaiians.

Statistical imputation applied to genome-wide array data is the most cost-effective approach to complete the catalog of genetic variation in a study population. However, imputed genotypes in underrepresented populations incur greater inaccuracies due to ascertainment bias and a lack of representation among reference individuals, further contributing to the obstacles to study these populations. Here we examined the consequences due to the lack of representation by genotyping in a large number of self-reported Native Hawaiians (N = 3693) a functionally important, Polynesian-specific variant in the CREBRF gene, rs373863828. We found the derived allele was significantly associated with several adiposity traits with large effects (e.g. ~ 1.28 kg/m2 per allele in body mass index as the most significant; P = 7.5 × 10-5), consistent with the original findings in Samoans. Due to the current absence of Polynesian representation in publicly accessible reference sequences, rs373863828 or its proxies could not be tested through imputation using these existing resources. Moreover, the association signals at the entire CREBRF locus could not be captured by alternative approaches, such as admixture mapping. In contrast, highly accurate imputation can be achieved even if a small number (<200) of internally constructed Polynesian reference individuals were available; this would increase sample size and improve the statistical evidence of associations. Taken together, our results suggest the alarming possibility that lack of representation in reference panels could inhibit discovery of functionally important loci such as CREBRF. Yet, they could be easily detected and prioritized with improved representation of diverse populations in sequencing studies.

Exome chip analyses identify genes affecting mortality after HLA-matched unrelated-donor blood and marrow transplantation.

Although survival outcomes have significantly improved, up to 40% of patients die within 1 year of HLA-matched unrelated-donor blood and marrow transplantation (BMT). To identify non-HLA genetic contributors to mortality after BMT, we performed the first exome-wide association study in the DISCOVeRY-BMT cohorts using the Illumina HumanExome BeadChip. This study includes 2473 patients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome and 2221 10/10 HLA-matched donors treated from 2000 to 2011. Single-variant and gene-level analyses were performed on overall survival (OS), transplantation-related mortality (TRM), and disease-related mortality (DRM). Genotype mismatches between recipients and donors in a rare nonsynonymous variant of testis-expressed gene TEX38 significantly increased risk of TRM, which was more dramatic when either the recipient or donor was female. Using the SKAT-O test to evaluate gene-level effects, variant genotypes of OR51D1 in recipients were significantly associated with OS and TRM. In donors, 4 (ALPP, EMID1, SLC44A5, LRP1), 1 (HHAT), and 2 genes (LYZL4, NT5E) were significantly associated with OS, TRM, and DRM, respectively. Inspection of NT5E crystal structures showed 4 of the associated variants affected the enzyme structure and likely decreased the catalytic efficiency of the enzyme. Further confirmation of these findings and additional functional studies may provide individualized risk prediction and prognosis, as well as alternative donor selection strategies.

Replication and validation of genetic polymorphisms associated with survival after allogeneic blood or marrow transplant.

Multiple candidate gene-association studies of non-HLA single-nucleotide polymorphisms (SNPs) and outcomes after blood or marrow transplant (BMT) have been conducted. We identified 70 publications reporting 45 SNPs in 36 genes significantly associated with disease-related mortality, progression-free survival, transplant-related mortality, and/or overall survival after BMT. Replication and validation of these SNP associations were performed using DISCOVeRY-BMT (Determining the Influence of Susceptibility COnveying Variants Related to one-Year mortality after BMT), a well-powered genome-wide association study consisting of 2 cohorts, totaling 2888 BMT recipients with acute myeloid leukemia, acute lymphoblastic leukemia, or myelodysplastic syndrome, and their HLA-matched unrelated donors, reported to the Center for International Blood and Marrow Transplant Research. Gene-based tests were used to assess the aggregate effect of SNPs on outcome. None of the previously reported significant SNPs replicated at P < .05 in DISCOVeRY-BMT. Validation analyses showed association with one previously reported donor SNP at P < .05 and survival; more associations would be anticipated by chance alone. No gene-based tests were significant at P < .05. Functional annotation with publicly available data shows these candidate SNPs most likely do not have biochemical function; only 13% of candidate SNPs correlate with gene expression or are predicted to impact transcription factor binding. Of these, half do not impact the candidate gene of interest; the other half correlate with expression of multiple genes. These findings emphasize the peril of pursing candidate approaches and the importance of adequately powered tests of unbiased genome-wide associations with BMT clinical outcomes given the ultimate goal of improving patient outcomes.

GWAS meta-analysis of 16 852 women identifies new susceptibility locus for endometrial cancer.

Endometrial cancer is the most common gynecological malignancy in the developed world. Although there is evidence of genetic predisposition to the disease, most of the genetic risk remains unexplained. We present the meta-analysis results of four genome-wide association studies (4907 cases and 11 945 controls total) in women of European ancestry. We describe one new locus reaching genome-wide significance (P < 5 × 10 -8) at 6p22.3 (rs1740828; P = 2.29 × 10 -8, OR = 1.20), providing evidence of an additional region of interest for genetic susceptibility to endometrial cancer.

Characterization of large structural genetic mosaicism in human autosomes.

Analyses of genome-wide association study (GWAS) data have revealed that detectable genetic mosaicism involving large (>2 Mb) structural autosomal alterations occurs in a fraction of individuals. We present results for a set of 24,849 genotyped individuals (total GWAS set II [TGSII]) in whom 341 large autosomal abnormalities were observed in 168 (0.68%) individuals. Merging data from the new TGSII set with data from two prior reports (the Gene-Environment Association Studies and the total GWAS set I) generated a large dataset of 127,179 individuals; we then conducted a meta-analysis to investigate the patterns of detectable autosomal mosaicism (n = 1,315 events in 925 [0.73%] individuals). Restricting to events >2 Mb in size, we observed an increase in event frequency as event size decreased. The combined results underscore that the rate of detectable mosaicism increases with age (p value = 5.5 × 10(-31)) and is higher in men (p value = 0.002) but lower in participants of African ancestry (p value = 0.003). In a subset of 47 individuals from whom serial samples were collected up to 6 years apart, complex changes were noted over time and showed an overall increase in the proportion of mosaic cells as age increased. Our large combined sample allowed for a unique ability to characterize detectable genetic mosaicism involving large structural events and strengthens the emerging evidence of non-random erosion of the genome in the aging population.

Replication of associations between genetic polymorphisms and chronic graft-versus-host disease.

Previous studies have identified single-nucleotide polymorphisms (SNPs) associated with the risk of chronic graft-versus-host disease (GVHD) after allogeneic hematopoietic cell transplantation. The current study determined whether these associations could be replicated in large cohorts of donors and recipients. Each SNP was tested with cohorts of patients having the same donor type (HLA-matched related, unrelated, or both) reported in the original publication, and testing was limited to the same genome (recipient or donor) and genetic model (dominant, recessive, or allelic) reported in the original study. The 21 SNPs reported in this study represent 19 genes, and the analysis encompassed 22 SNP association tests. The hazard ratio (HR) point estimates and risk ratio point estimates corresponding to odds ratios in previous studies consistently fall outside the 95% confidence intervals of HR estimates in the current study. Despite the large size of the cohorts available for the current study, the 95% confidence intervals for most HRs did not exclude 1.0. Three SNPs representing CTLA4, HPSE, and IL1R1 showed evidence of association with the risk of chronic GVHD in unrelated donor-recipient pairs from 1 cohort, but none of these associations was replicated when tested in unrelated donor-recipient pairs from an independent cohort. Two SNPs representing CCR6 and FGFR1OP showed possible associations with the risk of chronic GVHD in related donor-recipient pairs but not in unrelated donor-recipient pairs. These results remain to be tested for replication in other cohorts of related donor-recipient pairs.

Fine Mapping and Identification of BMI Loci in African Americans.

Genome-wide association studies (GWASs) primarily performed in European-ancestry (EA) populations have identified numerous loci associated with body mass index (BMI). However, it is still unclear whether these GWAS loci can be generalized to other ethnic groups, such as African Americans (AAs). Furthermore, the putative functional variant or variants in these loci mostly remain under investigation. The overall lower linkage disequilibrium in AA compared to EA populations provides the opportunity to narrow in or fine-map these BMI-related loci. Therefore, we used the Metabochip to densely genotype and evaluate 21 BMI GWAS loci identified in EA studies in 29,151 AAs from the Population Architecture using Genomics and Epidemiology (PAGE) study. Eight of the 21 loci (SEC16B, TMEM18, ETV5, GNPDA2, TFAP2B, BDNF, FTO, and MC4R) were found to be associated with BMI in AAs at 5.8 × 10(-5). Within seven out of these eight loci, we found that, on average, a substantially smaller number of variants was correlated (r(2) > 0.5) with the most significant SNP in AA than in EA populations (16 versus 55). Conditional analyses revealed GNPDA2 harboring a potential additional independent signal. Moreover, Metabochip-wide discovery analyses revealed two BMI-related loci, BRE (rs116612809, p = 3.6 × 10(-8)) and DHX34 (rs4802349, p = 1.2 × 10(-7)), which were significant when adjustment was made for the total number of SNPs tested across the chip. These results demonstrate that fine mapping in AAs is a powerful approach for both narrowing in on the underlying causal variants in known loci and discovering BMI-related loci.

Genome-wide testing of putative functional exonic variants in relationship with breast and prostate cancer risk in a multiethnic population.

Rare variation in protein coding sequence is poorly captured by GWAS arrays and has been hypothesized to contribute to disease heritability. Using the Illumina HumanExome SNP array, we successfully genotyped 191,032 common and rare non-synonymous, splice site, or nonsense variants in a multiethnic sample of 2,984 breast cancer cases, 4,376 prostate cancer cases, and 7,545 controls. In breast cancer, the strongest associations included either SNPs in or gene burden scores for genes LDLRAD1, SLC19A1, FGFBP3, CASP5, MMAB, SLC16A6, and INS-IGF2. In prostate cancer, one of the most associated SNPs was in the gene GPRC6A (rs2274911, Pro91Ser, OR = 0.88, P = 1.3 × 10(-5)) near to a known risk locus for prostate cancer; other suggestive associations were noted in genes such as F13A1, ANXA4, MANSC1, and GP6. For both breast and prostate cancer, several of the most significant associations involving SNPs or gene burden scores (sum of minor alleles) were noted in genes previously reported to be associated with a cancer-related phenotype. However, only one of the associations (rs145889899 in LDLRAD1, p = 2.5 × 10(-7) only seen in African Americans) for overall breast or prostate cancer risk was statistically significant after correcting for multiple comparisons. In addition to breast and prostate cancer, other cancer-related traits were examined (body mass index, PSA level, and alcohol drinking) with a number of known and potentially novel associations described. In general, these findings do not support there being many protein coding variants of moderate to high risk for breast and prostate cancer with odds ratios over a range that is probably required for protein coding variation to play a truly outstanding role in risk heritability. Very large sample sizes will be required to better define the role of rare and less penetrant coding variation in prostate and breast cancer disease genetics.

Multiple regions within 8q24 independently affect risk for prostate cancer.

After the recent discovery that common genetic variation in 8q24 influences inherited risk of prostate cancer, we genotyped 2,973 SNPs in up to 7,518 men with and without prostate cancer from five populations. We identified seven risk variants, five of them previously undescribed, spanning 430 kb and each independently predicting risk for prostate cancer (P = 7.9 x 10(-19) for the strongest association, and P < 1.5 x 10(-4) for five of the variants, after controlling for each of the others). The variants define common genotypes that span a more than fivefold range of susceptibility to cancer in some populations. None of the prostate cancer risk variants aligns to a known gene or alters the coding sequence of an encoded protein.

Genome-wide association study of endometrial cancer in E2C2.

Endometrial cancer (EC), a neoplasm of the uterine epithelial lining, is the most common gynecological malignancy in developed countries and the fourth most common cancer among US women. Women with a family history of EC have an increased risk for the disease, suggesting that inherited genetic factors play a role. We conducted a two-stage genome-wide association study of Type I EC. Stage 1 included 5,472 women (2,695 cases and 2,777 controls) of European ancestry from seven studies. We selected independent single-nucleotide polymorphisms (SNPs) that displayed the most significant associations with EC in Stage 1 for replication among 17,948 women (4,382 cases and 13,566 controls) in a multiethnic population (African America, Asian, Latina, Hawaiian and European ancestry), from nine studies. Although no novel variants reached genome-wide significance, we replicated previously identified associations with genetic markers near the HNF1B locus. Our findings suggest that larger studies with specific tumor classification are necessary to identify novel genetic polymorphisms associated with EC susceptibility.

Characterizing genetic risk at known prostate cancer susceptibility loci in African Americans.

GWAS of prostate cancer have been remarkably successful in revealing common genetic variants and novel biological pathways that are linked with its etiology. A more complete understanding of inherited susceptibility to prostate cancer in the general population will come from continuing such discovery efforts and from testing known risk alleles in diverse racial and ethnic groups. In this large study of prostate cancer in African American men (3,425 prostate cancer cases and 3,290 controls), we tested 49 risk variants located in 28 genomic regions identified through GWAS in men of European and Asian descent, and we replicated associations (at p≤0.05) with roughly half of these markers. Through fine-mapping, we identified nearby markers in many regions that better define associations in African Americans. At 8q24, we found 9 variants (p≤6×10(-4)) that best capture risk of prostate cancer in African Americans, many of which are more common in men of African than European descent. The markers found to be associated with risk at each locus improved risk modeling in African Americans (per allele OR = 1.17) over the alleles reported in the original GWAS (OR = 1.08). In summary, in this detailed analysis of the prostate cancer risk loci reported from GWAS, we have validated and improved upon markers of risk in some regions that better define the association with prostate cancer in African Americans. Our findings with variants at 8q24 also reinforce the importance of this region as a major risk locus for prostate cancer in men of African ancestry.

A genome-wide association study of breast cancer in women of African ancestry.

Genome-wide association studies (GWAS) in diverse populations are needed to reveal variants that are more common and/or limited to defined populations. We conducted a GWAS of breast cancer in women of African ancestry, with genotyping of >1,000,000 SNPs in 3,153 African American cases and 2,831 controls, and replication testing of the top 66 associations in an additional 3,607 breast cancer cases and 11,330 controls of African ancestry. Two of the 66 SNPs replicated (p < 0.05) in stage 2, which reached statistical significance levels of 10(-6) and 10(-5) in the stage 1 and 2 combined analysis (rs4322600 at chromosome 14q31: OR = 1.18, p = 4.3 × 10(-6); rs10510333 at chromosome 3p26: OR = 1.15, p = 1.5 × 10(-5)). These suggestive risk loci have not been identified in previous GWAS in other populations and will need to be examined in additional samples. Identification of novel risk variants for breast cancer in women of African ancestry will demand testing of a substantially larger set of markers from stage 1 in a larger replication sample.

Genome-wide association study of abdominal MRI-measured visceral fat: The multiethnic cohort adiposity phenotype study.

Few studies have explored the genetic underpinnings of intra-abdominal visceral fat deposition, which varies substantially by sex and race/ethnicity. Among 1,787 participants in the Multiethnic Cohort (MEC)-Adiposity Phenotype Study (MEC-APS), we conducted a genome-wide association study (GWAS) of the percent visceral adiposity tissue (VAT) area out of the overall abdominal area, averaged across L1-L5 (%VAT), measured by abdominal magnetic resonance imaging (MRI). A genome-wide significant signal was found on chromosome 2q14.3 in the sex-combined GWAS (lead variant rs79837492: Beta per effect allele = -4.76; P = 2.62 × 10-8) and in the male-only GWAS (lead variant rs2968545: (Beta = -6.50; P = 1.09 × 10-9), and one suggestive variant was found at 13q12.11 in the female-only GWAS (rs79926925: Beta = 6.95; P = 8.15 × 10-8). The negatively associated variants were most common in European Americans (T allele of rs79837492; 5%) and African Americans (C allele of rs2968545; 5%) and not observed in Japanese Americans, whereas the positively associated variant was most common in Japanese Americans (C allele of rs79926925, 5%), which was all consistent with the racial/ethnic %VAT differences. In a validation step among UK Biobank participants (N = 23,699 of mainly British and Irish ancestry) with MRI-based VAT volume, both rs79837492 (Beta = -0.026, P = 0.019) and rs2968545 (Beta = -0.028, P = 0.010) were significantly associated in men only (n = 11,524). In the MEC-APS, the association between rs79926925 and plasma sex hormone binding globulin levels reached statistical significance in females, but not in males, with adjustment for total adiposity (Beta = -0.24; P = 0.028), on the log scale. Rs79837492 and rs2968545 are located in intron 5 of CNTNAP5, and rs79926925, in an intergenic region between GJB6 and CRYL1. These novel findings differing by sex and racial/ethnic group warrant replication in additional diverse studies with direct visceral fat measurements.

Shared graft-versus-leukemia minor histocompatibility antigens in DISCOVeRY-BMT.

T-cell responses to minor histocompatibility antigens (mHAs) mediate graft-versus-leukemia (GVL) effects and graft-versus-host disease (GVHD) in allogeneic hematopoietic cell transplantation. Therapies that boost T-cell responses improve allogeneic hematopoietic cell transplant (alloHCT) efficacy but are limited by concurrent increases in the incidence and severity of GVHD. mHAs with expression restricted to hematopoietic tissue (GVL mHAs) are attractive targets for driving GVL without causing GVHD. Prior work to identify mHAs has focused on a small set of mHAs or population-level single-nucleotide polymorphism-association studies. We report the discovery of a large set of novel GVL mHAs based on predicted immunogenicity, tissue expression, and degree of sharing among donor-recipient pairs (DRPs) in the DISCOVeRY-BMT data set of 3231 alloHCT DRPs. The total number of predicted mHAs varied by HLA allele, and the total number and number of each class of mHA significantly differed by recipient genomic ancestry group. From the pool of predicted mHAs, we identified the smallest sets of GVL mHAs needed to cover 100% of DRPs with a given HLA allele. We used mass spectrometry to search for high-population frequency mHAs for 3 common HLA alleles. We validated 24 predicted novel GVL mHAs that are found cumulatively within 98.8%, 60.7%, and 78.9% of DRPs within DISCOVeRY-BMT that express HLA-A∗02:01, HLA-B∗35:01, and HLA-C∗07:02, respectively. We confirmed the immunogenicity of an example novel mHA via T-cell coculture with peptide-pulsed dendritic cells. This work demonstrates that the identification of shared mHAs is a feasible and promising technique for expanding mHA-targeting immunotherapeutics.

Germline Sequencing Analysis to Inform Clinical Gene Panel Testing for Aggressive Prostate Cancer.

Importance: Germline gene panel testing is recommended for men with advanced prostate cancer (PCa) or a family history of cancer. While evidence is limited for some genes currently included in panel testing, gene panels are also likely to be incomplete and missing genes that influence PCa risk and aggressive disease. Objective: To identify genes associated with aggressive PCa. Design, Setting, and Participants: A 2-stage exome sequencing case-only genetic association study was conducted including men of European ancestry from 18 international studies. Data analysis was performed from January 2021 to March 2023. Participants were 9185 men with aggressive PCa (including 6033 who died of PCa and 2397 with confirmed metastasis) and 8361 men with nonaggressive PCa. Exposure: Sequencing data were evaluated exome-wide and in a focused investigation of 29 DNA repair pathway and cancer susceptibility genes, many of which are included on gene panels. Main Outcomes and Measures: The primary study outcomes were aggressive (category T4 or both T3 and Gleason score ≥8 tumors, metastatic PCa, or PCa death) vs nonaggressive PCa (category T1 or T2 and Gleason score ≤6 tumors without known recurrence), and metastatic vs nonaggressive PCa. Results: A total of 17 546 men of European ancestry were included in the analyses; mean (SD) age at diagnosis was 65.1 (9.2) years in patients with aggressive PCa and 63.7 (8.0) years in those with nonaggressive disease. The strongest evidence of association with aggressive or metastatic PCa was noted for rare deleterious variants in known PCa risk genes BRCA2 and ATM (P ≤ 1.9 × 10-6), followed by NBN (P = 1.7 × 10-4). This study found nominal evidence (P < .05) of association with rare deleterious variants in MSH2, XRCC2, and MRE11A. Five other genes had evidence of greater risk (OR≥2) but carrier frequency differences between aggressive and nonaggressive PCa were not statistically significant: TP53, RAD51D, BARD1, GEN1, and SLX4. Deleterious variants in these 11 candidate genes were carried by 2.3% of patients with nonaggressive, 5.6% with aggressive, and 7.0% with metastatic PCa. Conclusions and Relevance: The findings of this study provide further support for DNA repair and cancer susceptibility genes to better inform disease management in men with PCa and for extending testing to men with nonaggressive disease, as men carrying deleterious alleles in these genes are likely to develop more advanced disease.

Genome-Wide Association Analyses Identify Variants in IRF4 Associated With Acute Myeloid Leukemia and Myelodysplastic Syndrome Susceptibility.

The role of common genetic variation in susceptibility to acute myeloid leukemia (AML), and myelodysplastic syndrome (MDS), a group of rare clonal hematologic disorders characterized by dysplastic hematopoiesis and high mortality, remains unclear. We performed AML and MDS genome-wide association studies (GWAS) in the DISCOVeRY-BMT cohorts (2,309 cases and 2,814 controls). Association analysis based on subsets (ASSET) was used to conduct a summary statistics SNP-based analysis of MDS and AML subtypes. For each AML and MDS case and control we used PrediXcan to estimate the component of gene expression determined by their genetic profile and correlate this imputed gene expression level with risk of developing disease in a transcriptome-wide association study (TWAS). ASSET identified an increased risk for de novo AML and MDS (OR = 1.38, 95% CI, 1.26-1.51, Pmeta = 2.8 × 10-12) in patients carrying the T allele at s12203592 in Interferon Regulatory Factor 4 (IRF4), a transcription factor which regulates myeloid and lymphoid hematopoietic differentiation. Our TWAS analyses showed increased IRF4 gene expression is associated with increased risk of de novo AML and MDS (OR = 3.90, 95% CI, 2.36-6.44, Pmeta = 1.0 × 10-7). The identification of IRF4 by both GWAS and TWAS contributes valuable insight on the role of genetic variation in AML and MDS susceptibility.

Genome-wide association study of pancreatic fat: The Multiethnic Cohort Adiposity Phenotype Study.

Several studies have found associations between higher pancreatic fat content and adverse health outcomes, such as diabetes and the metabolic syndrome, but investigations into the genetic contributions to pancreatic fat are limited. This genome-wide association study, comprised of 804 participants with MRI-assessed pancreatic fat measurements, was conducted in the ethnically diverse Multiethnic Cohort-Adiposity Phenotype Study (MEC-APS). Two genetic variants reaching genome-wide significance, rs73449607 on chromosome 13q21.2 (Beta = -0.67, P = 4.50x10-8) and rs7996760 on chromosome 6q14 (Beta = -0.90, P = 4.91x10-8) were associated with percent pancreatic fat on the log scale. Rs73449607 was most common in the African American population (13%) and rs79967607 was most common in the European American population (6%). Rs73449607 was also associated with lower risk of type 2 diabetes (OR = 0.95, 95% CI = 0.89-1.00, P = 0.047) in the Population Architecture Genomics and Epidemiology (PAGE) Study and the DIAbetes Genetics Replication and Meta-analysis (DIAGRAM), which included substantial numbers of non-European ancestry participants (53,102 cases and 193,679 controls). Rs73449607 is located in an intergenic region between GSX1 and PLUTO, and rs79967607 is in intron 1 of EPM2A. PLUTO, a lncRNA, regulates transcription of an adjacent gene, PDX1, that controls beta-cell function in the mature pancreas, and EPM2A encodes the protein laforin, which plays a critical role in regulating glycogen production. If validated, these variants may suggest a genetic component for pancreatic fat and a common etiologic link between pancreatic fat and type 2 diabetes.

Novel genetic variants associated with mortality after unrelated donor allogeneic hematopoietic cell transplantation.

BACKGROUND: Identification of non-human leukocyte antigen (HLA) genetic risk factors could improve survival after allogeneic blood or marrow transplant (BMT) through matching at additional loci or individualizing risk prediction. We hypothesized that non-HLA loci contributed significantly to 1-year overall survival (OS), disease related mortality (DRM) or transplant related mortality (TRM) after unrelated donor (URD)BMT. METHODS: We performed a genome-wide association study (GWAS) in 2,887 acute myeloid leukemia (AML), myelodysplastic syndrome (MDS) and acute lymphoblastic leukemia (ALL) patients and their ≥8/8 HLA-matched URDs comprising two independent cohorts treated from 2000-2011. FINDINGS: Using meta-analyses of both cohorts, genome-wide significant associations (p < 5 × 10-8) were identified in: recipient genomes with OS at MBNL1 (rs9990017, HR = 1.4, 95% CI 1.24-1.56, p = 3.3 × 10-8) and donor-recipient genotype mismatch with OS at LINC02774 (rs10927108, HR = 1.34, 95% CI 1.21-1.48, p = 2.0 × 10-8); donor genomes with DRM at PCNX4 (rs79076914, HR = 1.7, 95% CI 1.41-2.05, p = 3.15 × 10-8), LINC01194 (rs79498125, HR = 1.86, 95% CI 1.49-2.31, p = 2.84 × 10-8), ARID5B (rs2167710, HR = 1.5, 95% CI 1.31-1.73, p = 6.9 × 10-9) and CT49 (rs32250, HR = 1.44, 95% CI1.26-1.64, p = 2.6 × 10-8); recipient genomes at PILRB with TRM (rs141591562, HR = 2.33, 95% CI 1.74-3.12, p = 1.26 × 10-8) and donor-recipient genotype mismatch between EPGN and MTHF2DL with TRM (rs75868097, HR = 2.66, 95% CI 1.92-3.58, p = 4.6 × 10-9). Results publicly available at https://fuma.ctglab.nl/browse. INTERPRETATION: These data provide the first evidence that non-HLA common genetic variation at novel loci with biochemical function significantly impacts 1-year URD-BMT survival. Our findings have implications for donor selection, could guide treatment strategies and provide individualized risk prediction after future validation and functional studies. FUNDING: This project was funded by grants from the National Institutes of Health, USA.