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.

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 (ORq2 vs.q1 = 1.94, 95% CI, 0.86-4.39; ORq3 vs.q1 = 3.46, 95% CI, 1.57-7.63; ORq4 vs.q1 = 8.64, 95% CI, 3.92-19.02; and Ptrend < .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.

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.

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.

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.

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.