Novel MicroRNA-Regulated Transcript Networks Are Associated with Chemotherapy Response in Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) is a highly lethal gynecologic cancer, in part due to resistance to platinum-based chemotherapy reported among 20% of patients. This study aims to generate novel hypotheses of the biological mechanisms underlying chemotherapy resistance, which remain poorly understood. Differential expression analyses of mRNA- and microRNA-sequencing data from HGSOC patients of The Cancer Genome Atlas identified 21 microRNAs associated with angiogenesis and 196 mRNAs enriched for adaptive immunity and translation. Coexpression network analysis identified three microRNA networks associated with chemotherapy response enriched for lipoprotein transport and oncogenic pathways, as well as two mRNA networks enriched for ubiquitination and lipid metabolism. These network modules were replicated in two independent ovarian cancer cohorts. Moreover, integrative analyses of the mRNA/microRNA sequencing and single-nucleotide polymorphisms (SNPs) revealed potential regulation of significant mRNA transcripts by microRNAs and SNPs (expression quantitative trait loci). Thus, we report novel transcriptional networks and biological pathways associated with resistance to platinum-based chemotherapy in HGSOC patients. These results expand our understanding of the effector networks and regulators of chemotherapy response, which will help to improve the management of ovarian cancer.

Investigating the relationship between melatonin patterns and methylation in circadian genes among day shift and night shift workers.

OBJECTIVES: Mechanisms underlying the carcinogenicity of night shift work remain uncertain. One compelling yet understudied cancer mechanism may involve altered DNA methylation in circadian genes due to melatonin secretion patterns. The objective of this study was to explore the relationship between melatonin secretion patterns and circadian gene methylation among day and night shift workers. METHODS: Female healthcare employees (n=38 day workers, n=36 night shift workers) for whom we had urinary 6-sulfatoxymelatonin secretion data from a previous study were recontacted. New blood samples were collected and used to measure methylation levels at 1150 CpG loci across 22 circadian genes using the Illumina Infinium MethylationEPIC beadchip. Linear regression was used to examine the association between melatonin (acrophase and mesor) and M values for each CpG site (false discovery rate, q=0.2), while testing for effect modification by shift work status. RESULTS: Among night shift workers, a higher mesor (24 hours of mean production of melatonin) was associated with increased methylation in the body of RORA (q=0.02) and decreased methylation in the putative promoter region of MTNR1A (q=0.03). Later acrophase (ie, time of peak concentration) was associated with increased methylation in the putative promoter region of MTNR1A (q=0.20) and decreased methylation in the body of PER3 (q=0.20). No associations were identified among day workers. CONCLUSIONS: In conclusion, patterns in melatonin secretion were associated with differential circadian gene methylation among night shift workers. Melatonin and alteration of DNA methylation in circadian genes may be one pathway towards increased cancer risk, although larger-scale studies examining multiple time points are needed.

DNA methylation of circadian genes and markers of cardiometabolic risk in female hospital workers: An exploratory study.

Night shift work has been linked to increased risk of cardiovascular disease (CVD); however, the underlying mechanisms remain unclear. A compelling yet understudied mechanism involves differential DNA methylation of circadian genes. To investigate the relevance of this mechanism, we conducted an exploratory cross-sectional study of 74 female hospital personnel (38 day workers, 36 night shift workers). Sociodemographic, lifestyle, and health characteristics as well as shift work status and history were determined through self-report. Fasting blood samples were collected to measure markers of cardiometabolic risk and DNA was extracted to measure DNA methylation of 1150 cytosine-guanine (CpG) sites across 22 circadian genes. Associations between methylation levels at individual CpG sites (ß-values) and markers of cardiometabolic risk were analyzed while considering effect modification by shift work status. The false discovery rate was applied to account for multiple comparisons (q ≤ 0.20). Two CpG sites [cg06758649 (CRY1) and cg06899802 (CSNK1A1)] were differentially associated with waist circumference and body mass index by shift work status, and eight CpG sites [cg26103512 (CSNK1D), cg03941313 (CSNK1E), cg18217763 (CSNK1E), cg16682686 (DEC1), cg12061096 (RORA), cg10133825 (RORA), cg19652148 (RORA), and cg22904654 (RORA)] were differentially associated with LDL cholesterol concentration by shift work status (all q ≤ 0.20). Our findings suggest that the relationship between DNA methylation of circadian genes and cardiometabolic risk differs by day and night shift worker status, which may contribute to mechanisms of increased risk of CVD observed among night shift workers.

Exploring the impact of night shift work on methylation of circadian genes.

Night shift work is associated with increased breast cancer risk, but the molecular mechanisms are not well-understood. The objective of this study was to explore the relationship between night shift work parameters (current status, duration/years, and intensity) and methylation in circadian genes as a potential mechanism underlying the carcinogenic effects of night shift work. A cross-sectional study was conducted among 74 female healthcare employees (n = 38 day workers, n = 36 night shift workers). The Illumina Infinium MethylationEPIC beadchip was applied to DNA extracted from blood samples to measure methylation using a candidate gene approach at 1150 CpG loci across 22 circadian genes. Linear regression models were used to examine the association between night shift work parameters and continuous methylation measurements (ß-values) for each CpG site. The false-discovery rate (q = 0.2) was used to account for multiple comparisons. Compared to day workers, current night shift workers demonstrated hypermethylation in the 5'UTR region of CSNK1E (q = 0.15). Individuals that worked night shifts for ≥10 years exhibited hypomethylation in the gene body of NR1D1 (q = 0.08) compared to those that worked <10 years. Hypermethylation in the gene body of ARNTL was also apparent in those who worked ≥3 consecutive night shifts a week (q = 0.18). These findings suggest that night shift work is associated with differential methylation in core circadian genes, including CSNK1E, NR1D1 and ARNTL. Future, larger-scale studies with long-term follow-up and detailed night shift work assessment are needed to confirm and expand on these findings.

Gene networks and expression quantitative trait loci associated with adjuvant chemotherapy response in high-grade serous ovarian cancer.

BACKGROUND: A major impediment in the treatment of ovarian cancer is the relapse of chemotherapy-resistant tumors, which occurs in approximately 25% of patients. A better understanding of the biological mechanisms underlying chemotherapy resistance will improve treatment efficacy through genetic testing and novel therapies. METHODS: Using data from high-grade serous ovarian carcinoma (HGSOC) patients in the Cancer Genome Atlas (TCGA), we classified those who remained progression-free for 12 months following platinum-taxane combination chemotherapy as "chemo-sensitive" (N = 160) and those who had recurrence within 6 months as "chemo-resistant" (N = 110). Univariate and multivariate analysis of expression microarray data were used to identify differentially expressed genes and co-expression gene networks associated with chemotherapy response. Moreover, we integrated genomics data to determine expression quantitative trait loci (eQTL). RESULTS: Differential expression of the Valosin-containing protein (VCP) gene and five co-expression gene networks were significantly associated with chemotherapy response in HGSOC. VCP and the most significant co-expression network module contribute to protein processing in the endoplasmic reticulum, which has been implicated in chemotherapy response. Both univariate and multivariate analysis findings were successfully replicated in an independent ovarian cancer cohort. Furthermore, we identified 192 cis-eQTLs associated with the expression of network genes and 4 cis-eQTLs associated with BRCA2 expression. CONCLUSION: This study implicates both known and novel genes as well as biological processes underlying response to platinum-taxane-based chemotherapy among HGSOC patients.

Human milk fatty acid composition is associated with dietary, genetic, sociodemographic, and environmental factors in the CHILD Cohort Study.

BACKGROUND: Fatty acids are a vital component of human milk. They influence infant neurodevelopment and immune function, and they provide ∼50% of milk's energy content. OBJECTIVES: The objectives of this study were to characterize the composition of human milk fatty acids in a large Canadian birth cohort and identify factors influencing their variability. METHODS: In a subset of the CHILD cohort (n = 1094), we analyzed milk fatty acids at 3-4 mo postpartum using GLC. Individual and total SFAs, MUFAs, and n-3 and n-6 PUFAs were analyzed using SD scores and principal component analysis (PCA). Maternal diet, sociodemographic, health, and environmental factors were self-reported. Single-nucleotide polymorphisms were assessed in the fatty acid desaturase 1 (FADS1-rs174556) and 2 (FADS2-rs174575) genes. RESULTS: Fatty acid profiles were variable, with individual fatty acid proportions varying from 2- to >30-fold between women. Using PCA, we identified 4 milk fatty acid patterns: "MUFA and low SFA," "high n-6 PUFA," "high n-3 PUFA," and "high medium-chain fatty acids." In multivariable-adjusted analyses, fish oil supplementation and fatty cold water fish intake were positively associated with DHA and the "high n-3 PUFA" pattern. Mothers carrying the minor allele of FADS1-rs174556 had lower proportions of arachidonic acid (ARA; 20:4n-6). Independent of selected dietary variables and genetic variants, Asian ethnicity was associated with higher linoleic acid (18:2n-6) and total n-3 PUFAs. Ethnic differences in ARA were explained by FADS1 genotype. Maternal obesity was independently associated with higher total SFAs, the "high medium-chain fatty acid" pattern, and lower total MUFAs. Lactation stage, season, study site, and maternal education were also independently associated with some milk fatty acids. No associations were observed for maternal age, parity, delivery mode, or infant sex. CONCLUSIONS: This study provides unique insights about the "normal" variation in the composition of human milk fatty acids and the contributing dietary, genetic, sociodemographic, health, and environmental factors. Further research is required to assess implications for infant health.

Integrated Analysis of Human Milk Microbiota With Oligosaccharides and Fatty Acids in the CHILD Cohort.

Background: Human milk contains many bioactive components that are typically studied in isolation, including bacteria. We performed an integrated analysis of human milk oligosaccharides and fatty acids to explore their associations with milk microbiota. Methods: We studied a sub-sample of 393 mothers in the CHILD birth cohort. Milk was collected at 3-4 months postpartum. Microbiota was analyzed by 16S rRNA gene V4 sequencing. Oligosaccharides and fatty acids were analyzed by rapid high-throughput high performance and gas liquid chromatography, respectively. Dimension reduction was performed with principal component analysis for oligosaccharides and fatty acids. Center log-ratio transformation was applied to all three components. Associations between components were assessed using Spearman rank correlation, network visualization, multivariable linear regression, redundancy analysis, and structural equation modeling. P-values were adjusted for multiple comparisons. Key covariates were considered, including fucosyltransferase-2 (FUT2) secretor status of mother and infant, method of feeding (direct breastfeeding or pumped breast milk), and maternal fish oil supplement use. Results: Overall, correlations were strongest between milk components of the same type. For example, FUT2-dependent HMOs were positively correlated with each other, and Staphylococcus was negatively correlated with other core taxa. Some associations were also observed between components of different types. Using redundancy analysis and structural equation modeling, the overall milk fatty acid profile was significantly associated with milk microbiota composition. In addition, some individual fatty acids [22:6n3 (docosahexaenoic acid), 22:5n3, 20:5n3, 17:0, 18:0] and oligosaccharides (fucosyl-lacto-N-hexaose, lacto-N-hexaose, lacto-N-fucopentaose I) were associated with microbiota α diversity, while others (C18:0, 3'-sialyllactose, disialyl-lacto-N-tetraose) were associated with overall microbiota composition. Only a few significant associations between individual HMOs and microbiota were observed; notably, among mothers using breast pumps, Bifidobacterium prevalence was associated with lower abundances of disialyl-lacto-N-hexaose. Additionally, among non-secretor mothers, Staphylococcus was positively correlated with sialylated HMOs. Conclusion: Using multiple approaches to integrate and analyse milk microbiota, oligosaccharides, and fatty acids, we observed several associations between different milk components and microbiota, some of which were modified by secretor status and/or breastfeeding practices. Additional research is needed to further validate and mechanistically characterize these associations and determine their relevance to infant gut and respiratory microbiota development and health.

Whole genome sequencing identifies high-impact variants in well-known pharmacogenomic genes.

More than 1100 genetic loci have been correlated with drug response outcomes but disproportionately few have been translated into clinical practice. One explanation for the low rate of clinical implementation is that the majority of associated variants may be in linkage disequilibrium (LD) with the causal variants, which are often elusive. This study aims to identify and characterize likely causal variants within well-established pharmacogenomic genes using next-generation sequencing data from the 1000 Genomes Project. We identified 69,319 genetic variations within 160 pharmacogenomic genes, of which 8207 variants are in strong LD (r2>0.8) with known pharmacogenomic variants. Of the latter, eight are coding or structural variants predicted to have high impact, with 19 additional missense variants that are predicted to have moderate impact. In conclusion, we identified putatively functional variants within known pharmacogenomics loci that could account for the association signals and represent the missing causative variants underlying drug response phenotypes.

Identification of deleterious and regulatory genomic variations in known asthma loci.

BACKGROUND: Candidate gene and genome-wide association studies have identified hundreds of asthma risk loci. The majority of associated variants, however, are not known to have any biological function and are believed to represent markers rather than true causative mutations. We hypothesized that many of these associated markers are in linkage disequilibrium (LD) with the elusive causative variants. METHODS: We compiled a comprehensive list of 449 asthma-associated variants previously reported in candidate gene and genome-wide association studies. Next, we identified all sequence variants located within the 305 unique genes using whole-genome sequencing data from the 1000 Genomes Project. Then, we calculated the LD between known asthma variants and the sequence variants within each gene. LD variants identified were then annotated to determine those that are potentially deleterious and/or functional (i.e. coding or regulatory effects on the encoded transcript or protein). RESULTS: We identified 10,130 variants in LD (r2 > 0.6) with known asthma variants. Annotations of these LD variants revealed that several have potentially deleterious effects including frameshift, alternate splice site, stop-lost, and missense. Moreover, 24 of the LD variants have been reported to regulate gene expression as expression quantitative trait loci (eQTLs). CONCLUSIONS: This study is proof of concept that many of the genetic loci previously associated with complex diseases such as asthma are not causative but represent markers of disease, which are in LD with the elusive causative variants. We hereby report a number of potentially deleterious and regulatory variants that are in LD with the reported asthma loci. These reported LD variants could account for the original association signals with asthma and represent the true causative mutations at these loci.

Stress and Bronchodilator Response in Children with Asthma.

RATIONALE: Stress is associated with asthma morbidity in Puerto Ricans (PRs), who have reduced bronchodilator response (BDR). OBJECTIVES: To examine whether stress and/or a gene regulating anxiety (ADCYAP1R1) is associated with BDR in PR and non-PR children with asthma. METHODS: This was a cross-sectional study of stress and BDR (percent change in FEV1 after BD) in 234 PRs ages 9-14 years with asthma. We assessed child stress using the Checklist of Children's Distress Symptoms, and maternal stress using the Perceived Stress Scale. Replication analyses were conducted in two cohorts. Polymorphisms in ADCYAP1R1 were genotyped in our study and six replication studies. Multivariable models of stress and BDR were adjusted for age, sex, income, environmental tobacco smoke, and use of inhaled corticosteroids. MEASUREMENTS AND MAIN RESULTS: High child stress was associated with reduced BDR in three cohorts. PR children who were highly stressed (upper quartile, Checklist of Children's Distress Symptoms) and whose mothers had high stress (upper quartile, Perceived Stress Scale) had a BDR that was 10.2% (95% confidence interval, 6.1-14.2%) lower than children who had neither high stress nor a highly stressed mother. A polymorphism in ADCYAP1R1 (rs34548976) was associated with reduced BDR. This single-nucleotide polymorphism is associated with reduced expression of the gene for the ß2-adrenergic receptor (ADRB2) in CD4(+) lymphocytes of subjects with asthma, and it affects brain connectivity of the amygdala and the insula (a biomarker of anxiety). CONCLUSIONS: High child stress and an ADCYAP1R1 single-nucleotide polymorphism are associated with reduced BDR in children with asthma. This is likely caused by down-regulation of ADRB2 in highly stressed children.

Pharmacogenomics: novel loci identification via integrating gene differential analysis and eQTL analysis.

Nearly one-half of asthmatic patients do not respond to the most commonly prescribed controller therapy, inhaled corticosteroids (ICS). We conducted an expression quantitative trait loci (eQTL) analysis using >300 expression microarrays (from 117 lymphoblastoid cell lines) in corticosteroid (dexamethasone) treated and untreated cells derived from asthmatic subjects in the Childhood Asthma Management Program (CAMP) clinical trial. We then tested the associations of eQTL with longitudinal change in airway responsiveness to methacholine (LnPC20) on ICS. We identified 2484 cis-eQTL affecting 767 genes following dexamethasone treatment. A significant over-representation of lnPC20-associated cis-eQTL [190 single-nucleotide polymorphisms (SNPs)] among differentially expressed genes (odds ratio = 1.76, 95% confidence interval: 1.35-2.29) was noted in CAMP Caucasians. Forty-six of these 190 clinical associations were replicated in CAMP African Americans, including seven SNPs near six genes meeting criteria for genome-wide significance (P < 2 × 10(-7)). Notably, the majority of genome-wide findings would not have been uncovered via analysis of untreated samples. These results indicate that identifying eQTL after relevant environmental perturbation enables identification of true pharmacogenetic variants.

Genetic predictors associated with improvement of asthma symptoms in response to inhaled corticosteroids.

BACKGROUND: To date, genome-wide association studies (GWASs) of inhaled corticosteroid (ICS) response in asthmatic patients have focused primarily on lung function and exacerbations. OBJECTIVE: We hypothesized that GWAS analysis could identify novel genetic markers predicting a symptomatic response to ICSs. METHODS: We analyzed differences in asthma symptoms in response to ICSs in 124 white children from the Childhood Asthma Management Program (CAMP) trial using scores from diary cards. Of the 440,862 single nucleotide polymorphisms (SNPs) analyzed, the top 100 ranked SNPs were pursued for replication initially in subjects from the pediatric Childhood Asthma Research and Education trials (77 white children) and then in subjects from the adult Asthma Clinical Research Network (110 white adults) and Leukotriene Modifier or Corticosteroid or Corticosteroid-Salmeterol trials (110 white adults). RESULTS: The lowest P value for GWAS analysis in the CAMP trial was 8.94 × 10(-8) (rs2388639). Of the 60 SNPs available in the Childhood Asthma Research and Education Network trials, rs1558726 (combined P = 1.02 × 10(-5)), rs2388639 (combined P = 8.56 × 10(-9)), and rs10044254 (combined P = 9.16 × 10(-8)) independently replicated. However, these 3 SNPs were not additionally replicated in the adult asthmatic patients of the remaining trials. rs10044254 lies in the intronic region of F-box and leucine-rich repeat protein 7 (FBXL7) and is associated with decreased expression in immortalized B cells derived from CAMP participants. CONCLUSIONS: We have identified a novel SNP, rs10044254, associated with both decreased expression of FBXL7 and improved symptomatic response to ICSs in 2 independent pediatric cohorts. Our results suggest that there might be a specific genetic mechanism regulating symptomatic response to ICSs in children that does not carry over to adults.

Genomewide association between GLCCI1 and response to glucocorticoid therapy in asthma.

BACKGROUND: The response to treatment for asthma is characterized by wide interindividual variability, with a significant number of patients who have no response. We hypothesized that a genomewide association study would reveal novel pharmacogenetic determinants of the response to inhaled glucocorticoids. METHODS: We analyzed a small number of statistically powerful variants selected on the basis of a family-based screening algorithm from among 534,290 single-nucleotide polymorphisms (SNPs) to determine changes in lung function in response to inhaled glucocorticoids. A significant, replicated association was found, and we characterized its functional effects. RESULTS: We identified a significant pharmacogenetic association at SNP rs37972, replicated in four independent populations totaling 935 persons (P=0.0007), which maps to the glucocorticoid-induced transcript 1 gene (GLCCI1) and is in complete linkage disequilibrium (i.e., perfectly correlated) with rs37973. Both rs37972 and rs37973 are associated with decrements in GLCCI1 expression. In isolated cell systems, the rs37973 variant is associated with significantly decreased luciferase reporter activity. Pooled data from treatment trials indicate reduced lung function in response to inhaled glucocorticoids in subjects with the variant allele (P=0.0007 for pooled data). Overall, the mean (±SE) increase in forced expiratory volume in 1 second in the treated subjects who were homozygous for the mutant rs37973 allele was only about one third of that seen in similarly treated subjects who were homozygous for the wild-type allele (3.2±1.6% vs. 9.4±1.1%), and their risk of a poor response was significantly higher (odds ratio, 2.36; 95% confidence interval, 1.27 to 4.41), with genotype accounting for about 6.6% of overall inhaled glucocorticoid response variability. CONCLUSIONS: A functional GLCCI1 variant is associated with substantial decrements in the response to inhaled glucocorticoids in patients with asthma.

A functional XPNPEP2 promoter haplotype leads to reduced plasma aminopeptidase P and increased risk of ACE inhibitor-induced angioedema.

Angiotensin I-converting enzyme inhibitors (ACEi) are widely used antihypertensive agents that are associated with a potentially life-threatening reaction, ACEi-angioedema. Impaired metabolism of bradykinin and des-Arg(9) -bradykinin by aminopeptidase P (APP) is a key contributor to ACEi-angioedema. This study aimed to characterize the genetic regulation of the XPNPEP2 gene and identify the genetic factors contributing to variance in plasma APP activity and ACEi-angioedema. Additive genetic factors accounted for 47.3% of variance in plasma APP activity in healthy individuals. Nested deletion analysis identified the minimal promoter (-338 bp to -147 bp) and an enhancer region (-2,502 bp to -2,238 bp). Three polymorphisms (c.-2399C>A, c.-1612G>T, and c.-393G>A) were significantly associated with plasma APP activity. Haplotype ATG was significantly associated with reduced reporter gene activity and with reduced plasma APP activity. The c.-2399C>A polymorphism was located in an enhancer region and was predicted to differentially bind hepatic nuclear factor 4 (HNF4). Over expression of HNF4 increased the activation of haplotype ATG compared with haplotype CGG. In a case control study of subjects with a history of ACEi-angioedema, haplotype ATG was significantly associated with ACEi-angioedema (OR 4.87 [1.78-13.35] P = 0.002). The ATG haplotype is functional and contributes to ACEi-angioedema through a reduction in APP.

Global analysis of the impact of environmental perturbation on cis-regulation of gene expression.

Genetic variants altering cis-regulation of normal gene expression (cis-eQTLs) have been extensively mapped in human cells and tissues, but the extent by which controlled, environmental perturbation influences cis-eQTLs is unclear. We carried out large-scale induction experiments using primary human bone cells derived from unrelated donors of Swedish origin treated with 18 different stimuli (7 treatments and 2 controls, each assessed at 2 time points). The treatments with the largest impact on the transcriptome, verified on two independent expression arrays, included BMP-2 (t = 2h), dexamethasone (DEX) (t = 24 h), and PGE2 (t = 24 h). Using these treatments and control, we performed expression profiling for 18,144 RefSeq transcripts on biological replicates of the complete study cohort of 113 individuals (n(total) = 782) and combined it with genome-wide SNP-genotyping data in order to map treatment-specific cis-eQTLs (defined as SNPs located within the gene ± 250 kb). We found that 93% of cis-eQTLs at 1% FDR were observed in at least one additional treatment, and in fact, on average, only 1.4% of the cis-eQTLs were considered as treatment-specific at high confidence. The relative invariability of cis-regulation following perturbation was reiterated independently by genome-wide allelic expression tests where only a small proportion of variance could be attributed to treatment. Treatment-specific cis-regulatory effects were, however, 2- to 6-fold more abundant among differently expressed genes upon treatment. We further followed-up and validated the DEX-specific cis-regulation of the MYO6 and TNC loci and found top cis-regulatory variants located 180 kb and 250 kb upstream of the transcription start sites, respectively. Our results suggest that, as opposed to tissue-specificity of cis-eQTLs, the interactions between cellular environment and cis-variants are relatively rare (∼1.5%), but that detection of such specific interactions can be achieved by a combination of functional genomic approaches as described here.

Regulatory haplotypes in ARG1 are associated with altered bronchodilator response.

RATIONALE: ß2-agonists, the most common treatment for asthma, have a wide interindividual variability in response, which is partially attributed to genetic factors. We previously identified single nucleotide polymorphisms in the arginase 1 (ARG1) gene, which are associated with ß2-agonist bronchodilator response (BDR). OBJECTIVES: To identify cis-acting haplotypes in the ARG1 locus that are associated with BDR in patients with asthma and regulate gene expression in vitro. METHODS: We resequenced ARG1 in 96 individuals and identified three common, 5' haplotypes (denoted 1, 2, and 3). A haplotype-based association analysis of BDR was performed in three independent, adult asthma drug trial populations. Next, each haplotype was cloned into vectors containing a luciferase reporter gene and transfected into human airway epithelial cells (BEAS-2B) to ascertain its effect on gene expression. MEASUREMENTS AND MAIN RESULTS: BDR varied by haplotype in each of the three populations with asthma. Individuals with haplotype 1 were more likely to have higher BDR, compared to those with haplotypes 2 and 3, which is supported by odds ratios of 1.25 (95% confidence interval, 1.03-1.71) and 2.18 (95% confidence interval, 1.34-2.52), respectively. Luciferase expression was 50% greater in cells transfected with haplotype 1 compared to haplotypes 2 and 3. CONCLUSIONS: The identified ARG1 haplotypes seem to alter BDR and differentially regulate gene expression with a concordance of decreased BDR and reporter activity from haplotypes 2 and 3. These findings may facilitate pharmacogenetic tests to predict individuals who may benefit from other therapeutic agents in addition to ß(2)-agonists for optimal asthma management. Clinical trial registered with www.clinicaltrials.gov (NCT00156819, NCT00046644, and NCT00073840).

Thymic stromal lymphopoietin gene promoter polymorphisms are associated with susceptibility to bronchial asthma.

Thymic stromal lymphopoietin (TSLP) triggers dendritic cell--mediated T helper (Th) 2 inflammatory responses. A single-nucleotide polymorphism (SNP), rs3806933, in the promoter region of the TSLP gene creates a binding site for the transcription factor activating protein (AP)-1. The variant enhances AP-1 binding to the regulatory element, and increases the promoter--reporter activity of TSLP in response to polyinosinic-polycytidylic acid (poly[I:C]) stimulation in normal human bronchial epithelium (NHBE). We investigated whether polymorphisms including the SNP rs3806933 could affect the susceptibility to and clinical phenotypes of bronchial asthma. We selected three representative (i.e., Tag) SNPs and conducted association studies of the TSLP gene, using two independent populations (639 patients with childhood atopic asthma and 838 control subjects, and 641 patients with adult asthma and 376 control subjects, respectively). We further examined the effects of corticosteroids and a long-acting ß(2)-agonist (salmeterol) on the expression levels of the TSLP gene in response to poly(I:C) in NHBE. We found that the promoter polymorphisms rs3806933 and rs2289276 were significantly associated with disease susceptibility in both childhood atopic and adult asthma. The functional SNP rs3806933 was associated with asthma (meta-analysis, P = 0.000056; odds ratio, 1.29; 95% confidence interval, 1.14-1.47). A genotype of rs2289278 was correlated with pulmonary function. Moreover, the induction of TSLP mRNA and protein expression induced by poly(I:C) in NHBE was synergistically impaired by a corticosteroid and salmeterol. TSLP variants are significantly associated with bronchial asthma and pulmonary function. Thus, TSLP may serve as a therapeutic target molecule for combination therapy.

Pharmacogenetics of asthma therapy.

Asthma is a chronic disorder causing inflammation and reversible airway obstruction that affects approximately 300 million individuals worldwide. The incidence of asthma has nearly doubled in the past three decades resulting in higher rates of morbidity, mortality and health care costs. Despite the availability of several classes of asthma medications such as beta-agonists, leukotriene modifiers and corticosteroids, up to 50% of asthmatics do not benefit from one or more of these drugs. Studies have shown that asthma and drug response phenotypes such as forced expiratory volume in one second (FEV(1)) are heritable traits, indicating a genetic component of variable response to asthma drugs. This review summarizes the findings of pharmacogenetic investigations on the three main classes of asthma medications. In addition, the limitations of these genetic studies are discussed and future research avenues are proposed to identify novel genetic factors. Although numerous genes have been associated with variable response to common asthma drugs, results are often contradictory across different studies, and remain to be confirmed in larger replication cohorts. Nevertheless, literature in asthma pharmacogenetics demonstrates that genetic variants influence response to asthma treatments and may be used for predictive testing prior to drug administration to avoid adverse reactions and increase drug efficacy.

Predicting response to short-acting bronchodilator medication using Bayesian networks.

AIMS: Bronchodilator response tests measure the effect of beta(2)-agonists, the most commonly used short-acting reliever drugs for asthma. We sought to relate candidate gene SNP data with bronchodilator response and measure the predictive accuracy of a model constructed with genetic variants. MATERIALS & METHODS: Bayesian networks, multivariate models that are able to account for simultaneous associations and interactions among variables, were used to create a predictive model of bronchodilator response using candidate gene SNP data from 308 Childhood Asthma Management Program Caucasian subjects. RESULTS: The model found that 15 SNPs in 15 genes predict bronchodilator response with fair accuracy, as established by a fivefold cross-validation area under the receiver-operating characteristic curve of 0.75 (standard error: 0.03). CONCLUSION: Bayesian networks are an attractive approach to analyze large-scale pharmacogenetic SNP data because of their ability to automatically learn complex models that can be used for the prediction and discovery of novel biological hypotheses.