Molecular Analysis of NOTCH2 in Patients with Primary Open-Angle Glaucoma.

BACKGROUND AND PURPOSE: Transgenic mice overexpressing Notch2 in the uvea exhibit a hyperplastic ciliary body leading to increased IOP and glaucoma. The aim of this study was to investigate the possible presence of NOTCH2 variants in patients with primary open-angle glaucoma (POAG). METHODS: We screened DNA samples from 130 patients with POAG for NOTCH2 variants by denaturing high-performance liquid chromatography after PCR amplification and validated our data by direct Sanger sequencing. RESULTS: No mutations were observed in the coding regions of NOTCH2 or in the splice sites. 19 known SNPs (single nucleotide polymorphisms) were detected. An SNP located in intron 24, c.[4005+45A>G], was seen in 28.5% of the patients (37/130 patients). As this SNP is reported to have a minor allele frequency of 7% in the 1000 genomes database, it could be associated with POAG. However, we evaluated its frequency in an ethnic-matched control group of 96 subjects unaffected by POAG and observed a frequency of 29%, indicating that it was not related to POAG. CONCLUSION: NOTCH2 seemed to be a good candidate for POAG as it is expressed in the anterior segment in the human eye. However, mutational analysis did not show any causative mutation. This study also shows that proper ethnic-matched control groups are essential in association studies and that values given in databases are sometimes misleading.

Influence of MCHR2 and MCHR2-AS1 Genetic Polymorphisms on Body Mass Index in Psychiatric Patients and In Population-Based Subjects with Present or Past Atypical Depression.

Obesity development during psychotropic treatments represents a major health issue in psychiatry. Melanin-concentrating hormone receptor 2 (MCHR2) is a central receptor involved in energy homeostasis. MCHR2 shares its promoter region with MCHR2-AS1, a long antisense non-coding RNA. The aim of this study was to determine whether tagging single nucleotide polymorphisms (tSNPs) of MCHR2 and MCHR2-AS1 are associated with the body mass index (BMI) in the psychiatric and in the general population. The influence of MCHR2 and MCHR2-AS1 tSNPs on BMI was firstly investigated in a discovery psychiatric sample (n1 = 474). Positive results were tested for replication in two other psychiatric samples (n2 = 164, n3 = 178) and in two population-based samples (CoLaus, n4 = 5409; GIANT, n5 = 113809). In the discovery sample, TT carriers of rs7754794C>T had 1.08 kg/m2 (p = 0.04) lower BMI as compared to C-allele carriers. This observation was replicated in an independent psychiatric sample (-2.18 kg/m2; p = 0.009). The association of rs7754794C>T and BMI seemed stronger in subjects younger than 45 years (median of age). In the population-based sample, a moderate association was observed (-0.17 kg/m2; p = 0.02) among younger individuals (<45y). Interestingly, this association was totally driven by patients meeting lifetime criteria for atypical depression, i.e. major depressive episodes characterized by symptoms such as an increased appetite. Indeed, patients with atypical depression carrying rs7754794-TT had 1.17 kg/m2 (p = 0.04) lower BMI values as compared to C-allele carriers, the effect being stronger in younger individuals (-2.50 kg/m2; p = 0.03; interaction between rs7754794 and age: p-value = 0.08). This study provides new insights on the possible influence of MCHR2 and/or MCHR2-AS1 on obesity in psychiatric patients and on the pathophysiology of atypical depression.

Fast and Rigorous Computation of Gene and Pathway Scores from SNP-Based Summary Statistics.

Integrating single nucleotide polymorphism (SNP) p-values from genome-wide association studies (GWAS) across genes and pathways is a strategy to improve statistical power and gain biological insight. Here, we present Pascal (Pathway scoring algorithm), a powerful tool for computing gene and pathway scores from SNP-phenotype association summary statistics. For gene score computation, we implemented analytic and efficient numerical solutions to calculate test statistics. We examined in particular the sum and the maximum of chi-squared statistics, which measure the strongest and the average association signals per gene, respectively. For pathway scoring, we use a modified Fisher method, which offers not only significant power improvement over more traditional enrichment strategies, but also eliminates the problem of arbitrary threshold selection inherent in any binary membership based pathway enrichment approach. We demonstrate the marked increase in power by analyzing summary statistics from dozens of large meta-studies for various traits. Our extensive testing indicates that our method not only excels in rigorous type I error control, but also results in more biologically meaningful discoveries.

Prevalence and determinants of periodic limb movements in the general population.

OBJECTIVE: Periodic limb movements during sleep (PLMS) are sleep phenomena characterized by periodic episodes of repetitive stereotyped limb movements. The aim of this study was to describe the prevalence and determinants of PLMS in a middle to older aged general population. METHODS: Data from 2,162 subjects (51.2% women, mean age = 58.4 ± 11.1 years) participating in a population-based study (HypnoLaus, Lausanne, Switzerland) were collected. Assessments included laboratory tests, sociodemographic data, personal and treatment history, and full polysomnography at home. PLMS index (PLMSI) was determined, and PLMSI > 15/h was considered as significant. RESULTS: Prevalence of PLMSI > 15/h was 28.6% (31.3% in men, 26% in women). Compared to subjects with PLMSI ≤ 15/h, subjects with PLMSI > 15/h were older (p < 0.001), were predominantly males (p = 0.007), had a higher proportion of restless legs syndrome (RLS; p < 0.001), had a higher body mass index (p = 0.001), and had a lower mean glomerular filtration rate (p < 0.001). Subjects with PLMSI > 15/h also had a higher prevalence of diabetes, hypertension, and beta-blocker or hypnotic treatments. The prevalence of antidepressant use was higher, but not statistically significant (p = 0.07). Single nucleotide polymorphisms (SNPs) within BTBD9 (rs3923809), TOX3 (rs3104788), and MEIS1 (rs2300478) genes were significantly associated with PLSMI > 15/h. Conversely, mean hemoglobin and ferritin levels were similar in both groups. In the multivariate analysis, age, male gender, antidepressant intake, RLS, and rs3923809, rs3104788, and rs2300478 SNPs were independently associated with PLMSI > 15/h. INTERPRETATION: PLMS are highly prevalent in our middle-aged European population. Age, male gender, RLS, antidepressant treatment, and specific BTBD9, TOX3, and MEIS1 SNP distribution are independent predictors of PLMSI > 15/h. ANN NEUROL 2016;79:464-474.

Genomic evidence for adaptive inversion clines in Drosophila melanogaster.

Clines in chromosomal inversion polymorphisms-presumably driven by climatic gradients-are common but there is surprisingly little evidence for selection acting on them. Here we address this long-standing issue in Drosophila melanogaster by using diagnostic single nucleotide polymorphism (SNP) markers to estimate inversion frequencies from 28 whole-genome Pool-seq samples collected from 10 populations along the North American east coast. Inversions In(3L)P, In(3R)Mo, and In(3R)Payne showed clear latitudinal clines, and for In(2L)t, In(2R)NS, and In(3R)Payne the steepness of the clinal slopes changed between summer and fall. Consistent with an effect of seasonality on inversion frequencies, we detected small but stable seasonal fluctuations of In(2R)NS and In(3R)Payne in a temperate Pennsylvanian population over 4 years. In support of spatially varying selection, we observed that the cline in In(3R)Payne has remained stable for >40 years and that the frequencies of In(2L)t and In(3R)Payne are strongly correlated with climatic factors that vary latitudinally, independent of population structure. To test whether these patterns are adaptive, we compared the amount of genetic differentiation of inversions versus neutral SNPs and found that the clines in In(2L)t and In(3R)Payne are maintained nonneutrally and independent of admixture. We also identified numerous clinal inversion-associated SNPs, many of which exhibit parallel differentiation along the Australian cline and reside in genes known to affect fitness-related traits. Together, our results provide strong evidence that inversion clines are maintained by spatially-and perhaps also temporally-varying selection. We interpret our data in light of current hypotheses about how inversions are established and maintained.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

Genetic variation in IL28B is associated with chronic hepatitis C and treatment failure: a genome-wide association study.

BACKGROUND & AIMS: Hepatitis C virus (HCV) induces chronic infection in 50% to 80% of infected persons; approximately 50% of these do not respond to therapy. We performed a genome-wide association study to screen for host genetic determinants of HCV persistence and response to therapy. METHODS: The analysis included 1362 individuals: 1015 with chronic hepatitis C and 347 who spontaneously cleared the virus (448 were coinfected with human immunodeficiency virus [HIV]). Responses to pegylated interferon alfa and ribavirin were assessed in 465 individuals. Associations between more than 500,000 single nucleotide polymorphisms (SNPs) and outcomes were assessed by multivariate logistic regression. RESULTS: Chronic hepatitis C was associated with SNPs in the IL28B locus, which encodes the antiviral cytokine interferon lambda. The rs8099917 minor allele was associated with progression to chronic HCV infection (odds ratio [OR], 2.31; 95% confidence interval [CI], 1.74-3.06; P = 6.07 x 10(-9)). The association was observed in HCV mono-infected (OR, 2.49; 95% CI, 1.64-3.79; P = 1.96 x 10(-5)) and HCV/HIV coinfected individuals (OR, 2.16; 95% CI, 1.47-3.18; P = 8.24 x 10(-5)). rs8099917 was also associated with failure to respond to therapy (OR, 5.19; 95% CI, 2.90-9.30; P = 3.11 x 10(-8)), with the strongest effects in patients with HCV genotype 1 or 4. This risk allele was identified in 24% of individuals with spontaneous HCV clearance, 32% of chronically infected patients who responded to therapy, and 58% who did not respond (P = 3.2 x 10(-10)). Resequencing of IL28B identified distinct haplotypes that were associated with the clinical phenotype. CONCLUSIONS: The association of the IL28B locus with natural and treatment-associated control of HCV indicates the importance of innate immunity and interferon lambda in the pathogenesis of HCV infection.

β-Arrestin2 influences the response to methadone in opioid-dependent patients.

β-Arrestin2 (ARRB2) is a component of the G-protein-coupled receptor complex and is involved in μ-opioid and dopamine D(2) receptor signaling, two central processes in methadone signal transduction. We analyzed 238 patients in methadone maintenance treatment (MMT) and identified a haplotype block (rs34230287, rs3786047, rs1045280 and rs2036657) spanning almost the entire ARRB2 locus. Although none of these single nucleotide polymorphisms (SNPs) leads to a change in amino-acid sequence, we found that for all the SNPs analyzed, with exception of rs34230287, homozygosity for the variant allele confers a nonresponding phenotype (n=73; rs1045280C and rs2036657G: OR=3.1, 95% CI=1.5-6.3, P=0.004; rs3786047A: OR=2.5, 95% CI=1.2-5.1, P=0.02) also illustrated by a 12-fold shorter period of negative urine screening (P=0.01). The ARRB2 genotype may thus contribute to the interindividual variability in the response to MMT and help to predict response to treatment.

The P450 oxidoreductase genotype is associated with CYP3A activity in vivo as measured by the midazolam phenotyping test.

CYP3A4, CYP3A5 and CYP3A7 are hepatic enzymes that metabolize about 50% of drugs on the market, with a large overlap in their specificities. We investigated the genetic bases that contribute to the variation of CYP3A activity. We phenotyped 251 individuals from two independent studies (182 patients treated with methadone and 69 patients with clozapine) for CYP3A activity using the midazolam phenotyping test and genotyped them for CYP3A4, CYP3A5, and CYP3A7 genetic variants, including the single nucleotide polymorphism (SNP) rs4646437C>T in intron 7 of CYP3A4. Owing to the fact that CYP enzymes require electron transfer through the P450 oxidoreductase (POR), and functional impairment has been shown for the POR*28 SNP, this polymorphism was also analysed. We show that CYP3A4, CYP3A5 and CYP3A7 genotypes, including the SNP rs4646437C>T, do not reflect the inter-individual variability of CYP3A activity (P>0.1). In contrast, POR*28 TT genotype presents a 1.6-fold increase in CYP3A activity compared with POR*28C carriers (n = 182, P = 0.004). This finding was replicated in the second independent dataset (n = 69, P = 0.04). The SNP POR*28 seems to be a better genetic marker of the variability of total CYP3A activity in vivo than CYP3A4, CYP3A5 and CYP3A7 genetic variants.

DQB1 locus alone explains most of the risk and protection in narcolepsy with cataplexy in Europe.

STUDY OBJECTIVE: Prior research has identified five common genetic variants associated with narcolepsy with cataplexy in Caucasian patients. To replicate and/or extend these findings, we have tested HLA-DQB1, the previously identified 5 variants, and 10 other potential variants in a large European sample of narcolepsy with cataplexy subjects. DESIGN: Retrospective case-control study. SETTING: A recent study showed that over 76% of significant genome-wide association variants lie within DNase I hypersensitive sites (DHSs). From our previous GWAS, we identified 30 single nucleotide polymorphisms (SNPs) with P < 10(-4) mapping to DHSs. Ten SNPs tagging these sites, HLADQB1, and all previously reported SNPs significantly associated with narcolepsy were tested for replication. PATIENTS AND PARTICIPANTS: For GWAS, 1,261 narcolepsy patients and 1,422 HLA-DQB1*06:02-matched controls were included. For HLA study, 1,218 patients and 3,541 controls were included. MEASUREMENTS AND RESULTS: None of the top variants within DHSs were replicated. Out of the five previously reported SNPs, only rs2858884 within the HLA region (P < 2x10(-9)) and rs1154155 within the TRA locus (P < 2x10(-8)) replicated. DQB1 typing confirmed that DQB1*06:02 confers an extraordinary risk (odds ratio 251). Four protective alleles (DQB1*06:03, odds ratio 0.17, DQB1*05:01, odds ratio 0.56, DQB1*06:09 odds ratio 0.21, DQB1*02 odds ratio 0.76) were also identified. CONCLUSION: An overwhelming portion of genetic risk for narcolepsy with cataplexy is found at DQB1 locus. Since DQB1*06:02 positive subjects are at 251-fold increase in risk for narcolepsy, and all recent cases of narcolepsy after H1N1 vaccination are positive for this allele, DQB1 genotyping may be relevant to public health policy.

Mutation screening of the glutamate cysteine ligase modifier (GCLM) gene in patients with schizophrenia.

BACKGROUND: Experimental evidences show that glutathione and its rate-limiting synthesizing enzyme, the glutamate-cysteine ligase (GCL), are involved in the pathogenesis of schizophrenia. Furthermore, genetic association has been previously reported between two single nucleotide polymorphisms lying in noncoding regions of glutamate cysteine ligase modifier (GCLM) gene, which specifies for the modifier subunit of GCL and schizophrenia. OBJECTIVE: We wanted to investigate the presence of GCLM true functional mutations, likely in linkage disequilibrium with the previously identified single nucleotide polymorphism alleles, in the same set of cases that allowed the detection of the original association signal. METHODS: We screened all the coding regions of GCLM and their intronic flanking vicinities in 353 patients with schizophrenia by direct DNA sequencing. RESULTS: Ten sequence variations were identified, five of which were not previously described. None of these DNA changes was within the GCLM coding sequence and in-silico analysis failed to indicate functional impairment induced by these variations. Furthermore, screening of normal controls and downstream statistical analyses revealed no significant relationship of any of these DNA variants with schizophrenia. CONCLUSION: It is unlikely that functional mutations in the GCLM gene could play a major role in genetic predisposition to schizophrenia and further studies will be required to assess its etiological function in the disease.

Repeated unidirectional introgression towards Populus balsamifera in contact zones of exotic and native poplars.

As the evolutionary significance of hybridization is largely dictated by its extent beyond the first generation, we broadly surveyed patterns of introgression across a sympatric zone of two native poplars (Populus balsamifera, Populus deltoides) in Quebec, Canada within which European exotic Populus nigra and its hybrids have been extensively planted since the 1800s. Single nucleotide polymorphisms (SNPs) that appeared fixed within each species were characterized by DNA-sequencing pools of pure individuals. Thirty-five of these diagnostic SNPs were employed in a high-throughput assay that genotyped 635 trees of different age classes, sampled from 15 sites with various degrees of anthropogenic disturbance. The degree of admixture within sampled trees was then assessed through Bayesian clustering of genotypes. Hybrids were present in seven of the populations, with 2.4% of all sampled trees showing spontaneous admixture. Sites with hybrids were significantly more disturbed than pure stands, while hybrids comprised both immature juveniles and trees of reproductive age. All three possible F1s were detected. Advanced-generation hybrids were consistently biased towards P. balsamifera regardless of whether hybridization had occurred with P. deltoides or P. nigra. Gene exchange between P. deltoides and P. nigra was not detected beyond the F1 generation; however, detection of a trihybrid demonstrates that even this apparent reproductive isolation does not necessarily result in an evolutionary dead end. Collectively, results demonstrate the natural fertility of hybrid poplars and suggest that introduced genes could potentially affect the genetic integrity of native trees, similar to that arising from introgression between natives.

Genetic architecture of ambulatory blood pressure in the general population: insights from cardiovascular gene-centric array.

Genetic determinants of blood pressure are poorly defined. We undertook a large-scale, gene-centric analysis to identify loci and pathways associated with ambulatory systolic and diastolic blood pressure. We measured 24-hour ambulatory blood pressure in 2020 individuals from 520 white European nuclear families (the Genetic Regulation of Arterial Pressure of Humans in the Community Study) and genotyped their DNA using the Illumina HumanCVD BeadChip array, which contains ≈50 000 single nucleotide polymorphisms in >2000 cardiovascular candidate loci. We found a strong association between rs13306560 polymorphism in the promoter region of MTHFR and CLCN6 and mean 24-hour diastolic blood pressure; each minor allele copy of rs13306560 was associated with 2.6 mm Hg lower mean 24-hour diastolic blood pressure (P=1.2×10(-8)). rs13306560 was also associated with clinic diastolic blood pressure in a combined analysis of 8129 subjects from the Genetic Regulation of Arterial Pressure of Humans in the Community Study, the CoLaus Study, and the Silesian Cardiovascular Study (P=5.4×10(-6)). Additional analysis of associations between variants in gene ontology-defined pathways and mean 24-hour blood pressure in the Genetic Regulation of Arterial Pressure of Humans in the Community Study showed that cell survival control signaling cascades could play a role in blood pressure regulation. There was also a significant overrepresentation of rare variants (minor allele frequency: <0.05) among polymorphisms showing at least nominal association with mean 24-hour blood pressure indicating that a considerable proportion of its heritability may be explained by uncommon alleles. Through a large-scale gene-centric analysis of ambulatory blood pressure, we identified an association of a novel variant at the MTHFR/CLNC6 locus with diastolic blood pressure and provided new insights into the genetic architecture of blood pressure.

Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.