As part of an international consortium, case-parent trios were collected for a genome wide association study of isolated, non-syndromic oral clefts, including cleft lip (CL), cleft palate (CP) and cleft lip and palate (CLP). Non-syndromic oral clefts have a complex and heterogeneous etiology. Risk is influenced by genes, environmental factors, and differs markedly by gender. Family based association tests (FBAT) were used on 14,486 SNPs spanning the X chromosome, stratified by type of cleft and racial group. Significant results even after multiple comparisons correction were obtained for the Duchene’s muscular dystrophy (DMD) gene, the largest single gene in the human genome, among CL/P trios (both CL and CLP combined). When stratified into groups of European and Asian ancestry, stronger signals were obtained for Asians. Although conventional sliding window haplotype analysis showed no increase in significance, analysis selected combinations of the 25 most significant SNPs in DMD identified four SNPs together that attained genome-wide significance among Asian CL/P trios, raising the possibility of interaction between distant SNPs within DMD.
oral clefts; case-parent trios; X-linked; family-based association; DMD
Motivation: Family-based designs are regaining popularity for genomic sequencing studies because they provide a way to test cosegregation with disease of variants that are too rare in the population to be tested individually in a conventional case–control study.
Results: Where only a few affected subjects per family are sequenced, the probability that any variant would be shared by all affected relatives—given it occurred in any one family member—provides evidence against the null hypothesis of a complete absence of linkage and association. A P-value can be obtained as the sum of the probabilities of sharing events as (or more) extreme in one or more families. We generalize an existing closed-form expression for exact sharing probabilities to more than two relatives per family. When pedigree founders are related, we show that an approximation of sharing probabilities based on empirical estimates of kinship among founders obtained from genome-wide marker data is accurate for low levels of kinship. We also propose a more generally applicable approach based on Monte Carlo simulations. We applied this method to a study of 55 multiplex families with apparent non-syndromic forms of oral clefts from four distinct populations, with whole exome sequences available for two or three affected members per family. The rare single nucleotide variant rs149253049 in ADAMTS9 shared by affected relatives in three Indian families achieved significance after correcting for multiple comparisons (p=2×10−6).
Availability and implementation: Source code and binaries of the R package RVsharing are freely available for download at http://cran.r-project.org/web/packages/RVsharing/index.html.
firstname.lastname@example.org or email@example.com
Supplementary data are available at Bioinformatics online.
Genome-wide association studies have identified numerous single nucleotide polymorphisms (SNPs) as associated with colorectal cancer (CRC) risk in populations of European descent. However, their utility for predicting risk to CRC in Asians remains unknown. A case-cohort study (random sub-cohort N = 1,685) from the Korean Cancer Prevention Study-II (KCPS-II) (N = 145,842) was used. Twenty-three SNPs identified in previous 47 studies were genotyped on the KCPS-II sub-cohort members. A genetic risk score (GRS) was calculated by summing the number of risk alleles over all SNPs. Prediction models with or without GRS were evaluated in terms of the area under the receiver operating characteristic curve (AUROC) and the continuous net reclassification index (NRI).
Seven of 23 SNPs showed significant association with CRC and rectal cancer in Koreans, but not with colon cancer alone. AUROCs (95% CI) for traditional risk score (TRS) alone and TRS plus GRS were 0.73 (0.69–0.78) and 0.74 (0.70–0.78) for CRC, and 0.71 (0.65–0.77) and 0.74 (0.68–0.79) for rectal cancer, respectively. The NRI (95% CI) for a prediction model with GRS compared to the model with TRS alone was 0.17 (-0.05-0.37) for CRC and 0.41 (0.10–0.68) for rectal cancer alone.
Our results indicate genetic variants may be useful for predicting risk to CRC in the Koreans, especially risk for rectal cancer alone. Moreover, this study suggests effective prediction models for colon and rectal cancer should be developed separately.
Electronic supplementary material
The online version of this article (doi:10.1186/s12863-015-0207-y) contains supplementary material, which is available to authorized users.
Single nucleotide polymorphisms; Gene-traditional risk score; Colorectal cancer
Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence around 1 in 700 live births. The Runt-related transcription factor 2 (RUNX2) gene has been suggested as a candidate gene for CL/P based largely on mouse models; however, no human studies have focused on RUNX2 as a risk factor for CL/P. This study examines the association between markers in RUNX2 and isolated, nonsyndromic CL/P using a case-parent trio design, while considering parent-of-origin effects. Case-parent trios from four populations (77 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 24 single nucleotide polymorphisms (SNPs) in the RUNX2 gene. We performed the transmission disequilibrium test on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test and the parent-of-origin likelihood ratio test (PO-LRT). When all trios were combined, the transmission asymmetry test revealed a block of 11 SNPs showing excess maternal transmission significant at the P < 0.01 level, plus one SNP (rs1934328) showing excess paternal transmission (P = 0.002). For the 11 SNPs showing excess maternal transmission, odds ratios of being transmitted to the case from the mother ranged between 3.00 and 4.00. The parent-of-origin likelihood ratio tests for equality of maternal and paternal transmission were significant for three individual SNPs (rs910586, rs2819861, and rs1934328). Thus, RUNX2 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.
RUNX2; oral cleft; maternal transmission effects; parent-of-origin
Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of 1 in 700 live births. The paired box (PAX) genes have been suggested as candidate genes for CL/P based largely on mouse models; however, few human studies have focused on this gene family. This study tests for association between markers in four PAX genes and CL/P using a case-parent trio design considering parent-of-origin effects. Trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 34 single nucleotide polymorphisms (SNPs) in the PAX3, PAX6, PAX7, and PAX9 genes. We performed the transmission disequilibrium test (TDT) on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test (TAT) and the parent-of- origin likelihood ratio test (PO-LRT). TDT analysis showed one SNP (rs766325) in PAX7 yielding evidence of linkage and association when parent-of-origin was not considered, with an OR(transmission) = 1.62 (P = 0.003), and five SNPs in PAX6 (including two pairs in near perfect linkage disequilibrium). TAT analysis of all trios revealed two SNPs in PAX7 and four SNPs in PAX3 showing significant excess maternal transmission. For these six SNPs, the maternal OR(transmission) ranged between 1.74 and 2.40, and PO-LRT was also significant (P-values = 0.035–0.012). When this analysis was limited to trios with male cases, SNPs in PAX7 showed higher maternal OR(transmission) and greater significance. PAX genes may influence the risk of CL/P through maternal effects, possibly imprinting, which seems to be stronger among male cases.
PAX genes; oral cleft; maternal transmission effects; parent-of-origin
Isolated cleft palate is among the most common human birth defects. The TCOF1 gene has been suggested as a candidate gene for cleft palate based on animal models. This study tests for association between markers in TCOF1 and isolated, nonsyndromic cleft palate using a case-parent trio design considering parent-of-origin effects. Case-parent trios from three populations (comprising a total of 81 case-parent trios) were genotyped for single nucleotide polymorphisms (SNPs) in the TCOF1 gene. We used the transmission disequilibrium test and the transmission asymmetry test on individual SNPs. When all trios were combined, the odds ratio for transmission of the minor allele, OR(transmission), was significant for SNP rs15251 (OR = 2.88, P = 0.007), as well as rs2255796 and rs2569062 (OR = 2.08, P = 0.03; OR = 2.43, P = 0.041; respectively) when parent of origin was not considered. The transmission asymmetry test also revealed one SNP (rs15251) showing excess maternal transmission significant at the P = 0.005 level (OR = 6.50). Parent-of-origin effects were assessed using the parent-of-origin likelihood ratio test on both SNPs and haplotypes. While the parent-of-origin likelihood ratio test was only marginally significant for this SNP (P = 0.136), analysis of haplotypes of rs2255796 and rs15251 suggested excess maternal transmission. Therefore, these data suggest TCOF1 may influence risk of cleft palate through a parent-of-origin effect.
TCOF1; oral cleft palate; maternal transmission effects; parent-of-origin
The genetic risk factors for susceptibility to chronic obstructive
pulmonary disease (COPD) are still largely unknown. Additional genetic
variants are likely to be identified by genome-wide association studies in
larger cohorts or specific subgroups.
Genome-wide association analysis in COPDGene (non-Hispanic whites and
African-Americans) was combined with existing data from the ECLIPSE,
NETT/NAS, and GenKOLS (Norway) studies. Analyses were performed both using
all moderate-to-severe cases and the subset of severe cases. Top loci not
previously described as genome-wide significant were genotyped in the ICGN
study, and results combined in a joint meta-analysis.
Analysis of a total of 6,633 moderate-to-severe cases and 5,704
controls confirmed association at three known loci:
CHRNA3/CHRNA5/IREB2, FAM13A, and HHIP
(10−12 < P < 10−14),
and also showed significant evidence of association at a novel locus near
RIN3 (overall P, including ICGN =
5•4×10−9). In the severe COPD analysis
(n=3,497), the effects at two of three previously described loci were
significantly stronger; we also identified two additional loci previously
reported to affect gene expression of MMP12 and
TGFB2 (overall P = 2•6x10−9
and 8•3×10−9). RIN3 and
TGFB2 expression levels were reduced in a set of Lung
Tissue Research Consortium COPD lung tissue samples compared with
In a genome-wide study of COPD, we confirmed associations at three
known loci and found additional genome-wide significant associations with
moderate-to-severe COPD near RIN3 and with severe COPD near
MMP12 and TGFB2. Genetic variants,
apart from alpha-1 antitrypsin deficiency, increase the risk of COPD. Our
analysis of severe COPD suggests additional genetic variants may be
identified by focusing on this subgroup.
National Heart, Lung, and Blood Institute; the COPD Foundation
through contributions from AstraZeneca, Boehringer Ingelheim, Novartis, and
Sepracor; GlaxoSmithKline; Centers for Medicare and Medicaid Services;
Agency for Healthcare Research and Quality; US Department of Veterans
COPD patients have a great burden of comorbidity. However, it is not well established whether this is due to shared risk factors such as smoking, if they impact patients exercise capacity and quality of life, or whether there are racial disparities in their impact on COPD.
We analyzed data from 10,192 current and ex-smokers with (cases) and without COPD (controls) from the COPDGene® cohort to establish risk for COPD comorbidities adjusted for pertinent covariates. In adjusted models, we examined comorbidities prevalence and impact in African-Americans (AA) and Non-Hispanic Whites (NHW).
Comorbidities are more common in COPD compared to those with normal spirometry (controls), and the risk persists after adjustments for covariates including pack-years smoked. After adjustment for confounders, eight conditions were independently associated with worse exercise capacity, quality of life and dyspnea. There were racial disparities in the impact of comorbidities on exercise capacity, dyspnea and quality of life, presence of osteoarthritis and gastroesophageal reflux disease having a greater negative impact on all three outcomes in AAs than NHWs (p<0.05 for all interaction terms).
Individuals with COPD have a higher risk for comorbidities than controls, an important finding shown for the first time comprehensively after accounting for confounders. Individual comorbidities are associated with worse exercise capacity, quality of life, and dyspnea, in African-Americans compared to non-Hispanic Whites.
COPD; Comorbidities; Race
Rationale: Pulmonary emphysema overlaps partially with spirometrically defined chronic obstructive pulmonary disease and is heritable, with moderately high familial clustering.
Objectives: To complete a genome-wide association study (GWAS) for the percentage of emphysema-like lung on computed tomography in the Multi-Ethnic Study of Atherosclerosis (MESA) Lung/SNP Health Association Resource (SHARe) Study, a large, population-based cohort in the United States.
Methods: We determined percent emphysema and upper-lower lobe ratio in emphysema defined by lung regions less than −950 HU on cardiac scans. Genetic analyses were reported combined across four race/ethnic groups: non-Hispanic white (n = 2,587), African American (n = 2,510), Hispanic (n = 2,113), and Chinese (n = 704) and stratified by race and ethnicity.
Measurements and Main Results: Among 7,914 participants, we identified regions at genome-wide significance for percent emphysema in or near SNRPF (rs7957346; P = 2.2 × 10−8) and PPT2 (rs10947233; P = 3.2 × 10−8), both of which replicated in an additional 6,023 individuals of European ancestry. Both single-nucleotide polymorphisms were previously implicated as genes influencing lung function, and analyses including lung function revealed independent associations for percent emphysema. Among Hispanics, we identified a genetic locus for upper-lower lobe ratio near the α-mannosidase–related gene MAN2B1 (rs10411619; P = 1.1 × 10−9; minor allele frequency [MAF], 4.4%). Among Chinese, we identified single-nucleotide polymorphisms associated with upper-lower lobe ratio near DHX15 (rs7698250; P = 1.8 × 10−10; MAF, 2.7%) and MGAT5B (rs7221059; P = 2.7 × 10−8; MAF, 2.6%), which acts on α-linked mannose. Among African Americans, a locus near a third α-mannosidase–related gene, MAN1C1 (rs12130495; P = 9.9 × 10−6; MAF, 13.3%) was associated with percent emphysema.
Conclusions: Our results suggest that some genes previously identified as influencing lung function are independently associated with emphysema rather than lung function, and that genes related to α-mannosidase may influence risk of emphysema.
emphysema; computed tomography; multiethnic; cohort study; genetic association
ATP-binding cassette (ABC) proteins in the placenta regulate fetal exposure to xenobiotics. We hypothesized that functional polymorphisms in ABC genes influence risk for non-syndromic oral clefts (NSOC). Both family-based and case–control studies were undertaken to evaluate the association of nine potentially functional single-nucleotide polymorphisms within four ABC genes with risk of NSOC. Peripheral blood DNA from a total of 150 NSOC case-parent trios from Singapore and Taiwan were genotyped, as was cord blood DNA from 189 normal Chinese neonates used as controls. In trios, significant association was observed between the ABCB1 single-nucleotide polymorphisms and NSOC (P<0.05). Only ABCB1 rs1128503 retained significant association after Bonferroni correction (odds ratio (OR)=2.04; 95% confidence interval (CI)=1.42–2.98), while rs2032582 and rs1045642 showed nominal significance. Association with rs1128503 was replicated in a case–control analysis comparing NSOC probands with controls (OR=1.58; 95% CI=1.12–2.23). A comparison between the mothers of probands and controls showed no evidence of association, suggesting NSOC risk is determined by fetal and not maternal ABCB1 genotype. The two studies produced a combined OR of 1.79 (95% CI=1.38–2.30). The T-allele at rs1128503 was associated with higher risk. This study thus provides evidence that potentially functional polymorphisms in fetal ABCB1 modulate risk for NSOC, presumably through suboptimal exclusion of xenobiotics at the fetal–maternal interface.
cleft palate; cleft lip; ATP-binding cassette transporters; single-nucleotide polymorphism; disease susceptibility; placenta
Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common disorder with complex etiology. The Bone Morphogenetic Protein 4 gene (BMP4) has been considered a prime candidate gene with evidence accumulated from animal experimental studies, human linkage studies, as well as candidate gene association studies. The aim of the current study is to test for linkage and association between BMP4 and NSCL/P that could be missed in genome-wide association studies (GWAS) when genotypic (G) main effects alone were considered.
We performed the analysis considering G and interactions with multiple maternal environmental exposures using additive conditional logistic regression models in 895 Asian and 681 European complete NSCL/P trios. Single nucleotide polymorphisms (SNPs) that passed the quality control criteria among 122 genotyped and 25 imputed single nucleotide variants in and around the gene were used in analysis. Selected maternal environmental exposures during 3 months prior to and through the first trimester of pregnancy included any personal tobacco smoking, any environmental tobacco smoke in home, work place or any nearby places, any alcohol consumption and any use of multivitamin supplements. A novel significant association held for rs7156227 among Asian NSCL/P and non-syndromic cleft lip and palate (NSCLP) trios after Bonferroni correction which was not seen when G main effects alone were considered in either allelic or genotypic transmission disequilibrium tests. Odds ratios for carrying one copy of the minor allele without maternal exposure to any of the four environmental exposures were 0.58 (95%CI = 0.44, 0.75) and 0.54 (95%CI = 0.40, 0.73) for Asian NSCL/P and NSCLP trios, respectively. The Bonferroni P values corrected for the total number of 117 tested SNPs were 0.0051 (asymptotic P = 4.39*10−5) and 0.0065 (asymptotic P = 5.54*10−5), accordingly. In European trios, no significant association was seen for any SNPs after Bonferroni corrections for the total number of 120 tested SNPs.
Our findings add evidence from GWAS to support the role of BMP4 in susceptibility to NSCL/P originally identified in linkage and candidate gene association studies.
Preserved Ratio Impaired Spirometry (PRISm), defined as a reduced FEV1 in the setting of a preserved FEV1/FVC ratio, is highly prevalent and is associated with increased respiratory symptoms, systemic inflammation, and mortality. Studies investigating quantitative chest tomographic features, genetic associations, and subtypes in PRISm subjects have not been reported.
Data from current and former smokers enrolled in COPDGene (n = 10,192), an observational, cross-sectional study which recruited subjects aged 45–80 with ≥10 pack years of smoking, were analyzed. To identify epidemiological and radiographic predictors of PRISm, we performed univariate and multivariate analyses comparing PRISm subjects both to control subjects with normal spirometry and to subjects with COPD. To investigate common genetic predictors of PRISm, we performed a genome-wide association study (GWAS). To explore potential subgroups within PRISm, we performed unsupervised k-means clustering.
The prevalence of PRISm in COPDGene is 12.3%. Increased dyspnea, reduced 6-minute walk distance, increased percent emphysema and decreased total lung capacity, as well as increased segmental bronchial wall area percentage were significant predictors (p-value <0.05) of PRISm status when compared to control subjects in multivariate models. Although no common genetic variants were identified on GWAS testing, a significant association with Klinefelter’s syndrome (47XXY) was observed (p-value < 0.001). Subgroups identified through k-means clustering include a putative “COPD-subtype”, “Restrictive-subtype”, and a highly symptomatic “Metabolic-subtype”.
PRISm subjects are clinically and genetically heterogeneous. Future investigations into the pathophysiological mechanisms behind and potential treatment options for subgroups within PRISm are warranted.
Clinicaltrials.gov Identifier: NCT000608764.
Electronic supplementary material
The online version of this article (doi:10.1186/s12931-014-0089-y) contains supplementary material, which is available to authorized users.
Spirometry; Restriction; Lung diseases; Smoking
Allergic rhinitis is a common disease whose genetic basis is incompletely explained. We report an integrated genomic analysis of allergic rhinitis.
We performed genome wide association studies (GWAS) of allergic rhinitis in 5633 ethnically diverse North American subjects. Next, we profiled gene expression in disease-relevant tissue (peripheral blood CD4+ lymphocytes) collected from subjects who had been genotyped. We then integrated the GWAS and gene expression data using expression single nucleotide (eSNP), coexpression network, and pathway approaches to identify the biologic relevance of our GWAS.
GWAS revealed ethnicity-specific findings, with 4 genome-wide significant loci among Latinos and 1 genome-wide significant locus in the GWAS meta-analysis across ethnic groups. To identify biologic context for these results, we constructed a coexpression network to define modules of genes with similar patterns of CD4+ gene expression (coexpression modules) that could serve as constructs of broader gene expression. 6 of the 22 GWAS loci with P-value ≤ 1x10−6 tagged one particular coexpression module (4.0-fold enrichment, P-value 0.0029), and this module also had the greatest enrichment (3.4-fold enrichment, P-value 2.6 × 10−24) for allergic rhinitis-associated eSNPs (genetic variants associated with both gene expression and allergic rhinitis). The integrated GWAS, coexpression network, and eSNP results therefore supported this coexpression module as an allergic rhinitis module. Pathway analysis revealed that the module was enriched for mitochondrial pathways (8.6-fold enrichment, P-value 4.5 × 10−72).
Our results highlight mitochondrial pathways as a target for further investigation of allergic rhinitis mechanism and treatment. Our integrated approach can be applied to provide biologic context for GWAS of other diseases.
Genome-wide association study; Allergic rhinitis; Coexpression network; Expression single-nucleotide polymorphism; Coexpression module; Pathway; Mitochondria; Hay fever; Allergy
The 2011 Global Strategy for the Diagnosis, Management, and Prevention of COPD (GOLD) consensus report uses symptoms, exacerbation history and FEV1% to define four categories: A, low symptoms/low risk; B, high symptoms/low risk; C, low symptoms/high risk; and D, high symptoms/high risk where risk refers to exacerbations, hospitalization and death. Our objective was to determine (1) the influence of symptom instrument on category membership and (2) prospective exacerbation risk by category.
4,484 COPD subjects from COPDGene were analyzed. All subjects had smoking history ≥ 10 pack-years and FEV1/FVC<0·7. Categories were defined using the modified Medical Research Council Dyspnea [mMRC] (0–1 versus ≥ 2) and the Saint George’s Respiratory Questionnaire [SGRQ] (≥25 versus <25 as a surrogate for the COPD Assessment Test ≥ 10 versus <10) in addition to COPD exacerbations in the prior year (<2 versus ≥ 2), and FEV1% predicted (≥50 versus <50).
Category assignment using mMRC versus SGRQ were similar but not identical. Using the mMRC, category assignments were 34% A, 21% B, 8% C and 38% D and for SGRQ were 29% A, 25% B, 5% C and 41% D (kappa=0·77). Significant heterogeneity in exacerbation rates (exacerbations/person-year) were seen particularly within the D group, depending on the risk factor that determined category assignment (lung function (0·89), prior exacerbation history (1·34) or both (1·86), p<0·001.
The GOLD classification emphasizes the importance of symptoms and exacerbation risk in assessing COPD severity. The choice of symptom measure influences category assignment. The relative number of subjects with low symptoms and high risk for exacerbations (category C) is low. Differences in exacerbation rates for subjects in the highest risk category D were seen depending on whether risk was based on lung function, exacerbation history or both.
Characterization of genetic admixture of populations in the Americas and the Caribbean is of interest for anthropological, epidemiological, and historical reasons. Asthma has a higher prevalence and is more severe in populations with a high African component. Association of African ancestry with asthma has been demonstrated. We estimated admixture proportions of samples from six trihybrid populations of African descent and determined the relationship between African ancestry and asthma and total serum IgE levels (tIgE). We genotyped 237 ancestry informative markers in asthmatics and nonasthmatic controls from Barbados (190/277), Jamaica (177/529), Brazil (40/220), Colombia (508/625), African Americans from New York (207/171), and African Americans from Baltimore/Washington, D.C. (625/757). We estimated individual ancestries and evaluated genetic stratification using Structure and principal component analysis. Association of African ancestry and asthma and tIgE was evaluated by regression analysis. Mean SD African ancestry ranged from 0.76 ± 0.10 among Barbadians to 0.33 ± 0.13 in Colombians. The European component varied from 0.14 ± 0.05 among Jamaicans and Barbadians to 0.26 ± 0.08 among Colombians. African ancestry was associated with risk for asthma in Colombians (odds ratio (OR) = 4.5, P = 0.001) Brazilians (OR = 136.5, P = 0.003), and African Americans of New York (OR: 4.7; P = 0.040). African ancestry was also associated with higher tIgE levels among Colombians (β = 1.3, P = 0.04), Barbadians (β = 3.8, P = 0.03), and Brazilians (β = 1.6, P = 0.03). Our findings indicate that African ancestry can account for, at least in part, the association between asthma and its associated trait, tIgE levels.
African; asthma; ancestry
Hedgehog Interacting Protein (HHIP) was implicated in chronic obstructive pulmonary disease (COPD) by genome-wide association studies (GWAS). However, it remains unclear how HHIP contributes to COPD pathogenesis. To identify genes regulated by HHIP, we performed gene expression microarray analysis in a human bronchial epithelial cell line (Beas-2B) stably infected with HHIP shRNAs. HHIP silencing led to differential expression of 296 genes; enrichment for variants nominally associated with COPD was found. Eighteen of the differentially expressed genes were validated by real-time PCR in Beas-2B cells. Seven of 11 validated genes tested in human COPD and control lung tissues demonstrated significant gene expression differences. Functional annotation indicated enrichment for extracellular matrix and cell growth genes. Network modeling demonstrated that the extracellular matrix and cell proliferation genes influenced by HHIP tended to be interconnected. Thus, we identified potential HHIP targets in human bronchial epithelial cells that may contribute to COPD pathogenesis.
Hedgehog interacting protein (HHIP); Gene expression profiling; COPD (Chronic obstructive pulmonary disease); extracellular matrix (ECM); network modeling
Chronic mucus hypersecretion (CMH) is associated with an increased frequency of respiratory infections, excess lung function decline, and increased hospitalisation and mortality rates in the general population. It is associated with smoking, but it is unknown why only a minority of smokers develops CMH. A plausible explanation for this phenomenon is a predisposing genetic constitution. Therefore, we performed a genome wide association (GWA) study of CMH in Caucasian populations.
GWA analysis was performed in the NELSON-study using the Illumina 610 array, followed by replication and meta-analysis in 11 additional cohorts. In total 2,704 subjects with, and 7,624 subjects without CMH were included, all current or former heavy smokers (≥20 pack-years). Additional studies were performed to test the functional relevance of the most significant single nucleotide polymorphism (SNP).
A strong association with CMH, consistent across all cohorts, was observed with rs6577641 (p = 4.25×10−6, OR = 1.17), located in intron 9 of the special AT-rich sequence-binding protein 1 locus (SATB1) on chromosome 3. The risk allele (G) was associated with higher mRNA expression of SATB1 (4.3×10−9) in lung tissue. Presence of CMH was associated with increased SATB1 mRNA expression in bronchial biopsies from COPD patients. SATB1 expression was induced during differentiation of primary human bronchial epithelial cells in culture.
Our findings, that SNP rs6577641 is associated with CMH in multiple cohorts and is a cis-eQTL for SATB1, together with our additional observation that SATB1 expression increases during epithelial differentiation provide suggestive evidence that SATB1 is a gene that affects CMH.
Cigarette smoking is the major environmental risk factor for chronic obstructive pulmonary disease (COPD). Genome-wide association studies have provided compelling associations for three loci with COPD. In this study, we aimed to estimate direct, i.e., independent from smoking, and indirect effects of those loci on COPD development using mediation analysis. We included a total of 3,424 COPD cases and 1,872 unaffected controls with data on two smoking-related phenotypes: lifetime average smoking intensity and cumulative exposure to tobacco smoke (pack years). Our analysis revealed that effects of two linked variants (rs1051730 and rs8034191) in the AGPHD1/CHRNA3 cluster on COPD development are significantly, yet not entirely, mediated by the smoking-related phenotypes. Approximately 30 % of the total effect of variants in the AGPHD1/CHRNA3 cluster on COPD development was mediated by pack years. Simultaneous analysis of modestly (r2 = 0.21) linked markers in CHRNA3 and IREB2 revealed that an even larger (~42 %) proportion of the total effect of the CHRNA3 locus on COPD was mediated by pack years after adjustment for an IREB2 single nucleotide polymorphism. This study confirms the existence of direct effects of the AGPHD1/CHRNA3, IREB2, FAM13A and HHIP loci on COPD development. While the association of the AGPHD1/CHRNA3 locus with COPD is significantly mediated by smoking-related phenotypes, IREB2 appears to affect COPD independently of smoking.
Neuronal nicotinic acetylcholine receptor (nAChR) genes (CHRNA5/CHRNA3/CHRNB4) have been reproducibly associated with nicotine dependence, smoking behaviors, and lung cancer risk. Of the few reports that have focused on early smoking behaviors, association results have been mixed. This meta-analysis examines early smoking phenotypes and SNPs in the gene cluster to determine: (1) whether the most robust association signal in this region (rs16969968) for other smoking behaviors is also associated with early behaviors, and/or (2) if additional statistically independent signals are important in early smoking. We focused on two phenotypes: age of tobacco initiation (AOI) and age of first regular tobacco use (AOS). This study included 56,034 subjects (41 groups) spanning nine countries and evaluated five SNPs including rs1948, rs16969968, rs578776, rs588765, and rs684513. Each dataset was analyzed using a centrally generated script. Meta-analyses were conducted from summary statistics. AOS yielded significant associations with SNPs rs578776 (beta = 0.02, P = 0.004), rs1948 (beta = 0.023, P = 0.018), and rs684513 (beta = 0.032, P = 0.017), indicating protective effects. There were no significant associations for the AOI phenotype. Importantly, rs16969968, the most replicated signal in this region for nicotine dependence, cigarettes per day, and cotinine levels, was not associated with AOI (P = 0.59) or AOS (P = 0.92). These results provide important insight into the complexity of smoking behavior phenotypes, and suggest that association signals in the CHRNA5/A3/B4 gene cluster affecting early smoking behaviors may be different from those affecting the mature nicotine dependence phenotype.
CHRNA5; CHRNA3; CHRNB4; meta-analysis; nicotine; smoke
Copy number variants (CNVs) may play an important part in the development of common birth defects such as oral clefts, and individual patients with multiple birth defects (including clefts) have been shown to carry small and large chromosomal deletions. In this paper we investigate de novo deletions defined as DNA segments missing in an oral cleft proband but present in both unaffected parents. We compare de novo deletion frequencies in children of European ancestry with an isolated, non-syndromic oral cleft to frequencies in children of European ancestry from randomly sampled trios.
We identified a genome-wide significant 62 kilo base (kb) non-coding region on chromosome 7p14.1 where de novo deletions occur more frequently among oral cleft cases than controls. We also observed wider de novo deletions among cleft lip and palate (CLP) cases than seen among cleft palate (CP) and cleft lip (CL) cases.
This study presents a region where de novo deletions appear to be involved in the etiology of oral clefts, although the underlying biological mechanisms are still unknown. Larger de novo deletions are more likely to interfere with normal craniofacial development and may result in more severe clefts. Study protocol and sample DNA source can severely affect estimates of de novo deletion frequencies. Follow-up studies are needed to further validate these findings and to potentially identify additional structural variants underlying oral clefts.
Oral clefts; DNA copy numbers; de novo deletions; Case-parent trios
Nonsyndromic cleft palate (CP) is one of the most common human birth defects and both genetic and environmental risk factors contribute to its etiology. We conducted a genome-wide association study (GWAS) using 550 CP case-parent trios ascertained in an international consortium. Stratified analysis among trios with different ancestries was performed to test for GxE interactions with common maternal exposures using conditional logistic regression models. While no single nucleotide polymorphism (SNP) achieved genome-wide significance when considered alone, markers in SLC2A9 and the neighboring WDR1 on chromosome 4p16.1 gave suggestive evidence of gene-environment interaction with environmental tobacco smoke (ETS) among 259 Asian trios when the models included a term for GxE interaction. Multiple SNPs in these two genes were associated with increased risk of nonsyndromic CP if the mother was exposed to ETS during the peri-conceptual period (3 months prior to conception through the first trimester). When maternal ETS was considered, fifteen of 135 SNPs mapping to SLC2A9 and 9 of 59 SNPs in WDR1 gave P values approaching genome-wide significance (10−6
Isolated, non-syndromic cleft lip with or without cleft palate (iCL±P) is a common human congenital malformation with a complex and heterogeneous etiology. Genes coding for fibroblast growth factors and their receptors (FGF/FGFR genes) are excellent candidate genes.
We tested single nucleotide polymorphic (SNP) markers in 10 FGF/FGFR genes (including FGFBP1, FGF2, FGF10, FGF18, FGFR1, FGFR2, FGF19, FGF4, FGF3, and FGF9) for genotypic effects, interactions with one another, and with common maternal environmental exposures in 221 Asian and 76 Maryland case-parent trios ascertained through a child with iCL±P.
Both FGFR1 and FGF19 yielded evidence of linkage and association in the transmission disequilibrium test, confirming previous evidence. Haplotypes of three SNPs in FGFR1 were nominally significant among Asian trios. Estimated ORs for individual SNPs and haplotypes of multiple markers in FGF19 ranged between1.31-1.87. We also found suggestive evidence of maternal genotypic effects for markers in FGF2 and FGF10 among Asian trios. Tests for gene-environment (GxE) interaction between markers in FGFR2 and maternal smoking or multivitamin supplementation yielded significant evidence of GxE interaction separately. Tests of gene-gene (GxG) interaction using Cordell's method yielded significant evidence between SNPs in FGF9 and FGF18, which was confirmed in an independent sample of trios from an international consortium.
Our results suggest several genes in the FGF/FGFR family may influence risk to iCL±P through distinct biological mechanisms.
FGF/FGFR; oral clefts; maternal effects; gene-environment interaction; gene-gene interaction
Accelerated lung function decline is a key COPD phenotype; however its genetic control remains largely unknown.
We performed a genome-wide association study using the Illumina Human660W-Quad v.1_A BeadChip. Generalized estimation equations were used to assess genetic contributions to lung function decline over a 5-year period in 4,048 European-American Lung Health Study participants with largely mild COPD. Genotype imputation was performed using reference HapMap II data. To validate regions meeting genome-wide significance, replication of top SNPs was attempted in independent cohorts. Three genes (TMEM26, ANK3 and FOXA1) within the regions of interest were selected for tissue expression studies using immunohistochemistry.
Measurements and Main Results
Two intergenic SNPs (rs10761570, rs7911302) on chromosome 10 and one SNP on chromosome 14 (rs177852) met genome-wide significance after Bonferroni. Further support for the chromosome 10 region was obtained by imputation, the most significantly associated imputed SNPs (rs10761571, rs7896712) being flanked by observed markers rs10761570 and rs7911302. Results were not replicated in four general population cohorts or a smaller cohort of subjects with moderate to severe COPD; however, we show novel expression of genes near regions of significantly associated SNPS, including TMEM26 and FOXA1 in airway epithelium and lung parenchyma, and ANK3 in alveolar macrophages. Levels of expression were associated with lung function and COPD status.
We identified two novel regions associated with lung function decline in mild COPD. Genes within these regions were expressed in relevant lung cells and their expression related to airflow limitation suggesting they may represent novel candidate genes for COPD susceptibility.
COPD; lung function decline; GWAS; genome wide association; genes; polymorphisms
Rationale: A genome-wide association study (GWAS) for circulating chronic obstructive pulmonary disease (COPD) biomarkers could identify genetic determinants of biomarker levels and COPD susceptibility.
Objectives: To identify genetic variants of circulating protein biomarkers and novel genetic determinants of COPD.
Methods: GWAS was performed for two pneumoproteins, Clara cell secretory protein (CC16) and surfactant protein D (SP-D), and five systemic inflammatory markers (C-reactive protein, fibrinogen, IL-6, IL-8, and tumor necrosis factor-α) in 1,951 subjects with COPD. For genome-wide significant single nucleotide polymorphisms (SNPs) (P < 1 × 10−8), association with COPD susceptibility was tested in 2,939 cases with COPD and 1,380 smoking control subjects. The association of candidate SNPs with mRNA expression in induced sputum was also elucidated.
Measurements and Main Results: Genome-wide significant susceptibility loci affecting biomarker levels were found only for the two pneumoproteins. Two discrete loci affecting CC16, one region near the CC16 coding gene (SCGB1A1) on chromosome 11 and another locus approximately 25 Mb away from SCGB1A1, were identified, whereas multiple SNPs on chromosomes 6 and 16, in addition to SNPs near SFTPD, had genome-wide significant associations with SP-D levels. Several SNPs affecting circulating CC16 levels were significantly associated with sputum mRNA expression of SCGB1A1 (P = 0.009–0.03). Several SNPs highly associated with CC16 or SP-D levels were nominally associated with COPD in a collaborative GWAS (P = 0.001–0.049), although these COPD associations were not replicated in two additional cohorts.
Conclusions: Distant genetic loci and biomarker-coding genes affect circulating levels of COPD-related pneumoproteins. A subset of these protein quantitative trait loci may influence their gene expression in the lung and/or COPD susceptibility.
Clinical trial registered with www.clinicaltrials.gov (NCT 00292552).
biomarker; chronic obstructive pulmonary disease; genome-wide association study
Results 1-25 (74)
This will clear all selections from your clipboard. Do you wish proceed?
Clipboard is full! Please remove an item and try again.