A Genome-Wide Association Study for Diabetic Retinopathy in a Japanese Population: Potential Association with a Long Intergenic Non-Coding RNA

Elucidation of the genetic susceptibility factors for diabetic retinopathy (DR) is important to gain insight into the pathogenesis of DR, and may help to define genetic risk factors for this condition. In the present study, we conducted a three-stage genome-wide association study (GWAS) to identify DR susceptibility loci in Japanese patients, which comprised a total of 837 type 2 diabetes patients with DR (cases) and 1,149 without DR (controls). From the stage 1 genome-wide scan of 446 subjects (205 cases and 241 controls) on 614,216 SNPs, 249 SNPs were selected for the stage 2 replication in 623 subjects (335 cases and 288 controls). Eight SNPs were further followed up in a stage 3 study of 297 cases and 620 controls. The top signal from the present association analysis was rs9362054 in an intron of RP1-90L14.1 showing borderline genome-wide significance (P_met = 1.4×10⁻⁷, meta-analysis of stage 1 and stage 2, allele model). RP1-90L14.1 is a long intergenic non-coding RNA (lincRNA) adjacent to KIAA1009/QN1/CEP162 gene; CEP162 plays a critical role in ciliary transition zone formation before ciliogenesis. The present study raises the possibility that the dysregulation of ciliary-associated genes plays a role in susceptibility to DR.     

Genome-Wide Association Study of Pancreatic Cancer in Japanese Population

Pancreatic cancer shows very poor prognosis and is the fifth leading cause of cancer death in Japan. Previous studies indicated some genetic factors contributing to the development and progression of pancreatic cancer; however, there are limited reports for common genetic variants to be associated with this disease, especially in the Asian population. We have conducted a genome-wide association study (GWAS) using 991 invasive pancreatic ductal adenocarcinoma cases and 5,209 controls, and identified three loci showing significant association (P-value<5×10⁻⁷) with susceptibility to pancreatic cancer. The SNPs that showed significant association carried estimated odds ratios of 1.29, 1.32, and 3.73 with 95% confidence intervals of 1.17–1.43, 1.19–1.47, and 2.24–6.21; P-value of 3.30×10⁻⁷, 3.30×10⁻⁷, and 4.41×10⁻⁷; located on chromosomes 6p25.3, 12p11.21 and 7q36.2, respectively. These associated SNPs are located within linkage disequilibrium blocks containing genes that have been implicated some roles in the oncogenesis of pancreatic cancer.     

Genome-Wide Association Study of Peripheral Arterial Disease in a Japanese Population

Characteristics of peripheral arterial disease (PAD) are the occlusion or stenosis of multiple vessel sites caused mainly by atherosclerosis and chronic lower limb ischemia. To identify PAD susceptible loci, we conducted a genome-wide association study (GWAS) with 785 cases and 3,383 controls in a Japanese population using 431,666 single nucleotide polymorphisms (SNP). After staged analyses including a total of 3,164 cases and 20,134 controls, we identified 3 novel PAD susceptibility loci at IPO5/RAP2A, EDNRA and HDAC9 with genome wide significance (combined P = 6.8 x 10⁻¹⁴, 5.3 x 10⁻⁹ and 8.8 x 10⁻⁸, respectively). Fine-mapping at the IPO5/RAP2A locus revealed that rs9584669 conferred risk of PAD. Luciferase assay showed that the risk allele at this locus reduced expression levels of IPO5. To our knowledge, these are the first genetic risk factors for PAD.     

Genome-Wide Association Study of Breast Cancer in the Japanese Population

Breast cancer is the most common malignancy among women in worldwide including Japan. Several studies have identified common genetic variants to be associated with the risk of breast cancer. Due to the complex linkage disequilibrium structure and various environmental exposures in different populations, it is essential to identify variants associated with breast cancer in each population, which subsequently facilitate the better understanding of mammary carcinogenesis. In this study, we conducted a genome-wide association study (GWAS) as well as whole-genome imputation with 2,642 cases and 2,099 unaffected female controls. We further examined 13 suggestive loci (P<1.0×10⁻⁵) using an independent sample set of 2,885 cases and 3,395 controls and successfully validated two previously-reported loci, rs2981578 (combined P-value of 1.31×10⁻¹², OR = 1.23; 95% CI = 1.16–.30) on chromosome 10q26 (FGFR2), rs3803662 (combined P-value of 2.79×10⁻¹¹, OR = 1.21; 95% CI = 1.15–.28) and rs12922061 (combined P-value of 3.97×10⁻¹⁰, OR = 1.23; 95% CI = 1.15–.31) on chromosome 16q12 (TOX3-LOC643714). Weighted genetic risk score on the basis of three significantly associated variants and two previously reported breast cancer associated loci in East Asian population revealed that individuals who carry the most risk alleles in category 5 have 2.2 times higher risk of developing breast cancer in the Japanese population than those who carry the least risk alleles in reference category 1. Although we could not identify additional loci associated with breast cancer, our study utilized one of the largest sample sizes reported to date, and provided genetic status that represent the Japanese population. Further local and international collaborative study is essential to identify additional genetic variants that could lead to a better, accurate prediction for breast cancer.     

Genome-Wide Association Study for Levels of Total Serum IgE Identifies HLA-C in a Japanese Population

Most of the previously reported loci for total immunoglobulin E (IgE) levels are related to Th2 cell-dependent pathways. We undertook a genome-wide association study (GWAS) to identify genetic loci responsible for IgE regulation. A total of 479,940 single nucleotide polymorphisms (SNPs) were tested for association with total serum IgE levels in 1180 Japanese adults. Fine-mapping with SNP imputation demonstrated 6 candidate regions: the PYHIN1/IFI16, MHC classes I and II, LEMD2, GRAMD1B, and chr13∶60576338 regions. Replication of these candidate loci in each region was assessed in 2 independent Japanese cohorts (n = 1110 and 1364, respectively). SNP rs3130941 in the HLA-C region was consistently associated with total IgE levels in 3 independent populations, and the meta-analysis yielded genome-wide significance (P = 1.07×10⁻¹⁰). Using our GWAS results, we also assessed the reproducibility of previously reported gene associations with total IgE levels. Nine of 32 candidate genes identified by a literature search were associated with total IgE levels after correction for multiple testing. Our findings demonstrate that SNPs in the HLA-C region are strongly associated with total serum IgE levels in the Japanese population and that some of the previously reported genetic associations are replicated across ethnic groups.     

A Genome-Wide Association Study Identifies Potential Susceptibility Loci for Hirschsprung Disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.     

Association of VEGF Gene Polymorphisms with Diabetic Retinopathy: A Meta-Analysis

Background

Studies on the association of vascular endothelial growth factor (VEGF) gene -460T/C and -2578C/A polymorphisms with diabetic retinopathy (DR) have reported conflicting results. The aim of the present study was to assess the association by using meta-analysis.

Methods

A systematic search of electronic databases (PubMed, EMBASE, Elsevier Science Direct, ISI Web of Science, CBM, CNKI and VIP) was carried out until Sept 18, 2013. The pooled odds ratios (ORs) and their corresponding 95% confidence intervals (CIs) were used to assess the strength of the association.

Results

Eleven studies (-460T/C: 6 studies including 932 cases and 722 controls; -2578C/A: 6 studies including 1,071 cases and 1,137 controls) were involved in this meta-analysis. Significant association was found for -460T/C polymorphism (C versus T: OR=1.48, 95%CI=1.07–2.05, P=0.02; TC+CC versus TT: OR=1.78, 95%CI=1.02–3.12, P=0.04; CC versus TT+TC: OR=1.76, 95%CI=1.10–2.81, P=0.02), but not for -2578C/A polymorphism (P>0.05). Similar results were found in the subgroup analysis.

Conclusions

This meta-analysis demonstrates that DR is associated with VEGF gene -460T/C polymorphism, but not -2578C/A polymorphism. Further case-control studies based on larger sample size are still needed, especially for -2578C/A polymorphism.     

Genome-Wide Association Study for Wool Production Traits in a Chinese Merino Sheep Population

Genome-wide association studies (GWAS) provide a powerful approach for identifying quantitative trait loci without prior knowledge of location or function. To identify loci associated with wool production traits, we performed a genome-wide association study on a total of 765 Chinese Merino sheep (JunKen type) genotyped with 50 K single nucleotide polymorphisms (SNPs). In the present study, five wool production traits were examined: fiber diameter, fiber diameter coefficient of variation, fineness dispersion, staple length and crimp. We detected 28 genome-wide significant SNPs for fiber diameter, fiber diameter coefficient of variation, fineness dispersion, and crimp trait in the Chinese Merino sheep. About 43% of the significant SNP markers were located within known or predicted genes, including YWHAZ, KRTCAP3, TSPEAR, PIK3R4, KIF16B, PTPN3, GPRC5A, DDX47, TCF9, TPTE2, EPHA5 and NBEA genes. Our results not only confirm the results of previous reports, but also provide a suite of novel SNP markers and candidate genes associated with wool traits. Our findings will be useful for exploring the genetic control of wool traits in sheep.     

Genome-Wide Association Study for Autism Spectrum Disorder in Taiwanese Han Population

Background

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with strong genetic components. Several recent genome-wide association (GWA) studies in Caucasian samples have reported a number of gene regions and loci correlated with the risk of ASD—albeit with very little consensus across studies.

Methods

A two-stage GWA study was employed to identify common genetic variants for ASD in the Taiwanese Han population. The discovery stage included 315 patients with ASD and 1,115 healthy controls, using the Affymetrix SNP array 6.0 platform for genotyping. Several gene regions were then selected for fine-mapping and top markers were examined in extended samples. Single marker, haplotype, gene-based, and pathway analyses were conducted for associations.

Results

Seven SNPs had p-values ranging from 3.4~9.9*10⁻⁶, but none reached the genome-wide significant level. Five of them were mapped to three known genes (OR2M4, STYK1, and MNT) with significant empirical gene-based p-values in OR2M4 (p = 3.4*10⁻⁵) and MNT (p = 0.0008). Results of the fine-mapping study showed single-marker associations in the GLIS1 (rs12082358 and rs12080993) and NAALADL2 (rs3914502 and rs2222447) genes, and gene-based associations for the OR2M3-OR2T5 (olfactory receptor genes, p = 0.02), and GLIPR1/KRR1 gene regions (p = 0.015). Pathway analyses revealed important pathways for ASD, such as olfactory and G protein–coupled receptors signaling pathways.

Conclusions

We reported Taiwanese Han specific susceptibility genes and variants for ASD. However, further replication in other Asian populations is warranted to validate our findings. Investigation in the biological functions of our reported genetic variants might also allow for better understanding on the underlying pathogenesis of autism.     

Long Non-coding RNAs: Pivotal Epigenetic Regulators in Diabetic Retinopathy

Diabetic retinopathy (DR) is a severe complication of diabetes; however, its mechanism is not fully understood. Evidence has recently revealed that long non-coding RNAs (lncRNAs) are abnormally expressed in DR, and lncRNAs may function as pivotal regulators. LncRNAs are able to modulate gene expression at the epigenetic level by acting as scaffolds of histone modification complexes and sponges of binding with microRNAs (miRNAs). LncRNAs are believed to be important epigenetic regulators, which may become beneficial in the diagnosis and therapy of DR. However, the mechanisms of lncRNAs in DR are still unclear. In this review, we summarize the possible functions and mechanisms of lncRNAs in epigenetic regulation to target genes in the progression of DR.     

Genome-Wide Association Study Identifies Variants in PMS1 Associated with Serum Ferritin in a Chinese Population

Only a small proportion of genetic variation in serum ferritin has been explained by variant genetic studies, and genome-wide association study (GWAS) for serum ferritin has not been investigated widely in Chinese population. We aimed at exploring the novel genetic susceptibility to serum ferritin, and performed this two stage GWAS in a healthy Chinese population of 3,495 men aged 20–69 y, including 1,999 unrelated subjects in the first stage and 1,496 independent individuals in the second stage. Serum ferritin was measured with electrochemiluminescence immunoassay, and DNA samples were collected for genotyping. A total of 1,940,243 SNPs were tested by using multivariate linear regression analysis. After adjusting for population stratification, age and BMI, the rs5742933 located in the 5′UTR region of PMS1 gene on chromosome 2 was the most significantly associated with ferritin concentrations (P-combined = 2.329×10⁻¹⁰) (β = −0.11, 95% CI: −0.14, −0.07). Moreover, this marker was about 200kb away from the candidate gene SLC40A1 which is responsible for iron export. PMS1 gene was the novel genetic susceptibility to serum ferritin in Chinese males and its relation to SLC40A1 needs further study.     

Inverse Association between High Blood 25-Hydroxyvitamin D Levels and Diabetic Retinopathy in a Representative Korean Population

Purpose

To investigate the association between 25-hydroxyvitamin D and diabetic retinopathy (DR).

Methods

A population-based cross-sectional study using a nation-wide, systemically stratified, multistage, clustered sampling method included a total of 18,363 subjects aged ≥40 years who participated in the Korean National Health and Nutrition Examination Survey during 2008–2012. All participants participated in standardized interviews, blood 25-hydroxyvitamin D level evaluations, and comprehensive ophthalmic examinations. Seven standard retinal fundus photographs were obtained from both eyes after pupil dilatation. DR was graded according to the modified Airlie House classification system.

Results

The blood 25-hydroxyvitamin D levels were 19.2 ng/mL in men and 17.9 ng/mL in women. After adjusting for potential confounders, including age, sex, diabetes duration, hemoglobin A1c levels, and hypertension, the odds ratios (OR) for any DR and proliferative DR among men decreased significantly in the highest blood 25-hydroxyvitamin D level quintile relative to the lowest quintile (OR, 0.37; 95% confidence interval [CI], 0.18–0.76; P for trend  = 0.004 and OR, 0.15; 95% CI, 0.03–0.83; P for trend  = 0.043).

Conclusions

This study provides the first epidemiologic findings of the inverse relationships of blood 25-hydroxyvitamin D levels with any DR and proliferative DR only in men. Considering anti-angiogenic and anti-fibrotic action of vitamin D, further studies including longitudinal and interventional analysis are warranted.     

Pathway Analysis for Genome-Wide Association Study of Lung Cancer in Han Chinese Population

Genome-wide association studies (GWAS) have identified a number of genetic variants associated with lung cancer risk. However, these loci explain only a small fraction of lung cancer hereditability and other variants with weak effect may be lost in the GWAS approach due to the stringent significance level after multiple comparison correction. In this study, in order to identify important pathways involving the lung carcinogenesis, we performed a two-stage pathway analysis in GWAS of lung cancer in Han Chinese using gene set enrichment analysis (GSEA) method. Predefined pathways by BioCarta and KEGG databases were systematically evaluated on Nanjing study (Discovery stage: 1,473 cases and 1,962 controls) and the suggestive pathways were further to be validated in Beijing study (Replication stage: 858 cases and 1,115 controls). We found that four pathways (achPathway, metPathway, At1rPathway and rac1Pathway) were consistently significant in both studies and the P values for combined dataset were 0.012, 0.010, 0.022 and 0.005 respectively. These results were stable after sensitivity analysis based on gene definition and gene overlaps between pathways. These findings may provide new insights into the etiology of lung cancer.     

Replication Study in a Japanese Population of Six Susceptibility Loci for Type 2 Diabetes Originally Identified by a Transethnic Meta-Analysis of Genome-Wide Association Studies

AimWe performed a replication study in a Japanese population to evaluate the association between type 2 diabetes and six susceptibility loci (TMEM154, SSR1, FAF1, POU5F1, ARL15, and MPHOSPH9) originally identified by a transethnic meta-analysis of genome-wide association studies (GWAS) in 2014.MethodsWe genotyped 7,620 Japanese participants (5,817 type 2 diabetes patients and 1,803 controls) for each of the single nucleotide polymorphisms (SNPs) using a multiplex polymerase chain reaction invader assay. The association of each SNP locus with the disease was evaluated using logistic regression analysis.ResultsOf the six SNPs examined in this study, four (rs6813195 near TMEM154, rs17106184 in FAF1, rs3130501 in POU5F1 and rs4275659 near MPHOSPH9) had the same direction of effect as in the original reports, but two (rs9505118 in SSR1 and rs702634 in ARL15) had the opposite direction of effect. Among these loci, rs3130501 and rs4275659 were nominally associated with type 2 diabetes (rs3130501; p = 0.017, odds ratio [OR] = 1.113, 95% confidence interval [CI] 1.019–1.215, rs4275659; p = 0.012, OR = 1.127, 95% CI 1.026–1.238, adjusted for sex, age and body mass index), but we did not observe a significant association with type 2 diabetes for any of the six evaluated SNP loci in our Japanese population.ConclusionsOur results indicate that effects of the six SNP loci identified in the transethnic GWAS meta-analysis are not major among the Japanese, although SNPs in POU5F1 and MPHOSPH9 loci may have some effect on susceptibility to type 2 diabetes in this population.     

Association of Complement C5 Gene Polymorphisms with Proliferative Diabetic Retinopathy of Type 2 Diabetes in a Chinese Han Population

Purpose

To investigate the association of C5 SNPs with proliferative diabetic retinopathy (PDR) of type 2 diabetes (T2D).

Methods

A total of four C5 SNPs including rs2269067, rs7040033, rs1017119 and rs7027797 were genotyped in 400 PDR patients with T2D (cases) and 600 non- proliferative diabetic retinopathy PDR (NPDR) with T2D patients (controls) by using PCR-RFLP method. mRNA expression was examined by real-time PCR. Cytokine production was detected by ELISA.

Results

The frequency of GG genotype of C5 rs2269067 was significantly increased in cases compared with controls (Pc = 3.4×10⁻⁵, OR = 1.87). And C5 mRNA expression was significantly increased in rs2269067 GG cases as compared with CG or CC cases (P = 0.003, P = 0.001, respectively). Moreover, the production of IL-6 was significantly increased in rs2269067 GG cases compared to CG cases or CC cases (P = 0.002, P = 0.001, respectively).

Conclusions

C5 rs2269067 GG genotype confers risk for PDR of T2D in Chinese han population and is associated with an elevated C5 mRNA expression and an increased IL-6 production.