Mitochondrial Genome Analysis of Primary Open Angle Glaucoma Patients

Primary open angle glaucoma (POAG) is a multi-factorial optic disc neuropathy characterized by accelerating damage of the retinal ganglion cells and atrophy of the optic nerve head. The vulnerability of the optic nerve damage leading to POAG has been postulated to result from oxidative stress and mitochondrial dysfunction. In this study, we investigated the possible involvement of the mitochondrial genomic variants in 101 patients and 71 controls by direct sequencing of the entire mitochondrial genome. The number of variable positions in the mtDNA with respect to the revised Cambridge Reference Sequence (rCRS), have been designated “Segregating Sites”. The segregating sites present only in the patients or controls have been designated “Unique Segregating Sites (USS)”. The population mutation rate (θ = 4N_eμ) as estimated by Watterson’s θ (θ_w), considering only the USS, was significantly higher among the patients (p = 9.8×10⁻¹⁵) compared to controls. The difference in θ_w and the number of USS were more pronounced when restricted to the coding region (p<1.31×10⁻²¹ and p = 0.006607, respectively). Further analysis of the region revealed non-synonymous variations were significantly higher in Complex I among the patients (p = 0.0053). Similar trends were retained when USS was considered only within complex I (frequency 0.49 vs 0.31 with p<0.0001 and mutation rate p-value <1.49×10⁻⁴³) and ND5 within its gene cluster (frequency 0.47 vs 0.23 with p<0.0001 and mutation rate p-value <4.42×10⁻⁴⁷). ND5 is involved in the proton pumping mechanism. Incidentally, glaucomatous trabecular meshwork cells have been reported to be more sensitive to inhibition of complex I activity. Thus mutations in ND5, expected to inhibit complex I activity, could lead to generation of oxidative stress and favor glaucomatous condition.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

Genome-Wide Association Study Identifies a Novel Canine Glaucoma Locus

Glaucoma is an optic neuropathy and one of the leading causes of blindness. Its hereditary forms are classified into primary closed-angle (PCAG), primary open-angle (POAG) and primary congenital glaucoma (PCG). Although many loci have been mapped in human, only a few genes have been identified that are associated with the development of glaucoma and the genetic basis of the disease remains poorly understood. Glaucoma has also been described in many dog breeds, including Dandie Dinmont Terriers (DDT) in which it is a late-onset (>7 years) disease. We designed clinical and genetic studies to better define the clinical features of glaucoma in the DDT and to identify the genetic cause. Clinical diagnosis was based on ophthalmic examinations of the affected dogs and 18 additionally investigated unaffected DDTs. We collected DNA from over 400 DTTs and a genome wide association study was performed in a cohort of 23 affected and 23 controls, followed by a fine mapping, a replication study and candidate gene sequencing. The clinical study suggested that ocular abnormalities including abnormal iridocorneal angles and pectinate ligament dysplasia are common (50% and 72%, respectively) in the breed and the disease resembles human PCAG. The genetic study identified a novel 9.5 Mb locus on canine chromosome 8 including the 1.6 Mb best associated region (p = 1.63×10⁻¹⁰, OR = 32 for homozygosity). Mutation screening in five candidate genes did not reveal any causative variants. This study indicates that although ocular abnormalities are common in DDTs, the genetic risk for glaucoma is conferred by a novel locus on CFA8. The canine locus shares synteny to a region in human chromosome 14q, which harbors several loci associated with POAG and PCG. Our study reveals a new locus for canine glaucoma and ongoing molecular studies will likely help to understand the genetic etiology of the disease.     

Variant near ADAMTS9 Known to Associate with Type 2 Diabetes Is Related to Insulin Resistance in Offspring of Type 2 Diabetes Patients—EUGENE2 Study

Backround

A meta-analysis combining results from three genome-wide association studies and followed by large-scale replication identified six novel type 2 diabetes loci. Subsequent studies of the effect of these variants on estimates of the beta-cell function and insulin sensitivity have been inconclusive. We examined these variants located in or near the JAZF1 (rs864745), THADA (rs7578597), TSPAN8 (rs7961581), ADAMTS9 (rs4607103), NOTCH2 (rs10923931) and the CDC123/CAMK1D (rs12779790) genes for associations with measures of pancreatic beta-cell function and insulin sensitivity.

Methodology/Results

Oral and intravenous glucose stimulated insulin release (n = 849) and insulin sensitivity (n = 596) estimated from a hyperinsulinemic euglycemic clamp were measured in non-diabetic offspring of type 2 diabetic patients from five European populations. Assuming an additive genetic model the diabetes-associated major C-allele of rs4607103 near ADAMTS9 associated with reduced insulin-stimulated glucose uptake (p = 0.002) during a hyperinsulinemic euglycemic clamp. However, following intravenous and oral administration of glucose serum insulin release was increased in individuals with the C-allele (p = 0.003 and p = 0.01, respectively). A meta-analyse combining clamp and IVGTT data from a total of 905 non-diabetic individuals showed that the C-risk allele associated with decreased insulin sensitivity (p = 0.003) and increased insulin release (p = 0.002). The major T-allele of the intronic JAZF1 rs864745 conferring increased diabetes risk was associated with increased 2^nd phase serum insulin release during an IVGTT (p = 0.03), and an increased fasting serum insulin level (p = 0.001). The remaining variants did not show any associations with insulin response, insulin sensitivity or any other measured quantitative traits.

Conclusion

The present studies suggest that the diabetogenic impact of the C-allele of rs4607103 near ADAMTS9 may in part be mediated through decreased insulin sensitivity of peripheral tissues.     

Physical Activity Attenuates the Genetic Predisposition to Obesity in 20,000 Men and Women from EPIC-Norfolk Prospective Population Study

Background

We have previously shown that multiple genetic loci identified by genome-wide association studies (GWAS) increase the susceptibility to obesity in a cumulative manner. It is, however, not known whether and to what extent this genetic susceptibility may be attenuated by a physically active lifestyle. We aimed to assess the influence of a physically active lifestyle on the genetic predisposition to obesity in a large population-based study.

Methods and Findings

We genotyped 12 SNPs in obesity-susceptibility loci in a population-based sample of 20,430 individuals (aged 39–79 y) from the European Prospective Investigation of Cancer (EPIC)-Norfolk cohort with an average follow-up period of 3.6 y. A genetic predisposition score was calculated for each individual by adding the body mass index (BMI)-increasing alleles across the 12 SNPs. Physical activity was assessed using a self-administered questionnaire. Linear and logistic regression models were used to examine main effects of the genetic predisposition score and its interaction with physical activity on BMI/obesity risk and BMI change over time, assuming an additive effect for each additional BMI-increasing allele carried. Each additional BMI-increasing allele was associated with 0.154 (standard error [SE] 0.012) kg/m² (p = 6.73×10⁻³⁷) increase in BMI (equivalent to 445 g in body weight for a person 1.70 m tall). This association was significantly (p _interaction = 0.005) more pronounced in inactive people (0.205 [SE 0.024] kg/m² [p = 3.62×10⁻¹⁸; 592 g in weight]) than in active people (0.131 [SE 0.014] kg/m² [p = 7.97×10⁻²¹; 379 g in weight]). Similarly, each additional BMI-increasing allele increased the risk of obesity 1.116-fold (95% confidence interval [CI] 1.093–1.139, p = 3.37×10⁻²⁶) in the whole population, but significantly (p _interaction = 0.015) more in inactive individuals (odds ratio [OR] = 1.158 [95% CI 1.118–1.199; p = 1.93×10⁻¹⁶]) than in active individuals (OR = 1.095 (95% CI 1.068–1.123; p = 1.15×10⁻¹²]). Consistent with the cross-sectional observations, physical activity modified the association between the genetic predisposition score and change in BMI during follow-up (p _interaction = 0.028).

Conclusions

Our study shows that living a physically active lifestyle is associated with a 40% reduction in the genetic predisposition to common obesity, as estimated by the number of risk alleles carried for any of the 12 recently GWAS-identified loci. Please see later in the article for the Editors'' Summary     

Human Intelligence and Polymorphisms in the DNA Methyltransferase Genes Involved in Epigenetic Marking

Epigenetic mechanisms have been implicated in syndromes associated with mental impairment but little is known about the role of epigenetics in determining the normal variation in human intelligence. We measured polymorphisms in four DNA methyltransferases (DNMT1, DNMT3A, DNMT3B and DNMT3L) involved in epigenetic marking and related these to childhood and adult general intelligence in a population (n = 1542) consisting of two Scottish cohorts born in 1936 and residing in Lothian (n = 1075) or Aberdeen (n = 467). All subjects had taken the same test of intelligence at age 11yrs. The Lothian cohort took the test again at age 70yrs. The minor T allele of DNMT3L SNP 11330C>T (rs7354779) allele was associated with a higher standardised childhood intelligence score; greatest effect in the dominant analysis but also significant in the additive model (coefficient = 1.40_additive; 95%CI 0.22,2.59; p = 0.020 and 1.99_dominant; 95%CI 0.55,3.43; p = 0.007). The DNMT3L C allele was associated with an increased risk of being below average intelligence (OR 1.25_additive; 95%CI 1.05,1.51; p = 0.011 and OR 1.37_dominant; 95%CI 1.11,1.68; p = 0.003), and being in the lowest 40^th (p_additive = 0.009; p_dominant = 0.002) and lowest 30^th (p_additive = 0.004; p_dominant = 0.002) centiles for intelligence. After Bonferroni correction for the number variants tested the link between DNMT3L 11330C>T and childhood intelligence remained significant by linear regression and centile analysis; only the additive regression model was borderline significant. Adult intelligence was similarly linked to the DNMT3L variant but this analysis was limited by the numbers studied and nature of the test and the association was not significant after Bonferroni correction. We believe that the role of epigenetics in the normal variation in human intelligence merits further study and that this novel finding should be tested in other cohorts.     

The p53 Codon 72 PRO/PRO Genotype May Be Associated with Initial Central Visual Field Defects in Caucasians with Primary Open Angle Glaucoma

Background

Loss of vision in glaucoma is due to apoptotic retinal ganglion cell loss. While p53 modulates apoptosis, gene association studies between p53 variants and glaucoma have been inconsistent. In this study we evaluate the association between a p53 variant functionally known to influence apoptosis (codon 72 Pro/Arg) and the subset of primary open angle glaucoma (POAG) patients with early loss of central visual field.

Methods

Genotypes for the p53 codon 72 polymorphism (Pro/Arg) were obtained for 264 POAG patients and 400 controls from the U.S. and in replication studies for 308 POAG patients and 178 controls from Australia (GIST). The glaucoma patients were divided into two groups according to location of initial visual field defect (either paracentral or peripheral). All cases and controls were Caucasian with European ancestry.

Results

The p53-PRO/PRO genotype was more frequent in the U.S. POAG patients with early visual field defects in the paracentral regions compared with those in the peripheral regions or control group (p = 2.7×10⁻⁵). We replicated this finding in the GIST cohort (p  = 7.3×10⁻³, and in the pooled sample (p = 6.6×10⁻⁷) and in a meta-analysis of both the US and GIST datasets (1.3×10⁻⁶, OR 2.17 (1.58–2.98 for the PRO allele).

Conclusions

These results suggest that the p53 codon 72 PRO/PRO genotype is potentially associated with early paracentral visual field defects in primary open-angle glaucoma patients.     

SLC22A1-ABCB1 Haplotype Profiles Predict Imatinib Pharmacokinetics in Asian Patients with Chronic Myeloid Leukemia

Objective

This study aimed to explore the influence of SLC22A1, PXR, ABCG2, ABCB1 and CYP3A5*3 genetic polymorphisms on imatinib mesylate (IM) pharmacokinetics in Asian patients with chronic myeloid leukemia (CML).

Patients and Methods

Healthy subjects belonging to three Asian populations (Chinese, Malay, Indian; n = 70 each) and CML patients (n = 38) were enrolled in a prospective pharmacogenetics study. Imatinib trough (C_0h) and clearance (CL) were determined in the patients at steady state. Haplowalk method was applied to infer the haplotypes and generalized linear model (GLM) to estimate haplotypic effects on IM pharmacokinetics. Association of haplotype copy numbers with IM pharmacokinetics was defined by Mann-Whitney U test.

Results

Global haplotype score statistics revealed a SLC22A1 sub-haplotypic region encompassing three polymorphisms (rs3798168, rs628031 and IVS7+850C>T), to be significantly associated with IM clearance (p = 0.013). Haplotype-specific GLM estimated that the haplotypes AGT and CGC were both associated with 22% decrease in clearance compared to CAC [CL (*10⁻² L/hr/mg): CAC vs AGT: 4.03 vs 3.16, p = 0.017; CAC vs CGC: 4.03 vs 3.15, p = 0.017]. Patients harboring 2 copies of AGT or CGC haplotypes had 33.4% lower clearance and 50% higher C_0h than patients carrying 0 or 1 copy [CL (*10⁻² L/hr/mg): 2.19 vs 3.29, p = 0.026; C_0h (*10⁻⁶ 1/ml): 4.76 vs 3.17, p = 0.013, respectively]. Further subgroup analysis revealed SLC22A1 and ABCB1 haplotypic combinations to be significantly associated with clearance and C_0h (p = 0.002 and 0.009, respectively).

Conclusion

This exploratory study suggests that SLC22A1-ABCB1 haplotypes may influence IM pharmacokinetics in Asian CML patients.     

Glaucomatous Patterns in Frequency Doubling Technology (FDT) Perimetry Data Identified by Unsupervised Machine Learning Classifiers

Purpose

The variational Bayesian independent component analysis-mixture model (VIM), an unsupervised machine-learning classifier, was used to automatically separate Matrix Frequency Doubling Technology (FDT) perimetry data into clusters of healthy and glaucomatous eyes, and to identify axes representing statistically independent patterns of defect in the glaucoma clusters.

Methods

FDT measurements were obtained from 1,190 eyes with normal FDT results and 786 eyes with abnormal FDT results from the UCSD-based Diagnostic Innovations in Glaucoma Study (DIGS) and African Descent and Glaucoma Evaluation Study (ADAGES). For all eyes, VIM input was 52 threshold test points from the 24-2 test pattern, plus age.

Results

FDT mean deviation was −1.00 dB (S.D. = 2.80 dB) and −5.57 dB (S.D. = 5.09 dB) in FDT-normal eyes and FDT-abnormal eyes, respectively (p<0.001). VIM identified meaningful clusters of FDT data and positioned a set of statistically independent axes through the mean of each cluster. The optimal VIM model separated the FDT fields into 3 clusters. Cluster N contained primarily normal fields (1109/1190, specificity 93.1%) and clusters G₁ and G₂ combined, contained primarily abnormal fields (651/786, sensitivity 82.8%). For clusters G₁ and G₂ the optimal number of axes were 2 and 5, respectively. Patterns automatically generated along axes within the glaucoma clusters were similar to those known to be indicative of glaucoma. Fields located farther from the normal mean on each glaucoma axis showed increasing field defect severity.

Conclusions

VIM successfully separated FDT fields from healthy and glaucoma eyes without a priori information about class membership, and identified familiar glaucomatous patterns of loss.     

An Intronic SINE Insertion in FAM161A that Causes Exon-Skipping Is Associated with Progressive Retinal Atrophy in Tibetan Spaniels and Tibetan Terriers

Progressive retinal atrophy (PRA) in dogs is characterised by the degeneration of the photoreceptor cells of the retina, resulting in vision loss and eventually complete blindness. The condition affects more than 100 dog breeds and is known to be genetically heterogeneous between breeds. Around 19 mutations have now been identified that are associated with PRA in around 49 breeds, but for the majority of breeds the mutation(s) responsible have yet to be identified. Using genome-wide association with 22 Tibetan Spaniel PRA cases and 10 controls, we identified a novel PRA locus, PRA3, on CFA10 (p_raw = 2.01×10⁻⁵, p_genome = 0.014), where a 3.8 Mb region was homozygous within 12 cases. Using targeted next generation sequencing, a short interspersed nuclear element insertion was identified near a splice acceptor site in an intron of a provocative gene, FAM161A. Analysis of mRNA from an affected dog revealed that the SINE causes exon skipping, resulting in a frame shift, leading to a downstream premature termination codon and possibly a truncated protein product. This mutation segregates with the disease in 22 out of 35 cases tested (63%). Of the PRA controls, none are homozygous for the mutation, 15% carry the mutation and 85% are homozygous wildtype. This mutation was also identified in Tibetan Terriers, although our results indicate that PRA is genetically heterogeneous in both Tibetan Spaniels and Tibetan Terriers.     

Genome-Wide Association of Body Fat Distribution in African Ancestry Populations Suggests New Loci

Central obesity, measured by waist circumference (WC) or waist-hip ratio (WHR), is a marker of body fat distribution. Although obesity disproportionately affects minority populations, few studies have conducted genome-wide association study (GWAS) of fat distribution among those of predominantly African ancestry (AA). We performed GWAS of WC and WHR, adjusted and unadjusted for BMI, in up to 33,591 and 27,350 AA individuals, respectively. We identified loci associated with fat distribution in AA individuals using meta-analyses of GWA results for WC and WHR (stage 1). Overall, 25 SNPs with single genomic control (GC)-corrected p-values<5.0×10⁻⁶ were followed-up (stage 2) in AA with WC and with WHR. Additionally, we interrogated genomic regions of previously identified European ancestry (EA) WHR loci among AA. In joint analysis of association results including both Stage 1 and 2 cohorts, 2 SNPs demonstrated association, rs2075064 at LHX2, p = 2.24×10⁻⁸ for WC-adjusted-for-BMI, and rs6931262 at RREB1, p = 2.48×10⁻⁸ for WHR-adjusted-for-BMI. However, neither signal was genome-wide significant after double GC-correction (LHX2: p = 6.5×10⁻⁸; RREB1: p = 5.7×10⁻⁸). Six of fourteen previously reported loci for waist in EA populations were significant (p<0.05 divided by the number of independent SNPs within the region) in AA studied here (TBX15-WARS2, GRB14, ADAMTS9, LY86, RSPO3, ITPR2-SSPN). Further, we observed associations with metabolic traits: rs13389219 at GRB14 associated with HDL-cholesterol, triglycerides, and fasting insulin, and rs13060013 at ADAMTS9 with HDL-cholesterol and fasting insulin. Finally, we observed nominal evidence for sexual dimorphism, with stronger results in AA women at the GRB14 locus (p for interaction = 0.02). In conclusion, we identified two suggestive loci associated with fat distribution in AA populations in addition to confirming 6 loci previously identified in populations of EA. These findings reinforce the concept that there are fat distribution loci that are independent of generalized adiposity.     

Genome-Wide Analysis in German Shepherd Dogs Reveals Association of a Locus on CFA 27 with Atopic Dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1×10⁻⁵) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n _cases = 91, n _controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0×10⁻⁹), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (p_raw = 3.1×10⁻⁷, p_genome = 0.03). The total associated region was defined as a ∼1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ∼209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.     

Discovery and Functional Annotation of SIX6 Variants in Primary Open-Angle Glaucoma

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10⁻¹¹). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10⁻¹⁰) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.     

Timing and Intensity of Light Correlate with Body Weight in Adults

Light exposure can influence sleep and circadian timing, both of which have been shown to influence weight regulation. The goal of this study was to evaluate the relationship between ambient light, sleep and body mass index. Participants included 54 individuals (26 males, mean age 30.6, SD = 11.7 years). Light levels, sleep midpoint and duration were measured with wrist actigraphy (Actiwatch-L) for 7 days. BMI was derived from self-reported height and weight. Caloric intake was determined from 7 days of food logs. For each participant, light and activity data were output in 2 minute epochs, smoothed using a 5 point (10 minute) moving average and then aggregated over 24 hours. The mean light timing above 500 lux (MLiT⁵⁰⁰) was defined as the average clock time of all aggregated data points above 500 lux. MLiT⁵⁰⁰ was positively correlated with BMI (r = 0.51, p<0.001), and midpoint of sleep (r = 0.47, p<0.01). In a multivariable linear regression model including MLiT⁵⁰⁰ and midpoint of sleep, MLiT⁵⁰⁰ was a significant predictor of BMI (B = 1.26 SE = 0.34, β = 0.53 p = 0.001, r ² _Δ = 0.22). Adjusting for covariates, MLiT⁵⁰⁰ remained an independent predictor of BMI (B = 1.28 SE = 0.36, β = 0.54, p = 0.002, r ² _Δ = 0.20). The full model accounted for 34.7% of the variance in BMI (p = 0.01). Exposure to moderate levels of light at biologically appropriate times can influence weight, independent of sleep timing and duration.     

Hepatocyte Growth Factor Genetic Variations and Primary Angle-Closure Glaucoma in the Han Chinese Population

Purpose

The aim of this study is to examine whether or not hepatocyte growth factor (HGF) genetic variations are associated with susceptibility to primary angle-closure glaucoma (PACG) in the Han Chinese population.

Methods

Three single-nucleotide polymorphisms (SNPs)–rs5745718, rs17427817, and rs3735520–in the HGF gene were genotyped in 238 adult patients with PACG and 287 age-, sex-, and ethnically matched healthy controls by using a polymerase chain reaction restriction fragment length polymorphism assay. Data was analyzed by χ² analysis.

Results

The three tested analyzed polymorphisms in the HGF gene were in Hardy-Weinberg equilibrium, in all the subjects. The frequencies of the genotype and allele of rs5745718 and rs1742817 in the HGF gene were significantly different between the PACG patients and the controls. On one hand, the frequencies of the CC genotype and C allele of rs5745718 were significantly decreased in PACG patients compared with controls (P_c = 1.40×10⁻³; P_c = 3.21×10⁻⁴, respectively); however, on the other hand, significantly decreased frequencies of the GG genotype and the G allele of rs17427817 were observed in PACG patients compared with the controls (P_c = 0.006,; P_c = 6.06×10⁻⁴, respectively). A comparison of the distributions of the genotypes and alleles of rs3735520 showed no statistically significant differences between the PACG patients and the controls (p_c>0.05). The haplotype analysis results showed that the CGC haplotype frequency was significantly decreased in the patients with PACG compared with the controls (p_c<0.001). No difference was detected between the patients and the controls with regard to the other haplotypes.

Conclusions

Our study suggests that rs5745718 and rs17427817 are associated with a decreased risk of PACG in the Chinese Han population. The CGC haplotype was demonstrated to possibly play a protective role against PACG in this population.     

Highly Significant Association between Two Common Single Nucleotide Polymorphisms in CORIN Gene and Preeclampsia in Caucasian Women

Preeclampsia is a frequent medical complication during pregnancy. Corin, a serine protease which activates pro-atrial natriuretic peptide, has recently been shown to be involved in the pathophysiology of preeclampsia. The aim of this study was to search for CORIN gene variations and their association to preeclampsia in Caucasian and African women. Our study population was composed of 571 pregnant women (295 with preeclampsia and 276 normotensive controls) matched for maternal and gestational age, and ethnic origin. The 22 exons of the CORIN gene were sequenced in a discovery sample (n = 260), where 31 single nucleotide polymorphisms were identified. In a replication sample (n = 311), 4 single nucleotide polymorphisms were tested. Two minor alleles (C for rs2271036 and G for rs2271037) were significantly associated to preeclampsia. Adjusted odds ratios [95% confidence interval] were 2.5 [1.2–3.8] (p = 0.007) and 2.3 [1.5–3.5] (p = 1.3×10⁻⁴), respectively. These associations were ethnic-specific, as only found in the Caucasian of subjects (odds ratio = 3.5 [1.8–6.6], p = 1.1×10⁻⁴; odds ratio = 3.1 [1.7–5.8], p = 2.1×10⁻⁴, for each single nucleotide polymorphism, respectively). The two single nucleotide polymorphisms are in almost perfect linkage disequilibrium (r² = 0.93). No specific association was found with severe preeclampsia, early-onset preeclampsia nor fetal growth retardation. In conclusion, this is the first report of a highly significant association between these two single nucleotide polymorphisms in CORIN gene and preeclampsia. Our findings further support the probability of a critical role of corin in preeclamspia pathophysiology at the uteroplacental interface.     

Point Mutations Associated with Organophosphate and Carbamate Resistance in Chinese Strains of Culex pipiens quinquefasciatus (Diptera: Culicidae)

Acetylcholinesterase resistance has been well documented in many insects, including several mosquito species. We tested the resistance of five wild, Chinese strains of the mosquito Culex pipiens quinquefasciatus to two kinds of pesticides, dichlorvos and propoxur. An acetylcholinesterase gene (ace1) was cloned and sequenced from a pooled sample of mosquitoes from these five strains and the amino acids of five positions were found to vary (V185M, G247S, A328S, A391T, and T682A). Analysis of the correlation between mutation frequencies and resistance levels (LC₅₀) suggests that two point mutations, G247S (r² = 0.732, P = 0.065) and A328S (r² = 0.891, P = 0.016), are associated with resistance to propoxur but not to dichlorvos. Although the V185M mutation was not associated with either dichlorvos or propoxur resistance, its RS genotype frequency was correlated with propoxur resistance (r² = 0.815, P = 0.036). And the HWE test showed the A328S mutation is linked with V185M, also with G247S mutation. This suggested that these three mutations may contribute synergistically to propoxur resistance. The T682A mutation was negatively correlated with propoxur (r² = 0.788, P = 0.045) resistance. Knowledge of these mutations may help design strategies for managing pesticide resistance in wild mosquito populations.     

Impaired Chemosensitivity of Mouse Dorsal Raphe Serotonergic Neurons Overexpressing Serotonin 1A (Htr1a) Receptors

Background

Serotonergic system participates in a wide range of physiological processes and behaviors, but its role is generally considered as modulatory and noncrucial, especially concerning life-sustaining functions. We recently created a transgenic mouse line in which a functional deficit in serotonin homeostasis due to excessive serotonin autoinhibition was produced by inducing serotonin 1A receptor (Htr1a) overexpression selectively in serotonergic neurons (Htr1a raphe-overexpressing or Htr1a^RO mice). Htr1a^RO mice exhibit episodes of autonomic dysregulation, cardiovascular crises and death, resembling those of sudden infant death syndrome (SIDS) and revealing a life-supporting role of serotonergic system in autonomic control. Since midbrain serotonergic neurons are chemosensitive and are implicated in arousal we hypothesized that their chemosensitivity might be impaired in Htr1a^RO mice.

Principal findings

Loose-seal cell-attached recordings in brainstem slices revealed that serotonergic neurons in dorsal raphe nucleus of Htr1a^RO mice have dramatically reduced responses to hypercapnic challenge as compared with control littermates. In control mice, application of 9% CO₂ produced an increase in firing rate of serotonergic neurons (0.260±0.041 Hz, n = 20, p = 0.0001) and application of 3% CO₂ decreased their firing rate (−0.142±0.025 Hz, n = 17, p = 0.0008). In contrast, in Htr1a^RO mice, firing rate of serotonergic neurons was not significantly changed by 9% CO₂ (0.021±0.034 Hz, n = 16, p = 0.49) and by 3% CO₂ (0.012±0.046 Hz, n = 12, p = 0.97).

Conclusions

Our findings support the hypothesis that chemosensitivity of midbrain serotonergic neurons provides a physiological mechanism for arousal responses to life-threatening episodes of hypercapnia and that functional impairment, such as excessive autoinhibition, of midbrain serotonergic neuron responses to hypercapnia may contribute to sudden death.     

A Large Scale Gene-Centric Association Study of Lung Function in Newly-Hired Female Cotton Textile Workers with Endotoxin Exposure

Background

Occupational exposure to endotoxin is associated with decrements in pulmonary function, but how much variation in this association is explained by genetic variants is not well understood.

Objective

We aimed to identify single nucleotide polymorphisms (SNPs) that are associated with the rate of forced expiratory volume in one second (FEV₁) decline by a large scale genetic association study in newly-hired healthy young female cotton textile workers.

Methods

DNA samples were genotyped using the Illumina Human CVD BeadChip. Change rate in FEV₁ was modeled as a function of each SNP genotype in linear regression model with covariate adjustment. We controlled the type 1 error in study-wide level by permutation method. The false discovery rate (FDR) and the family-wise error rate (FWER) were set to be 0.10 and 0.15 respectively.

Results

Two SNPs were found to be significant (P<6.29×10⁻⁵), including rs1910047 (P = 3.07×10⁻⁵, FDR = 0.0778) and rs9469089 (P = 6.19×10⁻⁵, FDR = 0.0967), as well as other eight suggestive (P<5×10⁻⁴) associated SNPs. Gene-gene and gene-environment interactions were also observed, such as rs1910047 and rs1049970 (P = 0.0418, FDR = 0.0895); rs9469089 and age (P = 0.0161, FDR = 0.0264). Genetic risk score analysis showed that the more risk loci the subjects carried, the larger the rate of FEV₁ decline occurred (P _trend = 3.01×10⁻¹⁸). However, the association was different among age subgroups (P = 7.11×10⁻⁶) and endotoxin subgroups (P = 1.08×10⁻²). Functional network analysis illustrates potential biological connections of all interacted genes.

Conclusions

Genetic variants together with environmental factors interact to affect the rate of FEV₁ decline in cotton textile workers.