Implications of Ocular Confounding Factors for Aqueous Humor Proteomic and Metabolomic Analyses in Retinal Diseases.

Purpose: To assess the impact of ocular confounding factors on aqueous humor (AH) proteomic and metabolomic analyses for retinal disease characterization. Methods: This study recruited 138 subjects (eyes): 102 with neovascular age-related macular degeneration (nAMD), 18 with diabetic macular edema (DME), and 18 with cataract (control group). AH samples underwent analysis using Olink Target 96 proteomics and Metabolon's metabolomics platform Data analysis included correlation, differential abundance, and gene-set analysis. Results: In total, 756 proteins and 408 metabolites were quantified in AH. Total AH protein concentration was notably higher in nAMD (3.2-fold) and DME (4.1-fold) compared to controls. Pseudophakic eyes showed higher total AH protein concentrations than phakic eyes (e.g., 1.6-fold in nAMD) and a specific protein signature indicative of matrix remodeling. Unexpectedly, pupil-dilating drugs containing phenylephrine/tropicamide increased several AH proteins, notably interleukin-6 (5.4-fold in nAMD). Correcting for these factors revealed functionally relevant protein correlation clusters and disease-relevant, differentially abundant proteins across the groups. Metabolomics analysis, for which the relevance of confounder adjustment was less apparent, suggested insufficiently controlled diabetes and chronic hyperglycemia in the DME group. Conclusions: AH protein concentration, pseudophakia, and pupil dilation with phenylephrine/tropicamide are important confounding factors for AH protein analyses. When these factors are considered, AH analyses can more clearly reveal disease-relevant factors. Translational Relevance: Considering AH protein concentration, lens status, and phenylephrine/tropicamide administration as confounders is crucial for accurate interpretation of AH protein data.

Independent genome-wide scans identify a chromosome 18 quantitative-trait locus influencing dyslexia.

Developmental dyslexia is defined as a specific and significant impairment in reading ability that cannot be explained by deficits in intelligence, learning opportunity, motivation or sensory acuity. It is one of the most frequently diagnosed disorders in childhood, representing a major educational and social problem. It is well established that dyslexia is a significantly heritable trait with a neurobiological basis. The etiological mechanisms remain elusive, however, despite being the focus of intensive multidisciplinary research. All attempts to map quantitative-trait loci (QTLs) influencing dyslexia susceptibility have targeted specific chromosomal regions, so that inferences regarding genetic etiology have been made on the basis of very limited information. Here we present the first two complete QTL-based genome-wide scans for this trait, in large samples of families from the United Kingdom and United States. Using single-point analysis, linkage to marker D18S53 was independently identified as being one of the most significant results of the genome in each scan (P< or =0.0004 for single word-reading ability in each family sample). Multipoint analysis gave increased evidence of 18p11.2 linkage for single-word reading, yielding top empirical P values of 0.00001 (UK) and 0.0004 (US). Measures related to phonological and orthographic processing also showed linkage at this locus. We replicated linkage to 18p11.2 in a third independent sample of families (from the UK), in which the strongest evidence came from a phoneme-awareness measure (most significant P value=0.00004). A combined analysis of all UK families confirmed that this newly discovered 18p QTL is probably a general risk factor for dyslexia, influencing several reading-related processes. This is the first report of QTL-based genome-wide scanning for a human cognitive trait.

Whole Genome Sequencing Identifies Novel Common and Low-Frequency Variants Associated With Age-Related Macular Degeneration.

Purpose: To identify associations of common, low-frequency, and rare variants with advanced age-related macular degeneration (AMD) using whole genome sequencing (WGS). Methods: WGS data were obtained for 2123 advanced AMD patients (participants of clinical trials for advanced AMD) and 2704 controls (participants of clinical trials for asthma [N = 2518] and Alzheimer's disease [N = 186]), and joint genotype calling was performed, followed by quality control of the dataset. Single variant association analyses were performed for all identified common, low-frequency, and rare variants. Gene-based tests were executed for rare and low-frequency variants using SKAT-O and three groups of variants based on putative impact information: (1) all variants, (2) modifier impact variants, and (3) high- and moderate-impact variants. To ascertain independence of the identified associations from previously reported AMD and asthma loci, conditional analyses were performed. Results: Previously identified AMD variants at the CFH, ARMS2/HTRA1, APOE, and C3 loci were associated with AMD at a genome-wide significance level. We identified new single variant associations for common variants near the PARK7 gene and in the long non-coding RNA AC103876.1, and for a rare variant near the TENM3 gene. In addition, gene-based association analyses identified a burden of modifier variants in eight intergenic and gene-spanning regions and of high- and moderate-impact variants in the C3, CFHR5, SLC16A8, and CFI genes. Conclusions: We describe the largest WGS study in AMD to date. We confirmed previously identified associations and identified several novel associations that are worth exploring in further follow-up studies.

Correlation of Aqueous, Vitreous, and Serum Protein Levels in Patients With Retinal Diseases.

Purpose: To further establish aqueous humor (AH) as a clinically suitable source of protein biomarkers in retinal diseases by evaluating the correlation of a large panel of proteins between AH, vitreous humor (VH), and serum (SE). Methods: We enrolled 60 subjects (eyes) with various non-infectious retinal diseases. AH, VH, and SE proteins were analyzed using the Olink Target 96 platform (1196 protein assays in total). We compared these three matrices in terms of quantification overlap, principal component analysis, and correlation. Results: In the AH, VH, and SE samples, 841, 917, and 1133 proteins, respectively, were consistently quantified above the limit of detection in more than 30% of patients. AH and VH shared 812 of these proteins. AH and VH samples overlapped along principal component 1, but SE samples were distinct. We identified 490 proteins with significant (false discovery rate [FDR]-adjusted P < 0.05) and relevant correlations (correlation coefficient > 0.5) between AH and VH, compared to only 33 and 40 proteins for VH and SE and for AH and SE, respectively. Conclusions: Due to a close correlation between protein concentrations in the AH and VH and a clear difference from the SE, AH has the potential to serve as a substitute for VH and may hold significance in identifying protein biomarkers and novel targets related to retinal diseases. Translational Relevance: This study further supports AH as a clinically suitable source of protein biomarkers in retinal diseases. In addition, the identified AH and VH correlations can inform the selection of protein biomarker candidates in future translational research.

Longitudinal Evaluation of Visual Function Impairments in Early and Intermediate Age-Related Macular Degeneration Patients.

Purpose: To evaluate visual function (VF) changes in early and intermediate age-related macular degeneration (eAMD and iAMD) over 24 months. Design: Prospective, observational natural history study. Participants: Participants were enrolled at the Duke Eye Center. Methods: A total of 101 subjects (33 with eAMD, 47 with iAMD, and 21 normal controls) were recruited. Visual function (VF) tests included best-corrected visual acuity (BCVA), low- luminance visual acuity (LLVA), microperimetry (MP), cone contrast tests (CCTs), and dark adaptation (DA). Mixed-effect model repeated measures based on absolute values and change from baseline identified VF tests differentiating AMD from controls and revealing longitudinal VF decline when controlling for covariates (baseline value, age, coronary artery disease, dry eye, and phakic status). Nine AMD genetic risk variants, combinations of these (genetic burden score), reticular pseudodrusen (RPD), and hyperreflective foci (HRF) were tested as predictors of diagnosis and VF performance. Main Outcome Measures: Longitudinal changes in VF metrics over 24 months. Results: A total of 70 subjects completed the 2-year visit (22 with eAMD, 31 with iAMD, and 17 controls). Percent reduced threshold (PRT) on MP and CCT red significantly distinguished iAMD versus controls after 12 and 24 months, respectively. Cone contrast test red, PRT, and absolute threshold (AT) on MP showed significant longitudinal deterioration of VF in iAMD versus baseline at 12 months and onward, however, with a reduced rate of worsening. The DA data confirmed a preexisting functional deficit in iAMD at baseline and revealed an increasing proportion of poorly performing iAMD subjects in DA over the study period. None of the other VF measures showed consistent significant changes among the normal, early, and intermediate groups or over time. The genetic burden score was significantly associated with AMD diagnosis (relative risk for iAMD = 1.64, P < 0.01) and DA (r = 0.42, P = 0.00005). Reticular pseudodrusen and HRF showed moderate associations with DA and weak to moderate associations with MP variables. Conclusions: In iAMD, MP variables, CCT red, and DA revealed slow and nonlinear functional decline over 24 months. A structure-function relationship in eAMD and iAMD stages was demonstrated among HRF, RPD, and DA, possibly modified by genetic risk factors. These structural and functional features represent potential end points for clinical trials in iAMD.

Genetic association of complement receptor 1 polymorphism rs3818361 in Alzheimer's disease.

Complement receptor 1 gene polymorphism rs3818361 was recently shown to increase the risk of Alzheimer's disease (AD). We performed an independent replication study of this genetic variant in 2,470 individuals from Spain. By applying an allelic model, we observed a trend toward an association between this marker and late-onset AD susceptibility in our case-control study (odds ratio = 1.114, 95% confidence interval: 0.958-1.296, P = .16). Meta-analysis of available studies (n = 31,771 individuals), including previous studies and public genome-wide association study resources (Alzheimer's Disease Neuroimaging Initiative, Translational Genomics Research Institute, and Multi-site Collaborative Study for Genotype-Phenotype Associations in Alzheimer's Disease), strongly supports the effect of rs3818361 (odds ratio = 1.180, 95% confidence interval: 1.113-1.252, P < 2.99E-8) and suggests the existence of between-study heterogeneity (P < .05). We concluded that the complement receptor 1 gene may contribute to AD risk, although its effect size could be smaller than previously estimated.

"Does replication groups scoring reduce false positive rate in SNP interaction discovery? Response".

BACKGROUND: The genomewide evaluation of genetic epistasis is a computationally demanding task, and a current challenge in Genetics. HFCC (Hypothesis-Free Clinical Cloning) is one of the methods that have been suggested for genomewide epistasis analysis. In order to perform an exhaustive search of epistasis, HFCC has implemented several tools and data filters, such as the use of multiple replication groups, and direction of effect and control filters. A recent article has claimed that the use of multiple replication groups (as implemented in HFCC) does not reduce the false positive rate, and we hereby try to clarify these issues. RESULTS/DISCUSSION: HFCC uses, as an analysis strategy, the possibility of replicating findings in multiple replication groups, in order to select a liberal subset of preliminary results that are above a statistical criterion and consistent in direction of effect. We show that the use of replication groups and the direction filter reduces the false positive rate of a study, although at the expense of lowering the overall power of the study. A post-hoc analysis of these selected signals in the combined sample could then be performed to select the most promising results. CONCLUSION: Replication of results in independent samples is generally used in scientific studies to establish credibility in a finding. Nonetheless, the combined analysis of several datasets is known to be a preferable and more powerful strategy for the selection of top signals. HFCC is a flexible and complete analysis tool, and one of its analysis options combines these two strategies: a preliminary multiple replication group analysis to eliminate inconsistent false positive results, and a post-hoc combined-group analysis to select the top signals.

Genetic structure of the Spanish population.

BACKGROUND: Genetic admixture is a common caveat for genetic association analysis. Therefore, it is important to characterize the genetic structure of the population under study to control for this kind of potential bias. RESULTS: In this study we have sampled over 800 unrelated individuals from the population of Spain, and have genotyped them with a genome-wide coverage. We have carried out linkage disequilibrium, haplotype, population structure and copy-number variation (CNV) analyses, and have compared these estimates of the Spanish population with existing data from similar efforts. CONCLUSIONS: In general, the Spanish population is similar to the Western and Northern Europeans, but has a more diverse haplotypic structure. Moreover, the Spanish population is also largely homogeneous within itself, although patterns of micro-structure may be able to predict locations of origin from distant regions. Finally, we also present the first characterization of a CNV map of the Spanish population. These results and original data are made available to the scientific community.

WWOX gene is associated with HDL cholesterol and triglyceride levels.

BACKGROUND: Altered lipid profile, and in particular low HDL and high triglyceride (TG) plasma levels, are within the major determinants of cardiovascular diseases. The identification of quantitative trait loci (QTL) affecting these lipid levels is a relevant issue for predictive purposes. The WWOX gene has been recently associated with HDL levels. This gene is located at chromosome 16q23, a region previously linked to familial combined hyperlipidemia (FCHL) and HDL. Our objective is to perform a genetic association analysis at the WWOX gene region with HDL, TG and TG/HDL ratio. METHODS: A quantitative association analysis performed in 801 individuals selected from the Spanish general population. RESULTS: For HDL levels, two regions of intron 8 display clustering of positive signals (p < 0.05) but none of them was associated in the haplotypic analysis (0.07 ≤ p ≤ 0.165). For TG levels not only intron 8 but also a 27 kb region spanning from the promoter region to intron 4 are associated in this study. For the TG/HDL genetic association analysis, positive signals are coincident with those of the isolated traits. Interestingly, haplotypic analysis at the 5' region showed that variation in this region modified both HDL and TG levels, especially the latter (p = 0.003). CONCLUSIONS: Our results suggest that WWOX is a QTL for both TG and HDL.

An entropy test for single-locus genetic association analysis.

BACKGROUND: The etiology of complex diseases is due to the combination of genetic and environmental factors, usually many of them, and each with a small effect. The identification of these small-effect contributing factors is still a demanding task. Clearly, there is a need for more powerful tests of genetic association, and especially for the identification of rare effects RESULTS: We introduce a new genetic association test based on symbolic dynamics and symbolic entropy. Using a freely available software, we have applied this entropy test, and a conventional test, to simulated and real datasets, to illustrate the method and estimate type I error and power. We have also compared this new entropy test to the Fisher exact test for assessment of association with low-frequency SNPs. The entropy test is generally more powerful than the conventional test, and can be significantly more powerful when the genotypic test is applied to low allele-frequency markers. We have also shown that both the Fisher and Entropy methods are optimal to test for association with low-frequency SNPs (MAF around 1-5%), and both are conservative for very rare SNPs (MAF<1%) CONCLUSIONS: We have developed a new, simple, consistent and powerful test to detect genetic association of biallelic/SNP markers in case-control data, by using symbolic dynamics and symbolic entropy as a measure of gene dependence. We also provide a standard asymptotic distribution of this test statistic. Given that the test is based on entropy measures, it avoids smoothed nonparametric estimation. The entropy test is generally as good or even more powerful than the conventional and Fisher tests. Furthermore, the entropy test is more computationally efficient than the Fisher's Exact test, especially for large number of markers. Therefore, this entropy-based test has the advantage of being optimal for most SNPs, regardless of their allele frequency (Minor Allele Frequency (MAF) between 1-50%). This property is quite beneficial, since many researchers tend to discard low allele-frequency SNPs from their analysis. Now they can apply the same statistical test of association to all SNPs in a single analysis., which can be especially helpful to detect rare effects.

Genomewide linkage scan reveals novel loci modifying age of onset of Huntington's disease in the Venezuelan HD kindreds.

The age of onset of Huntington's disease (HD) is inversely correlated with the CAG length in the HD gene. The CAG repeat length accounts for 70% of the variability in HD age of onset. However, 90% of individuals worldwide with expanded alleles possess between 40 and 50 CAG repeat lengths in their HD gene. For these people, the size of their repeat only determines 44% of the variability in their age of onset. Once the effect of the CAG repeat has been accounted for, the residual variance in age of onset is a heritable trait. Targeted candidate gene studies and a genome scan have suggested some loci as potential modifiers of the age of onset of HD. We analyzed the large Venezuelan kindreds in which the HD gene was originally identified. These kindreds offer greater analytic power than standard sib-pair designs. We developed novel pedigree-member selection procedures to maximize power. Using a 5,858-single-nucleotide-polymorphism marker panel, we performed a genomewide linkage analysis. We discovered two novel loci on chromosome 2. Chromosome 2p25 (logarithm of the odds ratio (LOD)=4.29) and 2q35 (LOD=3.39) may contain genes that modify age of onset. A third linkage peak on chromosome 6q22 (LOD=2.48) may confirm the most promising locus from a previous genome scan. Two other candidate loci are suggestive on chromosome 5 (LOD=3.31 at 5p14 and LOD=3.14 at 5q32). All these regions harbor candidate genes that are potential HD modifier genes. Finding these modifier genes can reveal accessible and promising new therapeutic pathways and targets to ameliorate and cure HD.

A method for detecting epistasis in genome-wide studies using case-control multi-locus association analysis.

BACKGROUND: The difficulty in elucidating the genetic basis of complex diseases roots in the many factors that can affect the development of a disease. Some of these genetic effects may interact in complex ways, proving undetectable by current single-locus methodology. RESULTS: We have developed an analysis tool called Hypothesis Free Clinical Cloning (HFCC) to search for genome-wide epistasis in a case-control design. HFCC combines a relatively fast computing algorithm for genome-wide epistasis detection, with the flexibility to test a variety of different epistatic models in multi-locus combinations. HFCC has good power to detect multi-locus interactions simulated under a variety of genetic models and noise conditions. Most importantly, HFCC can accomplish exhaustive genome-wide epistasis search with large datasets as demonstrated with a 400,000 SNP set typed on a cohort of Parkinson's disease patients and controls. CONCLUSION: With the current availability of genetic studies with large numbers of individuals and genetic markers, HFCC can have a great impact in the identification of epistatic effects that escape the standard single-locus association analyses.

Longitudinal Study of Language and Speech of Twins at 4 and 6 Years: Twinning Effects Decrease, Zygosity Effects Disappear, and Heritability Increases.

Purpose: This study investigates the heritability of language, speech, and nonverbal cognitive development of twins at 4 and 6 years of age. Possible confounding effects of twinning and zygosity, evident at 2 years, were investigated among other possible predictors of outcomes. Method: The population-based twin sample included 627 twin pairs and 1 twin without a co-twin (197 monozygotic and 431 dizygotic), 610 boys and 645 girls, 1,255 children in total. Nine phenotypes from the same comprehensive direct behavioral assessment protocol were investigated at 4 and 6 years of age. Twinning effects were estimated for each phenotype at each age using general linear mixed models using maximum likelihood. Results: Twinning effects decreased from 4 to 6 years; zygosity effects disappeared by 6 years. Heritability increased from 4 to 6 years across all 9 phenotypes, and the heritability estimates were higher than reported previously, in the range of .44-.92 at 6 years. The highest estimate, .92, was for the clinical grammar marker. Conclusions: Across multiple dimensions of speech, language, and nonverbal cognition, heritability estimates are robust. A finiteness marker of grammar shows the highest inherited influences in this early period of children's language acquisition.

A mega-analysis of expression quantitative trait loci (eQTL) provides insight into the regulatory architecture of gene expression variation in liver.

Genome-wide association studies (GWAS) have identified numerous genetic variants in the human genome associated with diseases and traits. Nevertheless, for most loci the causative variant is still unknown. Expression quantitative trait loci (eQTL) in disease relevant tissues is an excellent approach to correlate genetic association with gene expression. While liver is the primary site of gene transcription for two pathways relevant to age-related macular degeneration (AMD), namely the complement system and cholesterol metabolism, we explored the contribution of AMD associated variants to modulate liver gene expression. We extracted publicly available data and computed the largest eQTL data set for liver tissue to date. Genotypes and expression data from all studies underwent rigorous quality control. Subsequently, Matrix eQTL was used to identify significant local eQTL. In total, liver samples from 588 individuals revealed 202,489 significant eQTL variants affecting 1,959 genes (Q-Value < 0.001). In addition, a further 101 independent eQTL signals were identified in 93 of the 1,959 eQTL genes. Importantly, our results independently reinforce the notion that high density lipoprotein metabolism plays a role in AMD pathogenesis. Taken together, our study generated a first comprehensive map reflecting the genetic regulatory landscape of gene expression in liver.

Combined Plasma and Cerebrospinal Fluid Signature for the Prediction of Midterm Progression From Mild Cognitive Impairment to Alzheimer Disease.

Importance: A reliable method of detecting Alzheimer disease (AD) in its prodromal state is needed for patient stratification in clinical trials or for personalizing existing or potential upcoming therapies. Current cerebrospinal fluid (CSF)- or imaging-based single biomarkers for AD offer reliable identification of patients with underlying AD but insufficient prediction of the rate of AD progression. Objective: To optimize prediction of progression from mild cognitive impairment (MCI) to AD dementia by combining information from diverse patient variables. Design, Setting, and Participants: This cohort study from the Alzheimer Disease Neuroimaging Initiative (ADNI) enrolled 928 patients with MCI at baseline and 249 selected variables available in the ADNI data set. Variables included clinical and demographic data, cognitive scores, magnetic resonance imaging-based brain volumetric data, the apolipoprotein E (APOE) and translocase of outer mitochondrial membrane 40 homolog (TOMM40) genotypes, and analyte levels measured in the CSF and plasma. Data were collected in July 2012 and analyzed from July 1, 2012, to June 1, 2015. Main Outcomes and Measures: Progression from MCI to AD within 1 to 6 years. To determine whether combinations of markers could predict progression from MCI to AD within 1 to 6 years, the elastic net algorithm was used in an iterative resampling of a training- and test-based variable selection and modeling approach. Results: Among the 928 patients with MCI in the ADNI database, 94 had 224 of the required variables available for the modeling. The results showed the contributions of age, Clinical Dementia Rating Sum of Boxes composite test score, hippocampal volume, and multiple plasma and CSF factors in modeling progression to AD. A combination of apolipoprotein A-II and cortisol levels in plasma and fibroblast growth factor 4, heart-type fatty acid binding protein, calcitonin, and tumor necrosis factor-related apoptosis-inducing ligand receptor 3 (TRAIL-R3) in CSF allowed for reliable prediction of disease status 3 years from the time of sample collection (80% classification accuracy, 88% sensitivity, and 70% specificity). Conclusions and Relevance: These study findings suggest that a combination of markers measured in plasma and CSF, distinct from ß-amyloid and tau, could prove useful in predicting midterm progression from MCI to AD dementia. Such a large-scale, multivariable-based analytical approach could be applied to other similar large data sets involving AD and beyond.

Quantitative trait locus for reading disability on chromosome 6p is pleiotropic for attention-deficit/hyperactivity disorder.

Comorbidity is pervasive among both adult and child psychiatric disorders; however, the etiological mechanisms underlying the majority of comorbidities are unknown. This study used genetic linkage analysis to assess the etiology of comorbidity between reading disability (RD) and attention-deficit hyperactivity disorder (ADHD), two common childhood disorders that frequently co-occur. Sibling pairs (N = 85) were ascertained initially because at least one individual in each pair exhibited a history of reading difficulties. Univariate linkage analyses in sibling pairs selected for ADHD from within this RD-ascertained sample suggested that a quantitative trait locus (QTL) on chromosome 6p is a susceptibility locus for ADHD. Because this QTL is in the same region as a well-replicated QTL for reading disability, subsequent bivariate analyses were conducted to test if this QTL contributed to comorbidity between the two disorders. Analyses of data from sib pairs selected for reading deficits revealed suggestive bivariate linkage for ADHD and three measures of reading difficulty, indicating that comorbidity between RD and ADHD may be due at least in part to pleiotropic effects of a QTL on chromosome 6p.

Late language emergence in 24-month-old twins: heritable and increased risk for late language emergence in twins.

PURPOSE: This study investigated the etiology of late language emergence (LLE) in 24-month-old twins, considering possible twinning, zygosity, gender, and heritability effects for vocabulary and grammar phenotypes. METHOD: A population-based sample of 473 twin pairs participated. Multilevel modeling estimated means and variances of vocabulary and grammar phenotypes, controlling for familiality. Heritability was estimated with DeFries-Fulker regression and variance components models to determine effects of heritability, shared environment, and nonshared environment. RESULTS: Twins had lower average language scores than norms for single-born children, with lower average performance for monozygotic than dizygotic twins and for boys than girls, although gender and zygosity did not interact. Gender did not predict LLE. Significant heritability was detected for vocabulary (0.26) and grammar phenotypes (0.52 and 0.43 for boys and girls, respectively) in the full sample and in the sample selected for LLE (0.42 and 0.44). LLE and the appearance of Word Combinations were also significantly heritable (0.22-0.23). CONCLUSIONS: The findings revealed an increased likelihood of LLE in twin toddlers compared with single-born children that is modulated by zygosity and gender differences. Heritability estimates are consistent with previous research for vocabulary and add further suggestion of heritable differences in early grammar acquisition.

A colorectal cancer susceptibility new variant at 4q26 in the Spanish population identified by genome-wide association analysis.

BACKGROUND: Non-hereditary colorectal cancer (CRC) is a complex disorder resulting from the combination of genetic and non-genetic factors. Genome-wide association studies (GWAS) are useful for identifying such genetic susceptibility factors. However, the single loci so far associated with CRC only represent a fraction of the genetic risk for CRC development in the general population. Therefore, many other genetic risk variants alone and in combination must still remain to be discovered. The aim of this work was to search for genetic risk factors for CRC, by performing single-locus and two-locus GWAS in the Spanish population. RESULTS: A total of 801 controls and 500 CRC cases were included in the discovery GWAS dataset. 77 single nucleotide polymorphisms (SNP)s from single-locus and 243 SNPs from two-locus association analyses were selected for replication in 423 additional CRC cases and 1382 controls. In the meta-analysis, one SNP, rs3987 at 4q26, reached GWAS significant p-value (p = 4.02×10(-8)), and one SNP pair, rs1100508 CG and rs8111948 AA, showed a trend for two-locus association (p = 4.35×10(-11)). Additionally, our GWAS confirmed the previously reported association with CRC of five SNPs located at 3q36.2 (rs10936599), 8q24 (rs10505477), 8q24.21(rs6983267), 11q13.4 (rs3824999) and 14q22.2 (rs4444235). CONCLUSIONS: Our GWAS for CRC patients from Spain confirmed some previously reported associations for CRC and yielded a novel candidate risk SNP, located at 4q26. Epistasis analyses also yielded several novel candidate susceptibility pairs that need to be validated in independent analyses.

Genetic study of neurexin and neuroligin genes in Alzheimer's disease.

The interaction between neurexins and neuroligins promotes the formation of functional synaptic structures. Recently, it has been reported that neurexins and neuroligins are proteolytically processed by presenilins at synapses. Based on this interaction and the role of presenilins in familial Alzheimer's disease (AD), we hypothesized that dysfunction of the neuroligin-neurexin pathway might be associated with AD. To explore this hypothesis, we carried out a meta-analysis of five genome-wide association studies (GWAS) comprising 1, 256 SNPs in the NRXN1, NRXN2, NRXN3, and NLGN1 genes (3,009 cases and 3,006 control individuals). We identified a marker in the NRXN3 gene (rs17757879) that showed a consistent protective effect in all GWAS, however, the statistical significance obtained did not resist multiple testing corrections (OR = 0.851, p = 0.002). Nonetheless, gender analysis revealed that this effect was restricted to males. A combined meta-analysis of the former five GWAS together with a replication Spanish sample consisting of 1,785 cases and 1,634 controls confirmed this observation (rs17757879, OR = 0.742, 95% CI = 0.632-0.872, p = 0.00028, final meta-analysis). We conclude that NRXN3 might have a role in susceptibility to AD in males.

Exploratory analysis of seven Alzheimer's disease genes: disease progression.

The relationships between genome wide association study-identified and replicated genetic variants associated with Alzheimer's disease (AD) risk and disease progression or therapeutic responses in AD patients are almost unexplored. Seven hundred and one AD patients with at least 3 different cognitive evaluations and genotypic information for APOE and 6 genome wide association study-significant single-nucleotide polymorphisms were selected for this study. Mean differences in Global Deterioration Score and Mini Mental State Examination (MMSE) were evaluated using nonparametric tests, general linear model and mixed models for repeated measurements. Each chart was also reviewed for evidence of treatment with any cholinesterase inhibitor, memantine, or both. Relationships between therapeutic protocols, genetic markers, and progression were explored using stratified analysis looking for specific effects on progression in each therapeutic category separately. Neither calculation rendered a Bonferroni-corrected statistically significant difference in any genetic marker. Mixed model results suggested differences in the average point in MMSE test for patients carrying PICALM GA or AA genotype compared with GG carriers at the end of the follow-up (MMSE mean difference = -0.57; 95% confidence interval, -1.145 to 0.009; p = 0.047). This observation remained unaltered after covariate adjustments although it did not achieve predefined multiple testing significance threshold. The PICALM single-nucleotide polymorphism also displayed a significant effect protecting against rapid progression during pharmacogenetic assays although its observed effect displayed heterogeneity among AD therapeutic protocols (p = 0.039). None of the studied genetic markers were convincingly linked to AD progression or drug response. However, by using different statistical approaches, the PICALM rs3851179 marker displayed consistent but weak effects on disease progression phenotypes.