COPD; Genetics; Association analysis; Consortium
The pulmonary function measures of forced expiratory volume in one second (FEV1) and its ratio to forced vital capacity (FVC) are used in the diagnosis and monitoring of lung diseases and predict cardiovascular mortality in the general population. Genome wide association studies (GWAS) have identified numerous loci associated with FEV1 and FEV1/FVC but the causal variants remain uncertain. We hypothesized that novel or rare variants poorly tagged by GWAS may explain the significant associations between FEV1/FVC and two genes: ADAM19 and HTR4.
Methods and Results
We sequenced ADAM19 and its promoter region along with the approximately 21 kb portion of HTR4 harboring GWAS SNPs for pulmonary function and analyzed associations with FEV1/FVC among 3,983 participants of European ancestry from Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE). Meta-analysis of common variants in each region identified statistically significant associations (316 tests, P < 1.58×10−4) with FEV1/FVC for 14 ADAM19 SNPs and 24 HTR4 SNPs. After conditioning on the sentinel GWAS hit in each gene [ADAM19 rs1422795, minor allele frequency (MAF)=0.33 and HTR4 rs11168048, MAF=0.40] one SNP remained statistically significant (ADAM19 rs13155908, MAF = 0.12, P = 1.56×10−4). Analysis of rare variants (MAF < 1%) using Sequence Kernel Association Test did not identify associations with either region.
Sequencing identified one common variant associated with FEV1/FVC independently of the sentinel ADAM19 GWAS hit and supports the original HTR4 GWAS findings. Rare variants do not appear to underlie GWAS associations with pulmonary function for common variants in ADAM19 and HTR4.
genetic polymorphism; lung; population studies; DNA sequencing; Genome Wide Association Study
Forced vital capacity (FVC), a spirometric measure of pulmonary function, reflects lung volume and is used to diagnose and monitor lung diseases. We performed genome-wide association study meta-analysis of FVC in 52,253 individuals from 26 studies and followed up the top associations in 32,917 additional individuals of European ancestry. We found six new regions associated at genome-wide significance (P < 5 × 10−8) with FVC in or near EFEMP1, BMP6, MIR-129-2/HSD17B12, PRDM11, WWOX, and KCNJ2. Two (GSTCD and PTCH1) loci previously associated with spirometric measures were related to FVC. Newly implicated regions were followed-up in samples of African American, Korean, Chinese, and Hispanic individuals. We detected transcripts for all six newly implicated genes in human lung tissue. The new loci may inform mechanisms involved in lung development and pathogenesis of restrictive lung disease.
We examined common variants in the fatty acid binding protein 4 gene (FABP4) and plasma levels of FABP4 in adults aged 65 and older from the Cardiovascular Health Study. We genotyped rs16909187, rs1054135, rs16909192, rs10808846, rs7018409, rs2290201, and rs6992708 and measured circulating FABP4 levels among 3190 European Americans and 660 African Americans. Among European Americans, the minor alleles of six single nucleotide polymorphisms (SNP) were associated with lower FABP4 levels (all p ≤ 0.01). Among African Americans, the SNP with the lowest minor allele frequency was associated with lower FABP4 levels (p = 0.015). The C-A haplotype of rs16909192 and rs2290201 was associated with lower FABP4 levels in both European Americans (frequency = 16 %; p = 0.001) and African Americans (frequency = 8 %; p = 0.04). The haplotype combined a SNP in the first intron with one in the 3′untranslated region. However, the alleles associated with lower FABP4 levels were associated with higher fasting glucose in meta-analyses from the MAGIC consortium. These results demonstrate associations of common SNP and haplotypes in the FABP4 gene with lower plasma FABP4 but higher fasting glucose levels.
Fatty acid binding proteins; Metabolism; Genetics
Genome-wide association studies (GWAS) have identified numerous loci influencing cross-sectional lung function, but less is known about genes influencing longitudinal change in lung function.
We performed GWAS of the rate of change in forced expiratory volume in the first second (FEV1) in 14 longitudinal, population-based cohort studies comprising 27,249 adults of European ancestry using linear mixed effects model and combined cohort-specific results using fixed effect meta-analysis to identify novel genetic loci associated with longitudinal change in lung function. Gene expression analyses were subsequently performed for identified genetic loci. As a secondary aim, we estimated the mean rate of decline in FEV1 by smoking pattern, irrespective of genotypes, across these 14 studies using meta-analysis.
The overall meta-analysis produced suggestive evidence for association at the novel IL16/STARD5/TMC3 locus on chromosome 15 (P = 5.71 × 10-7). In addition, meta-analysis using the five cohorts with ≥3 FEV1 measurements per participant identified the novel ME3 locus on chromosome 11 (P = 2.18 × 10-8) at genome-wide significance. Neither locus was associated with FEV1 decline in two additional cohort studies. We confirmed gene expression of IL16, STARD5, and ME3 in multiple lung tissues. Publicly available microarray data confirmed differential expression of all three genes in lung samples from COPD patients compared with controls. Irrespective of genotypes, the combined estimate for FEV1 decline was 26.9, 29.2 and 35.7 mL/year in never, former, and persistent smokers, respectively.
In this large-scale GWAS, we identified two novel genetic loci in association with the rate of change in FEV1 that harbor candidate genes with biologically plausible functional links to lung function.
Chronic obstructive pulmonary disease (COPD) is linked to cardiovascular disease; however, there are few studies on the associations of cardiovascular genes with COPD.
We assessed the association of lung function with 2,100 genes selected for cardiovascular diseases among 20,077 European-Americans and 6,900 African-Americans. We performed replication of significant loci in the other racial group and an independent consortium of Europeans, tested the associations of significant loci with percent emphysema, and examined gene expression in an independent sample. We then tested the association of a related lipid biomarker with FEV1/FVC and percent emphysema.
We identified one new polymorphism for FEV1/FVC (rs805301) in European-Americans (p=1.3×10−6) and a second (rs707974) in the combined European-American and African-American analysis (p=1.38×10−7). Both SNPs flank the gene for apolipoprotein M (apoM), a component of HDL. Both replicated in an independent cohort. SNPs in a second gene related to apoM and HDL, PCSK9, were associated with FEV1/FVC among African-Americans. rs707974 was associated with percent emphysema among European-Americans and African-Americans, and APOM expression was related to FEV1/FVC and percent emphysema. Higher HDL levels were associated with lower FEV1/FVC and greater percent emphysema.
These findings suggest a novel role for the APOM/HDL pathway in the pathogenesis of COPD and emphysema.
Apolipoproteins; Cholesterol; Percent Emphysema; Polymorphism, Single Nucleotide; Pulmonary Disease, Chronic Obstructive
Intake of marine- based omega-3 fatty acids [eicosapentaenoic acids (EPA), docosapentaenoic acid (DPA), and docosahexaenoic acids (DHA)] is recommended to prevent coronary heart disease (CHD). Stearidonic acid (SDA), a plant- based omega-3 fatty acid (FA), is a precursor of EPA and may be more readily converted to EPA than alpha-linolenic acid (ALA). While transgenic soybean might supply SDA at low costs, it is unclear whether SDA is associated with CHD risk. Furthermore, associations of other omega-3 FAs with CHD risk remain inconsistent. This ancillary study examined the association of red blood cell SDA as well as other omega-3 FAs with the risk of CHD. In a prospective nested case-control study of the Physicians’ Health Study, we randomly selected 1,000 pairs of incident CHD with matching controls. Red blood cell FAs were measured using gas chromatography. We used conditional logistic regression to estimate relative risks.
Mean age was 68.7±8.7 y. In a multivariable model controlling for matching factors and established CHD risk factors, odds ratio for CHD for each standard deviation increase of log-SDA was 1.03 (95% CI: 0.90, 1.18). Corresponding values for log ALA and log-marine omega-3 FAs were 1.04 (95% CI: 0.94, 1.16) and 0.97 (95% CI: 0.88, 1.07), respectively.
In conclusion, our data did not show an association between red blood cell SDA, ALA, or marine omega-3 FAs and the risk of CHD in male physicians.
Stearidonic acid (SDA); omega-3 fatty acids; alpha-linolenic acid (ALA); coronary heart disease
Vitamin D is associated with lung health in epidemiologic studies, but mechanisms mediating observed associations are poorly understood. This study explores mechanisms for an effect of vitamin D in lung through an in vivo gene expression study, an expression quantitative trait loci (eQTL) analysis in lung tissue, and a population-based cohort study of sequence variants.
Microarray analysis investigated the association of gene expression in small airway epithelial cells with serum 25(OH)D in adult non-smokers. Sequence variants in candidate genes identified by the microarray were investigated in a lung tissue eQTL database, and also in relation to cross-sectional pulmonary function in the Health, Aging, and Body Composition (Health ABC) study, stratified by race, with replication in the Framingham Heart Study (FHS).
13 candidate genes had significant differences in expression by serum 25(OH)D (nominal p < 0.05), and a genome-wide significant eQTL association was detected for SGPP2. In Health ABC, SGPP2 SNPs were associated with FEV1 in both European- and African-Americans, and the gene-level association was replicated in European-American FHS participants. SNPs in 5 additional candidate genes (DAPK1, FSTL1, KAL1, KCNS3, and RSAD2) were associated with FEV1 in Health ABC participants.
SGPP2, a sphingosine-1-phosphate phosphatase, is a novel vitamin D-responsive gene associated with lung function. The identified associations will need to be followed up in further studies.
Vitamin D; Airflow obstruction; FEV1; SGPP2; FEV1/FVC
Limited data are available on the association between adiponectin and incident heart failure. In the current ancillary study to the Physicians' Health Study, we used a prospective nested-case control design to examine whether plasma adiponectin concentration was related to the risk of heart failure. We selected 787 incident heart failure cases and 787 matched controls for the current analysis. Each control was selected using a risk set sampling technique at the time of the occurrence of the index case and matched on year of birth, age at blood collection, and race. Adiponectin was measured using ELISA. Heart failure occurrence was self-reported in annual follow-up questionnaire. Validation of self-reported heart failure in this cohort has been published. The mean age was 58.7 years. In a conditional logistic regression adjusting for age, race, time of blood collection, year of birth, hypertension, atrial fibrillation, smoking, alcohol intake, and exercise, estimates of the relative risk (95% confidence interval) were 1.0 (ref), 0.74 (0.53–1.04), 0.67 (0.48–0.94), 0.70 (0.50–0.99), and 0.92 (0.65–1.30) from the lowest to the highest quintile of adiponectin, respectively, p for quadratic trend 0.004. Additional adjustment for potential mediating factors including diabetes, C-reactive protein, and body mass index led to the attenuation of the estimate of effect [1.0 (ref), 0.81 (0.57–1.15), 0.75 (0.53–1.06), 0.83 (0.58–1.18), and 1.26 (0.87–1.81) across consecutive quintiles of adiponectin]. Our data are consistent with a J-shaped association between total adiponectin and the risk of heart failure among US male physicians.
Adiponectin; epidemiology; heart failure; risk factors
Genome-wide association studies (GWAS) have identified determinants of chronic obstructive pulmonary disease, asthma and lung function level, however none addressed decline in lung function.
We conducted the first GWAS on age-related decline in forced expiratory volume in the first second (FEV1) and in its ratio to forced vital capacity (FVC) stratified a priori by asthma status.
Discovery cohorts included adults of European ancestry (1441 asthmatics, 2677 non-asthmatics; Epidemiological Study on the Genetics and Environment of Asthma (EGEA); Swiss Cohort Study on Air Pollution And Lung And Heart Disease In Adults (SAPALDIA); European Community Respiratory Health Survey (ECRHS)). The associations of FEV1 and FEV1/FVC decline with 2.5 million single nucleotide polymorphisms (SNPs) were estimated. Thirty loci were followed-up by in silico replication (1160 asthmatics, 10858 non-asthmatics: Atherosclerosis Risk in Communities (ARIC); Framingham Heart Study (FHS); British 1958 Birth Cohort (B58C); Dutch asthma study).
Main signals identified differed between asthmatics and non-asthmatics. None of the SNPs reached genome-wide significance. The association between the height related gene DLEU7 and FEV1 decline suggested for non-asthmatics in the discovery phase was replicated (discovery P=4.8×10−6; replication P=0.03) and additional sensitivity analyses point to a relation to growth. The top ranking signal, TUSC3, associated with FEV1/FVC decline in asthmatics (P=5.3×10−8) did not replicate. SNPs previously associated with cross-sectional lung function were not prominently associated with decline.
Genetic heterogeneity of lung function may be extensive. Our results suggest that genetic determinants of longitudinal and cross-sectional lung function differ and vary by asthma status.
Asthma; cohort studies; genome-wide association; lung function decline; heterogeneity
Rationale: Genome-wide association studies (GWAS) have identified loci influencing lung function, but fewer genes influencing chronic obstructive pulmonary disease (COPD) are known.
Objectives: Perform meta-analyses of GWAS for airflow obstruction, a key pathophysiologic characteristic of COPD assessed by spirometry, in population-based cohorts examining all participants, ever smokers, never smokers, asthma-free participants, and more severe cases.
Methods: Fifteen cohorts were studied for discovery (3,368 affected; 29,507 unaffected), and a population-based family study and a meta-analysis of case-control studies were used for replication and regional follow-up (3,837 cases; 4,479 control subjects). Airflow obstruction was defined as FEV1 and its ratio to FVC (FEV1/FVC) both less than their respective lower limits of normal as determined by published reference equations.
Measurements and Main Results: The discovery meta-analyses identified one region on chromosome 15q25.1 meeting genome-wide significance in ever smokers that includes AGPHD1, IREB2, and CHRNA5/CHRNA3 genes. The region was also modestly associated among never smokers. Gene expression studies confirmed the presence of CHRNA5/3 in lung, airway smooth muscle, and bronchial epithelial cells. A single-nucleotide polymorphism in HTR4, a gene previously related to FEV1/FVC, achieved genome-wide statistical significance in combined meta-analysis. Top single-nucleotide polymorphisms in ADAM19, RARB, PPAP2B, and ADAMTS19 were nominally replicated in the COPD meta-analysis.
Conclusions: These results suggest an important role for the CHRNA5/3 region as a genetic risk factor for airflow obstruction that may be independent of smoking and implicate the HTR4 gene in the etiology of airflow obstruction.
chronic obstructive pulmonary disease; single-nucleotide polymorphism; genes
Atopy and plasma IgE concentration are genetically complex traits, and the specific genetic risk factors that lead to IgE dysregulation and clinical atopy are an area of active investigation.
To ascertain the genetic risk factors which lead to IgE dysregulation.
A genome wide association study (GWAS) was performed in 6,819 participants from the Framingham Heart Study (FHS). Seventy of the top SNPs were selected based on p-values and linkage disequilibrium among neighboring SNPs and evaluated in a meta-analysis with five independent populations from the KORA, B58C, and CAMP cohorts.
Thirteen SNPs located in the region of three genes, FCER1A, STAT6, and IL-13, were found to have genome-wide significance in the FHS GWAS. The most significant SNPs from the three regions were rs2251746 (FCER1A, p-value 2.11×10-12), rs1059513 (STAT6, p-value 2.87×10-08), and rs1295686 (IL-13, p-value 3.55×10-08). Four additional gene regions - HLA-G, HLA-DQA2, HLA-A, and DARC - reached genome-wide statistical significance in meta-analysis combining FHS and replication cohorts, although the DARC association did not appear independent of SNPs in the nearby FCER1A gene.
This GWAS of the FHS has identified genetic loci in HLA genes that may have a role in the pathogenesis of IgE dysregulation and atopy. It also confirmed the association of known susceptibility loci, FCER1A, STAT6, and IL-13, for the dysregulation of total IgE.
total IgE; atopy; asthma; GWAS
Genome-wide association studies have identified numerous genetic loci for spirometic measures of pulmonary function, forced expiratory volume in one second (FEV1), and its ratio to forced vital capacity (FEV1/FVC). Given that cigarette smoking adversely affects pulmonary function, we conducted genome-wide joint meta-analyses (JMA) of single nucleotide polymorphism (SNP) and SNP-by-smoking (ever-smoking or pack-years) associations on FEV1 and FEV1/FVC across 19 studies (total N = 50,047). We identified three novel loci not previously associated with pulmonary function. SNPs in or near DNER (smallest PJMA = 5.00×10−11), HLA-DQB1 and HLA-DQA2 (smallest PJMA = 4.35×10−9), and KCNJ2 and SOX9 (smallest PJMA = 1.28×10−8) were associated with FEV1/FVC or FEV1 in meta-analysis models including SNP main effects, smoking main effects, and SNP-by-smoking (ever-smoking or pack-years) interaction. The HLA region has been widely implicated for autoimmune and lung phenotypes, unlike the other novel loci, which have not been widely implicated. We evaluated DNER, KCNJ2, and SOX9 and found them to be expressed in human lung tissue. DNER and SOX9 further showed evidence of differential expression in human airway epithelium in smokers compared to non-smokers. Our findings demonstrated that joint testing of SNP and SNP-by-environment interaction identified novel loci associated with complex traits that are missed when considering only the genetic main effects.
Measures of pulmonary function provide important clinical tools for evaluating lung disease and its progression. Genome-wide association studies have identified numerous genetic risk factors for pulmonary function but have not considered interaction with cigarette smoking, which has consistently been shown to adversely impact pulmonary function. In over 50,000 study participants of European descent, we applied a recently developed joint meta-analysis method to simultaneously test associations of gene and gene-by-smoking interactions in relation to two major clinical measures of pulmonary function. Using this joint method to incorporate genetic main effects plus gene-by-smoking interaction, we identified three novel gene regions not previously related to pulmonary function: (1) DNER, (2) HLA-DQB1 and HLA-DQA2, and (3) KCNJ2 and SOX9. Expression analyses in human lung tissue from ours or prior studies indicate that these regions contain genes that are plausibly involved in pulmonary function. This work highlights the utility of employing novel methods for incorporating environmental interaction in genome-wide association studies to identify novel genetic regions.
Mutations in the leucine-rich repeat kinase 2 gene (LRRK2), located at 12q12, are the most common known genetic causes of Parkinson’s disease (PD). Studies of LRRK2 mutation carriers have shown incomplete and age-dependent penetrance and previous studies have suggested that inherited susceptibility factors may modify the penetrance of LRRK2 mutations.
Genomewide linkage to age of onset of LRRK2-related PD was evaluated in a sample of 113 LRRK2 mutation carriers from 64 families using single nucleotide polymorphism data from the Illumina HumanCNV370 genotyping array. Association between onset age and SNPs located under suggestive linkage peaks was also evaluated.
The top LOD-score for onset age (LOD-score=2.43) was located in the chromosome 1q32.1 region. Moderate linkage to onset was also identified at 16q12.1 (LOD-score=1.58). Examination of single nucleotide polymorphism association to PD onset under the linkage peaks revealed no statistically significant SNP associations.
The two novel genomic regions identified may harbor modifiers of LRRK2-related PD onset age or penetrance and further study of these regions may provide important insight into LRRK2-related PD.
Parkinson’s Disease; LRRK2; Linkage
Pulmonary function measures reflect respiratory health and predict mortality, and are used in the diagnosis of chronic obstructive pulmonary disease (COPD). We tested genome-wide association with the forced expiratory volume in 1 second (FEV1) and the ratio of FEV1 to forced vital capacity (FVC) in 48,201 individuals of European ancestry, with follow-up of top associations in up to an additional 46,411 individuals. We identified new regions showing association (combined P<5×10−8) with pulmonary function, in or near MFAP2, TGFB2, HDAC4, RARB, MECOM (EVI1), SPATA9, ARMC2, NCR3, ZKSCAN3, CDC123, C10orf11, LRP1, CCDC38, MMP15, CFDP1, and KCNE2. Identification of these 16 new loci may provide insight into the molecular mechanisms regulating pulmonary function and into molecular targets for future therapy to alleviate reduced lung function.
Although family history is a well-established risk factor for Parkinson's disease (PD), fewer than 5% of PD cases can be attributed to known genetic mutations. The etiology for the remainder of PD cases is unclear; however, neuronal accumulation of the protein α-synuclein is common to nearly all patients, implicating pathways that influence α-synuclein in PD pathogenesis. We report a genome-wide significant association (P = 3.97 × 10−8) between a polymorphism, rs1564282, in the cyclin-G-associated kinase (GAK) gene and increased PD risk, with a meta-analysis odds ratio of 1.48. This association result is based on the meta-analysis of three publicly available PD case–control genome-wide association study and genotyping from a new, independent Italian cohort. Microarray expression analysis of post-mortem frontal cortex from PD and control brains demonstrates a significant association between rs1564282 and higher α-synuclein expression, a known cause of early onset PD. Functional knockdown of GAK in cell culture causes a significant increase in toxicity when α-synuclein is over-expressed. Furthermore, knockdown of GAK in rat primary neurons expressing the A53T mutation of α-synuclein, a well-established model for PD, decreases cell viability. These observations provide evidence that GAK is associated with PD risk and suggest that GAK and α-synuclein interact in a pathway involved in PD pathogenesis. The GAK protein, a serine/threonine kinase, belongs to a family of proteins commonly targeted for drug development. This, combined with GAK's observed relationship to the levels of α-synuclein expression and toxicity, suggests that the protein is an attractive therapeutic target for the treatment of PD.
Duplications and triplications of the α-synuclein (SNCA) gene increase risk for PD, suggesting increased expression levels of the gene to be associated with increased PD risk. However, past SNCA expression studies in brain tissue report inconsistent results. We examined expression of the full-length SNCA transcript (140 amino acid protein isoform), as well as total SNCA mRNA levels in 165 frontal cortex samples (101 PD, 64 control) using quantitative real-time polymerase chain reaction. Additionally, we evaluated the relationship of eight SNPs in both 5′ and 3′ regions of SNCA with the gene expression levels. The association between postmortem interval (PMI) and SNCA expression was different for PD and control samples: SNCA expression decreased with increasing PMI in cases, while staying relatively constant in controls. For short PMI, SNCA expression was increased in PD relative to control samples, whereas for long PMI, SNCA expression in PD was decreased relative to control samples.
Like most complex phenotypes, exceptional longevity is thought to reflect a combined influence of environmental (e.g., lifestyle choices, where we live) and genetic factors. To explore the genetic contribution, we undertook a genome-wide association study of exceptional longevity in 801 centenarians (median age at death 104 years) and 914 genetically matched healthy controls. Using these data, we built a genetic model that includes 281 single nucleotide polymorphisms (SNPs) and discriminated between cases and controls of the discovery set with 89% sensitivity and specificity, and with 58% specificity and 60% sensitivity in an independent cohort of 341 controls and 253 genetically matched nonagenarians and centenarians (median age 100 years). Consistent with the hypothesis that the genetic contribution is largest with the oldest ages, the sensitivity of the model increased in the independent cohort with older and older ages (71% to classify subjects with an age at death>102 and 85% to classify subjects with an age at death>105). For further validation, we applied the model to an additional, unmatched 60 centenarians (median age 107 years) resulting in 78% sensitivity, and 2863 unmatched controls with 61% specificity. The 281 SNPs include the SNP rs2075650 in TOMM40/APOE that reached irrefutable genome wide significance (posterior probability of association = 1) and replicated in the independent cohort. Removal of this SNP from the model reduced the accuracy by only 1%. Further in-silico analysis suggests that 90% of centenarians can be grouped into clusters characterized by different “genetic signatures” of varying predictive values for exceptional longevity. The correlation between 3 signatures and 3 different life spans was replicated in the combined replication sets. The different signatures may help dissect this complex phenotype into sub-phenotypes of exceptional longevity.
African Americans have been understudied in genome wide association studies of diabetes and related traits. In the current study, we examined the joint association of single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) with fasting insulin and an index of insulin resistance (HOMA-IR) in the HyperGEN study, a family based study with proband ascertainment for hypertension. This analysis is restricted to 1,040 African Americans without diabetes. We generated allele specific CNV genotypes at 872,243 autosomal loci using Birdsuite, a freely available multi-stage program. Joint tests of association for SNPs and CNVs were performed using linear mixed models adjusting for covariates and familial relationships. Our results highlight SNPs associated with fasting insulin and HOMA-IR (rs6576507 and rs8026527, 3.7*10−7≤P≤1.1*10−5) near ATPase, class V, type 10A (ATP10A), and the L Type voltage dependent calcium channel (CACNA1D, rs1401492, P≤5.2*10−6). ATP10A belongs to a family of aminophospholipid-transporting ATPases and has been associated with type 2 diabetes in mice. CACNA1D has been linked to pancreatic beta cell generation in mice. The two most significant copy variable markers (rs10277702 and rs361367; P<2.0*10−4) were in the beta variable region of the T-cell receptor gene (TCRVB). Human and mouse TCR has been shown to mimic insulin and its receptor and could contribute to insulin resistance. Our findings differ from genome wide association studies of fasting insulin and other diabetes related traits in European populations, highlighting the continued need to investigate unique genetic influences for understudied populations such as African Americans.
Copy number variants (CNVs) are known to cause Mendelian forms of Parkinson disease (PD), most notably in SNCA and PARK2. PARK2 has a recessive mode of inheritance; however, recent evidence demonstrates that a single CNV in PARK2 (but not a single missense mutation) may increase risk for PD. We recently performed a genome-wide association study for PD that excluded individuals known to have either a LRRK2 mutation or two PARK2 mutations. Data from the Illumina370Duo arrays were re-clustered using only white individuals with high quality intensity data, and CNV calls were made using two algorithms, PennCNV and QuantiSNP. After quality assessment, the final sample included 816 cases and 856 controls. Results varied between the two CNV calling algorithms for many regions, including the PARK2 locus (genome-wide p = 0.04 for PennCNV and p = 0.13 for QuantiSNP). However, there was consistent evidence with both algorithms for two novel genes, USP32 and DOCK5 (empirical, genome-wide p-values<0.001). PARK2 CNVs tended to be larger, and all instances that were molecularly tested were validated. In contrast, the CNVs in both novel loci were smaller and failed to replicate using real-time PCR, MLPA, and gel electrophoresis. The DOCK5 variation is more akin to a VNTR than a typical CNV and the association is likely caused by artifact due to DNA source. DNA for all the cases was derived from whole blood, while the DNA for all controls was derived from lymphoblast cell lines. The USP32 locus contains many SNPs with low minor allele frequency leading to a loss of heterozygosity that may have been spuriously interpreted by the CNV calling algorithms as support for a deletion. Thus, only the CNVs within the PARK2 locus could be molecularly validated and associated with PD susceptibility.
Allergic rhinitis (AR) affects up to 80% of children with asthma and increases asthma severity. Thymic stromal lymphopoietin (TSLP) is a key mediator of allergic inflammation. The role of the TSLP gene (TSLP) in the pathogenesis of AR has not been studied.
To test for associations between variants in TSLP, TSLP-related genes, and AR in children with asthma.
We genotyped 15 single nucleotide polymorphisms (SNPs) in TSLP, OX40L, IL7R, and RXRα in three independent cohorts: 592 asthmatic Costa Rican children and their parents, 422 nuclear families of North American children with asthma, and 239 Swedish children with asthma. We tested for associations between these SNPs and AR. As we previously reported sex-specific effects for TSLP, we performed overall and sex-stratified analyses. We additionally performed secondary analyses for gene-by-gene interactions.
Across the three cohorts, the T allele of TSLP SNP rs1837253 was undertransmitted in boys with AR and asthma as compared to boys with asthma alone. The SNP was associated with reduced odds for AR (odds ratios ranging from 0.56 to 0.63, with corresponding Fisher's combined P value of 1.2 × 10-4). Our findings were significant after accounting for multiple comparisons. SNPs in OX40L, IL7R, and RXRα were not consistently associated with AR in children with asthma. There were nominally significant interactions between gene pairs.
TSLP SNP rs1837253 is associated with reduced odds for AR in boys with asthma. Our findings support a role for TSLP in the pathogenesis of AR in children with asthma.
Single nucleotide polymorphisms (SNPs) in thymic stromal lymphopoietin (TSLP) have been associated with IgE (in girls) and asthma (in general). We sought to determine whether TSLP SNPs are associated with asthma in a sex-specific fashion.
We conducted regular and sex-stratified analyses of association between SNPs in TSLP and asthma in families of asthmatic children in Costa Rica. Significant findings were replicated in white and African-American participants in the Childhood Asthma Management Program, in African Americans in the Genomic Research on Asthma in the African Diaspora study, in whites and Hispanics in the Children’s Health Study, and in whites in the Framingham Heart Study (FHS).
Two SNPs in TSLP (rs1837253 and rs2289276) were significantly associated with a reduced risk of asthma in combined analyses of all cohorts (p values of 2×10−5 and 1×10−5, respectively). In a sex-stratified analysis, the T allele of rs1837253 was significantly associated with a reduced risk of asthma in males only (p= 3×10−6). Alternately, the T allele of rs2289276 was significantly associated with a reduced risk of asthma in females only (p= 2×10−4). Findings for rs2289276 were consistent in all cohorts except the FHS.
TSLP variants are associated with asthma in a sex-specific fashion.
asthma; genetic association; sex-specific; thymic stromal lymphopoietin; TSLP
Asthma is a chronic respiratory disease whose genetic basis has been explored for over two decades, most recently via genome-wide association studies. We sought to find asthma-susceptibility variants by using probands from a single population in both family-based and case-control association designs.
We used probands from the Childhood Asthma Management Program (CAMP) in two primary genome-wide association study designs: (1) probands were combined with publicly available population controls in a case-control design, and (2) probands and their parents were used in a family-based design. We followed a two-stage replication process utilizing three independent populations to validate our primary findings.
We found that single nucleotide polymorphisms with similar case-control and family-based association results were more likely to replicate in the independent populations, than those with the smallest p-values in either the case-control or family-based design alone. The single nucleotide polymorphism that showed the strongest evidence for association to asthma was rs17572584, which replicated in 2/3 independent populations with an overall p-value among replication populations of 3.5E-05. This variant is near a gene that encodes an enzyme that has been implicated to act coordinately with modulators of Th2 cell differentiation and is expressed in human lung.
Our results suggest that using probands from family-based studies in case-control designs, and combining results of both family-based and case-control approaches, may be a way to augment our ability to find SNPs associated with asthma and other complex diseases.
Measurements of lung function by spirometry are heritable traits that reflect respiratory health and predict morbidity and mortality. We meta-analyzed genome-wide association studies for two clinically important measures, forced expiratory volume in the first second (FEV1) and its ratio to forced vital capacity (FEV1/FVC), an indicator of airflow obstruction. This meta-analysis included 20,890 participants of European ancestry from four CHARGE consortium studies: Atherosclerosis Risk in Communities (ARIC), Cardiovascular Health Study (CHS), Framingham Heart Study (FHS), and Rotterdam Study (RS). We identified eight loci associated with FEV1/FVC (HHIP, GPR126, ADAM19, AGER-PPT2, FAM13A, PTCH1, PID1, and HTR4) and one locus associated with FEV1 (INTS12-GSTCD-NPNT) at or near genome-wide significance (P<5×10−8) in CHARGE; all but 3 loci (FAM13A, PTCH1, and PID1) replicated with the SpiroMeta consortium. Our findings of novel loci influencing pulmonary function may offer insights into chronic lung disease pathogenesis.