Andrade, Mariza de

Abstract Abstract 1148 Background: Factor V Leiden (F5rs6025) has been suggested as a disease-susceptibility SNP for deep vein thrombosis (DVT) but not PE. Given the significantly reduced survival after PE, identifying individuals at risk for PE is important for targeting prophylaxis and improving survival. Objectives: To identify PE disease-susceptibility genes. Methods: We performed in silicoGWAS analyses of VTE cases using genotype data imputed to ∼2.5 million SNPs from adults with objectively-diagnosed PE ± DVT (n=745), and DVT only (n=758). We tested these SNPs for an association with PE ± DVT after adjusting for age, sex, stroke and/or myocardial infarction, and USA state of residence using unconditional logistic regression. Results: Of all SNPs, three were most strongly associated with PE ± DVT: one at BDH1 on chromosome (chr) 3q29 (rs1309873, OR=1.56, p=4.02E-07) and two at SLC39A10 on chr 2q32.3 (rs7578216, OR=2.55, p=4.57E-07; and rs4589778, OR=2.49, p=7.52E-07). Individually adjusting for F5 rs6025 (chr 1), F2 rs1799963 (prothrombin G20210A; chr 11) and ABO rs8176719 (ABO blood type O allele; chr 9) did not materially affect these associations. None of these three SNPs associated with PE were intragenic. BDH1 encodes for 3-hydroxybutyrate dehydrogenase 1, a key mitochondrial enzyme in fatty acid catabolism, and SLC39A10encodes for solute carrier family 39, member 10, an important protein for zinc transport. Conclusions: Chromosomes 3q29 (BDH1) and 2q32.2 (SLC39A10) may harbor PE disease-susceptibility gene regions (genes). These findings require independent confirmation. If confirmed, additional studies addressing the biologic mechanisms through which these genes and/or gene regions operate to increase PE risk also will be required. Disclosures: No relevant conflicts of interest to declare.

Abstract Abstract 150 Background: While bivariate interactions between Factor V Leiden, prothrombin G20201A and hereditary deficiency of antithrombin, protein C and protein S compound VTE risk, whether other gene-gene interactions are associated with VTE is largely unknown. Objective: To test gene-gene interactions for an association with VTE. Methods: Cases (n=1004) were Mayo Clinic European-American patients of non-Hispanic ancestry with objectively-diagnosed VTE in the absence of active cancer, venous catheter or antiphospholipid antibodies. Controls (n=1066) were Mayo Clinic outpatients without VTE who were frequency-matched on case age, gender, race, MI/stroke status and state of residence. We selected candidate genes relevant to the anticoagulant, procoagulant, fibrinolytic and innate immunity pathways, focusing on platelet, monocyte, neutrophil and endothelial cell agonists, receptors, ligands, signal transduction and adhesion molecules, granule contents and effectors; plasma proteases and inhibitors; matrix metalloproteases; inflammatory cytokines and receptors; estrogen, progesterone and androgen receptors; co-regulators and enzymes related to estrogen metabolism; important enzymes for catechol, homocysteine, thromboxane A2 and prostacyclin biosynthesis and metabolism; and HMG CoA reductase. For these genes (n=780), we selected all non-synonymous coding single nucleotide polymorphisms (SNPs) with minor allele frequency ≥0.5%; the remaining SNPs were selected using an LD tagging algorithm (Carlson et al. AJHG 2004); 558 ancestry-informative markers were also included (total n=13,237 SNPs). Leukocyte genomic DNA was genotyped using a custom Illumina Infinium iSelect chemistry and platform, including appropriate controls. For these analyses, all pairwise SNP-SNP interactions for 12,313 SNPs were tested using logistic regression. In addition, we tested all Factor V Leiden -, prothrombin G20210A - and ABO non-O blood group - SNP interactions. Results: The mean ± SD case and control ages were 55.3 ± 16.4 and 56.5 ± 15.9 years, respectively, and 51% were female. Analysis of ancestry-informative markers gave no evidence of population stratification. Among almost 76 million pairwise interactions tested, 504, 44 and 8 SNP-SNP interactions reached p-values ≤ 1E-5, 1E-6, and 1E-7 respectively; none reached the Bonferroni statistical significance threshold (≤ 6.6E-10). The gene-gene pairwise interactions with p-values ≤ 1E-7 included: Leptin receptor - Estrogen receptor 1, α1 Antichymotrypsin - β2 glycoprotein 1, CRADD(caspace) - B1 Bradykinin receptor, ILF10 - Prolactin receptor and GFRA2(TGFβ neurotrophic factor receptor) - Bactericidal permeability-increasing protein. Gene-gene pairwise interactions with the lowest p-value for Factor V Leiden -, prothrombin G20210A - and ABO non-group O - were Protein kinase C beta (2.1E-5), von Willebrand Factor (7.6E-4) and Phospholipase C gamma-2 (6.7E-4), respectively. Conclusion: Pairwise gene-gene interaction analyses suggest potential associations between VTE and novel genes and gene pathways. These potential associations require confirmation in future replication studies. Disclosures: No relevant conflicts of interest to declare.

In this work related to familial aggregation of familial venous thromboembolism, the investigators report genomic and transcriptomic association of 16 novel susceptibility loci for venous thromboembolism.

Abstract Lung cancer (LC) is the primary cause of cancer-related deaths in the United States. Whereas smoking and other environmental factors strongly increase LC risk, multiple genetic variants also contribute to risk in smokers. Furthermore, among smokers, some families have been identified with an abnormally high prevalence of LC, suggesting that unknown genetic factors can greatly increase LC risk in smokers. The typically short survival after LC diagnosis impedes collection of detailed genotypic information on any single large family pedigree, impairing the identification of putative high-risk factors. Therefore, the Genetic Epidemiology of LC Consortium has collected epidemiological and genetic data from a number of families with high numbers of LC cases from eight different sites across the US. In this study, we have obtained whole exome sequences (WES) from 290 members of 28 families, including 66 LC cases. We used a gene-based approach to allow for the possibility that different families may contain different variants of the same gene. Variants were filtered for i) allele frequency, ii) functional effect using combined annotation-dependent depletion (CADD), and, iii) affecting a gene with either a known or suspected role in cancer. We further selected variants based on their segregation in family pedigrees in a pattern consistent with a large effect on LC risk. Candidate LC risk genes were then identified as those represented in at least two families by the same or different variants. We further culled the list of genes by requiring the presence of at least one rare, functional variant enriched in the WES of 1060 cases relative to 899 controls from the Transdisciplinary Research on Cancer of the Lung consortium. This analysis narrowed our results to two genes, one being E2A, a member of the E family of bHLH transcription factors. Whereas loss-of-function mutations in E2A drive lymphoid cancers, the E2A protein also participates in an oncogenic heterodimer with TWIST1 that promotes the epithelial-mesenchymal transition and is implicated in multiple cancer types. Furthermore, the E2A/TWIST1 heterodimer is the primary TWIST1-containing complex implicated in oncogenesis, and silencing of E2A in KRAS-mutant non-small cell lung cancer (NSCLC) cell lines results in oncogene-stimulated senescence and apoptosis. Our data identified three distinct E2A variants present in all 10 sequenced LC cases in the 5 families in which those variants are found. Finally, two of these E2A variants are located only 57 nucleotides from each other, immediately adjacent to sequences encoding a transcription activation domain, suggesting that both variants alter the same specific protein function. These data identify E2A variants as likely high risk factors for LC in smokers and validate our general approach for identifying genetic factors with a large impact on LC risk. Citation Format: Claudio W. Pikielny, Anthony M. Musolf, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, yafang Li, Ming You, Richard Wilson, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Ambrose I. Granizo-Mackenzie, Yanhong Liu, Ramaswamy Govindan, James McKay, Rayjean Hung, John K. Field, David C. Christiani, Joan E. Bailey-Wilson, Christopher I. Amos. Familial studies identify variants in the E2A transcription factor as putative risk factors for lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr LB-053.

Abstract Lung cancer (LC) kills more people in the United States each year than any other cancer. While it is well known that a variety of environmental factors (particularly tobacco smoke) strongly increase the risk of LC, there are multiple associated genetic variants with small contributions to risk. High aggregation of LC within rare individual families suggest that there are high-risk genetic variants as well. However, these genetic risk factors for LC are under studied due to the rapid fatality of LC. We studied 28 highly aggregated extended high-risk familial lung cancer (HRFLC) families collected from eight different sites across the US. Whole exome sequencing was performed on 290 individuals from these families to identify potential risk variants for HRFLC using genetic linkage analysis. Quality control was performed on the sequence data, filtering on parameters such as read depth, genotype quality, missingness, and Mendelian inconsistencies. Identity-by-descent (IBD) values were also calculated to verify correct familial relationships. Quality control procedures left approximately 400,000 SNVs and indels for analysis.Parametric two-point linkage analysis was performed assuming an autosomal dominant mode of inheritance. Disease allele frequency was set to 1% with a penetrance of 80% for carriers and 1% phenocopy rate. While we did not identify any genome-wide significant variants across the 28 families, multiple suggestive variants were identified. The largest cluster of suggestive variants was located at 14q32 in the CATSPERB gene. Given the likely locus heterogeneity in LC (combined with the lack of power in some families), it is not surprising that none of the variants were significant across the families; looking at the individual family results proved more informative. Long haplotypes linked to LC risk were identified in multiple families. These long runs of positive linkages, which have little to no negative linkage evidence across them, are characteristic of true linkage signals in these types of analysis. Two of the most interesting linked regions were at 7p36.1 and 4q21.23-28.23. The 7p signal (observed in a single family) was genome-wide suggestive and located within the SSPO gene. SSPO has been implicated in breast and skin cancer (melanoma); it is a novel lung cancer signal. The 4q linkage (again observed in a single family) covers a large chunk of 4q and contains multiple potential candidate genes, however, the best candidate gene is PTPN13, a gene implicated in lung cancer but never in familial lung cancer. We are currently evaluating the individual family results of several other pedigrees and plan to perform additional analyses to confirm these linkages. Citation Format: Anthony Musolf, Haiming Sun, Bilal A. Moiz, Claudio W. Pikielny, Mariza de Andrade, Diptasri Mandal, Colette Gaba, Ping Yang, Yafang Li, Ming You, Richard K. Wilson, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Ramaswamy Govindan, Joan E. Bailey-Wilson. Familial lung cancer exhibits multiple novel linked haplotypes within pedigrees [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 4176.

Abstract Objective: We studied the long-term risk of depressive and anxiety symptoms in women who underwent bilateral oophorectomy before menopause. Design: We conducted a cohort study among all women residing in Olmsted County, MN, who underwent bilateral oophorectomy before the onset of menopause for a noncancer indication from 1950 through 1987. Each member of the bilateral oophorectomy cohort was matched by age with a referent woman from the same population who had not undergone an oophorectomy. In total, we studied 666 women with bilateral oophorectomy and 673 referent women. Women were followed for a median of 24 years, and depressive and anxiety symptoms were assessed using a structured questionnaire via a direct or proxy telephone interview performed from 2001 through 2006. Results: Women who underwent bilateral oophorectomy before the onset of menopause had an increased risk of depressive symptoms diagnosed by a physician (hazard ratio = 1.54, 95% CI: 1.04-2.26, adjusted for age, education, and type of interview) and of anxiety symptoms (adjusted hazard ratio = 2.29, 95% CI: 1.33-3.95) compared with referent women. The findings remained consistent after excluding depressive or anxiety symptoms that first occurred within 10 years after oophorectomy. The associations were greater with younger age at oophorectomy but did not vary across indications for the oophorectomy. In addition, treatment with estrogen to age 50 years in women who underwent bilateral oophorectomy at younger ages did not modify the risk. Conclusions: Bilateral oophorectomy performed before the onset of menopause is associated with an increased long-term risk of depressive and anxiety symptoms.

Abstract Lung cancer is the deadliest cancer in the United States, contributing approximately 25% of all cancer deaths. Lung cancer risk is well-documented to increase in response to environmental factors, particularly tobacco smoking. As one might expect for a complex trait, there is also a significant genetic risk component in lung cancer, though it is not very well studied. We obtained whole exome sequencing (WES) genotype data from Washington University on 204 subjects from 25 extended families. These individuals were recruited from families with a history of lung cancer and previous analyses showed these 25 families be informative. The purpose of this study is to identify potential risk variants for lung cancer by performing genetic linkage analysis. Quality control was performed on the sequence data, filtering on parameters such as depth (less than 10), genotype quality (less than 10), missingness, and Mendelian inconsistencies. Identity-by-descent (IBD) values were also calculated to verify correct familial relationships. Quality control procedures left approximately 500,000 SNVs and indels for analysis. We performed two-point parametric linkage analysis assuming an autosomal dominant mode of inheritance with a disease allele frequency of 1%, a 10% penetrance for carriers and a 1% penetrance for non-carriers. Two discrete sets of linkage analyses were performed. One was a variant-based analysis, which evaluated linkage between the phenotype and individual SNVs or indels. The second was gene-based analysis, which created a multi-allelic pseudomarker corresponding to a gene from haplotypes of rare variants (minor allele frequency <= 0.01) located within that particular gene. Two-point linkage analysis was then performed on the pseudomarkers. While the variant-based analysis did not identify any genome-wide significant results, several were identified by the gene-based analysis. The highest HLOD scores were both greater than the genome-wide significance level of HLOD = 3.3 and were located on two regions of chromosome 1q: 1q42.13-43 and 1q21.2-21.1. In both regions, the significant HLOD scores clustered around known cancer genes. At the 1q42.13-43 region, signals centered on the cancer-implicated genes OBSCN and RYR2 (also a known mesothelioma gene), while at 1q21.2-21.1 the signals centered on five genes in the neuroblastoma breakpoint family (NBPF) genes, a cluster of recently duplicated genes. NBPF genes have previously been implicated in a variety of different cancers including lung cancer. We are currently performing additional analyses to corroborate the significant results, as well as examining the individual families to determine which families are driving the significant signals. Citation Format: Anthony M. Musolf, Haiming Sun, Bilal A. Moiz, Diptasri Mandal, Mariza de Andrade, Colette Gaba, Ping Yang, Yafang Li, Ming You, Elena Y. Kupert, Marshall W. Anderson, Ann G. Schwartz, Susan M. Pinney, Christopher I. Amos, Joan E. Bailey-Wilson. Whole exome sequencing identifies significantly linked regions on multiple chromosomes in families with a history of lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3276.

Abstract Objective: Response to menopausal hormone therapy (MHT) shows individual variation. SLCO1B1 encodes the OATP1B1 transporter expressed in the liver that transports many endogenous substances, including estrone sulfate, from the blood into hepatocytes. This study evaluated the relationship between genetic variation in SLCO1B1 and response to MHT in women enrolled in the Kronos Early Estrogen Prevention Study (KEEPS) at Mayo Clinic, Rochester, MN. Methods: KEEPS participants were randomized to oral conjugated equine estrogen (n = 33, oCEE), transdermal 17β-estradiol (n = 33, tE2), or placebo (n = 34) for 48 months. Menopausal symptoms (hot flashes, night sweats, insomnia, palpitations) were self-reported before treatment and at 48 months. Estrone (E1), E2, and sulfated conjugates (E1S, E2S) were measured using high-performance liquid chromatography-tandem mass spectrometry. SLCO1B1 rs4149056 (c.521T>C, p.Val174Ala) was genotyped using a TaqMan assay. Results: After adjusting for treatment, there was a significant association between the SLCO1B1 rs4149056 TT genotype (encoding normal function transporter) and lower E1S, E1S/E1, and E2S (P = 0.032, 0.010, and 0.008, respectively) compared with women who were heterozygous (TC) or homozygous (CC) for the reduced function allele. The interactions between genotype, treatment, and E2S concentration were stronger in women assigned to tE2 (P = 0.013) than the women taking oCEE (P = 0.056). Among women assigned to active treatment, women with the CT genotype showed a significantly greater decrease in night sweats (P = 0.041) than those with the TT genotype. Conclusions: Individual variation in sulfated estrogens is explained, in part, by genetic variation in SLCO1B1. Bioavailability of sulfated estrogens may contribute to relief of night sweats.