Rare variants in CAPN2 increase risk for isolated hypoplastic left heart syndrome

Hypoplastic left heart syndrome (HLHS) is a severe congenital heart defect (CHD) characterized by hypoplasia of the left ventricle and aorta along with stenosis or atresia of the aortic and mitral valves. HLHS represents only ∼4%-8% of all CHDs but accounts for ∼25% of deaths. HLHS is an isolated defect (i.e., iHLHS) in 70% of families, the vast majority of which are simplex. Despite intense investigation, the genetic basis of iHLHS remains largely unknown. We performed exome sequencing on 331 families with iHLHS aggregated from four independent cohorts. A Mendelian-model-based analysis demonstrated that iHLHS was not due to single, large-effect alleles in genes previously reported to underlie iHLHS or CHD in >90% of families in this cohort. Gene-based association testing identified increased risk for iHLHS associated with variation in CAPN2 (p = 1.8 × 10-5), encoding a protein involved in functional adhesion. Functional validation studies in a vertebrate animal model (Xenopus laevis) confirmed CAPN2 is essential for cardiac ventricle morphogenesis and that in vivo loss of calpain function causes hypoplastic ventricle phenotypes and suggest that human CAPN2707C>T and CAPN21112C>T variants, each found in multiple individuals with iHLHS, are hypomorphic alleles. Collectively, our findings show that iHLHS is typically not a Mendelian condition, demonstrate that CAPN2 variants increase risk of iHLHS, and identify a novel pathway involved in HLHS pathogenesis.

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.

Exome sequencing identifies variants in infants with sacral agenesis

BACKGROUND

Sacral agenesis (SA) consists of partial or complete absence of the caudal end of the spine and often presents with additional birth defects. Several studies have examined gene variants for syndromic forms of SA, but only one has examined exomes of children with non-syndromic SA.

METHODS

Using buccal cell specimens from families of children with non-syndromic SA, exomes of 28 child-parent trios (eight with and 20 without a maternal diagnosis of pregestational diabetes) and two child-father duos (neither with diagnosis of maternal pregestational diabetes) were exome sequenced.

RESULTS

Three children had heterozygous missense variants in ID1 (Inhibitor of DNA Binding 1), with CADD scores >20 (top 1% of deleterious variants in the genome); two children inherited the variant from their fathers and one from the child's mother. Rare missense variants were also detected in PDZD2 (PDZ Domain Containing 2; N = 1) and SPTBN5 (Spectrin Beta, Non-erythrocytic 5; N = 2), two genes previously suggested to be associated with SA etiology. Examination of variants with autosomal recessive and X-linked recessive inheritance identified five and two missense variants, respectively. Compound heterozygous variants were identified in several genes. In addition, 12 de novo variants were identified, all in different genes in different children.

CONCLUSIONS

To our knowledge, this is the first study reporting a possible association between ID1 and non-syndromic SA. Although maternal pregestational diabetes has been strongly associated with SA, the missense variants in ID1 identified in two of three children were paternally inherited. These findings add to the knowledge of gene variants associated with non-syndromic SA and provide data for future studies.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

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.

Examining Individual and Synergistic Contributions of PTSD and Genetics to Blood Pressure: A Trans-Ethnic Meta-Analysis

Growing research suggests that posttraumatic stress disorder (PTSD) may be a risk factor for poor cardiovascular health, and yet our understanding of who might be at greatest risk of adverse cardiovascular outcomes after trauma is limited. In this study, we conducted the first examination of the individual and synergistic contributions of PTSD symptoms and blood pressure genetics to continuous blood pressure levels. We harnessed the power of the Psychiatric Genomics Consortium-PTSD Physical Health Working Group and investigated these associations across 11 studies of 72,224 trauma-exposed individuals of European (n = 70,870) and African (n = 1,354) ancestry. Genetic contributions to blood pressure were modeled via polygenic scores (PGS) for systolic blood pressure (SBP) and diastolic blood pressure (DBP) that were derived from a prior trans-ethnic blood pressure genome-wide association study (GWAS). Results of trans-ethnic meta-analyses revealed significant main effects of the PGS on blood pressure levels [SBP: β = 2.83, standard error (SE) = 0.06, p < 1E-20; DBP: β = 1.32, SE = 0.04, p < 1E-20]. Significant main effects of PTSD symptoms were also detected for SBP and DBP in trans-ethnic meta-analyses, though there was significant heterogeneity in these results. When including data from the largest contributing study - United Kingdom Biobank - PTSD symptoms were negatively associated with SBP levels (β = -1.46, SE = 0.44, p = 9.8E-4) and positively associated with DBP levels (β = 0.70, SE = 0.26, p = 8.1E-3). However, when excluding the United Kingdom Biobank cohort in trans-ethnic meta-analyses, there was a nominally significant positive association between PTSD symptoms and SBP levels (β = 2.81, SE = 1.13, p = 0.01); no significant association was observed for DBP (β = 0.43, SE = 0.78, p = 0.58). Blood pressure PGS did not significantly moderate the associations between PTSD symptoms and blood pressure levels in meta-analyses. Additional research is needed to better understand the extent to which PTSD is associated with high blood pressure and how genetic as well as contextual factors may play a role in influencing cardiovascular risk.

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.

International meta-analysis of PTSD genome-wide association studies identifies sex- and ancestry-specific genetic risk loci

The risk of posttraumatic stress disorder (PTSD) following trauma is heritable, but robust common variants have yet to be identified. In a multi-ethnic cohort including over 30,000 PTSD cases and 170,000 controls we conduct a genome-wide association study of PTSD. We demonstrate SNP-based heritability estimates of 5-20%, varying by sex. Three genome-wide significant loci are identified, 2 in European and 1 in African-ancestry analyses. Analyses stratified by sex implicate 3 additional loci in men. Along with other novel genes and non-coding RNAs, a Parkinson's disease gene involved in dopamine regulation, PARK2, is associated with PTSD. Finally, we demonstrate that polygenic risk for PTSD is significantly predictive of re-experiencing symptoms in the Million Veteran Program dataset, although specific loci did not replicate. These results demonstrate the role of genetic variation in the biology of risk for PTSD and highlight the necessity of conducting sex-stratified analyses and expanding GWAS beyond European ancestry populations.

Genome-wide identification of DNA methylation QTLs in whole blood highlights pathways for cardiovascular disease

Identifying methylation quantitative trait loci (meQTLs) and integrating them with disease-associated variants from genome-wide association studies (GWAS) may illuminate functional mechanisms underlying genetic variant-disease associations. Here, we perform GWAS of >415 thousand CpG methylation sites in whole blood from 4170 individuals and map 4.7 million cis- and 630 thousand trans-meQTL variants targeting >120 thousand CpGs. Independent replication is performed in 1347 participants from two studies. By linking cis-meQTL variants with GWAS results for cardiovascular disease (CVD) traits, we identify 92 putatively causal CpGs for CVD traits by Mendelian randomization analysis. Further integrating gene expression data reveals evidence of cis CpG-transcript pairs causally linked to CVD. In addition, we identify 22 trans-meQTL hotspots each targeting more than 30 CpGs and find that trans-meQTL hotspots appear to act in cis on expression of nearby transcriptional regulatory genes. Our findings provide a powerful meQTL resource and shed light on DNA methylation involvement in human diseases.

Differentially methylated CpGs can inform on disease mechanisms and be useful as biomarkers. Here, the authors perform GWAS for DNA methylation in whole blood, cis- and trans-meQTL mapping, followed by Mendelian randomization analysis that links meQTLs with cardiovascular diseases.

Exome sequencing of family trios from the National Birth Defects Prevention Study: Tapping into a rich resource of genetic and environmental data

BACKGROUND

The National Birth Defects Prevention Study (NBDPS) is a multisite, population-based, case-control study of genetic and nongenetic risk factors for major structural birth defects. Eligible women had a pregnancy affected by a birth defect or a liveborn child without a birth defect between 1997 and 2011. They were invited to complete a telephone interview to collect pregnancy exposure data and were mailed buccal cell collection kits to collect specimens from themselves, their child (if living), and their child's father. Over 23,000 families representing more than 30 major structural birth defects provided DNA specimens.

METHODS

To evaluate their utility for exome sequencing (ES), specimens from 20 children with colonic atresia were studied. Evaluations were conducted on specimens collected using cytobrushes stored and transported in open versus closed packaging, on native genomic DNA (gDNA) versus whole genome amplified (WGA) products and on a library preparation protocol adapted to low amounts of DNA.

RESULTS

The DNA extracted from brushes in open packaging yielded higher quality sequence data than DNA from brushes in closed packaging. Quality metrics of sequenced gDNA were consistently higher than metrics from corresponding WGA products and were consistently high when using a low input protocol.

CONCLUSIONS

This proof-of-principle study established conditions under which ES can be applied to NBDPS specimens. Successful sequencing of exomes from well-characterized NBDPS families indicated that this unique collection can be used to investigate the roles of genetic variation and gene-environment interaction effects in birth defect etiologies, providing a valuable resource for birth defect researchers.

Powerful and Efficient Strategies for Genetic Association Testing of Symptom and Questionnaire Data in Psychiatric Genetic Studies

Genetic studies of psychiatric disorders often deal with phenotypes that are not directly measurable. Instead, researchers rely on multivariate symptom data from questionnaires and surveys like the PTSD Symptom Scale (PSS) and Beck Depression Inventory (BDI) to indirectly assess a latent phenotype of interest. Researchers subsequently collapse such multivariate questionnaire data into a univariate outcome to represent a surrogate for the latent phenotype. However, when a causal variant is only associated with a subset of collapsed symptoms, the effect will be challenging to detect using the univariate outcome. We describe a more powerful strategy for genetic association testing in this situation that jointly analyzes the original multivariate symptom data collectively using a statistical framework that compares similarity in multivariate symptom-scale data from questionnaires to similarity in common genetic variants across a gene. We use simulated data to demonstrate this strategy provides substantially increased power over standard approaches that collapse questionnaire data into a single surrogate outcome. We also illustrate our approach using GWAS data from the Grady Trauma Project and identify genes associated with BDI not identified using standard univariate techniques. The approach is computationally efficient, scales to genome-wide studies, and is applicable to correlated symptom data of arbitrary dimension.

Risk of gastroschisis with maternal genitourinary infections: the US National birth defects prevention study 1997-2011

OBJECTIVE

To assess the association between occurrence and timing of maternal self-reported genitourinary tract infection (urinary tract infections [UTIs] and/or sexually transmitted infection [STI]) and risk for gastroschisis in the offspring.

DESIGN

Population-based case-control study.

SETTING

National Birth Defects Prevention Study, a multisite study in the USA.

PARTICIPANTS

Mothers of 1366 gastroschisis cases and 11 238 healthy controls.

MAIN OUTCOME MEASURES

Crude and adjusted ORs (aORs) with 95% CIs.

RESULTS

Genitourinary infections were frequent in case (19.3%) and control women (9.9%) during the periconceptional period (defined as 3 months prior to 3 months after conception). UTI and/or STI in the periconceptional period were associated with similarly increased risks for gastroschisis (aOR 1.5, 95% CI 1.3 to 1.8; aOR 1.6, 95% CI 1.2 to 2.3, respectively). The risk was increased with a UTI before (aOR 2.5; 95% CI 1.4 to 4.5) or after (aOR 1.7; 95% CI 1.1 to 2.6) conception only among women ≥25 years of age. The risk was highest among women <20 years of age with an STI before conception (aOR 3.6; 95% CI 1.5 to 8.4) and in women ≥25 years of age, the risk was similar for before (aOR 2.9; 95% CI 1.0 to 8.5) and after (aOR 2.8; 95% CI 1.3 to 6.1) conception. A specific STI pathogen was reported in 89.3% (50/56) of cases and 84.3% (162/191) of controls with Chlamydia trachomatis the most common (25/50 cases, 50%; 58/162 controls, 36%) and highest among women <20 years of age (16/25 cases, 64%; 22/33 controls, 67%).

CONCLUSIONS

UTI and/or STI were associated with an increased risk for gastroschisis, with the strength of the association varying by maternal age and timing of infection.

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.

Cross-disorder risk gene CACNA1C differentially modulates susceptibility to psychiatric disorders during development and adulthood

Single-nucleotide polymorphisms (SNPs) in CACNA1C, the α1C subunit of the voltage-gated L-type calcium channel Cav1.2, rank among the most consistent and replicable genetics findings in psychiatry and have been associated with schizophrenia, bipolar disorder and major depression. However, genetic variants of complex diseases often only confer a marginal increase in disease risk, which is additionally influenced by the environment. Here we show that embryonic deletion of Cacna1c in forebrain glutamatergic neurons promotes the manifestation of endophenotypes related to psychiatric disorders including cognitive decline, impaired synaptic plasticity, reduced sociability, hyperactivity and increased anxiety. Additional analyses revealed that depletion of Cacna1c during embryonic development also increases the susceptibility to chronic stress, which suggest that Cav1.2 interacts with the environment to shape disease vulnerability. Remarkably, this was not observed when Cacna1c was deleted in glutamatergic neurons during adulthood, where the later deletion even improved cognitive flexibility, strengthened synaptic plasticity and induced stress resilience. In a parallel gene × environment design in humans, we additionally demonstrate that SNPs in CACNA1C significantly interact with adverse life events to alter the risk to develop symptoms of psychiatric disorders. Overall, our results further validate Cacna1c as a cross-disorder risk gene in mice and humans, and additionally suggest a differential role for Cav1.2 during development and adulthood in shaping cognition, sociability, emotional behavior and stress susceptibility. This may prompt the consideration for pharmacological manipulation of Cav1.2 in neuropsychiatric disorders with developmental and/or stress-related origins.

Largest GWAS of PTSD (N=20 070) yields genetic overlap with schizophrenia and sex differences in heritability

The Psychiatric Genomics Consortium-Posttraumatic Stress Disorder group (PGC-PTSD) combined genome-wide case-control molecular genetic data across 11 multiethnic studies to quantify PTSD heritability, to examine potential shared genetic risk with schizophrenia, bipolar disorder, and major depressive disorder and to identify risk loci for PTSD. Examining 20 730 individuals, we report a molecular genetics-based heritability estimate (h2SNP) for European-American females of 29% that is similar to h2SNP for schizophrenia and is substantially higher than h2SNP in European-American males (estimate not distinguishable from zero). We found strong evidence of overlapping genetic risk between PTSD and schizophrenia along with more modest evidence of overlap with bipolar and major depressive disorder. No single-nucleotide polymorphisms (SNPs) exceeded genome-wide significance in the transethnic (overall) meta-analysis and we do not replicate previously reported associations. Still, SNP-level summary statistics made available here afford the best-available molecular genetic index of PTSD-for both European- and African-American individuals-and can be used in polygenic risk prediction and genetic correlation studies of diverse phenotypes. Publication of summary statistics for ∼10 000 African Americans contributes to the broader goal of increased ancestral diversity in genomic data resources. In sum, the results demonstrate genetic influences on the development of PTSD, identify shared genetic risk between PTSD and other psychiatric disorders and highlight the importance of multiethnic/racial samples. As has been the case with schizophrenia and other complex genetic disorders, larger sample sizes are needed to identify specific risk loci.

Estrogen-dependent association of HDAC4 with fear in female mice and women with PTSD

Women are at increased risk of developing post-traumatic stress disorder (PTSD) following a traumatic event. Recent studies suggest that this may be mediated, in part, by circulating estrogen levels. This study evaluated the hypothesis that individual variation in response to estrogen levels contributes to fear regulation and PTSD risk in women. We evaluated DNA methylation from blood of female participants in the Grady Trauma Project and found that serum estradiol levels associates with DNA methylation across the genome. For genes expressed in blood, we examined the association between each CpG site and PTSD diagnosis using linear models that adjusted for cell proportions and age. After multiple test correction, PTSD associated with methylation of CpG sites in the HDAC4 gene, which encodes histone deacetylase 4, and is involved in long-term memory formation and behavior. DNA methylation of HDAC4 CpG sites were tagged by a nearby single-nucleotide polymorphism (rs7570903), which also associated with HDAC4 expression, fear-potentiated startle and resting-state functional connectivity of the amygdala in traumatized humans. Using auditory Pavlovian fear conditioning in a rodent model, we examined the regulation of Hdac4 in the amygdala of ovariectomized (OVX) female mice. Hdac4 messenger RNA levels were higher in the amygdala 2 h after tone-shock presentations, compared with OVX-homecage control females. In naturally cycling females, tone-shock presentations increased Hdac4 expression relative to homecage controls for metestrous (low estrogen) but not the proestrous (high estrogen) group. Together, these results support an estrogenic influence of HDAC4 regulation and expression that may contribute to PTSD in women.

A putative causal relationship between genetically determined female body shape and posttraumatic stress disorder

Background

The nature and underlying mechanisms of the observed increased vulnerability to posttraumatic stress disorder (PTSD) in women are unclear.

Methods

We investigated the genetic overlap of PTSD with anthropometric traits and reproductive behaviors and functions in women. The analysis was conducted using female-specific summary statistics from large genome-wide association studies (GWAS) and a cohort of 3577 European American women (966 PTSD cases and 2611 trauma-exposed controls). We applied a high-resolution polygenic score approach and Mendelian randomization analysis to investigate genetic correlations and causal relationships.

Results

We observed an inverse association of PTSD with genetically determined anthropometric traits related to body shape, independent of body mass index (BMI). The top association was related to BMI-adjusted waist circumference (WC_adj; R = –0.079, P < 0.001, Q = 0.011). We estimated a relative decrease of 64.6% (95% confidence interval = 27.5–82.7) in the risk of PTSD per 1-SD increase in WC_adj. MR-Egger regression intercept analysis showed no evidence of pleiotropic effects in this association (P_pleiotropy = 0.979). We also observed associations of genetically determined WC_adj with age at first sexual intercourse and number of sexual partners (P = 0.013 and P < 0.001, respectively).

Conclusions

There is a putative causal relationship between genetically determined female body shape and PTSD, which could be mediated by evolutionary mechanisms involved in human sexual behaviors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0491-4) contains supplementary material, which is available to authorized users.

Genome-wide association study of positive emotion identifies a genetic variant and a role for microRNAs

Positive affect denotes a state of pleasurable engagement with the environment eliciting positive emotion such as contentment, enthusiasm or happiness. Positive affect is associated with favorable psychological, physical and economic outcomes in many longitudinal studies. With a heritability of ⩽64%, positive affect is substantially influenced by genetic factors; however, our understanding of genetic pathways underlying individual differences in positive affect is still limited. Here, through a genome-wide association study of positive affect in African-American participants, we identify a single-nucleotide polymorphism, rs322931, significantly associated with positive affect at P<5 × 10^-8, and replicate this association in another cohort. Furthermore, we show that the minor allele of rs322931 predicts expression of microRNAs miR-181a and miR-181b in human brain and blood, greater nucleus accumbens reactivity to positive emotional stimuli and enhanced fear inhibition. Prior studies have suggested that miR-181a is part of the reward neurocircuitry. Taken together, we identify a novel genetic variant for further elucidation of genetic underpinning of positive affect that mediates positive emotionality potentially via the nucleus accumbens and miR-181.

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.

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.

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.

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.

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.

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.

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.

Association of CRP genetic variation and CRP level with elevated PTSD symptoms and physiological responses in a civilian population with high levels of trauma

OBJECTIVE

Increased systemic inflammation is associated with stress-related psychopathology. Specifically, levels of the proinflammatory marker C-reactive protein (CRP) are elevated in individuals with posttraumatic stress disorder (PTSD). Furthermore, single-nucleotide polymorphisms (SNPs) in the CRP gene are associated with CRP level, risk for cardiovascular disease, and obesity. The authors examined whether polymorphisms within the CRP gene and increased CRP levels are associated with PTSD symptoms and fear physiology in a civilian population with high levels of trauma.

METHOD

Cross-sectional data and DNA samples were collected from 2,698 individuals recruited from an inner-city public hospital that serves a primarily African American, low-socioeconomic-status population. A subgroup of 187 participants participated in further interviews, testing, and physiological measures; of these, 135 were assessed using the fear-potentiated startle paradigm to assess fear-related phenotypes of PTSD.

RESULTS

One SNP within the CRP gene, rs1130864, was significantly associated with increased PTSD symptoms (N=2,692), including "being overly alert" as the most significant individual symptom (N=2,698). Additionally, CRP genotype was associated with the odds of PTSD diagnosis (N=2,692). This SNP was also associated with increased CRP level (N=137), and high CRP levels (>3 mg/L) were positively associated with PTSD symptoms (N=187) and fear-potentiated startle to a safety signal (N=135).

CONCLUSIONS

Together, these data indicate that genetic variability in the CRP gene is associated with serum CRP level and PTSD symptom severity, including that of hyperarousal symptoms. Elevated CRP levels were also associated with exacerbated fear-related psychophysiology and PTSD symptom ratings and diagnosis. These findings suggest a potential mechanism by which an increased proinflammatory state may lead to heightened PTSD symptoms.

Correcting systematic inflation in genetic association tests that consider interaction effects: application to a genome-wide association study of posttraumatic stress disorder

IMPORTANCE

Genetic association studies of psychiatric outcomes often consider interactions with environmental exposures and, in particular, apply tests that jointly consider gene and gene-environment interaction effects for analysis. Using a genome-wide association study (GWAS) of posttraumatic stress disorder (PTSD), we report that heteroscedasticity (defined as variability in outcome that differs by the value of the environmental exposure) can invalidate traditional joint tests of gene and gene-environment interaction.

OBJECTIVES

To identify the cause of bias in traditional joint tests of gene and gene-environment interaction in a PTSD GWAS and determine whether proposed robust joint tests are insensitive to this problem.

DESIGN, SETTING, AND PARTICIPANTS

The PTSD GWAS data set consisted of 3359 individuals (978 men and 2381 women) from the Grady Trauma Project (GTP), a cohort study from Atlanta, Georgia. The GTP performed genome-wide genotyping of participants and collected environmental exposures using the Childhood Trauma Questionnaire and Trauma Experiences Inventory.

MAIN OUTCOMES AND MEASURES

We performed joint interaction testing of the Beck Depression Inventory and modified PTSD Symptom Scale in the GTP GWAS. We assessed systematic bias in our interaction analyses using quantile-quantile plots and genome-wide inflation factors.

RESULTS

Application of the traditional joint interaction test to the GTP GWAS yielded systematic inflation across different outcomes and environmental exposures (inflation-factor estimates ranging from 1.07 to 1.21), whereas application of the robust joint test to the same data set yielded no such inflation (inflation-factor estimates ranging from 1.01 to 1.02). Simulated data further revealed that the robust joint test is valid in different heteroscedasticity models, whereas the traditional joint test is invalid. The robust joint test also has power similar to the traditional joint test when heteroscedasticity is not an issue.

CONCLUSIONS AND RELEVANCE

We believe the robust joint test should be used in candidate-gene studies and GWASs of psychiatric outcomes that consider environmental interactions. To make the procedure useful for applied investigators, we created a software tool that can be called from the popular PLINK package for analysis.

Interaction of the ADRB2 gene polymorphism with childhood trauma in predicting adult symptoms of posttraumatic stress disorder

IMPORTANCE

Posttraumatic stress disorder (PTSD), while highly prevalent (7.6% over a lifetime), develops only in a subset of trauma-exposed individuals. Genetic risk factors in interaction with trauma exposure have been implicated in PTSD vulnerability.

OBJECTIVE

To examine the association of 3755 candidate gene single-nucleotide polymorphisms with PTSD development in interaction with a history of childhood trauma.

DESIGN, SETTING, AND PARTICIPANTS

Genetic association study in an Ohio National Guard longitudinal cohort (n = 810) of predominantly male soldiers of European ancestry, with replication in an independent Grady Trauma Project (Atlanta, Georgia) cohort (n = 2083) of predominantly female African American civilians.

MAIN OUTCOMES AND MEASURES

Continuous measures of PTSD severity, with a modified (interview) PTSD checklist in the discovery cohort and the PTSD Symptom Scale in the replication cohort.

RESULTS

Controlling for the level of lifetime adult trauma exposure, we identified the novel association of a single-nucleotide polymorphism within the promoter region of the ADRB2 (Online Mendelian Inheritance in Man 109690) gene with PTSD symptoms in interaction with childhood trauma (rs2400707, P = 1.02 × 10-5, significant after correction for multiple comparisons). The rs2400707 A allele was associated with relative resilience to childhood adversity. An rs2400707 × childhood trauma interaction predicting adult PTSD symptoms was replicated in the independent predominantly female African American cohort.

CONCLUSIONS AND RELEVANCE

Altered adrenergic and noradrenergic function has been long believed to have a key etiologic role in PTSD development; however, direct evidence of this link has been missing. The rs2400707 polymorphism has been linked to function of the adrenergic system, but, to our knowledge, this is the first study to date linking the ADRB2 gene to PTSD or any psychiatric disorders. These findings have important implications for PTSD etiology, chronic pain, and stress-related comorbidity, as well as for both primary prevention and treatment strategies.

Follow-up and extension of a prior genome-wide association study of posttraumatic stress disorder: gene × environment associations and structural magnetic resonance imaging in a highly traumatized African-American civilian population

Follow-up and Extension of a Prior Genome-wide Association Study of Posttraumatic Stress Disorder: Gene × Environment Associations and Structural Magnetic Resonance Imaging in a Highly Traumatized African-American Civilian Population

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.

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.

Epigenomic association analysis identifies smoking-related DNA methylation sites in African Americans

Cigarette smoking is an environmental risk factor for many chronic diseases, and disease risk can often be managed by smoking control. Smoking can induce cellular and molecular changes, including epigenetic modification, but the short- and long-term epigenetic modifications caused by cigarette smoking at the gene level have not been well understood. Recent studies have identified smoking-related DNA methylation (DNAm) sites in Caucasians. To determine whether the same DNAm sites associate with smoking in African Americans, and to identify novel smoking-related DNAm sites, we conducted a methylome-wide association study of cigarette smoking using a discovery sample of 972 African Americans, and a replication sample of 239 African Americans with two array-based methods. Among 15 DNAm sites significantly associated with smoking after correction for multiple testing in our discovery sample, 5 DNAm sites are replicated in an independent cohort, and 14 sites in the replication sample have effects in the same direction as in the discovery sample. The top two smoking-related DNAm sites in F2RL3 (factor II receptor-like 3) and GPR15 (G-protein-coupled receptor 15) observed in African Americans are consistent with previous findings in Caucasians. The associations between the replicated DNAm sites and smoking remain significant after adjusting for genetic background. Despite the distinct genetic background between African Americans and Caucasians, the DNAm from the two ethnic groups shares common associations with cigarette smoking, which suggests a common molecular mechanism of epigenetic modification influenced by environmental exposure.

Sex dependent influence of a functional polymorphism in steroid 5-α-reductase type 2 (SRD5A2) on post-traumatic stress symptoms

A non-synonymous, single nucleotide polymorphism (SNP) in the gene coding for steroid 5-α-reductase type 2 (SRD5A2) is associated with reduced conversion of testosterone to dihydrotestosterone (DHT). Because SRD5A2 participates in the regulation of testosterone and cortisol metabolism, hormones shown to be dysregulated in patients with PTSD, we examined whether the V89L variant (rs523349) influences risk for post-traumatic stress disorder (PTSD). Study participants (N = 1,443) were traumatized African-American patients of low socioeconomic status with high rates of lifetime trauma exposure recruited from the primary care clinics of a large, urban hospital. PTSD symptoms were measured with the post-traumatic stress symptom scale (PSS). Subjects were genotyped for the V89L variant (rs523349) of SRD5A2. We initially found a significant sex-dependent effect of genotype in male but not female subjects on symptoms. Associations with PTSD symptoms were confirmed using a separate internal replication sample with identical methods of data analysis, followed by pooled analysis of the combined samples (N = 1,443, sex × genotype interaction P < 0.002; males: n = 536, P < 0.001). These data support the hypothesis that functional variation within SRD5A2 influences, in a sex-specific way, the severity of post-traumatic stress symptoms and risk for diagnosis of PTSD.

ADCYAP1R1 genotype associates with post-traumatic stress symptoms in highly traumatized African-American females

Pituitary adenylate cyclase-activating polypeptide (PACAP) and its receptor (PAC1) play a critical role in biological processes that mediate stress response and have been implicated in psychological outcome following trauma. Our previous work [Ressler et al. (2011); Nature 470:492-497] demonstrated that a variant, rs2267735, in the gene encoding PAC1 (ADCYAP1R1) is associated with post-traumatic stress disorder (PTSD) in a primarily African-American cohort of highly traumatized females. We sought to extend and replicate our previous finding in a similarly trauma-exposed, replicate sample of 1,160 African-American adult male and female patients. Self-reported psychiatric measures were collected, and DNA was obtained for genetic analysis. Using linear regression models to test for association with PTSD symptom severity under an additive (allelic) model, we found a genotype × trauma interaction in females (P < 0.001), but not males (P > 0.1); however, there was no main effect of genotype as in our previous study. The observed interaction suggests a genetic association that increases with the degree of trauma exposure in females only. This interaction remained significant in females, but not males, after controlling for age (P < 0.001), income (P < 0.01), past substance abuse (P < 0.001), depression severity (P = 0.02), or child abuse (P < 0.0005), and all five combined (P = 0.01). No significant effects of genotype (or interactions) were found when modeling depression severity when controlling for comorbid PTSD symptom severity (P > 0.1), demonstrating the relative specificity of this variant for PTSD symptoms. A meta-analysis with the previously reported African-American samples revealed a strong association between PTSD symptom severity and the interaction between trauma and genotype in females (N = 1424, P < 0.0001).

Dysregulated biomarkers induce distinct pathways in preterm birth

OBJECTIVE

To document racial disparity in biomarker concentrations in maternal/fetal plasma and amniotic fluid between African Americans and European Americans with spontaneous preterm birth (PTB; cases) and normal term birth (controls), and their contribution to distinct pathophysiological pathways of PTB.

DESIGN

Nested case-control study.

SETTING

The Perinatal Research Center, Nashville, Tennessee, USA.

SAMPLE

Maternal and fetal plasma and amniotic fluid samples were collected from 105 cases (59 African American and 46 European American) and 86 controls (40 African American and 46 European American).

METHODS

Thirty-six biomarkers were analysed using the protein microarray approach.

MAIN OUTCOME MEASURES

Differences in biomarker concentrations between cases and controls of different races in maternal, fetal and intra-amniotic compartments, and the risk of PTB. Dysregulated biomarker-induced PTB pathways associated with PTB in each race were determined using ingenuity pathway analysis (IPA).

RESULTS

Racial disparity was observed in biomarker concentrations in each compartment between cases and controls: amniotic fluid, IL8 and MIP1α differed between case and controls in European Americans, whereas ANGPT2, Eotaxin, ICAM-1, IL-1β, IL1RA, RANTES and TNFα differed between case and controls in African Americans. In both races the FAS ligand, MCP-3 and TNFR-I differed between cases and controls. For fetal plasma, ANGPT2, Eotaxin, FGF basic, ICAM-1, IGF-I, IL10, IL-1β, IL2, IP10 KGF, MCP-3, MIP1α, PDGF-BB, TGFα, TGFβ1, TIMP1, TNFα, TNFR-I, TNFR-II and VEGF differed between cases and controls in European Americans, whereas only MMP7 differed between cases and controls in African Americans. IL-8 differed between cases and controls in both races. For maternal plasma, IL1RA, MMP7 and VEGF differed between cases and controls in European Americans, whereas ANGPT2, FGF basic, IL-1β, IL5, IL6R, KGF, MCP-3, MIP1α, TIMP1 and TNFα differed between cases and controls in African Americans. ANG, IL8 and TNFR-I differed between cases and controls in both races.

CONCLUSIONS

We conclude that: (1) biomarker concentrations in maternal, fetal and intra-amniotic compartments differ between cases and controls; (2) there is racial disparity in the biomarker profile in each of the compartments; (3) substantial numbers of dysregulated fetal plasma biomarkers contribute to PTB in European Americans, whereas maternal plasma biomarkers contribute to PTB in African Americans; and (4) both inflammation and haematological functions are associated with PTB in European Americans, but maternal proinflammatory changes dominate PTB in African Americans. Biomarker analyses document racial disparity and the distinct pathophysiological contributions from different compartments that can determine pregnancy outcome.

Socially modulated cell proliferation is independent of gonadal steroid hormones in the brain of the adult green treefrog (Hyla cinerea)

Gonadal steroid hormones have been shown to influence adult neurogenesis in addition to their well-defined role in regulating social behavior. Adult neurogenesis consists of several processes including cell proliferation, which can be studied via 5-bromo-2'-deoxyuridine (BrdU) labeling. In a previous study we found that social stimulation altered both cell proliferation and levels of circulating gonadal steroids, leaving the issue of cause/effect unclear. In this study, we sought to determine whether socially modulated BrdU-labeling depends on gonadal hormone changes. We investigated this using a gonadectomy-implant paradigm and by exposing male and female green treefrogs (Hyla cinerea) to their conspecific chorus or control stimuli (i.e. random tones). Our results indicate that socially modulated cell proliferation occurred independently of gonadal hormone levels; furthermore, neither androgens in males nor estrogen in females increased cell proliferation in the preoptic area (POA) and infundibular hypothalamus, brain regions involved in endocrine regulation and acoustic communication. In fact, elevated estrogen levels decreased cell proliferation in those brain regions in the implanted female. In male frogs, evoked calling behavior was positively correlated with BrdU-labeling in the POA; however, statistical analysis showed that this behavior did not mediate socially induced cell proliferation. These results show that the social modulation of cell proliferation can occur without gonadal hormone involvement in either male or female adult anuran amphibians, and confirms that it is independent of a behavioral response in males.

Sex-specific modulation of cell proliferation by socially relevant stimuli in the adult green treefrog brain (Hyla cinerea)

Social experience plays an important role in regulating the neural, physiological and hormonal changes that accompany the expression of reproductive behavior in vertebrates. This suite of functions is sexually dimorphic, with different neural control areas preeminent in males and females. In anuran amphibians, social experience comes in the form of acoustic communication, which is central to their reproductive behavior. We sought to determine whether acoustic cues regulate cell proliferation in the brain of adult green treefrogs (Hyla cinerea). Our results show that both male and female treefrogs that heard a conspecific chorus during the breeding season exhibited increased brain cell proliferation compared to animals that heard random tones. Increased cell proliferation, as assessed by the number of 5-bromo-2'-deoxyuridine-immunoreactive (BrdU+) cells, were found near the ventricles of acoustically sensitive brain regions such as the preoptic area (POA) and the infundibular hypothalamus (IF). Sex differences emerged in the location of this socially modulated cell proliferation: increases occurred primarily in the male POA and the female IF. In addition, gonadal steroid hormones might have played a role in the social modulation of cell proliferation: by statistically control- ling for hormone level, we revealed that androgens might influence socially induced increases in BrdU+ cells in the male POA, but estrogen did not contribute to socially induced increases in the female IF. These results indicate that the reception of social cues increases cell proliferation in brain regions mediating sexual behavior and endocrine regulation, and moreover that social modulation of cell proliferation occurs in a sexually differentiated fashion.

Regional distribution and migration of proliferating cell populations in the adult brain of Hyla cinerea (Anura, Amphibia)

We examined the distribution of adult cell proliferation throughout the brain of an anuran amphibian using 5-bromo-2'-deoxyuridine (BrdU). BrdU, a thymidine analog, is a commonly used cellular marker that is incorporated into actively dividing progenitor cells. Adult green treefrogs, Hyla cinerea, received injections of BrdU and were sacrificed 2 h, 2 days, 2 weeks, or 30 days later. Immunohistochemistry revealed BrdU-immunopositive (BrdU+) cells to be distributed in ventricular zones throughout the brain. The heaviest concentrations of cells were located in the telencephalon, primarily in the ventrolateral region of the lateral ventricles, and the ventricles of olfactory bulbs. Numerous BrdU+ cells were located around the preoptic and hypothalamic recesses and few around the third ventricle in the diencephalon. Proceeding caudally towards the midbrain, there was a marked decrease in BrdU labeling and few BrdU+ cells were found in the hindbrain. Consistent with previous studies in ectothermic vertebrates, BrdU+ cells were found predominantly in the ventricular zone (VZ) and immediately adjacent to the VZ; at later time points (i.e., 30 days), the cells appeared to have migrated into parenchymal regions. The extent of cellular proliferation in anurans is similar to that of fishes and reptiles and thus is more widespread compared to mammals.

Environmental Enrichment Alters the Behavioral Profile of Ratsnakes (Elaphe)

This study investigated the effects of environmentally enriched and standard laboratory housing conditions on behavioral performance in 16 subadult ratsnakes (Elaphe obsoleta) using a split-clutch design. In a problem-solving task, snakes housed in enriched environments (EC) exhibited shorter latencies to the goal hole as compared to snakes housed in standard conditions (SC). In an open field task, EC snakes tended to habituate more quickly than SC snakes with repeated testing. A feeding task did not reveal any significant differences between EC and SC snakes. A discriminant function analysis correctly assigned all snakes to their appropriate housing treatment groups, based on the responses in each of the behavioral tasks. This suggests that each group had a distinct behavioral profile; that is, EC snakes were more behaviorally adaptive than SC snakes. This study demonstrated that housing conditions can affect the behavior of captive snakes and produce improvements in behavior similar to those seen in mammalian enrichment studies.

Multiple pathways of neuroprotection against oxidative stress and excitotoxic injury in immature primary hippocampal neurons

In the immature brain hydrogen peroxide accumulates after excitotoxic hypoxia-ischemia and is neurotoxic. Immature hippocampal neurons were exposed to N-methyl-D-aspartate (NMDA), a glutamate agonist, and hydrogen peroxide (H(2)O(2)) and the effects of free radical scavenging and transition metal chelation on neurotoxicity were studied. alpha-Phenyl-N-tert.-butylnitrone (PBN), a known superoxide scavenger, attenuated both H(2)O(2) and NMDA mediated toxicity. Treatment with desferrioxamine (DFX), an iron chelator, at the time of exposure to H(2)O(2) was ineffective, but pretreatment was protective. DFX also protected against NMDA toxicity. TPEN, a metal chelator with higher affinities for a broad spectrum of transition metal ions, also protected against H(2)O(2) toxicity but was ineffective against NMDA induced toxicity. These data suggest that during exposure to free radical and glutamate agonists, the presence of iron and other free metal ions contribute to neuronal cell death. In the immature nervous system this neuronal injury can be attenuated by free radical scavengers and metal chelators.