Genome-wide association study identifies four pan-ancestry loci for suicidal ideation in the Million Veteran Program.

Suicidal ideation (SI) often precedes and predicts suicide attempt and death, is the most common suicidal phenotype and is over-represented in veterans. The genetic architecture of SI in the absence of suicide attempt (SA) is unknown, yet believed to have distinct and overlapping risk with other suicidal behaviors. We performed the first GWAS of SI without SA in the Million Veteran Program (MVP), identifying 99,814 SI cases from electronic health records without a history of SA or suicide death (SD) and 512,567 controls without SI, SA or SD. GWAS was performed separately in the four largest ancestry groups, controlling for sex, age and genetic substructure. Ancestry-specific results were combined via meta-analysis to identify pan-ancestry loci. Four genome-wide significant (GWS) loci were identified in the pan-ancestry meta-analysis with loci on chromosomes 6 and 9 associated with suicide attempt in an independent sample. Pan-ancestry gene-based analysis identified GWS associations with DRD2, DCC, FBXL19, BCL7C, CTF1, ANNK1, and EXD3. Gene-set analysis implicated synaptic and startle response pathways (q's<0.05). European ancestry (EA) analysis identified GWS loci on chromosomes 6 and 9, as well as GWS gene associations in EXD3, DRD2, and DCC. No other ancestry-specific GWS results were identified, underscoring the need to increase representation of diverse individuals. The genetic correlation of SI and SA within MVP was high (rG = 0.87; p = 1.09e-50), as well as with post-traumatic stress disorder (PTSD; rG = 0.78; p = 1.98e-95) and major depressive disorder (MDD; rG = 0.78; p = 8.33e-83). Conditional analysis on PTSD and MDD attenuated most pan-ancestry and EA GWS signals for SI without SA to nominal significance, with the exception of EXD3 which remained GWS. Our novel findings support a polygenic and complex architecture for SI without SA which is largely shared with SA and overlaps with psychiatric conditions frequently comorbid with suicidal behaviors.

Nontargeted Plasma Proteomic Analysis of Renal Disease and Pulmonary Hypertension in Patients with Sickle Cell Disease.

Sickle cell disease (SCD) is characterized by red blood cell sickling, vaso-occlusion, hemolytic anemia, damage to multiple organ systems, and, as a result, shortened life expectancy. Sickle cell disease nephropathy (SCDN) and pulmonary hypertension (pHTN) are common and frequently co-occurring complications of SCD; both are associated with markedly accelerated mortality. To identify candidate circulating biomarkers of SCDN and pHTN, we used mass spectrometry to quantify the relative abundance of >1000 proteins in plasma samples from 189 adults with SCD from the Outcome Modifying Genes in SCD (OMG-SCD) cohort (ProteomeXchange identifier PXD048716). Forty-four proteins were differentially abundant in SCDN, most significantly cystatin-C and collagen α-1(XVIII) chain (COIA1), and 55 proteins were dysregulated in patients with SCDN and pHTN, most significantly insulin-like growth factor-binding protein 6 (IBP6). Network analysis identified a module of 133 coregulated proteins significantly associated with SCDN, that was enriched for extracellular matrix proteins, insulin-like growth factor binding proteins, cell adhesion proteins, EGF-like calcium binding proteins, and several cadherin family members. Collectively, these data provide a comprehensive understanding of plasma protein changes in SCDN and pHTN which validate numerous studies of chronic kidney disease and suggest shared profiles of protein disruption in kidney dysfunction and pHTN among SCD patients.

Rare and de novo coding variants in chromodomain genes in Chiari I malformation.

Chiari I malformation (CM1), the displacement of the cerebellum through the foramen magnum into the spinal canal, is one of the most common pediatric neurological conditions. Individuals with CM1 can present with neurological symptoms, including severe headaches and sensory or motor deficits, often as a consequence of brainstem compression or syringomyelia (SM). We conducted whole-exome sequencing (WES) on 668 CM1 probands and 232 family members and performed gene-burden and de novo enrichment analyses. A significant enrichment of rare and de novo non-synonymous variants in chromodomain (CHD) genes was observed among individuals with CM1 (combined p = 2.4 × 10-10), including 3 de novo loss-of-function variants in CHD8 (LOF enrichment p = 1.9 × 10-10) and a significant burden of rare transmitted variants in CHD3 (p = 1.8 × 10-6). Overall, individuals with CM1 were found to have significantly increased head circumference (p = 2.6 × 10-9), with many harboring CHD rare variants having macrocephaly. Finally, haploinsufficiency for chd8 in zebrafish led to macrocephaly and posterior hindbrain displacement reminiscent of CM1. These results implicate chromodomain genes and excessive brain growth in CM1 pathogenesis.

Variation and impact of polygenic hematologic traits in monogenic sickle cell disease.

Several of the complications observed in sickle cell disease (SCD) are influenced by variation in hematologic traits (HT), such as fetal hemoglobin (HbF) level and neutrophil count. Previous large-scale genome-wide association studies carried out in largely healthy individuals have identified thousands of variants associated with HT, which have then been used to develop multi-ancestry polygenic trait scores (PTS). Here, we tested whether these PTS associate with HT in SCD patients and if they can improve statistical models associated with SCD-related complications. In 2,056 SCD patients, we found that the PTS predicted less HT variance than in non-SCD individuals of African ancestry. This was particularly striking at the Duffy/DARC locus, where we observed an epistatic interaction between the SCD genotype and the Duffy null variant (rs2814778) that led to a two-fold weaker effect on neutrophil count. PTS for these HT which are measured as part of routine practice were not associated with complications in SCD. In contrast, we found that a simple PTS for HbF that includes only six variants explained a large fraction of the phenotypic variation (20.5-27.1%), associated with acute chest syndrome and stroke risk, and improved the statistical modeling of the vaso-occlusive crisis rate. Using Mendelian randomization, we found that increasing HbF by 4.8% reduces stroke risk by 39% (P=0.0006). Taken together, our results highlight the importance of validating PTS in large diseased populations before proposing their implementation in the context of precision medicine initiatives.

A genome-wide association study of suicide attempts in the million veterans program identifies evidence of pan-ancestry and ancestry-specific risk loci.

To identify pan-ancestry and ancestry-specific loci associated with attempting suicide among veterans, we conducted a genome-wide association study (GWAS) of suicide attempts within a large, multi-ancestry cohort of U.S. veterans enrolled in the Million Veterans Program (MVP). Cases were defined as veterans with a documented history of suicide attempts in the electronic health record (EHR; N = 14,089) and controls were defined as veterans with no documented history of suicidal thoughts or behaviors in the EHR (N = 395,064). GWAS was performed separately in each ancestry group, controlling for sex, age and genetic substructure. Pan-ancestry risk loci were identified through meta-analysis and included two genome-wide significant loci on chromosomes 20 (p = 3.64 × 10-9) and 1 (p = 3.69 × 10-8). A strong pan-ancestry signal at the Dopamine Receptor D2 locus (p = 1.77 × 10-7) was also identified and subsequently replicated in a large, independent international civilian cohort (p = 7.97 × 10-4). Additionally, ancestry-specific genome-wide significant loci were also detected in African-Americans, European-Americans, Asian-Americans, and Hispanic-Americans. Pathway analyses suggested over-representation of many biological pathways with high clinical significance, including oxytocin signaling, glutamatergic synapse, cortisol synthesis and secretion, dopaminergic synapse, and circadian rhythm. These findings confirm that the genetic architecture underlying suicide attempt risk is complex and includes both pan-ancestry and ancestry-specific risk loci. Moreover, pathway analyses suggested many commonly impacted biological pathways that could inform development of improved therapeutics for suicide prevention.

Potential causal role of l-glutamine in sickle cell disease painful crises: A Mendelian randomization analysis.

In a recent clinical trial, the metabolite l-glutamine was shown to reduce painful crises in sickle cell disease (SCD) patients. To support this observation and identify other metabolites implicated in SCD clinical heterogeneity, we profiled 129 metabolites in the plasma of 705 SCD patients. We tested correlations between metabolite levels and six SCD-related complications (painful crises, cholecystectomy, retinopathy, leg ulcer, priapism, aseptic necrosis) or estimated glomerular filtration rate (eGFR), and used Mendelian randomization (MR) to assess causality. We found a potential causal relationship between l-glutamine levels and painful crises (N = 1278, odds ratio (OR) [95% confidence interval] = 0.68 [0.52-0.89], P = 0.0048). In two smaller SCD cohorts (N = 299 and 406), the protective effect of l-glutamine was observed (OR = 0.82 [0.50-1.34]), although the MR result was not significant (P = 0.44). We identified 66 significant correlations between the levels of other metabolites and SCD-related complications or eGFR. We tested these correlations for causality using MR analyses and found no significant causal relationship. The baseline levels of quinolinic acid were associated with prospectively ascertained survival in SCD patients, and this effect was dependent on eGFR. Metabolomics provide a promising approach to prioritize small molecules that may serve as biomarkers or drug targets in SCD.

Rapid decline in estimated glomerular filtration rate in sickle cell anemia: results of a multicenter pooled analysis.

Not available.

Clinical and metabolomic risk factors associated with rapid renal function decline in sickle cell disease.

Sickle cell disease (SCD) nephropathy and lower estimated glomerular filtration rate (eGFR) are risk factors for early mortality. Furthermore, rate of eGFR decline predicts progression to end-stage renal disease in many clinical settings. However, factors predicting renal function decline in SCD are poorly documented. Using clinical, laboratory, genetic, and metabolomic data, we evaluated predictors of renal function decline in a longitudinal cohort of 288 adults (mean age 33.0 years). In 193 subjects with 5-year follow-up data, mean rate of eGFR decline was 2.35 mL/min/1.73 m2 /year, nearly twice that of African American adults overall. Hyperfiltration was prevalent at baseline (61.1%), and 36.8% of subjects experienced rapid eGFR decline (≥3 mL/min/1.73 m2 /year). Severe Hb genotype; proteinuria; higher platelet and reticulocyte counts, and systolic BP; and lower Hb level and BMI were associated with rapid decline. A risk scoring system was created using these 7 variables and was highly predictive of rapid eGFR decline, with odds of rapid decline increasing 1.635-fold for every point increment (P < 0.0001). Rapid eGFR decline was also associated with higher organ system severity score and peak creatinine. Additionally, two metabolites (asymmetric dimethylarginine and quinolinic acid) were associated with rapid decline. Further investigation into longitudinal SCD nephropathy (SCDN) trajectory, early markers of SCDN, and tools for risk stratification should inform interventional studies targeted to slowing GFR decline and improving SCD outcomes.

In vivo Modeling Implicates APOL1 in Nephropathy: Evidence for Dominant Negative Effects and Epistasis under Anemic Stress.

African Americans have a disproportionate risk for developing nephropathy. This disparity has been attributed to coding variants (G1 and G2) in apolipoprotein L1 (APOL1); however, there is little functional evidence supporting the role of this protein in renal function. Here, we combined genetics and in vivo modeling to examine the role of apol1 in glomerular development and pronephric filtration and to test the pathogenic potential of APOL1 G1 and G2. Translational suppression or CRISPR/Cas9 genome editing of apol1 in zebrafish embryos results in podocyte loss and glomerular filtration defects. Complementation of apol1 morphants with wild-type human APOL1 mRNA rescues these defects. However, the APOL1 G1 risk allele does not ameliorate defects caused by apol1 suppression and the pathogenicity is conferred by the cis effect of both individual variants of the G1 risk haplotype (I384M/S342G). In vivo complementation studies of the G2 risk allele also indicate that the variant is deleterious to protein function. Moreover, APOL1 G2, but not G1, expression alone promotes developmental kidney defects, suggesting a possible dominant-negative effect of the altered protein. In sickle cell disease (SCD) patients, we reported previously a genetic interaction between APOL1 and MYH9. Testing this interaction in vivo by co-suppressing both transcripts yielded no additive effects. However, upon genetic or chemical induction of anemia, we observed a significantly exacerbated nephropathy phenotype. Furthermore, concordant with the genetic interaction observed in SCD patients, APOL1 G2 reduces myh9 expression in vivo, suggesting a possible interaction between the altered APOL1 and myh9. Our data indicate a critical role for APOL1 in renal function that is compromised by nephropathy-risk encoding variants. Moreover, our interaction studies indicate that the MYH9 locus is also relevant to the phenotype in a stressed microenvironment and suggest that consideration of the context-dependent functions of both proteins will be required to develop therapeutic paradigms.

The association of single-nucleotide polymorphisms in the oxytocin receptor and G protein-coupled receptor kinase 6 (GRK6) genes with oxytocin dosing requirements and labor outcomes.

BACKGROUND: Oxytocin is a potent uterotonic agent that is widely used for induction and augmentation of labor. Oxytocin has a narrow therapeutic index and the optimal dosing for any individual woman varies widely. OBJECTIVE: The objective of this study was to determine whether genetic variation in the oxytocin receptor (OXTR) or in the gene encoding G protein-coupled receptor kinase 6 (GRK6), which regulates desensitization of the oxytocin receptor, could explain variation in oxytocin dosing and labor outcomes among women being induced near term. STUDY DESIGN: Pregnant women with a singleton gestation residing in Durham County, NC, were prospectively enrolled as part of the Healthy Pregnancy, Healthy Baby cohort study. Those women undergoing an induction of labor at 36 weeks or greater were genotyped for 18 haplotype-tagging single-nucleotide polymorphisms in OXTR and 7 haplotype-tagging single-nucleotide polymorphisms in GRK6 using TaqMan assays. Linear regression was used to examine the relationship between maternal genotype and maximal oxytocin infusion rate, total oxytocin dose received, and duration of labor. Logistic regression was used to test for the association of maternal genotype with mode of delivery. For each outcome, backward selection techniques were utilized to control for important confounding variables and additive genetic models were used. Race/ethnicity was included in all models because of differences in allele frequencies across populations, and Bonferroni correction for multiple testing was used. RESULTS: DNA was available from 482 women undergoing induction of labor at 36 weeks or greater. Eighteen haplotype-tagging single-nucleotide polymorphisms within OXTR and 7 haplotype-tagging single-nucleotide polymorphisms within GRK6 were examined. Five single-nucleotide polymorphisms in OXTR showed nominal significance with maximal infusion rate of oxytocin, and two single-nucleotide polymorphisms in OXTR were associated with total oxytocin dose received. One single-nucleotide polymorphism in OXTR and two single-nucleotide polymorphisms in GRK6 were associated with duration of labor, one of which met the multiple testing threshold (P = .0014, rs2731664 [GRK6], mean duration of labor, 17.7 hours vs 20.2 hours vs 23.5 hours for AA, AC, and CC genotypes, respectively). Three single-nucleotide polymorphisms, two in OXTR and one in GRK6, showed nominal significance with mode of delivery. CONCLUSION: Genetic variation in OXTR and GRK6 is associated with the amount of oxytocin required as well as the duration of labor and risk for cesarean delivery among women undergoing induction of labor near term. With further research, pharmacogenomic approaches may potentially be utilized to develop personalized treatment to improve safety and efficacy outcomes among women undergoing induction of labor.

Discovery and functional annotation of SIX6 variants in primary open-angle glaucoma.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10(-11)). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10(-10)) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

Epigenome-wide association of PTSD from heterogeneous cohorts with a common multi-site analysis pipeline.

Compelling evidence suggests that epigenetic mechanisms such as DNA methylation play a role in stress regulation and in the etiologic basis of stress related disorders such as Post traumatic Stress Disorder (PTSD). Here we describe the purpose and methods of an international consortium that was developed to study the role of epigenetics in PTSD. Inspired by the approach used in the Psychiatric Genomics Consortium, we brought together investigators representing seven cohorts with a collective sample size of N = 1147 that included detailed information on trauma exposure, PTSD symptoms, and genome-wide DNA methylation data. The objective of this consortium is to increase the analytical sample size by pooling data and combining expertise so that DNA methylation patterns associated with PTSD can be identified. Several quality control and analytical pipelines were evaluated for their control of genomic inflation and technical artifacts with a joint analysis procedure established to derive comparable data over the cohorts for meta-analysis. We propose methods to deal with ancestry population stratification and type I error inflation and discuss the advantages and disadvantages of applying robust error estimates. To evaluate our pipeline, we report results from an epigenome-wide association study (EWAS) of age, which is a well-characterized phenotype with known epigenetic associations. Overall, while EWAS are highly complex and subject to similar challenges as genome-wide association studies (GWAS), we demonstrate that an epigenetic meta-analysis with a relatively modest sample size can be well-powered to identify epigenetic associations. Our pipeline can be used as a framework for consortium efforts for EWAS.

Joint eQTL assessment of whole blood and dura mater tissue from individuals with Chiari type I malformation.

BACKGROUND: Expression quantitative trait loci (eQTL) play an important role in the regulation of gene expression. Gene expression levels and eQTLs are expected to vary from tissue to tissue, and therefore multi-tissue analyses are necessary to fully understand complex genetic conditions in humans. Dura mater tissue likely interacts with cranial bone growth and thus may play a role in the etiology of Chiari Type I Malformation (CMI) and related conditions, but it is often inaccessible and its gene expression has not been well studied. A genetic basis to CMI has been established; however, the specific genetic risk factors are not well characterized. RESULTS: We present an assessment of eQTLs for whole blood and dura mater tissue from individuals with CMI. A joint-tissue analysis identified 239 eQTLs in either dura or blood, with 79% of these eQTLs shared by both tissues. Several identified eQTLs were novel and these implicate genes involved in bone development (IPO8, XYLT1, and PRKAR1A), and ribosomal pathways related to marrow and bone dysfunction, as potential candidates in the development of CMI. CONCLUSIONS: Despite strong overall heterogeneity in expression levels between blood and dura, the majority of cis-eQTLs are shared by both tissues. The power to detect shared eQTLs was improved by using an integrative statistical approach. The identified tissue-specific and shared eQTLs provide new insight into the genetic basis for CMI and related conditions.

Hepatic gene expression profiles differentiate presymptomatic patients with mild versus severe nonalcoholic fatty liver disease.

UNLABELLED: Clinicians rely upon the severity of liver fibrosis to segregate patients with well-compensated nonalcoholic fatty liver disease (NAFLD) into subpopulations at high- versus low-risk for eventual liver-related morbidity and mortality. We compared hepatic gene expression profiles in high- and low-risk NAFLD patients to identify processes that distinguish the two groups and hence might be novel biomarkers or treatment targets. Microarray analysis was used to characterize gene expression in percutaneous liver biopsies from low-risk, "mild" NAFLD patients (fibrosis stage 0-1; n = 40) and high-risk, "severe" NAFLD patients (fibrosis stage 3-4; n = 32). Findings were validated in a second, independent cohort and confirmed by real-time polymerase chain reaction and immunohistochemistry (IHC). As a group, patients at risk for bad NAFLD outcomes had significantly worse liver injury and more advanced fibrosis (severe NAFLD) than clinically indistinguishable NAFLD patients with a good prognosis (mild NAFLD). A 64-gene profile reproducibly differentiated severe NAFLD from mild NAFLD, and a 20-gene subset within this profile correlated with NAFLD severity, independent of other factors known to influence NAFLD progression. Multiple genes involved with tissue repair/regeneration and certain metabolism-related genes were induced in severe NAFLD. Ingenuity Pathway Analysis and IHC confirmed deregulation of metabolic and regenerative pathways in severe NAFLD and revealed overlap among the gene expression patterns of severe NAFLD, cardiovascular disease, and cancer. CONCLUSION: By demonstrating specific metabolic and repair pathways that are differentially activated in livers with severe NAFLD, gene profiling identified novel targets that can be exploited to improve diagnosis and treatment of patients who are at greatest risk for NAFLD-related morbidity and mortality.

Gene-centric association study of acute chest syndrome and painful crisis in sickle cell disease patients.

Patients with sickle cell disease (SCD) present with a wide range of clinical complications. Understanding this clinical heterogeneity offers the prospects to tailor the right treatments to the right patients and also guide the development of novel therapies. Several environmental (eg, nutrition) and nonenvironmental (eg, fetal hemoglobin levels, α-thalassemia status) factors are known to modify SCD severity. To find new genetic modifiers of SCD severity, we performed a gene-centric association study in 1514 African American participants from the Cooperative Study of Sickle Cell Disease (CSSCD) for acute chest syndrome (ACS) and painful crisis. From the initial results, we selected 36 single nucleotide polymorphism (SNPs) and genotyped them for replication in 387 independent patients from the CSSCD, 318 SCD patients recruited at Georgia Health Sciences University, and 449 patients from the Duke SCD cohort. In the combined analysis, an association between ACS and rs6141803 reached array-wide significance (P = 4.1 × 10(-7)). This SNP is located 8.2 kilobases upstream of COMMD7, a gene highly expressed in the lung that interacts with nuclear factor-κB signaling. Our results provide new leads to gaining a better understanding of clinical variability in SCD, a "simple" monogenic disease.

Alloimmunization in sickle cell disease: changing antibody specificities and association with chronic pain and decreased survival.

BACKGROUND: Alloimmunization remains a significant complication of transfusion and has been associated with multiple factors, including inflammation, an important pathophysiologic mechanism in sickle cell disease (SCD). We explored whether alloimmunization is associated with disease severity in SCD. STUDY DESIGN AND METHODS: Adult SCD patients were enrolled in a study of outcome-modifying genes in SCD. Historical records of patients with SCD at two participating institutions were reviewed for data on antigen phenotype and alloimmunization. Differences in demographic, clinical, and laboratory findings; end-organ damage; and overall disease severity were then compared between alloimmunized and nonalloimmunized patients. RESULTS: Of 319 patients, 87 (27%) were alloimmunized. Alloantibody specificities differed from those previously described, especially due to the significantly higher frequency of anti-S. Although alloimmunization was not associated with frequency of vasoocclusive episodes, a higher percentage of alloimmunized patients had chronic pain, as defined by daily use of short-acting narcotics (p = 0.006), long-acting narcotics (p = 0.013), or both (p = 0.03). Additionally, alloimmunized patients had poorer survival (hazard ratio, 1.92; p = 0.01) and were more likely to have avascular necrosis (p = 0.024), end-organ damage (p = 0.049), and red blood cell autoantibodies (p < 0.001), even after controlling for the effects of age, sex, and hemoglobin diagnosis. Alloimmunization was not associated with other SCD-related complications, such as acute chest syndrome or stroke. CONCLUSION: Alloimmunization in SCD may be associated with chronic pain, risk of end-organ damage, and shorter survival. These novel findings suggest new directions for the investigation of immune response-mediated pathways common to alloimmunization and chronic pain.

EFFECT OF THE APOE ε4 ALLELE AND COMBAT EXPOSURE ON PTSD AMONG IRAQ/AFGHANISTAN-ERA VETERANS.

BACKGROUND: The apolipoprotein E (APOE) ε4 allele has been implicated in a range of neuropsychiatric conditions. The present research examined if the ε4 allele of the APOE gene moderated the effect of combat exposure on posttraumatic stress disorder (PTSD) among Iraq/Afghanistan-era veterans. METHOD: Participants included 765 non-Hispanic White (NHW) and 859 non-Hispanic Black (NHB) Iraq/Afghanistan-era veterans. A structured interview established psychiatric diagnoses. Combat exposure and PTSD symptom severity were assessed via self-report. RESULTS: The most common lifetime diagnoses were depression (39.2%), PTSD (38.4%), and alcohol dependence (24.38%). After correcting for multiple comparisons, no significant effects were observed on any of the outcomes among the NHW sample; however, within the NHB sample, significant gene × environment (G × E) interactions were observed for lifetime PTSD (P = .0029) and PTSD symptom severity (P = .0009). In each case, the APOE ε4 allele had no effect on the outcomes when combat exposure was low; however, when combat exposure was high, an additive effect was observed such that ε4 homozygotes exposed to high levels of combat reported the highest rates of PTSD (92%) and the worst symptom severity scores on the Davidson Trauma Scale (M = 79.5). CONCLUSIONS: Although preliminary, these findings suggest that the APOE ε4 allele, in conjunction with exposure to high levels of combat exposure, may increase veterans' risk for developing PTSD.