A Common Variant in MIR182 Is Associated With Primary Open-Angle Glaucoma in the NEIGHBORHOOD Consortium

Purpose

Noncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG).

Methods

Using the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR-182 expression in AH between five HTG cases and five controls.

Results

Only rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] = 1.23, 95% confidence interval [CI]: 1.11–1.42, P = 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR = 1.26, 95% CI: 1.08–1.47, P = 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P = 0.03) without controlling for medication treatment.

Conclusions

Our integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

Evaluating genetic susceptibility to Staphylococcus aureus bacteremia in African Americans using admixture mapping

The incidence of Staphylococcus aureus bacteremia (SAB) is significantly higher in African American (AA) than in European-descended populations. We used admixture mapping (AM) to test the hypothesis that genomic variations with different frequencies in European and African ancestral genomes influence susceptibility to SAB in AAs. A total of 565 adult AAs (390 cases with SAB; 175 age-matched controls) were genotyped for AM analysis. A case-only admixture score and a mixed χ²(1df) score (MIX) to jointly evaluate both single-nucleotide polymorphism (SNP) and admixture association (P<5.00e-08) were computed using MIXSCORE. In addition, a permutation scheme was implemented to derive multiplicity adjusted P-values (genome-wide 0.05 significance threshold: P<9.46e-05). After empirical multiplicity adjustment, one region on chromosome 6 (52 SNPs, P=4.56e-05) in the HLA class II region was found to exhibit a genome-wide statistically significant increase in European ancestry. This region encodes genes involved in HLA-mediated immune response and these results provide additional evidence for genetic variation influencing HLA-mediated immunity, modulating susceptibility to SAB.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Generation of disease-specific autopsy-confirmed iPSCs lines from postmortem isolated Peripheral Blood Mononuclear Cells

Understanding the molecular mechanisms that underlie neurodegenerative disorders has been hampered by a lack of readily available model systems that replicate the complexity of the human disease. Recent advances in stem cell technology have facilitated the derivation of patient-specific stem cells from a variety of differentiated cell types. These induced pluripotent stem cells (iPSCs) are attractive disease models since they can be grown and differentiated to produce large numbers of disease-relevant cell types. However, most iPSC lines are derived in advance of, and without the benefit of, neuropathological confirmation of the donor - the gold standard for many disease classifications and measurement of disease severity. While others have reported the generation of autopsy-confirmed iPSC lines from patient explants, these methods require outgrowth of cadaver tissue, which require additional time and is often only successful ∼50% of the time. Here we report the rapid generation of autopsy-confirmed iPSC lines from peripheral blood mononuclear cells (PBMCs) drawn postmortem. Since this approach doesn't require the propagation of previously frozen cadaver tissue, iPSC can be rapidly and efficiently produced from patients with autopsy-confirmed pathology. These matched iPSC-derived patient-specific neurons and postmortem brain tissue will support studies of specific mechanisms that drive the pathogenesis of neurodegenerative diseases.

Heritability of Choroidal Thickness in the Amish

PURPOSE

To evaluate the heritability of choroidal thickness and its relationship to age-related macular degeneration (AMD).

DESIGN

Cohort study.

PARTICIPANTS

Six hundred eighty-nine individuals from Amish families with early or intermediate AMD.

METHODS

Ocular coherence tomography was used to quantify choroidal thickness, and fundus photography was used to classify eyes into categories using a modified Clinical Age-Related Maculopathy Staging (CARMS) system. Repeatability and heritability of choroidal thickness and its phenotypic and genetic correlations with the AMD phenotype (CARMS category) were estimated using a generalized linear mixed model (GLMM) approach that accounted for relatedness, repeated measures (left and right eyes), and the effects of age, gender, and refraction.

MAIN OUTCOME MEASURES

Heritability of choroidal thickness and its phenotypic and genetic correlation with the AMD phenotype (CARMS category).

RESULTS

Phenotypic correlation between choroidal thickness and CARMS category was moderate (Spearman's rank correlation, r_s = -0.24; n = 1313 eyes) and significant (GLMM posterior mean, -4.27; 95% credible interval [CI], -7.88 to -0.79; P = 0.02) after controlling for relatedness, age, gender, and refraction. Eyes with advanced AMD had thinner choroids than eyes without AMD (posterior mean, -73.8; 95% CI, -94.7 to -54.6; P < 0.001; n = 1178 eyes). Choroidal thickness was highly repeatable within individuals (repeatability, 0.78; 95% CI, 0.68 to 0.89) and moderately heritable (heritability, 0.40; 95% CI, 0.14 to 0.51), but did not show significant genetic correlation with CARMS category, although the effect size was moderate (genetic correlation, -0.18; 95% CI, -0.49 to 0.16). Choroidal thickness also varied with age, gender, and refraction. The CARMS category showed moderate heritability (heritability, 0.49; 95% CI, 0.26 to 0.72).

CONCLUSIONS

We quantify the heritability of choroidal thickness for the first time, highlighting a heritable, quantitative trait that is measurable in all individuals regardless of AMD affection status, and moderately phenotypically correlated with AMD severity. Choroidal thickness therefore may capture variation not captured by the CARMS system. However, because the genetic correlation between choroidal thickness and AMD severity was not significant in our data set, genes associated with the 2 traits may not overlap substantially. Future studies should therefore test for genetic variation associated with choroidal thickness to determine the overlap in genetic basis with AMD.

Assessing the Association of Mitochondrial Genetic Variation With Primary Open-Angle Glaucoma Using Gene-Set Analyses

Purpose

Recent studies indicate that mitochondrial proteins may contribute to the pathogenesis of primary open-angle glaucoma (POAG). In this study, we examined the association between POAG and common variations in gene-encoding mitochondrial proteins.

Methods

We examined genetic data from 3430 POAG cases and 3108 controls derived from the combination of the GLAUGEN and NEIGHBOR studies. We constructed biological-system coherent mitochondrial nuclear-encoded protein gene-sets by intersecting the MitoCarta database with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We examined the mitochondrial gene-sets for association with POAG and with normal-tension glaucoma (NTG) and high-tension glaucoma (HTG) subsets using Pathway Analysis by Randomization Incorporating Structure.

Results

We identified 22 KEGG pathways with significant mitochondrial protein-encoding gene enrichment, belonging to six general biological classes. Among the pathway classes, mitochondrial lipid metabolism was associated with POAG overall (P = 0.013) and with NTG (P = 0.0006), and mitochondrial carbohydrate metabolism was associated with NTG (P = 0.030). Examining the individual KEGG pathway mitochondrial gene-sets, fatty acid elongation and synthesis and degradation of ketone bodies, both lipid metabolism pathways, were significantly associated with POAG (P = 0.005 and P = 0.002, respectively) and NTG (P = 0.0004 and P < 0.0001, respectively). Butanoate metabolism, a carbohydrate metabolism pathway, was significantly associated with POAG (P = 0.004), NTG (P = 0.001), and HTG (P = 0.010).

Conclusions

We present an effective approach for assessing the contributions of mitochondrial genetic variation to open-angle glaucoma. Our findings support a role for mitochondria in POAG pathogenesis and specifically point to lipid and carbohydrate metabolism pathways as being important.

Preparedness of Educational Speech-Language Pathologists to Provide Services to Students with Traumatic Brain Injury

This study used a questionnaire to survey South Carolina's educational speech-language pathologists (SLPs)to determine their perceived knowledge about TBI and their preparedness to provide services to students with traumatic brain injury (TBI). Of the 768 questionnaires distributed, 220 (29%)were returned from 43 of 46 counties. Respondents had received degrees from 72 university and college training programs in 27 states. Results indicated that the majority of respondents did not express a high level of comfort regarding preparation for pediatric TBI rehabilitation, but SLPs with graduate training or national certification were more likely to feel prepared than other SLPs. The need for a stronger knowledge base in pediatric TBI was reflected in the mean number of correct responses (15.19)on a 20 item true-false test of general TBI knowledge. The survey found that interaction among respondents, referring rehabilitation teams, and school multidisciplinary teams was not high. The respondents indicated a desire for more education about TBI through direct inservice training. A list of references useful for preparing SLPs to provide services for students with TBI is included.

Whole exome sequencing of extreme age-related macular degeneration phenotypes

PURPOSE

Demographic, environmental, and genetic risk factors for age-related macular degeneration (AMD) have been identified; however, a substantial portion of the variance in AMD disease risk and heritability remains unexplained. To identify AMD risk variants and generate hypotheses for future studies, we performed whole exome sequencing for 75 individuals whose phenotype was not well predicted by their genotype at known risk loci. We hypothesized that these phenotypically extreme individuals were more likely to carry rare risk or protective variants with large effect sizes.

METHODS

A genetic risk score was calculated in a case-control set of 864 individuals (467 AMD cases, 397 controls) based on 19 common (≥1% minor allele frequency, MAF) single nucleotide variants previously associated with the risk of advanced AMD in a large meta-analysis of advanced cases and controls. We then selected for sequencing 39 cases with bilateral choroidal neovascularization with the lowest genetic risk scores to detect risk variants and 36 unaffected controls with the highest genetic risk score to detect protective variants. After minimizing the influence of 19 common genetic risk loci on case-control status, we targeted single variants of large effect and the aggregate effect of weaker variants within genes and pathways. Single variant tests were conducted on all variants, while gene-based and pathway analyses were conducted on three subsets of data: 1) rare (≤1% MAF in the European population) stop, splice, or damaging missense variants, 2) all rare variants, and 3) all variants. All analyses controlled for the effects of age and sex.

RESULTS

No variant, gene, or pathway outside regions known to be associated with risk for advanced AMD reached genome-wide significance. However, we identified several variants with substantial differences in allele frequency between cases and controls with strong additive effects on affection status after controlling for age and sex. Protective effects trending toward significance were detected at two loci identified in single-variant analyses: an intronic variant in FBLN7 (the gene encoding fibulin 7) and at three variants near pyridoxal (pyridoxine, vitamin B6) kinase (PDXK). Aggregate rare-variant analyses suggested evidence for association at ASRGL1, a gene previously linked to photoreceptor cell death, and at BSDC1. In known AMD loci we also identified 29 novel or rare damaging missense or stop/splice variants in our sample of cases and controls.

CONCLUSIONS

Identified variants and genes may highlight regions important in the pathogenesis of AMD and are key targets for replication.

DNA variants in CACNA1C modify Parkinson disease risk only when vitamin D level is deficient

Objective:

To evaluate the association between the genetic variants in CACNA1C, which encodes the α1 subunit of the L-type voltage-sensitive calcium channel (LVSCC) and Parkinson disease (PD) while accounting for interactions with vitamin D concentration.

Methods:

Two independent case-control data sets (478 cases and 431 controls; 482 cases and 412 controls) were used. Joint effects of single nucleotide polymorphisms (SNPs) and SNP-vitamin D interaction were analyzed by comparing models containing vitamin D deficiency, SNP genotypes, SNP-vitamin D interaction, and covariates to a restricted model with only vitamin D deficiency and covariates. Meta-analysis was used to combine the joint effects in the 2 data sets. Analysis was stratified by vitamin D deficiency to demonstrate the pattern of SNP-vitamin D interaction.

Results:

Vitamin D deficiency was associated with PD in both data sets (odds ratio [OR] = 1.9–2.7, p ≤ 0.009). SNP rs34621387 demonstrated a significant joint effect (meta-analysis, p = 7.5 × 10⁻⁵; Bonferroni corrected, p = 0.02). The G allele at rs34621387 is associated with PD in vitamin D-deficient individuals in both data sets (OR = 2.0–2.1, confidence interval = 1.3–3.5, p = 0.002) but is not associated with PD in vitamin D–nondeficient individuals (p > 0.8 in both data sets).

Conclusions:

Previous studies suggest that vitamin D deficiency is associated with PD and sustained opening of LVSCC contributes to the selective vulnerability of dopaminergic neurons in PD. Our data demonstrate that the association between genetic variations in CACNA1C and PD depends on vitamin D deficiency, providing one potential mechanism underlying the association between vitamin D deficiency and PD.

Genetic Association Analysis of Drusen Progression

PURPOSE

Age-related macular degeneration is a common form of vision loss affecting older adults. The etiology of AMD is multifactorial and is influenced by environmental and genetic risk factors. In this study, we examine how 19 common risk variants contribute to drusen progression, a hallmark of AMD pathogenesis.

METHODS

Exome chip data was made available through the International AMD Genomics Consortium (IAMDGC). Drusen quantification was carried out with color fundus photographs using an automated drusen detection and quantification algorithm. A genetic risk score (GRS) was calculated per subject by summing risk allele counts at 19 common genetic risk variants weighted by their respective effect sizes. Pathway analysis of drusen progression was carried out with the software package Pathway Analysis by Randomization Incorporating Structure.

RESULTS

We observed significant correlation with drusen baseline area and the GRS in the age-related eye disease study (AREDS) dataset (ρ = 0.175, P = 0.006). Measures of association were not statistically significant between drusen progression and the GRS (P = 0.54). Pathway analysis revealed the cell adhesion molecules pathway as the most highly significant pathway associated with drusen progression (corrected P = 0.02).

CONCLUSIONS

In this study, we explored the potential influence of known common AMD genetic risk factors on drusen progression. Our results from the GRS analysis showed association of increasing genetic burden (from 19 AMD associated loci) to baseline drusen load but not drusen progression in the AREDS dataset while pathway analysis suggests additional genetic contributors to AMD risk.

The Application of Genetic Risk Scores in Age-Related Macular Degeneration: A Review

Age-related macular degeneration (AMD), a highly prevalent and impactful disease of aging, is inarguably influenced by complex interactions between genetic and environmental factors. Various risk scores have been tested that assess measurable genetic and environmental contributions to disease. We herein summarize and review the ability and utility of these numerous models for prediction of AMD and suggest additional risk factors to be incorporated into clinically useful predictive models of AMD.

A Locus at 5q33.3 Confers Resistance to Tuberculosis in Highly Susceptible Individuals

Immunosuppression resulting from HIV infection increases the risk of progression to active tuberculosis (TB) both in individuals newly exposed to Mycobacterium tuberculosis (MTB) and in those with latent infections. We hypothesized that HIV-positive individuals who do not develop TB, despite living in areas where it is hyperendemic, provide a model of natural resistance. We performed a genome-wide association study of TB resistance by using 581 HIV-positive Ugandans and Tanzanians enrolled in prospective cohort studies of TB; 267 of these individuals developed active TB, and 314 did not. A common variant, rs4921437 at 5q33.3, was significantly associated with TB (odds ratio = 0.37, p = 2.11 × 10(-8)). This variant lies within a genomic region that includes IL12B and is embedded in an H3K27Ac histone mark. The locus also displays consistent patterns of linkage disequilibrium across African populations and has signals of strong selection in populations from equatorial Africa. Along with prior studies demonstrating that therapy with IL-12 (the cytokine encoded in part by IL12B, associated with longer survival following MTB infection in mice deficient in CD4 T cells), our results suggest that this pathway might be an excellent target for the development of new modalities for treating TB, especially for HIV-positive individuals. Our results also indicate that studying extreme disease resistance in the face of extensive exposure can increase the power to detect associations in complex infectious disease.

Polymorphisms in HLA Class II Genes Are Associated With Susceptibility to Staphylococcus aureus Infection in a White Population

BACKGROUND

Staphylococcus aureus can cause life-threatening infections. Human susceptibility to S. aureus infection may be influenced by host genetic variation.

METHODS

A genome-wide association study (GWAS) in a large health plan-based cohort included biologic specimens from 4701 culture-confirmed S. aureus cases and 45 344 matched controls; 584 535 single-nucleotide polymorphisms (SNPs) were genotyped on an array specific to individuals of European ancestry. Coverage was increased by imputation of >25 million common SNPs, using the 1000 Genomes Reference panel. In addition, human leukocyte antigen (HLA) serotypes were also imputed.

RESULTS

Logistic regression analysis, performed under the assumption of an additive genetic model, revealed several imputed SNPs (eg, rs115231074: odds ratio [OR], 1.22 [P = 1.3 × 10(-10)]; rs35079132: OR, 1.24 [P = 3.8 × 10(-8)]) achieving genome-wide significance on chromosome 6 in the HLA class II region. One adjacent genotyped SNP was nearly genome-wide significant (rs4321864: OR, 1.13; P = 8.8 × 10(-8)). These polymorphisms are located near the genes encoding HLA-DRA and HLA-DRB1. Results of further logistic regression analysis, in which the most significant GWAS SNPs were conditioned on HLA-DRB1*04 serotype, showed additional support for the strength of association between HLA class II genetic variants and S. aureus infection.

CONCLUSIONS

Our study results are the first reported evidence of human genetic susceptibility to S. aureus infection.

Linkage of familial essential tremor to chromosome 5q35

BACKGROUND

Essential tremor is a neurological condition characterized by tremor during voluntary movement. To date, 3 loci linked to familial essential tremor have been identified.

METHODS

We examined 48 essential tremor patients in 5 large essential tremor pedigrees in our data set for genetic linkage using an Affymetrix Axiom array. Linkage analysis was performed using an affecteds-only dominant model in SIMWALK2. To incorporate all genotype information, GERMLINE was used to identify genome segments shared identical-by-descent in pairs of affecteds. Exome sequencing was performed in pedigrees showing evidence of linkage.

RESULTS

For one family, chromosomes 5 and 18 showed genome-wide significant linkage to essential tremor. Shared segment analysis excluded the 18p11 candidate region and reduced the 5q35 region by 1 megabase. Exome sequencing did not identify a potential causative variant in this region.

CONCLUSION

A locus on chromosome 5 is linked to essential tremor. Further research is needed to identify a causative variant. © 2016 International Parkinson and Movement Disorder Society.

Overlap between Parkinson disease and Alzheimer disease in ABCA7 functional variants

OBJECTIVE

Given their reported function in phagocytosis and clearance of protein aggregates in Alzheimer disease (AD), we hypothesized that variants in ATP-binding cassette transporter A7 (ABCA7) might be involved in Parkinson disease (PD).

METHODS

ABCA7 variants were identified using whole-exome sequencing (WES) on 396 unrelated patients with PD and 222 healthy controls. In addition, we used the publicly available WES data from the Parkinson's Progression Markers Initiative (444 patients and 153 healthy controls) as a second, independent data set.

RESULTS

We observed a higher frequency of loss-of-function (LOF) variants and rare putative highly functional variants (Combined Annotation Dependent Depletion [CADD] >20) in clinically diagnosed patients with PD than in healthy controls in both data sets. Overall, we identified LOF variants in 11 patients and 1 healthy control (odds ratio [OR] 4.94, Fisher exact p = 0.07). Four of these variants have been previously implicated in AD risk (p.E709AfsX86, p.W1214X, p.L1403RfsX7, and rs113809142). In addition, rare variants with CADD >20 were observed in 19 patients vs 3 healthy controls (OR 2.85, Fisher exact p = 0.06).

CONCLUSION

The presence of ABCA7 LOF variants in clinically defined PD suggests that they might be risk factors for neurodegeneration in general, especially those variants hallmarked by protein aggregation. More studies will be needed to evaluate the overall impact of this transporter in neurodegenerative disease.

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

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

Estimating cumulative pathway effects on risk for age-related macular degeneration using mixed linear models

Background

Age-related macular degeneration (AMD) is the leading cause of irreversible visual loss in the elderly in developed countries and typically affects more than 10 % of individuals over age 80. AMD has a large genetic component, with heritability estimated to be between 45 % and 70 %. Numerous variants have been identified and implicate various molecular mechanisms and pathways for AMD pathogenesis but those variants only explain a portion of AMD’s heritability. The goal of our study was to estimate the cumulative genetic contribution of common variants on AMD risk for multiple pathways related to the etiology of AMD, including angiogenesis, antioxidant activity, apoptotic signaling, complement activation, inflammatory response, response to nicotine, oxidative phosphorylation, and the tricarboxylic acid cycle. While these mechanisms have been associated with AMD in literature, the overall extent of the contribution to AMD risk for each is unknown.

Methods

In a case–control dataset with 1,813 individuals genotyped for over 600,000 SNPs we used Genome-wide Complex Trait Analysis (GCTA) to estimate the proportion of AMD risk explained by SNPs in genes associated with each pathway. SNPs within a 50 kb region flanking each gene were also assessed, as well as more distant, putatively regulatory SNPs, based on DNaseI hypersensitivity data from ocular tissue in the ENCODE project.

Results

We found that 19 previously associated AMD risk SNPs contributed to 13.3 % of the risk for AMD in our dataset, while the remaining genotyped SNPs contributed to 36.7 % of AMD risk. Adjusting for the 19 risk SNPs, the complement activation and inflammatory response pathways still explained a statistically significant proportion of additional risk for AMD (9.8 % and 17.9 %, respectively), with other pathways showing no significant effects (0.3 % – 4.4 %).

Discussion

Our results show that SNPs associated with complement activation and inflammation significantly contribute to AMD risk, separately from the risk explained by the 19 known risk SNPs. We found that SNPs within 50 kb regions flanking genes explained additional risk beyond genic SNPs, suggesting a potential regulatory role, but that more distant SNPs explained less than 0.5 % additional risk for each pathway.

Conclusions

From these analyses we find that the impact of complement SNPs on risk for AMD extends beyond the established genome-wide significant SNPs.

Electronic supplementary material

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

PARK10 is a major locus for sporadic neuropathologically confirmed Parkinson disease

OBJECTIVE

To minimize pathologic heterogeneity in genetic studies of Parkinson disease (PD), the Autopsy-Confirmed Parkinson Disease Genetics Consortium conducted a genome-wide association study using both patients with neuropathologically confirmed PD and controls.

METHODS

Four hundred eighty-four cases and 1,145 controls met neuropathologic diagnostic criteria, were genotyped, and then imputed to 3,922,209 variants for genome-wide association study analysis.

RESULTS

A small region on chromosome 1 was strongly associated with PD (rs10788972; p = 6.2 × 10(-8)). The association peak lies within and very close to the maximum linkage peaks of 2 prior positive linkage studies defining the PARK10 locus. We demonstrate that rs10788972 is in strong linkage disequilibrium with rs914722, the single nucleotide polymorphism defining the PARK10 haplotype previously shown to be significantly associated with age at onset in PD. The region containing the PARK10 locus was significantly reduced from 10.6 megabases to 100 kilobases and contains 4 known genes: TCEANC2, TMEM59, miR-4781, and LDLRAD1.

CONCLUSIONS

We confirm the association of a PARK10 haplotype with the risk of developing idiopathic PD. Furthermore, we significantly reduce the size of the PARK10 region. None of the candidate genes in the new PARK10 region have been previously implicated in the biology of PD, suggesting new areas of potential research. This study strongly suggests that reducing pathologic heterogeneity may enhance the application of genetic association studies to PD.

Vitamin D from different sources is inversely associated with Parkinson disease

An inverse association between Parkinson disease (PD) and total vitamin D levels has been reported, but whether vitamin D from different sources, that is, 25(OH)D2 (from diet and supplements) and 25(OH)D3 (mainly from sunlight exposure), all contribute to the association is unknown. Plasma total 25(OH)D, 25(OH)D2, and 25(OH)D3 levels were measured by liquid chromatography-tandem mass spectrometry in PD patients (n = 478) and controls (n = 431). Total 25(OH)D was categorized by clinical insufficiency or deficiency; 25(OH)D2 and 25(OH)D3 were analyzed in quartiles. Vitamin D deficiency (total 25[OH]D < 20 ng/mL) and vitamin D insufficiency (total 25[OH]D < 30 ng/mL) are associated with PD risk (odds ratio [OR] = 2.6 [deficiency] and 2.1 [insufficiency]; P < 0.0001), adjusting for age, sex, and sampling season. Both 25(OH)D2 and 25(OH)D3 levels are inversely associated with PD (P(trend) < 0.0001). The association between 25(OH)D2 and PD risk is largely confined to individuals with low 25(OH)D3 levels (P(trend) = 0.0008 and 0.12 in individuals with 25[OH]D3 < 20 ng/mL and 25[OH]D3 ≥ 20 ng/mL, respectively). Our data confirm the association between vitamin D deficiency and PD, and for the first time demonstrate an inverse association of 25(OH)D2 with PD. Given that 25(OH)D2 concentration is independent of sunlight exposure, this new finding suggests that the inverse association between vitamin D levels and PD is not simply attributable to lack of sunlight exposure in PD patients with impaired mobility. The current study, however, cannot exclude the possibility that gastrointestinal dysfunction, a non-motor PD symptom, contributes to the lower vitamin D2 levels in PD patients.

DNA copy number variants of known glaucoma genes in relation to primary open-angle glaucoma

PURPOSE

We examined the role of DNA copy number variants (CNVs) of known glaucoma genes in relation to primary open angle glaucoma (POAG).

METHODS

Our study included DNA samples from two studies (NEIGHBOR and GLAUGEN). All the samples were genotyped with the Illumina Human660W_Quad_v1 BeadChip. After removing non-blood-derived and amplified DNA samples, we applied quality control steps based on the mean Log R Ratio and the mean B allele frequency. Subsequently, data from 3057 DNA samples (1599 cases and 1458 controls) were analyzed with PennCNV software. We defined CNVs as those ≥5 kilobases (kb) in size and interrogated by ≥5 consecutive probes. We further limited our investigation to CNVs in known POAG-related genes, including CDKN2B-AS1, TMCO1, SIX1/SIX6, CAV1/CAV2, the LRP12-ZFPM2 region, GAS7, ATOH7, FNDC3B, CYP1B1, MYOC, OPTN, WDR36, SRBD1, TBK1, and GALC.

RESULTS

Genomic duplications of CDKN2B-AS1 and TMCO1 were each found in a single case. Two cases carried duplications in the GAS7 region. Genomic deletions of SIX6 and ATOH7 were each identified in one case. One case carried a TBK1 deletion and another case carried a TBK1 duplication. No controls had duplications or deletions in these six genes. A single control had a duplication in the MYOC region. Deletions of GALC were observed in five cases and two controls.

CONCLUSIONS

The CNV analysis of a large set of cases and controls revealed the presence of rare CNVs in known POAG susceptibility genes. Our data suggest that these rare CNVs may contribute to POAG pathogenesis and merit functional evaluation.

The relationship between obsessive-compulsive symptoms and PARKIN genotype: The CORE-PD study

BACKGROUND

Few studies have systematically investigated the association between PARKIN genotype and psychiatric co-morbidities of Parkison's disease (PD). PARKIN-associated PD is characterized by severe nigral dopaminergic neuronal loss, a finding that may have implications for behaviors rooted in dopaminergic circuits such as obsessive-compulsive symptoms (OCS).

METHODS

The Schedule of Compulsions and Obsessions Patient Inventory (SCOPI) was administered to 104 patients with early-onset PD and 257 asymptomatic first-degree relatives. Carriers of one and two PARKIN mutations were compared with noncarriers.

RESULTS

Among patients, carriers scored lower than noncarriers in adjusted models (one-mutation: 13.9 point difference, P = 0.03; two-mutation: 24.1, P = 0.001), where lower scores indicate less OCS. Among asymptomatic relatives, a trend toward the opposite was seen: mutation carriers scored higher than noncarriers (one mutation, P = 0.05; two mutations, P = 0.13).

CONCLUSIONS

First, a significant association was found between PARKIN mutation status and obsessive-compulsive symptom level in both PD and asymptomatic patients, suggesting that OCS might represent an early non-motor dopamine-dependent feature. Second, irrespective of disease status, heterozygotes were significantly different from noncarriers, suggesting that PARKIN heterozygosity may contribute to phenotype. © 2014 International Parkinson and Movement Disorder Society.

Rare variant APOC3 R19X is associated with cardio-protective profiles in a diverse population-based survey as part of the Epidemiologic Architecture for Genes Linked to Environment Study

BACKGROUND

A founder mutation was recently discovered and described as conferring favorable lipid profiles and reduced subclinical atherosclerotic disease in a Pennsylvania Amish population. Preliminary data have suggested that this null mutation APOC3 R19X (rs76353203) is rare in the general population.

METHODS AND RESULTS

To better describe the frequency and lipid profile in the general population, we as part of the Population Architecture using Genomics and Epidemiology I Study and the Epidemiological Architecture for Genes Linked to Environment Study genotyped rs76353203 in 1113 Amish participants from Ohio and Indiana and 19 613 participants from the National Health and Nutrition Examination Surveys (NHANES III, 1999 to 2002, and 2007 to 2008). We found no carriers among the Ohio and Indiana Amish. Of the 19 613 NHANES participants, we identified 31 participants carrying the 19X allele, for an overall allele frequency of 0.08%. Among fasting adults, the 19X allele was associated with lower triglycerides (n=7603; β=-71.20; P=0.007) and higher high-density lipoprotein cholesterol (n=8891; β=15.65; P=0.0002) and, although not significant, lower low-density lipoprotein cholesterol (n=6502; β= -4.85; P=0.68) after adjustment for age, sex, and race/ethnicity. On average, 19X allele participants had approximately half the triglyceride levels (geometric means, 51.3 to 69.7 versus 134.6 to 141.3 mg/dL), >20% higher high-density lipoprotein cholesterol levels (geometric means, 56.8 to 74.4 versus 50.38 to 53.36 mg/dL), and lower low-density lipoprotein cholesterol levels (geometric means, 104.5 to 128.6 versus 116.1 to 125.7 mg/dL) compared with noncarrier participants.

CONCLUSIONS

These data demonstrate that APOC3 19X exists in the general US population in multiple racial/ethnic groups and is associated with cardio-protective lipid profiles.

Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson's disease

We conducted a meta-analysis of Parkinson's disease genome-wide association studies using a common set of 7,893,274 variants across 13,708 cases and 95,282 controls. Twenty-six loci were identified as having genome-wide significant association; these and 6 additional previously reported loci were then tested in an independent set of 5,353 cases and 5,551 controls. Of the 32 tested SNPs, 24 replicated, including 6 newly identified loci. Conditional analyses within loci showed that four loci, including GBA, GAK-DGKQ, SNCA and the HLA region, contain a secondary independent risk variant. In total, we identified and replicated 28 independent risk variants for Parkinson's disease across 24 loci. Although the effect of each individual locus was small, risk profile analysis showed substantial cumulative risk in a comparison of the highest and lowest quintiles of genetic risk (odds ratio (OR) = 3.31, 95% confidence interval (CI) = 2.55-4.30; P = 2 × 10(-16)). We also show six risk loci associated with proximal gene expression or DNA methylation.

Hypothesis-independent pathway analysis implicates GABA and acetyl-CoA metabolism in primary open-angle glaucoma and normal-pressure glaucoma

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p < 0.001. In addition, the human leukocyte antigen (HLA) gene family was significantly associated with POAG (p < 0.001). In the POAG subset with normal-pressure glaucoma (NPG), the butanoate metabolism pathway was also significantly associated (p < 0.001) as well as the MAPK and Hedgehog signaling pathways (hsa04010 and hsa04340), glycosaminoglycan biosynthesis-heparan sulfate pathway (hsa00534) and the phenylalanine, tyrosine and tryptophan biosynthesis pathway (hsa0400). The butanoate metabolism pathway overall, and specifically the aspects of the pathway that contribute to GABA and acetyl-CoA metabolism, was the only pathway significantly associated with both POAG and NPG. Collectively these results implicate GABA and acetyl-CoA metabolism in glaucoma pathogenesis, and suggest new potential therapeutic targets.

Rare complement factor H variant associated with age-related macular degeneration in the Amish

PURPOSE

Age-related macular degeneration is the leading cause of blindness among the adult population in the developed world. To further the understanding of this disease, we have studied the genetically isolated Amish population of Ohio and Indiana.

METHODS

Cumulative genetic risk scores were calculated using the 19 known allelic associations. Exome sequencing was performed in three members of a small Amish family with AMD who lacked the common risk alleles in complement factor H (CFH) and ARMS2/HTRA1. Follow-up genotyping and association analysis was performed in a cohort of 973 Amish individuals, including 95 with self-reported AMD.

RESULTS

The cumulative genetic risk score analysis generated a mean genetic risk score of 1.12 (95% confidence interval [CI]: 1.10, 1.13) in the Amish controls and 1.18 (95% CI: 1.13, 1.22) in the Amish cases. This mean difference in genetic risk scores is statistically significant (P = 0.0042). Exome sequencing identified a rare variant (P503A) in CFH. Association analysis in the remainder of the Amish sample revealed that the P503A variant is significantly associated with AMD (P = 9.27 × 10(-13)). Variant P503A was absent when evaluated in a cohort of 791 elderly non-Amish controls, and 1456 non-Amish cases.

CONCLUSIONS

Data from the cumulative genetic risk score analysis suggests that the variants reported by the AMDGene consortium account for a smaller genetic burden of disease in the Amish compared with the non-Amish Caucasian population. Using exome sequencing data, we identified a novel missense mutation that is shared among a densely affected nuclear Amish family and located in a gene that has been previously implicated in AMD risk.

Epiregulin (EREG) and human V-ATPase (TCIRG1): genetic variation, ethnicity and pulmonary tuberculosis susceptibility in Guinea-Bissau and The Gambia

We analyzed two West African samples (Guinea-Bissau: n=289 cases and 322 controls; The Gambia: n=240 cases and 248 controls) to evaluate single-nucleotide polymorphisms (SNPs) in Epiregulin (EREG) and V-ATPase (T-cell immune regulator 1 (TCIRG1)) using single and multilocus analyses to determine whether previously described associations with pulmonary tuberculosis (PTB) in Vietnamese and Italians would replicate in African populations. We did not detect any significant single locus or haplotype associations in either sample. We also performed exploratory pairwise interaction analyses using Visualization of Statistical Epistasis Networks (ViSEN), a novel method to detect only interactions among multiple variables, to elucidate possible interaction effects between SNPs and demographic factors. Although we found no strong evidence of marginal effects, there were several significant pairwise interactions that were identified in either the Guinea-Bissau or the Gambian samples, two of which replicated across populations. Our results indicate that the effects of EREG and TCIRG1 variants on PTB susceptibility, to the extent that they exist, are dependent on gene-gene interactions in West African populations as detected with ViSEN. In addition, epistatic effects are likely to be influenced by inter- and intra-population differences in genetic or environmental context and/or the mycobacterial lineages causing disease.

Dusp3 and Psme3 are associated with murine susceptibility to Staphylococcus aureus infection and human sepsis

Using A/J mice, which are susceptible to Staphylococcus aureus, we sought to identify genetic determinants of susceptibility to S. aureus, and evaluate their function with regard to S. aureus infection. One QTL region on chromosome 11 containing 422 genes was found to be significantly associated with susceptibility to S. aureus infection. Of these 422 genes, whole genome transcription profiling identified five genes (Dcaf7, Dusp3, Fam134c, Psme3, and Slc4a1) that were significantly differentially expressed in a) S. aureus -infected susceptible (A/J) vs. resistant (C57BL/6J) mice and b) humans with S. aureus blood stream infection vs. healthy subjects. Three of these genes (Dcaf7, Dusp3, and Psme3) were down-regulated in susceptible vs. resistant mice at both pre- and post-infection time points by qPCR. siRNA-mediated knockdown of Dusp3 and Psme3 induced significant increases of cytokine production in S. aureus-challenged RAW264.7 macrophages and bone marrow derived macrophages (BMDMs) through enhancing NF-κB signaling activity. Similar increases in cytokine production and NF-κB activity were also seen in BMDMs from CSS11 (C57BL/6J background with chromosome 11 from A/J), but not C57BL/6J. These findings suggest that Dusp3 and Psme3 contribute to S. aureus infection susceptibility in A/J mice and play a role in human S. aureus infection.

Role of autophagy genetic variants for the risk of Candida infections

Candida albicans can cause candidemia in neutropenic and critically ill patients and oropharyngeal candidiasis in human immunodeficiency virus (HIV)-positive patients with low CD4(+) counts. Because all patients at risk do not develop Candida infections, it is possible that a patient's genetic background might play a role in his or her susceptibility to infection. Autophagy mediates pathogen clearance and modulation of inflammation. Our aim was to assess the effect of genetic variations in the ATG16L1 and IRGM autophagy genes on the susceptibility of patients with candidemia and oropharyngeal candidiasis. We assessed genetic variations in the ATG16L1 and IRGM genes in a cohort of candidemia patients of both African and European origin. In addition, we evaluated the effect of these polymorphisms on the susceptibility to oropharyngeal candidiasis of an HIV-positive cohort from Tanzania. Functional studies have been performed to assess the effect of the ATG16L1 and IRGM genetic variants on both in vitro and in vivo cytokine production. The results indicate that ATG16L1 variants modulate production of tumor necrosis factor-alpha, but not other cytokines, while no effects were seen in the presence of IRGM polymorphisms. In addition, no significant associations between the single-nucleotide polymorphisms in the ATG16L1 and IRGM genetic variants and the incidence of candidemia or oropharyngeal candidiasis were identified. Despite moderate effects on the modulation of proinflammatory cytokine production, genetic variation in the autophagy genes ATG16L1 and IRGM has a minor impact on the susceptibility to both mucosal and systemic Candida infections.

Absence of C9ORF72 expanded or intermediate repeats in autopsy-confirmed Parkinson's disease

BACKGROUND

We have reported that intermediate repeat lengths of the C9ORF72 repeat are a risk factor for Parkinson's disease (PD) in a clinically diagnosed data set. Because 10% to 25% of clinically diagnosed PD have different diagnoses upon autopsy, we hypothesized that this may reflect phenotypic heterogeneity or concomitant pathology of other neurodegenerative disorders.

METHODS

We screened 488 autopsy-confirmed PD cases for expansion haplotype tag rs3849942T. In 196 identified haplotype carriers, the C9ORF72 repeat was genotyped using the repeat-primed polymerase chain reaction assay.

RESULTS

No larger (intermediate or expanded) repeats were found in these autopsy-confirmed PD samples. This absence of larger repeats is significantly different from the frequency in clinically diagnosed datasets (P = 0.002).

CONCLUSIONS

Our results suggest that expanded or intermediate C9ORF72 repeats in clinically diagnosed PD or parkinsonism might be an indication of heterogeneity in clinically diagnosed PD cases. Further studies are needed to elucidate the potential contribution of the C9ORF72 repeat to autopsy-confirmed PD.

A genome-wide association study of variants associated with acquisition of Staphylococcus aureus bacteremia in a healthcare setting

Background

Humans vary in their susceptibility to acquiring Staphylococcus aureus infection, and research suggests that there is a genetic basis for this variability. Several recent genome-wide association studies (GWAS) have identified variants that may affect susceptibility to infectious diseases, demonstrating the potential value of GWAS in this arena.

Methods

We conducted a GWAS to identify common variants associated with acquisition of S. aureus bacteremia (SAB) resulting from healthcare contact. We performed a logistic regression analysis to compare patients with healthcare contact who developed SAB (361 cases) to patients with healthcare contact in the same hospital who did not develop SAB (699 controls), testing 542,410 SNPs and adjusting for age (by decade), sex, and 6 significant principal components from our EIGENSTRAT analysis. Additionally, we evaluated the joint effect of the host and pathogen genomes in association with severity of SAB infection via logistic regression, including an interaction of host SNP with bacterial genotype, and adjusting for age (by decade), sex, the 6 significant principal components, and dialysis status. Bonferroni corrections were applied in both analyses to control for multiple comparisons.

Results

Ours is the first study that has attempted to evaluate the entire human genome for variants potentially involved in the acquisition or severity of SAB. Although this study identified no common variant of large effect size to have genome-wide significance for association with either the risk of acquiring SAB or severity of SAB, the variant (rs2043436) most significantly associated with severity of infection is located in a biologically plausible candidate gene (CDON, a member of the immunoglobulin family) and may warrant further study.

Conclusions

The genetic architecture underlying SAB is likely to be complex. Future investigations using larger samples, narrowed phenotypes, and advances in both genotyping and analytical methodologies will be important tools for identifying causative variants for this common and serious cause of healthcare-associated infection.

Cognitive and motor function in long-duration PARKIN-associated Parkinson disease

IMPORTANCE

Data on the long-term cognitive outcomes of patients with PARKIN-associated Parkinson disease (PD) are unknown but may be useful when counseling these patients.

OBJECTIVE

Among patients with early-onset PD of long duration, we assessed cognitive and motor performances, comparing homozygotes and compound heterozygotes who carry 2 PARKIN mutations with noncarriers.

DESIGN, SETTING, AND PARTICIPANTS

Cross-sectional study of 44 participants at 17 different movement disorder centers who were in the Consortium on Risk for Early-Onset PD study with a duration of PD greater than the median duration (>14 years): 4 homozygotes and 17 compound heterozygotes (hereafter referred to as carriers) and 23 noncarriers.

MAIN OUTCOMES AND MEASURES

Unified Parkinson Disease Rating Scale Part III (UPDRS-III) and Clinical Dementia Rating scores and neuropsychological performance. Linear regression models were applied to assess the association between PARKIN mutation status and cognitive domain scores and UPDRS-III scores. Models were adjusted for age, education, disease duration, language, and levodopa equivalent daily dose.

RESULTS

Carriers had an earlier age at onset of PD (P < .001) and were younger (P = .004) at time of examination than noncarriers. They performed better than noncarriers on the Mini-Mental State Examination (P = .010) and were more likely to receive lower scores on the Clinical Dementia Rating (P = .003). In multivariate analyses, carriers performed better than noncarriers on the UPDRS-III (P = .02) and on tests of attention (P = .03), memory (P = .03), and visuospatial (P = .02) cognitive domains.

CONCLUSIONS AND RELEVANCE

In cross-sectional analyses, carriers demonstrated better cognitive and motor performance than did noncarriers with long disease duration, suggesting slower disease progression. A longitudinal follow-up study is required to confirm these findings.

Autophagy is redundant for the host defense against systemic Candida albicans infections

Autophagy has been demonstrated to play an important role in the immunity against intracellular pathogens, but very little is known about its role in the host defense against fungal pathogens such as Candida albicans. Therefore, the role of autophagy for the host defense against C. albicans was assessed by complementary approaches using mice defective in autophagy, as well as immunological and genetic studies in humans. Although C. albicans induced LC3-II formation in macrophages, myeloid cell-specific ATG7(-/-) mice with defects in autophagy did not display an increased susceptibility to disseminated candidiasis. In in vitro experiments in human blood mononuclear cells, blocking autophagy modulated cytokine production induced by lipopolysaccharide, but not by C. albicans. Furthermore, autophagy modulation in human monocytes did not influence the phagocytosis and killing of C. albicans. Finally, 18 single-nucleotide polymorphisms in 13 autophagy genes were not associated with susceptibility to candidemia or clinical outcome of disease in a large cohort of patients, and there was no correlation between these genetic variants and cytokine production in either candidemia patients or healthy controls. Based on these complementary in vitro and in vivo studies, it can be concluded that autophagy is redundant for the host response against systemic infections with C. albicans.

Genetic factors in nonsmokers with age-related macular degeneration revealed through genome-wide gene-environment interaction analysis

Relatively little is known about the interaction between genes and environment in the complex etiology of age-related macular degeneration (AMD). This study aimed to identify novel factors associated with AMD by analyzing gene-smoking interactions in a genome-wide association study of 1207 AMD cases and 686 controls of Caucasian background with genotype data on 668,238 single nucleotide polymorphisms (SNPs) after quality control. Participants' history of smoking at least 100 cigarettes lifetime was determined by a self-administered questionnaire. SNP associations modeled the effect of the minor allele additively on AMD using logistic regression, with adjustment for age, sex, and ever/never smoking. Joint effects of SNPs and smoking were examined comparing a null model containing only age, sex, and smoking against an extended model including genotypic and interaction terms. Genome-wide significant main effects were detected at three known AMD loci: CFH (P = 7.51×10(-30) ), ARMS2 (P = 1.94×10(-23) ), and RDBP/CFB/C2 (P = 4.37×10(-10) ), while joint effects analysis revealed three genomic regions with P < 10(-5) . Analyses stratified by smoking found genetic associations largely restricted to nonsmokers, with one notable exception: the chromosome 18q22.1 intergenic SNP rs17073641 (between SERPINB8 and CDH7), more strongly associated in nonsmokers (OR = 0.57, P = 2.73 × 10(-5) ), with an inverse association among smokers (OR = 1.42, P = 0.00228), suggesting that smoking modifies the effect of some genetic polymorphisms on AMD risk.

Linkage and association of successful aging to the 6q25 region in large Amish kindreds

Successful aging (SA) is a multidimensional phenotype involving living to older age with high physical function, preserved cognition, and continued social engagement. Several domains underlying SA are heritable, and identifying health-promoting polymorphisms and their interactions with the environment could provide important information regarding the health of older adults. In the present study, we examined 263 cognitively intact Amish individuals age 80 and older (74 SA and 189 "normally aged") all of whom are part of a single 13-generation pedigree. A genome-wide association study of 630,309 autosomal single nucleotide polymorphisms (SNPs) was performed and analyzed for linkage using multipoint analyses and for association using the modified quasi-likelihood score test. There was evidence for linkage on 6q25-27 near the fragile site FRA6E region with a dominant model maximum multipoint heterogeneity LOD score = 3.2. The 1-LOD-down support interval for this linkage contained one SNP for which there was regionally significant evidence of association (rs205990, p = 2.36 × 10(-5)). This marker survived interval-wide Bonferroni correction for multiple testing and was located between the genes QKI and PDE10A. Other areas of chromosome 6q25-q27 (including the FRA6E region) contained several SNPs associated with SA (minimum p = 2.89 × 10(-6)). These findings suggest potentially novel genes in the 6q25-q27 region linked and associated with SA in the Amish; however, these findings should be verified in an independent replication cohort.

Estrogen pathway polymorphisms in relation to primary open angle glaucoma: an analysis accounting for gender from the United States

PURPOSE

Circulating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender.

METHODS

We included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP <22 mmHg (normal pressure glaucoma [NPG]) at diagnosis. We conducted these analyses for each gender separately and then jointly in men and women.

RESULTS

Among women, the estrogen SNP pathway was associated with POAG overall (permuted p=0.006) and HPG (permuted p<0.001) but not NPG (permuted p=0.09). Interestingly, there was no relation between the estrogen SNP pathway and POAG when men were considered alone (permuted p>0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01).

CONCLUSIONS

The estrogen SNP pathway was associated with POAG among women.

C9ORF72 intermediate repeat copies are a significant risk factor for Parkinson disease

We set out to determine whether expansions in the C9ORF72 repeat found in amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) families are associated with Parkinson disease (PD). We determined the repeat size in a total of 889 clinically ascertained patients (including PD and essential tremor plus Parkinsonism (ETP)) and 1144 controls using a repeat-primed PCR assay. We found that large C9ORF72 repeat expansions (>30 repeats) were not contributing to PD risk. However, PD and ETP cases had a significant increase in intermediate (>20 to 30+) repeat copies compared to controls. Overall, 14 cases (13 PD, 1 ETP) and three controls had >20 repeat copies (Fisher's exact test p = 0.002). Further, seven cases and no controls had >23 repeat copies (p = 0.003). Our results suggest that intermediate copy numbers of the C9ORF72 repeat contribute to risk for PD and ETP. This also suggests that PD, ALS and FTD share some pathophysiological mechanisms of disease. Further studies are needed to elucidate the contribution of the C9ORF72 repeat in the overall PD population and to determine whether other common genetic risk factors exist between these neurodegenerative disorders.

Variants at chromosome 10q26 locus and the expression of HTRA1 in the retina

Variations in a locus at chromosome 10q26 are strongly associated with the risk of age-related macular degeneration (AMD). The most significantly associated haplotype includes a nonsynonymous SNP rs10490924 in the exon 1 of ARMS2 and rs11200638 in the promoter region of HTRA1. It is under debate which gene(s), ARMS2, HTRA1 or some other genes are functionally responsible for the genetic association. To verify whether the associated variants correlate with a higher HTRA1 expression level as previously reported, HTRA1 mRNA and protein were measured in a larger human retina-RPE-choroid samples (n = 82). Results show there is no significant change of HTRA1 mRNA level among genotypes at rs11200638, rs10490924 or an indel variant of ARMS2. Furthermore, two AMD-associated synonymous SNPs rs1049331 and rs2293870 in HTRA1 exon 1 do not change its protein level either. These results suggest that the AMD-associated variants in the chromosome 10q26 locus do not significantly affect the expression of HTRA1.

Parkinson disease loci in the mid-western Amish

Previous evidence has shown that Parkinson disease (PD) has a heritable component, but only a small proportion of the total genetic contribution to PD has been identified. Genetic heterogeneity complicates the verification of proposed PD genes and the identification of new PD susceptibility genes. Our approach to overcome the problem of heterogeneity is to study a population isolate, the mid-western Amish communities of Indiana and Ohio. We performed genome-wide association and linkage analyses on 798 individuals (31 with PD), who are part of a 4,998 member pedigree. Through these analyses, we identified a region on chromosome 5q31.3 that shows evidence of association (p value < 1 × 10(-4)) and linkage (multipoint HLOD = 3.77). We also found further evidence of linkage on chromosomes 6 and 10 (multipoint HLOD 4.02 and 4.35 respectively). These data suggest that locus heterogeneity, even within the Amish, may be more extensive than previously appreciated.

Functional genomics identifies type I interferon pathway as central for host defense against Candida albicans

Candida albicans is the most common human fungal pathogen causing mucosal and systemic infections. However, human antifungal immunity remains poorly defined. Here, by integrating transcriptional analysis and functional genomics, we identified Candida-specific host defense mechanisms in humans. Candida induced significant expression of genes from the type I interferon (IFN) pathway in human peripheral blood mononuclear cells. This unexpectedly prominent role of type I IFN pathway in anti-Candida host defense was supported by additional evidence. Polymorphisms in type I IFN genes modulated Candida-induced cytokine production and were correlated with susceptibility to systemic candidiasis. In in-vitro experiments, type I IFNs skewed Candida-induced inflammation from a Th17-response toward a Th1-response. Patients with chronic mucocutaneaous candidiasis displayed defective expression of genes in the type I IFN pathway. These findings indicate that the type I IFN pathway is a main signature of Candida-induced inflammation and plays a crucial role in anti-Candida host defense in humans.

High-resolution survey in familial Parkinson disease genes reveals multiple independent copy number variation events in PARK2

A high density comparative genomic hybridization array was designed to evaluate CNVs in the genomic region of six familial PD genes in 181 PD cases and 67 controls. No CNV was found in PARK7, ATP13A2, PINK1, and LRRK2. Intronic-only CNVs were found in SNCA and PARK2 but were not associated with PD risk. A whole-gene duplication of SNCA was found in one case. The allele frequency of PARK2 exonic CNV is significantly higher in cases than in controls (P = 0.02), higher in early-onset (AAO ≤ 40) than in late-onset cases (P = 0.001), and higher in familial than in sporadic cases (P = 0.005). Except for single exon 2 duplications, all PARK2 exonic CNVs have different breakpoints, even when the same exon(s) were involved. In conclusion, except for SNCA and PARK2, CNVs are not a major contributing mechanism for the familial PD genes examined. The majority of PARK2 exonic CNVs are not recurrent.

Evaluating power and type 1 error in large pedigree analyses of binary traits

Studying population isolates with large, complex pedigrees has many advantages for discovering genetic susceptibility loci; however, statistical analyses can be computationally challenging. Allelic association tests need to be corrected for relatedness among study participants, and linkage analyses require subdividing and simplifying the pedigree structures. We have extended GenomeSIMLA to simulate SNP data in complex pedigree structures based on an Amish pedigree to generate the same structure and distribution of sampled individuals. We evaluated type 1 error rates when no disease SNP was simulated and power when disease SNPs with recessive, additive, and dominant modes of inheritance and odds ratios of 1.1, 1.5, 2.0, and 5.0 were simulated. We generated subpedigrees with a maximum bit-size of 24 using PedCut and performed two-point and multipoint linkage using Merlin. We also ran MQLS on the subpedigrees and unified pedigree. We saw no inflation of type 1 error when running MQLS on either the whole pedigrees or the sub-pedigrees, and we saw low type 1 error for two-point and multipoint linkage. Power was reduced when running MQLS on the subpedigrees versus the whole pedigree, and power was low for two-point and multipoint linkage analyses of the subpedigrees. These data suggest that MQLS has appropriate type 1 error rates in our Amish pedigree structure, and while type 1 error does not seem to be affected when dividing the pedigree prior to linkage analysis, power to detect linkage is diminished when the pedigree is divided.

Seven new loci associated with age-related macular degeneration

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P < 5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P < 5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Whole exome sequencing of rare variants in EIF4G1 and VPS35 in Parkinson disease

OBJECTIVE

Recently, vacuolar protein sorting 35 (VPS35) and eukaryotic translation initiation factor 4 gamma 1 (EIF4G1) have been identified as 2 causal Parkinson disease (PD) genes. We used whole exome sequencing for rapid, parallel analysis of variations in these 2 genes.

METHODS

We performed whole exome sequencing in 213 patients with PD and 272 control individuals. Those rare variants (RVs) with <5% frequency in the exome variant server database and our own control data were considered for analysis. We performed joint gene-based tests for association using RVASSOC and SKAT (Sequence Kernel Association Test) as well as single-variant test statistics.

RESULTS

We identified 3 novel VPS35 variations that changed the coded amino acid (nonsynonymous) in 3 cases. Two variations were in multiplex families and neither segregated with PD. In EIF4G1, we identified 11 (9 nonsynonymous and 2 small indels) RVs including the reported pathogenic mutation p.R1205H, which segregated in all affected members of a large family, but also in 1 unaffected 86-year-old family member. Two additional RVs were found in isolated patients only. Whereas initial association studies suggested an association (p = 0.04) with all RVs in EIF4G1, subsequent testing in a second dataset for the driving variant (p.F1461) suggested no association between RVs in the gene and PD.

CONCLUSIONS

We confirm that the specific EIF4G1 variation p.R1205H seems to be a strong PD risk factor, but is nonpenetrant in at least one 86-year-old. A few other select RVs in both genes could not be ruled out as causal. However, there was no evidence for an overall contribution of genetic variability in VPS35 or EIF4G1 to PD development in our dataset.

Genome-wide association and linkage study in the Amish detects a novel candidate late-onset Alzheimer disease gene

To identify novel late-onset Alzheimer disease (LOAD) risk genes, we have analysed Amish populations of Ohio and Indiana. We performed genome-wide SNP linkage and association studies on 798 individuals (109 with LOAD). We tested association using the Modified Quasi-Likelihood Score test and also performed two-point and multipoint linkage analyses. We found that LOAD was significantly associated with APOE (P= 9.0 × 10-6) in all our ascertainment regions except for the Adams County, Indiana, community (P= 0.55). Genome-wide, the most strongly associated SNP was rs12361953 (P= 7.92 × 10-7). A very strong, genome-wide significant multipoint peak [recessive heterogeneity multipoint LOD (HLOD) = 6.14, dominant HLOD = 6.05] was detected on 2p12. Three additional loci with multipoint HLOD scores >3 were detected on 3q26, 9q31 and 18p11. Converging linkage and association results, the most significantly associated SNP under the 2p12 peak was at rs2974151 (P= 1.29 × 10-4). This SNP is located in CTNNA2, which encodes catenin alpha 2, a neuronal-specific catenin known to have function in the developing brain. These results identify CTNNA2 as a novel candidate LOAD gene, and implicate three other regions of the genome as novel LOAD loci. These results underscore the utility of using family-based linkage and association analyses in isolated populations to identify novel loci for traits with complex genetic architecture.

Population-based case-control association studies

This unit provides an overview of the design and analysis of population-based case-control studies of genetic risk factors for complex disease. Considerations specific to genetic studies are emphasized. The unit reviews basic study designs differentiating case-control studies from others, presents different genetic association strategies (candidate gene, genome-wide association, and high-throughput sequencing), introduces basic methods of statistical analysis for case-control data and approaches to combining case-control studies, and discusses measures of association and impact. Admixed populations, controlling for confounding (including population stratification), consideration of multiple loci and environmental risk factors, and complementary analyses of haplotypes, genes, and pathways are briefly discussed. Readers are referred to basic texts on epidemiology for more details on general conduct of case-control studies.

Retinoblastoma treatment: impact of the glycolytic inhibitor 2-deoxy-d-glucose on molecular genomics expression in LH(BETA)T(AG) retinal tumors

PURPOSE

The purpose of this study was to evaluate the effect of 2-deoxy-D-glucose (2-DG) on the spatial distribution of the genetic expression of key elements involved in angiogenesis, hypoxia, cellular metabolism, and apoptosis in LH(BETA)T(AG) retinal tumors.

METHODS

The right eye of each LH(BETA)T(AG) transgenic mouse (n = 24) was treated with either two or six subconjunctival injections of 2-DG (500 mg/kg) or saline control at 16 weeks of age. A gene expression array analysis was performed on five different intratumoral regions (apex, center, base, anterior-lateral, and posterior-lateral) using Affymetrix GeneChip Mouse Gene 1.0 ST arrays. To test for treatment effects of each probe within each region, a two-way analysis of variance was used.

RESULTS

Significant differences between treatment groups (ie, 0, 2, and 6 injections) were found as well as differences among the five retinal tumor regions evaluated (P < 0.01). More than 100 genes were observed to be dysregulated by ≥2-fold difference in expression between the three treatment groups, and their dysregulation varied across the five regions assayed. Several genes involved in pathways important for tumor cell growth (ie, angiogenesis, hypoxia, cellular metabolism, and apoptosis) were identified.

CONCLUSIONS

2-DG was found to significantly alter the gene expression in LH(BETA)T(AG) retinal tumor cells according to their location within the tumor as well as the treatment schedule. 2-DG's effects on genetic expression found here correlate with previous reported results on varied processes involved in its in vitro and in vivo activity in inhibiting tumor cell growth.

Human genetic susceptibility to Candida infections

Infections with Candida spp. have different manifestations in humans, ranging from mucosal to bloodstream and deep-seated disseminated infections. Immunocompromised patients have increased susceptibility to these types of infections, due to reduced capacity to elicit effective innate or adaptive immunity. In addition, rare and common genetic variants in the human genome have been identified that influence susceptibility to Candida infections. Genetic determinants of primary immunodeficiencies leading to chronic mucocutaneous candidiasis have been reported, and polymorphisms in genes that are known to be involved in anti-Candida host defense are associated with increased susceptibility to systemic infection. These findings have greatly increased our understanding of pathways important for anti-Candida defense in humans, and patterns of prevalence of Candida infections. In addition, these pathways may offer novel therapeutic targets for treatment. This review provides an overview of the current insights in genetic susceptibility to Candida infections and their consequences for the immune response against Candida.