Genome Wide Association Study of Neuropathic Ocular Pain

Purpose

To conduct a genome-wide association study (GWAS) of individuals with neuropathic ocular pain (NOP) symptoms to identify genomic variants that may predispose to NOP development.

Design

Prospective study of individuals with NOP.

Participants

Three hundred twenty-nine patients recruited from the Miami Veterans Affairs eye clinic.

Methods

The Neuropathic Pain Symptom Inventory modified for the eye (NPSI-Eye) was completed to calculate a NPSI-Eye-Sub-Score (summed ratings of burning and wind sensitivity) as an indicator of NOP severity. A GWAS was performed for the NPSI-Eye-Sub-Score with a significance threshold of P < 5 × 10-8. A gene-based analysis was performed using the multimarker analysis of genomic annotation software (in the functional mapping and annotation of GWAS online platform). The 13 865 778 single nucleotide polymorphisms (SNPs) from our GWAS analysis were mapped to 10 834 protein coding genes, and significant genes were run through gene set enrichment analysis.

Main Outcome Measures

Identification of SNPs and protein products that may be associated with the development of NOP.

Results

One hundred seventy-one SNPs reached a threshold of P < 10-5, of which 10 SNPs reached the suggestive level of significance of P < 5 × 10-7 and 1 SNP met our genome-wide significance threshold of P < 5 × 10-8. This lead SNP, rs140293404 (P = 1.23 × 10-8), is an intronic variant found within gene ENSG00000287251 coding for transcript ENST00000662732.1. Rs140293404 is in linkage disequilibrium with exon variant rs7926353 (r2 > 0.8) within ENSG00000279046 coding for transcript ENST00000624288.1. The most significant genes from gene-based tests were matrix metalloproteinase-19 (MMP19) (P = 1.12 × 10-5), zinc finger RNA-binding motif and serine/arginine rich-1 (ZRSR1) (P = 1.48 × 10-4), CTC-487M23.8 (P = 1.79 × 10-4), receptor expression-enhancing protein-5 (REEP5) (P = 2.36 × 10-4), and signal recognition particle-19 (SRP19) (P = 2.56 × 10-4). From gene set enrichment analysis, the sensory perception (false discovery rate = 6.57 × 10-3) and olfactory signaling (false discovery rate = 1.63 × 10-2) pathways were enriched with the most significant genes.

Conclusions

Our GWAS revealed genes with protein products that may impact sensory perception, lending biological plausibility to a role for SNPs identified by our GWAS in the development of NOP. A better understanding of the biological relevance of these genes and pathways in the pathophysiology associated with NOP may facilitate future novel mechanism-based treatments.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Founder population-specific weights yield improvements in performance of polygenic risk scores for Alzheimer disease in the Midwestern Amish

Alzheimer disease (AD) is the most common type of dementia and is estimated to affect 6 million Americans. Risk for AD is multifactorial, including both genetic and environmental risk factors. AD genomic research has generally focused on identification of risk variants. Using this information, polygenic risk scores (PRSs) can be calculated to quantify an individual's relative disease risk due to genetic factors. The Amish are a founder population descended from German and Swiss Anabaptist immigrants. They experienced a genetic bottleneck after arrival in the United States, making their genetic architecture different from the broader European ancestry population. Prior work has demonstrated the lack of transferability of PRSs across populations. Here, we compared the performance of PRSs derived from genome-wide association studies (GWASs) of Amish individuals to those derived from a large European ancestry GWAS. Participants were screened for cognitive impairment with further evaluation for AD. Genotype data were imputed after collection via Illumina genotyping arrays. The Amish individuals were split into two groups based on the primary site of recruitment. For each group, GWAS was conducted with account for relatedness and adjustment for covariates. PRSs were then calculated using weights from the other Amish group. PRS models were evaluated with and without covariates. The Amish-derived PRSs distinguished between dementia status better than the European-derived PRS in our Amish populations and demonstrated performance improvements despite a smaller training sample size. This work highlighted considerations for AD PRS usage in populations that cannot be adequately described by basic race/ethnicity or ancestry classifications.

Multi-ancestry genome-wide meta-analysis of 56,241 individuals identifies LRRC4C, LHX5-AS1 and nominates ancestry-specific loci PTPRK , GRB14 , and KIAA0825 as novel risk loci for Alzheimer's disease: the Alzheimer's Disease Genetics Consortium

Limited ancestral diversity has impaired our ability to detect risk variants more prevalent in non-European ancestry groups in genome-wide association studies (GWAS). We constructed and analyzed a multi-ancestry GWAS dataset in the Alzheimer's Disease (AD) Genetics Consortium (ADGC) to test for novel shared and ancestry-specific AD susceptibility loci and evaluate underlying genetic architecture in 37,382 non-Hispanic White (NHW), 6,728 African American, 8,899 Hispanic (HIS), and 3,232 East Asian individuals, performing within-ancestry fixed-effects meta-analysis followed by a cross-ancestry random-effects meta-analysis. We identified 13 loci with cross-ancestry associations including known loci at/near CR1 , BIN1 , TREM2 , CD2AP , PTK2B , CLU , SHARPIN , MS4A6A , PICALM , ABCA7 , APOE and two novel loci not previously reported at 11p12 ( LRRC4C ) and 12q24.13 ( LHX5-AS1 ). Reflecting the power of diverse ancestry in GWAS, we observed the SHARPIN locus using 7.1% the sample size of the original discovering single-ancestry GWAS (n=788,989). We additionally identified three GWS ancestry-specific loci at/near ( PTPRK ( P =2.4×10 -8 ) and GRB14 ( P =1.7×10 -8 ) in HIS), and KIAA0825 ( P =2.9×10 -8 in NHW). Pathway analysis implicated multiple amyloid regulation pathways (strongest with P adjusted =1.6×10 -4 ) and the classical complement pathway ( P adjusted =1.3×10 -3 ). Genes at/near our novel loci have known roles in neuronal development ( LRRC4C, LHX5-AS1 , and PTPRK ) and insulin receptor activity regulation ( GRB14 ). These findings provide compelling support for using traditionally-underrepresented populations for gene discovery, even with smaller sample sizes.

Visuospatial and Verbal Memory Differences in Amish Individuals With Alzheimer Disease and Related Dementias

BACKGROUND

Verbal and visuospatial memory impairments are common to Alzheimer disease and Related Dementias (ADRD), but the patterns of decline in these domains may reflect genetic and lifestyle influences. The latter may be pertinent to populations such as the Amish who have unique lifestyle experiences.

METHODS

Our data set included 420 Amish and 401 CERAD individuals. Sex-adjusted, age-adjusted, and education-adjusted Z-scores were calculated for the recall portions of the Constructional Praxis Delay (CPD) and Word List Delay (WLD). ANOVAs were then used to examine the main and interaction effects of cohort (Amish, CERAD), cognitive status (case, control), and sex on CPD and WLD Z-scores.

RESULTS

The Amish performed better on the CPD than the CERAD cohort. In addition, the difference between cases and controls on the CPD and WLD were smaller in the Amish and Amish female cases performed better on the WLD than the CERAD female cases.

DISCUSSION

The Amish performed better on the CPD task, and ADRD-related declines in CPD and WLD were less severe in the Amish. In addition, Amish females with ADRD may have preferential preservation of WLD. This study provides evidence that the Amish exhibit distinct patterns of verbal and visuospatial memory loss associated with aging and ADRD.

Identification of novel genes for age-at-onset of Alzheimer's disease by combining quantitative and survival trait analyses

INTRODUCTION

Our understanding of the genetic predisposition for age-at-onset (AAO) of Alzheimer's disease (AD) is limited. Here, we sought to identify genes modifying AAO and examined whether any have sex-specific effects.

METHODS

Genome-wide association analysis were performed on imputed genetic data of 9219 AD cases and 10,345 controls from 20 cohorts of the Alzheimer's Disease Genetics Consortium. AAO was modeled from cases directly and as a survival outcome.

RESULTS

We identified 11 genome-wide significant loci (P < 5 × 10-8 ), including six known AD-risk genes and five novel loci, UMAD1, LUZP2, ARFGEF2, DSCAM, and 4q25, affecting AAO of AD. Additionally, 39 suggestive loci showed strong association. Twelve loci showed sex-specific effects on AAO including CD300LG and MLX/TUBG2 for females and MIR4445 for males.

DISCUSSION

Genes that influence AAO of AD are excellent therapeutic targets for delaying onset of AD. Several loci identified include genes with promising functional implications for AD.

Genetic variants in the SHISA6 gene are associated with delayed cognitive impairment in two family datasets

INTRODUCTION

Studies of cognitive impairment (CI) in Amish communities have identified sibships containing CI and cognitively unimpaired (CU) individuals. We hypothesize that CU individuals may carry protective alleles delaying age at onset (AAO) of CI.

METHODS

A total of 1522 individuals screened for CI were genotyped. The outcome studied was AAO for CI individuals or age at last normal exam for CU individuals. Cox mixed-effects models examined association between age and single nucleotide variants (SNVs).

RESULTS

Three SNVs were significantly associated (P < 5 × 10-8 ) with AAO on chromosomes 6 (rs14538074; hazard ratio [HR] = 3.35), 9 (rs534551495; HR = 2.82), and 17 (rs146729640; HR = 6.38). The chromosome 17 association was replicated in the independent National Institute on Aging Genetics Initiative for Late-Onset Alzheimer's Disease dataset.

DISCUSSION

The replicated genome-wide significant association with AAO on chromosome 17 is located in the SHISA6 gene, which is involved in post-synaptic transmission in the hippocampus and is a biologically plausible candidate gene for Alzheimer's disease.

The genetic architecture of Alzheimer disease risk in the Ohio and Indiana Amish

Alzheimer disease (AD) is the most common type of dementia and is currently estimated to affect 6.2 million Americans. It ranks as the sixth leading cause of death in the United States, and the proportion of deaths due to AD has been increasing since 2000, while the proportion of many other leading causes of deaths have decreased or remained constant. The risk for AD is multifactorial, including genetic and environmental risk factors. Although APOE ε4 remains the largest genetic risk factor for AD, more than 26 other loci have been associated with AD risk. Here, we recruited Amish adults from Ohio and Indiana to investigate AD risk and protective genetic effects. As a founder population that typically practices endogamy, variants that are rare in the general population may be of a higher frequency in the Amish population. Since the Amish have a slightly lower incidence and later age of onset of disease, they represent an excellent and unique population for research on protective genetic variants. We compared AD risk in the Amish and to a non-Amish population through APOE genotype, a non-APOE genetic risk score of genome-wide significant variants, and a non-APOE polygenic risk score considering all of the variants. Our results highlight the lesser relative impact of APOE and differing genetic architecture of AD risk in the Amish compared to a non-Amish, general European ancestry population.

Association Between Polymorphisms in DNA Damage Repair Genes and Radiation Therapy-Induced Early Adverse Skin Reactions in a Breast Cancer Population: A Polygenic Risk Score Approach

PURPOSE

Genetic variations in DNA damage repair (DDR) genes may influence radiation therapy (RT)-induced acute normal tissue toxicity in patients with breast cancer. Identifying an individual or multiple single-nucleotide polymorphisms (SNPs) associated with RT-induced early adverse skin reactions (EASR) is critical for precision medicine in radiation oncology.

METHODS AND MATERIALS

At the completion of RT, EASR was assessed using the Oncology Nursing Society scale (0-6) in 416 patients with breast cancer, and Oncology Nursing Society score ≥4 was considered RT-induced EASR. PLINK set-based tests and subsequent individual SNP association analyses were conducted to identify genes and SNPs associated with EASR among the 53 DDR genes and 1968 SNPs. A weighted polygenic risk score (PRS) model was constructed to ascertain the association between the joint effect of risk alleles and EASR.

RESULTS

The study population consisted of 264 Hispanic whites, 86 blacks or African Americans, 55 non-Hispanic whites, and 11 others. A total of 115 patients (27.6%) developed EASR. Five genes (ATM, CHEK1, ERCC2, RAD51C, and TGFB1) were significantly associated with RT-induced EASR. Nine SNPs within these 5 genes were further identified: ATM rs61915066, CHEK1 rs11220184, RAD51C rs302877, rs405684, TBFB1 rs4803455, rs2241714, and ERCC2 rs60152947, rs10404465, rs1799786. In a multivariable-adjusted PRS model, patients in a higher quartile of PRS were more likely to develop EASR compared with patients in the lowest quartile (OR_{q2 vs.q1} = 1.94, 95% CI, 0.86-4.39; OR_{q3 vs.q1} = 3.46, 95% CI, 1.57-7.63; OR_{q4 vs.q1} = 8.64, 95% CI, 3.92-19.02; and P_trend < .0001).

CONCLUSIONS

We newly identified the associations between 9 SNPs in ATM, CHEK1, RAD51C, TGFB1, and ERCC2 and RT-induced EASR. PRS modeling showed its potential in identifying populations at risk. Multiple SNPs in DDR genes may jointly contribute to interindividual variation in RT-induced EASR. Validation in an independent external cohort is required to determine the clinical significance of these predictive biomarkers.

Genome-wide enriched pathway analysis of acute post-radiotherapy pain in breast cancer patients: a prospective cohort study

BACKGROUND

Adjuvant radiotherapy (RT) can increase the risk of developing pain; however, the molecular mechanisms of RT-related pain remain unclear. The current study aimed to identify susceptibility loci and enriched pathways for clinically relevant acute post-RT pain, defined as having moderate to severe pain (pain score ≥ 4) at the completion of RT.

METHODS

We conducted a genome-wide association study (GWAS) with 1,344,832 single-nucleotide polymorphisms (SNPs), a gene-based analysis using PLINK set-based tests of 19,621 genes, and a functional enrichment analysis of a gene list of 875 genes with p < 0.05 using NIH DAVID functional annotation module with KEGG pathways and GO terms (n = 380) among 1112 breast cancer patients.

RESULTS

About 29% of patients reported acute post-RT pain. None of SNPs nor genes reached genome-wide significant level. Four SNPs showed suggestive associations with post-RT pain; rs16970540 in RFFL or near the LIG3 gene (p = 1.7 × 10-6), rs4584690, and rs7335912 in ABCC4/MPR4 gene (p = 5.5 × 10-6 and p = 7.8 × 10-6, respectively), and rs73633565 in EGFL6 gene (p = 8.1 × 10-6). Gene-based analysis suggested the potential involvement of neurotransmitters, olfactory receptors, and cytochrome P450 in post-RT pain, whereas functional analysis showed glucuronidation (FDR-adjusted p value = 9.46 × 10-7) and olfactory receptor activities (FDR-adjusted p value = 0.032) as the most significantly enriched biological features.

CONCLUSIONS

This is the first GWAS suggesting that post-RT pain is a complex polygenic trait influenced by many biological processes and functions such as glucuronidation and olfactory receptor activities. If validated in larger populations, the results can provide biological targets for pain management to improve cancer patients' quality of life. Additionally, these genes can be further tested as predictive biomarkers for personalized pain management.

PLINK: Key Functions for Data Analysis

Genetic data analysis of large numbers of single nucleotide variants (SNVs), including genome-wide association studies (GWAS), exome chips, and whole exome (WES) or whole-genome (WGS) sequencing data, requires well defined processing steps. As a result, several freely available analytic toolkits have been developed to streamline these processes. Among these, PLINK is the most comprehensive in terms of its quality control and analytic modules, although its focus remains on SNVs. PLINK fulfills two analytic needs-aiding the process of performing quality control (QC) on large data sets and providing basic statistical tools to analyze the variants in genetic models. The current version of PLINK (v1.90b) has incorporated several sophisticated statistical modeling features, such as those that were introduced by GCTA (genome-wide complex trait analysis), including mixed-model association analysis and cluster-based algorithms. Although PLINK is diverse in its applicability to data management and analysis, in some instances, other available tools offer more optimal options. Here we provide a practical overview of major PLINK features with respect to QC, data management, and association mapping, along with learned shortcuts and limitations to be considered. In cases where PLINK features are limited, we provide alternative approaches using additional freely available pipelines. © 2018 by John Wiley & Sons, Inc.

Properties of global- and local-ancestry adjustments in genetic association tests in admixed populations

Population substructure can lead to confounding in tests for genetic association, and failure to adjust properly can result in spurious findings. Here we address this issue of confounding by considering the impact of global ancestry (average ancestry across the genome) and local ancestry (ancestry at a specific chromosomal location) on regression parameters and relative power in ancestry-adjusted and -unadjusted models. We examine theoretical expectations under different scenarios for population substructure; applying different regression models, verifying and generalizing using simulations, and exploring the findings in real-world admixed populations. We show that admixture does not lead to confounding when the trait locus is tested directly in a single admixed population. However, if there is more complex population structure or a marker locus in linkage disequilibrium (LD) with the trait locus is tested, both global and local ancestry can be confounders. Additionally, we show the genotype parameters of adjusted and unadjusted models all provide tests for LD between the marker and trait locus, but in different contexts. The local ancestry adjusted model tests for LD in the ancestral populations, while tests using the unadjusted and the global ancestry adjusted models depend on LD in the admixed population(s), which may be enriched due to different ancestral allele frequencies. Practically, this implies that global-ancestry adjustment should be used for screening, but local-ancestry adjustment may better inform fine mapping and provide better effect estimates at trait loci.

Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants

Background

Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD.

Methods

We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD.

Results

We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1.

Conclusions

We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-015-0034-z) contains supplementary material, which is available to authorized users.

Exome sequencing of extended families with autism reveals genes shared across neurodevelopmental and neuropsychiatric disorders

BACKGROUND

Autism spectrum disorders (ASDs) comprise a range of neurodevelopmental conditions of varying severity, characterized by marked qualitative difficulties in social relatedness, communication, and behavior. Despite overwhelming evidence of high heritability, results from genetic studies to date show that ASD etiology is extremely heterogeneous and only a fraction of autism genes have been discovered.

METHODS

To help unravel this genetic complexity, we performed whole exome sequencing on 100 ASD individuals from 40 families with multiple distantly related affected individuals. All families contained a minimum of one pair of ASD cousins. Each individual was captured with the Agilent SureSelect Human All Exon kit, sequenced on the Illumina Hiseq 2000, and the resulting data processed and annotated with Burrows-Wheeler Aligner (BWA), Genome Analysis Toolkit (GATK), and SeattleSeq. Genotyping information on each family was utilized in order to determine genomic regions that were identical by descent (IBD). Variants identified by exome sequencing which occurred in IBD regions and present in all affected individuals within each family were then evaluated to determine which may potentially be disease related. Nucleotide alterations that were novel and rare (minor allele frequency, MAF, less than 0.05) and predicted to be detrimental, either by altering amino acids or splicing patterns, were prioritized.

RESULTS

We identified numerous potentially damaging, ASD associated risk variants in genes previously unrelated to autism. A subset of these genes has been implicated in other neurobehavioral disorders including depression (SLIT3), epilepsy (CLCN2, PRICKLE1), intellectual disability (AP4M1), schizophrenia (WDR60), and Tourette syndrome (OFCC1). Additional alterations were found in previously reported autism candidate genes, including three genes with alterations in multiple families (CEP290, CSMD1, FAT1, and STXBP5). Compiling a list of ASD candidate genes from the literature, we determined that variants occurred in ASD candidate genes 1.65 times more frequently than in random genes captured by exome sequencing (P = 8.55 × 10-5).

CONCLUSIONS

By studying these unique pedigrees, we have identified novel DNA variations related to ASD, demonstrated that exome sequencing in extended families is a powerful tool for ASD candidate gene discovery, and provided further evidence of an underlying genetic component to a wide range of neurodevelopmental and neuropsychiatric diseases.

A common MUC5B promoter polymorphism and pulmonary fibrosis

BACKGROUND

The mutations that have been implicated in pulmonary fibrosis account for only a small proportion of the population risk.

METHODS

Using a genomewide linkage scan, we detected linkage between idiopathic interstitial pneumonia and a 3.4-Mb region of chromosome 11p15 in 82 families. We then evaluated genetic variation in this region in gel-forming mucin genes expressed in the lung among 83 subjects with familial interstitial pneumonia, 492 subjects with idiopathic pulmonary fibrosis, and 322 controls. MUC5B expression was assessed in lung tissue.

RESULTS

Linkage and fine mapping were used to identify a region of interest on the p-terminus of chromosome 11 that included gel-forming mucin genes. The minor-allele of the single-nucleotide polymorphism (SNP) rs35705950, located 3 kb upstream of the MUC5B transcription start site, was present at a frequency of 34% among subjects with familial interstitial pneumonia, 38% among subjects with idiopathic pulmonary fibrosis, and 9% among controls (allelic association with familial interstitial pneumonia, P=1.2×10(-15); allelic association with idiopathic pulmonary fibrosis, P=2.5×10(-37)). The odds ratios for disease among subjects who were heterozygous and those who were homozygous for the minor allele of this SNP were 6.8 (95% confidence interval [CI], 3.9 to 12.0) and 20.8 (95% CI, 3.8 to 113.7), respectively, for familial interstitial pneumonia and 9.0 (95% CI, 6.2 to 13.1) and 21.8 (95% CI, 5.1 to 93.5), respectively, for idiopathic pulmonary fibrosis. MUC5B expression in the lung was 14.1 times as high in subjects who had idiopathic pulmonary fibrosis as in those who did not (P<0.001). The variant allele of rs35705950 was associated with up-regulation in MUC5B expression in the lung in unaffected subjects (expression was 37.4 times as high as in unaffected subjects homozygous for the wild-type allele, P<0.001). MUC5B protein was expressed in lesions of idiopathic pulmonary fibrosis.

CONCLUSIONS

A common polymorphism in the promoter of MUC5B is associated with familial interstitial pneumonia and idiopathic pulmonary fibrosis. Our findings suggest that dysregulated MUC5B expression in the lung may be involved in the pathogenesis of pulmonary fibrosis. (Funded by the National Heart, Lung, and Blood Institute and others.).

A modifier locus on chromosome 5 contributes to L1 cell adhesion molecule X-linked hydrocephalus in mice

Humans with L1 cell adhesion molecule (L1CAM) mutations exhibit X-linked hydrocephalus, as well as other severe neurological disorders. L1-6D mutant mice, which are homozygous for a deletion that removes the sixth immunoglobulin-like domain of L1cam, seldom display hydrocephalus on the 129/Sv background. However, the same L1-6D mutation produces severe hydrocephalus on the C57BL/6J background. To begin to understand how L1cam deficiencies result in hydrocephalus and to identify modifier loci that contribute to X-linked hydrocephalus by genetically interacting with L1cam, we conducted a genome-wide scan on F2 L1-6D mice, bred from L1-6D 129S2/SvPasCrlf and C57BL/6J mice. Linkage studies, utilizing chi-square tests and quantitative trait loci mapping techniques, were performed. Candidate modifier loci were further investigated in an extension study. Linkage was confirmed for a locus on chromosome 5, which we named L1cam hydrocephalus modifier 1 (L1hydro1), p = 4.04 X 10(-11).

A genome-wide association study of autism reveals a common novel risk locus at 5p14.1

Although autism is one of the most heritable neuropsychiatric disorders, its underlying genetic architecture has largely eluded description. To comprehensively examine the hypothesis that common variation is important in autism, we performed a genome-wide association study (GWAS) using a discovery dataset of 438 autistic Caucasian families and the Illumina Human 1M beadchip. 96 single nucleotide polymorphisms (SNPs) demonstrated strong association with autism risk (p-value < 0.0001). The validation of the top 96 SNPs was performed using an independent dataset of 487 Caucasian autism families genotyped on the 550K Illumina BeadChip. A novel region on chromosome 5p14.1 showed significance in both the discovery and validation datasets. Joint analysis of all SNPs in this region identified 8 SNPs having improved p-values (3.24E-04 to 3.40E-06) than in either dataset alone. Our findings demonstrate that in addition to multiple rare variations, part of the complex genetic architecture of autism involves common variation.

Association of AGTR1 with 18-month treatment outcome in late-life depression

OBJECTIVE

Converging lines of evidence implicate vascular factors in late-life depression, and argue that late-life depression is a distinct entity among the mood disorders. The A1166C polymorphism in the angiotensin II receptor, vascular type 1 (AGTR1) gene has been associated with a range of vascular diseases. This study investigated the association of AGTR1 genotype on 18-month treatment outcome in late-life depression.

METHODS

In a large, prospective cohort study, patients with late-life depression received individualized treatment using a standardized algorithm. The authors genotyped participants at the AGTR1 A1166C single nucleotide polymorphism (SNP) using standardized methodology, then used survival analysis to estimate the impact of A1166C and demographic variables on time to remission during 18 months of follow-up.

RESULTS

The hazard ratio for AGTR1 homozygous C/C status was 0.37. The A1166C SNP showed evidence for genotypic and allelic association in a comparison of remitted and unremitted/censored subjects.

CONCLUSION

Consistent with its association with numerous vascular disorders, AGTR1 is associated with treatment outcome in late-life depression. Further studies are needed to replicate this finding, and to investigate the impact of other genetic markers of vascular disease on late-life depression outcome.

Phenotypic definition of Chiari type I malformation coupled with high-density SNP genome screen shows significant evidence for linkage to regions on chromosomes 9 and 15

Chiari type I malformation (CMI; OMIM 118420) is narrowly defined when the tonsils of the cerebellum extend below the foramen magnum, leading to a variety of neurological symptoms. It is widely thought that a small posterior fossa (PF) volume, relative to the total cranial volume leads to a cramped cerebellum and herniation of the tonsils into the top of the spinal column. In a collection of magnetic resonance imagings (MRIs) from affected individuals and their family members, we measured correlations between ten cranial morphologies and estimated their heritability in these families. Correlations between bones delineating the PF and significant heritability of PF volume (0.955, P = 0.003) support the cramped PF theory and a genetic basis for this condition. In a collection of 23 families with 71 affected individuals, we performed a genome wide linkage screen of over 10,000 SNPs across the genome to identify regions of linkage to CMI. Two-point LOD scores on chromosome 15 reached 3.3 and multipoint scores in this region identified a 13 cM region with LOD scores over 1 (15q21.1-22.3). This region contains a biologically plausible gene for CMI, fibrillin-1, which is a major gene in Marfan syndrome and has been linked to Shprintzen-Goldberg syndrome, of which CMI is a distinguishing characteristic. Multipoint LOD scores on chromosome 9 maximized at 3.05, identifying a 40 cM region with LOD scores over 1 (9q21.33-33.1) and a tighter region with multipoint LOD scores over 2 that was only 8.5 cM. This linkage evidence supports a genetic role in Chiari malformation and justifies further exploration with fine mapping and investigation of candidate genes in these regions.

High-density single nucleotide polymorphism screen in a large multiplex neural tube defect family refines linkage to loci at 7p21.1-pter and 2q33.1-q35

BACKGROUND

Neural tube defects (NTDs) are considered complex, with both genetic and environmental factors implicated. To date, no major causative genes have been identified in humans despite several investigations. The first genomewide screen in NTDs demonstrated evidence of linkage to chromosomes 7 and 10. This screen included 44 multiplex families and consisted of 402 microsatellite markers spaced approximately 10 cM apart. Further investigation of the genomic screen data identified a single large multiplex family, pedigree 8776, as primarily driving the linkage results on chromosome 7.

METHODS

To investigate this family more thoroughly, a high-density single nucleotide polymorphism (SNP) screen was performed. Two-point and multipoint linkage analyses were performed using both parametric and nonparametric methods.

RESULTS

For both the microsatellite and SNP markers, linkage analysis suggested the involvement of a locus or loci proximal to the telomeric regions of chromosomes 2q and 7p, with both regions generating a LOD* score of 3.0 using a nonparametric identity by descent relative sharing method.

CONCLUSIONS

The regions with the strongest evidence for linkage map proximal to the telomeres on these two chromosomes. In addition to mutations and/or variants in a major gene, these loci may harbor a microdeletion and/or translocation; potentially, polygenic factors may also be involved. This single family may be promising for narrowing the search for NTD susceptibility genes.

Neural tube defects and folate pathway genes: family-based association tests of gene-gene and gene-environment interactions

BACKGROUND

Folate metabolism pathway genes have been examined for association with neural tube defects (NTDs) because folic acid supplementation reduces the risk of this debilitating birth defect. Most studies addressed these genes individually, often with different populations providing conflicting results.

OBJECTIVES

Our study evaluates several folate pathway genes for association with human NTDs, incorporating an environmental cofactor: maternal folate supplementation.

METHODS

In 304 Caucasian American NTD families with myelomeningocele or anencephaly, we examined 28 polymorphisms in 11 genes: folate receptor 1, folate receptor 2, solute carrier family 19 member 1, transcobalamin II, methylenetetrahydrofolate dehydrogenase 1, serine hydroxymethyl-transferase 1, 5,10-methylenetetrahydrofolate reductase (MTHFR), 5-methyltetrahydrofolate-homo-cysteine methyltransferase, 5-methyltetrahydrofolate-homocysteine methyltransferase reductase, betaine-homocysteine methyltransferase (BHMT), and cystathionine-beta-synthase.

RESULTS

Only single nucleotide polymorphisms (SNPs) in BHMT were significantly associated in the overall data set; this significance was strongest when mothers took folate-containing nutritional supplements before conception. The BHMT SNP rs3733890 was more significant when the data were stratified by preferential transmission of the MTHFR rs1801133 thermolabile T allele from parent to offspring. Other SNPs in folate pathway genes were marginally significant in some analyses when stratified by maternal supplementation, MTHFR, or BHMT allele transmission.

CONCLUSIONS

BHMT rs3733890 is significantly associated in our data set, whereas MTHFR rs1801133 is not a major risk factor. Further investigation of folate and methionine cycle genes will require extensive SNP genotyping and/or resequencing to identify novel variants, inclusion of environmental factors, and investigation of gene-gene interactions in large data sets.

A genome-wide linkage analysis of dementia in the Amish

Susceptibility genes for Alzheimer's disease are proving to be highly challenging to detect and verify. Population heterogeneity may be a significant confounding factor contributing to this difficulty. To increase the power for disease susceptibility gene detection, we conducted a genome-wide genetic linkage screen using individuals from the relatively isolated, genetically homogeneous, Amish population. Our genome linkage analysis used a 407-microsatellite-marker map (average density 7 cM) to search for autosomal genes linked to dementia in five Amish families from four Midwestern U.S. counties. Our highest two-point lod score (3.01) was observed at marker D4S1548 on chromosome 4q31. Five other regions (10q22, 3q28, 11p13, 4q28, 19p13) also demonstrated suggestive linkage with markers having two-point lod scores >2.0. While two of these regions are novel (4q31 and 11p13), the other regions lie close to regions identified in previous genome scans in other populations. Our results identify regions of the genome that may harbor genes involved in a subset of dementia patients, in particular the North American Amish community.

Clinical and pathologic features of familial interstitial pneumonia

RATIONALE

Several lines of evidence suggest that genetic factors and environmental exposures play a role in the development of pulmonary fibrosis.

OBJECTIVES

We evaluated families with 2 or more cases of idiopathic interstitial pneumonia among first-degree family members (familial interstitial pneumonia, or FIP), and identified 111 families with FIP having 309 affected and 360 unaffected individuals.

METHODS

The presence of probable or definite FIP was based on medical record review in 28 cases (9.1%); clinical history, diffusing capacity of carbon monoxide (DL(CO)), and chest X-ray in 16 cases (5.2%); clinical history, DL(CO), and high-resolution computed tomography chest scan in 191 cases (61.8%); clinical history and surgical lung biopsy in 56 cases (18.1%); and clinical history and autopsy in 18 cases (5.8%).

RESULTS

Older age (68.3 vs. 53.1; p < 0.0001), male sex (55.7 vs. 37.2%; p < 0.0001), and having ever smoked cigarettes (67.3 vs. 34.1%; p < 0.0001) were associated with the development of FIP. After controlling for age and sex, having ever smoked cigarettes remained strongly associated with the development of FIP (odds ratio(adj), 3.6; 95% confidence interval, 1.3-9.8). Evidence of aggregation of disease was highly significant (p < 0.001) among sibling pairs, and 20 pedigrees demonstrated vertical transmission, consistent with autosomal dominant inheritance. Forty-five percent of pedigrees demonstrated phenotypic heterogeneity, with some pedigrees demonstrating several subtypes of idiopathic interstitial pneumonia occurring within the same families.

CONCLUSIONS

These findings suggest that FIP may be caused by an interaction between a specific environmental exposure and a gene (or genes) that predisposes to the development of several subtypes of idiopathic interstitial pneumonia.

SNPs in the neural cell adhesion molecule 1 gene (NCAM1) may be associated with human neural tube defects

Neural tube defects (NTDs) are common birth defects, occurring in approximately 1/1,000 births; both genetic and environmental factors are implicated. To date, no major genetic risk factors have been identified. Throughout development, cell adhesion molecules are strongly implicated in cell-cell interactions, and may play a role in the formation and closure of the neural tube. To evaluate the role of neural cell adhesion molecule 1 (NCAM1) in risk of human NTDs, we screened for novel single-nucleotide polymorphisms (SNPs) within the gene. Eleven SNPs across NCAM1 were genotyped using TaqMan. We utilized a family-based approach to evaluate evidence for association and/or linkage disequilibrium. We evaluated American Caucasian simplex lumbosacral myelomeningocele families (n=132 families) using the family based association test (FBAT) and the pedigree disequilibrium test (PDT). Association analysis revealed a significant association between risk for NTDs and intronic SNP rs2298526 using both the FBAT test (P=0.0018) and the PDT (P=0.0025). Using the HBAT version of the FBAT to look for haplotype association, all pairwise comparisons with SNP rs2298526 were also significant. A replication study set, consisting of 72 additional families showed no significant association; however, the overall trend for overtransmission of the less common allele of SNP rs2298526 remained significant in the combined sample set. In addition, we analyzed the expression pattern of the NCAM1 protein in human embryos, and while NCAM1 is not expressed within the neural tube at the time of closure, it is expressed in the surrounding and later in differentiated neurons of the CNS. These results suggest variations in NCAM1 may influence risk for human NTDs.

Life after the screen: making sense of many P-values

A multiple analytic approach may be useful for analyzing complex traits since different methods extract both similar and distinct, but complementary pieces of information from genome screen data on extended pedigrees. We examined the usefulness of combining p-values both across methods and across adjacent markers, taking into account the observed correlation structure among these p-values. To this end, we employed the recently proposed truncated product method [Zaykin et al., Genet Epidemiol, in press]. It appears that this approach is helpful for visualizing priority regions for follow-up analysis and reducing the number of false-positive linkage signals.

SNPing away at complex diseases: analysis of single-nucleotide polymorphisms around APOE in Alzheimer disease

There has been great interest in the prospects of using single-nucleotide polymorphisms (SNPs) in the search for complex disease genes, and several initiatives devoted to the identification and mapping of SNPs throughout the human genome are currently underway. However, actual data investigating the use of SNPs for identification of complex disease genes are scarce. To begin to look at issues surrounding the use of SNPs in complex disease studies, we have initiated a collaborative SNP mapping study around APOE, the well-established susceptibility gene for late-onset Alzheimer disease (AD). Sixty SNPs in a 1.5-Mb region surrounding APOE were genotyped in samples of unrelated cases of AD, in controls, and in families with AD. Standard tests were conducted to look for association of SNP alleles with AD, in cases and controls. We also used family-based association analyses, including recently developed methods to look for haplotype association. Evidence of association (P Collapse

Key Words Collapse

MESH Headings

• Age of Onset
• Alleles
• Alzheimer Disease/epidemiology
• Alzheimer Disease/genetics
• Apolipoproteins E/genetics
• Case-Control Studies
• Chromosome Mapping/methods
• Gene Frequency/genetics
• Genetic Predisposition to Disease/genetics
• Genotype
• Haplotypes/genetics
• Humans
• Linkage Disequilibrium/genetics
• Lod Score
• Middle Aged
• Models, Genetic
• Polymorphism, Single Nucleotide/genetics

Collapse

Grants

• 5 R01 NS31153-07 NINDS NIH HHS
• R01 NS031153 NINDS NIH HHS
• 2 P01 NS26630-11A1 NINDS NIH HHS
• 5 P50 AG05128-16 NIA NIH HHS
• Wellcome Trust
• P01 NS026630 NINDS NIH HHS
• P50 AG005128 NIA NIH HHS
• U24 AG021886 NIA NIH HHS

Collapse

A genetic linkage map of the baboon (Papio hamadryas) genome based on human microsatellite polymorphisms

A first-generation genetic linkage map of the baboon (Papio hamadryas) genome was developed for use in biomedical and evolutionary genetics. Pedigreed baboons (n = 694) were selected from the breeding colony maintained by the Southwest Foundation for Biomedical Research. To facilitate comparison with the human genome, the baboon linkage map consists primarily of human microsatellite loci amplified using published human PCR primers. Genotypes for 325 human microsatellites and 6 novel baboon microsatellites were used in linkage analyses performed with the MultiMap expert system. The resulting sex-averaged meiotic recombination map covers all 20 baboon autosomes, with average spacing among loci of 7.2 cM. Direct comparison among homologous (orthologous) loci reveals that, for 7 human autosomes, locus order is conserved between humans and baboons. For the other 15 autosomes, one or more rearrangements distinguish the two genomes. The total centimorgan distances among homologous markers are 28.0% longer in the human genome than in the baboon, suggesting that rates of recombination may be higher in humans. This baboon linkage map is the first reported for any nonhuman primate species and creates opportunities for mapping quantitative trait loci in baboons, as well as for comparative evolutionary analyses of genome structure.

Prior segregation analysis and the power to detect linkage

Complex parametric segregation and linkage analysis was performed on the simulated quantitative trait Q1 for all 200 replicates of the nuclear families data set. The segregation analysis inferred a major gene in 46% of the replicates. Among all replicates, including those that rejected a major gene, the power to detect suggestive linkage to any of three loci affecting Q1 was 0.600 and the false positive rate was 0.002. Among the replicates where a major gene was found, the power to detect suggestive linkage was 0.652 and the false positive rate was also 0.002. Thus, for purposes of linkage to this complex trait, a prior segregation then linkage analysis approach located a gene in 30% of all replicates, whereas a linkage only approach located a gene in 60% of all replicates.

Identifying influential individuals in linkage analysis: application to a quantitative trait locus detected in the COGA data

Once linkage is detected to a quantitative trait locus (QTL), the next step towards localizing the gene involved may be to identify those families, or individuals, in whom the putative mutations are segregating. In this paper, we describe a jackknife procedure for identifying individuals (and families) who contribute disproportionately to the linkage. Following initial detection of linkage to a QTL, the strategy involves sequentially removing each individual (or each family) from the analysis and recomputing the lod score associated with the linked region using data from all remaining subjects (or families). This procedure can be used to determine if particular observations have substantial impact on evidence for linkage. Identification of such observations may provide insights for further efforts to localize the QTL.

Two major loci control variation in beta-lipoprotein cholesterol and response to dietary fat and cholesterol in baboons

We explored the genetic control of cholesterolemic responses to dietary cholesterol and fat in 575 pedigreed baboons. We measured cholesterol in beta-lipoproteins (low density lipoprotein cholesterol [LDLC]) in blood drawn from baboons while they were consuming a baseline (low in cholesterol and fat) diet, a high-saturated fat (lard) diet, and a high-cholesterol, high-saturated fat diet. In addition to baseline levels (LDLC(Base)), we analyzed two variables for diet response: LDLC(RF), which represents the LDLC response to increasing dietary fat (ie, high-fat diet minus baseline), and LDLC(RC), which represents the LDLC response to increasing dietary cholesterol level (ie, high-cholesterol, high-fat diet minus high-fat diet). Heritabilities (h2) of the 3 traits were 0.59 for LDLC(Base), 0.14 for LDLC(RF), and 0.59 for LDLC(RC). In addition, LDLC(Base) and LDLC(RC) had a significant genetic correlation (ie, rhoG=0.54), suggesting that 1 or more genes exert pleiotropic effects on the 2 traits. Segregation analyses detected a single major locus that accounted for nearly all genetic variation in LDLC(RC) and some genetic variation in LDLC(Base) and LDLC(RF) and confirmed the presence of a different major locus that influences LDLC(Base) alone. Preliminary linkage analyses indicated that neither locus was linked to the LDL receptor gene, a likely candidate locus for LDLC. Detection of these major loci with large effects on the LDLC response to dietary cholesterol in a nonhuman primate offers hope of detecting and ultimately identifying similar loci that determine LDLC variation in human populations.

Characterization of a composite gradient gel for the electrophoretic separation of lipoproteins

We describe a protocol for making a new type of gradient gel, the Composite gradient gel, that was designed to resolve plasma lipoproteins using nondenaturing gradient gel electrophoresis. The new gel format allows analysis both of high density lipoproteins (HDLs) and low density lipoproteins (LDLs) on the same gel. The gel gave highly repeatable (r2 = 0.999) size estimates. We compared lipoprotein phenotypes determined from the new gradient gel with those obtained using specialized HDL and LDL gradient gels. The comparisons indicated that the Composite gel gave lipoprotein particle size estimates for HDLs and LDLs that were virtually identical to those obtained, respectively, from the specialized HDL and LDL gradient gels. We measured median diameters, which reflect the distributions of absorbance, for LDLs and for HDLs and found that the Composite gel gave lipoprotein size distributions that were virtually identical to those measured using the specialized LDL and HDL gels. Finally, comparison of fractional absorbance for six lipoprotein size intervals obtained from the Composite and specialized gels revealed a close correlation (r2 = 0.828). Thus, it appears that both LDL and HDL size phenotypes may be evaluated simultaneously using a single gradient gel format.