Multivariate investigation of aging in mouse models expressing the Alzheimer's protective APOE2 allele: integrating cognitive metrics, brain imaging, and blood transcriptomics

APOE allelic variation is critical in brain aging and Alzheimer's disease (AD). The APOE2 allele associated with cognitive resilience and neuroprotection against AD remains understudied. We employed a multipronged approach to characterize the transition from middle to old age in mice with APOE2 allele, using behavioral assessments, image-derived morphometry and diffusion metrics, structural connectomics, and blood transcriptomics. We used sparse multiple canonical correlation analyses (SMCCA) for integrative modeling, and graph neural network predictions. Our results revealed brain sub-networks associated with biological traits, cognitive markers, and gene expression. The cingulate cortex emerged as a critical region, demonstrating age-associated atrophy and diffusion changes, with higher fractional anisotropy in males and middle-aged subjects. Somatosensory and olfactory regions were consistently highlighted, indicating age-related atrophy and sex differences. The hippocampus exhibited significant volumetric changes with age, with differences between males and females in CA3 and CA1 regions. SMCCA underscored changes in the cingulate cortex, somatosensory cortex, olfactory regions, and hippocampus in relation to cognition and blood-based gene expression. Our integrative modeling in aging APOE2 carriers revealed a central role for changes in gene pathways involved in localization and the negative regulation of cellular processes. Our results support an important role of the immune system and response to stress. This integrative approach offers novel insights into the complex interplay among brain connectivity, aging, and sex. Our study provides a foundation for understanding the impact of APOE2 allele on brain aging, the potential for detecting associated changes in blood markers, and revealing novel therapeutic intervention targets.

Genomic Structural Equation Modeling Reveals Latent Phenotypes in the Human Cortex with Distinct Genetic Architecture

Genetic contributions to human cortical structure manifest pervasive pleiotropy. This pleiotropy may be harnessed to identify unique genetically-informed parcellations of the cortex that are neurobiologically distinct from functional, cytoarchitectural, or other cortical parcellation schemes. We investigated genetic pleiotropy by applying genomic structural equation modeling (SEM) to map the genetic architecture of cortical surface area (SA) and cortical thickness (CT) for the 34 brain regions recently reported in the ENIGMA cortical GWAS. Genomic SEM uses the empirical genetic covariance estimated from GWAS summary statistics with LD score regression (LDSC) to discover factors underlying genetic covariance, which we are denoting genetically informed brain networks (GIBNs). Genomic SEM can fit a multivariate GWAS from summary statistics for each of the GIBNs, which can subsequently be used for LD score regression (LDSC). We found the best-fitting model of cortical SA identified 6 GIBNs and CT identified 4 GIBNs. The multivariate GWASs of these GIBNs identified 74 genome-wide significant (GWS) loci (p<5×10-8), including many previously implicated in neuroimaging phenotypes, behavioral traits, and psychiatric conditions. LDSC of GIBN GWASs found that SA-derived GIBNs had a positive genetic correlation with bipolar disorder (BPD), and cannabis use disorder, indicating genetic predisposition to a larger SA in the specific GIBN is associated with greater genetic risk of these disorders. A negative genetic correlation was observed with attention deficit hyperactivity disorder (ADHD), major depressive disorder (MDD), and insomnia, indicating genetic predisposition to a larger SA in the specific GIBN is associated with lower genetic risk of these disorders. CT GIBNs displayed a negative genetic correlation with alcohol dependence. Jointly modeling the genetic architecture of complex traits and investigating multivariate genetic links across phenotypes offers a new vantage point for mapping the cortex into genetically informed networks.

Polarized Desmosome and Hemidesmosome Shedding via Exosomes is an Early Indicator of Outer Blood-Retina Barrier Dysfunction

The retinal pigmented epithelium (RPE) constitutes the outer blood-retinal barrier, enables photoreceptor function of the eye, and is constantly exposed to oxidative stress. As such, dysfunction of the RPE underlies pathology leading to development of age-related macular degeneration (AMD), the leading cause of vision loss among the elderly in industrialized nations. A major responsibility of the RPE is to process photoreceptor outer segments, which relies on the proper functioning of its endocytic pathways and endosomal trafficking. Exosomes and other extracellular vesicles from RPE are an essential part of these pathways and may be early indicators of cellular stress. To test the role of exosomes that may underlie the early stages of AMD, we used a polarized primary RPE cell culture model under chronic subtoxic oxidative stress. Unbiased proteomic analyses of highly purified basolateral exosomes from oxidatively stressed RPE cultures revealed changes in proteins involved in epithelial barrier integrity. There were also significant changes in proteins accumulating in the basal-side sub-RPE extracellular matrix during oxidative stress, that could be prevented with an inhibitor of exosome release. Thus, chronic subtoxic oxidative stress in primary RPE cultures induces changes in exosome content, including basal-side specific desmosome and hemidesmosome shedding via exosomes. These findings provide novel biomarkers of early cellular dysfunction and opportunity for therapeutic intervention in age-related retinal diseases, (e.g., AMD) and broadly from blood-CNS barriers in other neurodegenerative diseases.

Genetic regulation of fetal hemoglobin across global populations

Human genetic variation has enabled the identification of several key regulators of fetal-to-adult hemoglobin switching, including BCL11A, resulting in therapeutic advances. However, despite the progress made, limited further insights have been obtained to provide a fuller accounting of how genetic variation contributes to the global mechanisms of fetal hemoglobin (HbF) gene regulation. Here, we have conducted a multi-ancestry genome-wide association study of 28,279 individuals from several cohorts spanning 5 continents to define the architecture of human genetic variation impacting HbF. We have identified a total of 178 conditionally independent genome-wide significant or suggestive variants across 14 genomic windows. Importantly, these new data enable us to better define the mechanisms by which HbF switching occurs in vivo. We conduct targeted perturbations to define BACH2 as a new genetically-nominated regulator of hemoglobin switching. We define putative causal variants and underlying mechanisms at the well-studied BCL11A and HBS1L-MYB loci, illuminating the complex variant-driven regulation present at these loci. We additionally show how rare large-effect deletions in the HBB locus can interact with polygenic variation to influence HbF levels. Our study paves the way for the next generation of therapies to more effectively induce HbF in sickle cell disease and β-thalassemia.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Evaluating Associations between Average Pain Intensity and Genetic Variation in People with Sickle Cell Disease: An Exploratory Study

BACKGROUND

Pain is one of the most common and deleterious symptoms experienced by individuals with sickle cell disease (SCD). There is a paucity of studies identifying potential genetic mechanisms of pain in this population.

AIM

Examine associations between 11 functional single nucleotide polymorphisms in 9 candidate genes with reports of average pain intensity in individuals with sickle cell disease.

METHOD

Cross-sectional analyses were performed on data and blood samples collected through the Duke SCD Implementation Consortium Registry. Participants were asked to rate their pain "on the average" using an 11-point numeric rating scale (0 = no pain; 10 = pain as bad as you can imagine). We genotyped 11 single nucleotide polymorphisms in 9 pain-related genes using TaqMan® Genotyping Assays. Associations between each polymorphism and reports of average pain were evaluated.

RESULTS

The 86 participants (mean age: 28.7 years; 64% female) included in this study reported moderate pain on average (Mean = 4, Standard Deviation = 2.4). ICAM1 rs1799969 was the only genetic polymorphism that was significantly associated with pain (p = .01). Individuals with one or more minor alleles had lower average pain (Mean = 1.25, Standard Deviation = 1.50) than individuals without a minor allele (Mean = 4.13, Standard Deviation = 2.25). The effect size for ICAM1 rs1799969 was 1.30, which is considered large. The effect sizes for all other single nucleotide polymorphisms ranged from small to medium (range: 0-0.3).

CONCLUSIONS

Our findings provide preliminary evidence that the minor allele in ICAM1 rs1799969 had protective effects against experiencing more severe pain in sickle cell disease.

Structural variation across 138,134 samples in the TOPMed consortium

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Whole genome sequencing identifies candidate genes for familial essential tremor and reveals biological pathways implicated in essential tremor aetiology

Background

Essential tremor (ET), one of the most common neurological disorders, has a phenotypically heterogeneous presentation characterized by bilateral kinetic tremor of the arms and, in some patients, tremor involving other body regions (e.g., head, voice). Genetic studies suggest that ET is genetically heterogeneous.

Methods

We analyzed whole genome sequence data (WGS) generated on 104 multi-generational white families with European ancestry affected by ET. Genome-wide parametric linkage and association scans were analyzed using adjusted logistic regression models through the application of the Pseudomarker software. To investigate the additional contribution of rare variants in familial ET, we also performed an aggregate variant non-parametric linkage (NPL) analysis using the collapsed haplotype method implemented in CHP-NPL software.

Findings

Parametric linkage analysis of common variants identified several loci with significant evidence of linkage (HLOD ≥3.6). Among the gene regions within the strongest ET linkage peaks were BTC (4q13.3, HLOD=4.53), N6AMT1 (21q21.3, HLOD=4.31), PCDH9 (13q21.32, HLOD=4.21), EYA1 (8q13.3, HLOD=4.04), RBFOX1 (16p13.3, HLOD=4.02), MAPT (17q21.31, HLOD=3.99) and SCARB2 (4q21.1, HLOD=3.65). CHP-NPL analysis identified fifteen additional genes with evidence of significant linkage (LOD ≥3.8). These genes include TUBB2A, VPS33B, STEAP1B, SPINK5, ZRANB1, TBC1D3C, PDPR, NPY4R, ETS2, ZNF736, SPATA21, ARL17A, PZP, BLK and CCDC94. In one ET family contributing to the linkage peak on chromosome 16p13.3, we identified a likely pathogenic heterozygous canonical splice acceptor variant in exon 2 of RBFOX1 (ENST00000547372; c.4-2A>G), that co-segregated with the ET phenotype in the family.

Interpretation

Linkage and association analyses of WGS identified several novel ET candidate genes, which are implicated in four major pathways that include 1) the epidermal growth factor receptor-phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha-AKT serine/threonine kinase 1 (EGFR-PI3K-AKT) and Mitogen-activated protein Kinase 1 (ERK) pathways, 2) Reactive oxygen species (ROS) and DNA repair, 3) gamma-aminobutyric acid-ergic (GABAergic) system and 4) RNA binding and regulation of RNA processes. Our study provides evidence for a possible overlap in the genetic architecture of ET, neurological disease, cancer and aging. The genes and pathways identified can be prioritized in future genetic and functional studies.

Funding

National Institutes of Health, NINDS, NS073872 (USA) and NIA AG058131(USA).

Alcohol use and alcohol use disorder differ in their genetic relationships with PTSD: A genomic structural equation modelling approach

PURPOSE

Posttraumatic Stress Disorder (PTSD) is associated with increased alcohol use and alcohol use disorder (AUD), which are all moderately heritable. Studies suggest the genetic association between PTSD and alcohol use differs from that of PTSD and AUD, but further analysis is needed.

BASIC PROCEDURES

We used genomic Structural Equation Modeling (genomicSEM) to analyze summary statistics from large-scale genome-wide association studies (GWAS) of European Ancestry participants to investigate the genetic relationships between PTSD (both diagnosis and re-experiencing symptom severity) and a range of alcohol use and AUD phenotypes.

MAIN FINDINGS

When we differentiated genetic factors for alcohol use and AUD we observed improved model fit relative to models with all alcohol-related indicators loading onto a single factor. The genetic correlations (rG) of PTSD were quite discrepant for the alcohol use and AUD factors. This was true when modeled as a three-correlated-factor model (PTSD-AUD rG:.36, p < .001; PTSD-alcohol use rG: -0.17, p < .001) and as a Bifactor model, in which the common and unique portions of alcohol phenotypes were pulled out into an AUD-specific factor (rG with PTSD:.40, p < .001), AU-specific factor (rG with PTSD: -0.57, p < .001), and a common alcohol factor (rG with PTSD:.16, NS).

PRINCIPAL CONCLUSIONS

These results indicate the genetic architecture of alcohol use and AUD are differentially associated with PTSD. When the portions of variance unique to alcohol use and AUD are extracted, their genetic associations with PTSD vary substantially, suggesting different genetic architectures of alcohol phenotypes in people with PTSD.

Clonal hematopoiesis in sickle cell disease

BACKGROUND

Curative gene therapies for sickle cell disease (SCD) are currently undergoing clinical evaluation. The occurrence of myeloid malignancies in these trials has prompted safety concerns. Individuals with SCD are predisposed to myeloid malignancies, but the underlying causes remain undefined. Clonal hematopoiesis (CH) is a pre-malignant condition that also confers significant predisposition to myeloid cancers. While it has been speculated that CH may play a role in SCD-associated cancer predisposition, limited data addressing this issue have been reported.

METHODS

Here, we leveraged 74,190 whole genome sequences to robustly study CH in SCD. Somatic mutation calling methods were used to assess CH in all samples and comparisons between individuals with and without SCD were performed.

RESULTS

While we had sufficient power to detect a greater than 2-fold increased rate of CH, we found no detectable variation in rate or clone properties between individuals affected by SCD and controls. The rate of CH in individuals with SCD was unaltered by hydroxyurea use.

CONCLUSIONS

We did not observe an increased risk for acquiring detectable CH in SCD, at least as measured by whole genome sequencing. These results should help guide ongoing efforts and further studies that seek to better define the risk factors underlying myeloid malignancy predisposition in SCD and help ensure that curative therapies can be more safely applied.

FUNDING

Funding was provided by the New York Stem Cell Foundation and National Institutes of Health. The funders had no role in study design or reporting.

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation

Common genetic risk for neuropsychiatric disorders is enriched in regulatory elements active during cortical neurogenesis. However, it remains poorly understood as to how these variants influence gene regulation. To model the functional impact of common genetic variation on the noncoding genome during human cortical development, we performed the assay for transposase accessible chromatin using sequencing (ATAC-seq) and analyzed chromatin accessibility quantitative trait loci (QTL) in cultured human neural progenitor cells and their differentiated neuronal progeny from 87 donors. We identified significant genetic effects on 988/1,839 neuron/progenitor regulatory elements, with highly cell-type and temporally specific effects. A subset (roughly 30%) of chromatin accessibility-QTL were also associated with changes in gene expression. Motif-disrupting alleles of transcriptional activators generally led to decreases in chromatin accessibility, whereas motif-disrupting alleles of repressors led to increases in chromatin accessibility. By integrating cell-type-specific chromatin accessibility-QTL and brain-relevant genome-wide association data, we were able to fine-map and identify regulatory mechanisms underlying noncoding neuropsychiatric disorder risk loci.

Sex and Menopause Modify the Effect of Single Nucleotide Polymorphism Genotypes on Fibrosis in NAFLD

The development of fibrosis in nonalcoholic fatty liver disease (NAFLD) is influenced by genetics, sex, and menopausal status, but whether genetic susceptibility to fibrosis is influenced by sex and reproductive status is unclear. Our aim was to identify metabolism-related single nucleotide polymorphisms (SNPs), whose effect on NAFLD fibrosis is significantly modified by sex and menopausal status. We performed a cross-sectional, proof-of-concept study of 616 patients in the Duke NAFLD Clinical Database and Biorepository. The primary outcome was nonalcoholic steatohepatitis-Clinical Research Network (NASH-CRN) fibrosis stage. Menopause status was self-reported; age 51 years was used as a surrogate for menopause in patients with missing menopause data. The Metabochip was used to obtain 98,359 SNP genotypes in known metabolic pathway genes for each patient. We used additive genetic models to characterize sex and menopause-specific effects of SNP genotypes on NAFLD fibrosis stage. In the main effects analysis, none of the SNPs were associated with fibrosis at P < 0.05 after correcting for multiple comparisons. Twenty-five SNPs significantly interacted with sex/menopause to affect fibrosis stage (interaction P < 0.0001). After removal of loci in linkage disequilibrium, 10 independent loci were identified. Six were in the following genes: KCNIP4 (potassium voltage-gated channel interacting protein 4), PSORS1C1 (psoriasis susceptibility 1 candidate 1), KLHL8 (Kelch-like family member 8), GLRA1 (glycine receptor alpha 1), NOTCH2 (notch receptor 2), and PRKCH (protein kinase C eta), and four SNPs were intergenic. In stratified models, four SNPs were significant in premenopausal and postmenopausal women, three only in postmenopausal women, two in men and postmenopausal women, and one only in premenopausal women. Conclusion: We identified 10 loci with a significant sex/menopause interaction with respect to fibrosis. None of these SNPs were significant in all sex/menopause groups, suggesting modulation of genetic susceptibility to fibrosis by sex and menopause status. Future studies of genetic predictors of NAFLD progression should account for sex and menopause.

Human myelomeningocele risk and ultra-rare deleterious variants in genes associated with cilium, WNT-signaling, ECM, cytoskeleton and cell migration

Myelomeningocele (MMC) affects one in 1000 newborns annually worldwide and each surviving child faces tremendous lifetime medical and caregiving burdens. Both genetic and environmental factors contribute to disease risk but the mechanism is unclear. This study examined 506 MMC subjects for ultra-rare deleterious variants (URDVs, absent in gnomAD v2.1.1 controls that have Combined Annotation Dependent Depletion score ≥ 20) in candidate genes either known to cause abnormal neural tube closure in animals or previously associated with human MMC in the current study cohort. Approximately 70% of the study subjects carried one to nine URDVs among 302 candidate genes. Half of the study subjects carried heterozygous URDVs in multiple genes involved in the structure and/or function of cilium, cytoskeleton, extracellular matrix, WNT signaling, and/or cell migration. Another 20% of the study subjects carried heterozygous URDVs in candidate genes associated with gene transcription regulation, folate metabolism, or glucose metabolism. Presence of URDVs in the candidate genes involving these biological function groups may elevate the risk of developing myelomeningocele in the study cohort.

Diversity of variant alleles encoding Kidd, Duffy, and Kell antigens in individuals with sickle cell disease using whole genome sequencing data from the NHLBI TOPMed Program

BACKGROUND

Genetic variants in the SLC14A1, ACKR1, and KEL genes, which encode Kidd, Duffy, and Kell red blood cell antigens, respectively, may result in weakened expression of antigens or a null phenotype. These variants are of particular interest to individuals with sickle cell disease (SCD), who frequently undergo chronic transfusion therapy with antigen-matched units. The goal was to describe the diversity and the frequency of variants in SLC14A1, ACKR1, and KEL genes among individuals with SCD using whole genome sequencing (WGS) data.

STUDY DESIGN AND METHODS

Two large SCD cohorts were studied: the Recipient Epidemiology and Donor Evaluation Study III (REDS-III) (n = 2634) and the Outcome Modifying Gene in SCD (OMG) (n = 640). Most of the studied individuals were of mixed origin. WGS was performed as part of the National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine (TOPMed) program.

RESULTS

In SLC14A1, variants included four encoding a weak Jka phenotype and five null alleles (JKnull ). JKA*01N.09 was the most common JKnull . One possible JKnull mutation was novel: c.812G>T. In ACKR1, identified variants included two that predicted Fyx (FY*X) and one corresponding to the c.-67T>C GATA mutation. The c.-67T>C mutation was associated with FY*A (FY*01N.01) in four participants. FY*X was identified in 49 individuals. In KEL, identified variants included three null alleles (KEL*02N.17, KEL*02N.26, and KEL*02N.04) and one allele predicting Kmod phenotype, all in heterozygosity.

CONCLUSIONS

We described the diversity and distribution of SLC14A1, ACKR1, and KEL variants in two large SCD cohorts, comprising mostly individuals of mixed ancestry. This information may be useful for planning the transfusion support of patients with SCD.

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

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

Burden of rare deleterious variants in WNT signaling genes among 511 myelomeningocele patients

Genes in the noncanonical WNT signaling pathway controlling planar cell polarity have been linked to the neural tube defect myelomeningocele. We hypothesized that some genes in the WNT signaling network have a higher mutational burden in myelomeningocele subjects than in reference subjects in gnomAD. Exome sequencing data from 511 myelomeningocele subjects was obtained in-house and data from 29,940 ethnically matched subjects was provided by version 2 of the publicly available Genome Aggregation Database. To compare mutational burden, we collapsed rare deleterious variants across each of 523 human WNT signaling genes in case and reference populations. Ten WNT signaling genes were disrupted with a higher mutational burden among Mexican American myelomeningocele subjects compared to reference subjects (Fishers exact test, P ≤ 0.05) and seven different genes were disrupted among individuals of European ancestry compared to reference subjects. Gene ontology enrichment analyses indicate that genes disrupted only in the Mexican American population play a role in planar cell polarity whereas genes identified in both populations are important for the regulation of canonical WNT signaling. In summary, evidence for WNT signaling genes that may contribute to myelomeningocele in humans is presented and discussed.

Identification of novel candidate risk genes for myelomeningocele within the glucose homeostasis/oxidative stress and folate/one-carbon metabolism networks

Background

Neural tube defects (NTDs) are the second most common complex birth defect, yet, our understanding of the genetic contribution to their development remains incomplete. Two environmental factors associated with NTDs are Folate and One Carbon Metabolism (FOCM) and Glucose Homeostasis and Oxidative Stress (GHOS). Utilizing next‐generation sequencing of a large patient cohort, we identify novel candidate genes in these two networks to provide insights into NTD mechanisms.

Methods

Exome sequencing (ES) was performed in 511 patients, born with myelomeningocele, divided between European American and Mexican American ethnicities. Healthy control data from the Genome Aggregation database were ethnically matched and used as controls. Rare, high fidelity, nonsynonymous predicted damaging missense, nonsense, or canonical splice site variants in independently generated candidate gene lists for FOCM and GHOS were identified. We used a gene‐based collapsing approach to quantify mutational burden in case and controls, with the control cohort estimated using cumulative allele frequencies assuming Hardy–Weinberg equilibrium.

Results

We identified 45 of 837 genes in the FOCM network and 22 of 568 genes in the GHOS network as possible NTD risk genes with p < 0.05. No nominally significant risk genes were shared between ethnicities. Using a novel approach to mutational burden we identify 55 novel NTD risk associations.

Conclusions

We provide a means of utilizing large publicly available sequencing datasets as controls for sequencing projects examining rare disease. This approach confirmed existing risk genes for myelomeningocele and identified possible novel risk genes. Lastly, it suggests possible distinct genetic etiologies for this malformation between different ethnicities.

The genetic architecture of the human cerebral cortex

The cerebral cortex underlies our complex cognitive capabilities, yet little is known about the specific genetic loci that influence human cortical structure. To identify genetic variants that affect cortical structure, we conducted a genome-wide association meta-analysis of brain magnetic resonance imaging data from 51,665 individuals. We analyzed the surface area and average thickness of the whole cortex and 34 regions with known functional specializations. We identified 199 significant loci and found significant enrichment for loci influencing total surface area within regulatory elements that are active during prenatal cortical development, supporting the radial unit hypothesis. Loci that affect regional surface area cluster near genes in Wnt signaling pathways, which influence progenitor expansion and areal identity. Variation in cortical structure is genetically correlated with cognitive function, Parkinson's disease, insomnia, depression, neuroticism, and attention deficit hyperactivity disorder.

The TFAP2A-IRF6-GRHL3 genetic pathway is conserved in neurulation

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect.

Development of Common Data Elements for Use in Chiari Malformation Type I Clinical Research: An NIH/NINDS Project

The management of Chiari I malformation (CMI) is controversial because treatment methods vary and treatment decisions rest on incomplete understanding of its complex symptom patterns, etiologies, and natural history. Validity of studies that attempt to compare treatment of CMI has been limited because of variable terminology and methods used to describe study subjects. The goal of this project was to standardize terminology and methods by developing a comprehensive set of Common Data Elements (CDEs), data definitions, case report forms (CRFs), and outcome measure recommendations for use in CMI clinical research, as part of the CDE project at the National Institute of Neurological Disorders and Stroke (NINDS) of the US National Institutes of Health. A working group, comprising over 30 experts, developed and identified CDEs, template CRFs, data dictionaries, and guidelines to aid investigators starting and conducting CMI clinical research studies. The recommendations were compiled, internally reviewed, and posted online for external public comment. In October 2016, version 1.0 of the CMI CDE recommendations became available on the NINDS CDE website. The recommendations span these domains: Core Demographics/Epidemiology; Presentation/Symptoms; Co-Morbidities/Genetics; Imaging; Treatment; and Outcome. Widespread use of CDEs could facilitate CMI clinical research trial design, data sharing, retrospective analyses, and consistent data sharing between CMI investigators around the world. Updating of CDEs will be necessary to keep them relevant and applicable to evolving research goals for understanding CMI and its treatment.

Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia

Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.

Testosterone Pathway Genetic Polymorphisms in Relation to Primary Open-Angle Glaucoma: An Analysis in Two Large Datasets

Purpose

Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk.

Methods

We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]).

Results

In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations.

Conclusions

Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

APOE ε4 associated with preserved executive function performance and maintenance of temporal and cingulate brain volumes in younger adults

The APOE ε4 allele is associated with cognitive deficits and brain atrophy in older adults, but studies in younger adults are mixed. We examined APOE genotype effects on cognition and brain structure in younger adults and whether genotype effects differed by age and with presence of depression. 157 adults (32 % ε4 carriers, 46 % depressed) between 20 and 50 years of age completed neuropsychological testing, 131 of which also completed 3 T cranial MRI. We did not observe a direct effect of APOE genotype on cognitive performance or structural MRI measures. A significant genotype by age interaction was observed for executive function, where age had less of an effect on executive function in ε4 carriers. Similar interactions were observed for the entorhinal cortex, rostral and caudal anterior cingulate cortex and parahippocampal gyrus, where the effect of age on regional volumes was reduced in ε4 carriers. There were no significant interactions between APOE genotype and depression diagnosis. The ε4 allele benefits younger adults by allowing them to maintain executive function performance and volumes of cingulate and temporal cortex regions with aging, at least through age fifty years.

Genetic correlations between intraocular pressure, blood pressure and primary open-angle glaucoma: a multi-cohort analysis

Primary open-angle glaucoma (POAG) is the most common chronic optic neuropathy worldwide. Epidemiological studies show a robust positive relation between intraocular pressure (IOP) and POAG and modest positive association between IOP and blood pressure (BP), while the relation between BP and POAG is controversial. The International Glaucoma Genetics Consortium (n=27 558), the International Consortium on Blood Pressure (n=69 395), and the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (n=37 333), represent genome-wide data sets for IOP, BP traits and POAG, respectively. We formed genome-wide significant variant panels for IOP and diastolic BP and found a strong relation with POAG (odds ratio and 95% confidence interval: 1.18 (1.14-1.21), P=1.8 × 10^-27) for the former trait but no association for the latter (P=0.93). Next, we used linkage disequilibrium (LD) score regression, to provide genome-wide estimates of correlation between traits without the need for additional phenotyping. We also compared our genome-wide estimate of heritability between IOP and BP to an estimate based solely on direct measures of these traits in the Erasmus Rucphen Family (ERF; n=2519) study using Sequential Oligogenic Linkage Analysis Routines (SOLAR). LD score regression revealed high genetic correlation between IOP and POAG (48.5%, P=2.1 × 10^-5); however, genetic correlation between IOP and diastolic BP (P=0.86) and between diastolic BP and POAG (P=0.42) were negligible. Using SOLAR in the ERF study, we confirmed the minimal heritability between IOP and diastolic BP (P=0.63). Overall, IOP shares genetic basis with POAG, whereas BP has limited shared genetic correlation with IOP or POAG.

Age at natural menopause genetic risk score in relation to age at natural menopause and primary open-angle glaucoma in a US-based sample

OBJECTIVE

Several attributes of female reproductive history, including age at natural menopause (ANM), have been related to primary open-angle glaucoma (POAG). We assembled 18 previously reported common genetic variants that predict ANM to determine their association with ANM or POAG.

METHODS

Using data from the Nurses' Health Study (7,143 women), we validated the ANM weighted genetic risk score in relation to self-reported ANM. Subsequently, to assess the relation with POAG, we used data from 2,160 female POAG cases and 29,110 controls in the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (NEIGHBORHOOD), which consists of 8 datasets with imputed genotypes to 5.6+ million markers. Associations with POAG were assessed in each dataset, and site-specific results were meta-analyzed using the inverse weighted variance method.

RESULTS

The genetic risk score was associated with self-reported ANM (P = 2.2 × 10) and predicted 4.8% of the variance in ANM. The ANM genetic risk score was not associated with POAG (Odds Ratio (OR) = 1.002; 95% Confidence Interval (CI): 0.998, 1.007; P = 0.28). No single genetic variant in the panel achieved nominal association with POAG (P ≥0.20). Compared to the middle 80 percent, there was also no association with the lowest 10 percentile or highest 90 percentile of genetic risk score with POAG (OR = 0.75; 95% CI: 0.47, 1.21; P = 0.23 and OR = 1.10; 95% CI: 0.72, 1.69; P = 0.65, respectively).

CONCLUSIONS

A genetic risk score predicting 4.8% of ANM variation was not related to POAG; thus, genetic determinants of ANM are unlikely to explain the previously reported association between the two phenotypes.

Transcriptome analysis of adult and fetal trabecular meshwork, cornea, and ciliary body tissues by RNA sequencing

PURPOSE

To characterize the transcriptional landscape of human adult and fetal trabecular meshwork (TM), cornea, and ciliary body (CB) tissues, and to evaluate the expression level of candidate genes selected from genetic association studies of primary-open angle glaucoma, central corneal thickness, intraocular pressure, vertical cup to disc ratio, and optic nerve parameters.

METHODS

Deep RNA sequencing was performed on the selected human tissues. Transcriptome analyses were performed to 1) characterize the total number of expressed genes, 2) identify the most highly expressed genes, 3) estimate the number of novel transcripts, and 4) evaluate the expression of candidate genes in each tissue. Finally, a differential gene expression analysis was conducted to compare the adult and fetal ocular tissues.

RESULTS

There was an average of 12,362 protein coding genes and 3725 novel transcripts expressed in each tissue. The top most expressed genes in each tissue included SPARC (fetal cornea and TM), APOD (adult TM), CLU (adult cornea), and PTGDS (adult and fetal CB). Twenty-nine candidate genes selected from genetic association studies primarily showed high expression levels in the trabecular meshwork and cornea. Comparison of adult and fetal samples identified 2012 and 1261 differentially expressed protein-coding genes within the cornea and trabecular meshwork, respectively.

CONCLUSIONS

This study has provided an unbiased glimpse into the transcriptome of three essential anterior ocular tissues, resulting in the development of several novel hypotheses. These data can be used in the future to better guide ocular research questions.

Genome-wide association analysis identifies TXNRD2, ATXN2 and FOXC1 as susceptibility loci for primary open-angle glaucoma

Primary open angle glaucoma (POAG) is a leading cause of blindness world-wide. To identify new susceptibility loci, we meta-analyzed GWAS results from 8 independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significant SNPs in two Australian studies (1,252 cases and 2,592 controls), 3 European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of top SNPs identified three novel loci: rs35934224[T] within TXNRD2 (odds ratio (OR) = 0.78, P = 4.05×10⁻¹¹ encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] within ATXN2 (OR = 1.17, P = 8.73×10⁻¹⁰), and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76×10⁻¹⁰). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest novel targets for preventative therapies.

TMEM231, mutated in orofaciodigital and Meckel syndromes, organizes the ciliary transition zone

TMEM231, a functional component of the MKS complex at the ciliary transition zone, is mutated in orofaciodigital syndrome type 3 and Meckel syndrome.

The Meckel syndrome (MKS) complex functions at the transition zone, located between the basal body and axoneme, to regulate the localization of ciliary membrane proteins. We investigated the role of Tmem231, a two-pass transmembrane protein, in MKS complex formation and function. Consistent with a role in transition zone function, mutation of mouse Tmem231 disrupts the localization of proteins including Arl13b and Inpp5e to cilia, resulting in phenotypes characteristic of MKS such as polydactyly and kidney cysts. Tmem231 and B9d1 are essential for each other and other complex components such as Mks1 to localize to the transition zone. As in mouse, the Caenorhabditis elegans orthologue of Tmem231 localizes to and controls transition zone formation and function, suggesting an evolutionarily conserved role for Tmem231. We identified TMEM231 mutations in orofaciodigital syndrome type 3 (OFD3) and MKS patients that compromise transition zone function. Thus, Tmem231 is critical for organizing the MKS complex and controlling ciliary composition, defects in which cause OFD3 and MKS.

Pregnancy continuation and organizational religious activity following prenatal diagnosis of a lethal fetal defect are associated with improved psychological outcome

OBJECTIVE

The aim of the article is to examine the psychological impact, specifically symptoms of grief, post-traumatic stress and depression, in women and men who either terminated or continued a pregnancy following prenatal diagnosis of a lethal fetal defect.

METHOD

This project investigated a diagnostically homogeneous group composed of 158 women and 109 men who lost a pregnancy to anencephaly, a lethal neural tube defect. Participants completed the Perinatal Grief Scale, Impact of Event Scale - Revised and Beck Depression Inventory-II, which measure symptoms of grief, post-traumatic stress and depression, respectively. Demographics, religiosity and pregnancy choices were also collected. Gender-specific analysis of variance was performed for instrument total scores and subscales.

RESULTS

Women who terminated reported significantly more despair (p = 0.02), avoidance (p = 0.008) and depression (p = 0.04) than women who continued the pregnancy. Organizational religious activity was associated with a reduction in grief (Perinatal Grief Scale subscales) in both women (p = 0.02, p = 0.04 and p = 0.03) and men (p = 0.047).

CONCLUSION

There appears to be a psychological benefit to women to continue the pregnancy following a lethal fetal diagnosis. Following a lethal fetal diagnosis, the risks and benefits, including psychological effects, of termination and continuation of pregnancy should be discussed in detail with an effort to be as nondirective as possible.

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

BACKGROUND

The apolipoprotein E (APOE) ε4 allele has been implicated in a range of neuropsychiatric conditions. The present research examined if the ε4 allele of the APOE gene moderated the effect of combat exposure on posttraumatic stress disorder (PTSD) among Iraq/Afghanistan-era veterans.

METHOD

Participants included 765 non-Hispanic White (NHW) and 859 non-Hispanic Black (NHB) Iraq/Afghanistan-era veterans. A structured interview established psychiatric diagnoses. Combat exposure and PTSD symptom severity were assessed via self-report.

RESULTS

The most common lifetime diagnoses were depression (39.2%), PTSD (38.4%), and alcohol dependence (24.38%). After correcting for multiple comparisons, no significant effects were observed on any of the outcomes among the NHW sample; however, within the NHB sample, significant gene × environment (G × E) interactions were observed for lifetime PTSD (P = .0029) and PTSD symptom severity (P = .0009). In each case, the APOE ε4 allele had no effect on the outcomes when combat exposure was low; however, when combat exposure was high, an additive effect was observed such that ε4 homozygotes exposed to high levels of combat reported the highest rates of PTSD (92%) and the worst symptom severity scores on the Davidson Trauma Scale (M = 79.5).

CONCLUSIONS

Although preliminary, these findings suggest that the APOE ε4 allele, in conjunction with exposure to high levels of combat exposure, may increase veterans' risk for developing PTSD.

Angiogenic, neurotrophic, and inflammatory system SNPs moderate the association between birth weight and ADHD symptom severity

Low birth weight is associated with increased risk for Attention-Deficit/Hyperactivity Disorder (ADHD); however, the etiological underpinnings of this relationship remain unclear. This study investigated if genetic variants in angiogenic, dopaminergic, neurotrophic, kynurenine, and cytokine-related biological pathways moderate the relationship between birth weight and ADHD symptom severity. A total of 398 youth from two multi-site, family-based studies of ADHD were included in the analysis. The sample consisted of 360 ADHD probands, 21 affected siblings, and 17 unaffected siblings. A set of 164 SNPs from 31 candidate genes, representing five biological pathways, were included in our analyses. Birth weight and gestational age data were collected from a state birth registry, medical records, and parent report. Generalized Estimating Equations tested for main effects and interactions between individual SNPs and birth weight centile in predicting ADHD symptom severity. SNPs within neurotrophic (NTRK3) and cytokine genes (CNTFR) were associated with ADHD inattentive symptom severity. There was no main effect of birth weight centile on ADHD symptom severity. SNPs within angiogenic (NRP1 & NRP2), neurotrophic (NTRK1 & NTRK3), cytokine (IL16 & S100B), and kynurenine (CCBL1 & CCBL2) genes moderate the association between birth weight centile and ADHD symptom severity. The SNP main effects and SNP × birth weight centile interactions remained significant after adjusting for multiple testing. Genetic variability in angiogenic, neurotrophic, and inflammatory systems may moderate the association between restricted prenatal growth, a proxy for an adverse prenatal environment, and risk to develop ADHD.

Association of gene variants of the renin-angiotensin system with accelerated hippocampal volume loss and cognitive decline in old age

OBJECTIVE

Genetic factors confer risk for neuropsychiatric phenotypes, but the polygenic etiology of these phenotypes makes identification of genetic culprits challenging. An approach to this challenge is to examine the effects of genetic variation on relevant endophenotypes, such as hippocampal volume loss. A smaller hippocampus is associated with gene variants of the renin-angiotensin system (RAS), a system implicated in vascular disease. However, no studies to date have investigated longitudinally the effects of genetic variation of RAS on the hippocampus.

METHOD

The authors examined the effects of polymorphisms of AGTR1, the gene encoding angiotensin-II type 1 receptor of RAS, on longitudinal hippocampal volumes of older adults. In all, 138 older adults (age ≥60 years) were followed for an average of about 4 years. The participants underwent repeated structural MRI and comprehensive neurocognitive testing, and they were genotyped for four AGTR1 single-nucleotide polymorphisms (SNPs) with low pairwise linkage disequilibrium values and apolipoprotein E (APOE) genotype.

RESULTS

Genetic variants at three AGTR1 SNPs (rs2638363, rs1492103, and rs2675511) were independently associated with accelerated hippocampal volume loss over the 4-year follow-up period in the right but not left hemisphere. Intriguingly, these AGTR1 risk alleles also predicted worse episodic memory performance but were not related to other cognitive measures. Two risk variants (rs2638363 and rs12721331) interacted with the APOE4 allele to accelerate right hippocampal volume loss.

CONCLUSIONS

Risk genetic variants of the RAS may accelerate memory decline in older adults, an effect that may be conferred by accelerated hippocampal volume loss. Molecules involved in this system may hold promise as early therapeutic targets for late-life neuropsychiatric disorders.

Heavy metals, organic solvents, and multiple sclerosis: An exploratory look at gene-environment interactions

Exposure to heavy metals and organic solvents are potential etiologic factors for multiple sclerosis (MS), but their interaction with MS-associated genes is under-studied. The authors explored the relationship between environmental exposure to lead, mercury, and solvents and 58 single-nucleotide polymorphisms (SNPs) in MS-associated genes. Data from a population-based case-control study of 217 prevalent MS cases and 496 age-, race-, gender-, and geographically matched controls were used to fit conditional logistic regression models of the association between the chemical, gene, and MS, adjusting for education and ancestry. MS cases were more likely than controls to report lead (odds ratio [OR] = 2.03; 95% confidence interval [CI]: 1.07, 3.86) and mercury exposure (OR = 2.06; 95% CI: 1.08, 3.91). Findings of potential gene-environment interactions between SNPs in TNF-α, TNF-β, TCA-β, VDR, MBP, and APOE, and lead, mercury, or solvents should be considered cautiously due to limited sample size.

Association between the oxytocin receptor (OXTR) gene and mesolimbic responses to rewards

BACKGROUND

There has been significant progress in identifying genes that confer risk for autism spectrum disorders (ASDs). However, the heterogeneity of symptom presentation in ASDs impedes the detection of ASD risk genes. One approach to understanding genetic influences on ASD symptom expression is to evaluate relations between variants of ASD candidate genes and neural endophenotypes in unaffected samples. Allelic variations in the oxytocin receptor (OXTR) gene confer small but significant risk for ASDs for which the underlying mechanisms may involve associations between variability in oxytocin signaling pathways and neural response to rewards. The purpose of this preliminary study was to investigate the influence of allelic variability in the OXTR gene on neural responses to monetary rewards in healthy adults using functional magnetic resonance imaging (fMRI).

METHODS

The moderating effects of three single nucleotide polymorphisms (SNPs) (rs1042778, rs2268493 and rs237887) of the OXTR gene on mesolimbic responses to rewards were evaluated using a monetary incentive delay fMRI task.

RESULTS

T homozygotes of the rs2268493 SNP demonstrated relatively decreased activation in mesolimbic reward circuitry (including the nucleus accumbens, amygdala, insula, thalamus and prefrontal cortical regions) during the anticipation of rewards but not during the outcome phase of the task. Allelic variation of the rs1042778 and rs237887 SNPs did not moderate mesolimbic activation during either reward anticipation or outcomes.

CONCLUSIONS

This preliminary study suggests that the OXTR SNP rs2268493, which has been previously identified as an ASD risk gene, moderates mesolimbic responses during reward anticipation. Given previous findings of decreased mesolimbic activation during reward anticipation in ASD, the present results suggest that OXTR may confer ASD risk via influences on the neural systems that support reward anticipation.

Stac3 is a component of the excitation-contraction coupling machinery and mutated in Native American myopathy

Excitation-contraction coupling, the process that regulates contractions by skeletal muscles, transduces changes in membrane voltage by activating release of Ca(2+) from internal stores to initiate muscle contraction. Defects in excitation-contraction coupling are associated with muscle diseases. Here we identify Stac3 as a novel component of the excitation-contraction coupling machinery. Using a zebrafish genetic screen, we generate a locomotor mutation that is mapped to stac3. We provide electrophysiological, Ca(2+) imaging, immunocytochemical and biochemical evidence that Stac3 participates in excitation-contraction coupling in muscles. Furthermore, we reveal that a mutation in human STAC3 is the genetic basis of the debilitating Native American myopathy (NAM). Analysis of NAM stac3 in zebrafish shows that the NAM mutation decreases excitation-contraction coupling. These findings enhance our understanding of both excitation-contraction coupling and the pathology of myopathies.

Relationship between methylome and transcriptome in patients with nonalcoholic fatty liver disease

BACKGROUND & AIMS

Cirrhosis and liver cancer are potential outcomes of advanced nonalcoholic fatty liver disease (NAFLD). It is not clear what factors determine whether patients will develop advanced or mild NAFLD, limiting noninvasive diagnosis and treatment before clinical sequelae emerge. We investigated whether DNA methylation profiles can distinguish patients with mild disease from those with advanced NAFLD, and how these patterns are functionally related to hepatic gene expression.

METHODS

We collected frozen liver biopsies and clinical data from patients with biopsy-proven NAFLD (56 in the discovery cohort and 34 in the replication cohort). Samples were divided into groups based on histologic severity of fibrosis: F0-1 (mild) and F3-4 (advanced). DNA methylation profiles were determined and coupled with gene expression data from the same biopsies; differential methylation was validated in subsets of the discovery and replication cohorts. We then analyzed interactions between the methylome and transcriptome.

RESULTS

Clinical features did not differ between patients known to have mild or advanced fibrosis based on biopsy analysis. There were 69,247 differentially methylated CpG sites (76% hypomethylated, 24% hypermethylated) in patients with advanced vs mild NAFLD (P < .05). Methylation at fibroblast growth factor receptor 2, methionine adenosyl methyltransferase 1A, and caspase 1 was validated by bisulfite pyrosequencing and the findings were reproduced in the replication cohort. Methylation correlated with gene transcript levels for 7% of differentially methylated CpG sites, indicating that differential methylation contributes to differences in expression. In samples with advanced NAFLD, many tissue repair genes were hypomethylated and overexpressed, and genes in certain metabolic pathways, including 1-carbon metabolism, were hypermethylated and underexpressed.

CONCLUSIONS

Functionally relevant differences in methylation can distinguish patients with advanced vs mild NAFLD. Altered methylation of genes that regulate processes such as steatohepatitis, fibrosis, and carcinogenesis indicate the role of DNA methylation in progression of NAFLD.

Cadmium exposure and the epigenome: Exposure-associated patterns of DNA methylation in leukocytes from mother-baby pairs

Cadmium (Cd) is prevalent in the environment yet understudied as a developmental toxicant. Cd partially crosses the placental barrier from mother to fetus and is linked to detrimental effects in newborns. Here we examine the relationship between levels of Cd during pregnancy and 5-methylcytosine (5mC) levels in leukocyte DNA collected from 17 mother-newborn pairs. The methylation of cytosines is an epigenetic mechanism known to impact transcriptional signaling and influence health endpoints. A methylated cytosine-guanine (CpG) island recovery assay was used to assess over 4.6 million sites spanning 16,421 CpG islands. Exposure to Cd was classified for each mother-newborn pair according to maternal blood levels and compared with levels of cotinine. Subsets of genes were identified that showed altered DNA methylation levels in their promoter regions in fetal DNA associated with levels of Cd (n = 61), cotinine (n = 366), or both (n = 30). Likewise, in maternal DNA, differentially methylated genes were identified that were associated with Cd (n = 92) or cotinine (n = 134) levels. While the gene sets were largely distinct between maternal and fetal DNA, functional similarities at the biological pathway level were identified including an enrichment of genes that encode for proteins that control transcriptional regulation and apoptosis. Furthermore, conserved DNA motifs with sequence similarity to specific transcription factor binding sites were identified within the CpG islands of the gene sets. This study provides evidence for distinct patterns of DNA methylation or "footprints" in fetal and maternal DNA associated with exposure to Cd.

Genome-wide association study and meta-analysis of intraocular pressure

Elevated intraocular pressure (IOP) is a major risk factor for glaucoma and is influenced by genetic and environmental factors. Recent genome-wide association studies (GWAS) reported associations with IOP at TMCO1 and GAS7, and with primary open-angle glaucoma (POAG) at CDKN2B-AS1, CAV1/CAV2, and SIX1/SIX6. To identify novel genetic variants and replicate the published findings, we performed GWAS and meta-analysis of IOP in >6,000 subjects of European ancestry collected in three datasets: the NEI Glaucoma Human genetics collaBORation, GLAUcoma Genes and ENvironment study, and a subset of the Age-related Macular Degeneration-Michigan, Mayo, AREDS and Pennsylvania study. While no signal achieved genome-wide significance in individual datasets, a meta-analysis identified significant associations with IOP at TMCO1 (rs7518099-G, p = 8.0 × 10(-8)). Focused analyses of five loci previously reported for IOP and/or POAG, i.e., TMCO1, CDKN2B-AS1, GAS7, CAV1/CAV2, and SIX1/SIX6, revealed associations with IOP that were largely consistent across our three datasets, and replicated the previously reported associations in both effect size and direction. These results confirm the involvement of common variants in multiple genomic regions in regulating IOP and/or glaucoma risk.

Interaction of HLA-DRB1*1501 and TNF-Alpha in a Population-based Case-control Study of Multiple Sclerosis

This study was conducted to determine whether single nucleotide polymorphisms (SNPs) in nine genes (human leukocyte antigen (HLA), T cell receptor beta (TCA receptor β), tumor necrosis factor α (TNF α), tumor necrosis factor β (TNF β), apolipoprotein E (APOE), interleukin 7 receptor alpha chain (IL7RA) interleukin 2 receptor alpha chain (IL2RA) myelin basic protein (MBP) and vitamin D receptor (VDR)) associated with multiple sclerosis (MS) could be replicated in a population-based sample, and to determine if these associations are modified by presence of HLA DRB1*1501. DNA was available from 722 individuals (223 with MS and 499 controls) who participated in a population-based case-control study. Cases and controls were matched on ancestry, age, gender and geographic area. HLA DRB1*1501 risk allele (T) was confirmed in this population using a genotypic test, controlling for multiple comparisons. Examining the effect of each SNP in the presence or absence of the HLA DRB1*1501 risk allele identified significant associations with TNF α -1031 (rs1799964) among those without the HLA risk allele. No additional interactions were significant in a cases-only analysis. Our results indicate that an interaction between SNPs in TNF α and HLA DRB1*1501 may influence the risk of developing MS.

Negative life stress and longitudinal hippocampal volume changes in older adults with and without depression

Major depressive disorder is associated with smaller hippocampal volumes but the mechanisms underlying this relationship are unclear. To examine the effect of environmental influences, we examined the relationship between self-reported stressors and two-year change in hippocampal volume. Seventy elderly nondepressed subjects and eighty-nine elderly depressed subjects were followed for two years. The number of negative stressful life events (nSLE), perceived stress levels, and cranial MRI were obtained at baseline and at the two-year assessment. For secondary analyses, subjects provided blood for 5-HTTLPR polymorphism genotyping. After controlling for covariates including presence or absence of depression, greater numbers of baseline nSLEs were significantly associated with greater baseline hippocampal volumes bilaterally. Greater numbers of baseline nSLEs were also associated with reduction in hippocampal volume over two years in the right but not the left hemisphere. Neither perceived stress levels nor changes in stress measures were significantly associated with hippocampal volume measures. However, in secondary analyses, we found that increases in perceived stress over time was associated with volume reduction of the left hippocampus, but only in 5-HTTLPR L/L homozygotes. Our findings suggest different short- and long-term effects of negative life stressors on hippocampal volumes in older adults. These effects appear independent on the presence or absence of depression. Furthermore, these effects may be moderated by genetic polymorphisms in key neurotransmitter systems. These novel findings have important implications for understanding environmental influences on brain aging.

Outcome and life satisfaction of adults with myelomeningocele

BACKGROUND

Myelomeningocele (MMC) commonly causes impairments in body structure and functions as well as cognitive disabilities that can have an adverse effect on adult life. Improved medical care has resulted in increased numbers of individuals with MMC surviving to adulthood, however little is known about the impact of MMC on the lives of adults age 25 years or older.

OBJECTIVE

To gain a better understanding of outcomes in education, employment, relationships, reproduction and life satisfaction of adults with MMC.

METHODS

A primarily quantitative multiple-choice questionnaire designed to capture outcomes in education, employment, relationships and reproduction, along with a previously validated life satisfaction checklist (LiSat-11), was completed by adults with MMC. Relationships between demographic variables, outcomes and life satisfaction were determined using cross tabulation analysis, logistic regression and linear regression.

RESULTS

Ninety adults with MMC, age 25-85 years (median age 32), reported a diverse range of outcomes in education, employment, relationships and reproduction. The most consistent variable associated with difficulty attaining adult milestones was hydrocephalus, the presence of which reduced the likelihood of living independently (p ≤ 0.001), having a partner (p = 0.003) and reproducing (p ≤ 0.001), but did not contribute to reduced life satisfaction.

CONCLUSIONS

Adults with MMC, especially those without hydrocephalus, can obtain gainful employment, live independently, form partner relationships and have children, and these achievements contribute to life satisfaction. While MMC does not affect overall reported life satisfaction for adults, attention should be paid to specific domains with less reported satisfaction.

Genome-wide association analyses identify multiple loci associated with central corneal thickness and keratoconus

Central corneal thickness (CCT) is associated with eye conditions including keratoconus and glaucoma. We performed a meta-analysis on >20,000 individuals in European and Asian populations that identified 16 new loci associated with CCT at genome-wide significance (P < 5 × 10(-8)). We further showed that 2 CCT-associated loci, FOXO1 and FNDC3B, conferred relatively large risks for keratoconus in 2 cohorts with 874 cases and 6,085 controls (rs2721051 near FOXO1 had odds ratio (OR) = 1.62, 95% confidence interval (CI) = 1.4-1.88, P = 2.7 × 10(-10), and rs4894535 in FNDC3B had OR = 1.47, 95% CI = 1.29-1.68, P = 4.9 × 10(-9)). FNDC3B was also associated with primary open-angle glaucoma (P = 5.6 × 10(-4); tested in 3 cohorts with 2,979 cases and 7,399 controls). Further analyses implicate the collagen and extracellular matrix pathways in the regulation of CCT.

Design, methodological issues and participation in a multiple sclerosis case-control study

OBJECTIVES

This study was conducted to determine whether the risk of developing multiple sclerosis (MS) was associated with certain environmental exposures or genetic factors previously reported to influence MS risk. This paper describes the methodological issues, study design and characteristics of the study population.

MATERIALS AND METHODS

Individuals with definite MS were identified from a prevalence study conducted in three geographic areas. The target number of cases was not reached, so an additional study area was added. Identifying clinic controls was inefficient, so controls were recruited using random digit dialing. All study participants completed a detailed questionnaire regarding environmental exposures using computer-assisted telephone interviewing, and blood was collected for genetic analysis.

RESULTS

In total, 276 cases and 590 controls participated, but participation rates were low, ranging from 28.4% to 38.9%. Only one-third (33.6%) of individuals identified in the prevalence study agreed to participate in the case-control study. Cases were more likely to be non-Hispanic white and older than their source populations as identified in the preceding prevalence study (P < 0.05). Most participants provided a blood sample for genotyping (91%; n = 789).

CONCLUSIONS

Epidemiological studies play a key role in identifying genetic and environmental factors that are associated with complex diseases like MS. Methodological issues arise in every study, and investigators need to be able to detect, respond to and correct problems in a timely and scientifically valid manner.

Bipolar disorder, brain-derived neurotrophic factor (BDNF) Val66Met polymorphism and brain morphology

In this study of the effect of bipolar status and presence of BDNF Val66Met polymorphism on differences in regional brain volumes, we hypothesized based on previous studies that 1) bipolar subjects will have smaller regional brain volumes than healthy controls; 2) BDNF Met66 allele carriers within the same population are likely to have smaller regional brain volumes as compared to Val66 homozygyotes. In our Caucasian sample of 166 bipolar subjects and 64 gender-matched healthy controls, we found significant decreases in total (p = 0.005) and regional gray matter volumes in bipolar patients compared to healthy controls, more pronounced in the inferior and posterior parts of the brain, together with a concomitant increase in total CSF (p = 0.012) particularly in the lateral ventricles (p = 0.023). However, there was no difference in white matter volumes noted by other studies. Furthermore we did not find significant differences in other brain regions that have been reported by other authors. Nor did we find a significant effect of BDNF on these measurements.

Genome-wide analysis of central corneal thickness in primary open-angle glaucoma cases in the NEIGHBOR and GLAUGEN consortia

PURPOSE

To investigate the effects of central corneal thickness (CCT)-associated variants on primary open-angle glaucoma (POAG) risk using single nucleotide polymorphisms (SNP) data from the Glaucoma Genes and Environment (GLAUGEN) and National Eye Institute (NEI) Glaucoma Human Genetics Collaboration (NEIGHBOR) consortia.

METHODS

A replication analysis of previously reported CCT SNPs was performed in a CCT dataset (n = 1117) and these SNPs were then tested for association with POAG using a larger POAG dataset (n = 6470). Then a CCT genome-wide association study (GWAS) was performed. Top SNPs from this analysis were selected and tested for association with POAG. cDNA libraries from fetal and adult brain and ocular tissue samples were generated and used for candidate gene expression analysis.

RESULTS

Association with one of 20 previously published CCT SNPs was replicated: rs12447690, near the ZNF469 gene (P = 0.001; β = -5.08 μm/allele). None of these SNPs were significantly associated with POAG. In the CCT GWAS, no SNPs reached genome-wide significance. After testing 50 candidate SNPs for association with POAG, one SNP was identified, rs7481514 within the neurotrimin (NTM) gene, that was significantly associated with POAG in a low-tension subset (P = 0.00099; Odds Ratio [OR] = 1.28). Additionally, SNPs in the CNTNAP4 gene showed suggestive association with POAG (top SNP = rs1428758; P = 0.018; OR = 0.84). NTM and CNTNAP4 were shown to be expressed in ocular tissues.

CONCLUSIONS

The results suggest previously reported CCT loci are not significantly associated with POAG susceptibility. By performing a quantitative analysis of CCT and a subsequent analysis of POAG, SNPs in two cell adhesion molecules, NTM and CNTNAP4, were identified and may increase POAG susceptibility in a subset of cases.