Shadoo (Sprn) and prion disease incubation time in mice

Prion diseases are transmissible neurodegenerative disorders of mammalian species and include scrapie, bovine spongiform encephalopathy (BSE), and variant Creutzfeldt-Jakob disease (vCJD). The prion protein (PrP) plays a key role in the disease, with coding polymorphism in both human and mouse influencing disease susceptibility and incubation time, respectively. Other genes are also thought to be important and a plausible candidate is Sprn, which encodes the PrP-like protein Shadoo (Sho). Sho is expressed in the adult central nervous system and exhibits neuroprotective activity reminiscent of PrP in an in vitro assay. To investigate the role of Sprn in prion disease incubation time we sequenced the open reading frame (ORF) in a diverse panel of mice and saw little variation except in strains derived from wild-trapped mice. Sequencing the untranslated regions revealed polymorphisms that allowed us to carry out an association study of incubation period in the Northport heterogeneous stock of mice inoculated with Chandler/RML prions. We also examined the expression level of Sprn mRNA in the brains of normal and prion-infected mice and saw no correlation with either genotype or incubation time. We therefore conclude that Sprn does not play a major role in prion disease incubation time in these strains of mice.     

Investigation of mcp1 as a quantitative trait gene for prion disease incubation time in mouse

The genetic basis of prion disease incubation time is principally determined by polymorphisms in the prion protein gene, Prnp. However, it is now known that other genetic factors are important. Several quantitative trait loci (QTL) have been identified across the genome including a broad region of linkage on Mmu11. Monocyte chemoattractant protein 1 (MCP-1) maps to this region and has been associated with microglial activation and reduced survival in the ME7 mouse scrapie model of prion disease. We have identified 10 polymorphisms, 3 of which are nonsynonomous, in Mcp1 between "long" (CAST) and "short" (SJL or NZW) incubation-time mouse strains. Crosses between these strains and Mcp1(-/-) mice inoculated with the Chandler/RML mouse scrapie prion strain formed the basis of a quantitative complementation test. In these models loss of Mcp1 did not show an increase in incubation time suggesting that the effects of Mcp1 may be specific to the ME7 prion strain and that Mcp1 does not contribute to the QTL described on Mmu11.     

Possible Single-Nucleotide Polymorphism Loci Associated with Systemic Sclerosis Susceptibility: A Genetic Association Study in a Chinese Han Population

Objective

The aim of this study was to confirm the association of RHOB and FAM167A-BLK gene polymorphisms with susceptibility to systemic sclerosis (SSc) in a Chinese Han population.

Methods

A total of 248 SSc patients and 251 healthy controls of Chinese Han ethnicity, which visited the department of dermatology of Peking Union Medical College Hospital, were included in the study. Six selected single nucleotide polymorphisms (SNPs) in the RHOB and FAM167A-BLK regions were selected as markers and were genotyped using a MassARRAY system, which is based on the matrix-assisted laser desorption/ionization time of flight mass spectrometry technique.

Results

Three SNPs in the coding regions of the RHOB and FAM167A-BLK genes displayed an association with SSc: (1) rs1062292T, which is a newly discovered SNP in the RHOB gene (P = 0.03, odds ratio [OR] = 1.62, 95% confidence interval (CI) = 1.05–2.50), (2) rs2736340T (P = 0.03, OR = 1.39, 95%CI = 1.03–1.85), and (3) rs13277113A (P = 0.04, OR = 1.34, 95%CI = 1.01–1.76), both in the FAM167A-BLK gene. Our results support previous findings that vaiants in the RHOB and FAM167A-BLK genes may be associated with susceptibility to SSc. However, the loci of the SNPs in RHOB region that displayed an association with SSc are quite different from the loci which were identified in studies of Caucasian populations.

Conclusion

Our results confirm that RHOB and FAM167A-BLK polymorphisms exist in Chinese Han SSc patients. Therefore, variants of the RHOB and FAM167A-BLK genes are promising genetic markers for SSc.     

Oxidative Stress Induced Age Dependent Meibomian Gland Dysfunction in Cu,Zn-Superoxide Dismutase-1 (Sod1) Knockout Mice

Purpose

The purpose of our study was to investigate alterations in the meibomian gland (MG) in Cu, Zn-Superoxide Dismutase-1 knockout (Sod1 ^−/−) mouse.

Methods

Tear function tests [Break up time (BUT) and cotton thread] and ocular vital staining test were performed on Sod1 ^−/− male mice (n = 24) aged 10 and 50 weeks, and age and sex matched wild–type (+/+) mice (n = 25). Tear and serum samples were collected at sacrifice for inflammatory cytokine assays. MG specimens underwent Hematoxylin and Eosin staining, Mallory staining for fibrosis, Oil Red O lipid staining, TUNEL staining, immunohistochemistry stainings for 4HNE, 8-OHdG and CD45. Transmission electron microscopic examination (TEM) was also performed.

Results

Corneal vital staining scores in the Sod1 ^−/− mice were significantly higher compared with the wild type mice throughout the follow-up. Tear and serum IL-6 and TNF-α levels also showed significant elevations in the 10 to 50 week Sod1 ^−/− mice. Oil Red O staining showed an accumulation of large lipid droplets in the Sod1 ^−/− mice at 50 weeks. Immunohistochemistry revealed both increased TUNEL and oxidative stress marker stainings of the MG acinar epithelium in the Sod1 ^−/− mice compared to the wild type mice. Immunohistochemistry staining for CD45 showed increasing inflammatory cell infiltrates from 10 to 50 weeks in the Sod1 ^−/− mice compared to the wild type mice. TEM revealed prominent mitochondrial changes in 50 week Sod1 ^−/− mice.

Conclusions

Our results suggest that reactive oxygen species might play a vital role in the pathogensis of meibomian gland dysfunction. The Sod1 ^−/− mouse appears to be a promising model for the study of reactive oxygen species associated MG alterations.     

SORL1 Is Genetically Associated with Late-Onset Alzheimer’s Disease in Japanese,Koreans and Caucasians

To discover susceptibility genes of late-onset Alzheimer’s disease (LOAD), we conducted a 3-stage genome-wide association study (GWAS) using three populations: Japanese from the Japanese Genetic Consortium for Alzheimer Disease (JGSCAD), Koreans, and Caucasians from the Alzheimer Disease Genetic Consortium (ADGC). In Stage 1, we evaluated data for 5,877,918 genotyped and imputed SNPs in Japanese cases (n = 1,008) and controls (n = 1,016). Genome-wide significance was observed with 12 SNPs in the APOE region. Seven SNPs from other distinct regions with p-values <2×10⁻⁵ were genotyped in a second Japanese sample (885 cases, 985 controls), and evidence of association was confirmed for one SORL1 SNP (rs3781834, P = 7.33×10⁻⁷ in the combined sample). Subsequent analysis combining results for several SORL1 SNPs in the Japanese, Korean (339 cases, 1,129 controls) and Caucasians (11,840 AD cases, 10,931 controls) revealed genome wide significance with rs11218343 (P = 1.77×10⁻⁹) and rs3781834 (P = 1.04×10⁻⁸). SNPs in previously established AD loci in Caucasians showed strong evidence of association in Japanese including rs3851179 near PICALM (P = 1.71×10⁻⁵) and rs744373 near BIN1 (P = 1.39×10⁻⁴). The associated allele for each of these SNPs was the same as in Caucasians. These data demonstrate for the first time genome-wide significance of LOAD with SORL1 and confirm the role of other known loci for LOAD in Japanese. Our study highlights the importance of examining associations in multiple ethnic populations.     

Case-Control Approach to Identify Plasmodium falciparum Polymorphisms Associated with Severe Malaria

Background

Studies to identify phenotypically-associated polymorphisms in the Plasmodium falciparum 23 Mb genome will require a dense array of marker loci. It was considered promising to undertake initial allelic association studies to prospect for virulence polymorphisms in Thailand, as the low endemicity would allow higher levels of linkage disequilibrium (LD) than would exist in more highly endemic areas.

Methodology/Principal Findings

Assessment of LD was first made with 11 microsatellite loci widely dispersed in the parasite genome, and 16 microsatellite loci covering a ∼140 kb region of chromosome 2 (an arbitrarily representative non-telomeric part of the genome), in a sample of 100 P. falciparum isolates. The dispersed loci showed minimal LD (Index of Association, I^S_A = 0.013, P = 0.10), while those on chromosome 2 showed significant LD values mostly between loci <5 kb apart. A disease association study was then performed comparing parasites in 113 severe malaria cases and 245 mild malaria controls. Genotyping was performed on almost all polymorphisms in the binding domains of three erythrocyte binding antigens (eba175, eba140 and eba181), and repeat sequence polymorphisms ∼2 kb apart in each of three reticulocyte binding homologues (Rh1, Rh2a/b, and Rh4). Differences between cases and controls were seen for (i) codons 388-90 in eba175, and (ii) a repeat sequence centred on Rh1 codon 667.

Conclusions/Significance

Allelic association studies on P. falciparum require dense genotypic markers, even in a population of only moderate endemicity that has more extensive LD than highly endemic populations. Disease-associated polymorphisms in the eba175 and Rh1 genes encode differences in the middle of previously characterised erythrocyte binding domains, marking these for further investigation.     

DNA Methyl Transferase (DNMT) Gene Polymorphisms Could Be a Primary Event in Epigenetic Susceptibility to Schizophrenia

DNA methylation has been implicated in the etiopathology of various complex disorders. DNA methyltransferases are involved in maintaining and establishing new methylation patterns. The aim of the present study was to investigate the inherent genetic variations within DNA methyltransferase genes in predisposing to susceptibility to schizophrenia. We screened for polymorphisms in DNA methyltransferases, DNMT1, DNMT3A, DNMT3B and DNMT3L in 330 schizophrenia patients and 302 healthy controls for association with Schizophrenia in south Indian population. These polymorphisms were also tested for subgroup analysis with patient''s gender, age of onset and family history. DNMT1 rs2114724 (genotype P = .004, allele P = 0.022) and rs2228611 (genotype P = 0.004, allele P = 0.022) were found to be significantly associated at genotypic and allelic level with Schizophrenia in South Indian population. DNMT3B rs2424932 genotype (P = 0.023) and allele (P = 0.0063) increased the risk of developing schizophrenia in males but not in females. DNMT3B rs1569686 (genotype P = 0.027, allele P = 0.033) was found to be associated with early onset of schizophrenia and also with family history and early onset (genotype P = 0.009). DNMT3L rs2070565 (genotype P = 0.007, allele P = 0.0026) confers an increased risk of developing schizophrenia at an early age in individuals with family history. In-silico prediction indicated functional relevance of these SNPs in regulating the gene. These observations might be crucial in addressing and understanding the genetic control of methylation level differences from ethnic viewpoint. Functional significance of genotype variations within the DNMTs indeed suggest that the genetic nature of methyltransferases should be considered while addressing epigenetic events mediated by methylation in Schizophrenia.     

Effects of Polymorphisms -1112C/T and +2044A/G in Interleukin-13 Gene on Asthma Risk: A Meta-Analysis

Background

Associations between interleukin-13 (IL-13) polymorphisms and asthma risk remained controversial and ambiguous. Therefore, we performed a meta-analysis to assess the associations between IL-13 polymorphisms and asthma susceptibility.

Methods

Pubmed, EMBASE, Chinese National Knowledge Infrastructure (CNKI) and Wangfang databases were searched. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to calculate the strength of association in the random-effects model.

Results

Thirty-four studies were included in this meta-analysis. The results indicated that IL13 -1112C/T polymorphism was significantly associated with asthma risk (OR = 1.20, 95% CI 1.08–1.34, P = 0.0009) in a dominant genetic model. When stratifying for race, IL13 -1112C/T polymorphism exhibited increased asthma risk in Caucasians (OR = 1.30, 95% CI 1.09–1.55, P = 0.003), while no significant association was found in Asians and African Americans. In the subgroup analysis based on atopic status, significant association was observed in atopic patients (OR = 1.25, 95% CI 1.07–1.45, P = 0.004) but not in the non-atopic patients. In addition, a significant association between IL13+2044A/G polymorphism and asthma risk was observed (OR = 1.18, 95% CI 1.08–1.28, P = 0.0002). In the subgroup analysis by ethnicity, there were significant associations between IL13+2044A/G polymorphism and asthma risk in Asians (OR = 1.19, 95% CI 1.04–1.36, P = 0.01) and Caucasians (OR = 1.22, 95% CI 1.06–1.40, P = 0.005) but not in African Americans. In the subgroup analysis stratified by atopic status, a marginal significant association was found in atopic patients (OR = 1.12, 95% CI 1.00–1.26, P = 0.05).

Conclusions

This meta-analysis suggested that the IL13 -1112C/T and +2044A/G polymorphisms were risk factors for asthma.     

Common Genetic Variation Near the Phospholamban Gene Is Associated with Cardiac Repolarisation: Meta-Analysis of Three Genome-Wide Association Studies

To identify loci affecting the electrocardiographic QT interval, a measure of cardiac repolarisation associated with risk of ventricular arrhythmias and sudden cardiac death, we conducted a meta-analysis of three genome-wide association studies (GWAS) including 3,558 subjects from the TwinsUK and BRIGHT cohorts in the UK and the DCCT/EDIC cohort from North America. Five loci were significantly associated with QT interval at P<1×10⁻⁶. To validate these findings we performed an in silico comparison with data from two QT consortia: QTSCD (n = 15,842) and QTGEN (n = 13,685). Analysis confirmed the association between common variants near NOS1AP (P = 1.4×10⁻⁸³) and the phospholamban (PLN) gene (P = 1.9×10⁻²⁹). The most associated SNP near NOS1AP (rs12143842) explains 0.82% variance; the SNP near PLN (rs11153730) explains 0.74% variance of QT interval duration. We found no evidence for interaction between these two SNPs (P = 0.99). PLN is a key regulator of cardiac diastolic function and is involved in regulating intracellular calcium cycling, it has only recently been identified as a susceptibility locus for QT interval. These data offer further mechanistic insights into genetic influence on the QT interval which may predispose to life threatening arrhythmias and sudden cardiac death.     

The Joint Effect of hOGG1, APE1, and ADPRT Polymorphisms and Cooking Oil Fumes on the Risk of Lung Adenocarcinoma in Chinese Non-Smoking Females

Background

The human 8-oxoguanine DNA glycosylase 1 (hOGG1), apurinic/apyrimidinic endonuclease 1 (APE1), and adenosine diphosphate ribosyl transferase (ADPRT) genes play an important role in the DNA base excision repair pathway. Single nucleotide polymorphisms (SNPs) in critical genes are suspected to be associated with the risk of lung cancer. This study aimed to identify the association between the polymorphisms of hOGG1 Ser326Cys, APE1 Asp148Glu, and ADPRT Val762Ala, and the risk of lung adenocarcinoma in the non-smoking female population, and investigated the interaction between genetic polymorphisms and environmental exposure in lung adenocarcinoma.

Methods

We performed a hospital-based case-control study, including 410 lung adenocarcinoma patients and 410 cancer-free hospital control subjects who were matched for age. Each case and control was interviewed to collect information by well-trained interviewers. A total of 10 ml of venous blood was collected for genotype testing. Three polymorphisms were analyzed by the polymerase chain reaction-restriction fragment length polymorphism technique.

Results

We found that individuals who were homozygous for the variant hOGG1 326Cys/Cys showed a significantly increased risk of lung adenocarcinoma (OR = 1.54; 95% CI: 1.01–2.36; P = 0.045). When the combined effect of variant alleles was analyzed, we found an increased OR of 1.89 (95% CI: 1.24–2.88, P = 0.003) for lung adenocarcinoma individuals with more than one homozygous variant allele. In stratified analyses, we found that the OR for the gene-environment interaction between Ser/Cys and Cys/Cys genotypes of hOGG1 codon 326 and cooking oil fumes for the risk of lung adenocarcinoma was 1.37 (95% CI: 0.77–2.44; P = 0.279) and 2.79 (95% CI: 1.50–5.18; P = 0.001), respectively.

Conclusions

The hOGG1 Ser326Cys polymorphism might be associated with the risk of lung adenocarcinoma in Chinese non-smoking females. Furthermore, there is a significant gene-environment association between cooking oil fumes and hOGG1 326 Cys/Cys genotype in lung adenocarcinoma among female non-smokers.     

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.     

Sequence Analysis of Six Blood Pressure Candidate Regions in 4,178 Individuals: The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Targeted Sequencing Study

Background

Genome-wide association studies (GWAS) identified multiple loci for blood pressure (BP) and hypertension. Six genes – ATP2B1, CACNB2, CYP17A1, JAG1, PLEKHA7, and SH2B3 – were evaluated for sequence variation with large effects on systolic blood pressure (SBP), diastolic blood pressure (DBP), pulse pressure (PP), and mean arterial pressure (MAP).

Methods and Results

Targeted genomic sequence was determined in 4,178 European ancestry participants from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium. Common variants (≥50 minor allele copies) were evaluated individually and rare variants (minor allele frequency, MAF≤1%) were aggregated by locus. 464 common variants were identified across the 6 genes. An upstream CYP17A1 variant, rs11191416 (MAF = 0.09), was the most significant association for SBP (P = 0.0005); however the association was attenuated (P = 0.0469) after conditioning on the GWAS index single nucleotide polymorphism (SNP). A PLEKHA7 intronic variant was the strongest DBP association (rs12806040, MAF = 0.007, P = 0.0006) and was not in LD (r² = 0.01) with the GWAS SNP. A CACNB2 intronic SNP, rs1571787, was the most significant association with PP (MAF = 0.27, P = 0.0003), but was not independent from the GWAS SNP (r² = 0.34). Three variants (rs6163 and rs743572 in the CYP17A1 region and rs112467382 in PLEKHA7) were associated with BP traits (P<0.001). Rare variation, aggregately assessed in the 6 regions, was not significantly associated with BP measures.

Conclusion

Six targeted gene regions, previously identified by GWAS, did not harbor novel variation with large effects on BP in this sample.     

Genome-Wide Association Study in BRCA1 Mutation Carriers Identifies Novel Loci Associated with Breast and Ovarian Cancer Risk

BRCA1-associated breast and ovarian cancer risks can be modified by common genetic variants. To identify further cancer risk-modifying loci, we performed a multi-stage GWAS of 11,705 BRCA1 carriers (of whom 5,920 were diagnosed with breast and 1,839 were diagnosed with ovarian cancer), with a further replication in an additional sample of 2,646 BRCA1 carriers. We identified a novel breast cancer risk modifier locus at 1q32 for BRCA1 carriers (rs2290854, P = 2.7×10⁻⁸, HR = 1.14, 95% CI: 1.09–1.20). In addition, we identified two novel ovarian cancer risk modifier loci: 17q21.31 (rs17631303, P = 1.4×10⁻⁸, HR = 1.27, 95% CI: 1.17–1.38) and 4q32.3 (rs4691139, P = 3.4×10⁻⁸, HR = 1.20, 95% CI: 1.17–1.38). The 4q32.3 locus was not associated with ovarian cancer risk in the general population or BRCA2 carriers, suggesting a BRCA1-specific association. The 17q21.31 locus was also associated with ovarian cancer risk in 8,211 BRCA2 carriers (P = 2×10⁻⁴). These loci may lead to an improved understanding of the etiology of breast and ovarian tumors in BRCA1 carriers. Based on the joint distribution of the known BRCA1 breast cancer risk-modifying loci, we estimated that the breast cancer lifetime risks for the 5% of BRCA1 carriers at lowest risk are 28%–50% compared to 81%–100% for the 5% at highest risk. Similarly, based on the known ovarian cancer risk-modifying loci, the 5% of BRCA1 carriers at lowest risk have an estimated lifetime risk of developing ovarian cancer of 28% or lower, whereas the 5% at highest risk will have a risk of 63% or higher. Such differences in risk may have important implications for risk prediction and clinical management for BRCA1 carriers.     

Functional Polymorphisms of FAS and FASL Gene and Risk of Breast Cancer – Pilot Study of 134 Cases

Fas/Fas ligand (FasL) system is one of the key apoptotic signaling entities in the extrinsic apoptotic pathway. De-regulation of this pathway, i.e. by mutations may prevent the immune system from the removal of newly-formed tumor cells, and thus lead to tumor formation. The present study investigated the association between −1377 G/A (rs2234767) and −670 A/G (rs1800682) polymorphisms in Fas as well as single nucleotide polymorphisms INV2nt −124 A/G (rs5030772) and −844 C/T (rs763110) in FasL in a sample of Iranian patients with breast cancer. This case-control study was done on 134 breast cancer patients and 152 normal women. Genomic DNA was extracted from whole blood samples. The polymorphisms were determined by using tetra-ARMS-PCR method. There was no significant difference in the genotype distribution of FAS rs2234767 polymorphism between cases and controls. FAS rs1800682, FASL rs5030772, and FASL rs763110 genotypes showed significant associations with an increasing risk of breast cancer (odds ratio OR = 3.18, P = 0.019; OR = 5.08, P = 0.012; OR = 2.40, P = 0.024, respectively). In conclusion, FAS rs2234767 was not associated with breast cancer risk. Though, FAS rs1800682, FASL rs5030772, and FASL rs763110 polymorphisms were associated with the risk of breast cancer in the examined population.     

Association of Mitochondrial DNA Polymerase γ Gene POLG1 Polymorphisms with Parkinsonism in Chinese Populations

Background

Mitochondrial DNA polymerase gamma (POLG1) mutations were associated with levodopa-responsive Parkinsonism. POLG1 gene contains a number of common nonsynonymous SNPs and intronic regulatory SNPs which may have functional consequences. It is of great interest to discover polymorphisms variants associated with Parkinson''s disease (PD), both in isolation and in combination with specific SNPs.

Materials and Methods

We conducted a case-control study and genotyped twenty SNPs and poly-Q polymorphisms of POLG1 gene including in 344 Chinese sporadic PD patients and 154 healthy controls. All the polymorphisms of POLG1 we found in this study were sequenced by PCR products with dye terminator methods using an ABI-3100 sequencer. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) for association between twenty POLG1 SNPs and PD were calculated using the program Haploview.

Principal Results

We provided evidence for strong association of four intronic SNPs of the POLG1 gene (new report: c.2070-12T>A and rs2307439: c.2070-64G>A in intron 11, P = 0.00011, OR = 1.727; rs2302084: c.3105-11T>C and rs2246900: c.3105-36A>G in intron 19, P = 0.00031, OR = 1.648) with PD. However, we did not identify any significant association between ten exonic SNPs of POLG1 and PD. Linkage disequilibrium analysis indicated that c.2070-12T>A and c.2070-64G>A could be parsed into one block as Haplotype 1 as well as c.3105-11T>C and c.3105-36A>G in Haplotype 2. In addition, case and control study on association of POLG1 CAG repeat (poly-Q) alleles with PD showed a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD. Although intronic SNPs associated with PD didn''t influence POLG1 mRNA alternative splicing, there was a strong association of c.2070-12T>A and c.2070-64G>A with decreased POLG1 mRNA level and protein levels.

Conclusions

Our findings indicate that POLG1 may play a role in the pathogenesis of PD in Chinese populations.     

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.