A genome-wide association study provides evidence for association of chromosome 8p23 (MYP10) and 10q21.1 (MYP15) with high myopia in the French population.

PURPOSE. Myopia is a complex trait affected by both genetic and environmental factors. High myopia is associated with increased risk of sight-threatening eye disorders such as retinal detachment. The purpose of this genome-wide association study was to identify susceptibility genes contributing to high myopia in the French population. METHODS. High myopic cases were genotyped using Affymetrix SNP 6.0 chips and population controls were selected from the GABRIEL French dataset in which samples were genotyped by Illumina Human610 quad array. The association study was conducted using 152,234 single nucleotide polymorphisms that were present on both manufacturers' chips in 192 high myopic cases and 1064 controls to identify associated regions. Imputation was performed on peak regions. RESULTS. Associations were found at known myopia locus MYP10 on chromosome 8p23 and MYP15 on chromosome 10q21.1. Rs189798 (8p23) and rs10825992 (10q21.1) showed the strongest associations in these regions (P=6.32x10-7 and P=2.17x10-5, respectively). The imputed results at 8p23 showed 2 peaks of interest. The first spanned 30kb including rs189798 between MIR4660 and PPP1R3B with the most significant association at rs17155227 (P=1.07x10-10). The second novel peak was 4kb in length, encompassing MIR124-1 and the MSRA gene, with the strongest association at rs55864141 (P=1.30x10-7). The peak of imputed data at 10q21.1 was 70kb in length between ZWINT and MIR3924, with rs3107503 having the lowest P value (P=1.54x10-7). CONCLUSION. We provide evidence for the association of MYP10 at 8p23 and MYP15 at 10p21.1 with high myopia in the French population and refine these regions of association.

Conservation genetics of Ireland's sole population of the River water crowfoot (Ranunculus fluitans Lam.)

Populations of many freshwater species are becoming increasingly threatened as a result of a wide range of anthropogenically mediated factors. In the present study, we wanted to assess levels and patterns of genetic diversity in Ireland's sole population of the River water crowfoot (Ranunculus fluitans), which is restricted to a 12 km stretch of a single river, to assist the formation of conservation strategies. Analysis using amplified fragment length polymorphism (AFLP) indicated comparable levels of genetic diversity to those exhibited by a more extensive population of the species in England, and revealed no evidence of clonal reproduction. Allele-specific PCR analysis of five nuclear single nucleotide polymorphisms (SNPs) indicated no evidence of hybridization with its more abundant congener Ranunculus penicillatus, despite previous anecdotal reports of the occurrence of hybrids. Although the population currently exhibits healthy levels of genetic diversity and is not at risk of genetic assimilation via hybridization with R. penicillatus, it still remains vulnerable to other factors such as stochastic events and invasive species. © 2013 Elsevier B.V. All rights reserved.

Association Analysis of Dyslipidemia-Related Genes in Diabetic Nephropathy

Type 1 diabetes (T1D) increases risk of the development of microvascular complications and cardiovascular disease (CVD). Dyslipidemia is a common risk factor in the pathogenesis of both CVD and diabetic nephropathy (DN), with CVD identified as the primary cause of death in patients with DN. In light of this commonality, we assessed single nucleotide polymorphisms (SNPs) in thirty-seven key genetic loci previously associated with dyslipidemia in a T1D cohort using a casecontrol design. SNPs (n = 53) were genotyped using Sequenom in 1467 individuals with T1D (718 cases with proteinuric nephropathy and 749 controls without nephropathy i.e. normal albumin excretion). Cases and controls were white and recruited from the UK and Ireland. Association analyses were performed using PLINK to compare allele frequencies in cases and controls. In a sensitivity analysis, samples from control individuals with reduced renal function (estimated glomerular filtration rate,60 ml/min/1.73 m2) were excluded. Correction for multiple testing was performed by permutation testing. A total of 1394 samples passed quality control filters. Following regression analysis adjusted by collection center, gender, duration of diabetes, and average HbA1c, two SNPs were significantly associated with DN. rs4420638 in the APOC1 region (odds ratio [OR] = 1.51; confidence intervals [CI]: 1.19–1.91; P = 0.001) and rs1532624 in CETP (OR = 0.82; CI: 0.69–0.99; P = 0.034); rs4420638 was also significantly associated in a sensitivity analysis (P = 0.016) together with rs7679 (P = 0.027). However, no association was significant following correction for multiple testing. Subgroup analysis of end-stage renal disease status failed to reveal any association. Our results suggest common variants associated with dyslipidemia are not strongly associated with DN in T1D among white individuals. Our findings, cannot entirely exclude these key genes which are central to the process of dyslipidemia, from involvement in DN pathogenesis as our study had limited power to detect variants of small effect size. Analysis in larger independent cohorts is required.

Extensive local-scale gene flow and long-term population stability in the intertidal mollusc Katharina tunicata (Mollusca: Polyplacophora)

The dispersal capabilities of intertidal organisms may represent a key factor to their survival in the face of global warming, as species that cannot adapt to the various effects of climate change will have to migrate to track suitable habitat. Although species with pelagic larval phases might be expected to have a greater capacity for dispersal than those with benthic larvae, interspecies comparisons have shown that this is not always the case. Consequently, population genetic approaches are being increasingly used to gain insights into dispersal through studying patterns of gene flow. In the present study, we used nuclear single-nucleotide polymorphisms (SNPs) and mitochondrial DNA (mtDNA) sequencing to elucidate fine-scale patterns of genetic variation between populations of the Black Katy Chiton, Katharina tunicata, separated by 15-150 km in south-west Vancouver Island. Both the nuclear and mitochondrial data sets revealed no genetic differentiation between the populations studied, and an isolation-with-migration analysis indicated extensive local-scale gene flow, suggesting an absence of barriers to dispersal. Population demographic analysis also revealed long-term population stability through previous periods of climate change associated with the Pleistocene glaciations. Together, the findings of the present study suggest that this high potential for dispersal may allow K. tunicata to respond to current global warming by tracking suitable habitat, consistent with its long-term demographic stability through previous changes in the Earth's climate. (C) 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106, 589597.

Investigation of Genetic Variation in Scavenger Receptor Class B, Member 1 (SCARB1) and Association with Serum Carotenoids

Objective: To investigate association of scavenger receptor class B, member 1 (SCARB1) genetic variants with serum carotenoid levels of lutein (L) and zeaxanthin (Z) and macular pigment optical density (MPOD).
Design: A cross-sectional study of healthy adults aged 20 to 70.
Participants: We recruited 302 participants after local advertisement.
Methods: We measured MPOD by customized heterochromatic flicker photometry. Fasting blood samples were taken for serum L and Z measurement by high-performance liquid chromatography and lipoprotein analysis by spectrophotometric assay. Forty-seven single nucleotide polymorphisms (SNPs) across SCARB1 were genotyped using Sequenom technology. Association analyses were performed using PLINK to compare allele and haplotype means, with adjustment for potential confounding and correction for multiple comparisons by permutation testing. Replication analysis was performed in the TwinsUK and Carotenoids in Age-Related Eye Disease Study (CAREDS) cohorts.
Main Outcome Measures: Odds ratios for MPOD area, serum L and Z concentrations associated with genetic variations in SCARB1 and interactions between SCARB1 and gender.
Results: After multiple regression analysis with adjustment for age, body mass index, gender, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, smoking, and dietary L and Z levels, 5 SNPs were significantly associated with serum L concentration and 1 SNP with MPOD (P<0.01). Only the association between rs11057841 and serum L withstood correction for multiple comparisons by permutation testing (P<0.01) and replicated in the TwinsUK cohort (P = 0.014). Independent replication was also observed in the CAREDS cohort with rs10846744 (P = 2×10-4), an SNP in high linkage disequilibrium with rs11057841 (r2 = 0.93). No interactions by gender were found. Haplotype analysis revealed no stronger association than obtained with single SNP analyses.
Conclusions: Our study has identified association between rs11057841 and serum L concentration (24% increase per T allele) in healthy subjects, independent of potential confounding factors. Our data supports further evaluation of the role for SCARB1 in the transport of macular pigment and the possible modulation of age-related macular degeneration risk through combating the effects of oxidative stress within the retina.
Financial Disclosure(s): Proprietary or commercial disclosures may be found after the references. Ophthalmology 2013;120:1632–1640 © 2013 by the American Academy of Ophthalmology.

Trait directed de novo population transcriptome dissects genetic regulation of a balanced polymorphism in phosphorus nutrition/arsenate tolerance in a wild grass Holcus lanatus

The aim of this study was to characterize the transcriptome of a balanced polymorphism, under the regulation of a single gene, for phosphate fertilizer responsiveness/arsenate toler- ance in wild grass Holcus lanatus genotypes screened from the same habitat.

De novo transcriptome sequencing, RNAseq (RNA sequencing) and single nucleotide poly- morphism (SNP) calling were conducted on RNA extracted from H.lanatus. Roche 454 sequencing data were assembled into c. 22 000 isotigs, and paired-end Illumina reads for phosphorus-starved (P) and phosphorus-treated (P+) genovars of tolerant (T) and nontoler- ant (N) phenotypes were mapped to this reference transcriptome.

Heatmaps of the gene expression data showed strong clustering of each P+/P treated genovar, as well as clustering by N/T phenotype. Statistical analysis identified 87 isotigs to be significantly differentially expressed between N and T phenotypes and 258 between P+ and P treated plants. SNPs and transcript expression that systematically differed between N and T phenotypes had regulatory function, namely proteases, kinases and ribonuclear RNA- binding protein and transposable elements.

A single gene for arsenate tolerance led to distinct phenotype transcriptomes and SNP pro- files, with large differences in upstream post-translational and post-transcriptional regulatory genes rather than in genes directly involved in P nutrition transport and metabolism per se.

Association Study of 167 Candidate Genes for Schizophrenia Selected by a Multi-Domain Evidence-Based Prioritization Algorithm and Neurodevelopmental Hypothesis

Integrating evidence from multiple domains is useful in prioritizing disease candidate genes for subsequent testing. We ranked all known human genes (n = 3819) under linkage peaks in the Irish Study of High-Density Schizophrenia Families using three different evidence domains: 1) a meta-analysis of microarray gene expression results using the Stanley Brain collection, 2) a schizophrenia protein-protein interaction network, and 3) a systematic literature search. Each gene was assigned a domain-specific p-value and ranked after evaluating the evidence within each domain. For comparison to this
ranking process, a large-scale candidate gene hypothesis was also tested by including genes with Gene Ontology terms related to neurodevelopment. Subsequently, genotypes of 3725 SNPs in 167 genes from a custom Illumina iSelect array were used to evaluate the top ranked vs. hypothesis selected genes. Seventy-three genes were both highly ranked and involved in neurodevelopment (category 1) while 42 and 52 genes were exclusive to neurodevelopment (category 2) or highly ranked (category 3), respectively. The most significant associations were observed in genes PRKG1, PRKCE, and CNTN4 but no individual SNPs were significant after correction for multiple testing. Comparison of the approaches showed an excess of significant tests using the hypothesis-driven neurodevelopment category. Random selection of similar sized genes from two independent genome-wide association studies (GWAS) of schizophrenia showed the excess was unlikely by chance. In a further meta-analysis of three GWAS datasets, four candidate SNPs reached nominal significance. Although gene ranking using integrated sources of prior information did not enrich for significant results in the current experiment, gene selection using an a priori hypothesis (neurodevelopment) was superior to random selection. As such, further development of gene ranking strategies using more carefully selected sources of information is warranted.

Neutral Genomic Microevolution of a Recently Emerged Pathogen, Salmonella enterica Serovar Agona

Salmonella enterica serovar Agona has caused multiple food-borne outbreaks of gastroenteritis since it was first isolated in 1952. We analyzed the genomes of 73 isolates from global sources, comparing five distinct outbreaks with sporadic infections as well as food contamination and the environment. Agona consists of three lineages with minimal mutational diversity: only 846 single nucleotide polymorphisms (SNPs) have accumulated in the non-repetitive, core genome since Agona evolved in 1932 and subsequently underwent a major population expansion in the 1960s. Homologous recombination with other serovars of S. enterica imported 42 recombinational tracts (360 kb) in 5/143 nodes within the genealogy, which resulted in 3,164 additional SNPs. In contrast to this paucity of genetic diversity, Agona is highly diverse according to pulsed-field gel electrophoresis (PFGE), which is used to assign isolates to outbreaks. PFGE diversity reflects a highly dynamic accessory genome associated with the gain or loss (indels) of 51 bacteriophages, 10 plasmids, and 6 integrative conjugational elements (ICE/IMEs), but did not correlate uniquely with outbreaks. Unlike the core genome, indels occurred repeatedly in independent nodes (homoplasies), resulting in inaccurate PFGE genealogies. The accessory genome contained only few cargo genes relevant to infection, other than antibiotic resistance. Thus, most of the genetic diversity within this recently emerged pathogen reflects changes in the accessory genome, or is due to recombination, but these changes seemed to reflect neutral processes rather than Darwinian selection. Each outbreak was caused by an independent clade, without universal, outbreak-associated genomic features, and none of the variable genes in the pan-genome seemed to be associated with an ability to cause outbreaks. © 2013 Achtman et al

Effect of interleukin-6 polymorphisms on human longevity:a systematic review and meta-analysis

Several studies have assessed changes in frequency of -174 interleukin (IL)-6 single nucleotide polymorphism (SNP) with age. If IL-6 tracks with disability and age-related diseases, then there should be reduction, in the oldest old, of the frequency of homozygous GG subjects, who produce higher IL-6 levels. However, discordant results have been obtained. To explore the relationship between this polymorphism and longevity, we analyzed individual data on long-living subjects and controls from eight case-control studies conducted in Europeans, using meta-analysis. There was no significant difference in the IL-6 genotype between the oldest old and controls (Odds Ratio [OR]=0.96; 95% C.I.: 0.77-1.20; p=0.71), but there was significant between-study heterogeneity (I2=55.5%). In a subgroup analyses when male centenarians from the three Italian studies were included, the frequency of the IL-6 -174 GG genotype was significantly lower than the other genotypes (OR=0.49; 95% C.I.: 0.31-0.80; p=0.004), with no evidence of heterogeneity (I2=0%). Our data supports a negative association between the GG genotype of IL-6 SNP and longevity in Italian centenarians, with males who carry the genotype being two times less likely to reach extreme old age compared with subjects carrying CC or CG genotypes. These findings were not replicated in other European groups suggesting a possible interaction between genetics, sex and environment in reaching longevity.

Five authors reply

We read with interest the comments offered by Drs. Hughes and Bradley (1) on our systematic review (2). Four single nucleotide polymorphisms (SNPs), rs9332739 and rs547154 in the complement component 2 gene (C2) and rs4151667 and rs641153 in the complement factor B gene (CFB), were pooled. Hughes and Bradley point out that we omitted the most common variant, rs12614. In fact, rs12614 is in high linkage disequilibrium (LD) with rs641153, which was included, and the major allele of both of these SNPs is in the range of 90% (population code, CEU, in the International HapMap Project (http://hapmap.ncbi.nlm.nih.gov/)). Moreover, our review was initiated in September 2010, at which point only 4 studies had published associations with rs12614, whereas 14 studies (n = 11,378) were available for rs641153. While it is true that both SNPs are better analyzed as a haplotype, these data were simply not available for pooling.
Hughes and Bradley also point out that we obtained and pooled new data that were not previously published. While it is recommended that contact with authors be completed as part of a comprehensive meta-analysis, we acknowledge that these additional data were not previously published and peer reviewed and, hence, do not have the same level of transparency. However, given that sample collections often increase over time and that the instrumentation for genotyping is continually improving, we thought that it would be advantageous to use the most recent information; this is a subjective decision.
We also agree that the allele frequencies given by Kaur et al. (3) were exactly opposite to those expected and were suggestive of strand flipping. However, we specifically queried this with the lead author on 2 separate occasions and were assured it was not.
Hughes and Bradley do make an interesting suggestion that SNPs in high LD should be used as a gauge of genotyping quality in HuGE reviews. This is an interesting idea but difficult to put into practice as the r2 parameter they propose as a measure of LD has some unusual properties. Although r2 is a measure of LD, it is also linked to the allele frequency; even small differences in allele frequencies between 2 linked SNPs can reduce the r2 dramatically. Wray (4) explored these effects and found that, at a baseline allele frequency of 10%, even a difference in allele frequency between 2 SNPs as small as 2% can drop the r2 value below 0.8. This degree of allele frequency difference is consistent with what could be expected for sampling error. Furthermore, when we look at 2 linked dialleleic SNPs, giving 4 possible haplotypes, the absence of 1 haplotype dramatically reduces r2, despite the 2 loci being in high LD as measured by D'. In fact, this is the situation for rs12614 and rs641153, where the low frequency of 1 haplotype means that the r2 is 0.01 but the D' is 1.
Hughes and Bradley also suggest consideration of genotype call rate restrictions as an inclusion criterion for metaanalysis. This would be more appropriate when focusing on genetic variants per se, as considered within the context of a genome-wide association study or other specific genetic analysis where large numbers of SNPs are evaluated (5).
The concerns raised by Hughes and Bradley reflect the limited ability of a meta-analysis based on summary data to tease out inconsistencies best identified at the individual level. We agree that SNPs in LD should be evaluated, but this will not necessarily be straightforward. A move to make genetic data sets publicly available, as in the Database of Genotypes and Phenotypes (http://www.ncbi.nlm.nih.gov/ gap), is a step in the right direction for greater transparency.

Molecular Validation of the Schizophrenia Spectrum

Background: Early descriptive work and controlled family and adoption studies support the hypothesis that a range of personality and nonschizophrenic psychotic disorders aggregate in families of schizophrenic probands. Can we validate, using molecular polygene scores from genome-wide association studies (GWAS), this schizophrenia spectrum? Methods: The predictive value of polygenic findings reported by the Psychiatric GWAS Consortium (PGC) was applied to 4 groups of relatives from the Irish Study of High-Density Schizophrenia Families (ISHDSF; N = s) differing on their assignment within the schizophrenia spectrum. Genome-wide single nucleotide polymorphism data for affected and unaffected relatives were used to construct per-individual polygene risk scores based on the PGC stage-I results. We compared mean polygene scores in the ISHDSF with mean scores in ethnically matched population controls (N = 929). Results: The schizophrenia polygene score differed significantly across diagnostic categories and was highest in those with narrow schizophrenia spectrum, lowest in those with no psychiatric illness, and in-between in those classified in the intermediate, broad, and very broad schizophrenia spectrum. Relatives of all of these groups of affected subjects, including those with no diagnosis, had schizophrenia polygene scores significantly higher than the control sample. Conclusions: In the relatives of high-density families, the observed pattern of enrichment of molecular indices of schizophrenia risk suggests an underlying, continuous liability distribution and validates, using aggregate common risk alleles, a genetic basis for the schizophrenia spectrum disorders. In addition, as predicted by genetic theory, nonpsychotic members of multiply-affected schizophrenia families are significantly enriched for replicated, polygenic risk variants compared with the general population.

Genetic variation in the 6p22.3 gene DTNBP1, the human ortholog of the mouse dysbindin gene, is associated with schizophrenia

Prior evidence has supported the existence of multiple susceptibility genes for schizophrenia. Multipoint linkage analysis of the 270 Irish high-density pedigrees that we have studied, as well as results from several other samples, suggest that at least one such gene is located in region 6p24-21. In the present study, family-based association analysis of 36 simple sequence-length-polymorphism markers and of 17 SNP markers implicated two regions, separated by approximately 7 Mb. The first region, and the focus of this report, is 6p22.3. In this region, single-nucleotide polymorphisms within the 140-kb gene DTNBP1 (dystrobrevin-binding protein 1, or dysbindin) are strongly associated with schizophrenia. Uncorrected, empirical P values produced by the program TRANSMIT were significant (P

Identification of a high-risk haplotype for the dystrobrevin binding protein 1 (DTNBP1) gene in the Irish study of high-density schizophrenia families

A recent report showed significant associations between several SNPs in a previously unknown EST cluster with schizophrenia. (1). The cluster was identified as the human dystrobrevin binding protein 1 gene (DTNBP1) by sequence database comparisons and homology with mouse DTNBP1. (2). However, the linkage disequilibrium (LD) among the SNPs in DTNBP1 as well as the pattern of significant SNP-schizophrenia association was complex. This raised several questions such as the number of susceptibility alleles that may be involved and the size of the region where the actual disease mutation(s) could be located. To address these questions, we performed different single-marker tests on the 12 previously studied and 2 new SNPs in DTNBP1 that were re-scored using an improved procedure, and performed a variety of haplotype analyses. The sample consisted of 268 Irish multiplex families selected for high density of schizophrenia. Results suggested a simple structure where the LD in the target region could be explained by 6 haplotypes that together accounted for 96% of haplotype diversity in the whole sample. From these six, a single high-risk haplotype was identified that showed a significant association with schizophrenia and explained the pattern of significant findings in the analyses with individual markers. This haplotype was 30 kb long, had a large effect, could be measured with two tag SNPs only, had a frequency of 6% in our sample, seemed to be of relatively recent origin in evolutionary terms, and was equally distributed over Ireland. Implications of these findings for follow-up and replication studies are discussed.

Variants in the catechol-o-methyltransferase (COMT) gene are associated with schizophrenia in Irish high-density families

The enzyme catechol-o-methyltransferase (COMT) transfers a methyl group from adenosylmethionine to catecholamines including the neurotransmitters dopamine, epinephrine and norepinephrine. This methylation results in the degradation of catecholamines. The involvement of the COMT gene in the metabolic pathway of these neurotransmitters has made it an attractive candidate gene for many psychiatric disorders. In this article, we reported our study of association of COMT with schizophrenia in Irish families with a high density of schizophrenia. Three single nucleotide polymorphisms (SNPs) were genotyped for the 274 such families and within-family transmission disequilibrium tests were performed. SNP rs4680, which is the functional Val/Met polymorphism, showed modest association with the disease by the TRANSMIT, FBAT and PDT programs, while the other two SNPs were negative. These SNPs showed lower level of LDs with each other in the Irish subjects than in Ashkenazi Jews. Haplotype analysis indicated that a haplotype, haplotype A-G-A for SNPs rs737865-rs4680-rs165599, was preferentially transmitted to the affected subjects. This was different from the reported G-G-G haplotype found in Ashkenazi Jews, but both haplotypes shared the Val allele. We concluded that COMT gene is associated with schizophrenia and carries a small but significant risk to the susceptibility in the Irish subjects.