Bone marrow-derived epithelial cells and hair follicle stem cells contribute to development of chronic cutaneous neoplasms.

We used allogeneic bone marrow transplantation (BMT) and a mouse multistage cutaneous carcinogenesis model to probe recruitment of bone marrow-derived epithelial cells (BMDECs) in skin tumors initiated with the carcinogen, dimethylbenz[a]anthracene (DMBA), and promoted with 12-O-tetradecanolyphorbol-13-acetate (TPA). BMDECs clustered in the lesional epithelium, expressed cytokeratins, proliferated, and stratified. We detected cytokeratin induction in plastic-adherent bone marrow cells (BMCs) cultured in the presence of filter-separated keratinocytes (KCs) and bone morphogenetic protein 5 (BMP5). Lineage-depleted BMCs migrated towards High Mobility Group Box 1 (HMGB1) protein and epidermal KCs in ex vivo invasion assays. Naive female mice receiving BMTs from DMBA-treated donors developed benign and malignant lesions after TPA promotion alone. We conclude that BMDECs contribute to the development of papillomas and dysplasia, demonstrating a systemic contribution to these lesions. Furthermore, carcinogen-exposed BMCs can initiate benign and malignant lesions upon tumor promotion. Ultimately, these findings may suggest targets for treatment of non-melanoma skin cancers.

Investigation of previously implicated genetic variants in chronic tic disorders: a transmission disequilibrium test approach.

Genetic studies in Tourette syndrome (TS) are characterized by scattered and poorly replicated findings. We aimed to replicate findings from candidate gene and genome-wide association studies (GWAS). Our cohort included 465 probands with chronic tic disorder (93% TS) and both parents from 412 families (some probands were siblings). We assessed 75 single nucleotide polymorphisms (SNPs) in 465 parent-child trios; 117 additional SNPs in 211 trios; and 4 additional SNPs in 254 trios. We performed SNP and gene-based transmission disequilibrium tests and compared nominally significant SNP results with those from a large independent case-control cohort. After quality control 71 SNPs were available in 371 trios; 112 SNPs in 179 trios; and 3 SNPs in 192 trios. 17 were candidate SNPs implicated in TS and 2 were implicated in obsessive-compulsive disorder (OCD) or autism spectrum disorder (ASD); 142 were tagging SNPs from eight monoamine neurotransmitter-related genes (including dopamine and serotonin); 10 were top SNPs from TS GWAS; and 13 top SNPs from attention-deficit/hyperactivity disorder, OCD, or ASD GWAS. None of the SNPs or genes reached significance after adjustment for multiple testing. We observed nominal significance for the candidate SNPs rs3744161 (TBCD) and rs4565946 (TPH2) and for five tagging SNPs; none of these showed significance in the independent cohort. Also, SLC1A1 in our gene-based analysis and two TS GWAS SNPs showed nominal significance, rs11603305 (intergenic) and rs621942 (PICALM). We found no convincing support for previously implicated genetic polymorphisms. Targeted re-sequencing should fully appreciate the relevance of candidate genes.

An Analytic Solution to the Computation of Power and Sample Size for Genetic Association Studies under a Pleiotropic Mode of Inheritance.

Our motivation here is to calculate the power of 3 statistical tests used when there are genetic traits that operate under a pleiotropic mode of inheritance and when qualitative phenotypes are defined by use of thresholds for the multiple quantitative phenotypes. Specifically, we formulate a multivariate function that provides the probability that an individual has a vector of specific quantitative trait values conditional on having a risk locus genotype, and we apply thresholds to define qualitative phenotypes (affected, unaffected) and compute penetrances and conditional genotype frequencies based on the multivariate function. We extend the analytic power and minimum-sample-size-necessary (MSSN) formulas for 2 categorical data-based tests (genotype, linear trend test [LTT]) of genetic association to the pleiotropic model. We further compare the MSSN of the genotype test and the LTT with that of a multivariate ANOVA (Pillai). We approximate the MSSN for statistics by linear models using a factorial design and ANOVA. With ANOVA decomposition, we determine which factors most significantly change the power/MSSN for all statistics. Finally, we determine which test statistics have the smallest MSSN. In this work, MSSN calculations are for 2 traits (bivariate distributions) only (for illustrative purposes). We note that the calculations may be extended to address any number of traits. Our key findings are that the genotype test usually has lower MSSN requirements than the LTT. More inclusive thresholds (top/bottom 25% vs. top/bottom 10%) have higher sample size requirements. The Pillai test has a much larger MSSN than both the genotype test and the LTT, as a result of sample selection. With these formulas, researchers can specify how many subjects they must collect to localize genes for pleiotropic phenotypes.

A PAX1 enhancer locus is associated with susceptibility to idiopathic scoliosis in females.

Idiopathic scoliosis (IS) is a common paediatric musculoskeletal disease that displays a strong female bias. By performing a genome-wide association study (GWAS) of 3,102 individuals, we identify significant associations with 20p11.22 SNPs for females (P=6.89 × 10(-9)) but not males (P=0.71). This association with IS is also found in independent female cohorts from the United States of America and Japan (overall P=2.15 × 10(-10), OR=1.30 (rs6137473)). Unexpectedly, the 20p11.22 IS risk alleles were previously associated with protection from early-onset alopecia, another sexually dimorphic condition. The 174-kb associated locus is distal to PAX1, which encodes paired box 1, a transcription factor involved in spine development. We identify a sequence in the associated locus with enhancer activity in zebrafish somitic muscle and spinal cord, an activity that is abolished by IS-associated SNPs. We thus identify a sexually dimorphic IS susceptibility locus, and propose the first functionally defined candidate mutations in an enhancer that may regulate expression in specific spinal cells.

Mapping genes with longitudinal phenotypes via Bayesian posterior probabilities.

Most association studies focus on disease risk, with less attention paid to disease progression or severity. These phenotypes require longitudinal data. This paper presents a new method for analyzing longitudinal data to map genes in both population-based and family-based studies. Using simulated systolic blood pressure measurements obtained from Genetic Analysis Workshop 18, we cluster the phenotype data into trajectory subgroups. We then use the Bayesian posterior probability of being in the high subgroup as a quantitative trait in an association analysis with genotype data. This method maintains high power (>80%) in locating genes known to affect the simulated phenotype for most specified significance levels (α). We believe that this method can be useful to aid in the discovery of genes that affect severity or progression of disease.

Chromatin profiling reveals regulatory network shifts and a protective role for hepatocyte nuclear factor 4α during colitis.

Transcriptional regulatory mechanisms likely contribute to the etiology of inflammatory bowel disease (IBD), as genetic variants associated with the disease are disproportionately found at regulatory elements. However, the transcription factors regulating colonic inflammation are unclear. To identify these transcription factors, we mapped epigenomic changes in the colonic epithelium upon inflammation. Epigenetic marks at transcriptional regulatory elements responded dynamically to inflammation and indicated a shift in epithelial transcriptional factor networks. Active enhancer chromatin structure at regulatory regions bound by the transcription factor hepatocyte nuclear factor 4α (HNF4A) was reduced during colitis. In agreement, upon an inflammatory stimulus, HNF4A was downregulated and showed a reduced ability to bind chromatin. Genetic variants that confer a predisposition to IBD map to HNF4A binding sites in the human colon cell line CaCo2, suggesting impaired HNF4A binding could underlie genetic susceptibility to IBD. Despite reduced HNF4A binding during inflammation, a temporal knockout model revealed HNF4A still actively protects against inflammatory phenotypes and promotes immune regulatory gene expression in the inflamed colonic epithelium. These findings highlight the potential for HNF4A agonists as IBD therapeutics.

A meta-analysis identifies adolescent idiopathic scoliosis association with LBX1 locus in multiple ethnic groups.

BACKGROUND: Adolescent idiopathic scoliosis (AIS) is a common rotational deformity of the spine that presents in children worldwide, yet its etiology is poorly understood. Recent genome-wide association studies (GWAS) have identified a few candidate risk loci. One locus near the chromosome 10q24.31 LBX1 gene (OMIM #604255) was originally identified by a GWAS of Japanese subjects and replicated in additional Asian populations. To extend this result, and to create larger AIS cohorts for the purpose of large-scale meta-analyses in multiple ethnicities, we formed a collaborative group called the International Consortium for Scoliosis Genetics (ICSG). METHODS: Here, we report the first ICSG study, a meta-analysis of the LBX1 locus in six Asian and three non-Asian cohorts. RESULTS: We find significant evidence for association of this locus with AIS susceptibility in all nine cohorts. Results for seven cohorts containing both genders yielded P=1.22×10-43 for rs11190870, and P=2.94×10-48 for females in all nine cohorts. Comparing the regional haplotype structures for three populations, we refined the boundaries of association to a ∼25 kb block encompassing the LBX1 gene. The LBX1 protein, a homeobox transcription factor that is orthologous to the Drosophila ladybird late gene, is involved in proper migration of muscle precursor cells, specification of cardiac neural crest cells, and neuronal determination in developing neural tubes. CONCLUSIONS: Our results firmly establish the LBX1 region as the first major susceptibility locus for AIS in Asian and non-Hispanic white groups, and provide a platform for larger studies in additional ancestral groups.

Genetic variants in GPR126 are associated with adolescent idiopathic scoliosis.

Adolescent idiopathic scoliosis (AIS) is the most common pediatric skeletal disease. We previously reported a locus on chromosome 10q24.31 associated with AIS susceptibility in Japanese using a genome-wide association study (GWAS) consisting of 1,033 cases and 1,473 controls. To identify additional AIS-associated loci, we expanded the study by adding X-chromosome SNPs in the GWAS and increasing the size of the replication cohorts. Through a stepwise association study including 1,819 cases and 25,939 controls, we identified a new susceptibility locus on chromosome 6q24.1 in Japanese (P = 2.25 × 10(-10); odds ratio (OR) = 1.28). The most significantly associated SNP, rs6570507, was in GPR126 (encoding G protein-coupled receptor 126). Its association was replicated in Han Chinese and European-ancestry populations (combined P = 1.27 × 10(-14); OR = 1.27). GPR126 was highly expressed in cartilage, and the knockdown of gpr126 in zebrafish caused delayed ossification of the developing spine. Our results should provide insights into the etiology and pathogenesis of AIS.

A novel method for analyzing genetic association with longitudinal phenotypes.

Knowledge of genes influencing longitudinal patterns may offer information about predicting disease progression. We developed a systematic procedure for testing association between SNP genotypes and longitudinal phenotypes. We evaluated false positive rates and statistical power to localize genes for disease progression. We used genome-wide SNP data from the Framingham Heart Study. With longitudinal data from two real studies unrelated to Framingham, we estimated three trajectory curves from each study. We performed simulations by randomly selecting 500 individuals. In each simulation replicate, we assigned each individual to one of the three trajectory groups based on the underlying hypothesis (null or alternative), and generated corresponding longitudinal data. Individual Bayesian posterior probabilities (BPPs) for belonging to a specific trajectory curve were estimated. These BPPs were treated as a quantitative trait and tested (using the Wald test) for genome-wide association. Empirical false positive rates and power were calculated. Our method maintained the expected false positive rate for all simulation models. Also, our method achieved high empirical power for most simulations. Our work presents a method for disease progression gene mapping. This method is potentially clinically significant as it may allow doctors to predict disease progression based on genotype and determine treatment accordingly.

IL-18R1 and IL-18RAP SNPs may be associated with bronchopulmonary dysplasia in African-American infants.

INTRODUCTION: The genetic contribution to the development of bronchopulmonary dysplasia (BPD) in prematurely born infants is substantial, but information related to the specific genes involved is lacking. RESULTS: Genotype analysis revealed, after multiple comparisons correction, two significant single-nucleotide polymorphism (SNPs), rs3771150 (IL-18RAP) and rs3771171 (IL-18R1), in African Americans (AAs) with BPD (vs. AAs without BPD; q < 0.05). No associations with Caucasian (CA) BPD, AA or CA respiratory distress syndrome (RDS), or prematurity in either AAs or CAs were identified with these SNPs. Respective frequencies were 0.098 and 0.093 in infants without BPD and 0.38 for each SNP in infants with BPD. In the replication set (82 cases; 102 controls), the P values were 0.012 for rs3771150 and 0.07 for rs3771171. Combining P values using Fisher's method, overall P values were 8.31 × 10(-7) for rs3771150 and 6.33 × 10(-6) for rs3771171. DISCUSSION: We conclude that IL-18RAP and IL-18R1 SNPs identify AA infants at risk for BPD. These genes may contribute to AA BPD pathogenesis via inflammatory-mediated processes and require further study. METHODS: We conducted a case-control SNP association study of candidate genes (n = 601) or 6,324 SNPs in 1,091 prematurely born infants with gestational age <35 weeks, with or without neonatal lung disease including BPD. BPD was defined as a need for oxygen at 28 days.

TDT-HET: a new transmission disequilibrium test that incorporates locus heterogeneity into the analysis of family-based association data.

BACKGROUND: Locus heterogeneity is one of the most documented phenomena in genetics. To date, relatively little work had been done on the development of methods to address locus heterogeneity in genetic association analysis. Motivated by Zhou and Pan's work, we present a mixture model of linked and unlinked trios and develop a statistical method to estimate the probability that a heterozygous parent transmits the disease allele at a di-allelic locus, and the probability that any trio is in the linked group. The purpose here is the development of a test that extends the classic transmission disequilibrium test (TDT) to one that accounts for locus heterogeneity. RESULTS: Our simulations suggest that, for sufficiently large sample size (1000 trios) our method has good power to detect association even the proportion of unlinked trios is high (75%). While the median difference (TDT-HET empirical power - TDT empirical power) is approximately 0 for all MOI, there are parameter settings for which the power difference can be substantial. Our multi-locus simulations suggest that our method has good power to detect association as long as the markers are reasonably well-correlated and the genotype relative risk are larger. Results of both single-locus and multi-locus simulations suggest our method maintains the correct type I error rate.Finally, the TDT-HET statistic shows highly significant p-values for most of the idiopathic scoliosis candidate loci, and for some loci, the estimated proportion of unlinked trios approaches or exceeds 50%, suggesting the presence of locus heterogeneity. CONCLUSIONS: We have developed an extension of the TDT statistic (TDT-HET) that allows for locus heterogeneity among coded trios. Benefits of our method include: estimates of parameters in the presence of heterogeneity, and reasonable power even when the proportion of linked trios is small. Also, we have extended multi-locus methods to TDT-HET and have demonstrated that the empirical power may be high to detect linkage. Last, given that we obtain PPBs, we conjecture that the TDT-HET may be a useful method for correctly identifying linked trios. We anticipate that researchers will find this property increasingly useful as they apply next-generation sequencing data in family based studies.

Single-variant and multi-variant trend tests for genetic association with next-generation sequencing that are robust to sequencing error.

As with any new technology, next-generation sequencing (NGS) has potential advantages and potential challenges. One advantage is the identification of multiple causal variants for disease that might otherwise be missed by SNP-chip technology. One potential challenge is misclassification error (as with any emerging technology) and the issue of power loss due to multiple testing. Here, we develop an extension of the linear trend test for association that incorporates differential misclassification error and may be applied to any number of SNPs. We call the statistic the linear trend test allowing for error, applied to NGS, or LTTae,NGS. This statistic allows for differential misclassification. The observed data are phenotypes for unrelated cases and controls, coverage, and the number of putative causal variants for every individual at all SNPs. We simulate data considering multiple factors (disease mode of inheritance, genotype relative risk, causal variant frequency, sequence error rate in cases, sequence error rate in controls, number of loci, and others) and evaluate type I error rate and power for each vector of factor settings. We compare our results with two recently published NGS statistics. Also, we create a fictitious disease model based on downloaded 1000 Genomes data for 5 SNPs and 388 individuals, and apply our statistic to those data. We find that the LTTae,NGS maintains the correct type I error rate in all simulations (differential and non-differential error), while the other statistics show large inflation in type I error for lower coverage. Power for all three methods is approximately the same for all three statistics in the presence of non-differential error. Application of our statistic to the 1000 Genomes data suggests that, for the data downloaded, there is a 1.5% sequence misclassification rate over all SNPs. Finally, application of the multi-variant form of LTTae,NGS shows high power for a number of simulation settings, although it can have lower power than the corresponding single-variant simulation results, most probably due to our specification of multi-variant SNP correlation values. In conclusion, our LTTae,NGS addresses two key challenges with NGS disease studies; first, it allows for differential misclassification when computing the statistic; and second, it addresses the multiple-testing issue in that there is a multi-variant form of the statistic that has only one degree of freedom, and provides a single p value, no matter how many loci.

Genome-wide association studies of adolescent idiopathic scoliosis suggest candidate susceptibility genes.

Adolescent idiopathic scoliosis (AIS) is an unexplained and common spinal deformity seen in otherwise healthy children. Its pathophysiology is poorly understood despite intensive investigation. Although genetic underpinnings are clear, replicated susceptibility loci that could provide insight into etiology have not been forthcoming. To address these issues, we performed genome-wide association studies (GWAS) of ∼327 000 single nucleotide polymorphisms (SNPs) in 419 AIS families. We found strongest evidence of association with chromosome 3p26.3 SNPs in the proximity of the CHL1 gene (P < 8 × 10(-8) for rs1400180). We genotyped additional chromosome 3p26.3 SNPs and tested replication in two follow-up case-control cohorts, obtaining strongest results when all three cohorts were combined (rs10510181 odds ratio = 1.49, 95% confidence interval = 1.29-1.73, P = 2.58 × 10(-8)), but these were not confirmed in a separate GWAS. CHL1 is of interest, as it encodes an axon guidance protein related to Robo3. Mutations in the Robo3 protein cause horizontal gaze palsy with progressive scoliosis (HGPPS), a rare disease marked by severe scoliosis. Other top associations in our GWAS were with SNPs in the DSCAM gene encoding an axon guidance protein in the same structural class with Chl1 and Robo3. We additionally found AIS associations with loci in CNTNAP2, supporting a previous study linking this gene with AIS. Cntnap2 is also of functional interest, as it interacts directly with L1 and Robo class proteins and participates in axon pathfinding. Our results suggest the relevance of axon guidance pathways in AIS susceptibility, although these findings require further study, particularly given the apparent genetic heterogeneity in this disease.

A cost-effective statistical method to correct for differential genotype misclassification when performing case-control genetic association.

BACKGROUND/AIMS: There is a growing interest regarding the effect of differential misclassification on power and type I error rate in genome-wide association studies. We present an extension of a previously published test statistic: the likelihood ratio test allowing for errors (LRTAE). This test uses double-sample information on a subset of individuals to increase power for genetic association in the presence of nondifferential misclassification. METHODS: We extend the original LRTAE by allowing for differential genotype misclassification between case and control populations. We label this new statistic as LRT(D)A(M)E . We test the performance of this statistic with data simulated under differential misclassification specifications and two different types of genetic models: null and power. For simulations using the null model, we specify that there is no difference between case and control genotype frequencies before the introduction of errors. For simulations under power, we consider three modes of inheritance: dominant, multiplicative, and recessive. RESULTS: We show that the LRT(D)A(M)E , with p values computed using permutation, maintains a correct type I error rate under the null model after the introduction of differential genotyping errors. Also, we find that as little as 10 to 15% of double-sampled genotype data is needed to achieve this effect. Aside from a few situations (particularly recessive mode of inheritance simulations) the LRT(D)A(M)E version that calculates p values through permutation requires 15 to 20% double sampling to maintain an 80% power for a 0.05 significance level and approximately 20% double sampling for a 0.01 significance level.

Evidence for association of two variants of the nociceptin/orphanin FQ receptor gene OPRL1 with vulnerability to develop opiate addiction in Caucasians.

OBJECTIVES: The OPRL1 gene encodes the nociceptin/orphanin FQ receptor, which plays a role in regulating tolerance and behavioral responses to morphine. However, there is limited information on whether variants of OPRL1 are associated with vulnerability to develop opiate addiction. In this study, we examined five variants of OPRL1 and their role in determining vulnerability to develop opiate addiction. METHODS: We recruited 447 individuals: 271 former severe heroin addicts and 176 healthy controls. Using a 5'-fluorogenic exonuclease assay, we genotyped individuals at five variants in OPRL1. It was then determined whether there was a significant association of allele, genotype, or haplotype frequency with vulnerability to develop opiate addiction. RESULTS: When the cohort was stratified by ethnicity, we found that, in Caucasians but not in African-Americans or Hispanics, the allele frequency of rs6090041 and rs6090043 were significantly associated point-wise with opiate addiction (P=0.03 and 0.04, respectively). Of the haplotypes formed by these two variants, one haplotype was found to be associated with protection from developing opiate addiction in both African-Americans (point-wise P=0.04) and Caucasians (point-wise P=0.04), and another haplotype with vulnerability to develop opiate addiction in Caucasians only (P=0.020). CONCLUSION: This study provides evidence for an association of two variants of the OPRL1 gene, rs6090041 and rs6090043, with vulnerability to develop opiate addiction, suggesting a role for nociceptin/orphanin FQ receptor in the development of opiate addiction.

EDNRA variants associate with smooth muscle mRNA levels, cell proliferation rates, and cystic fibrosis pulmonary disease severity.

Airway inflammation and pulmonary disease are heterogeneous phenotypes in cystic fibrosis (CF) patients, even among patients with the same cystic fibrosis transmembrane conductance regulator (CFTR) genotype. Endothelin, a proinflammatory peptide and smooth muscle agonist, is increased in CF airways, potentially contributing to the pulmonary phenotype. Four cohorts of CF patients were screened for variants in endothelin pathway genes to determine whether any of these variants associated with pulmonary function. An initial cohort of 808 CF patients homozygous for the common CF mutation, DeltaF508, showed significant association for polymorphisms in the endothelin receptor A gene, EDNRA (P = 0.04), but not in the related endothelin genes (EDN1, EDN2, EDN3, or EDNRB) or NOS1, NOS2A, or NOS3. Variants within EDNRA were examined in three additional cohorts of CF patients, 238 patients from Seattle, WA, 303 from Ireland and the U.K., and 228 from Cleveland, OH, for a total of 1,577 CF patients. The three additional groups each demonstrated a significant association between EDNRA 3'-untranslated region (UTR) variant rs5335 and pulmonary function (P = 0.002). At the molecular level, single nucleotide primer extension assays suggest that the effect of the variants is quantitative. EDNRA mRNA levels from cultured primary tracheal smooth muscle cells are greater for the allele that appears to be deleterious to lung function than for the protective allele, suggesting a mechanism by which increased receptor function is harmful to the CF airway. Finally, cell proliferation studies using human airway smooth muscle cells demonstrated that cells homozygous for the deleterious allele proliferate at a faster rate than those homozygous for the protective allele.

Genomewide scan for linkage reveals evidence of several susceptibility loci for alopecia areata.

Alopecia areata (AA) is a genetically determined, immune-mediated disorder of the hair follicle that affects 1%-2% of the U.S. population. It is defined by a spectrum of severity that ranges from patchy localized hair loss on the scalp to the complete absence of hair everywhere on the body. In an effort to define the genetic basis of AA, we performed a genomewide search for linkage in 20 families with AA consisting of 102 affected and 118 unaffected individuals from the United States and Israel. Our analysis revealed evidence of at least four susceptibility loci on chromosomes 6, 10, 16 and 18, by use of several different statistical approaches. Fine-mapping analysis with additional families yielded a maximum multipoint LOD score of 3.93 on chromosome 18, a two-point affected sib pair (ASP) LOD score of 3.11 on chromosome 16, several ASP LOD scores >2.00 on chromosome 6q, and a haplotype-based relative risk LOD of 2.00 on chromosome 6p (in the major histocompatibility complex locus). Our findings confirm previous studies of association of the human leukocyte antigen locus with human AA, as well as the C3H-HeJ mouse model for AA. Interestingly, the major loci on chromosomes 16 and 18 coincide with loci for psoriasis reported elsewhere. These results suggest that these regions may harbor gene(s) involved in a number of different skin and hair disorders.

Density-based clustering in haplotype analysis for association mapping.

Clustering of related haplotypes in haplotype-based association mapping has the potential to improve power by reducing the degrees of freedom without sacrificing important information about the underlying genetic structure. We have modified a generalized linear model approach for association analysis by incorporating a density-based clustering algorithm to reduce the number of coefficients in the model. Using the GAW 15 Problem 3 simulated data, we show that our novel method can substantially enhance power to detect association with the binary rheumatoid arthritis (RA) phenotype at the HLA-DRB1 locus on chromosome 6. In contrast, clustering did not appreciably improve performance at locus D, perhaps a consequence of a rare susceptibility allele and of the overwhelming effect of HLA-DRB1/locus C, 5 cM distal. Optimization of parameters governing the clustering algorithm identified a set of parameters that delivered nearly ideal performance in a variety of situations. The cluster-based score test was valid over a wide range of haplotype diversity, and was robust to severe departures from Hardy-Weinberg equilibrium encountered near HLA-DRB1 in RA case-control samples.

Beta 2 adrenergic receptor polymorphisms in cystic fibrosis.

There has been a recent emphasis on identifying modifier genes that influence the severity of cystic fibrosis (CF) lung disease. The beta-2-adrenergic receptor is expressed on airway smooth muscle, is the target for inhaled beta agonists, and has several common polymorphisms in its gene, ADRB2. Polymorphisms changing glycine to arginine or glutamate to glutamine in codons 16 and 27, respectively, were associated with differences in clinical response to inhaled beta agonists in individuals with asthma. We compared acute airway responsiveness and 5-year decline in pulmonary function in CF patients with different ADRB2 genotypes. One hundred and six subjects performed spirometry before and after the administration of an inhaled bronchodilator, and had ADRB2 genotype determined for codons 16 and 27. Comparing the percent change in FEV(1) and FEF(25-75) continuously revealed differences in the degree of airway responsiveness to bronchodilator between ADRB2-genotyped groups. However, there was no significant relationship between the ADRB2 genotype at positions 16 and 27 and bronchodilator response when defined as 12% improvement in FEV(1). Five-year decline in percent predicted FEV(1) showed no association with ADRB2 genotype. These data are consistent with variants of the ADRB2 gene having different responses to bronchodilator, but the long-term effects, if any, are not apparent over a 5-year period.