Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Birth defects associated with paternal firefighting in the National Birth Defects Prevention Study.

BACKGROUND: Few studies have evaluated birth defects among children of firefighters. We investigated associations between birth defects and paternal work as a firefighter compared to work in non-firefighting and police officer occupations. METHODS: We analyzed 1997-2011 data from the multi-site case-control National Birth Defects Prevention Study. Cases included fetuses or infants with major structural birth defects and controls included a random sample of live-born infants without major birth defects. Mothers of infants self-reported information about parents' occupations held during pregnancy. We investigated associations between paternal firefighting and birth defect groups using logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (CIs). Referent groups included families reporting fathers working non-firefighting and police officer jobs. RESULTS: Occupational groups included 227 firefighters, 36,285 non-firefighters, and 433 police officers. Twenty-nine birth defects were analyzed. In adjusted analyses, fathers of children with total anomalous pulmonary venous return (TAPVR; OR = 3.1; 95% CI = 1.1-8.7), cleft palate (OR = 1.8; 95% CI = 1.0-3.3), cleft lip (OR = 2.2; 95% CI = 1.2-4.2), and transverse limb deficiency (OR = 2.2; 95% CI = 1.1-4.7) were more likely than fathers of controls to be firefighters, versus non-firefighters. In police-referent analyses, fathers of children with cleft palate were 2.4 times more likely to be firefighters than fathers of controls (95% CI = 1.1-5.4). CONCLUSIONS: Paternal firefighting may be associated with an elevated risk of birth defects in offspring. Additional studies are warranted to replicate these findings. Further research may contribute to a greater understanding of the reproductive health of firefighters and their families for guiding workplace practices.

Exome sequencing identifies genetic variants in anophthalmia and microphthalmia.

Anophthalmia and microphthalmia (A/M) are rare birth defects affecting up to 2 per 10,000 live births. These conditions are manifested by the absence of an eye or reduced eye volumes within the orbit leading to vision loss. Although clinical case series suggest a strong genetic component in A/M, few systematic investigations have been conducted on potential genetic contributions owing to low population prevalence. To overcome this challenge, we utilized DNA samples and data collected as part of the National Birth Defects Prevention Study (NBDPS). The NBDPS employed multi-center ascertainment of infants affected by A/M. We performed exome sequencing on 67 family trios and identified numerous genes affected by rare deleterious nonsense and missense variants in this cohort, including de novo variants. We identified 9 nonsense changes and 86 missense variants that are absent from the reference human population (Genome Aggregation Database), and we suggest that these are high priority candidate genes for A/M. We also performed literature curation, single cell transcriptome comparisons, and molecular pathway analysis on the candidate genes and performed protein structure modeling to determine the potential pathogenic variant consequences on PAX6 in this disease.

Multi-omic biomarker identification and validation for diagnosing warzone-related post-traumatic stress disorder.

Post-traumatic stress disorder (PTSD) impacts many veterans and active duty soldiers, but diagnosis can be problematic due to biases in self-disclosure of symptoms, stigma within military populations, and limitations identifying those at risk. Prior studies suggest that PTSD may be a systemic illness, affecting not just the brain, but the entire body. Therefore, disease signals likely span multiple biological domains, including genes, proteins, cells, tissues, and organism-level physiological changes. Identification of these signals could aid in diagnostics, treatment decision-making, and risk evaluation. In the search for PTSD diagnostic biomarkers, we ascertained over one million molecular, cellular, physiological, and clinical features from three cohorts of male veterans. In a discovery cohort of 83 warzone-related PTSD cases and 82 warzone-exposed controls, we identified a set of 343 candidate biomarkers. These candidate biomarkers were selected from an integrated approach using (1) data-driven methods, including Support Vector Machine with Recursive Feature Elimination and other standard or published methodologies, and (2) hypothesis-driven approaches, using previous genetic studies for polygenic risk, or other PTSD-related literature. After reassessment of ~30% of these participants, we refined this set of markers from 343 to 28, based on their performance and ability to track changes in phenotype over time. The final diagnostic panel of 28 features was validated in an independent cohort (26 cases, 26 controls) with good performance (AUC = 0.80, 81% accuracy, 85% sensitivity, and 77% specificity). The identification and validation of this diverse diagnostic panel represents a powerful and novel approach to improve accuracy and reduce bias in diagnosing combat-related PTSD.

Exome sequencing of child-parent trios with bladder exstrophy: Findings in 26 children.

Bladder exstrophy (BE) is a rare, lower ventral midline defect with the bladder and part of the urethra exposed. The etiology of BE is unknown but thought to be influenced by genetic variation with more recent studies suggesting a role for rare variants. As such, we conducted paired-end exome sequencing in 26 child/mother/father trios. Three children had rare (allele frequency ≤ 0.0001 in several public databases) inherited variants in TSPAN4, one with a loss-of-function variant and two with missense variants. Two children had loss-of-function variants in TUBE1. Four children had rare missense or nonsense variants (one per child) in WNT3, CRKL, MYH9, or LZTR1, genes previously associated with BE. We detected 17 de novo missense variants in 13 children and three de novo loss-of-function variants (AKR1C2, PRRX1, PPM1D) in three children (one per child). We also detected rare compound heterozygous loss-of-function variants in PLCH2 and CLEC4M and rare inherited missense or loss-of-function variants in additional genes applying autosomal recessive (three genes) and X-linked recessive inheritance models (13 genes). Variants in two genes identified may implicate disruption in cell migration (TUBE1) and adhesion (TSPAN4) processes, mechanisms proposed for BE, and provide additional evidence for rare variants in the development of this defect.

Infant Mortality Attributable to Birth Defects - United States, 2003-2017.

Birth defects are a leading cause of infant mortality in the United States, accounting for 20.6% of infant deaths in 2017 (1). Rates of infant mortality attributable to birth defects (IMBD) have generally declined since the 1970s (1-3). U.S. linked birth/infant death data from 2003-2017 were used to assess trends in IMBD. Overall, rates declined 10% during 2003-2017, but decreases varied by maternal and infant characteristics. During 2003-2017, IMBD rates decreased 4% for infants of Hispanic mothers, 11% for infants of non-Hispanic black (black) mothers, and 12% for infants of non-Hispanic white (white) mothers. In 2017, these rates were highest among infants of black mothers (13.3 per 10,000 live births) and were lowest among infants of white mothers (9.9). During 2003-2017, IMBD rates for infants who were born extremely preterm (20-27 completed gestational weeks), full term (39-40 weeks), and late term/postterm (41-44 weeks) declined 20%-29%; rates for moderate (32-33 weeks) and late preterm (34-36 weeks) infants increased 17%. Continued tracking of IMBD rates can help identify areas where efforts to reduce IMBD are needed, such as among infants born to black and Hispanic mothers and those born moderate and late preterm (32-36 weeks).

Maternal occupational exposure to solvents and gastroschisis in offspring - National Birth Defects Prevention Study 1997-2011.

OBJECTIVES: The aim of this study was to assess the association between maternal occupational exposure to solvents and gastroschisis in offspring. METHODS: We used data from the National Birth Defects Prevention Study, a large population-based case-control study of major birth defects conducted in 10 US states from 1997 to 2011. Infants with gastroschisis were ascertained by active birth defects surveillance systems. Control infants without major birth defects were selected from vital records or birth hospital records. Self-reported maternal occupational histories were collected by telephone interview. Industrial hygienists reviewed this information to estimate exposure to aromatic, chlorinated and petroleum-based solvents from 1 month before conception through the first trimester of pregnancy. Cumulative exposure to solvents was estimated for the same period accounting for estimated exposure intensity and frequency, job duration and hours worked per week. ORs and 95% CIs were estimated to assess the association between exposure to any solvents or solvent classes, and gastroschisis risk. RESULTS: Among 879 cases and 7817 controls, the overall prevalence of periconceptional solvent exposure was 7.3% and 7.4%, respectively. Exposure to any solvent versus no exposure to solvents was not associated with gastroschisis after adjusting for maternal age (OR 1.00, 95% CI 0.75 to 1.32), nor was an association noted for solvent classes. There was no exposure-response relationship between estimated cumulative solvent exposure and gastroschisis after adjusting for maternal age. CONCLUSION: Our study found no association between maternal occupational solvent exposure and gastroschisis in offspring. Further research is needed to understand risk factors for gastroschisis.

Kernel machine methods for integrative analysis of genome-wide methylation and genotyping studies.

Many large GWAS consortia are expanding to simultaneously examine the joint role of DNA methylation in addition to genotype in the same subjects. However, integrating information from both data types is challenging. In this paper, we propose a composite kernel machine regression model to test the joint epigenetic and genetic effect. Our approach works at the gene level, which allows for a common unit of analysis across different data types. The model compares the pairwise similarities in the phenotype to the pairwise similarities in the genotype and methylation values; and high correspondence is suggestive of association. A composite kernel is constructed to measure the similarities in the genotype and methylation values between pairs of samples. We demonstrate through simulations and real data applications that the proposed approach can correctly control type I error, and is more robust and powerful than using only the genotype or methylation data in detecting trait-associated genes. We applied our method to investigate the genetic and epigenetic regulation of gene expression in response to stressful life events using data that are collected from the Grady Trauma Project. Within the kernel machine testing framework, our methods allow for heterogeneity in effect sizes, nonlinear, and interactive effects, as well as rapid P-value computation.

Problematic alcohol use associates with sodium channel and clathrin linker 1 (SCLT1) in trauma-exposed populations.

Excessive alcohol use is extremely prevalent in the United States, particularly among trauma-exposed individuals. While several studies have examined genetic influences on alcohol use and related problems, this has not been studied in the context of trauma-exposed populations. We report results from a genome-wide association study of alcohol consumption and associated problems as measured by the alcohol use disorders identification test (AUDIT) in a trauma-exposed cohort. Results indicate a genome-wide significant association between total AUDIT score and rs1433375 [N = 1036, P = 2.61 × 10-8 (dominant model), P = 7.76 × 10-8 (additive model)], an intergenic single-nucleotide polymorphism located 323 kb upstream of the sodium channel and clathrin linker 1 (SCLT1) at 4q28. rs1433375 was also significant in a meta-analysis of two similar, but independent, cohorts (N = 1394, P = 0.0004), the Marine Resiliency Study and Systems Biology PTSD Biomarkers Consortium. Functional analysis indicated that rs1433375 was associated with SCLT1 gene expression and cortical-cerebellar functional connectivity measured via resting state functional magnetic resonance imaging. Together, findings suggest a role for sodium channel regulation and cerebellar functioning in alcohol use behavior. Identifying mechanisms underlying risk for problematic alcohol use in trauma-exposed populations is critical for future treatment and prevention efforts.

Prioritizing individual genetic variants after kernel machine testing using variable selection.

Kernel machine learning methods, such as the SNP-set kernel association test (SKAT), have been widely used to test associations between traits and genetic polymorphisms. In contrast to traditional single-SNP analysis methods, these methods are designed to examine the joint effect of a set of related SNPs (such as a group of SNPs within a gene or a pathway) and are able to identify sets of SNPs that are associated with the trait of interest. However, as with many multi-SNP testing approaches, kernel machine testing can draw conclusion only at the SNP-set level, and does not directly inform on which one(s) of the identified SNP set is actually driving the associations. A recently proposed procedure, KerNel Iterative Feature Extraction (KNIFE), provides a general framework for incorporating variable selection into kernel machine methods. In this article, we focus on quantitative traits and relatively common SNPs, and adapt the KNIFE procedure to genetic association studies and propose an approach to identify driver SNPs after the application of SKAT to gene set analysis. Our approach accommodates several kernels that are widely used in SNP analysis, such as the linear kernel and the Identity by State (IBS) kernel. The proposed approach provides practically useful utilities to prioritize SNPs, and fills the gap between SNP set analysis and biological functional studies. Both simulation studies and real data application are used to demonstrate the proposed approach.

Relationship between Toxoplasma gondii seropositivity and acoustic startle response in an inner-city population.

Toxoplasma gondii (TOXO) is a neuroinvasive protozoan parasite that induces the formation of persistent cysts in mammalian brains. It infects approximately 1.1million people in the United States annually. Latent TOXO infection is implicated in the etiology of psychiatric disorders, especially schizophrenia (SCZ), and has been correlated with modestly impaired cognition. The acoustic startle response (ASR) is a reflex seen in all mammals. It is mediated by a simple subcortical circuit, and provides an indicator of neural function. We previously reported the association of TOXO with slowed acoustic startle latency, an index of neural processing speed, in a sample of schizophrenia and healthy control subjects. The alterations in neurobiology with TOXO latent infection may not be specific to schizophrenia. Therefore we examined TOXO in relation to acoustic startle in an urban, predominately African American, population with mixed psychiatric diagnoses, and healthy controls. Physiological and diagnostic data along with blood samples were collected from 364 outpatients treated at an inner-city hospital. TOXO status was determined with an ELISA assay for TOXO-specific IgG. A discrete titer was calculated based on standard cut-points as an indicator of seropositivity, and the TOXO-specific IgG concentration served as serointensity. A series of linear regression models were used to assess the association of TOXO seropositivity and serointensity with ASR magnitude and latency in models adjusting for demographics and psychiatric diagnoses (PTSD, major depression, schizophrenia, psychosis, substance abuse). ASR magnitude was 11.5% higher in TOXO seropositive subjects compared to seronegative individuals (p=0.01). This effect was more pronounced in models with TOXO serointensity that adjusted for sociodemographic covariates (F=7.41, p=0.0068; F=10.05, p=0.0017), and remained significant when psychiatric diagnoses were stepped into the models. TOXO showed no association with startle latency (t=0.49, p=0.63) in an unadjusted model, nor was TOXO associated with latency in models that included demographic factors. After stepping in individual psychiatric disorders, we found a significant association of latency with a diagnosis of PTSD (F=5.15, p=0.024), but no other psychiatric diagnoses, such that subjects with PTSD had longer startle latency. The mechanism by which TOXO infection is associated with high startle magnitude is not known, but possible mechanisms include TOXO cyst burden in the brain, parasite recrudescence, or molecular mimicry of a host epitope by TOXO. Future studies will focus on the neurobiology underlying the effects of latent TOXO infection as a potential inroad to the development of novel treatment targets for psychiatric disease.

Association Between Infant Mortality Attributable to Birth Defects and Payment Source for Delivery - United States, 2011-2013.

Birth defects are a leading cause of infant mortality in the United States (1), accounting for approximately 20% of infant deaths. The rate of infant mortality attributable to birth defects (IMBD) in the United States in 2014 was 11.9 per 10,000 live births (1). Rates of IMBD differ by race/ethnicity (2), age group at death (2), and gestational age at birth (3). Insurance type is associated with survival among infants with congenital heart defects (CHD) (4). In 2003, a checkbox indicating principal payment source for delivery was added to the U.S. standard birth certificate (5). To assess IMBD by payment source for delivery, CDC analyzed linked U.S. birth/infant death data for 2011-2013 from states that adopted the 2003 revision of the birth certificate. The results indicated that IMBD rates for preterm (<37 weeks of gestation) and term (≥37 weeks) infants whose deliveries were covered by Medicaid were higher during the neonatal (<28 days) and postneonatal (≥28 days to <1 year) periods compared with infants whose deliveries were covered by private insurance. Similar differences in postneonatal mortality were observed for the three most common categories of birth defects listed as a cause of death: central nervous system (CNS) defects, CHD, and chromosomal abnormalities. Strategies to ensure quality of care and access to care might reduce the difference between deliveries covered by Medicaid and those covered by private insurance.

PACAP receptor gene polymorphism impacts fear responses in the amygdala and hippocampus.

We have recently found higher circulating levels of pituitary adenylate cyclase-activating polypeptide (PACAP) associated with posttraumatic stress disorder (PTSD) symptoms in a highly traumatized cohort of women but not men. Furthermore, a single nucleotide polymorphism in the PACAP receptor gene ADCYAP1R1, adenylate cyclase activating polypeptide 1 receptor type 1, was associated with individual differences in PTSD symptoms and psychophysiological markers of fear and anxiety. The current study outlines an investigation of individual differences in brain function associated with ADCYAP1R1 genotype. Forty-nine women who had experienced moderate to high levels of lifetime trauma participated in a functional MRI task involving passive viewing of threatening and neutral face stimuli. Analyses focused on the amygdala and hippocampus, regions that play central roles in the pathophysiology of PTSD and are known to have high densities of PACAP receptors. The risk genotype was associated with increased reactivity of the amygdala and hippocampus to threat stimuli and decreased functional connectivity between the amygdala and hippocampus. The findings indicate that the PACAP system modulates medial temporal lobe function in humans. Individual differences in ADCYAP1R1 genotype may contribute to dysregulated fear circuitry known to play a central role in PTSD and other anxiety disorders.

Kernel Approach for Modeling Interaction Effects in Genetic Association Studies of Complex Quantitative Traits.

The etiology of complex traits likely involves the effects of genetic and environmental factors, along with complicated interaction effects between them. Consequently, there has been interest in applying genetic association tests of complex traits that account for potential modification of the genetic effect in the presence of an environmental factor. One can perform such an analysis using a joint test of gene and gene-environment interaction. An optimal joint test would be one that remains powerful under a variety of models ranging from those of strong gene-environment interaction effect to those of little or no gene-environment interaction effect. To fill this demand, we have extended a kernel machine based approach for association mapping of multiple SNPs to consider joint tests of gene and gene-environment interaction. The kernel-based approach for joint testing is promising, because it incorporates linkage disequilibrium information from multiple SNPs simultaneously in analysis and permits flexible modeling of interaction effects. Using simulated data, we show that our kernel machine approach typically outperforms the traditional joint test under strong gene-environment interaction models and further outperforms the traditional main-effect association test under models of weak or no gene-environment interaction effects. We illustrate our test using genome-wide association data from the Grady Trauma Project, a cohort of highly traumatized, at-risk individuals, which has previously been investigated for interaction effects.

GENOME-WIDE ASSOCIATION STUDY (GWAS) AND GENOME-WIDE BY ENVIRONMENT INTERACTION STUDY (GWEIS) OF DEPRESSIVE SYMPTOMS IN AFRICAN AMERICAN AND HISPANIC/LATINA WOMEN.

BACKGROUND: Genome-wide association studies (GWAS) have made little progress in identifying variants linked to depression. We hypothesized that examining depressive symptoms and considering gene-environment interaction (GxE) might improve efficiency for gene discovery. We therefore conducted a GWAS and genome-wide by environment interaction study (GWEIS) of depressive symptoms. METHODS: Using data from the SHARe cohort of the Women's Health Initiative, comprising African Americans (n = 7,179) and Hispanics/Latinas (n = 3,138), we examined genetic main effects and GxE with stressful life events and social support. We also conducted a heritability analysis using genome-wide complex trait analysis (GCTA). Replication was attempted in four independent cohorts. RESULTS: No SNPs achieved genome-wide significance for main effects in either discovery sample. The top signals in African Americans were rs73531535 (located 20 kb from GPR139, P = 5.75 × 10(-8) ) and rs75407252 (intronic to CACNA2D3, P = 6.99 × 10(-7) ). In Hispanics/Latinas, the top signals were rs2532087 (located 27 kb from CD38, P = 2.44 × 10(-7) ) and rs4542757 (intronic to DCC, P = 7.31 × 10(-7) ). In the GEWIS with stressful life events, one interaction signal was genome-wide significant in African Americans (rs4652467; P = 4.10 × 10(-10) ; located 14 kb from CEP350). This interaction was not observed in a smaller replication cohort. Although heritability estimates for depressive symptoms and stressful life events were each less than 10%, they were strongly genetically correlated (rG = 0.95), suggesting that common variation underlying self-reported depressive symptoms and stressful life event exposure, though modest on their own, were highly overlapping in this sample. CONCLUSIONS: Our results underscore the need for larger samples, more GEWIS, and greater investigation into genetic and environmental determinants of depressive symptoms in minorities.

Oxytocin Receptor Genetic and Epigenetic Variations: Association With Child Abuse and Adult Psychiatric Symptoms.

Childhood abuse can alter biological systems and increase risk for adult psychopathology. Epigenetic mechanisms, alterations in DNA structure that regulate the gene expression, are a potential mechanism underlying this risk. While abuse associates with methylation of certain genes, particularly those in the stress response system, no study to date has evaluated abuse and methylation of the oxytocin receptor (OXTR). However, studies support a role for OXTR in the link between abuse and adverse adult outcomes, showing that abuse can confer greater risk for psychiatric symptoms in those with specific OXTR genotypes. This study therefore sought to (a) assess the role of epigenetics in the link between abuse and psychopathology and (b) begin to integrate the genetic and epigenetic literature by exploring associations between OXTR genotypes and DNA CpG methylation. Data on 18 OXTR CpG sites, 44 single nucleotide polymorphisms, childhood abuse, and adult depression and anxiety symptoms were assessed in 393 African American adults (age = 41 ± 12.8 years). Overall, 68% of genotypes were associated with methylation of nearby CpG sites, with a subset surviving multiple test correction. Child abuse associated with higher methylation of two CpG sites yet did not survive correction or serve as a mediator of psychopathology. However, abuse interacted with CpG methylation to predict psychopathology. These findings suggest a role for OXTR in understanding the influence of early environments on adult psychiatric symptoms.

Stress and Bronchodilator Response in Children with Asthma.

RATIONALE: Stress is associated with asthma morbidity in Puerto Ricans (PRs), who have reduced bronchodilator response (BDR). OBJECTIVES: To examine whether stress and/or a gene regulating anxiety (ADCYAP1R1) is associated with BDR in PR and non-PR children with asthma. METHODS: This was a cross-sectional study of stress and BDR (percent change in FEV1 after BD) in 234 PRs ages 9-14 years with asthma. We assessed child stress using the Checklist of Children's Distress Symptoms, and maternal stress using the Perceived Stress Scale. Replication analyses were conducted in two cohorts. Polymorphisms in ADCYAP1R1 were genotyped in our study and six replication studies. Multivariable models of stress and BDR were adjusted for age, sex, income, environmental tobacco smoke, and use of inhaled corticosteroids. MEASUREMENTS AND MAIN RESULTS: High child stress was associated with reduced BDR in three cohorts. PR children who were highly stressed (upper quartile, Checklist of Children's Distress Symptoms) and whose mothers had high stress (upper quartile, Perceived Stress Scale) had a BDR that was 10.2% (95% confidence interval, 6.1-14.2%) lower than children who had neither high stress nor a highly stressed mother. A polymorphism in ADCYAP1R1 (rs34548976) was associated with reduced BDR. This single-nucleotide polymorphism is associated with reduced expression of the gene for the ß2-adrenergic receptor (ADRB2) in CD4(+) lymphocytes of subjects with asthma, and it affects brain connectivity of the amygdala and the insula (a biomarker of anxiety). CONCLUSIONS: High child stress and an ADCYAP1R1 single-nucleotide polymorphism are associated with reduced BDR in children with asthma. This is likely caused by down-regulation of ADRB2 in highly stressed children.

Accounting for population stratification in DNA methylation studies.

DNA methylation is an important epigenetic mechanism that has been linked to complex diseases and is of great interest to researchers as a potential link between genome, environment, and disease. As the scale of DNA methylation association studies approaches that of genome-wide association studies, issues such as population stratification will need to be addressed. It is well-documented that failure to adjust for population stratification can lead to false positives in genetic association studies, but population stratification is often unaccounted for in DNA methylation studies. Here, we propose several approaches to correct for population stratification using principal components (PCs) from different subsets of genome-wide methylation data. We first illustrate the potential for confounding due to population stratification by demonstrating widespread associations between DNA methylation and race in 388 individuals (365 African American and 23 Caucasian). We subsequently evaluate the performance of our PC-based approaches and other methods in adjusting for confounding due to population stratification. Our simulations show that (1) all of the methods considered are effective at removing inflation due to population stratification, and (2) maximum power can be obtained with single-nucleotide polymorphism (SNP)-based PCs, followed by methylation-based PCs, which outperform both surrogate variable analysis and genomic control. Among our different approaches to computing methylation-based PCs, we find that PCs based on CpG sites chosen for their potential to proxy nearby SNPs can provide a powerful and computationally efficient approach to adjust for population stratification in DNA methylation studies when genome-wide SNP data are unavailable.

A genome-wide identified risk variant for PTSD is a methylation quantitative trait locus and confers decreased cortical activation to fearful faces.

Genetic factors appear to be highly relevant to predicting differential risk for the development of post-traumatic stress disorder (PTSD). In a discovery sample, we conducted a genome-wide association study (GWAS) for PTSD using a small military cohort (Systems Biology PTSD Biomarkers Consortium; SBPBC, N = 147) that was designed as a case-controlled sample of highly exposed, recently returning veterans with and without combat-related PTSD. A genome-wide significant single nucleotide polymorphism (SNP), rs717947, at chromosome 4p15 (N = 147, ß = 31.34, P = 1.28 × 10(-8) ) was found to associate with the gold-standard diagnostic measure for PTSD (the Clinician Administered PTSD Scale). We conducted replication and follow-up studies in an external sample, a larger urban community cohort (Grady Trauma Project, GTP, N = 2006), to determine the robustness and putative functionality of this risk variant. In the GTP replication sample, SNP rs717947 associated with PTSD diagnosis in females (N = 2006, P = 0.005), but not males. SNP rs717947 was also found to be a methylation quantitative trait locus (meQTL) in the GTP replication sample (N = 157, P = 0.002). Further, the risk allele of rs717947 was associated with decreased medial and dorsolateral cortical activation to fearful faces (N = 53, P < 0.05) in the GTP replication sample. These data identify a genome-wide significant polymorphism conferring risk for PTSD, which was associated with differential epigenetic regulation and with differential cortical responses to fear in a replication sample. These results may provide new insight into understanding genetic and epigenetic regulation of PTSD and intermediate phenotypes that contribute to this disorder.

Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type.

BACKGROUND: Individual genotypes at specific loci can result in different patterns of DNA methylation. These methylation quantitative trait loci (meQTLs) influence methylation across extended genomic regions and may underlie direct SNP associations or gene-environment interactions. We hypothesized that the detection of meQTLs varies with ancestral population, developmental stage, and tissue type. We explored this by analyzing seven datasets that varied by ancestry (African American vs. Caucasian), developmental stage (neonate vs. adult), and tissue type (blood vs. four regions of postmortem brain) with genome-wide DNA methylation and SNP data. We tested for meQTLs by constructing linear regression models of methylation levels at each CpG site on SNP genotypes within 50 kb under an additive model controlling for multiple tests. RESULTS: Most meQTLs mapped to intronic regions, although a limited number appeared to occur in synonymous or nonsynonymous coding SNPs. We saw significant overlap of meQTLs between ancestral groups, developmental stages, and tissue types, with the highest rates of overlap within the four brain regions. Compared with a random group of SNPs with comparable frequencies, meQTLs were more likely to be 1) represented among the most associated SNPs in the WTCCC bipolar disorder results and 2) located in microRNA binding sites. CONCLUSIONS: These data give us insight into how SNPs impact gene regulation and support the notion that peripheral blood may be a reliable correlate of physiological processes in other tissues.