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.

DNA Methylation Variation Is Identified in Monozygotic Twins Discordant for Non-syndromic Cleft Lip and Palate.

Non-syndromic cleft lip with or without cleft palate (NSCLP) is the most common craniofacial birth defect. The etiology of NSCLP is complex with multiple genes and environmental factors playing causal roles. Although studies have identified numerous genetic markers associated with NSCLP, the role of epigenetic variation remains relatively unexplored. Because of their identical DNA sequences, monozygotic (MZ) twins discordant for NSCLP are an ideal model for examining the potential contribution of DNA methylation to non-syndromic orofacial clefting. In this study, we compared the patterns of whole genome DNA methylation in six MZ twin pairs discordant for NSCLP. Differentially methylated positions (DMPs) and regions (DMRs) were identified in NSCLP candidate genes, including differential methylation in MAFB and ZEB2 in two independent MZ twin pairs. In addition to DNA methylation differences in NSCLP candidate genes, we found common differential methylation in genes belonging to the Hippo signaling pathway, implicating this mechanosensory pathway in the etiology of NSCLP. The results of this novel approach using MZ twins discordant for NSCLP suggests that differential methylation is one mechanism contributing to NSCLP, meriting future studies on the role of DNA methylation in familial and sporadic NSCLP.

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.

Early-Onset Alzheimer Disease and Candidate Risk Genes Involved in Endolysosomal Transport.

Importance: Mutations in APP, PSEN1, and PSEN2 lead to early-onset Alzheimer disease (EOAD) but account for only approximately 11% of EOAD overall, leaving most of the genetic risk for the most severe form of Alzheimer disease unexplained. This extreme phenotype likely harbors highly penetrant risk variants, making it primed for discovery of novel risk genes and pathways for AD. Objective: To search for rare variants contributing to the risk for EOAD. Design, Setting, and Participants: In this case-control study, whole-exome sequencing (WES) was performed in 51 non-Hispanic white (NHW) patients with EOAD (age at onset <65 years) and 19 Caribbean Hispanic families previously screened as negative for established APP, PSEN1, and PSEN2 causal variants. Participants were recruited from John P. Hussman Institute for Human Genomics, Case Western Reserve University, and Columbia University. Rare, deleterious, nonsynonymous, or loss-of-function variants were filtered to identify variants in known and suspected AD genes, variants in multiple unrelated NHW patients, variants present in 19 Hispanic EOAD WES families, and genes with variants in multiple unrelated NHW patients. These variants/genes were tested for association in an independent cohort of 1524 patients with EOAD, 7046 patients with late-onset AD (LOAD), and 7001 cognitively intact controls (age at examination, >65 years) from the Alzheimer's Disease Genetics Consortium. The study was conducted from January 21, 2013, to October 13, 2016. Main Outcomes and Measures: Alzheimer disease diagnosed according to standard National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer Disease and Related Disorders Association criteria. Association between Alzheimer disease and genetic variants and genes was measured using logistic regression and sequence kernel association test-optimal gene tests, respectively. Results: Of the 1524 NHW patients with EOAD, 765 (50.2%) were women and mean (SD) age was 60.0 (4.9) years; of the 7046 NHW patients with LOAD, 4171 (59.2%) were women and mean (SD) age was 77.4 (8.6) years; and of the 7001 NHW controls, 4215 (60.2%) were women and mean (SD) age was 77.4 (8.6) years. The gene PSD2, for which multiple unrelated NHW cases had rare missense variants, was significantly associated with EOAD (P = 2.05 × 10-6; Bonferroni-corrected P value [BP] = 1.3 × 10-3) and LOAD (P = 6.22 × 10-6; BP = 4.1 × 10-3). A missense variant in TCIRG1, present in a NHW patient and segregating in 3 cases of a Hispanic family, was more frequent in EOAD cases (odds ratio [OR], 2.13; 95% CI, 0.99-4.55; P = .06; BP = 0.413), and significantly associated with LOAD (OR, 2.23; 95% CI, 1.37-3.62; P = 7.2 × 10-4; BP = 5.0 × 10-3). A missense variant in the LOAD risk gene RIN3 showed suggestive evidence of association with EOAD after Bonferroni correction (OR, 4.56; 95% CI, 1.26-16.48; P = .02, BP = 0.091). In addition, a missense variant in RUFY1 identified in 2 NHW EOAD cases showed suggestive evidence of an association with EOAD as well (OR, 18.63; 95% CI, 1.62-213.45; P = .003; BP = 0.129). Conclusions and Relevance: The genes PSD2, TCIRG1, RIN3, and RUFY1 all may be involved in endolysosomal transport-a process known to be important to development of AD. Furthermore, this study identified shared risk genes between EOAD and LOAD similar to previously reported genes, such as SORL1, PSEN2, and TREM2.

Genetic Characterization and Influence on Inflammatory Bowel Disease Expression in a Diverse Hispanic South Florida Cohort.

OBJECTIVES: Hispanics represent an understudied inflammatory bowel disease (IBD) population. Prior studies examining genetic predisposition to IBD in Hispanics are limited. In this study, we examined whether European-derived IBD variants confer risk in Hispanics and their influence on IBD phenotype in Hispanics compared to non-Hispanic whites (NHW). METHODS: Self-identified Hispanics and NHWs with IBD were included. Hispanic controls were included for our genetic analyses. We performed single-variant testing at previously identified Crohn's disease (CD) and ulcerative colitis (UC) IBD variants in Hispanic cases and controls. These risk variants were used to compute individual genetic risk scores. Genetic risk scores and phenotype associations were compared between Hispanic and NHW. RESULTS: A total of 1,115 participants were included: 698 controls and 417 IBD patients (230 Hispanics). We found evidence of association within our Hispanic cohort at 22 IBD risk loci, with ~76% of the risk loci demonstrating over-representation of the European risk allele; these included loci corresponding to IL23R and NOD2 genes. CD genetic risk score for Hispanics (199.67) was similar to the score for NHW (200.33), P=0.51; the same was true in UC. Genetic risk scores did not predict IBD phenotype or complications in Hispanics or NHW except for a younger age of CD onset in Hispanics (P=0.04). CONCLUSIONS: This study highlights the fundamental importance of these loci in IBD pathogenesis including in our diverse Hispanic population. Future studies looking at non-genetic mechanisms of disease are needed to explain differences in age of presentation and phenotype between Hispanics and NHW.

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 follow-up study for left ventricular mass on chromosome 12p11 identifies potential candidate genes.

BACKGROUND: Left ventricular mass (LVM) is an important risk factor for cardiovascular disease. Previously we found evidence for linkage to chromosome 12p11 in Dominican families, with a significant increase in a subset of families with high average waist circumference (WC). In the present study, we use association analysis to further study the genetic effect on LVM. METHODS: Association analysis with LVM was done in the one LOD critical region of the linkage peak in an independent sample of 897 Caribbean Hispanics. Genotype data were available on 7085 SNPs from 23 to 53 MB on chromosome 12p11. Adjustment was made for vascular risk factors and population substructure using an additive genetic model. Subset analysis by WC was performed to test for a difference in genetic effects between the high and low WC subsets. RESULTS: In the overall analysis, the most significant association was found to rs10743465, downstream of the SOX5 gene (p = 1.27E-05). Also, 19 additional SNPs had nominal p < 0.001. In the subset analysis, the most significant difference in genetic effect between those with high and low WC occurred with rs1157480 (p = 1.37E-04 for the difference in ß coefficients), located upstream of TMTC1. Twelve additional SNPs in or near 6 genes had p < 0.001. CONCLUSIONS: The current study supports previously identified evidence by linkage for a genetic effect on LVM on chromosome 12p11 using association analysis in population-based Caribbean Hispanic cohort. SOX5 may play an important role in the regulation of LVM. An interaction of TMTC1 with abdominal obesity may contribute to phenotypic variation of LVM.

A dominant mutation in RPE65 identified by whole-exome sequencing causes retinitis pigmentosa with choroidal involvement.

Linkage testing using Affymetrix 6.0 SNP Arrays mapped the disease locus in TCD-G, an Irish family with autosomal dominant retinitis pigmentosa (adRP), to an 8.8 Mb region on 1p31. Of 50 known genes in the region, 11 candidates, including RPE65 and PDE4B, were sequenced using di-deoxy capillary electrophoresis. Simultaneously, a subset of family members was analyzed using Agilent SureSelect All Exome capture, followed by sequencing on an Illumina GAIIx platform. Candidate gene and exome sequencing resulted in the identification of an Asp477Gly mutation in exon 13 of the RPE65 gene tracking with the disease in TCD-G. All coding exons of genes not sequenced to sufficient depth by next generation sequencing were sequenced by di-deoxy sequencing. No other potential disease-causing variants were found to segregate with disease in TCD-G. The Asp477Gly mutation was not present in Irish controls, but was found in a second Irish family provisionally diagnosed with choroideremia, bringing the combined maximum two-point LOD score to 5.3. Mutations in RPE65 are a known cause of recessive Leber congenital amaurosis (LCA) and recessive RP, but no dominant mutations have been reported. Protein modeling suggests that the Asp477Gly mutation may destabilize protein folding, and mutant RPE65 protein migrates marginally faster on SDS-PAGE, compared with wild type. Gene therapy for LCA patients with RPE65 mutations has shown great promise, raising the possibility of related therapies for dominant-acting mutations in this gene.

Genetic loci for blood lipid levels identified by linkage and association analyses in Caribbean Hispanics.

To identify genetic loci influencing blood lipid levels in Caribbean Hispanics, we first conducted a genome-wide linkage scan in 1,211 subjects from 100 Dominican families on five lipid quantitative traits: total cholesterol (TC), low density lipoprotein cholesterol (LDL-C), high density lipoprotein cholesterol (HDL-C), triglycerides (TG), and LDL-C/HDL-C ratio. We then investigated the association between blood lipid levels and 21,361 single nucleotide polymorphisms (SNP) under the 1-logarithm of odds (LOD) unit down regions of linkage peaks in an independent community-based subcohort (N = 814, 42% Dominican) from the Northern Manhattan Study (NOMAS). We found significant linkage evidence for LDL-C/HDL-C on 7p12 (multipoint LOD = 3.91) and for TC on 16q23 (LOD = 3.35). In addition, we identified suggestive linkage evidence of LOD > 2.0 on 15q23 for TG, 16q23 for LDL-C, 19q12 for TC and LDL-C, and 20p12 for LDL-C. In the association analysis of the linkage peaks, we found that seven SNPs near FLJ45974 were associated with LDL-C/HDL-C with a nominal P < 3.5 × 10(-5), in addition to associations (P < 0.0001) for other lipid traits with SNPs in or near CDH13, SUMF2, TLE3, FAH, ARNT2, TSHZ3, ZNF343, RPL7AL2, and TMC3. Further studies are warranted to perform in-depth investigations of functional genetic variants in these regions.

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.).

Carotid plaque and candidate genes related to inflammation and endothelial function in Hispanics from northern Manhattan.

BACKGROUND AND PURPOSE: The genetic influence on carotid atherosclerotic plaque is mostly unknown. This study examines the association between carotid plaque and single nucleotide polymorphisms in selected genes implicated in inflammation and endothelial function. METHODS: A total of 43 genes (197 single nucleotide polymorphisms) involved in inflammation and endothelial function were interrogated in 287 Dominicans from the Northern Manhattan Study (mean age, 64±7 years; 58% women) who had undergone high-resolution B-mode ultrasound for examination of carotid plaque. Using an additive genetic model, multiple logistic regression analyses were conducted, a within-gene haplotype analysis was performed, and interactions between genes were examined. Results were validated in an independent set of 301 Dominicans. RESULTS: Carotid plaque was present in 143 (47%) participants. Nine genes had at least 1 single nucleotide polymorphism associated (P≤0.01) with carotid plaque phenotypes: TNF, NOS2A, IL6R, TNFSF4, PPARA, IL1A, TLR4, ITGA2, and HABP2. Single nucleotide polymorphisms in TNFSF4, PPARA, TLR4, ITGA2, and HABP2 were also implicated with the same carotid phenotype in the validation analysis. Haplotype analysis revealed an additional gene of interest, VCAM1. CONCLUSIONS: We report novel associations between variations in 10 genes involved in inflammation and endothelial function and carotid plaque phenotypes in a Dominican sample, with replication for 5 genes in an independent Dominican sample.